Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids.

Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.

[1]  B. Roux,et al.  Representation of Ion–Protein Interactions Using the Drude Polarizable Force-Field , 2015, The journal of physical chemistry. B.

[2]  Erik Lindahl,et al.  Molecular dynamics simulation of NMR relaxation rates and slow dynamics in lipid bilayers , 2001 .

[3]  Alexander D. MacKerell,et al.  Six-site polarizable model of water based on the classical Drude oscillator. , 2013, The Journal of chemical physics.

[4]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[5]  Benoît Roux,et al.  AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA. , 2013, Journal of chemical theory and computation.

[6]  Peter Pulay,et al.  Systematic AB Initio Gradient Calculation of Molecular Geometries, Force Constants, and Dipole Moment Derivatives , 1979 .

[7]  E. Dufourc,et al.  Restatement of order parameters in biomembranes: calculation of C-C bond order parameters from C-D quadrupolar splittings. , 1995, Biophysical journal.

[8]  Benoît Roux,et al.  A polarizable force field of dipalmitoylphosphatidylcholine based on the classical Drude model for molecular dynamics simulations of lipids. , 2013, The journal of physical chemistry. B.

[9]  Klaus Schulten,et al.  High-performance scalable molecular dynamics simulations of a polarizable force field based on classical Drude oscillators in NAMD. , 2011, The journal of physical chemistry letters.

[10]  Alexander D. MacKerell,et al.  Force Field for Peptides and Proteins based on the Classical Drude Oscillator. , 2013, Journal of chemical theory and computation.

[11]  Alexander D. MacKerell,et al.  An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications , 2016, Chemical reviews.

[12]  Siewert J. Marrink,et al.  Methodological issues in lipid bilayer simulations , 2003 .

[13]  J. Dilger,et al.  The dielectric constant of phospholipid bilayers and the permeability of membranes to ions. , 1979, Science.

[14]  Leo Radom,et al.  Harmonic Vibrational Frequencies: An Evaluation of Hartree−Fock, Møller−Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors , 1996 .

[15]  Alexander D. MacKerell,et al.  A simple polarizable model of water based on classical Drude oscillators , 2003 .

[16]  V. Parsegian,et al.  Hydration forces between phospholipid bilayers , 1989 .

[17]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[18]  A. D. McLean,et al.  Contracted Gaussian basis sets for molecular calculations. I. Second row atoms, Z=11–18 , 1980 .

[19]  J. Seelig,et al.  Bilayers of dipalmitoyl-3-sn-phosphatidylcholine. Conformational differences between the fatty acyl chains. , 1975, Biochimica et biophysica acta.

[20]  Alexander D. MacKerell,et al.  Recent Advances in Polarizable Force Fields for Macromolecules: Microsecond Simulations of Proteins Using the Classical Drude Oscillator Model , 2014, The journal of physical chemistry letters.

[21]  Alexander D. MacKerell,et al.  A polarizable model of water for molecular dynamics simulations of biomolecules , 2006 .

[22]  John Katsaras,et al.  Molecular structures of fluid phosphatidylethanolamine bilayers obtained from simulation-to-experiment comparisons and experimental scattering density profiles. , 2015, The journal of physical chemistry. B.

[23]  Jeremy Pencer,et al.  Lipid bilayer structure determined by the simultaneous analysis of neutron and X-ray scattering data. , 2008, Biophysical journal.

[24]  Alexander D. MacKerell,et al.  Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator. , 2005, Journal of chemical theory and computation.

[25]  Alexander D. MacKerell,et al.  CHARMM additive and polarizable force fields for biophysics and computer-aided drug design. , 2015, Biochimica et biophysica acta.

[26]  J. Sheats,et al.  A photochemical technique for measuring lateral diffusion of spin-labeled phospholipids in membranes. , 1978, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Alexander D. MacKerell,et al.  Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics. , 2017, Journal of chemical theory and computation.

[28]  J. Nagle,et al.  Temperature dependence of structure, bending rigidity, and bilayer interactions of dioleoylphosphatidylcholine bilayers. , 2008, Biophysical journal.

