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

Sintesi Tecnica: Energie Libere di Idratazione degli Alcani Lineari: Deviazioni Sistematiche nei Modelli d'Acqua Comuni e Loro Correzione

Enunciato del Problema
I campi di forza molecolari comuni sovrastimano sistematicamente le energie libere di idratazione ($\Delta G_{hyd}$) dei soluti idrofobici, portando a un effetto idrofobico esagerato nelle simulazioni. Questa discrepanza si manifesta come una sottostima delle concentrazioni micellari critiche per i tensioattivi e una sovrastima delle energie libere di adsorbimento interfacciale olio-acqua. Sebbene gli alcani costituiscano un sistema modello fondamentale per comprendere la solvatazione non polare, le parametrizzazioni standard — in particolare la combinazione del campo di forza TraPPE-UA per gli alcani con vari modelli d'acqua e le regole di miscelazione standard di Lorentz-Berthelot — non riescono a riprodurre accuratamente i dati sperimentali di solubilità. Studi precedenti hanno notato queste deviazioni, ma è mancata una correzione sistematica applicabile a un'ampia gamma di lunghezze della catena degli alcani (dal metano all'eicosano) e a molteplici modelli d'acqua.

Metodologia
Gli autori hanno calcolato le energie libere di idratazione per alcani lineari che variano dal metano ($C_1$) all'eicosano ($C_{20}H_{42}$) utilizzando la Perturbazione dell'Energia Libera (FEP) nell'insieme isotermo-isobaro. Lo studio ha impiegato il campo di forza united-atom TraPPE-UA per gli alcani combinato con quattro distinti modelli d'acqua:

• Modelli a tre siti: SPC/E e OPC3.
• Modelli a quattro siti: TIP4P/2005 e OPC.

Le simulazioni sono state condotte utilizzando LAMMPS con un potenziale di Lennard-Jones a nucleo morbido per evitare singolarità durante le trasformazioni alchemiche. Il parametro di accoppiamento $\lambda$ è stato variato in 40 finestre. Lo studio ha utilizzato le regole di miscelazione standard di Lorentz-Berthelot per i confronti iniziali e ha applicato correzioni di coda a lungo raggio per potenziali non spostati.

Per correggere la sovrastima sistematica, gli autori hanno utilizzato una strategia di riparametrizzazione basata sulla decomposizione dell'energia libera di idratazione:
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
dove $\Delta G_{cavity}$ è l'energia libera di creazione di una cavità delle dimensioni del soluto e $\Delta H_{att}$ rappresenta le interazioni attrattive soluto-solvente (approssimate come l'energia di interazione di Lennard-Jones, $\Delta U_{att}$). Utilizzando valori consolidati per $\Delta G_{cavity}$, gli autori hanno derivato il $\Delta H_{att}$ sperimentale target e hanno regolato il parametro di profondità del pozzo di Lennard-Jones alcani-acqua ($\epsilon$) per corrispondere a questi target a 300 K. Il parametro di distanza ($\sigma$) è stato mantenuto costante, poiché le energie libere di cavità sono altamente sensibili a $\sigma$.

Inoltre, lo studio ha valutato:

1. Il campo di forza GAFF (General Amber Force Field) combinato con i modelli d'acqua TIP4P/2005 e SPC/E.
2. Il modello HH-alkane precedentemente proposto (Ashbaugh et al.) esteso fino a $C_{20}$.
3. L'impatto dell'applicazione di un potenziale di Lennard-Jones spostato (dove il potenziale è azzerato al cutoff) rispetto ai potenziali non spostati con correzioni di coda.

Risultati Chiave

• Sovrastima Sistematica: Tutti i modelli d'acqua testati (SPC/E, OPC3, TIP4P/2005, OPC) combinati con TraPPE-UA hanno sovrastimato $\Delta G_{hyd}$ rispetto ai dati sperimentali e alle stime basate sul contributo di gruppo. La deviazione era più pronunciata per i modelli a tre siti e aumentava con la lunghezza della catena degli alcani.
• Successo della Riparametrizzazione: La regolazione della profondità del pozzo di Lennard-Jones alcani-acqua ($\epsilon$) di circa il 5% (aumentandola rispetto al valore di Lorentz-Berthelot) ha allineato con successo i risultati delle simulazioni ai valori sperimentali a 300 K.
  • I modelli riparametrizzati hanno mantenuto un eccellente accordo con i dati sperimentali in un intervallo di temperature da 290 a 350 K.
  • I modelli TIP4P/2005, SPC/E e OPC3 riparametrizzati hanno mostrato deviazioni quadratiche medie (RMSD) di 0,07, 0,10 e 0,06 $kT$, rispettivamente.
• Prestazioni di GAFF: Il campo di forza GAFF combinato con i modelli d'acqua TIP4P/2005 o SPC/E ha prodotto deviazioni più piccole dai valori sperimentali rispetto alla combinazione originale TraPPE-UA/TIP4P/2005, sebbene abbia ancora mostrato tendenze sistematiche (sovrastima per alcani piccoli con TIP4P/2005, sottostima per alcani grandi).
• Modello HH-Alkane: L'estensione del modello HH-alkane (originariamente validato solo fino al butano) all'eicosano ha confermato la sua capacità di riprodurre i dati sperimentali/basati sul contributo di gruppo fino a $C_{20}$.
• Spostamento del Potenziale: Lo studio ha dimostrato che lo spostamento del potenziale di Lennard-Jones del suo valore al cutoff sferico introduce deviazioni sistematiche in $\Delta G_{hyd}$. Sebbene l'aggiunta di correzioni di coda al potenziale spostato possa recuperare matematicamente i risultati non spostati, l'atto dello spostamento stesso esacerba gli errori se non corretto rigorosamente, evidenziando la sensibilità delle energie libere alchemiche agli schemi di troncamento del potenziale.

Significato e Affermazioni
Il documento afferma che la sovrastima sistematica dell'idratazione idrofobica nei campi di forza comuni è un fenomeno robusto che può essere corretto attraverso un aggiustamento in un solo passo del parametro di profondità del pozzo di Lennard-Jones alcani-acqua. Gli autori sottolineano che:

1. Trasferibilità della Correzione: L'aumento di circa il 5% di $\epsilon$ è coerente tra diversi modelli d'acqua (sia a 3 siti che a 4 siti), suggerendo un meccanismo di correzione generalizzabile per soluti non polari.
2. Limiti delle Regole di Miscelazione: I risultati evidenziano che i campi di forza parametrizzati sulla termodinamica dei componenti puri spesso performano male nelle miscele quando vengono applicate regole euristica di miscelazione (come quelle di Lorentz-Berthelot) senza una specifica calibrazione delle interazioni incrociate.
3. Insight Metodologico: Lo studio fornisce una via pratica per migliorare le previsioni delle energie libere di idratazione calibrando direttamente i parametri di interazione contro le energie libere di formazione della cavità e i dati sperimentali di solubilità, piuttosto che affidarsi esclusivamente alle regole di miscelazione standard.
4. Spostamenti del Potenziale: Il lavoro funge da nota cautelativa riguardo all'uso di potenziali spostati nei calcoli delle energie libere, dimostrando che tali spostamenti possono introdurre errori sistematici significativi nella termodinamica dell'idratazione.

Gli autori concludono che energie libere di idratazione accurate per specie non polari come gli alcani possono essere ottenute calibrando direttamente i parametri di interazione utilizzando le energie libere di formazione della cavità, migliorando così l'affidabilità delle simulazioni molecolari per l'autoassemblaggio idrofobico e i fenomeni di solvatazione.