Particle-based membrane model for mesoscopic simulation of cellular dynamics.

We present a simple and computationally efficient coarse-grained and solvent-free model for simulating lipid bilayer membranes. In order to be used in concert with particle-based reaction-diffusion simulations, the model is purely based on interacting and reacting particles, each representing a coarse patch of a lipid monolayer. Particle interactions include nearest-neighbor bond-stretching and angle-bending and are parameterized so as to reproduce the local membrane mechanics given by the Helfrich energy density over a range of relevant curvatures. In-plane fluidity is implemented with Monte Carlo bond-flipping moves. The physical accuracy of the model is verified by five tests: (i) Power spectrum analysis of equilibrium thermal undulations is used to verify that the particle-based representation correctly captures the dynamics predicted by the continuum model of fluid membranes. (ii) It is verified that the input bending stiffness, against which the potential parameters are optimized, is accurately recovered. (iii) Isothermal area compressibility modulus of the membrane is calculated and is shown to be tunable to reproduce available values for different lipid bilayers, independent of the bending rigidity. (iv) Simulation of two-dimensional shear flow under a gravity force is employed to measure the effective in-plane viscosity of the membrane model and show the possibility of modeling membranes with specified viscosities. (v) Interaction of the bilayer membrane with a spherical nanoparticle is modeled as a test case for large membrane deformations and budding involved in cellular processes such as endocytosis. The results are shown to coincide well with the predicted behavior of continuum models, and the membrane model successfully mimics the expected budding behavior. We expect our model to be of high practical usability for ultra coarse-grained molecular dynamics or particle-based reaction-diffusion simulations of biological systems.

[1]  H. Noguchi,et al.  Shape transitions of fluid vesicles and red blood cells in capillary flows. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Reinhard Lipowsky,et al.  Tubulation and aggregation of spherical nanoparticles adsorbed on vesicles. , 2012, Physical review letters.

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

[4]  Gerhard Hummer,et al.  Formation and Stability of Lipid Membrane Nanotubes. , 2017, ACS nano.

[5]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[6]  R. Pastor,et al.  Structure and elasticity of lipid membranes with genistein and daidzein bioflavinoids using X-ray scattering and MD simulations. , 2012, The journal of physical chemistry. B.

[7]  P. T. Wolde,et al.  Simulating biochemical networks at the particle level and in time and space: Green's function reaction dynamics. , 2005 .

[8]  M. Deserno,et al.  Determining the Gaussian curvature modulus of lipid membranes in simulations. , 2012, Biophysical journal.

[9]  J. Nagle,et al.  Comparing simulations of lipid bilayers to scattering data: the GROMOS 43A1-S3 force field. , 2013, The journal of physical chemistry. B.

[10]  Hongyan Yuan,et al.  One-particle-thick, solvent-free, coarse-grained model for biological and biomimetic fluid membranes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Feng Feng,et al.  Finite element modeling of lipid bilayer membranes , 2006, J. Comput. Phys..

[12]  Brian A. Camley,et al.  Lipid bilayer domain fluctuations as a probe of membrane viscosity. , 2010, Biophysical journal.

[13]  Gregory A Voth,et al.  A second generation mesoscopic lipid bilayer model: connections to field-theory descriptions of membranes and nonlocal hydrodynamics. , 2006, The Journal of chemical physics.

[14]  M. Tuckerman,et al.  A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal–isobaric ensemble , 2006 .

[15]  R. Lipowsky,et al.  Cooperative wrapping of nanoparticles by membrane tubes. , 2014, Soft matter.

[16]  Oliver F. Lange,et al.  Recognition Dynamics Up to Microseconds Revealed from an RDC-Derived Ubiquitin Ensemble in Solution , 2008, Science.

[17]  Rumiana Dimova,et al.  Recent developments in the field of bending rigidity measurements on membranes. , 2014, Advances in colloid and interface science.

[18]  P. R. ten Wolde,et al.  Combining molecular dynamics with mesoscopic Green's function reaction dynamics simulations. , 2015, The Journal of chemical physics.