[29]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[30]  P. Westerman,et al.  A model of orientational ordering in phosphatidylcholine bilayers based on conformational analysis of the glycerol backbone region. , 1985, Biophysical journal.

[31]  Alexander D. MacKerell,et al.  Polyunsaturated fatty acids in lipid bilayers: intrinsic and environmental contributions to their unique physical properties. , 2002, Journal of the American Chemical Society.

[32]  G. Gomila,et al.  Nanoscale measurement of the dielectric constant of supported lipid bilayers in aqueous solutions with electrostatic force microscopy. , 2013, Biophysical journal.

[33]  G. Stark.  Negative hydrophobic ions as transport-mediators for positive ions: evidence for a carrier mechanism. , 1980, Biochimica et biophysica acta.

[34]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[35]  Thomas B Woolf,et al.  Atomistic simulations of biologically realistic transmembrane potential gradients. , 2004, The Journal of chemical physics.

[36]  Alexander D. MacKerell,et al.  Chapter 1 Considerations for Lipid Force Field Development , 2008 .

[37]  J. Seelig,et al.  Conformation and motion of the choline head group in bilayers of dipalmitoyl-3-sn-phosphatidylcholine. , 1975, Biochemistry.

[38]  Taehoon Kim,et al.  CHARMM‐GUI: A web‐based graphical user interface for CHARMM , 2008, J. Comput. Chem..

[39]  Alan E Mark,et al.  Lipid Bilayers: The Effect of Force Field on Ordering and Dynamics. , 2012, Journal of chemical theory and computation.

[40]  J. Nagle,et al.  Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using X-ray scattering from oriented multilamellar arrays and from unilamellar vesicles. , 2005, Biophysical journal.

[41]  H. Scheidt,et al.  Diffusion of cholesterol and its precursors in lipid membranes studied by 1H pulsed field gradient magic angle spinning NMR. , 2005, Biophysical journal.

[42]  Structure and polymorphism of the hydrocarbon chains of lipids: a study of lecithin-water phases. , 1973 .

[43]  P. Devaux,et al.  Order parameters of unsaturated phospholipids in membranes and the effect of cholesterol: a 1H–13C solid-state NMR study at natural abundance , 2005, European Biophysics Journal.

[44]  T. Unruh,et al.  Molecular mechanism of long-range diffusion in phospholipid membranes studied by quasielastic neutron scattering. , 2010, Journal of the American Chemical Society.

[45]  Jeffery B. Klauda,et al.  Simulation-based methods for interpreting x-ray data from lipid bilayers. , 2006, Biophysical journal.

[46]  O. H. S. Ollila,et al.  Atomistic resolution structure and dynamics of lipid bilayers in simulations and experiments. , 2016, Biochimica et biophysica acta.

[47]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[48]  V A Parsegian,et al.  Membrane dipole potentials, hydration forces, and the ordering of water at membrane surfaces. , 1992, Biophysical journal.

[49]  H. A. Lorentz.  Ueber die Anwendung des Satzes vom Virial in der kinetischen Theorie der Gase , 1881 .

[50]  John F. Nagle,et al.  Structure of Fully Hydrated Fluid Phase Lipid Bilayers with Monounsaturated Chains , 2006, The Journal of Membrane Biology.

[51]  Arieh Warshel,et al.  Modeling electrostatic effects in proteins. , 2006, Biochimica et biophysica acta.

[52]  Alexander D. MacKerell,et al.  Polarizable empirical force field for alkanes based on the classical Drude oscillator model. , 2005, The journal of physical chemistry. B.

[53]  B. Roux,et al.  Simulation of Osmotic Pressure in Concentrated Aqueous Salt Solutions , 2010 .

[54]  J. Tabony,et al.  Quasielastic neutron scattering measurements of fast local translational diffusion of lipid molecules in phospholipid bilayers. , 1991, Biochimica et biophysica acta.

[55]  Alexander D. MacKerell,et al.  Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field. , 2010, Journal of chemical theory and computation.

[56]  Alexander D. MacKerell,et al.  Understanding the dielectric properties of liquid amides from a polarizable force field. , 2008, The journal of physical chemistry. B.

