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

Technische Zusammenfassung: Hydratationsfreie Energien linearer Alkane: Systematische Abweichungen in gängigen Wassermodellen und deren Korrektur

Problemstellung
Gängige molekulare Kraftfelder überschätzen die Hydratationsfreien Energien ($\Delta G_{hyd}$) hydrophober Soluaten systematisch, was in Simulationen zu einem übertriebenen hydrophoben Effekt führt. Diese Diskrepanz äußert sich in unterschätzten kritischen Mizellbildungskonzentrationen für Tenside und überschätzten freien Adsorptionsenergien an der Öl-Wasser-Grenzfläche. Während Alkane als fundamentales Modellsystem zum Verständnis der nichtpolaren Solvatation dienen, scheitern Standardparametrisierungen – insbesondere die Kombination des TraPPE-UA-Kraftfelds für Alkane mit verschiedenen Wassermodellen unter Verwendung der Standard-Lorentz-Berthelot-Mischungsregeln – daran, experimentelle Löslichkeitsdaten genau wiederzugeben. Bisherige Studien haben diese Abweichungen zwar bemerkt, doch eine systematische Korrektur, die über einen weiten Bereich von Alkankettenlängen (von Methan bis Eicosan) und für mehrere Wassermodelle anwendbar ist, fehlte.

Methodik
Die Autoren berechneten die Hydratationsfreien Energien für lineare Alkane von Methan ($C_1$) bis Eicosan ($C_{20}H_{42}$) unter Verwendung der Freien-Energie-Perturbation (FEP) im isotherm-isobaren Ensemble. Die Studie verwendete das TraPPE-UA-Vereinheitlichte-Atom-Kraftfeld für Alkane in Kombination mit vier verschiedenen Wassermodellen:

• Drei-Punkt-Modelle: SPC/E und OPC3.
• Vier-Punkt-Modelle: TIP4P/2005 und OPC.

Die Simulationen wurden mit LAMMPS unter Verwendung eines Soft-Core-Lennard-Jones-Potenzials durchgeführt, um Singularitäten während der alchemischen Transformationen zu vermeiden. Der Kopplungsparameter $\lambda$ wurde in 40 Fenstern variiert. Für den initialen Vergleich wurden die Standard-Lorentz-Berthelot-Mischungsregeln verwendet, und es wurden Langreichweitenschwanzkorrekturen für unverschobene Potenziale angewendet.

Um die systematische Überschätzung zu korrigieren, nutzten die Autoren eine Reparametrisierungsstrategie, die auf der Zerlegung der Hydratationsfreien Energie basiert:
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
wobei $\Delta G_{cavity}$ die freie Energie zur Erzeugung einer solutengroßen Kavität ist und $\Delta H_{att}$ die anziehenden Soluat-Lösungsmittel-Wechselwirkungen darstellt (approximiert als Lennard-Jones-Wechselwirkungsenergie, $\Delta U_{att}$). Durch die Nutzung etablierter Werte für $\Delta G_{cavity}$ leiteten die Autoren das Ziel-experimentelle $\Delta H_{att}$ ab und passten den Lennard-Jones-Topftiefen-Parameter ($\epsilon$) für Alkan-Wasser-Wechselwirkungen so an, dass diese Ziele bei 300 K erreicht wurden. Der Distanzparameter ($\sigma$) wurde konstant gehalten, da die Kavität-Freie Energie hochsensitiv auf $\sigma$ reagiert.

Zusätzlich evaluierte die Studie:

1. Das General Amber Force Field (GAFF) in Kombination mit den Wassermodellen TIP4P/2005 und SPC/E.
2. Das zuvor vorgeschlagene HH-Alkan-Modell (Ashbaugh et al.), erweitert auf $C_{20}$.
3. Den Einfluss der Anwendung eines verschobenen Lennard-Jones-Potenzials (wobei das Potenzial am Cutoff auf Null gesetzt wird) im Vergleich zu unverschobenen Potenzialen mit Schwanzkorrekturen.

Hauptergebnisse

