Hydration dependent studies of highly aligned multilayer lipid membranes by neutron scattering

Quasi-elastic neutron scattering reveals the relationship between the dynamical behavior of phospholipid headgroups and hydration water

The dynamics of hydration water (HW) in 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine (DMPE) was investigated by means of quasi-elastic neutron scattering (QENS) and compared with those observed in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). The headgroup dynamics of DMPE was investigated using a mixture of tail-deuterated DMPE and D2O, and the QENS profiles were interpreted as consisting of three modes. The fast mode comprised the rotation of hydrogen atoms in -NH3+ and -CH2- groups in the headgroup of DMPE, the medium-speed mode comprised fluctuations in the entire DMPE molecule, and the slow mode comprised fluctuations in the membrane. These interpretations were confirmed using molecular dynamics (MD) simulations. The HW dynamics analysis was performed on a tail-deuterated DMPE and H2O mixture. The QENS profiles were analyzed in terms of three modes: (1) a slow mode, identified as loosely bound HW in the DMPC membrane; (2) a medium-speed mode similar to free HW in the DMPC membrane; and (3) a fast mode, identified as rotational motion. The relaxation time for the fast mode was approximately six times shorter than that of rotational water in DMPC, consistent with the results of terahertz time-domain spectroscopy. The activation energy of medium-speed HW in DMPE differed from that of free HW in DMPC, suggesting the presence of different hydration states or hydrogen-bonded networks around the phosphocholine and phosphoethanolamine headgroups.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

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: 2. Modeling of local lipid dynamics at the sub-nanosecond timescale in phospholipid membranes

Biological membranes are generally formed by lipids and proteins. Often, the membrane properties are studied through model membranes formed by phospholipids only. They are molecules composed by a hydrophilic head group and hydrophobic tails, which can present a panoply of various motions, including small localized movements of a few atoms up to the diffusion of the whole lipid or collective motions of many of them. In the past, efforts were made to measure these motions experimentally by incoherent neutron scattering and to quantify them, but with upcoming modern neutron sources and instruments, such models can now be improved. In the present work, we expose a quantitative and exhaustive study of lipid dynamics on DMPC and DMPG membranes, using the Matryoshka model recently developed by our group. The model is confronted here to experimental data collected on two different membrane samples, at three temperatures and two instruments. Despite such complexity, the model describes reliably the data and permits to extract a series of parameters. The results compare also very well to other values found in the literature.

Multiscale lipid membrane dynamics as revealed by neutron spectroscopy

The plasma membrane is one of the principal structural components of the cell and, therefore, one of the key components of the cellular life. Because the membrane's dynamics links the membrane's structure and function, the complexity and the broad range of the membrane's motions are essential for the enormously diverse functionality of the cell membrane. Even for the main membrane component, the lipid bilayer, considered alone, the range and complexity of the lipid motions are remarkable. Spanning the time scale from sub-picosecond to minutes and hours, the lipid motion in a bilayer is challenging to study even when a broad array of dynamic measurement techniques is employed. Neutron scattering plays a special role among such dynamic measurement techniques, particularly, because it involves the energy transfers commensurate with the typical intra- and inter- molecular dynamics and the momentum transfers commensurate with intra- and inter-molecular distances. Thus, using neutron scattering-based techniques, the spatial and temporal information on the lipid motion can be obtained and analysed simultaneously. Protium vs. deuterium sensitivity and non-destructive character of the neutron probe add to the remarkable prowess of neutron scattering for elucidating the lipid dynamics. Herein we present an overview of the neutron scattering-based studies of lipid dynamics in model membranes, with a discussion of the direct relevance and implications to the real-life cell membranes. The latter are much more complex systems than simple model membranes, consisting of heterogeneous non-stationary domains composed of lipids, proteins, and other small molecules, such as carbohydrates. Yet many fundamental aspects of the membrane behavior and membrane interactions with other molecules can be understood from neutron scattering measurements of the model membranes. For example, such studies can provide a great deal of information on the interactions of antimicrobial compounds with the lipid matrix of a pathogen membrane, or the interactions of drug molecules with the plasma membrane. Finally, we briefly discuss the recently emerging field of neutron scattering membrane studies with a reach far beyond the model membrane systems.

Evaluation of sample cell materials for aqueous solutions used in quasi-elastic neutron scattering measurements

