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

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

Apatites in Gale Crater

ChemCam is an active remote sensing instrument suite that has operated successfully on MSL since landing Aug. 6th, 2012. It uses laser pulses to remove dust and to analyze rocks up to 7 m away. Laser-induced breakdown spectroscopy (LIBS) obtains emission spectra of materials ablated from the samples in electronically excited states. The intensities of the emission lines scale with the abundances of the related element. ChemCam is sensitive to most major rock-forming elements as well as to a set of minor and trace elements such as F, Cl, Li, P, Sr, Ba, and Rb. The measured chemical composition can then be used to infer the mineralogical composition of the ablated material. Here, we report a summary of inferred apatite detections along the MSL traverse at Gale Crater. We present the geologic settings of these findings and derive some interpretations about the formation conditions of apatite in time and space

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

Vibrational Spectra of a Mechanosensitive Channel

We report the simulated vibrational spectra of a mechanosensitive membrane channel in different gating states. Our results show that while linear absorption is insensitive to structural differences, linear dichroism and sum-frequency generation spectroscopies are sensitive to the orientation of the transmembrane helices, which is changing during the opening process. Linear dichroism cannot distinguish an intermediate structure from the closed structure, but sum-frequency generation can. In addition, we find that two-dimensional infrared spectroscopy can be used to distinguish all three investigated gating states of the mechanosensitive membrane channel.

Calibration of the MSL/ChemCam/LIBS Remote Sensing Composition Instrument

The ChemCam instrument suite on board the 2011 Mars Science Laboratory (MSL) Rover, Curiosity, will provide remote-sensing composition information for rock and soil samples within seven meters of the rover using a laser-induced breakdown spectroscopy (LIBS) system, and will provide context imaging with a resolution of 0.10 mradians using the remote micro-imager (RMI) camera. The high resolution is needed to image the small analysis footprint of the LIBS system, at 0.2-0.6 mm diameter. This fine scale analytical capability will enable remote probing of stratigraphic layers or other small features the size of "blueberries" or smaller. ChemCam is intended for rapid survey analyses within 7 m of the rover, with each measurement taking less than 6 minutes. Repeated laser pulses remove dust coatings and provide depth profiles through weathering layers, allowing detailed investigation of rock varnish features as well as analysis of the underlying pristine rock composition. The LIBS technique uses brief laser pulses greater than 10 MW/square mm to ablate and electrically excite material from the sample of interest. The plasma emits photons with wavelengths characteristic of the elements present in the material, permitting detection and quantification of nearly all elements, including the light elements H, Li, Be, B, C, N, O. ChemCam LIBS projects 14 mJ of 1067 nm photons on target and covers a spectral range of 240-850 nm with resolutions between 0.15 and 0.60 nm FWHM. The Nd:KGW laser is passively cooled and is tuned to provide maximum power output from -10 to 0 C, though it can operate at 20% degraded energy output at room temperature. Preliminary calibrations were carried out on the flight model (FM) in 2008. However, the detectors were replaced in 2009, and final calibrations occurred in April-June, 2010. This presentation describes the LIBS calibration and characterization procedures and results, and details plans for final analyses during rover system thermal testing, planned for early March

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

Pain as a global public health priority

<p>Abstract</p> <p>Background</p> <p>Pain is an enormous problem globally. Estimates suggest that 20% of adults suffer from pain globally and 10% are newly diagnosed with chronic pain each year. Nevertheless, the problem of pain has primarily been regarded as a medical problem, and has been little addressed by the field of public health.</p> <p>Discussion</p> <p>Despite the ubiquity of pain, whether acute, chronic or intermittent, public health scholars and practitioners have not addressed this issue as a public health problem. The importance of viewing pain through a public health lens allows one to understand pain as a multifaceted, interdisciplinary problem for which many of the causes are the social determinants of health. Addressing pain as a global public health issue will also aid in priority setting and formulating public health policy to address this problem, which, like most other chronic non-communicable diseases, is growing both in absolute numbers and in its inequitable distribution across the globe.</p> <p>Summary</p> <p>The prevalence, incidence, and vast social and health consequences of global pain requires that the public health community give due attention to this issue. Doing so will mean that health care providers and public health professionals will have a more comprehensive understanding of pain and the appropriate public health and social policy responses to this problem.</p