[57]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[58]  F. Allen.  The Cambridge Structural Database: a quarter of a million crystal structures and rising. , 2002, Acta crystallographica. Section B, Structural science.

[59]  Björn O. Roos,et al.  Second-order perturbation theory with a complete active space self-consistent field reference function , 1992 .

[60]  Benoît Roux,et al.  Modeling induced polarization with classical Drude oscillators: Theory and molecular dynamics simulation algorithm , 2003 .

[61]  R. Thurmond,et al.  Molecular areas of phospholipids as determined by 2H NMR spectroscopy. Comparison of phosphatidylethanolamines and phosphatidylcholines. , 1991, Biophysical journal.

[62]  Alexander D. MacKerell,et al.  Implementation of extended Lagrangian dynamics in GROMACS for polarizable simulations using the classical Drude oscillator model , 2015, J. Comput. Chem..

[63]  Ronald J Clarke,et al.  Hydrophobic ion hydration and the magnitude of the dipole potential. , 2002, Biophysical journal.

[64]  R. Clegg,et al.  Translational diffusion of lipids in liquid crystalline phase phosphatidylcholine multibilayers. A comparison of experiment with theory. , 1985, Biochemistry.

[65]  G. Shipley,et al.  Nature of the Thermal pretransition of synthetic phospholipids: dimyristolyl- and dipalmitoyllecithin. , 1976, Biochemistry.

[66]  Frank Bringezu,et al.  Mechanism of the lamellar/inverse hexagonal phase transition examined by high resolution x-ray diffraction. , 2003, Biophysical journal.

[67]  K. Gawrisch,et al.  Effect of Unsaturated Lipid Chains on Dimensions, Molecular Order and Hydration of Membranes , 2001 .

[68]  G. Bryant,et al.  Interactions between soluble sugars and POPC (1-palmitoyl-2-oleoylphosphatidylcholine) during dehydration: vitrification of sugars alters the phase behavior of the phospholipid. , 1994, Biochimica et biophysica acta.

[69]  D. Haydon,et al.  Ion transfer across lipid membranes in the presence of gramicidin A. II. The ion selectivity. , 1972, Biochimica et biophysica acta.

[70]  R. Pastor,et al.  Mechanical properties of lipid bilayers from molecular dynamics simulation. , 2015, Chemistry and physics of lipids.

[71]  J. Seelig,et al.  The dynamic structure of fatty acyl chains in a phospholipid bilayer measured by deuterium magnetic resonance. , 1974, Biochemistry.

[72]  Mu-Ping Nieh,et al.  Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. , 2011, Biochimica et biophysica acta.

[73]  Alexander D. MacKerell,et al.  All‐atom polarizable force field for DNA based on the classical drude oscillator model , 2014, J. Comput. Chem..

[74]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[75]  H. Jarrell,et al.  Acyl chain dynamics of phosphatidylethanolamines containing oleic acid and dihydrosterculic acid: 2H NMR relaxation studies. , 1985, Biochemistry.

[76]  Benoît Roux,et al.  On the statistical equivalence of restrained-ensemble simulations with the maximum entropy method. , 2013, The Journal of chemical physics.

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

[78]  Brad A. Bauer,et al.  Charge equilibration force fields for molecular dynamics simulations of lipids, bilayers, and integral membrane protein systems. , 2012, Biochimica et biophysica acta.

[79]  J. Seelig,et al.  Structure of Escherichia coli membranes. Glycerol auxotrophs as a tool for the analysis of the phospholipid head-group region by deuterium magentic resonance. , 1981, Biochemistry.

[80]  Benoît Roux,et al.  Restrained-ensemble molecular dynamics simulations based on distance histograms from double electron-electron resonance spectroscopy. , 2013, The journal of physical chemistry. B.

[81]  Alexander P. Lyubartsev,et al.  Derivation and Systematic Validation of a Refined All-Atom Force Field for Phosphatidylcholine Lipids , 2012, The journal of physical chemistry. B.

[82]  S. McLaughlin,et al.  The electrostatic properties of membranes. , 1989, Annual review of biophysics and biophysical chemistry.

[83]  J. Nagle,et al.  Structure of lipid bilayers. , 2000, Biochimica et biophysica acta.