For quasi-elastic neutron scattering (QENS) studies, sample cells made of pure or alloyed aluminium are frequently employed. Although the Al surface is protected by a passivating film, this film is not robust. Therefore, when the sample is an aqueous solution, chemical interactions between the Al surface and sample, promoted by corrosive entities such as chloride ions and certain conditions of pH, can compromise the integrity of the cell and interfere with the experiment. In this study, the corrosion susceptibilities of Al and its alloys were investigated by subjecting them to various treatments; the results were compared with those of other candidate materials with low chemical reactivity. This work showed that alloys with higher Al content and boehmite-coated surfaces are resistant to corrosion. In particular, for Al, the resistance is due to a reduction in the contact area achieved by reducing the surface roughness. QENS measurements of empty sample cells made of these materials revealed two results: (1) the profile of the cell fabricated with a copper-free Al alloy showed a minor dependence on the scattering vector magnitude Q and (2) reducing the real surface area of Al effectively suppresses its scattering intensity, while boehmite coating strengthens the scattering. Cells fabricated with Mo, Nb and single-crystal sapphire can be used as alternatives to Al because of their low scattering intensity and reduced dependence on Q.

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

Redefining the concept of hydration water near soft interfaces

Water determines the properties of biological systems. Therefore, understanding the nature of the mutual interaction between water and biosystems is of primary importance for a proper assessment of any biological activity, e.g., the efficacy of new drugs or vaccines. A convenient way to characterize the interactions between biosystems and water is to analyze their impact on water density and dynamics in the proximity of the interfaces. It is commonly accepted that water bulk density and dynamical properties are recovered at distances of the order of 1 nm away from the surface of biological systems. This notion leads to the definition of hydration or biological water as the nanoscopic layer of water covering the surface of biosystems and to the expectation that all the effects of the water-interface interaction are limited to this thin region. Here, we review some of our latest contributions, showing that phospholipid membranes affect the water dynamics, structural properties, and hydrogen bond network at a distance that is more than twice as large as the commonly evoked ∼1nm thick layer and of the order of 2.4 nm. Furthermore, we unveil that at a shorter distance ∼0.5nm from the membrane, instead, there is an additional interface between lipid-bound and unbound water. Bound water has a structural role in the stability of the membrane. Our results imply that the concept of hydration water should be revised or extended and pave the way to a deeper understanding of the mutual interactions between water and biological systems.

Understanding the Structure and Dynamics of Complex Biomembrane Interactions by Neutron Scattering Techniques

The membrane is one of the key structural materials of biology at the cellular level. Composed predominantly of a bilayer of lipids with embedded and bound proteins, it defines the boundaries of the cell and many organelles essential to life and therefore is involved in almost all biological processes. Membrane-specific interactions, such as drug binding to a membrane receptor or the interactions of an antimicrobial compound with the lipid matrix of a pathogen membrane, are of interest across the scientific disciplines. Herein we present a review, aimed at nonexperts, of the major neutron scattering techniques used in membrane studies: small-angle neutron scattering, neutron membrane diffraction, neutron reflectometry, quasielastic neutron scattering, and neutron spin echo. Neutron scattering techniques are well suited to studying biological membranes. The nondestructive nature of cold neutrons means that samples can be measured for long periods without fear of beam damage from ultraviolet, electron, or X-ray radiation, and neutron beams are highly penetrating, thus offering flexibility in samples and sample environments. Most important is the strong difference in neutron scattering lengths between the two most abundant forms of hydrogen, protium and deuterium. Changing the relative amounts of protium/deuterium in a sample allows the production of a series of neutron scattering data sets, enabling the observation of differing components within complex membrane architectures. This approach can be as simple as using the naturally occurring neutron contrast between different biomolecules to study components in a complex by changing the solution H₂O/D₂O ratio or as complex as selectively labeling individual components with hydrogen isotopes. This review presents an overview of each experimental technique with the neutron instrument configuration, related sample preparation and sample environment, and data analysis, highlighted by a special emphasis on using prominent neutron contrast to understand structure and dynamics. This review gives researchers a practical introduction to the often enigmatic suite of neutron beamlines, thereby lowering the barrier to taking advantage of these large-facility techniques to achieve new understandings of membranes and their interactions with other molecules.

Network Topology in Water Nanoconfined between Phospholipid Membranes

Water provides the driving force for the assembly and stability of many cellular components. Despite its impact on biological functions, a nanoscale understanding of the relationship between its structure and dynamics under soft confinement has remained elusive. As expected, water in contact with biological membranes recovers its bulk density and dynamics at ∼1 nm from phospholipid headgroups but surprisingly enhances its intermediate range order (IRO) over a distance, at least, twice as large. Here, we explore how the IRO is related to the water's hydrogen-bond network (HBN) and its coordination defects. We characterize the increased IRO by an alteration of the HBN up to more than eight coordination shells of hydration water. The HBN analysis emphasizes the existence of a bound-unbound water interface at ∼0.8 nm from the membrane. The unbound water has a distribution of defects intermediate between bound and bulk water, but with density and dynamics similar to bulk, while bound water has reduced thermal energy and many more HBN defects than low-temperature water. This observation could be fundamental for developing nanoscale models of biological interactions and for understanding how alteration of the water structure and topology, for example, due to changes in extracellular ions concentration, could affect diseases and signaling. More generally, it gives us a different perspective to study nanoconfined water.

