Energy controlled insertion of polar molecules in dense fluids.

We present a method to search low energy configurations of polar molecules in the complex potential energy surfaces associated with dense fluids. The search is done in the configurational space of the translational and rotational degrees of freedom of the molecule, combining steepest-descent and Newton-Raphson steps which embed information on the average sizes of the potential energy wells obtained from prior inspection of the liquid structure. We perform a molecular dynamics simulation of a liquid water shell which demonstrates that the method enables fast and energy-controlled water molecule insertion in aqueous environments. The algorithm finds low energy configurations of incoming water molecules around three orders of magnitude faster than direct random insertion. This method represents an important step towards dynamic simulations of open systems and it may also prove useful for energy-biased ensemble average calculations of the chemical potential.

[1]  P. V. Coveney,et al.  USHER: An algorithm for particle insertion in dense fluids , 2003 .

[2]  D. J. Adams,et al.  Grand canonical ensemble Monte Carlo for a Lennard-Jones fluid , 1975 .

[3]  B. Montgomery Pettitt,et al.  Dynamic simulations of water at constant chemical potential , 1992 .

[4]  K. Schulten,et al.  Electrostatic tuning of permeation and selectivity in aquaporin water channels. , 2003, Biophysical journal.

[5]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[6]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[7]  P. Coveney,et al.  Comparison of molecular dynamics with hybrid continuum-molecular dynamics for a single tethered polymer in a solvent. , 2004, The Journal of chemical physics.

[8]  Eirik Grude Flekkøy,et al.  Hybrid model for combined particle and continuum dynamics , 2000 .

[9]  S. A. Dodds,et al.  Chemical Physics , 1877, Nature.

[10]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[11]  K. Shing,et al.  The chemical potential in dense fluids and fluid mixtures via computer simulation , 1982 .

[12]  M. Levitt,et al.  Computer simulation of protein folding , 1975, Nature.

[13]  L. Scriven,et al.  Efficient molecular simulation of chemical potentials , 1989 .

[14]  H. Scheraga,et al.  Monte Carlo-minimization approach to the multiple-minima problem in protein folding. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[15]  David A. Kofke,et al.  Appropriate methods to combine forward and reverse free-energy perturbation averages , 2003 .

[16]  Gillian C. Lynch,et al.  Semi-grand canonical molecular dynamics simulation of bovine pancreatic trypsin inhibitor , 2000 .

[17]  J. Thornton,et al.  Modelling protein unfolding: A solvent insertion protocol , 1996 .

[18]  U. Bornscheuer,et al.  Protein Engineering , 2018, Methods in Molecular Biology.

[19]  A. Pohorille,et al.  Excess chemical potential of small solutes across water--membrane and water--hexane interfaces. , 1996, The Journal of chemical physics.

[20]  M. Saunders.  Stochastic exploration of molecular mechanics energy surfaces. Hunting for the global minimum , 1987 .

[21]  M. Mezei,et al.  Calculation of the Free Energy Profile of H2O, O2, CO, CO2, NO, and CHCl3 in a Lipid Bilayer with a Cavity Insertion Variant of the Widom Method , 2000 .

[22]  Mihaly Mezei,et al.  Grand-canonical ensemble Monte Carlo study of dense liquid Lennard-Jones, soft spheres and water , 1987 .

[23]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

[24]  Y. Guissani,et al.  A computer‐simulation study of hydrophobic hydration of rare gases and of methane. I. Thermodynamic and structural properties , 1991 .

[25]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[26]  P. Coveney,et al.  Continuum-particle hybrid coupling for mass, momentum, and energy transfers in unsteady fluid flow. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[28]  J M Thornton,et al.  Modelling protein unfolding: hen egg-white lysozyme. , 1997, Protein engineering.

[29]  S. Lowen.  The Biophysical Journal , 1960, Nature.