• Systematische Überschätzung: Alle getesteten Wassermodelle (SPC/E, OPC3, TIP4P/2005, OPC) in Kombination mit TraPPE-UA überschätzten $\Delta G_{hyd}$ im Vergleich zu experimentellen Daten und Gruppenträger-Schätzungen. Die Abweichung war bei Drei-Punkt-Modellen ausgeprägter und nahm mit der Alkankettenlänge zu.
• Erfolg der Reparametrisierung: Eine Anpassung der Lennard-Jones-Topftiefe ($\epsilon$) für Alkan-Wasser-Wechselwirkungen um etwa 5 % (Erhöhung relativ zum Lorentz-Berthelot-Wert) führte dazu, dass die Simulationsergebnisse bei 300 K erfolgreich mit experimentellen Werten übereinstimmten.
  • Die reparametrisierten Modelle zeigten über einen Temperaturbereich von 290–350 K eine hervorragende Übereinstimmung mit experimentellen Daten.
  • Die reparametrisierten TIP4P/2005-, SPC/E- und OPC3-Modelle wiesen quadratische mittlere Abweichungen (RMSDs) von 0,07, 0,10 bzw. 0,06 $kT$ auf.
• GAFF-Leistung: Das GAFF-Kraftfeld in Kombination mit TIP4P/2005 oder SPC/E-Wasser ergab kleinere Abweichungen von experimentellen Werten im Vergleich zur ursprünglichen Kombination TraPPE-UA/TIP4P/2005, zeigte jedoch weiterhin systematische Trends (Überschätzung für kleine Alkane mit TIP4P/2005, Unterschätzung für große Alkane).
• HH-Alkan-Modell: Die Erweiterung des HH-Alkan-Modells (ursprünglich nur bis Butan validiert) auf Eicosan bestätigte seine Fähigkeit, experimentelle/Gruppenträger-Daten bis $C_{20}$ wiederzugeben.
• Potenzialverschiebung: Die Studie zeigte, dass das Verschieben des Lennard-Jones-Potenzials um seinen Wert am sphärischen Cutoff systematische Abweichungen in $\Delta G_{hyd}$ einführt. Obwohl das Hinzufügen von Schwanzkorrekturen zum verschobenen Potenzial die Ergebnisse des unverschobenen Potenzials mathematisch wiederherstellen kann, verschärft der Akt des Verschiebens an sich die Fehler, wenn er nicht rigoros korrigiert wird, was die Sensitivität alchemischer Freier Energien gegenüber Potenzial-Abschnitts-Schemata unterstreicht.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die systematische Überschätzung der hydrophoben Hydratation in gängigen Kraftfeldern ein robustes Phänomen ist, das durch eine einstufige Anpassung des Lennard-Jones-Topftiefen-Parameters für Alkan-Wasser-Wechselwirkungen korrigiert werden kann. Die Autoren betonen Folgendes:

1. Übertragbarkeit der Korrektur: Die ~5%ige Erhöhung von $\epsilon$ ist über verschiedene Wassermodelle hinweg (sowohl 3-Punkt- als auch 4-Punkt-Modelle) konsistent, was auf einen generalisierbaren Korrekturmechanismus für nichtpolare Soluaten hindeutet.
2. Grenzen der Mischungsregeln: Die Ergebnisse verdeutlichen, dass Kraftfelder, die auf Thermodynamik reiner Komponenten parametrisiert wurden, oft schlecht in Mischungen abschneiden, wenn heuristische Mischungsregeln (wie Lorentz-Berthelot) ohne spezifische Kalibrierung der Kreuzwechselwirkungen angewendet werden.
3. Methodische Einsicht: Die Studie bietet einen praktischen Weg, um Vorhersagen der Hydratationsfreien Energie zu verbessern, indem Wechselwirkungsparameter direkt gegen Freie Energien der Kavitätsbildung und experimentelle Löslichkeitsdaten kalibriert werden, anstatt sich ausschließlich auf Standard-Mischungsregeln zu verlassen.
4. Potenzialverschiebungen: Die Arbeit dient als Warnhinweis bezüglich der Verwendung verschobener Potenziale in Freie-Energie-Berechnungen und zeigt, dass solche Verschiebungen signifikante systematische Fehler in der Hydratationsthermodynamik einführen können.

Die Autoren schließen, dass genaue Hydratationsfreie Energien für nichtpolare Spezies wie Alkane erreicht werden können, indem Wechselwirkungsparameter unter Verwendung von Freien Energien der Kavitätsbildung direkt kalibriert werden, wodurch die Zuverlässigkeit molekularer Simulationen für hydrophobe Selbstassemblierung und Solvatationsphänomene verbessert wird.