Dynamic coupling of a hydration layer to a fluid phospholipid membrane: intermittency and multiple time-scale relaxations

Understanding the coupling of a hydration layer and a lipid membrane is crucial to gaining access to membrane dynamics and understanding its functionality towards various biological processes.

Quantification of spatio-temporal scales of dynamical heterogeneity of water near lipid membranes above supercooling

A hydrated 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid membrane is investigated using an all atom molecular dynamics simulation at 308 K to determine the physical sources of universal slow relaxations of hydration layers and length-scale of the spatially heterogeneous dynamics. Continuously residing interface water (IW) molecules hydrogen bonded to different moieties of lipid heads in the membrane are identified. The non-Gaussian parameters of all classes of IW molecules show a cross-over from cage vibration to translational diffusion. A significant non-Gaussianity is observed for the IW molecules exhibiting large length correlations in translational van Hove functions. Two time-scales for the ballistic motions and hopping transitions are obtained from the self intermediate scattering functions of the IW molecules with an additional long relaxation, which disappears for bulk water. The long relaxation time-scales for the IW molecules obtained from the self intermediate scattering functions are in good accordance with the hydrogen bond relaxation time-scales irrespective of the nature of the chemical confinement and the confinement lifetime. Employing a block analysis approach, the length-scale of dynamical heterogeneities is captured from a transition from non-Gaussianity to Gaussianity in van Hove correlation functions of the IW molecules. The heterogeneity length-scale is comparable to the wave-length of the small and weak undulations of the membrane calculated by Fourier transforms of lipid tilts. This opens up a new avenue towards a possible correlation between heterogeneity length-scale and membrane curvature more significant for rippled membranes. Thus, our analyses provide a measure towards the spatio-temporal scale of dynamical heterogeneity of confined water near membranes.

Broadband dielectric spectroscopy from sub GHz to THz of hydrated lipid bilayer of DMPC

In order to study the dynamics of a phospholipid and its hydration water, we measured complex dielectric spectra of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) from the sub-GHz to the THz frequency region with varying the temperature and hydration level of the sample. Spectra obtained from a vector network analyzer and two terahertz time-domain spectrometers are adjusted, which enables us to analyze the dielectric spectra from the sub-GHz region to the THz region by a model function. We confirmed a fast relaxational mode in the sub-THz region, which was suggested by the previous work which only used a THz spectrometer.

Elastic compliance as a tool to understand Hofmeister ion specific effect in DMPC liposomes

Elastic compliance of DMPC liposomes with Hofmeister electrolytes: NaCl, Na₂SO₄, Na₂CO₃, NaNO₃, KCl and MgCl₂ studied using Quartz crystal microbalance with dissipation has been correlated with changes in their lamellar spacing from SAXS. The study suggests that hydration water of the different ions has an effect on the overall packing of the lipid bilayer that results as either a dehydrated liposome or where water smears the surface of the liposomes. Ratio of hydrogen bonded carbonyl and phosphate of polar region of the liposomes from ATR-FTIR spectroscopy, suggests that the polar groups are less hydrated due to the displacement of water by the electrolytes compared to pure DMPC and ordered in the sequence for cations as: K⁺ < Na⁺,Mg²⁺ and for anions as SO₄^2- < CO₃^2- < Cl^- < NO₃^-. These findings show the usefulness of Elastic compliance for structural studies of composite phospholipid bilayers, lipid-protein complexes and lipid systems of reduced dimensionalities.

Relaxation dynamics of saturated and unsaturated oriented lipid bilayers

We present experimental measurements and analysis of the dynamics and the phase behaviour of saturated DMPC and unsaturated DOPC oriented multi-lamellar bilayers. Elastic and inelastic neutron scattering were used to directly probe the dynamical processes of these membrane systems on time and length scales relevant to the internal and localized motion of lipid monomers. Mobility in this regime can be informative in elucidating the local interactions responsible for material properties of these fluid lipid systems. DMPC and DOPC are structurally similar in terms of their membrane hydrophobic thickness; however, they exhibit different mechanical properties in terms of both elastic compressibility and bending moduli. The analyses suggest that the constraint imposed by the double bonds in DOPC acyl chains restricts atomic motion in both liquid and gel phases compared to DMPC. We discuss applications of molecular dynamics to further elucidate the atomic details of the dynamical processes. Such an understanding may suggest how membrane properties can be tuned using a variety of different lipid species.

Deciphering interactions of ionic liquids with biomembrane

