Direct comparison of elastic incoherent neutron scattering experiments with molecular dynamics simulations of DMPC phase transitions

The dynamical Matryoshka model: 1. Incoherent neutron scattering functions for lipid dynamics in bilayers

Fluid lipid bilayers are the building blocks of biological membranes. Although there is a large amount of experimental data using incoherent quasi-elastic neutron scattering (QENS) techniques to study membranes, very little theoretical works have been developed to study the local dynamics of membranes. The main objective of this work is to build a theoretical framework to study and describe the local dynamics of lipids and derive analytical expressions of intermediate scattering functions (ISF) for QENS. As results, we developed the dynamical Matryoshka model which describes the local dynamics of lipid molecules in membrane layers as a nested hierarchical convolution of three motional processes: (i) individual motions described by the vibrational motions of H-atoms; (ii) internal motions including movements of the lipid backbone, head groups and tails, and (iii) molecule movements of the lipid molecule as a whole. The analytical expressions of the ISF associated with these movements are all derived. For use in analyzing the QENS experimental data, we also derived an analytical expression for the aggregate ISF of the Matryoshka model which involves an elastic term plus three inelastic terms of well-separated time scales and whose amplitudes and rates are functions of the lipid motions. And as an illustrative application, we used the aggregated ISF to analyze the experimental QENS data on a lipid sample of multilamellar bilayers of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine). It is clear from this analysis that the dynamical Matryoshka model describes very well the experimental data and allow extracting the dynamical parameters of the studied system.

The dynamical Matryoshka model: 3. Diffusive nature of the atomic motions contained in a new dynamical model for deciphering local lipid dynamics

