Using Open Data to Rapidly Benchmark Biomolecular Simulations : Phospholipid Conformational Dynamics

Molecular dynamics (MD) simulations are widely used to monitor time-resolved motions of biomacromolecules, although it often remains unknown how closely the conformational dynamics correspond to those occurring in real life. Here, we used a large set of open-access MD trajectories of phosphatidylcholine (PC) lipid bilayers to benchmark the conformational dynamics in several contemporary MD models (force fields) against nuclear magnetic resonance (NMR) data available in the literature: effective correlation times and spin-lattice relaxation rates. We found none of the tested MD models to fully reproduce the conformational dynamics. That said, the dynamics in CHARMM36 and Slipids are more realistic than in the Amber Lipid14, OPLS-based MacRog, and GROMOS-based Berger force fields, whose sampling of the glycerol backbone conformations is too slow. The performance of CHARMM36 persists when cholesterol is added to the bilayer, and when the hydration level is reduced. However, for conformational dynamics of the PC headgroup, both with and without cholesterol, Slipids provides the most realistic description because CHARMM36 overestimates the relative weight of similar to 1 ns processes in the headgroup dynamics. We stress that not a single new simulation was run for the present work. This demonstrates the worth of open-access MD trajectory databanks for the indispensable step of any serious MD study: benchmarking the available force fields. We believe this proof of principle will inspire other novel applications of MD trajectory databanks and thus aid in developing biomolecular MD simulations into a true computational microscope-not only for lipid membranes but for all biomacromolecular systems.Peer reviewe

Emerging Era of Biomolecular Membrane Simulations: Automated Physically-Justified Force Field Development and Quality-Evaluated Databanks

Peer reviewe

Headgroup Structure and Cation Binding in Phosphatidylserine Lipid Bilayers

Phosphatidylserine (PS) is a negatively charged lipid type commonly found in eukaryotic membranes, where it interacts with proteins via nonspecific electrostatic interactions as well as via specific binding. Moreover, in the presence of calcium ions, PS lipids can induce membrane fusion and phase separation. Molecular details of these phenomena remain poorly understood, partly because accurate models to interpret the experimental data have not been available. Here we gather a set of previously published experimental NMR data of C-H bond order parameter magnitudes, vertical bar S-CH vertical bar, for pure PS and mixed PS:PC (phosphatidylcholine) lipid bilayers and augment this data set by measuring the signs of S-CH in the PS headgroup using S-DROSS solid-state NMR spectroscopy. The augmented data set is then used to assess the accuracy of the PS headgroup structures in, and the cation binding to, PS-containing membranes in the most commonly used classical molecular dynamics (MD) force fields including CHARMM36, Lipidl7, MacRog, Slipids, GROMOS-CKP, Berger, and variants. We show large discrepancies between different force fields and that none of them reproduces the NMR data within experimental accuracy. However, the best MD models can detect the most essential differences between PC and PS headgroup structures. The cation binding affinity is not captured correctly by any of the PS force fields-an observation that is in line with our previous results for PC lipids. Moreover, the simulated response of the PS headgroup to bound ions can differ from experiments even qualitatively. The collected experimental data set and simulation results will pave the way for development of lipid force fields that correctly describe the biologically relevant negatively charged membranes and their interactions with ions. This work is part of the NMRlipids open collaboration project (nmrlipids.blogspot.fi).Peer reviewe

Chemistry specificity of DNA-polycation complex salt response: A simulation study of DNA, polylysine and polyethyleneimine

Peer reviewe

Ewald electrostatics for mixtures of point and continuous line charge

Peer reviewe

Polyelectrolyte Decomplexation via Addition of Salt: Charge Correlation Driven Zipper

We report the first atomic scale studies of polyelectrolyte decomplexation. The complex between DNA and polylysine is shown to destabilize and spontaneously open in a gradual, reversible zipper-like mechanism driven by an increase in solution salt concentration. Divalent CaCl₂ is significantly more effective than monovalent NaCl in destabilizing the complex due to charge correlations and water binding capability. The dissociation occurs accompanied by charge reversal in which charge correlations and ion binding chemistry play a key role. Our results are in agreement with experimental work on complex dissociation but in addition show the underlying microstructural correlations driving the behavior. Comparison of our full atomic level detail and dynamics results with theoretical works describing the PEs as charged, rigid rods reveals that although charge correlation involved theories provide qualitatively similar responses, considering also specific molecular chemistry and molecular level water contributions provides a more complete understanding of PE complex stability and dynamics. The findings may facilitate controlled release in gene delivery and more in general tuning of PE membrane permeability and mechanical characteristics through ionic strength

Interaction modes between asymmetrically and oppositely charged rods

The interaction of oppositely and asymmetrically charged rods in salt - a simple model of (bio)macromolecular assembly - is observed via simulation to exhibit two free energy minima, separated by a repulsive barrier. In contrast to similar minima in the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, the governing mechanism includes electrostatic attraction at large separation, osmotic repulsion at close range, and depletion attraction near contact. A model accounting for ion condensation and excluded volume is shown to be superior to a mean-field treatment in predicting the effect of charge asymmetry on the free-energy profile.Peer reviewe