Ionic liquids (ILs) are a special class of low-temperature (typically < 100 °C) molten salts, which have huge upsurge interest in the field of chemical synthesis, catalysis, electrochemistry, pharmacology, and biotechnology, mainly due to their highly tunable nature and exceptional properties. However, practical uses of ILs are restricted mainly due to their adverse actions on organisms. Understanding interactions of ILs with biomembrane is prerequisite to assimilate the actions of these ionic compounds on the organism. Here, we review different biophysical methods to characterize interactions between ILs and phospholipid membrane, a model biomembrane. All these studies indicate that ILs interact profoundly with the lipid bilayer and modulate the structure, microscopic dynamics, and phase behavior of the membrane, which could be the fundamental cause of the observed toxicity of ILs. Effects of ILs on the membrane are found to be strongly dependent on the lipophilicity of the IL and are found to increase with the alkyl chain length of IL. This can be correlated with the observed higher toxicity of IL with the longer alkyl chain length. These informations would be useful to tune the toxicity of IL which is required in designing environment-friendly nontoxic solvents of the so-called green chemistry for various practical applications.

Effects of ionic liquids on the nanoscopic dynamics and phase behaviour of a phosphatidylcholine membrane

Ionic liquids (ILs) are potential candidates for new antimicrobials due to their tunable antibacterial and antifungal properties that are required to keep pace with the growing challenge of bacterial resistance. To a great extent their antimicrobial actions are related to the interactions of ILs with cell membranes. Here, we report the effects of ILs on the nanoscopic dynamics and phase behaviour of a dimyristoylphosphatidylcholine (DMPC) membrane, a model cell membrane, as studied using neutron scattering techniques. Two prototypical imidazolium-based ILs 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM[BF4]) and 1-decyl-3-methylimidazolium tetrafluoroborate (DMIM[BF4]), which differ only in terms of the alkyl chain length of cations, have been used for the present study. Fixed Elastic Window Scan (FEWS) shows that the incorporation of ILs affects the phase behaviour of the phospholipid membrane significantly and the transition from a solid gel to a fluid phase shifts to lower temperature. This is found to be consistent with our differential scanning calorimetry measurements. DMIM[BF4], which has a longer alkyl chain cation, affects the phase behaviour more strongly in comparison to BMIM[BF4]. The pressure-area isotherms of the DMPC monolayer measured at the air-water interface show that in the presence of ILs, isotherms shift towards higher area-per lipid molecule. DMIM[BF4] is found to shift the isotherm to a greater extent compared to BMIM[BF4]. Quasielastic neutron scattering (QENS) data show that both ILs act as a plasticizer, which enhances the fluidity of the membrane. DMIM[BF4] is found to be a stronger plasticizing agent in comparison to BMIM[BF4] that has a cation with a shorter alkyl chain. The incorporation of DMIM[BF4] enhances not only the long range lateral motion but also the localised internal motion of the lipids. On the other hand, BMIM[BF4] acts weakly in comparison to DMIM[BF4] and mainly alters the localised internal motion of the lipids. Any subtle change in the dynamical properties of the membrane can profoundly affect the stability of the cell. Hence, the dominant effect of the IL with the longer chain length on the dynamics of the phospholipid membrane might be correlated with its cytotoxic activity. QENS data analysis has provided a quantitative description of the effects of the two imidazolium-based ILs on the dynamical and phase behaviour of the model cell membrane, which is essential for a detailed understanding of their action mechanism.

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.

Dynamical Behavior of Hydration Water Molecules between Phospholipid Membranes

The dynamical behavior of hydration water sandwiched between 1,2-dimyristyl-sn-glycero-3-phosphocholine (DMPC) bilayers was investigated by quasi-elastic neutron scattering (QENS) in the range between 275 and 316 K, where the main transition temperature of DMPC is interposed. The results revealed that the hydration water could be categorized into three types of water: (1) free water, whose dynamical behavior is slightly different from that of bulk water; (2) loosely bound water, whose dynamical behavior is 1 order of magnitude slower than that of the free water; and (3) tightly bound water, whose dynamical behavior is comparable with that of DMPC molecules. The number of loosely bound and tightly bound water molecules per DMPC molecule monotonically decreased and increased with decreasing temperature, respectively, and the sum of these water molecules remained constant. The number of free water molecules per DMPC molecule was constant in the measured temperature range.

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.

Inelastic and quasi-elastic neutron scattering spectrometers in J-PARC

J-PARC, Japan Proton Accelerator Research Complex provides short pulse proton beam at a repetition rate 25Hz and the maximum power is expected to be 1MW. Materials and Life Science Experimental Facility (MLF) has 23 neutron beam ports and 21 instruments have already been operated or under construction/commissioning. There are 6 inelastic/quasi-elastic neutron scattering spectrometers and the complementary use of these spectrometers will open new insight for life science. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Structural Interpretation of the Large Slowdown of Water Dynamics at Stacked Phospholipid Membranes for Decreasing Hydration Level: All-Atom Molecular Dynamics

