Hacia las verdaderas energías de apilamiento de bases de ADN: MP2, CCSD (T), y cálculos completos del conjunto de...

Área de la Ciencia:

• Química computacional es la química computacional.
• El modelado biomolecular.
• Química cuántica es la química cuántica.

Sus antecedentes:

• Las interacciones de apilamiento de nucleobases son fundamentales para la estructura y función del ADN y el ARN.
• La cuantificación precisa de estas energías de apilamiento es esencial para comprender la estabilidad e interacciones del material genético.
• Los métodos computacionales anteriores a menudo subestimaban las energías de apilamiento, lo que requería cálculos más precisos.

Objetivo del estudio:

• Para determinar con precisión las energías de apilamiento de los homodímeros de pirimidina, uracilo, citosina y guanina.
• Para evaluar el impacto de las correcciones de correlación de orden superior en los cálculos de la energía de apilamiento.
• Proporcionar valores confiables de energía de apilamiento para la calibración de campos de fuerza y otros métodos computacionales.

Principales métodos:

• Utilizó la teoría de perturbaciones de Møller-Plesset de segundo orden (MP2) y el cúmulo acoplado con cálculos individuales y dobles y triples perturbadores (CCSD ((T)).
• Empleó una serie de conjuntos de bases de valencia dividida, consistentes en la correlación y de funciones de bonos.
• Extrapolado al límite del conjunto de base completo (CBS) para los cálculos de MP2 (CBS MP2) y las correcciones CCSD (T) aplicadas.

Principales resultados:

• Las energías de apilamiento del conjunto de bases completas MP2 (CBS MP2) fueron de -5.3, -12.3 y -11.2 kcal/mol para los dímeros de pirimidina, uracilo y citosina, respectivamente.
• Las correcciones de correlación de orden superior (CCSD) refinaron estos valores a -3.4, -11.6 y -10.4 kcal/mol.
• La energía de apilamiento del dímero de guanina se estimó en -12.9 kcal/mol después de la corrección CCSD (T), con una estimación de CBS MP2 de -14.8 kcal/mol.

Conclusiones:

• Las energías de apilamiento calculadas proporcionan puntos de referencia altamente precisos para el apilamiento de bases de ADN y ARN.
• El método MP2/6-31G ((0.25) subestima sistemáticamente las energías de apilamiento, lo que requiere una corrección para la calibración del campo de fuerza.
• Los valores precisos de energía de apilamiento son cruciales para mejorar la confiabilidad de las simulaciones moleculares y comprender las interacciones biomoleculares.