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.

Prediction, Refinement, and Persistency of Transmembrane Helix Dimers in Lipid Bilayers Using Implicit and Explicit Solvent/lipid Representations: Microsecond Molecular Dynamics Simulations of ErbB1/B2 and EphA1.

Molecular Dynamics Simulations Of Forced Unbending Of Integrin Alpha(V)Beta(3)

Molecular dynamics simulations of forced unbending of integrin α(v)β₃.

Electronic Polarization Stabilizes Tertiary Structure Prediction of HP-36.

All-Atom Molecular Dynamics Elucidating Molecular Mechanisms of Single-Transmembrane Model Peptide Dimerization in a Lipid Bilayer

Polymorphism of the Epidermal Growth Factor Receptor Extracellular Ligand Binding Domain: the Dimer Interface Depends on Domain Stabilization

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Evaluating the Dynamics and Electrostatic Interactions of Folded Proteins in Implicit Solvents

Investigation on Structure and Dynamics of Heterodimeric SAM:SAM Complexes Using Micro-Second Molecular Dynamics Simulations

Balancing Force Field Protein–Lipid Interactions To Capture Transmembrane Helix–Helix Association

Free energy, rates, and mechanism of transmembrane dimerization in lipid bilayers from dynamically unbiased molecular dynamics simulations

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

Exploring The Structure And Dynamics Of All-Atom Models For The Plexin Transmembrane Receptor Bound To Gtpases And To Lipid Bilayer

Architecture and Membrane Interactions of the EGF Receptor

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

Improved Efficiency of Replica Exchange Simulations through Use of a Hybrid Explicit/Implicit Solvation Model

Dissipative Particle Dynamics Simulation of Protein Folding in Explicit and Implicit Solvents: Coarse-Grained Model for Atomic Resolution

Efficient Exploration of Membrane-Associated Phenomena at Atomic Resolution

Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

Structural Refinement of Proteins by Restrained Molecular Dynamics Simulations with Non-Interacting Molecular Fragments

Examination Of Membrane Fusion By Dissipative Particle Dynamics Simulation And Comparison With Continuum Elastic Models