Coarse grained modeling of self assembled DNA 3D structure using pragmatic soft ellipsoid contact potential

Resumen Técnico: Modelado de grano grueso de estructuras 3D de ADN autoensamblado mediante un potencial de contacto de elipsoide suave pragmático

Autores: Abhirup Das y Jayashree Saha (Universidad de Calcuta).

1. El Problema

El estudio de la dinámica y el autoensamblaje del ADN es fundamental para la biología y la nanotecnología. Sin embargo, los investigadores se enfrentan a un dilema computacional: los modelos atomísticos (como AMBER o CHARMM) son extremadamente precisos pero limitados a escalas de tiempo muy cortas, mientras que los modelos de continuo carecen de los detalles moleculares necesarios. Aunque existen modelos de grano grueso (Coarse-Grained, CG), muchos utilizan geometrías esféricas para las bases nitrogenadas, lo que simplifica excesivamente la naturaleza plana y anisotrópica de estas moléculas, perdiendo información crucial sobre la orientación y las interacciones específicas de emparejamiento.

2. Metodología

Los autores proponen un nuevo modelo de grano grueso basado en el Potencial de Contacto de Elipsoide (ECP). A diferencia de los modelos esféricos, este utiliza elipsoides para representar las bases, permitiendo que la orientación esté integrada matemáticamente en la forma del potencial.

El modelo se compone de seis términos de energía ($V$):

• Interacciones enlazadas ($V_{bend}$ y $V_{dihedral}$): Modelan la flexibilidad de la columna vertebral (backbone) de azúcar-fosfato mediante potenciales armónicos y de coseno para los ángulos de flexión y torsión. Se utiliza el algoritmo RATTLE para mantener las distancias de enlace constantes.
• Interacción de volumen excluido ($V_{exc}$): Utiliza un potencial de tipo Weeks-Chandler-Anderson (WCA) para evitar el solapamiento de las partículas.
• Electroestática de largo alcance ($V_{el}$): Implementada mediante la aproximación de Debye-Hückel para modelar el efecto de la solución iónica.
• Interacción con el solvente ($V_{solv}$): Un modelo de solvente implícito mediante un potencial de Morse que actúa sobre los átomos de azúcar para representar los efectos de hidratación sin calcular explícitamente las moléculas de agua.
• Potencial de contacto de elipsoide suave ($V_{ECP}$): El núcleo del modelo. Es una extensión del potencial de Gay-Berne que permite interacciones entre partículas biaxiales no idénticas. Se divide en dos componentes:
  1. Emparejamiento (Pairing): Interacción entre bases de la cadena opuesta.
  2. Apilamiento (Stacking): Interacción entre bases adyacentes en la misma cadena.

El modelo se simuló utilizando el lenguaje Fortran 90 en un ensamble NVT con un termostato estocástico de Langevin, realizando procesos de recocido (annealing) para observar la transición de cadena sencilla (ssDNA) a doble cadena (dsDNA).

3. Contribuciones Clave

• Extensión del potencial ECP: La principal innovación es la adaptación de un potencial diseñado originalmente para lípidos a la geometría compleja y específica del ADN.
• Tratamiento de la anisotropía: Al usar elipsoides, el modelo captura de manera natural la geometría plana de las bases y cómo su orientación afecta la energía de interacción.
• Separación de interacciones: La distinción matemática entre el apilamiento (stacking) y el emparejamiento (pairing) permite un modelado más fiel de la estabilidad de la doble hélice y evita configuraciones erróneas como los "hairpins".
• Eficiencia computacional: Logra un equilibrio entre el detalle estructural y la velocidad de simulación, permitiendo estudiar procesos de renaturalización.

4. Resultados Principales

• Transición de fase: El modelo muestra una transición de fase clara durante la hibridación. El Potencial de Fuerza Media (PMF) presenta un salto de primer orden en la temperatura crítica ($T^* \approx 0.48$).
• Capacidad calorífica: Se observó un pico significativo en la capacidad calorífica específica ($C_V$) alrededor de la temperatura de transición, confirmando el cambio de fase.
• Curva de fusión (Melting Curve): El modelo reproduce con gran precisión las curvas de fusión experimentales, permitiendo una escala de energía y temperatura comparable con la realidad biológica.
• Propiedades estructurales:
  • Las distribuciones de los ángulos de flexión y diedros convergen a valores de equilibrio estables.
  • El diámetro de la hélice estimado es consistente con los valores experimentales.
  • El avance por par de bases (rise per base pair) se calculó en aproximadamente $3.677\ \text{Å}$, muy cerca del valor experimental de $3.4\ \text{Å}$.

5. Significado y Conclusión

Este trabajo demuestra que un modelo de grano grueso basado en elipses suaves es una herramienta robusta y eficaz para estudiar la biofísica del ADN. La capacidad de capturar la anisotropía de las bases sin el costo computacional de los modelos atomísticos abre nuevas vías para investigar procesos complejos como la formación de surcos (grooves), las interacciones ADN-proteína y el comportamiento del ADN en entornos de diferentes concentraciones salinas. Es un avance significativo para el diseño de dispositivos de nanotecnología basados en ADN.