Accelerating Convergence in Molecular Dynamics Simulations of Solutes in Lipid Membranes by Conducting a Random Walk along the Bilayer Normal.

All molecular dynamics simulations are susceptible to sampling errors, which degrade the accuracy and precision of observed values. The statistical convergence of simulations containing atomistic lipid bilayers is limited by the slow relaxation of the lipid phase, which can exceed hundreds of nanoseconds. These long conformational autocorrelation times are exacerbated in the presence of charged solutes, which can induce significant distortions of the bilayer structure. Such long relaxation times represent hidden barriers that induce systematic sampling errors in simulations of solute insertion. To identify optimal methods for enhancing sampling efficiency, we quantitatively evaluate convergence rates using generalized ensemble sampling algorithms in calculations of the potential of mean force for the insertion of the ionic side chain analog of arginine in a lipid bilayer. Umbrella sampling (US) is used to restrain solute insertion depth along the bilayer normal, the order parameter commonly used in simulations of molecular solutes in lipid bilayers. When US simulations are modified to conduct random walks along the bilayer normal using a Hamiltonian exchange algorithm, systematic sampling errors are eliminated more rapidly and the rate of statistical convergence of the standard free energy of binding of the solute to the lipid bilayer is increased 3-fold. We compute the ratio of the replica flux transmitted across a defined region of the order parameter to the replica flux that entered that region in Hamiltonian exchange simulations. We show that this quantity, the transmission factor, identifies sampling barriers in degrees of freedom orthogonal to the order parameter. The transmission factor is used to estimate the depth-dependent conformational autocorrelation times of the simulation system, some of which exceed the simulation time, and thereby identify solute insertion depths that are prone to systematic sampling errors and estimate the lower bound of the amount of sampling that is required to resolve these sampling errors. Finally, we extend our simulations and verify that the conformational autocorrelation times estimated by the transmission factor accurately predict correlation times that exceed the simulation time scale-something that, to our knowledge, has never before been achieved.

[1]  Samira Hezaveh,et al.  Diffusion of 1,2-dimethoxyethane and 1,2-dimethoxypropane through phosphatidycholine bilayers: a molecular dynamics study. , 2012, The journal of physical chemistry. B.

[2]  Toby W Allen,et al.  Assessing atomistic and coarse-grained force fields for protein-lipid interactions: the formidable challenge of an ionizable side chain in a membrane. , 2008, The journal of physical chemistry. B.

[3]  Y. Kaznessis,et al.  Free energy profile of the interaction between a monomer or a dimer of protegrin-1 in a specific binding orientation and a model lipid bilayer. , 2010, The journal of physical chemistry. B.

[4]  Daniel M Zuckerman,et al.  Equilibrium sampling in biomolecular simulations. , 2011, Annual review of biophysics.

[5]  Berk Hess,et al.  P-LINCS:  A Parallel Linear Constraint Solver for Molecular Simulation. , 2008, Journal of chemical theory and computation.

[6]  Erik Lindahl,et al.  Position‐resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations , 2007, Proteins.

[7]  W F Drew Bennett,et al.  Statistical Convergence of Equilibrium Properties in Simulations of Molecular Solutes Embedded in Lipid Bilayers. , 2011, Journal of chemical theory and computation.

[8]  K. Berka,et al.  Positioning of antioxidant quercetin and its metabolites in lipid bilayer membranes: implication for their lipid-peroxidation inhibition. , 2012, The journal of physical chemistry. B.

[9]  Justin L MacCallum,et al.  Computer simulation of the distribution of hexane in a lipid bilayer: spatially resolved free energy, entropy, and enthalpy profiles. , 2006, Journal of the American Chemical Society.

[10]  Herman J. C. Berendsen,et al.  Permeation Process of Small Molecules across Lipid Membranes Studied by Molecular Dynamics Simulations , 1996 .

[11]  O. Berger,et al.  Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. , 1997, Biophysical journal.

[12]  J. Essex,et al.  Behaviour of small solutes and large drugs in a lipid bilayer from computer simulations. , 2005, Biochimica et biophysica acta.

[13]  B. Roux Lonely Arginine Seeks Friendly Environment , 2007, The Journal of general physiology.

