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.

Accuracy and Reproducibility of Lipari-Szabo Order Parameters From Molecular Dynamics

Learning Effective Molecular Models from Experimental Observables.

Order parameter prediction from molecular dynamics simulations in proteins

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Probing Static Disorder In Arrhenius Kinetics By Single-Molecule Force Spectroscopy

Definition of Systematic, Approximately Separable and Modular Internal Coordinates (SASMIC) for macromolecular simulation

Measuring Intrinsic Disorder and Tracking Conformational Transitions Using Rosetta ResidueDisorder

Improving the accuracy of pose prediction by incorporating symmetry-related molecules

An Approximate Method in Using Molecular Mechanics Simulations to Study Slow Protein Conformational Changes

Bayesian nonparametric analysis of residence times for protein-lipid interactions in Molecular Dynamics simulations

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

Condensed-phase molecular representation to link structure and thermodynamics in molecular dynamics

Accurate p K a Calculations in Proteins with Reactive Molecular Dynamics Provide Physical Insight Into the Electrostatic Origins of Their Values

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Impact of replicas and simulation length on in silico behaviors of a protein domain

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

Adaptive Resolution Force Probe Simulations: Coarse Graining in the Ideal Gas Approximation

Prediction of Nmr Order Parameters in Proteins Using Weighted Protein Contact-Number Model

Developing a molecular dynamics force field for both folded and disordered protein states