[19]  Reinhard Lipowsky,et al.  Equilibrium structure and lateral stress distribution of amphiphilic bilayers from dissipative particle dynamics simulations , 2002 .

[20]  Anand Srivastava,et al.  Linear aggregation of proteins on the membrane as a prelude to membrane remodeling , 2013, Proceedings of the National Academy of Sciences.

[21]  J. Nagle.  Introductory lecture: basic quantities in model biomembranes. , 2013, Faraday discussions.

[22]  M. Deserno,et al.  Determining the bending modulus of a lipid membrane by simulating buckling. , 2013, The Journal of chemical physics.

[23]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[24]  Krassimir D. Danov,et al.  Falling ball viscosimetry of giant vesicle membranes: Finite-size effects , 1999 .

[25]  Kurt Kremer,et al.  Tunable generic model for fluid bilayer membranes. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  G. Gompper,et al.  Shape and orientation matter for the cellular uptake of nonspherical particles. , 2014, Nano letters.

[27]  G. Gompper,et al.  Wrapping of ellipsoidal nano-particles by fluid membranes , 2013, 1303.5567.

[28]  Berend Smit,et al.  Mesoscopic models of biological membranes , 2006 .

[29]  M. Laradji,et al.  Implicit-solvent mesoscale model based on soft-core potentials for self-assembled lipid membranes. , 2008, The Journal of chemical physics.

[30]  U. Schwarz,et al.  Studying protein assembly with reversible Brownian dynamics of patchy particles. , 2014, The Journal of chemical physics.

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

[32]  W. Helfrich.  Elastic Properties of Lipid Bilayers: Theory and Possible Experiments , 1973, Zeitschrift fur Naturforschung. Teil C: Biochemie, Biophysik, Biologie, Virologie.

[33]  Reinhard Lipowsky,et al.  The computational route from bilayer membranes to vesicle fusion , 2006, Journal of physics. Condensed matter : an Institute of Physics journal.

[34]  F. Noé,et al.  Protein conformational plasticity and complex ligand-binding kinetics explored by atomistic simulations and Markov models , 2015, Nature Communications.

[35]  Jaime Agudo-Canalejo,et al.  Wrapping of nanoparticles by membranes. , 2014, Advances in colloid and interface science.

[36]  Michael R. Shirts,et al.  Dynamical reweighting: improved estimates of dynamical properties from simulations at multiple temperatures. , 2011, The Journal of chemical physics.

[37]  Siewert J Marrink,et al.  Lipids on the move: simulations of membrane pores, domains, stalks and curves. , 2009, Biochimica et biophysica acta.

[38]  F. Noé,et al.  Explicit spatiotemporal simulation of receptor-G protein coupling in rod cell disk membranes. , 2014, Biophysical journal.

[39]  E. Neher,et al.  Protein scaffolds in the coupling of synaptic exocytosis and endocytosis , 2011, Nature Reviews Neuroscience.

[40]  Frank Noé,et al.  Simulation tools for particle-based reaction-diffusion dynamics in continuous space , 2014, BMC biophysics.

[41]  Helgi I Ingólfsson,et al.  Dry Martini, a coarse-grained force field for lipid membrane simulations with implicit solvent. , 2015, Journal of chemical theory and computation.

[42]  Subra Suresh,et al.  Size‐Dependent Endocytosis of Nanoparticles , 2009, Advanced materials.

[43]  H. Noguchi.  Anisotropic surface tension of buckled fluid membranes. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  A. Mikhailov,et al.  Coarse-grain model for lipid bilayer self-assembly and dynamics: multiparticle collision description of the solvent. , 2012, The Journal of chemical physics.

[45]  Frank Noé,et al.  Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states. , 2014, The Journal of chemical physics.

[46]  R. Altman,et al.  Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways , 2013, Nature chemistry.

[47]  Gregory A Voth,et al.  New insights into BAR domain-induced membrane remodeling. , 2009, Biophysical journal.