In accompanying papers [Bicout et al., BioRxiv https://doi.org/10.1101/2021.09.21.461198 (2021); Cissé et al., BioRxiv https://doi.org/10.1101/2022.03.30.486370 (2022)], a new model called Matryoshka model has been proposed to describe the geometry of atomic motions in phospholipid molecules in bilayers and multilamellar vesicles based on their quasielastic neutron scattering (QENS) spectra. Here, in order to characterize the relaxational aspects of this model, the energy widths of the QENS spectra of the samples were analyzed first in a model-free way. The spectra were decomposed into three Lorentzian functions, which are classified as slow, intermediate, and fast motions depending on their widths. The analysis provides the diffusion coefficients, residence times, and geometrical parameters for the three classes of motions. The results corroborate the parameter values such as the amplitudes and the mobile fractions of atomic motions obtained by the application of the Matryoshka model to the same samples. Since the current analysis was carried out independently of the development of the Matryoshka model, the present results enhance the validity of the model. The model will serve as a powerful tool to decipher the dynamics of lipid molecules not only in model systems, but also in more complex systems such as mixtures of different kinds of lipids or natural cell membranes.

High temperature molecular motions within a model protomembrane architecture

Modern phospholipid membranes are known to be in a functional, physiological state, corresponding to the liquid crystalline phase, only under very precise external conditions. The phase is characterised by specific lipid motions, which seem mandatory to permit sufficient flexibility and stability for the membrane. It can be assumed that similar principles hold for proto-membranes at the origin of life although they were likely composed of simpler, single chain fatty acids and alcohols. In the present study we investigated molecular motions of four types of model membranes to shed light on the variations of dynamics and structure from low to high temperature as protocells might have existed close to hot vents. We find a clear hierarchy among the flexibilities of the samples, where some structural parameters seem to depend on the lipid type used while others do not.

Influence of hydration on segmental chain librations and dynamical transition in lipid bilayers

Continuous wave electron paramagnetic resonance spectroscopy of chain-labeled phospholipids is used to investigate the effects of hydration on the librational oscillations and the dynamical transition of phospholipid membranes in the low-temperature range 120-270 K. Bilayers of dipalmitoylphostatidiycholine (DPPC) spin-labeled at the first acyl chain segments and at the methyl ends and prepared at full, low, and very low hydration are considered. The segmental mean-square angular amplitudes of librations, 〈α²〉, are larger in the bilayer interior than at the polar/apolar interface and larger in the fully and low hydrated than in the very low hydrated membranes. For chain segments at the beginning of the hydrocarbon region, 〈α²〉-values are markedly restricted and temperature independent in DPPC with the lowest water content, whereas they increase with temperature in the low and fully hydrated bilayers, particularly at the highest temperatures. For chain segments at the chain termini, the librational amplitudes increase progressively, first slowly and then more rapidly with temperature in bilayers at any level of hydration. From the temperature dependence of the mean-square librational amplitude, the dynamical transition is detected around 240 K at the polar/apolar interface in fully and low hydrated DPPC and at around 225 K at the inner hydrocarbon region for bilayers at any hydration condition. At the dynamical transition the bilayers cross low energy barriers of activation energy in the range 10-20 kJ/mol. The results highlight biophysical properties of DPPC bilayers at low-temperature and provide evidence of the effects of the hydration on the dynamical transition in bilayers.

Astragaloside IV- and nesfatin-1-encapsulated phosphatidylserine liposomes conjugated with wheat germ agglutinin and leptin to activate anti-apoptotic pathway and block phosphorylated tau protein expression for Parkinson's disease treatment

Heap-up of α-synuclein (α-Syn) and its association with tau protein are esteemed to trigger the onset of Parkinson's disease (PD). The purpose of this study was to develop multi-functional liposomes incorporated with 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC), cholesterol, 1,2-dimyristoyl-sn-glycero-3-phosphocholine and phosphatidylserine (PS) to load astragaloside IV (AS-IV) and nestifin-1 (NF-1), followed by grafting with wheat germ agglutinin (WGA) and leptin (Lep) (WGA-Lep-AS-IV-NF-1-PS-liposomes) to protect dopaminergic neurons from apoptosis. Experimental results showed that increasing the mole percentage of DSPC and PS enhanced the particle size, particle stability and entrapment efficiency of AS-IV and NF-1, and reduced the drug releasing rate. Strong affinity of NF-1 to PS was evidenced by nuclear magnetic resonance spectroscopy. WGA-Lep-AS-IV-NF-1-PS-liposomes diminished transendothelial electrical resistance and improved the capacity of propidium iodide, AS-IV and NF-1 to penetrate the blood-brain barrier (BBB). Immunocytochemical staining exhibited the ability of functionalized liposomes to target Lep receptor and α-Syn in MPP⁺-insulted SH-SY5Y cells. Western blots revealed a substantial reduction of α-Syn and phosphorylated tau protein in the anti-oxidative pathway through interaction with PS. During the course of treatment with WGA-Lep-AS-IV-NF-1-PS-liposomes, the combined activity of AS-IV and NF-1 and recognition capability simultaneously decreased the expression of Bax, and increased the expressions of Bcl-2, tyrosine hydroxylase and dopamine transporter. The liposomes carrying AS-IV and NF-1 can rescue degenerated neurons and are a promising formulation to achieve better PD management.

Investigation of the Action of Peptides on Lipid Membranes

Calorimetric and incoherent neutron scattering methods were employed to investigate the action of magainin 2 and PGLa peptides on the phase behavior and molecular dynamics of lipids mimicking cytoplasmic membranes of Gram-negative bacteria. The impact of the peptides, tested individually and cooperatively by differential scanning calorimetry, presented a broadened peak, sometimes with a second shoulder, depicting the phase transition temperature around 21 °C. Neutron scattering revealed a small but significant variation of the membrane dynamics due to the peptides in both in-plane and out-of-plane directions. Although we did not find a clear hint for synergy in the interplay of the two peptides, the calorimetric and neutron data give compatible results in terms of a decrease of the enthalpy due to the presence of the peptides, which destabilize the membrane. The dynamics in the two directions was differentiated when the individual peptides were added to the membranes, but the impact was smaller when both peptides were added together.

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and X-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts (Bolmatov et al., 2020). Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager's regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

Differences between calcium rich and depleted alpha-lactalbumin investigated by molecular dynamics simulations and incoherent neutron scattering

We present a study comparing atomic motional amplitudes in calcium rich and depleted alpha-lactalbumin. The investigations were performed by elastic incoherent neutron scattering (EINS) and molecular dynamics (MD) simulations. As the variations were expected to be very small, three different hydration levels and timescales (instrumental resolutions) were measured. In addition, we used two models to extract the mean square displacements (MSDs) from the EINS data, one taking into account the motional heterogeneity of the MSD. At a timescale of several nanoseconds, small differences in the amplitudes between the calcium enriched and depleted alpha-lactalbumin are visible, whereas at lower timescales no changes can be concluded within the statistics. The results are compared to MD simulations at 280 and 300 K by extracting the MSDs of the trajectories in two separate ways: first by direct calculation, and second by a virtual neutron experiment using the same models as for the experimental data. We show that the simulated data give qualitatively similar results as the experimental data but quantitatively there are differences. Furthermore, the distribution of the MSDs in the simulations suggests that the inclusion of heterogeneity is reasonable for alpha-lactalbumin, but a bi-or trimodal approach may be sufficient.

Solution Structure and Conformational Flexibility in the Active State of the Orange Carotenoid Protein. Part II: Quasielastic Neutron Scattering

Orange carotenoid proteins (OCPs), which are protecting cyanobacterial light-harvesting antennae from photodamage, undergo a pronounced structural change upon light absorption. In addition, the active state is anticipated to boost a significantly higher molecular flexibility similar to a "molten globule" state. Here, we used quasielastic neutron scattering to directly characterize the vibrational and conformational molecular dynamics of OCP in its ground and active states, respectively, on the picosecond time scale. At a temperature of 100 K, we observe mainly (vibronic) inelastic features with peak energies at 5 and 6 meV (40 and 48 cm^-1, respectively). At physiological temperatures, however, two (Lorentzian) quasielastic components represent localized protein motions, that is, stochastic structural fluctuations of protein side chains between various conformational substates of the protein. Global diffusion of OCP is not observed on the given time scale. The slower Lorentzian component is affected by illumination and can be well-characterized by a jump-diffusion model. While the jump diffusion constant D is (2.82 ± 0.01) × 10^-5 cm²/s at 300 K in the ground state, it is increased by ∼20% to (3.48 ± 0.01) × 10^-5 cm²/s in the active state, revealing a strong enhancement of molecular mobility. The increased mobility is also reflected in the average atomic mean square displacement ⟨u²⟩; we determine a ⟨u²⟩ of 1.47 ± 0.05 Å in the ground state, but 1.86 ± 0.05 Å in the active state (at 300 K). This effect is assigned to two factors: (i) the elongated structure of the active state with two widely separated protein domains is characterized by a larger number of surface residues with a concomitantly higher degree of motional freedom and (ii) a larger number of hydration water molecules bound at the surface of the protein. We thus conclude that the active state of the orange carotenoid protein displays an enhanced conformational dynamics. The higher degree of flexibility may provide additional channels for nonradiative decay so that harmful excess energy can be more efficiently converted to heat.

The effect of pH on the dynamics of natural membranes

Pure phospholipids and membrane fragments from bacterial cells living under various conditions were studied against the influence of the surrounding acidity on the internal dynamics. For that we compared mean square displacements extracted from elastic incoherent neutron scattering data, measured both at low and at neutral pH, of the phospholipids 1,2-dimyristoyl-sn-glycero-3-phosphocholine and of samples from neutralophilic and acidophilic micro-organisms (some being hyperthermophilic and others mesophilic). The lipids showed a slight shift in the phase transition temperature of about 4 degrees under pH variation and became slightly more mobile at lower pH. The membrane fragments not used to extreme acidic conditions were significantly more sensitive to variations in the pH values, whereas the acidophilic and -tolerant samples were much less influenced by this parameter. They presented the higher softness at low pH, which was closer to their native condition. Such findings might be a hint for adaptation mechanisms to different acidity conditions.

Thermodynamics of lipid multi-lamellar vesicles in presence of sterols at high hydrostatic pressure

We compared the effect of cholesterol at different concentration on the phase behaviour of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) multilamellar vesicles. We used pressure perturbation differential scanning calorimetry (PPC) that studies a system on the whole by giving access to relevant thermodynamic quantities, and elastic incoherent neutron scattering (EINS) that probes local motions of a system at the atomic level by allowing extraction of dynamical parameters. PPC revealed that the volume expansion coefficient of DMPC and DMPC/Cholesterol samples with 13 and 25 mol% cholesterol is a linear function of the heat capacity measured by differential scanning calorimetry. Neutron backscattering spectroscopy showed that the mean square displacements of H atoms do exhibit an increase with temperature and a decrease under increasing pressure. Cholesterol added at concentrations of 25 and 50 mol% led to suppression of the main phase transition. Taking advantage of these results, the present study aims (i) to show that calorimetry and EINS using the Bicout and Zaccai model equally permit to get access to thermodynamic quantities characterizing pure DMPC and DMPC/cholesterol mixtures, thus directly confirming the theoretical method, and (ii) to validate our approach as function of temperature and of pressure, as both are equally important and complementary thermodynamic variables.

The Dynamical Transition of Lipid Multilamellar Bilayers as a Matter of Cooperativity

The present study is the application of a two-state model formerly developed by Bicout and Zaccai [ Bicout , D. J. and Zaccai , G. Biophys. J. 2001 , 80 ( 3 ), 1115 - 1123 ] to describe the dynamical transition exhibited in the atomic mean square displacements of biological samples in terms of dynamic and thermodynamic parameters. Data were obtained by elastic incoherent neutron scattering on 1,2-dimyristoyl-sn-glycero-3-phosphocholine lipid membranes in various hydration states and on one partially per-deuterated lipid membrane. Fitting the data with the model allowed investigating which parts of lipid molecules were mainly involved in the dynamical transition, heads, tails, or both. Clear differences were found between the fully protonated and partially deuterated membranes. These findings shed light on the question of what is the degree of dynamical cooperativity of the atoms during the transition. Whereas the level of hydration does not significantly affect it, as the dry, the intermediate dry, and fully hydrated membranes all undergo a rather broad transition, the transition of the lipid tails is much sharper and sets in at much lower temperature than that of the heads. Therefore, the dynamical cooperativity appears high among the particles in the tails. Moreover, the transition of the lipid tails has to be completed first before the one of the head groups starts.