[14]  M. Klein,et al.  Effects of anesthetics on the structure of a phospholipid bilayer: molecular dynamics investigation of halothane in the hydrated liquid crystal phase of dipalmitoylphosphatidylcholine. , 1998, Biophysical journal.

[15]  Terry R. Stouch,et al.  Orientation and Diffusion of a Drug Analog in Biomembranes: Molecular Dynamics Simulations , 1995 .

[16]  N. Abbott,et al.  Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations. , 2005, Biophysical journal.

[17]  Brad A. Bauer,et al.  Free energetics and the role of water in the permeation of methyl guanidinium across the bilayer-water interface: insights from molecular dynamics simulations using charge equilibration potentials. , 2013, The journal of physical chemistry. B.

[18]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[19]  Daniel M Zuckerman,et al.  A Second Look at Canonical Sampling of Biomolecules Using Replica Exchange Simulation [J. Chem. Theory Comput. 2, 1200-1202 (2006)]. , 2006, Journal of chemical theory and computation.

[20]  D. Tieleman,et al.  Thermodynamic analysis of the effect of cholesterol on dipalmitoylphosphatidylcholine lipid membranes. , 2009, Journal of the American Chemical Society.

[21]  D. Tieleman,et al.  Thermodynamics of flip-flop and desorption for a systematic series of phosphatidylcholine lipids , 2009 .

[22]  A steroid in a lipid bilayer: localization, orientation, and energetics. , 2008, Biophysical journal.

[23]  Q. Cui,et al.  Does arginine remain protonated in the lipid membrane? Insights from microscopic pKa calculations. , 2008, Biophysical journal.

[24]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[25]  Interaction of K+ with a Phospholipid Bilayer: A Molecular Dynamics Study , 1997 .

[26]  Terry R. Stouch,et al.  Mechanism of Solute Diffusion through Lipid Bilayer Membranes by Molecular Dynamics Simulation , 1995 .

[27]  Shin-Ho Chung,et al.  Effect of Gating Modifier Toxins on Membrane Thickness: Implications for Toxin Effect on Gramicidin and Mechanosensitive Channels , 2013, Toxins.

[28]  Justin L. MacCallum,et al.  Partitioning of Amino Acid Side Chains into Lipid Bilayers: Results from Computer Simulations and Comparison to Experiment , 2007, The Journal of general physiology.

[29]  E. Pai,et al.  An iris-like mechanism of pore dilation in the CorA magnesium transport system. , 2010, Biophysical journal.

[30]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[31]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[32]  W F Drew Bennett,et al.  Molecular view of cholesterol flip-flop and chemical potential in different membrane environments. , 2009, Journal of the American Chemical Society.

[33]  Alexander D. MacKerell,et al.  Molecular-level organization of saturated and polyunsaturated fatty acids in a phosphatidylcholine bilayer containing cholesterol. , 2004, Biochemistry.

[34]  M. Berkowitz,et al.  Interaction between amyloid-beta (1-42) peptide and phospholipid bilayers: a molecular dynamics study. , 2009, Biophysical journal.

[35]  Toby W Allen,et al.  On the thermodynamic stability of a charged arginine side chain in a transmembrane helix , 2007, Proceedings of the National Academy of Sciences.

[36]  Henning Henschel,et al.  How warfarin's structural diversity influences its phospholipid bilayer membrane permeation. , 2013, The journal of physical chemistry. B.

[37]  M. Berkowitz,et al.  Energetics of cholesterol transfer between lipid bilayers. , 2008, The journal of physical chemistry. B.

[38]  Ramanan Krishnamoorti,et al.  Partitioning of nonsteroidal antiinflammatory drugs in lipid membranes: a molecular dynamics simulation study. , 2010, Biophysical journal.

[39]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[40]  Peter G Bolhuis,et al.  Reordering hydrogen bonds using hamiltonian replica exchange enhances sampling of conformational changes in biomolecular systems. , 2009, The journal of physical chemistry. B.

[41]  Interaction of Tryptophan Analogs with POPC Lipid Bilayers Investigated by Molecular Dynamics Calculations , 2002 .

[42]  R. Pomès,et al.  Equilibrium exchange enhances the convergence rate of umbrella sampling , 2008 .

