Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

Osmotic stress studies of G-protein-coupled receptor rhodopsin activation

We summarize and critically review osmotic stress studies of the G-protein-coupled receptor rhodopsin. Although small amounts of structural water are present in these receptors, the effect of bulk water on their function remains uncertain. Studies of the influences of osmotic stress on the GPCR archetype rhodopsin have given insights into the functional role of water in receptor activation. Experimental work has discovered that osmolytes shift the metarhodopsin equilibrium after photoactivation, either to the active or inactive conformations according to their molar mass. At least 80 water molecules are found to enter rhodopsin in the transition to the photoreceptor active state. We infer that this movement of water is both necessary and sufficient for receptor activation. If the water influx is prevented, e.g., by large polymer osmolytes or by dehydration, then the receptor functional transition is back shifted. These findings imply a new paradigm in which rhodopsin becomes solvent swollen in the activation mechanism. Water thus acts as an allosteric modulator of function for rhodopsin-like receptors in lipid membranes.

A review on the measurement of the bending rigidity of lipid membranes

This review provides an overview of the latest developments in both experimental and simulation techniques used to assess the bending rigidity of lipid membranes. It places special emphasis on experimental methods that utilize model vesicles to manipulate lipid compositions and other experimental parameters to determine the bending rigidity of the membrane. It also describes two commonly used simulation methods for estimating bending rigidity. The impact of various factors on membrane bending rigidity is summarized, including cholesterol, lipids, salt concentration, surface charge, membrane phase state, peptides, proteins, and polyethylene glycol. These factors are shown to influence the bending rigidity, contributing to a better understanding of the biophysical properties of membranes and their role in biological processes. Furthermore, the review discusses future directions and potential advancements in this research field, highlighting areas where further investigation is required.

Location of dopamine in lipid bilayers and its relevance to neuromodulator function

Dopamine (DA) is a neurotransmitter that also acts as a neuromodulator, with both functions being essential to brain function. Here, we present the first experimental measurement of DA location in lipid bilayers using x-ray diffuse scattering, solid-state deuterium NMR, and electron paramagnetic resonance. We find that the association of DA with lipid headgroups as seen in electron density profiles leads to an increase of intermembrane repulsion most likely due to electrostatic charging. DA location in the lipid headgroup region also leads to an increase of the cross-sectional area per lipid without affecting the bending rigidity significantly. The order parameters measured by solid-state deuterium NMR decrease in the presence of DA for the acyl chains of PC and PS lipids, consistent with an increase in the area per lipid due to DA. Most importantly, these results support the hypothesis that three-dimensional diffusion of DA to target membranes could be followed by relatively more efficient two-dimensional diffusion to receptors within those membranes.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Phospholipid headgroups govern area per lipid and emergent elastic properties of bilayers

Phospholipid bilayers are liquid-crystalline materials whose intermolecular interactions at mesoscopic length scales have key roles in the emergence of membrane physical properties. Here we investigated the combined effects of phospholipid polar headgroups and acyl chains on biophysical functions of membranes with solid-state 2H NMR spectroscopy. We compared the structural and dynamic properties of phosphatidylethanolamine and phosphatidylcholine with perdeuterated acyl chains in the solid-ordered (so) and liquid-disordered (ld) phases. Our analysis of spectral lineshapes of 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE-d62) and 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphocholine (DPPC-d62) in the so (gel) phase indicated an all-trans rotating chain structure for both lipids. Greater segmental order parameters (SCD) were observed in the ld (liquid-crystalline) phase for DPPE-d62 than for DPPC-d62 membranes, while their mixtures had intermediate values irrespective of the deuterated lipid type. Our results suggest the SCD profiles of the acyl chains are governed by methylation of the headgroups and are averaged over the entire system. Variations in the acyl chain molecular dynamics were further investigated by spin-lattice (R1Z) and quadrupolar-order relaxation (R1Q) measurements. The two acyl-perdeuterated lipids showed distinct differences in relaxation behavior as a function of the order parameter. The R1Z rates had a square-law dependence on SCD, implying collective mesoscopic dynamics, with a higher bending rigidity for DPPE-d62 than for DPPC-d62 lipids. Remodeling of lipid average and dynamic properties by methylation of the headgroups thus provides a mechanism to control the actions of peptides and proteins in biomembranes.

Cholesterol Stiffening of Lipid Membranes

Biomembrane order, dynamics, and other essential physicochemical parameters are controlled by cholesterol, a major component of mammalian cell membranes. Although cholesterol is well known to exhibit a condensing effect on fluid lipid membranes, the extent of stiffening that occurs with different degrees of lipid acyl chain unsaturation remains an enigma. In this review, we show that cholesterol locally increases the bending rigidity of both unsaturated and saturated lipid membranes, suggesting there may be a length-scale dependence of the bending modulus. We review our published data that address the origin of the mechanical effects of cholesterol on unsaturated and polyunsaturated lipid membranes and their role in biomembrane functions. Through a combination of solid-state deuterium NMR spectroscopy and neutron spin-echo spectroscopy, we show that changes in molecular packing cause the universal effects of cholesterol on the membrane bending rigidity. Our findings have broad implications for the role of cholesterol in lipid–protein interactions as well as raft-like mixtures, drug delivery applications, and the effects of antimicrobial peptides on lipid membranes.

Hydration-mediated G-protein-coupled receptor activation

Although G-protein–coupled receptors (GPCRs) control vast physiological pathways, their activation remains chemically and physically enigmatic. Our osmotic stress studies of the visual receptor rhodopsin have redefined the standard model of GPCR signaling by revealing the essential role of bulk water. We show results consistent with a large number of water molecules flooding the rhodopsin interior during activation to stabilize the effector binding conformation. These results suggest a model of GPCR activation in which the receptor becomes solvent-swollen upon formation of the active state. We thus demonstrate the mechanism whereby water acts as a powerful allosteric modulator of a pharmacologically important membrane protein family.

The Rhodopsin family of G-protein–coupled receptors (GPCRs) comprises the targets of nearly a third of all pharmaceuticals. Despite structural water present in GPCR X-ray structures, the physiological relevance of these solvent molecules to rhodopsin signaling remains unknown. Here, we show experimental results consistent with the idea that rhodopsin activation in lipid membranes is coupled to bulk water movements into the protein. To quantify hydration changes, we measured reversible shifting of the metarhodopsin equilibrium due to osmotic stress using an extensive series of polyethylene glycol (PEG) osmolytes. We discovered clear evidence that light activation entails a large influx of bulk water (∼80–100 molecules) into the protein, giving insight into GPCR activation mechanisms. Various size polymer osmolytes directly control rhodopsin activation, in which large solutes are excluded from rhodopsin and dehydrate the protein, favoring the inactive state. In contrast, small osmolytes initially forward shift the activation equilibrium until a quantifiable saturation point is reached, similar to gain-of-function protein mutations. For the limit of increasing osmolyte size, a universal response of rhodopsin to osmotic stress is observed, suggesting it adopts a dynamic, hydrated sponge-like state upon photoactivation. Our results demand a rethinking of the role of water dynamics in modulating various intermediates in the GPCR energy landscape. We propose that besides bound water, an influx of bulk water plays a necessary role in establishing the active GPCR conformation that mediates signaling.

Protein resonance assignment by solid-state NMR based on 1H-detected 13C double-quantum spectroscopy at fast MAS

Solid-state NMR spectroscopy is a powerful technique to study insoluble and non-crystalline proteins and protein complexes at atomic resolution. The development of proton (¹H) detection at fast magic-angle spinning (MAS) has considerably increased the analytical capabilities of the technique, enabling the acquisition of ¹H-detected fingerprint experiments in few hours. Here an approach based on double-quantum (DQ) ¹³C spectroscopy, detected on ¹H, is proposed for fast MAS regime (> 60 kHz) to perform the sequential assignment of insoluble proteins of small size, without any specific deuteration requirement. By combining two three-dimensional ¹H detected experiments correlating a ¹³C DQ dimension respectively to its intra-residue and sequential ¹⁵ N-¹H pairs, a sequential walk through DQ (Ca + CO) resonance is obtained. The approach takes advantage of fast MAS to achieve an efficient sensitivity and the addition of a DQ dimension provides spectral features useful for the resonance assignment process.

Activation of the G-Protein-Coupled Receptor Rhodopsin by Water

Visual rhodopsin is an important archetype for G-protein-coupled receptors, which are membrane proteins implicated in cellular signal transduction. Herein, we show experimentally that approximately 80 water molecules flood rhodopsin upon light absorption to form a solvent-swollen active state. An influx of mobile water is necessary for activating the photoreceptor, and this finding is supported by molecular dynamics (MD) simulations. Combined force-based measurements involving osmotic and hydrostatic pressure indicate the expansion occurs by changes in cavity volumes, together with greater hydration in the active metarhodopsin-II state. Moreover, we discovered that binding and release of the C-terminal helix of transducin is coupled to hydration changes as may occur in visual signal amplification. Hydration-dehydration explains signaling by a dynamic allosteric mechanism, in which the soft membrane matter (lipids and water) has a pivotal role in the catalytic G-protein cycle.