Hydration water determines the stability and function of phospholipid membranes as well as the interaction of membranes with other molecules. Experiments and simulations have shown that water dynamics slows down dramatically as the hydration decreases, suggesting that the interfacial water that dominates the average dynamics at low hydration is slower than water away from the membrane. Here, based on all-atom molecular dynamics simulations, we provide an interpretation of the slowdown of interfacial water in terms of the structure and dynamics of water–water and water–lipid hydrogen bonds (HBs). We calculate the rotational and translational slowdown of the dynamics of water confined in stacked phospholipid membranes at different levels of hydration, from completely hydrated to poorly hydrated membranes. For all hydrations, we analyze the distribution of HBs and find that water–lipids HBs last longer than water–water HBs and that at low hydration most of the water is in the interior of the membrane. We also show that water–water HBs become more persistent as the hydration is lowered. We attribute this effect (i) to HBs between water molecules that form, in turn, persistent HBs with lipids; (ii) to the hindering of the H-bonding switching between water molecules due to the lower water density at the interface; and (iii) to the higher probability of water–lipid HBs as the hydration decreases. Our interpretation of the large dynamic slowdown in water under dehydration is potentially relevant in understanding membrane biophysics at different hydration levels.

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

Neutron scattering techniques have been employed to investigate 1,2-dimyristoyl-sn -glycero-3-phosphocholine (DMPC) membranes in the form of multilamellar vesicles (MLVs) and deposited, stacked multilamellar-bilayers (MLBs), covering transitions from the gel to the liquid phase. Neutron diffraction was used to characterise the samples in terms of transition temperatures, whereas elastic incoherent neutron scattering (EINS) demonstrates that the dynamics on the sub-macromolecular length-scale and pico- to nano-second time-scale are correlated with the structural transitions through a discontinuity in the observed elastic intensities and the derived mean square displacements. Molecular dynamics simulations have been performed in parallel focussing on the length-, time- and temperature-scales of the neutron experiments. They correctly reproduce the structural features of the main gel-liquid phase transition. Particular emphasis is placed on the dynamical amplitudes derived from experiment and simulations. Two methods are used to analyse the experimental data and mean square displacements. They agree within a factor of 2 irrespective of the probed time-scale, i.e. the instrument utilized. Mean square displacements computed from simulations show a comparable level of agreement with the experimental values, albeit, the best match with the two methods varies for the two instruments. Consequently, experiments and simulations together give a consistent picture of the structural and dynamical aspects of the main lipid transition and provide a basis for future, theoretical modelling of dynamics and phase behaviour in membranes. The need for more detailed analytical models is pointed out by the remaining variation of the dynamical amplitudes derived in two different ways from experiments on the one hand and simulations on the other.

Strong Static Magnetic Fields Increase the Gel Signal in Partially Hydrated DPPC/DMPC Membranes

It was recently reported that static magnetic fields increase lipid order in the hydrophobic membrane core of dehydrated native plant plasma membranes [Poinapen, Soft Matter 9:6804-6813, 2013]. As plasma membranes are multicomponent, highly complex structures, in order to elucidate the origin of this effect, we prepared model membranes consisting of a lipid species with low and high melting temperature. By controlling the temperature, bilayers coexisting of small gel and fluid domains were prepared as a basic model for the plasma membrane core. We studied molecular order in mixed lipid membranes made of dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) using neutron diffraction in the presence of strong static magnetic fields up to 3.5 T. The contribution of the hydrophobic membrane core was highlighted through deuterium labeling the lipid acyl chains. There was no observable effect on lipid organization in fluid or gel domains at high hydration of the membranes. However, lipid order was found to be enhanced at a reduced relative humidity of 43%: a magnetic field of 3.5 T led to an increase of the gel signal in the diffraction patterns of 5%. While all biological materials have weak diamagnetic properties, the corresponding energy is too small to compete against thermal disorder or viscous effects in the case of lipid molecules. We tentatively propose that the interaction between the fatty acid chains’ electric moment and the external magnetic field is driving the lipid tails in the hydrophobic membrane core into a better ordered state.

Multi-component modeling of quasielastic neutron scattering from phospholipid membranes

We investigated molecular motions in the 0.3-350 ps time range of D2O-hydrated bilayers of 1-palmitoyl-oleoyl-sn-glycero-phosphocholine and 1,2-dimyristoyl-sn-glycero-phosphocholine in the liquid phase by quasielastic neutron scattering. Model analysis of sets of spectra covering scale lengths from 4.8 to 30 Å revealed the presence of three types of motion taking place on well-separated time scales: (i) slow diffusion of the whole phospholipid molecules in a confined cylindrical region; (ii) conformational motion of the phospholipid chains; and (iii) fast uniaxial rotation of the hydrogen atoms around their carbon atoms. Based on theoretical models for the hydrogen dynamics in phospholipids, the spatial extent of these motions was analysed in detail and the results were compared with existing literature data. The complex dynamics of protons was described in terms of elemental dynamical processes involving different parts of the phospholipid chain on whose motions the hydrogen atoms ride.

