Predicting the Limit of Intramolecular Hydrogen Bonding with Classical Molecular Dynamics
Molecular dynamics simulations were used to sample the conformational space (in the Figure) of flexible molecules and quantify their propensity to form intramolecular H-bonds in a variety of environments. Simulations quantitatively recapitulate experimental observables and provide insight on molecular behaviour under conditions not accessible experimentally.
Colizzi F, Hospital A, Zivanovic S, Orozco M.
Angew Chem Int Ed Engl 2019 Mar; 58: 3759.
The energetics of intramolecular recognition processes are governed by the balance of pre-organization and flexibility, which is often difficult to measure and hard to predict. Using classical MD simulations, we predict and quantify the effective strength of intramolecular hydrogen bonds between donor and acceptor sites separated by a variable alkyl linker in several solvents and crowded solutions. The balance of entropic and enthalpic contributions poses a solvent-dependent limit to the occurrence of intramolecular H-bonding. Still, free energies show a constant offset among different solvents with, for example, a 13 kJ mol(-1) difference between water and chloroform. Molecular crowding shows little effect on the thermodynamic equilibrium, but induces variations on the H-bond kinetics. The results are in quantitative agreement with experiments in chloroform and showcase a general strategy to investigate molecular interactions in different environments, extending the limits of current experiments towards the prospective prediction of H-bond interactions in a variety of contexts.
PubMed: 30681249. Doi: 10.1002/anie.201810922.