[43]  M. Berkowitz,et al.  Binding and reorientation of melittin in a POPC bilayer: computer simulations. , 2012, Biochimica et biophysica acta.

[44]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[45]  H. Vogel,et al.  Distribution of pentachlorophenol in phospholipid bilayers: a molecular dynamics study. , 2004, Biophysical journal.

[46]  Fancui Meng,et al.  Drug permeability prediction using PMF method , 2013, Journal of Molecular Modeling.

[47]  A. Mark,et al.  Convergence and sampling efficiency in replica exchange simulations of peptide folding in explicit solvent. , 2007, The Journal of chemical physics.

[48]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[49]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[50]  R. Faller,et al.  Density imbalances and free energy of lipid transfer in supported lipid bilayers. , 2009, The Journal of chemical physics.

[51]  M S Sansom,et al.  Alamethicin helices in a bilayer and in solution: molecular dynamics simulations. , 1999, Biophysical journal.

[52]  Bernard R Brooks,et al.  Simulations of membranes and other interfacial systems using P2(1) and Pc periodic boundary conditions. , 2002, Biophysical journal.

[53]  Molecular dynamics study of water pores in a phospholipid bilayer , 2002 .

[54]  Wonpil Im,et al.  Transmembrane helix assembly by window exchange umbrella sampling. , 2012, Physical review letters.

[55]  D. Kofke,et al.  Selection of temperature intervals for parallel-tempering simulations. , 2005, The Journal of chemical physics.

[56]  B. Roux The calculation of the potential of mean force using computer simulations , 1995 .

[57]  P. Biggin,et al.  Distribution and dynamics of adamantanes in a lipid bilayer. , 2008, Biophysical journal.

[58]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[59]  Wonpil Im,et al.  Revisiting hydrophobic mismatch with free energy simulation studies of transmembrane helix tilt and rotation. , 2010, Biophysical journal.

[60]  G. Voth,et al.  Protons may leak through pure lipid bilayers via a concerted mechanism. , 2005, Biophysical journal.

[61]  Toby W Allen,et al.  The role of tryptophan side chains in membrane protein anchoring and hydrophobic mismatch. , 2013, Biochimica et biophysica acta.

[62]  B Honig,et al.  Free-energy determinants of alpha-helix insertion into lipid bilayers. , 1996, Biophysical journal.

[63]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[64]  Alan Grossfield,et al.  Convergence of molecular dynamics simulations of membrane proteins , 2007, Proteins.

[65]  J. Mccammon,et al.  Free Energy for the Permeation of Na+ and Cl− Ions and Their Ion-Pair through a Zwitterionic Dimyristoyl Phosphatidylcholine Lipid Bilayer by Umbrella Integration with Harmonic Fourier Beads , 2009, Journal of the American Chemical Society.

[66]  D Peter Tieleman,et al.  Lipids out of equilibrium: energetics of desorption and pore mediated flip-flop. , 2006, Journal of the American Chemical Society.

[67]  Benoît Roux,et al.  Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations. , 2010, Journal of chemical theory and computation.

[68]  W. Shinoda,et al.  Molecular Dynamics Study on the Effects of Chain Branching on the Physical Properties of Lipid Bilayers: 2. Permeability , 2004 .

[69]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[70]  Ching-Hsing Yu,et al.  SciNet: Lessons Learned from Building a Power-efficient Top-20 System and Data Centre , 2010 .

[71]  Christian Kandt,et al.  Setting up and running molecular dynamics simulations of membrane proteins. , 2007, Methods.

[72]  I. Andricioaei,et al.  Poly(amidoamine) dendrimers on lipid bilayers I: Free energy and conformation of binding. , 2008, The journal of physical chemistry. B.

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

[74]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[75]  K. Hristova,et al.  A Look at Arginine in Membranes , 2010, The Journal of Membrane Biology.

[76]  Ignacio J General A Note on the Standard State's Binding Free Energy. , 2010, Journal of chemical theory and computation.

[77]  Jon Beckwith,et al.  The role of charged amino acids in the localization of secreted and membrane proteins , 1990, Cell.

[78]  A. Pohorille,et al.  Mechanism of unassisted ion transport across membrane bilayers. , 1996, Journal of the American Chemical Society.

