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.

Faster but not sweeter: A model of E. coli Re-level lipopolysaccharide for Martini 3 and a Martini 2 version with accelerated kinetics

Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli

Developing a Coarse-Grained Model for Bacterial Cell Walls and Evaluating Mechanical Properties and Free Energy Barriers

Lipid Flipping, O-antigen Shielding, and Water Dynamics Revealed by 100-200 μs Coarse-Grained Simulations of E. coli Outer Membrane Proteins

Modeling of Specific Lipopolysaccharide Binding Sites on a Gram-Negative Porin

Structures, Dynamics, and Water Permeation Free Energy Across Bilayers of Lipid A and Its Analog Studied with Molecular Dynamics Simulation

Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Capsular Polysaccharides

Atomistic Scale Effects of Lipopolysaccharide Modifications on Bacterial Outer Membrane Defenses

Single-molecule dynamics show a transient lipopolysaccharide transport bridge

A minimalist model lipid system mimicking the biophysical properties of Escherichia coli membrane

Functional Analysis of the Protein Machinery Required for Transport of Lipopolysaccharide to the Outer Membrane of Escherichia coli

A computational model for lipid-anchored polysaccharide export by the outer membrane protein GfcD

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Bacterial extracellular vesicles: towards realistic models for bacterial membranes in molecular interaction studies by surface plasmon resonance

Recalibration of MARTINI-3 Parameters for Improved Interactions between Peripheral Proteins and Lipid Bilayers

Optimized Protein-Excipient Interactions in the Martini 3 Force Field

Molecular interactions, elastic properties, and nanostructure of Langmuir bacterial-lipid monolayers: Towards solving the mystery in bacterial membrane asymmetry

Martini 3 building blocks for Lipid Nanoparticle design

Outer Membranes of Polymyxin-Resistant Acinetobacter Baumannii with Phosphoethanolamine-Modified Lipid A and Lipopolysaccharide Loss Display Different Atomic-Scale Interactions with Polymyxins

Full-Atom Model of the Agonist LPS-Bound Toll-like Receptor 4 Dimer in a Membrane Environment

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