Abstract:Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Binary bilayer simulations for partitioning within membranes

Understanding the free-energy landscape of phase separation in lipid bilayers using molecular dynamics

Dynamic simulations of multicomponent lipid membranes over long length and time scales

System size effects on the free energy landscapes from molecular dynamics of phase-separating bilayers

Large-scale simulations of fluctuating biological membranes

Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

Computational and theoretical approaches for studies of a lipid recognition protein on biological membranes

Dynamics of phase separation in mixed lipid membranes between two bounding walls

Computer simulations of a heterogeneous membrane with enhanced sampling techniques

Kinetics of symmetry and asymmetry in a phase-separating bilayer membrane

Identifying the Onset of Phase Separation in Quaternary Lipid Bilayer Systems from Coarse-Grained Simulations

Experimental models for dynamic compartmentalization of biomolecules in liquid organelles: Reversible formation and partitioning in aqueous biphasic systems

Multiphase model membraneless organelles

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

Computational studies of biomembrane systems: Theoretical considerations, simulation models, and applications

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Characterization of domain formation in complex membranes

A computational study of lateral phase separation in biological membranes

Exploring Free Energy Landscapes for Protein Partitioning into Membrane Domains in All-Atom and Coarse-Grained Simulations

Partitioning of benzoic acid into DMPC and blood-brain barrier mimetic bilayers.

Setting Up All-Atom Molecular Dynamics Simulations to Study the Interactions of Peripheral Membrane Proteins with Model Lipid Bilayers