[84]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[85]  Helgi I. Ingólfsson,et al.  Lipid and Peptide Diffusion in Bilayers: The Saffman-Delbrück Model and Periodic Boundary Conditions. , 2017, The journal of physical chemistry. B.

[86]  J. Nagle,et al.  Structure and interactions of fully hydrated dioleoylphosphatidylcholine bilayers. , 1998, Biophysical journal.

[87]  Jejoong Yoo,et al.  Improved Parametrization of Li+, Na+, K+, and Mg2+ Ions for All-Atom Molecular Dynamics Simulations of Nucleic Acid Systems , 2012 .

[88]  Alexander D. MacKerell,et al.  CHARMM Drude Polarizable Force Field for Aldopentofuranoses and Methyl-aldopentofuranosides. , 2015, The journal of physical chemistry. B.

[89]  Jeffery B. Klauda,et al.  A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties. , 2014, Biochimica et biophysica acta.

[90]  Alexander D. MacKerell,et al.  Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator , 2014, The journal of physical chemistry. B.

[91]  Alexander D. MacKerell,et al.  Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. , 2010, The journal of physical chemistry. B.

[92]  J. Seelig,et al.  Molecular order in cis and trans unsaturated phospholipid bilayers. , 1978, Biochemistry.

[93]  G Pabst,et al.  Structural information from multilamellar liposomes at full hydration: full q-range fitting with high quality x-ray data. , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[94]  R. Templer,et al.  Estimations of lipid bilayer geometry in fluid lamellar phases. , 2000, Biochimica et biophysica acta.

[95]  Margaret E. Johnson,et al.  Current status of the AMOEBA polarizable force field. , 2010, The journal of physical chemistry. B.

[96]  Alexander D. MacKerell,et al.  Development of the CHARMM Force Field for Lipids. , 2011, The journal of physical chemistry letters.

[97]  Alexander D. MacKerell,et al.  Ab Initio Calculations on the Use of Helium and Neon as Probes of the van der Waals Surfaces of Molecules , 1996 .

[98]  Alexander P Lyubartsev,et al.  An Extension and Further Validation of an All-Atomistic Force Field for Biological Membranes. , 2012, Journal of chemical theory and computation.

[99]  Alexander D. MacKerell,et al.  Development of CHARMM polarizable force field for nucleic acid bases based on the classical Drude oscillator model. , 2011, The journal of physical chemistry. B.

[100]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[101]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

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

[103]  C. Berde,et al.  Gel to liquid-crystalline phase transitions of aqueous dispersions of polyunsaturated mixed-acid phosphatidylcholines. , 1983, Biochemistry.

[104]  S. Dodd,et al.  Area per lipid and acyl length distributions in fluid phosphatidylcholines determined by (2)H NMR spectroscopy. , 2000, Biophysical journal.

[105]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[106]  Richard W. Pastor,et al.  Constant surface tension simulations of lipid bilayers: The sensitivity of surface areas and compressibilities , 1999 .

[107]  A. Petrov,et al.  Optical detection of phase transitions in simple and mixed lipid-water phases. , 1982, Biochimica et biophysica acta.

[108]  E. Evans,et al.  Effect of chain length and unsaturation on elasticity of lipid bilayers. , 2000, Biophysical journal.

[109]  Fred J Sigworth,et al.  Using cryo-EM to measure the dipole potential of a lipid membrane , 2006, Proceedings of the National Academy of Sciences.

[110]  Alexander D. MacKerell,et al.  Many-body polarization effects and the membrane dipole potential. , 2009, Journal of the American Chemical Society.

[111]  M. Sundaralihgam.  DISCUSSION PAPER: MOLECULAR STRUCTURES AND CONFORMATIONS OF THE PHOSPHOLIPIDS AND SPHINGOMYELINS * , 1972, Annals of the New York Academy of Sciences.

[112]  Benoît Roux,et al.  Structural refinement from restrained-ensemble simulations based on EPR/DEER data: application to T4 lysozyme. , 2013, The journal of physical chemistry. B.

[113]  Edward Teller,et al.  Interaction of the van der Waals Type Between Three Atoms , 1943 .