Exploring the Free Energy Landscape of Solutes Embedded in Lipid Bilayers.

Free energy calculations are vital for our understanding of biological processes on an atomistic scale and can offer insight to various mechanisms. However, in some cases, degrees of freedom (DOFs) orthogonal to the reaction coordinate have high energy barriers and/or long equilibration times, which prohibit proper sampling. Here we identify these orthogonal DOFs when studying the transfer of a solute from water to a model membrane. Important DOFs are identified in bulk liquids of different dielectric nature with metadynamics simulations and are used as reaction coordinates for the translocation process, resulting in two- and three-dimensional space of reaction coordinates. The results are in good agreement with experiments and elucidate the pitfalls of using one-dimensional reaction coordinates. The calculations performed here offer the most detailed free energy landscape of solutes embedded in lipid bilayers to date and show that free energy calculations can be used to study complex membrane translocation phenomena.

[1]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[2]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

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

[4]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[5]  D. Volk,et al.  Insight into NSAID-induced membrane alterations, pathogenesis and therapeutics: characterization of interaction of NSAIDs with phosphatidylcholine. , 2012, Biochimica et biophysica acta.

[6]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[7]  Karel Berka,et al.  Membrane Position of Ibuprofen Agrees with Suggested Access Path Entrance to Cytochrome P450 2C9 Active Site , 2011, The journal of physical chemistry. A.

[8]  David L Mobley,et al.  Perspective: Alchemical free energy calculations for drug discovery. , 2012, The Journal of chemical physics.

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

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

[11]  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.

[12]  A. Laio,et al.  A novel approach to the investigation of passive molecular permeation through lipid bilayers from atomistic simulations. , 2012, The journal of physical chemistry. B.

[13]  A. Avdeef,et al.  pH-Metric logP 10. Determination of Liposomal Membrane-Water Partition Coefficients of lonizable Drugs , 1998, Pharmaceutical Research.

[14]  A. Laio,et al.  Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamics. , 2006, The journal of physical chemistry. B.

[15]  P.-L. Chau.  New insights into the molecular mechanisms of general anaesthetics , 2010, British journal of pharmacology.

[16]  M. Parrinello,et al.  Sampling protein motion and solvent effect during ligand binding , 2012, Proceedings of the National Academy of Sciences.

[17]  H. Scheidt,et al.  Interaction of local anesthetics with lipid bilayers investigated by ¹H MAS NMR spectroscopy. , 2012, Biochimica et biophysica acta.

[18]  David L. Mobley,et al.  Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations , 2010, J. Comput. Aided Mol. Des..

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

[20]  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.

[21]  David L Mobley,et al.  Small Molecule Solvation Free Energy: Enhanced Conformational Sampling Using Expanded Ensemble Molecular Dynamics Simulation. , 2011, Journal of chemical theory and computation.

[22]  David L Mobley,et al.  The Confine-and-Release Method: Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change. , 2007, Journal of chemical theory and computation.

[23]  Benoît Roux,et al.  Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer. , 2012, Journal of chemical theory and computation.

[24]  D. Tieleman,et al.  Atomistic Simulations of Wimley-White Pentapeptides: Sampling of Structure and Dynamics in Solution. , 2013, Journal of chemical theory and computation.

[25]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[26]  Eric Darve,et al.  Assessing the efficiency of free energy calculation methods. , 2004, The Journal of chemical physics.

[27]  A. Laio,et al.  Assessing the accuracy of metadynamics. , 2005, The journal of physical chemistry. B.

[28]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[29]  B. Müller,et al.  Drug Liposome Partitioning as a Tool for the Prediction of Human Passive Intestinal Absorption , 1999, Pharmaceutical Research.

[30]  William Sinko,et al.  w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations , 2013, Journal of chemical theory and computation.

[31]  J. Baber,et al.  Distribution of general anesthetics in phospholipid bilayers determined using 2H NMR and 1H-1H NOE spectroscopy. , 1995, Biochemistry.

[32]  Christophe Chipot,et al.  Standard binding free energies from computer simulations: What is the best strategy? , 2013, Journal of chemical theory and computation.

[33]  David L Mobley,et al.  Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. , 2009, Journal of chemical theory and computation.

[34]  M. Klein,et al.  Distribution of halothane in a dipalmitoylphosphatidylcholine bilayer from molecular dynamics calculations. , 2000, Biophysical journal.

[35]  M. Parrinello,et al.  Targeting biomolecular flexibility with metadynamics. , 2010, Current opinion in structural biology.

[36]  Angel E Garcia,et al.  Molecular dynamics simulations suggest a mechanism for translocation of the HIV-1 TAT peptide across lipid membranes , 2007, Proceedings of the National Academy of Sciences.

[37]  Chau Pl,et al.  New insights into the molecular mechanisms of general anaesthetics , 2010 .

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

[39]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[40]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

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

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

[43]  Robert V. Swift,et al.  Back to the Future: Can Physical Models of Passive Membrane Permeability Help Reduce Drug Candidate Attrition and Move Us Beyond QSPR? , 2013, Chemical biology & drug design.

[44]  Stephen H. White,et al.  Experimentally determined hydrophobicity scale for proteins at membrane interfaces , 1996, Nature Structural Biology.

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

[46]  A. Warshel.  Calculations of enzymatic reactions: calculations of pKa, proton transfer reactions, and general acid catalysis reactions in enzymes. , 1981, Biochemistry.

[47]  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.

[48]  R. Krishnamoorti,et al.  Structural association of nonsteroidal anti-inflammatory drugs with lipid membranes. , 2012, Journal of the American Chemical Society.

[49]  Wonpil Im,et al.  Two Dimensional Window Exchange Umbrella Sampling for Transmembrane Helix Assembly. , 2013, Journal of chemical theory and computation.