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.

[Cervical management in the N0 squamous cell carcinoma of the tongue].

Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types.

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

Evaluating polarizable biomembrane simulations against experiments

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field.

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

Lipid21: Complex Lipid Membrane Simulations with AMBER

Drude Polarizable Force Field Parametrization of Carboxylate and N-acetyl Amine Carbohydrate Derivatives.

PhyNEO: A Neural-Network-Enhanced Physics-Driven Force Field Development Workflow for Bulk Organic Molecule and Polymer Simulations.

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A Polarizable Valence Electron Density Based Force Field for High-Energy Interactions between Atoms and Molecules

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Physically motivated, robust, ab initio force fields for CO2 and N2.

CHARMM‐GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Balancing Group 1 Monoatomic Ion-Polar Compound Interactions in the Polarizable Drude Force Field: Application in Protein and Nucleic Acid Systems

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

Water-phospholipid interactions at the interface of lipid membranes: comparison of different force fields

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields