Laterally Resolved Small-Angle Scattering Intensity from Lipid Bilayer Simulations: An Exact and a Limited-Range Treatment

SANS Spectra with PLUMED: Implementation and Application to Metainference

Using small-angle scattering with either X-ray or neutron sources has become common in the investigation of soft-matter systems. These experiments provide information about the coarse shape of the scattered objects, but obtaining more-detailed information can usually only be achieved with the aid of molecular simulations. In this Application Note, we report the implementation of an extension in PLUMED to compute the small-angle neutron scattering (SANS), which can be used for data processing as well for enhanced sampling, in particular with the metainference method to bias simulations and sample structures with a resulting spectrum in agreement with an experimental reference. Our implementation includes a resolution function that can be used to smear the SANS intensities according to beamline error sources and is compatible with both all-atom and coarse-grained simulations. Scripts to aid in the calculation of the scattering lengths when the system is coarse-grained and to aid in preparing the inputs are provided. We illustrate the use of the implementation with metainference by performing coarse-grained simulations of beta-octylglucoside and dodecylphosphocholine micelles in water. With different software and different Hamiltonians, we show that the metainference SANS bias can drive micelles to be split and to change shapes to achieve a better agreement with the experimental reference.

Cinematic reflectometry using QIKR, the quite intense kinetics reflectometer

The Quite Intense Kinetics Reflectometer (QIKR) will be a general-purpose, horizontal-sample-surface neutron reflectometer. Reflectometers measure the proportion of an incident probe beam reflected from a surface as a function of wavevector (momentum) transfer to infer the distribution and composition of matter near an interface. The unique scattering properties of neutrons make this technique especially useful in the study of soft matter, biomaterials, and materials used in energy storage. Exploiting the increased brilliance of the Spallation Neutron Source Second Target Station, QIKR will collect specular and off-specular reflectivity data faster than the best existing such machines. It will often be possible to collect complete specular reflectivity curves using a single instrument setting, enabling "cinematic" operation, wherein the user turns on the instrument and "films" the sample. Samples in time-dependent environments (e.g., temperature, electrochemical, or undergoing chemical alteration) will be observed in real time, in favorable cases with frame rates as fast as 1 Hz. Cinematic data acquisition promises to make time-dependent measurements routine, with time resolution specified during post-experiment data analysis. This capability will be deployed to observe such processes as in situ polymer diffusion, battery electrode charge-discharge cycles, hysteresis loops, and membrane protein insertion into lipid layers.

The dynamic face of lipid membranes

Cell membranes - primarily composed of lipids, sterols, and proteins - form a dynamic interface between living cells and their environment. They act as a mechanical barrier around the cell while selectively facilitating material transport, signal transduction, and various other functions necessary for the cell viability. The complex functionality of cell membranes and the hierarchical motions and responses they exhibit demand a thorough understanding of the origin of different membrane dynamics and how they are influenced by molecular additives and environmental cues. These dynamic modes include single-molecule diffusion, thermal fluctuations, and large-scale membrane deformations, to name a few. This review highlights advances in investigating structure-driven dynamics associated with model cell membranes, with a particular focus on insights gained from neutron scattering and spectroscopy experiments. We discuss the uniqueness of neutron contrast variation and its remarkable potential in probing selective membrane structure and dynamics on spatial and temporal scales over which key biological functions occur. We also present a summary of current and future opportunities in synergistic combinations of neutron scattering with molecular dynamics (MD) simulations to gain further understanding of the molecular mechanisms underlying complex membrane functions.

A combined molecular/continuum-modeling approach to predict the small-angle neutron scattering of curved membranes

This paper develops a framework to compute the small-angle neutron scattering (SANS) from highly curved, dynamically fluctuating, and potentially inhomogeneous membranes. This method is needed to compute the scattering from nanometer-scale membrane domains that couple to curvature, as predicted by molecular modeling. The detailed neutron scattering length density of a small planar bilayer patch is readily available via molecular dynamics simulation. A mathematical, mechanical transformation of the planar scattering length density is developed to predict the scattering from curved bilayers. By simulating a fluctuating, curved, surface-continuum model, long time- and length-scales can be reached while, with the aid of the planar-to-curved transformation, the molecular features of the scattering length density can be retained. A test case for the method is developed by constructing a coarse-grained lipid vesicle following a protocol designed to relieve both the osmotic stress inside the vesicle and the lipid-number stress between the leaflets. A question was whether the hybrid model would be able to replicate the scattering from the highly deformed inner and outer leaflets of the small vesicle. Matching the scattering of the full (molecular vesicle) and hybrid (continuum vesicle) models indicated that the inner and outer leaflets of the full vesicle were expanded laterally, consistent with previous simulations of the Martini forcefield that showed thinning in small vesicles. The vesicle structure is inconsistent with a zero-tension leaflet deformed by a single set of elastic parameters, and the results show that this is evident in the scattering. The method can be applied to translate observations of any molecular model's neutron scattering length densities from small patches to large length and timescales.