Multi-lamellar organization of fully deuterated lipid extracts of yeast membranes

Neutron scattering studies on mimetic biomembranes are currently limited by the low availability of deuterated unsaturated lipid species. In the present work, results from the first neutron diffraction experiments on fully deuterated lipid extracts from the yeast Pichia pastoris are presented. The structural features of these fully deuterated lipid stacks are compared with those of their hydrogenous analogues and with other similar synthetic systems. The influence of temperature and humidity on the samples has been investigated by means of small momentum-transfer neutron diffraction. All of the lipid extracts investigated self-assemble into multi-lamellar stacks having different structural periodicities; the stacking distances are affected by temperature and humidity without altering the basic underlying arrangement. At high relative humidity the deuterated and hydrogenous samples are similar in their multi-lamellar arrangement, being characterized by two main periodicities of ∼75 and ∼110 Å reflecting the presence of a large number of polar phospholipid molecules. Larger differences are found at lower relative humidity, where hydrogenous lipids are characterized by a larger single lamellar structure than that observed in the deuterated samples. In both cases the heterogeneity in composition is reflected in a wide structural complexity. The different behaviour upon dehydration can be related to compositional differences in the molecular composition of the two samples, which is attributed to metabolic effects related to the use of perdeuterated growth media.

Structure-based engineering of a minimal porin reveals loop-independent channel closure

Porins, like outer membrane protein G (OmpG) of Escherichia coli, are ideal templates among ion channels for protein and chemical engineering because of their robustness and simple architecture. OmpG shows fast transitions between open and closed states, which were attributed to loop 6 (L6). As flickering limits single-channel-based applications, we pruned L6 by either 8 or 12 amino acids. While the open probabilities of both L6 variants resemble that of native OmpG, their gating frequencies were reduced by 63 and 81%, respectively. Using the 3.2 Å structure of the shorter L6 variant in the open state, we engineered a minimal porin (220 amino acids), where all remaining extramembranous loops were truncated. Unexpectedly, this minimized porin still exhibited gating, but it was 5-fold less frequent than in OmpG. The residual gating of the minimal pore is hence independent of L6 rearrangements and involves narrowing of the ion conductance pathway most probably driven by global stretching-flexing deformations of the membrane-embedded β-barrel.

High hydrostatic pressure effects investigated by neutron scattering on lipid multilamellar vesicles

The effects of high hydrostatic pressure on the structure and dynamics of model membrane systems were investigated using neutron scattering. Diffraction experiments show shifts of the pre- and main-phase transitions of multilamellar vesicles of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) to higher temperatures with increased pressure which are close to results observed previously by other techniques, namely (10.4 ± 1.0) K kbar(-1) and (20.0 ± 0.5) K kbar(-1) for the two transitions. Backscattering spectroscopy reveals that the mean square displacements in the liquid phase are about 10% smaller at 300 bar and about 20% smaller at 600 bar compared to atmospheric pressure, whereas in the gel phase below the main phase transition the mean square displacements show a smaller difference in the dynamics of the three pressure values within the studied pressure range.

Correlation between supercoiling and conformational motions of the bacterial flagellar filament

The bacterial flagellar filament is a very large macromolecular assembly of a single protein, flagellin. Various supercoiled states of the filament exist, which are formed by two structurally different conformations of flagellin in different ratios. We investigated the correlation between supercoiling of the protofilaments and molecular dynamics in the flagellar filament using quasielastic and elastic incoherent neutron scattering on the picosecond and nanosecond timescales. Thermal fluctuations in the straight L- and R-type filaments were measured and compared to the resting state of the wild-type filament. Amplitudes of motion on the picosecond timescale were found to be similar in the different conformational states. Mean-square displacements and protein resilience on the 0.1 ns timescale demonstrate that the L-type state is more flexible and less resilient than the R-type, whereas the wild-type state lies in between. Our results provide strong support that supercoiling of the protofilaments in the flagellar filament is determined by the strength of molecular forces in and between the flagellin subunits.

Structural and dynamical properties of reconstituted myelin sheaths in the presence of myelin proteins MBP and P2 studied by neutron scattering

