Lipid Models for United-Atom Molecular Dynamics Simulations of Proteins.

United-atom force fields for molecular dynamics (MD) simulations provide a higher computational efficiency, especially in lipid membrane simulations, with little sacrifice in accuracy, when compared to all-atom force fields. Excellent united-atom lipid models are available, but in combination with depreciated protein force fields. In this work, a united-atom model of the lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine has been built with standard parameters of the force field GROMOS96 53a6 that reproduces the experimental area per lipid of a lipid bilayer within 3% accuracy to a value of 0.623 ± 0.011 nm(2) without the assumption of a constant surface area or the inclusion of surface pressure. In addition, the lateral self-diffusion constant and deuterium order parameters of the acyl chains are in agreement with experimental data. Furthermore, models for 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) result in areas per lipid of 0.625 nm(2) (DMPC), 0.693 nm(2) (POPC), and 0.700 nm(2) (POPG) from 40 ns MD simulations. Experimental lateral self-diffusion coefficients are reproduced satisfactorily by the simulation. The lipid models can form the basis for molecular dynamics simulations of membrane proteins with current and future versions of united-atom protein force fields.

[1]  J. Nagle,et al.  Closer look at structure of fully hydrated fluid phase DPPC bilayers. , 2006, Biophysical journal.

[2]  A. Kukol,et al.  Systematic molecular dynamics searching in a lipid bilayer: application to the glycophorin A and oncogenic ErbB-2 transmembrane domains. , 2006, Journal of molecular graphics & modelling.

[3]  J. Nagle,et al.  Lipid bilayer structure. , 2000, Current opinion in structural biology.

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

[5]  R A Sayle,et al.  RASMOL: biomolecular graphics for all. , 1995, Trends in biochemical sciences.

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

[7]  I. Vattulainen,et al.  Influence of cis double-bond parametrization on lipid membrane properties: how seemingly insignificant details in force-field change even qualitative trends. , 2008, The Journal of chemical physics.

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

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

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

[11]  G. Terstappen,et al.  In silico research in drug discovery. , 2001, Trends in pharmacological sciences.

[12]  M. Sansom,et al.  Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer. , 2006, Biochemistry.

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

[14]  Christian Kandt,et al.  Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid–protein interactions, side chain transfer free energies and model proteins , 2006, Journal of physics. Condensed matter : an Institute of Physics journal.

[15]  V. Parsegian,et al.  Measurement of repulsive forces between charged phospholipid bilayers. , 1978, Biochemistry.

[16]  Xicheng Wang,et al.  Mechanics of Channel Gating of the Nicotinic Acetylcholine Receptor , 2008, PLoS Comput. Biol..

[17]  H. Berendsen,et al.  A molecular dynamics study of the pores formed by Escherichia coli OmpF porin in a fully hydrated palmitoyloleoylphosphatidylcholine bilayer. , 1998, Biophysical journal.

[18]  P. Biggin,et al.  Molecular dynamics simulations of membrane proteins. , 2008, Methods in molecular biology.

[19]  H. Bernstein.  Recent changes to RasMol, recombining the variants. , 2000, Trends in biochemical sciences.

[20]  M. Patra,et al.  Molecular dynamics simulations of lipid bilayers: major artifacts due to truncating electrostatic interactions. , 2003, Biophysical journal.

[21]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[22]  E. Jakobsson,et al.  Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane. , 1995, Biophysical journal.

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

[24]  H. Berendsen,et al.  Molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine bilayer with different macroscopic boundary conditions and parameters , 1996 .

[25]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[26]  A. Kukol.  Site‒specific IR spectroscopy and molecular modelling combined towards solving transmembrane protein structure , 2005 .

[27]  Roman G. Efremov,et al.  Spatial Structure of the Dimeric Transmembrane Domain of the Growth Factor Receptor ErbB2 Presumably Corresponding to the Receptor Active State* , 2008, Journal of Biological Chemistry.

[28]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[29]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

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

[32]  I. Vattulainen,et al.  Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions. , 2007, Biophysical journal.

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

[34]  T. Róg,et al.  Phosphatidylethanolamine-phosphatidylglycerol bilayer as a model of the inner bacterial membrane. , 2005, Biophysical journal.

[35]  Roland Faller,et al.  Under the influence of alcohol: the effect of ethanol and methanol on lipid bilayers. , 2006, Biophysical journal.

[36]  M. Klein,et al.  United-atom acyl chains for CHARMM phospholipids. , 2008, The journal of physical chemistry. B.

[37]  Alexander D. MacKerell,et al.  Comparison of protein force fields for molecular dynamics simulations. , 2008, Methods in molecular biology.

[38]  D. Harries,et al.  Structure and fluctuations of charged phosphatidylserine bilayers in the absence of salt. , 2004, Biophysical journal.

[39]  Robert Vácha,et al.  Biomolecular simulations of membranes: physical properties from different force fields. , 2008, The Journal of chemical physics.

[40]  D. Peter Tieleman,et al.  A consistent potential energy parameter set for lipids: dipalmitoylphosphatidylcholine as a benchmark of the GROMOS96 45A3 force field , 2003, European Biophysics Journal.

[41]  P. Hünenberger,et al.  Molecular dynamics simulation of lipid bilayers with GROMOS96: Application of surface tension , 2005 .

[42]  M S Sansom,et al.  An alamethicin channel in a lipid bilayer: molecular dynamics simulations. , 1999, Biophysical journal.

[43]  L. Tamm,et al.  Annexin IV reduces the rate of lateral lipid diffusion and changes the fluid phase structure of the lipid bilayer when it binds to negatively charged membranes in the presence of calcium. , 1994, Biochemistry.

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

[45]  F. Ashcroft,et al.  Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study. , 2007, Biochemistry.

[46]  D. Tieleman,et al.  Molecular Dynamics Simulation of a Polyunsaturated Lipid Bilayer Susceptible to Lipid Peroxidation , 2004 .

[47]  O. Edholm,et al.  The range and shielding of dipole-dipole interactions in phospholipid bilayers. , 2004, Biophysical journal.

[48]  G. Heijne,et al.  Genome‐wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms , 1998, Protein science : a publication of the Protein Society.

[49]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[50]  G. Lindblom,et al.  Influence of Cholesterol and Water Content on Phospholipid Lateral Diffusion in Bilayers , 2003 .

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

[52]  Alexander D. MacKerell,et al.  An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids , 2000 .

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

[54]  S. May,et al.  Modelling of proteins in membranes. , 2006, Chemistry and physics of lipids.

[55]  D. Tieleman,et al.  Molecular dynamics simulation of a palmitoyl-oleoyl phosphatidylserine bilayer with Na+ counterions and NaCl. , 2004, Biophysical journal.

[56]  B. Tidor.  Molecular dynamics simulations , 1997, Current Biology.

[57]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[58]  Ian R. Gould,et al.  Structure and dynamics of phospholipid bilayers using recently developed general all‐atom force fields , 2008, J. Comput. Chem..

[59]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[60]  Wilfred F. van Gunsteren,et al.  Parametrization of aliphatic CHn united atoms of GROMOS96 force field , 1998, J. Comput. Chem..

[61]  D. Elmore.  Molecular dynamics simulation of a phosphatidylglycerol membrane , 2006, FEBS letters.

[62]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules J. Am. Chem. Soc. 1995, 117, 5179−5197 , 1996 .

[63]  T. Bayerl,et al.  Modulation of cytochrome C coupling to anionic lipid monolayers by a change of the phase state: a combined neutron and infrared reflection study. , 2000, Biophysical journal.

[64]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .