Rotational decoupling between the hydrophilic and hydrophobic regions in lipid membranes.

Cells use homeostatic mechanisms to ensure an optimal composition of distinct types of lipids in cellular membranes. The hydrophilic region of biological lipid membranes is mainly composed of several types of phospholipid headgroups that interact with incoming molecules, nanoparticles, and viruses, whereas the hydrophobic region consists of a distribution of acyl chains and sterols affecting membrane fluidity/rigidity related properties and forming an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of phospholipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well-known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm in which phospholipid dipole headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.

Filling the Gap with Long n-Alkanes: Incorporation of C20 and C30 into Phospholipid Membranes.

Investigating how hydrophobic molecules mix with phospholipid bilayers and how they affect membrane properties is commonplace in biophysics. Despite this, a molecular-level empirical description of a membrane model as simple as a phospholipid bilayer with long linear hydrophobic chains incorporated is still missing. Here, we present an unprecedented molecular characterization of the incorporation of two long n-alkanes, n-eicosane (C20) and n-triacontane (C30) with 20 and 30 carbons, respectively, in phosphatidylcholine (PC) bilayers using a combination of experimental techniques (2H NMR, 31P NMR, 1H-13C dipolar recoupling solid-state NMR, X-ray scattering, and cryogenic electron microscopy) and atomistic molecular dynamics (MD) simulations. At low hydration, deuterated C20 and C30 yield 2H NMR spectra evidencing anisotropic-motion, which demonstrates their miscibility in PC membranes up to a critical alkane-to-acyl-chain volume fraction, ϕc. The acquired 2H NMR spectra of C20 and C30 have notably different lineshapes. At low alkane volume fractions below ϕc, CHARMM36 MD simulations predict such 2H NMR spectra qualitatively and thus enable an atomistic-level interpretation of the spectra. Above ϕc, the 2H NMR lineshapes become characteristic of motions in the intermediate-regime that, together with the MD simulation results, suggest the onset of immiscibility between the alkane molecules and the acyl chains. For all the systems investigated, the phospholipid molecular structure is unperturbed by the presence of the alkanes. However, at conditions of excess hydration and at surprisingly low alkane fractions below ϕc, a peak characteristic of isotropic motion is observed in both the 2H spectra of the alkanes and 31P spectra of the phospholipids, strongly indicating that the incorporation of the alkanes induces a reduction on the average radius of the lipid vesicles.

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

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.

Time-domain proton-detected local-field NMR for molecular structure determination in complex lipid membranes.

Proton-detected local-field (PDLF) NMR spectroscopy, using magic-angle spinning and dipolar recoupling, is presently the most powerful experimental technique for obtaining atomistic structural information from small molecules undergoing anisotropic motion. Common examples include peptides, drugs, or lipids in model membranes and molecules that form liquid crystals. The measurements on complex systems are however compromised by the larger number of transients required. Retaining sufficient spectral quality in the direct dimension requires that the indirect time-domain modulation becomes too short for yielding dipolar splittings in the frequency domain. In such cases, the dipolar couplings can be obtained by fitting the experimental data; however ideal models often fail to fit PDLF data properly due to effects of radiofrequency field (RF) spatial inhomogeneity. Here, we demonstrate that by accounting for RF spatial inhomogeneity in the modeling of R-symmetry-based PDLF NMR experiments, the fitting accuracy is improved, facilitating the analysis of the experimental data. In comparison to the analysis of dipolar splittings without any fitting procedure, the accurate modeling of PDLF measurements makes possible three important improvements: the use of shorter experiments that enable the investigation of samples with a higher level of complexity, the measurement of C-H bond order parameters with smaller magnitudes |SCH| and of smaller variations of |SCH| caused by perturbations of the system, and the determination of |SCH| values with small differences from distinct sites having the same chemical shift. The increase in fitting accuracy is demonstrated by comparison with 2H NMR quadrupolar echo experiments on mixtures of deuterated and non-deuterated dimyristoylphosphatidylcholine (DMPC) and with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) membranes. Accurate modeling of PDLF NMR experiments is highly useful for investigating complex membrane systems. This is exemplified by application of the proposed fitting procedure for the characterization of membranes composed of a brain lipid extract with many distinct lipid types.

An atlas for human brain myelin content throughout the adult life span.

Myelin water imaging is a quantitative neuroimaging technique that provides the myelin water fraction (MWF), a metric highly specific to myelin content, and the intra-/extra-cellular T2 (IET2), which is related to water and iron content. We coupled high-resolution data from 100 adults with gold-standard methodology to create an optimized anatomical brain template and accompanying MWF and IET2 atlases. We then used the MWF atlas to characterize how myelin content relates to demographic factors. In most brain regions, myelin content followed a quadratic pattern of increase during the third decade of life, plateau at a maximum around the fifth decade, then decrease during later decades. The ranking of mean myelin content between brain regions remained consistent across age groups. These openly available normative atlases can facilitate evaluation of myelin imaging results on an individual basis and elucidate the distribution of myelin content between brain regions and in the context of aging.