[48]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[49]  O. Edholm,et al.  Areas of molecules in membranes consisting of mixtures. , 2005, Biophysical journal.

[50]  Hiroshi Noguchi,et al.  Dynamics of fluid vesicles in shear flow: effect of membrane viscosity and thermal fluctuations. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[51]  L Donati,et al.  Girsanov reweighting for path ensembles and Markov state models. , 2017, The Journal of chemical physics.

[52]  Gregory A Voth,et al.  Bridging microscopic and mesoscopic simulations of lipid bilayers. , 2002, Biophysical journal.

[53]  D. Tieleman,et al.  Perspective on the Martini model. , 2013, Chemical Society reviews.

[54]  Satoru Yamamoto,et al.  Budding and fission dynamics of two-component vesicles , 2003 .

[55]  Frank Noé,et al.  Dynamical Organization of Syntaxin-1A at the Presynaptic Active Zone , 2015, PLoS Comput. Biol..

[56]  M. Parrinello,et al.  Kinetics of protein–ligand unbinding: Predicting pathways, rates, and rate-limiting steps , 2015, Proceedings of the National Academy of Sciences.

[57]  Frank Noé,et al.  ReaDDyMM: Fast interacting particle reaction-diffusion simulations using graphical processing units. , 2015, Biophysical journal.

[58]  P. Tarazona,et al.  A computer simulation approach to quantify the true area and true area compressibility modulus of biological membranes. , 2015, The Journal of chemical physics.

[59]  D. Marsh,et al.  Elastic curvature constants of lipid monolayers and bilayers. , 2006, Chemistry and physics of lipids.

[60]  Frank Noé,et al.  ReaDDy - A Software for Particle-Based Reaction-Diffusion Dynamics in Crowded Cellular Environments , 2013, PloS one.

[61]  Frank Noé,et al.  MHC class II complexes sample intermediate states along the peptide exchange pathway , 2016, Nature Communications.

[62]  Gerhard Gompper,et al.  Network models of fluid, hexatic and polymerized membranes , 1997 .

[63]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[64]  Gregory A Voth,et al.  Protein-mediated transformation of lipid vesicles into tubular networks. , 2013, Biophysical journal.

[65]  Hiroshi Noguchi,et al.  Membrane Simulation Models from Nanometer to Micrometer Scale , 2008, 0812.0055.

[66]  Solvent-free model for self-assembling fluid bilayer membranes: stabilization of the fluid phase based on broad attractive tail potentials. , 2005, The Journal of chemical physics.

[67]  Watt W. Webb,et al.  Large-scale fluid/fluid phase separation of proteins and lipids in giant plasma membrane vesicles , 2007, Proceedings of the National Academy of Sciences.

[68]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[69]  Anand Srivastava,et al.  A Hybrid Approach for Highly Coarse-grained Lipid Bilayer Models. , 2013, Journal of chemical theory and computation.

[70]  E. Lindahl,et al.  Mesoscopic undulations and thickness fluctuations in lipid bilayers from molecular dynamics simulations. , 2000, Biophysical journal.

[71]  Max Hoffmann,et al.  Oscillations of Min-proteins in micropatterned environments: a three-dimensional particle-based stochastic simulation approach. , 2014, Soft matter.

[72]  G. Gompper,et al.  Fluctuations of polymerized, fluid and hexatic membranes: Continuum models and simulations , 1997 .

[73]  T. Kohyama.  Simulations of flexible membranes using a coarse-grained particle-based model with spontaneous curvature variables , 2009 .

[74]  Frank Noé,et al.  Lipid-mediated PX-BAR domain recruitment couples local membrane constriction to endocytic vesicle fission , 2017, Nature Communications.

[75]  Markus Deserno,et al.  Elastic deformation of a fluid membrane upon colloid binding. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[76]  M. Parrinello,et al.  From metadynamics to dynamics. , 2013, Physical review letters.

