Structural, dynamic, and electrostatic properties of fully hydrated DMPC bilayers from molecular dynamics simulations accelerated with graphical processing units (GPUs)

We present results of molecular dynamics simulations of fully hydrated DMPC bilayers performed on graphics processing units (GPUs) using current state‐of‐the‐art non‐polarizable force fields and a local GPU‐enabled molecular dynamics code named FEN ZI. We treat the conditionally convergent electrostatic interaction energy exactly using the particle mesh Ewald method (PME) for solution of Poisson's Equation for the electrostatic potential under periodic boundary conditions. We discuss elements of our implementation of the PME algorithm on GPUs as well as pertinent performance issues. We proceed to show results of simulations of extended lipid bilayer systems using our program, FEN ZI. We performed simulations of DMPC bilayer systems consisting of 17,004, 68,484, and 273,936 atoms in explicit solvent. We present bilayer structural properties (atomic number densities, electron density profiles), deuterium order parameters (SCD), electrostatic properties (dipole potential, water dipole moments), and orientational properties of water. Predicted properties demonstrate excellent agreement with experiment and previous all‐atom molecular dynamics simulations. We observe no statistically significant differences in calculated structural or electrostatic properties for different system sizes, suggesting the small bilayer simulations (less than 100 lipid molecules) provide equivalent representation of structural and electrostatic properties associated with significantly larger systems (over 1000 lipid molecules). We stress that the three system size representations will have differences in other properties such as surface capillary wave dynamics or surface tension related effects that are not probed in the current study. The latter properties are inherently dependent on system size. This contribution suggests the suitability of applying emerging GPU technologies to studies of an important class of biological environments, that of lipid bilayers and their associated integral membrane proteins. We envision that this technology will push the boundaries of fully atomic‐resolution modeling of these biological systems, thus enabling unprecedented exploration of meso‐scale phenomena (mechanisms, kinetics, energetics) with atomic detail at commodity hardware prices. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011

[1]  Jeremy C. Smith,et al.  Mechanism and kinetics of peptide partitioning into membranes from all-atom simulations of thermostable peptides. , 2010, Journal of the American Chemical Society.

[2]  B. Roux,et al.  The ionization state and the conformation of Glu-71 in the KcsA K(+) channel. , 2002, Biophysical journal.

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

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

[5]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[6]  Baroni,et al.  Conjugate gradient minimization of the energy functional: A new method for electronic structure calculation. , 1989, Physical review. B, Condensed matter.

[7]  Michela Taufer,et al.  FENZI: GPU-Enabled Molecular Dynamics Simulations of Large Membrane Regions Based on the CHARMM Force Field and PME , 2011, 2011 IEEE International Symposium on Parallel and Distributed Processing Workshops and Phd Forum.

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

[9]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[10]  P. Clancy,et al.  A critical study of simulations of the Lennard-Jones liquid-vapor interface , 1992 .

[11]  D. Peter Tieleman,et al.  Modifying the OPLS‐AA force field to improve hydration free energies for several amino acid side chains using new atomic charges and an off‐plane charge model for aromatic residues , 2007, J. Comput. Chem..

[12]  D Peter Tieleman,et al.  Computer simulation of partitioning of ten pentapeptides Ace-WLXLL at the cyclohexane/water and phospholipid/water interfaces , 2005 .

[13]  Luca Monticelli,et al.  Effect of lipid peroxidation on the properties of lipid bilayers: a molecular dynamics study. , 2007, Biophysical journal.

[14]  S. Feller,et al.  Molecular dynamics simulations of lipid bilayers , 2000 .

[15]  Todd J. Martinez,et al.  Graphical Processing Units for Quantum Chemistry , 2008, Computing in Science & Engineering.

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

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

[18]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[19]  Sagar A. Pandit,et al.  Molecular dynamics simulation of a dipalmitoylphosphatidylcholine bilayer with NaCl. , 2003, Biophysical journal.

[20]  Eric Darve,et al.  N-Body simulation on GPUs , 2006, SC.

[21]  Sagar A. Pandit,et al.  Mixed bilayer containing dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylserine: lipid complexation, ion binding, and electrostatics. , 2003, Biophysical journal.

