Hydration Free Energies of Linear Alkanes: Systematic Deviations in Common Water Models and Their Correction

Resumen Técnico: Energías Libres de Hidratación de Alcanos Lineales: Desviaciones Sistemáticas en Modelos de Agua Comunes y su Corrección

Planteamiento del Problema
Los campos de fuerza moleculares comunes sobreestiman sistemáticamente las energías libres de hidratación ($\Delta G_{hyd}$) de solutos hidrofóbicos, lo que conduce a un efecto hidrofóbico exagerado en las simulaciones. Esta discrepancia se manifiesta como concentraciones micelares críticas subestimadas para tensioactivos y energías libres de adsorción en la interfaz aceite-agua sobreestimadas. Si bien los alcanos sirven como un sistema modelo fundamental para comprender la solvatación no polar, las parametrizaciones estándar —específicamente la combinación del campo de fuerza TraPPE-UA para alcanos y diversos modelos de agua con reglas de mezcla Lorentz–Berthelot estándar— no logran reproducir con precisión los datos experimentales de solubilidad. Estudios anteriores han notado estas desviaciones, pero ha faltado una corrección sistemática aplicable a una amplia gama de longitudes de cadena de alcano (desde metano hasta eicosano) y múltiples modelos de agua.

Metodología
Los autores calcularon las energías libres de hidratación para alcanos lineales que van desde el metano ($C_1$) hasta el eicosano ($C_{20}H_{42}$) utilizando Perturbación de Energía Libre (FEP) dentro del ensemble isotérmico-isobárico. El estudio empleó el campo de fuerza unido-atómico TraPPE-UA para alcanos combinado con cuatro modelos de agua distintos:

• Modelos de tres sitios: SPC/E y OPC3.
• Modelos de cuatro sitios: TIP4P/2005 y OPC.

Las simulaciones se realizaron utilizando LAMMPS con un potencial de Lennard-Jones de núcleo blando para evitar singularidades durante las transformaciones alquímicas. El parámetro de acoplamiento $\lambda$ se varió en 40 ventanas. El estudio utilizó reglas de mezcla Lorentz–Berthelot estándar para comparaciones iniciales y aplicó correcciones de cola de largo alcance para potenciales no desplazados.

Para corregir la sobreestimación sistemática, los autores utilizaron una estrategia de reparametrización basada en la descomposición de la energía libre de hidratación:
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
donde $\Delta G_{cavity}$ es la energía libre de crear una cavidad del tamaño del soluto y $\Delta H_{att}$ representa las interacciones atractivas soluto-solvente (aproximadas como la energía de interacción de Lennard-Jones, $\Delta U_{att}$). Al utilizar valores establecidos para $\Delta G_{cavity}$, los autores derivaron el $\Delta H_{att}$ experimental objetivo y ajustaron el parámetro de profundidad del pozo de Lennard-Jones para la interacción alcano-agua ($\epsilon$) para igualar estos objetivos a 300 K. El parámetro de distancia ($\sigma$) se mantuvo constante, ya que las energías libres de cavidad son altamente sensibles a $\sigma$.

Adicionalmente, el estudio evaluó:

1. El Campo de Fuerza General Amber (GAFF) combinado con los modelos de agua TIP4P/2005 y SPC/E.
2. El modelo HH-alkane propuesto anteriormente (Ashbaugh et al.) extendido hasta $C_{20}$.
3. El impacto de aplicar un potencial de Lennard-Jones desplazado (donde el potencial se anula en el corte) frente a potenciales no desplazados con correcciones de cola.

Resultados Clave

• Sobreestimación Sistemática: Todos los modelos de agua probados (SPC/E, OPC3, TIP4P/2005, OPC) combinados con TraPPE-UA sobreestimaron $\Delta G_{hyd}$ en relación con los datos experimentales y las estimaciones de contribución de grupos. La desviación fue más pronunciada para los modelos de tres sitios y aumentó con la longitud de la cadena de alcano.
• Éxito de la Reparametrización: Ajustar la profundidad del pozo de Lennard-Jones para la interacción alcano-agua ($\epsilon$) en aproximadamente un 5% (aumentándola en relación con el valor Lorentz–Berthelot) alineó con éxito los resultados de la simulación con los valores experimentales a 300 K.
  • Los modelos reparametrizados mantuvieron un excelente acuerdo con los datos experimentales en un rango de temperatura de 290–350 K.
  • Los modelos TIP4P/2005, SPC/E y OPC3 reparametrizados mostraron desviaciones cuadráticas medias (RMSD) de 0.07, 0.10 y 0.06 $kT$, respectivamente.
• Rendimiento de GAFF: El campo de fuerza GAFF combinado con agua TIP4P/2005 o SPC/E produjo desviaciones menores de los valores experimentales en comparación con la combinación original TraPPE-UA/TIP4P/2005, aunque aún exhibió tendencias sistemáticas (sobreestimación para alcanos pequeños con TIP4P/2005, subestimación para alcanos grandes).
• Modelo HH-Alcano: La extensión del modelo HH-alkane (originalmente validado solo hasta butano) hasta eicosano confirmó su capacidad para reproducir datos experimentales/de contribución de grupos hasta $C_{20}$.
• Desplazamiento de Potenciales: El estudio demostró que desplazar el potencial de Lennard-Jones por su valor en el corte esférico introduce desviaciones sistemáticas en $\Delta G_{hyd}$. Si bien añadir correcciones de cola al potencial desplazado puede recuperar matemáticamente los resultados no desplazados, el acto de desplazar en sí mismo exacerba los errores si no se corrige rigurosamente, destacando la sensibilidad de las energías libres alquímicas a los esquemas de truncamiento de potenciales.

Significado y Afirmaciones
El artículo afirma que la sobreestimación sistemática de la hidratación hidrofóbica en campos de fuerza comunes es un fenómeno robusto que puede corregirse mediante un ajuste de un solo paso del parámetro de profundidad del pozo de Lennard-Jones para la interacción alcano-agua. Los autores enfatizan que:

1. Transferibilidad de la Corrección: El aumento de ~5% en $\epsilon$ es consistente entre diferentes modelos de agua (tanto de 3 sitios como de 4 sitios), lo que sugiere un mecanismo de corrección generalizable para solutos no polares.
2. Limitaciones de las Reglas de Mezcla: Los resultados destacan que los campos de fuerza parametrizados en termodinámica de componentes puros a menudo funcionan mal en mezclas cuando se aplican reglas de mezcla heurísticas (como Lorentz–Berthelot) sin una calibración específica de las interacciones cruzadas.
3. Perspectiva Metodológica: El estudio proporciona una vía práctica para mejorar las predicciones de energía libre de hidratación mediante la calibración directa de los parámetros de interacción contra las energías libres de formación de cavidades y los datos experimentales de solubilidad, en lugar de depender exclusivamente de reglas de mezcla estándar.
4. Desplazamientos de Potencial: El trabajo sirve como una nota de precaución sobre el uso de potenciales desplazados en cálculos de energía libre, demostrando que tales desplazamientos pueden introducir errores sistemáticos significativos en la termodinámica de hidratación.

Los autores concluyen que se pueden lograr energías libres de hidratación precisas para especies no polares como los alcanos mediante la calibración directa de los parámetros de interacción utilizando energías libres de formación de cavidades, mejorando así la fiabilidad de las simulaciones moleculares para fenómenos de autoensamblaje hidrofóbico y solvatación.