[77]  Frank Noé,et al.  HTMD: High-Throughput Molecular Dynamics for Molecular Discovery. , 2016, Journal of chemical theory and computation.

[78]  Hao Wu,et al.  Multiensemble Markov models of molecular thermodynamics and kinetics , 2016, Proceedings of the National Academy of Sciences.

[79]  Sean X. Sun,et al.  Shape transitions in lipid membranes and protein mediated vesicle fusion and fission. , 2007, The Journal of chemical physics.

[80]  G. Gompper,et al.  Statistical mechanics of membranes: freezing, undulations, and topology fluctuations , 2000 .

[81]  P. Canham.  The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell. , 1970, Journal of theoretical biology.

[82]  Petra Schwille,et al.  Translational diffusion in lipid membranes beyond the Saffman-Delbruck approximation. , 2008, Biophysical journal.

[83]  Udo Seifert,et al.  Configurations of fluid membranes and vesicles , 1997 .

[84]  Gregory A Voth,et al.  Interfacing molecular dynamics and macro-scale simulations for lipid bilayer vesicles. , 2002, Biophysical journal.

[85]  A. Gorfe,et al.  Simulating POPC and POPC/POPG Bilayers: Conserved Packing and Altered Surface Reactivity. , 2010, Journal of chemical theory and computation.

[86]  O. Farago.  “Water-free” computer model for fluid bilayer membranes , 2003, cond-mat/0304203.

[87]  Markus Deserno,et al.  Mesoscopic membrane physics: concepts, simulations, and selected applications. , 2009, Macromolecular rapid communications.

[88]  Marissa G. Saunders,et al.  Coarse-graining methods for computational biology. , 2013, Annual review of biophysics.

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

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

[91]  B. Berne,et al.  How and when does an anticancer drug leave its binding site? , 2017, Science Advances.

[92]  Daan Frenkel,et al.  Modeling flexible amphiphilic bilayers: a solvent-free off-lattice Monte Carlo study. , 2005, The Journal of chemical physics.

[93]  Alan E. Mark,et al.  Effect of Undulations on Surface Tension in Simulated Bilayers , 2001 .

[94]  Gregory A. Voth,et al.  Membrane tension controls the assembly of curvature-generating proteins , 2015, Nature Communications.

[95]  A C Maggs,et al.  Computer simulations of self-assembled membranes. , 1991, Science.

[96]  Reinhard Lipowsky,et al.  Computer simulations of bilayer membranes - self-assembly and interfacial tension. , 1998 .

[97]  T. Frankel.  The Geometry of Physics , 1997 .

[98]  Friederike Schmid,et al.  A Simple Computer Model for Liquid Lipid Bilayers , 2005 .

[99]  B. Różycki,et al.  Spontaneous curvature of bilayer membranes from molecular simulations: asymmetric lipid densities and asymmetric adsorption. , 2015, The Journal of chemical physics.

[100]  E. Evans,et al.  Bending resistance and chemically induced moments in membrane bilayers. , 1974, Biophysical journal.

[101]  F. Noé,et al.  Complete protein–protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modelling , 2017, Nature Chemistry.

[102]  Pietro Cicuta,et al.  Diffusion of liquid domains in lipid bilayer membranes. , 2007, The journal of physical chemistry. B.

[103]  Frank Noé,et al.  Higher-order architecture of rhodopsin in intact photoreceptors and its implication for phototransduction kinetics. , 2015, Structure.

[104]  Gregory A Voth,et al.  Membrane remodeling from N-BAR domain interactions: insights from multi-scale simulation. , 2007, Biophysical journal.

[105]  Y. Sugita,et al.  Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm , 2016, eLife.

[106]  R. Erban,et al.  Stochastic modelling of reaction–diffusion processes: algorithms for bimolecular reactions , 2009, Physical biology.

[107]  Klaus Schulten,et al.  Structural mechanism of voltage-dependent gating in an isolated voltage-sensing domain , 2013, Nature Structural &Molecular Biology.