The myelin sheath is a tightly packed, multilayered membrane structure wrapped around selected nerve axons in the central and the peripheral nervous system. Because of its electrical insulation of the axons, which allows fast, saltatory nerve impulse conduction, myelin is crucial for the proper functioning of the vertebrate nervous system. A subset of myelin-specific proteins is well-defined, but their influence on membrane dynamics, i.e. myelin stability, has not yet been explored in detail. We investigated the structure and the dynamics of reconstituted myelin membranes on a pico- to nanosecond timescale, influenced by myelin basic protein (MBP) and myelin protein 2 (P2), using neutron diffraction and quasi-elastic neutron scattering. A model for the scattering function describing molecular lipid motions is suggested. Although dynamical properties are not affected significantly by MBP and P2 proteins, they act in a highly synergistic manner influencing the membrane structure.

Dynamics measured by neutron scattering correlates with the organization of bioenergetics complexes in natural membranes from hyperthermophile and mesophile bacteria

Various models on membrane structure and organization of proteins and complexes in natural membranes emerged during the last years. However, the lack of systematic dynamical studies to complement structural investigations hindered the establishment of a more complete picture of these systems. Elastic incoherent neutron scattering gives access to the dynamics on a molecular level and was applied to natural membranes extracted from the hyperthermophile Aquifex aeolicus and the mesophile Wolinella succinogenes bacteria. The results permitted to extract a hierarchy of dynamic flexibility and atomic resilience within the samples, which correlated with the organization of proteins in bioenergetics complexes and the functionality of the membranes.

Dielectric relaxation change of water upon phase transition of a lipid bilayer probed by terahertz time domain spectroscopy

We investigate the influence of the 1, 2-ditetradecanoyl-sn-glycero-3-phosphocholine lipid bilayer phases on the water reorientation dynamics with terahertz time domain spectroscopy. The phase of the lipids was controlled by the temperature in the range of 14-35 °C. During the gel-to-fluid phase transition, the hydration water ratio drastically changed from 0.3 to 0.6. The absorption coefficient of the hydration water increased with the temperature in the gel phase and then decreased in the fluid phase. The dielectric relaxation time of the lipid solution decreased initially but then increased after the phase transition. This indicates that the hydration water reorientation dynamics are restricted by lipids and that this phenomenon is pronounced in a biologically relevant fluid phase.

Dynamical features in cationic micelles of varied chain length

Chain length is one of the parameters controlling the structural arrangement of micelle monomers, such that one can tailor the monomers for different applications, but the effect of chain length on the dynamical behavior of micelles is unknown. In this article, we report a study on the effect of varying chain length on the dynamical behavior of alkyltrimethylammonium bromide (C(n)TAB) micelles (n = 10, 12, 14, and 16) using incoherent quasielastic neutron scattering (QENS). The data analysis clearly shows the presence of two distinct motions: global motion of whole micelles and faster internal motions of the C(n)TAB monomers. The global diffusion is Fickian in nature, whereas the internal motions can be described with a model that considers the motions of the headgroup and the hydrophobic alkyl chain separately. Methyl groups in the headgroup undergo 3-fold jump rotations, and the hydrogen atoms belonging to the alkyl chain undergo localized translational diffusion. The hydrogen atoms belonging to the alkyl chain are confined within spherical volumes that increase linearly along the C(n)TAB chain: the hydrogen atoms closer to the headgroup move within smaller spheres with lower diffusion coefficients than those farther from the headgroup. The main result is that, with increasing chain length, the dynamics of the C(n)TAB monomer is greatly affected: diffusion is reduced and occurs in smaller spheres, and residence times are increased. Global motion is also hindered with increased chain length.

Temperature-dependent vibrational and conformational dynamics of photosystem II membrane fragments from spinach investigated by elastic and inelastic neutron scattering

Vibrational and conformational protein dynamics of photosystem II (PS II) membrane fragments from spinach were investigated by elastic and inelastic incoherent neutron scattering (EINS and IINS). As to the EINS experiments, the average atomic mean square displacement values of PS II membrane fragments hydrated at a relative humidity of 57% exhibit a dynamical transition at ~230K. In contrast, the dynamical transition was absent at a relative humidity of 44%. These findings are in agreement with previous studies which reported a "freezing" of protein mobility due to dehydration (Pieper et al. (2008) Eur. Biophys. J. 37: 657-663) and its correlation with an inhibition of electron transfer from Q(A)(-) to Q(B) (Kaminskaya et al. (2003) Biochemistry 42, 8119-8132). IINS spectra of a sample hydrated at a relative humidity of 57% show a distinct Boson peak at ~7.5meV at 20K, which shifts towards lower energy values upon temperature increase to 250K. This unexpected effect is interpreted in terms of a "softening" of the protein matrix along with the onset of conformational protein dynamics as revealed by the EINS experiments. Information on the density of vibrational states of pigment-protein complexes is important for a realistic calculation of excitation energy transfer kinetics and spectral lineshapes and is often routinely obtained by optical line-narrowing spectroscopy at liquid helium temperature. The data presented here demonstrate that IINS is a valuable experimental tool in determining the density of vibrational states not only at cryogenic, but also at nearly physiological temperatures up to 250K. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Dynamic properties of photosystem II membranes at physiological temperatures characterized by elastic incoherent neutron scattering. Increased flexibility associated with the inactivation of the oxygen evolving complex

