Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular...
Transcript of Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular...
![Page 1: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/1.jpg)
Dinámica Molecular de ProteínasDinámica Molecular de ProteínasModelado y Simulación ComputacionalModelado y Simulación Computacional
Profesores: Eliana K. Asciutto & Ignacio J. General2do cuatrimestre 2017
Escuela de Ciencia y TecnologíaUNSAM
![Page 2: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/2.jpg)
Dinámica Molecular de ProteínasDinámica Molecular de ProteínasModelado y Simulación ComputacionalModelado y Simulación Computacional
Dinámica MolecularDinámica Molecular
![Page 3: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/3.jpg)
3
¿Qué es la Dinámica Molecular (DM)?
Dinámica Molecular – UNSAM – 2017
Es el estudio de los movimientos de las moléculas. Algunas preguntas que puede responder:
● ¿Cómo son sus fluctuaciones estructurales? ● ¿Existen correlaciones entre los movimientos de distintas
partes de la molécula? ● ¿Cómo es el paisaje de energía libre?● ¿Se liga a otra molécula?● ¿Cuáles son los efectos del solvente?
![Page 4: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/4.jpg)
4
¿Qué es la Dinámica Molecular (DM)?
Dinámica Molecular – UNSAM – 2017
Algunas aplicaciones de DM a proteínas:
● Refinamiento de estructuras obtenidas por RMN, rayos X, etc● Predicción de estructuras a partir de secuencias● Cambios conformacionales● Termodinámica (ligamiento)
Energía
Entropía
Glóbulo fundido
Estado nativo
Estado desplegado
By Thomas Splettstoesser (www.scistyle.com) (Own work) [CC BY-SA 3.0 (http://creativecommons.org/licenses/by-sa/3.0)], via Wikimedia Commons
![Page 5: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/5.jpg)
5
¿Que es la Dinámica Molecular (DM)?
Dinámica Molecular – UNSAM – 2017
¿Cómo hacer una DM en una computadora?
Este es un problema conocido de la física: problema de n-cuerposproblema de n-cuerpos
No se puede resolver en forma exacta para más de 2 cuerpos (las proteínas tienen miles de átomos), pero si se pueden plantear las ecuaciones exactas y resolverlas en forma numérica, aproximada, usando computadoras.
ProcedimientoProcedimiento (para dinámica clásica):➔ Definir las interacciones entre átomos:
potencial U (resortes, Coulomb, etc)➔ Definir un conjunto de posiciones y velocidades iniciales para
todas las partículas➔ Aplicar las leyes de Newton a cada partícula para calcular su
nuevas posiciones y velocidades
![Page 6: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/6.jpg)
6
Algoritmo básico de DM
Dinámica Molecular – UNSAM – 2017
Dada la energía potencial, U(x), del sistema a estudiar:
x1, v
1condiciones iniciales
F = -∂U/∂x fuerza a partir de U
a = F/m 2da ley de Newton
v2 = v
1+a Δt cinemática del movimiento acelerado
x2 = x
1+v
2 Δt+O(Δt2) cinemática del movimiento acelerado
F
v1
v2F→a→Δv→Δx
Δx
![Page 7: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/7.jpg)
7
Ventajas de DM
Dinámica Molecular – UNSAM – 2017
Ventajas del planteo de DM:Ventajas del planteo de DM:● Solo es necesario conocer las interacciones interatómicas ● No es necesario presuponer nada sobre los procesos
interatómicos
Ventajas en cuanto a resultados:Ventajas en cuanto a resultados:● Provee detalles a nivel microscópico
Puede descubrir nuevos procesos o mecanismos físico/biológicos
![Page 8: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/8.jpg)
8
Limitaciones de DM
Dinámica Molecular – UNSAM – 2017
Descripción correcta→mecánica cuántica (Schrödinger)Pero su resolución es computacionalmente muy costosa
¿En qué casos se pueden usar ecuaciones clásicas?Cuando la longitud de onda de los átomos es menor que las Cuando la longitud de onda de los átomos es menor que las distancias típicas del sistema (distancias típicas del sistema (λ<λ<dd))
Usando la energía térmica efectiva de partícula libre, E ~ π KBT :
E=p2
2m⇒ p=√2m E →
hλ
=√2 m E ⇒ λ=h
√2m E
de Broglie
λt=h
√2 mπ K B T
Relación λ vs E para una partícula
![Page 9: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/9.jpg)
9
Limitaciones de DM
Dinámica Molecular – UNSAM – 2017
Para T = 300 K, y con un d típico ~ 2Å:
¿Se puede aproximar al H para ser tratado clásicamente?➔ Métodos SHAKE/RATTLEMétodos SHAKE/RATTLE: fijan los enlaces de H
(estos métodos remueven grados de libertad de H; este es justamente el problema, ya que dichos GL no están activos cuánticamente, pero si clásicamente)
➔ SHAKE y RATTLE también van a permitir aumentar el paso de aumentar el paso de tiempo (tiempo (Δt, Δt, time steptime step)), ya que la vibración de H es el movimiento más rápido del sistema.
M [u] λt [Å]
H 1 1.0
C 12 0.3
K 39 0.2
Excepto los átomos más livianos (H, He,...), es una buena aproximación considerar a los demás como partículas puntuales clásicas.
![Page 10: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/10.jpg)
10
Limitaciones de DM
Dinámica Molecular – UNSAM – 2017
Escalas de tiempo:
fs
ps
ns
μs
ms
s
Vibración enlaces
Rotación cadenas laterales
Movimientos bisagra
Transiciones alostéricasPlegamiento
El Δt de MD debería ser menor que el menor de estos tiempos
Δt = 1 fso
Δt = 2 fs usando SHAKE
Ligamiento sustratos
Computadoras modernas (clusters, GPUs):
(Δt = 2 fs) × 500.000.000 pasos ~ 1 μμss
![Page 11: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/11.jpg)
11
Limitaciones de DM
Dinámica Molecular – UNSAM – 2017
Escalas de distancia:
Número de interacciones en un sistema = n×n
Típicamente, en sistemas de proteína + solvente:
~ 10.000 - 1.000.000 átomos~ celdas cúbicas de 40 - 200 Å (1Å = 0.1 nm)
Esto es lo que determina el costo de la simulación
0
1 nm
10 nm
100 nm
1 μm
10 μm
1 fs 1 ps 1 ns 1 μs
l
tAb initio
Continuo
MD atómica
ASICs, GPUs
![Page 12: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/12.jpg)
12
Dinámica Molecular – UNSAM – 2017
AM
BE
R b
enc
hm
arks
http
://a
mbe
rmd.
org/
gpus
/ben
chm
arks
.htm
#Ben
chm
arks
![Page 13: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/13.jpg)
13
Limitaciones de DM
Dinámica Molecular – UNSAM – 2017
Escalas de distancia:
Número de interacciones en un sistema = n×n
Típicamente, en sistemas de proteína + solvente:
~ 10.000 - 1.000.000 átomos
~ celdas cúbicas de 40 - 200 Å
El tamaño finito de la celda puede introducir perturbaciones (artifacts). Por eso conviene tener una distancia extra entre la proteína y la pared de la celda.
Esto es lo que determina el costo de la simulación
![Page 14: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/14.jpg)
14
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
Tipos de energía potencial (U)
Los sistemas físicos en equilibrio tienden a minimizar su energía
La conformación de mínima energía será el estado nativoLa conformación de mínima energía será el estado nativo
(en principio, ya que al haber temperatura—interacciones entrópicas—el sistema es corrido del equilibrio)
¿Y como es la expresión de la energía?
Interacciones entálpicas:➢ de átomos ligados➢ de átomos no ligados
Interacciones entrópicas:➢ energía térmica
![Page 15: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/15.jpg)
15
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
U = Ulineal + Uangular + Udiedro + Ucoulomb + Uvdw
covalente2 átomos covalente
3 átomos
covalente4 átomos
no covalenteCoulomb
no covalenteVan der Waals
![Page 16: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/16.jpg)
16
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
Φ0
Θ0
exclusión de Pauli(repulsión)
fuerzas de dispersión(atractivas)
(o Lennard-Jones)
![Page 17: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/17.jpg)
17
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
n átomos → n2 interacciones
Ejemplo) n =100.000 átomos → 10.000.000.000 interacciones por paso→ extremadamente costoso
¿¿Cómo reducir el costo?Cómo reducir el costo?
![Page 18: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/18.jpg)
18
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
¿¿Cómo reducir el costo?Cómo reducir el costo?
1)Dejar de contar interacciones con Δrij > Rcutoff
➔ El costo cae a ~nRcutoff3
➔ Pero esto no es físicamente correcto: las interacciones electrostáticas de-caen con la distancia, pero hay más→se compensan
2)Particle Mesh Ewald (PME) permite, usando transformaciones de Fourier, sumar todas la interacciones de lar-go alcance a cambio de imponer condiciones periódicas de contorno (red cristalina). El costo es nlog(n).
![Page 19: Dinámica Molecular de Proteínas Modelado y Simulación ... · 3 ¿Qué es la Dinámica Molecular (DM)? Dinámica Molecular – UNSAM – 2017 Es el estudio de los movimientos de](https://reader030.fdocument.pub/reader030/viewer/2022040207/5e11807e58652f24e8517253/html5/thumbnails/19.jpg)
19
Campos de fuerza (Force Fields, FF)
Dinámica Molecular – UNSAM – 2017
¿¿Cómo reducir el costo?Cómo reducir el costo?
3)SHAKE/RATTLE: enlaces de hidrógeno constantes:Δt: 1fs→2fs
4)Fuerzas electrostáticas de largo alcance varían lentamente: no es necesario evaluarlas en cada paso de tiempo, hacerlo cada 2 o 3 pasos.