[22]  M. Klein,et al.  Molecular dynamics simulations of self-assembly and nanotube formation by amphiphilic molecules in aqueous solution: a coarse-grain approach , 2007 .

[23]  W. Im,et al.  Theoretical and computational models of biological ion channels , 2004, Quarterly Reviews of Biophysics.

[24]  Alexander D. MacKerell,et al.  An ab initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer. , 2005, The journal of physical chemistry. B.

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

[26]  Serdar Kuyucak,et al.  Gramicidin A channel as a test ground for molecular dynamics force fields. , 2003, Biophysical journal.

[27]  Klaus Schulten,et al.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

[28]  Michele Parrinello,et al.  Energy Conservation in Adaptive Hybrid Atomistic/Coarse-Grain Molecular Dynamics. , 2007, Journal of chemical theory and computation.

[29]  John E. Stone,et al.  Fast analysis of molecular dynamics trajectories with graphics processing units - Radial distribution function histogramming , 2011, J. Comput. Phys..

[30]  Michela Taufer,et al.  Molecular dynamics simulations of aqueous ions at the liquid–vapor interface accelerated using graphics processors , 2011, J. Comput. Chem..

[31]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[32]  Benoît Roux,et al.  Theoretical and computational models of ion channels. , 2002, Current opinion in structural biology.

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

[34]  Vijay S. Pande,et al.  Efficient nonbonded interactions for molecular dynamics on a graphics processing unit , 2010, J. Comput. Chem..

[35]  Joshua A. Anderson,et al.  General purpose molecular dynamics simulations fully implemented on graphics processing units , 2008, J. Comput. Phys..

[36]  Kyle A. Beauchamp,et al.  Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1-39). , 2010, Journal of the American Chemical Society.

[37]  B. Roux,et al.  Energetics of ion conduction through the gramicidin channel , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Wilfred F. van Gunsteren,et al.  An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase , 2001, J. Comput. Chem..

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

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

[41]  Sagar A. Pandit,et al.  Aqueous solutions next to phospholipid membrane surfaces: insights from simulations. , 2006, Chemical reviews.

[42]  Michael L Klein,et al.  Probing Membrane Insertion Activity of Antimicrobial Polymers via Coarse-grain Molecular Dynamics. , 2006, Journal of chemical theory and computation.

[43]  H. Berendsen,et al.  Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations. , 1996, Biophysical journal.

[44]  Stephen D. Bond,et al.  The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics , 1999 .

[45]  Klaus Schulten,et al.  Multilevel summation of electrostatic potentials using graphics processing units , 2009, Parallel Comput..

[46]  Harry A. Stern,et al.  Calculation of the dielectric permittivity profile for a nonuniform system: Application to a lipid bilayer simulation , 2003 .

[47]  Klaus Schulten,et al.  GPU-accelerated molecular modeling coming of age. , 2010, Journal of molecular graphics & modelling.

[48]  Sagar A. Pandit,et al.  Complexation of phosphatidylcholine lipids with cholesterol. , 2004, Biophysical journal.

[49]  J. Czub,et al.  Comparative molecular dynamics study of lipid membranes containing cholesterol and ergosterol. , 2006, Biophysical journal.

[50]  Vijay S. Pande,et al.  Hard Data on Soft Errors: A Large-Scale Assessment of Real-World Error Rates in GPGPU , 2009, 2010 10th IEEE/ACM International Conference on Cluster, Cloud and Grid Computing.

[51]  Klaus Schulten,et al.  Immersive Molecular Visualization and Interactive Modeling with Commodity Hardware , 2010, ISVC.

[52]  Vijay S. Pande,et al.  Accelerating molecular dynamic simulation on graphics processing units , 2009, J. Comput. Chem..

[53]  Sandeep Patel,et al.  Molecular dynamics simulations of a DMPC bilayer using nonadditive interaction models. , 2009, Biophysical journal.

[54]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[55]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .

[56]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .

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

[58]  G. Voth,et al.  Flexible simple point-charge water model with improved liquid-state properties. , 2006, The Journal of chemical physics.

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

[60]  Koji Yasuda,et al.  Two‐electron integral evaluation on the graphics processor unit , 2008, J. Comput. Chem..