Elastic incoherent neutron scattering (EINS), a non-invasive technique which is capable of measuring the mean square displacement of atoms in the sample, has been widely used in biology for exploring the dynamics of proteins and lipid membranes but studies on photosynthetic systems are scarce. In this study we investigated the dynamic characteristics of Photosystem II (PSII) membrane fragments between 280 and 340 K, i.e., in the physiological temperature range and in the range of thermal denaturation of some of the protein complexes. The mean square displacement values revealed the presence of a hydration-sensitive transition in the sample between 310 and 320 K, suggesting that the oxygen evolving complex (OEC) plays an important role in the transition. Indeed, in samples in which the OEC had been removed by TRIS- or heat-treatments (323 and 333 K) no such transition was found. Further support on the main role of OEC in these reorganizations is provided by data obtained from differential scanning calorimetry experiments, showing marked differences between the untreated and TRIS-treated samples. In contrast, circular dichroism spectra exhibited only minor changes in the excitonic interactions below 323 K, showing that the molecular organization of the pigment-protein complexes remains essentially unaffected. Our data, along with earlier incoherent neutron scattering data on PSII membranes at cryogenic temperatures (Pieper et al., Biochemistry 46:11398-11409, 2007), demonstrate that this technique can be applied to characterize the dynamic features of PSII membranes, and can be used to investigate photosynthetic membranes under physiologically relevant experimental conditions.

High-throughput optical quality control of lipid multilayers fabricated by dip-pen nanolithography

Surface supported phospholipid multilayers are promising materials for nanotechnology because of their tendency to self-organize, their innate biocompatibility, the possibility to encapsulate other materials within the multilayers, and the ability to control the multilayer thickness between ∼ 2 and 100 nm during fabrication. Dip-pen nanolithography (DPN) is an atomic force microscopy (AFM) based fabrication method that allows high-throughput fabrication and integration of a variety of micro- and nanostructured materials including lipid multilayers, with areal throughputs on the scale of cm(2) min(-1). Although multilayer thickness is a critical feature that determines the functionality of the lipid multilayer structures (for instance as carriers for other materials as well as optical scattering properties), reliable height characterization by AFM is slow (on the order of µm(2) min(-1)) and a bottleneck in the lithographic process. Here we describe a novel optical method to reliably measure the height of fluorescent multilayers with thicknesses above 10 nm, and widths above the optical diffraction limit based on calibrating the fluorescence intensity using one-time AFM height measurements. This allows large surface areas to be rapidly and quantitatively characterized using a standard fluorescence microscope. Importantly, different pattern dimensions (0D dots, 1D lines or 2D squares) require different calibration parameters, indicating that shape influences the optical properties of the structured lipid multilayers. This method has general implications in the systematic and high-throughput optical characterization of nanostructure-function relationships.

Interplay between hydration water and headgroup dynamics in lipid bilayers

In this study, the interplay between water and lipid dynamics has been investigated by broadband dielectric spectroscopy and modulated differential scanning calorimetry (MDSC). The multilamellar lipid bilayer system 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) has been studied over a broad temperature range at three different water contents: about 3, 6, and 9 water molecules per lipid molecule. The results from the dielectric relaxation measurements show that at temperatures <250 K the lipid headgroup rotation is described by a super-Arrhenius temperature dependence at the lowest hydration level and by the Arrhenius law at the highest hydration level. This difference in the temperature dependence of the lipid headgroup rotation can be explained by the increasing interaction between the headgroups with decreasing water content, which causes their rotational motion to be more cooperative in character. The main water relaxation shows an anomalous dependence on the water content in the supercooled and glassy regime. In contrast to the general behavior of interfacial water, the water dynamics is fastest in the driest sample and its temperature dependence is best described by a super-Arrhenius temperature dependence. The best explanation for this anomalous behavior is that the water relaxation becomes more determined by fast local lipid motions than by the intrinsic water dynamics at low water contents. In support for this interpretation is the finding that the relaxation time of the main water process is faster than that in most other host systems at temperatures below 180 K. Thus, the dielectric relaxation data show clearly the strong interplay between water and lipid dynamics; the water influences the lipid dynamics and vice versa. In the MDSC data, we observe a weak enthalpy relaxation at 203 K for the driest sample and at 179 K for the most hydrated sample, attributed to the freezing-in of the lipid headgroup rotation observed in the dielectric data, since this motion reaches a time scale of about 100 s at about the same temperatures.