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

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields’ Descriptions of POPC–Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C–H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches

Integrative determination of the atomic structure of mutant huntingtin exon 1 fibrils from Huntington's disease

Neurodegeneration in Huntington's disease (HD) is accompanied by the aggregation of fragments of the mutant huntingtin protein, a biomarker of disease progression. A particular pathogenic role has been attributed to the aggregation-prone huntingtin exon 1 (HttEx1) fragment, whose polyglutamine (polyQ) segment is expanded. Unlike amyloid fibrils from Parkinson's and Alzheimer's diseases, the atomic-level structure of HttEx1 fibrils has remained unknown, limiting diagnostic and treatment efforts. We present and analyze the structure of fibrils formed by polyQ peptides and polyQ-expanded HttEx1. Atomic-resolution perspectives are enabled by an integrative analysis and unrestrained all-atom molecular dynamics (MD) simulations incorporating experimental data from electron microscopy (EM), solid-state NMR, and other techniques. Visualizing the HttEx1 subdomains in atomic detail helps explaining the biological properties of these protein aggregates, as well as paves the way for targeting them for detection and degradation.</p

3D Scaffolds of Polycaprolactone/Copper-Doped Bioactive Glass: Architecture Engineering with Additive Manufacturing and CellularAssessments in a Coculture of Bone Marrow Stem Cells and Endothelial Cells

The local delivery of Cu2+ from copper-doped bioactive glass (Cu-BaG) was combined with 3D printing of polycaprolactone (PCL) scaffolds for its potent angiogenic effect in bone tissue engineering. PCL and Cu-BaG were, respectively, dissolved and dispersed in acetone to formulate a moderately homogeneous ink. The PCL/Cu-BaG scaffolds were fabricated via direct ink writing into a cold ethanol bath. The architecture of the printed scaffolds, including strut diameter, strut spacing, and porosity, were investigated and characterized. The PCL/Cu-BaG scaffolds showed a Cu-BaG content-dependent mechanical property, as the compressive Young\u27s modulus ranged from 7 to 13 MPa at an apparent porosity of 60%. The ion dissolution behavior in simulated body fluid was evaluated, and the hydroxyapatite-like precipitation on the strut surface was confirmed. Furthermore, the cytocompatibility of the PCL/Cu-BaG scaffolds was assessed in human bone marrow stem cell (hBMSC) culture, and a dose-dependent cytotoxicity of Cu2+ was observed. Here, the PCL/BaG scaffold induced the higher expression of late osteogenic genes OSTEOCALCIN and DLX5 in comparison to the PCL scaffold. The doping of Cu2+ in BaG elicited higher expression of the early osteogenic marker gene RUNX2a but decreased the expression of late osteogenic marker genes OSTEOCALCIN and DLX5 in comparison to the PCL/BaG scaffold, demonstrating the suppressing effect of Cu2+ on osteogenic differentiation of hBMSCs. In a coculture of hBMSCs and human umbilical vein endothelial cells, both the PCL/BaG and PCL/Cu-BaG scaffolds stimulated the formation of a denser tubule network, compared to the PCL scaffold. Meanwhile, only slightly higher gene expression of vWF was observed with the PCL/Cu-BaG scaffold than with the PCL/BaG scaffold, indicating the potent angiogenic effect of the released Cu2+

Molecular electrometer and binding of cations to phospholipid bilayers

Despite the vast amount of experimental and theoretical studies on the binding affinity of cations -especially the biologically relevant Na+ and Ca2+ - for phospholipid bilayers, there is no consensus in the literature. Here we show that by interpreting changes in the choline headgroup order parameters according to the 'molecular electrometer' concept [Seelig et al., Biochemistry, 1987, 26, 7535], one can directly compare the ion binding affinities between simulations and experiments. Our findings strongly support the view that in contrast to Ca2+ and other multivalent ions, Na+ and other monovalent ions (except Li+) do not specifically bind to phosphatidylcholine lipid bilayers at sub-molar concentrations. However, the Na+ binding affinity was overestimated by several molecular dynamics simulation models, resulting in artificially positively charged bilayers and exaggerated structural effects in the lipid headgroups. While qualitatively correct headgroup order parameter response was observed with Ca2+ binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca2+: lipid stoichiometry or the induced atomistic resolution structural changes. All scientific contributions to this open collaboration work were made publicly, using nmrlipids. blogspot.fi as the main communication platform.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