1 / 31

Bioinformática estructural

Bioinformática estructural. Principios y aplicaciones de la Biofisicoquímica molecular. Dr. Marcelo A. Martí Grupo de Modelado Molecular DQIAyQF Grupo de Bioinformática Estructural Dto.QB Facultad de Ciencias Exactas y Naturales Universidad de Buenos Aires

havyn
Download Presentation

Bioinformática estructural

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Bioinformática estructural Principios y aplicaciones de la Biofisicoquímica molecular Dr. Marcelo A. Martí Grupo de Modelado Molecular DQIAyQF Grupo de Bioinformática Estructural Dto.QB Facultad de Ciencias Exactas y Naturales Universidad de Buenos Aires marcelo@qi.fcen.uba.ar, marcelomarti@yahoo.com

  2. ¿Que es la Bioinformática estructural? • Bioinformática: Desarrollo y aplicación de métodos que convierten datos biológicos en conocimiento sobre los sistemas • 1) Desarrollo: teorías y algoritmos • 2) Aplicaciones: análisis-predicción • 3) Organización de datos: Arrays, Mass Spectrometry, DNA seqs • Bioinformática estructural aplicación de 1,2 y 3 a la estructura de biomoléculas. 1 y 2 => química computacional 3 bases de datos estructurales

  3. Z Y X ¿Qué es la Estructura? “Conocer la estructura implica conocer la posición espacial relativa de todos los átomos de la proteína”

  4. Otras Representaciones: Resalta la estructura adoptada por el “Backbone”

  5. Otras Representaciones: Representación del potencial electroestático superficial Representación Volumétrica

  6. Métodos de Simulación ComputacionalEn Química y Biología Herramientas fisicoquímicas implementadas en programas 1) Modelos: recetas para evaluar superficie energía potencial (SEP). (potenciales clásicos, QM, etc) 2) Esquemas para extraer información de SEP. Monte Carlo, Dinámica Molecular.

  7. MODELOS CUANTICOS • Resuelven (o tratan!) ecuación Schrödinger: H = E , Energía es función de la configuración electrónica y de las posiciones nucleares • Existen diferentes niveles de teoría (HF, DFT, MP2, CI, Semiempiricos) • Indispensable para: procesos reactivos, estados excitados. • En principio: válido para cualquier sistema. • En la práctica: limitaciones computacionales en cuanto al tamaño del sistema.

  8. MODELOS CLASICOS • Modelos basados en la ideas clásicas: • (ley Coulomb, resortes, etc) E = f(N:x,y,z) • Distintos niveles aproximación según complejidad • función propuesta, (atomísticos, grano grueso). • En principio: no son válidos para procesos reactivos. • Parametrización es clave y requiere artesanía! • Son los más utilizados para Biomoléculas!!

  9. Campos de Fuerza Clásicos¿cómo es un campo de fuerza? Términos de unión Términos de NO unión

  10. E Q3 V(Q2)-V(Q1) Q1 Q2 Q ¿Como obtengo la información relevante? Supongamos que podemos crear MUCHAS conformaciones “R” y leer E(R) para cada una puedo crear una superficie de energía potencial (SEP) Propiedades del sístema

  11. ¿Cómo obtengo la SEP? Proyección temporal de las coordenadas integrando las ecuaciones de Newton => TRAYECTORIA x(t) => x(t + dt) => x(t + 2dt)….. x(t + dt) = x(t) + vdt + ½ (F(t)/m) dt2 F = E/x Alternativa: Muestreo de la SEP x Monte Carlo

  12. Monte Carlo • A partir de una configuración de partida, se generan configuraciones nuevas “al azar” • Cada nueva configuración se evalúa p = exp(-E)/kT), donde E, es la diferencia de energía entre ambas config. • Si p > 1, se acepta la movida. (E < 0) • Si p< 1 (E > 0), se compara p con otro número al azar entre 0 y 1, X. Si p < X, se acepta la movida.

  13. ¿Qué quiere decir aceptar movida? “A” propiedad del sistema y Riuna confromación • Si acepto movida, uso esa configuración para la suma. Si rechazo movida, vuelvo a poner en la suma la configuración anterior. • Actualmente MC se utiliza como algoritmo eficiente de búsqueda en el espacio configuracional

  14. ¿Qué se obtiene de una simulación? Confirmación de las Hipótesis sobre el comportamiento FQ de los sistemas estudiados Entender determinantes moleculares de un fenómeno conocido. Obtener valor agregado: poder lograr información no accesible experimentalmente

  15. Un poquito más de FQ: “E vs G” • Los métodos QM como los de MD nos permiten obtener la Energía “E” directamente pero • El comportamiento de las Biomoléculas en solución esta determinado por la Energía Libre (G o A) G = H -T S A = E -T S En fases condensadas H  E, y G  A G = RT ln (Keq)

  16. Ejemplo 1 Estudio Dinámico Estructural de la Neuroglobina Humana

  17. La Neuroglobina Grupo Hemo Función: desconocida Control redox ?

  18. -O2 +O2 Mecanismo de regulación: 5c Oxy 6c Equilibrio 6c-5c determina la afinidad!! ¿Qué hicimos? • Dinámica Ngb 6c, 5c, Oxy, 6cS-S, 5cS-S • Perfil de energía libre 5c(S-S) <=> 6c(S-S)

  19. Comparación estructural RMSD clustering: RMSD(Å): 6c SH Oxy 5c SH 6c S-S 5c S-S 6c vs 5c 1.96 6c vs Oxy 1.42 6c vs 5c (S-S) 0.73

  20. S-S SH ------ Equilibrio 6c-5c

  21. Explicación Molecular 6c SH, 5c SH, Oxy 6cSS ---, 5cSS ---

  22. Ejemplo 2 Mecanismos de regulación de la afinidad por oxigeno en hemoproteínas x Métodos QM-MM

  23. Ejemplo: Calculo de energía de unión oxígeno en hemoproteínas Eu = Eoxy-(Edeoxy-EO2) Koff Eu Métodos Híbridos QM/MM • Juntan lo bueno y lo malo de ambos mundos • Precisión y relatividad en el sitio activo (QM) + efecto realístico del entorno proteico (MM) • Estudio de reacciones Químicas en sistemas complejos

  24. HisE7 ValE11 TyrB10 GlnE7 HisF8 GlnE11 HisF8 Relación estructura-afinidad

  25. Combinando MD con QM/MM Sitio Activo de la hemoglobina de Cerebratulus lacteus • Afinidad por O2 moderada • La mutación de TyrB10 reduce la afinidad 3 veces • La mutación de ThrE11 aumenta la afinidad 1000 veces

  26. Dinámica Molecular CerHb 2 conformaciones para la Tyr B10 ΔE O2 24.4 kcal/mol ΔE O2 32.8 kcal/mol Perfil de energía libre de ínter-conversión

  27. Conclusión: Efectos dinámicos sobre la afinidad *Sólo postulando un equilibrio dinámico se explican los resultados para la wt y las dos mutantes

  28. Ejemplo 3 Catálisis Enzimática x Métodos QM-MM

  29. Reacciones enzimáticos (x métodos QM-MM) Sitio activo: QM proteina: MM

  30. Reaccion de la corismato mutasa Prefenato Corismato

  31. Resulado: E Kcal/mol DE‡= 13.8 DE‡= 5.3 D(DE‡) = 8.5 D(DE‡exp) = 8 Reactivo DE = -22.0 DE = -24.9 Producto

More Related