[79]  Régis Pomès,et al.  Calculation of absolute protein-ligand binding free energy using distributed replica sampling. , 2008, The Journal of chemical physics.

[80]  Anders Wallqvist,et al.  Free energy difference in indolicidin attraction to eukaryotic and prokaryotic model cell membranes. , 2012, The journal of physical chemistry. B.

[81]  Y. Okamoto Institute for Molecular Science,et al.  Replica-exchange simulated tempering method for simulations of frustrated systems , 2000 .

[82]  A. Pohorille,et al.  Interactions of anesthetics with the membrane-water interface. , 1996, Chemical physics.

[83]  Terry R Stouch,et al.  Functional group dependence of solute partitioning to various locations within a DOPC bilayer: a comparison of molecular dynamics simulations with experiment. , 2011, Journal of pharmaceutical sciences.

[84]  Chris Neale,et al.  Simulated Tempering Distributed Replica Sampling, Virtual Replica Exchange, and Other Generalized-Ensemble Methods for Conformational Sampling. , 2009, Journal of chemical theory and computation.

[85]  S. Whittington,et al.  Monte carlo study of the interacting self-avoiding walk model in three dimensions , 1996 .

[86]  Wei Yang,et al.  Simultaneous escaping of explicit and hidden free energy barriers: application of the orthogonal space random walk strategy in generalized ensemble based conformational sampling. , 2009, The Journal of chemical physics.

[87]  I. Vattulainen,et al.  Mechanism for translocation of fluoroquinolones across lipid membranes. , 2012, Biochimica et biophysica acta.

[88]  Wilfred F van Gunsteren,et al.  Molecular simulation as an aid to experimentalists. , 2008, Current opinion in structural biology.

[89]  Edward Lyman,et al.  A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation. , 2006, Journal of chemical theory and computation.

[90]  Arun Yethiraj,et al.  Sequence-dependent interaction of β-peptides with membranes. , 2010, The journal of physical chemistry. B.

[91]  G. Voth,et al.  Mechanisms of passive ion permeation through lipid bilayers: insights from simulations. , 2006, The journal of physical chemistry. B.

[92]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[93]  Sunhwan Jo,et al.  Cholesterol flip-flop: insights from free energy simulation studies. , 2010, The journal of physical chemistry. B.

[94]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[95]  S Gnanakaran,et al.  Taste of sugar at the membrane: thermodynamics and kinetics of the interaction of a disaccharide with lipid bilayers. , 2013, Biophysical journal.

[96]  Régis Pomès,et al.  Distributed Replica Sampling. , 2006, Journal of chemical theory and computation.

[97]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[98]  Replica-exchange molecular dynamics simulation of small peptide in water and in ethanol , 2005 .

[99]  F. Rao,et al.  Replica exchange molecular dynamics simulations of amyloid peptide aggregation. , 2004, The Journal of chemical physics.

[100]  K. Berka,et al.  Molecular insight into affinities of drugs and their metabolites to lipid bilayers. , 2013, The journal of physical chemistry. B.

[101]  M. Troyer,et al.  Optimized parallel tempering simulations of proteins. , 2006, The Journal of chemical physics.

[102]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[103]  Mihaly Mezei,et al.  Studies on free energy calculations. I. Thermodynamic integration using a polynomial path , 1993 .

[104]  Matthew A. Wyczalkowski,et al.  Satisfying the fluctuation theorem in free-energy calculations with Hamiltonian replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[105]  H. Nymeyer,et al.  Indole localization in lipid membranes revealed by molecular simulation. , 2006, Biophysical journal.

[106]  J. Beckwith,et al.  Positively charged amino acid residues can act as topogenic determinants in membrane proteins. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[107]  Is arginine charged in a membrane? , 2008, Biophysical journal.

[108]  H. Nymeyer,et al.  Folding is not required for bilayer insertion: Replica exchange simulations of an α‐helical peptide with an explicit lipid bilayer , 2004, Proteins.

[109]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[110]  Angel E García,et al.  Free energy of translocating an arginine-rich cell-penetrating peptide across a lipid bilayer suggests pore formation. , 2013, Biophysical journal.

[111]  Peter A. Kollman,et al.  The lag between the Hamiltonian and the system configuration in free energy perturbation calculations , 1989 .

[112]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[113]  K. Berka,et al.  Lipid bilayer membrane affinity rationalizes inhibition of lipid peroxidation by a natural lignan antioxidant. , 2013, The journal of physical chemistry. B.

[114]  Peter Westh,et al.  Binding of serotonin to lipid membranes. , 2013, Journal of the American Chemical Society.

[115]  Herman J. C. Berendsen,et al.  Simulation of Water Transport through a Lipid Membrane , 1994 .

[116]  D. van der Spoel,et al.  Large influence of cholesterol on solute partitioning into lipid membranes. , 2012, Journal of the American Chemical Society.

[117]  Alexander P Lyubartsev,et al.  Implicit inclusion of atomic polarization in modeling of partitioning between water and lipid bilayers. , 2013, Physical chemistry chemical physics : PCCP.

[118]  R. Elber,et al.  Unassisted transport of N-acetyl-L-tryptophanamide through membrane: experiment and simulation of kinetics. , 2012, The journal of physical chemistry. B.

[119]  Tod D Romo,et al.  Block Covariance Overlap Method and Convergence in Molecular Dynamics Simulation. , 2011, Journal of chemical theory and computation.

[120]  Ulrich H E Hansmann,et al.  Folding proteins by first-passage-times-optimized replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[121]  Danial S. Dashti,et al.  Computing Alchemical Free Energy Differences with Hamiltonian Replica Exchange Molecular Dynamics (H-REMD) Simulations. , 2011, Journal of chemical theory and computation.

[122]  J. Ulander,et al.  Permeation across hydrated DPPC lipid bilayers: simulation of the titrable amphiphilic drug valproic acid. , 2003, Biophysical journal.

[123]  R. Dalbey Positively charged residues are important determinants of membrane protein topology. , 1990, Trends in biochemical sciences.

[124]  S. D. de Leeuw,et al.  Rapid free energy calculation of peptide self-assembly by REMD umbrella sampling. , 2008, The journal of physical chemistry. B.

[125]  A. Violi,et al.  Molecular dynamics simulation study of a pulmonary surfactant film interacting with a carbonaceous nanoparticle. , 2008, Biophysical journal.

[126]  Toby W Allen,et al.  Potential of mean force and pKa profile calculation for a lipid membrane-exposed arginine side chain. , 2008, The journal of physical chemistry. B.

[127]  A. Gorfe,et al.  Thermodynamics of Peptide Insertion and Aggregation in a Lipid Bilayer , 2008, The journal of physical chemistry. B.

[128]  D. Tieleman,et al.  Molecular simulation of rapid translocation of cholesterol, diacylglycerol, and ceramide in model raft and nonraft membranes[S] , 2012, Journal of Lipid Research.

[129]  W F Drew Bennett,et al.  Distribution of amino acids in a lipid bilayer from computer simulations. , 2008, Biophysical journal.

[130]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[131]  S. Marrink,et al.  Location, tilt, and binding: a molecular dynamics study of voltage-sensitive dyes in biomembranes. , 2009, The journal of physical chemistry. B.

[132]  Robert I Cukier,et al.  Molecular dynamics of apo-adenylate kinase: a distance replica exchange method for the free energy of conformational fluctuations. , 2006, The journal of physical chemistry. B.

[133]  F. Rao,et al.  Replica exchange molecular dynamics simulations of reversible folding , 2003 .

[134]  Karel Berka,et al.  Convergence of Free Energy Profile of Coumarin in Lipid Bilayer , 2012, Journal of chemical theory and computation.

[135]  W F Drew Bennett,et al.  The Role of Atomic Polarization in the Thermodynamics of Chloroform Partitioning to Lipid Bilayers. , 2012, Journal of chemical theory and computation.

[136]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[137]  Jonathan W. Essex,et al.  The development of replica-exchange-based free-energy methods , 2003 .

[138]  Jacek Waluk,et al.  Partitioning and localization of environment-sensitive 2-(2'-pyridyl)- and 2-(2'-pyrimidyl)-indoles in lipid membranes: a joint refinement using fluorescence measurements and molecular dynamics simulations. , 2010, The journal of physical chemistry. B.