Anomalous and anisotropic nanoscale diffusion of hydration water molecules in fluid lipid membranes

Solid-State Molecular Protonics Devices of Solid-Supported Biological Membranes Reveal the Mechanism of Long-Range Lateral Proton Transport

Lateral proton transport (PT) on the surface of biological membranes is a fundamental biochemical process in the bioenergetics of living cells, but a lack of available experimental techniques has resulted in a limited understanding of its mechanism. Here, we present a molecular protonics experimental approach to investigate lateral PT across membranes by measuring long-range (70 μm) lateral proton conduction via a few layers of lipid bilayers in a solid-state-like environment, i.e., without having bulk water surrounding the membrane. This configuration enables us to focus on lateral proton conduction across the surface of the membrane while decoupling it from bulk water. Hence, by controlling the relative humidity of the environment, we can directly explore the role of water in the lateral PT process. We show that proton conduction is dependent on the number of water molecules and their structure and on membrane composition, where we explore the role of the headgroup, the tail saturation, the membrane phase, and membrane fluidity. The measured PT as a function of temperature shows an inverse temperature dependency, which we explain by the desorption and adsorption of water molecules into the solid membrane platform. We explain our findings by discussing the role of percolating hydrogen bonding within the membrane structure in a Grotthuss-like mechanism.

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.

Neutron scattering studies on dynamics of lipid membranes

Neutron scattering methods are powerful tools for the study of the structure and dynamics of lipid bilayers in length scales from sub Å to tens to hundreds nm and the time scales from sub ps to μs. These techniques also are nondestructive and, perhaps most importantly, require no additives to label samples. Because the neutron scattering intensities are very different for hydrogen- and deuterium-containing molecules, one can replace the hydrogen atoms in a molecule with deuterium to prepare on demand neutron scattering contrast without significantly altering the physical properties of the samples. Moreover, recent advances in neutron scattering techniques, membrane dynamics theories, analysis tools, and sample preparation technologies allow researchers to study various aspects of lipid bilayer dynamics. In this review, we focus on the dynamics of individual lipids and collective membrane dynamics as well as the dynamics of hydration water.

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.

Water dynamics in human cancer and non-cancer tissues

Normal-to-malignant transformation is a poorly understood process associated with cellular biomechanical properties. These are strongly dependent on the dynamical behaviour of water, known to play a fundamental role in normal cellular activity and in the maintenance of the three-dimensional architecture of the tissue and the functional state of biopolymers. In this study, quasi-elastic neutron scattering was used to probe the dynamical behaviour of water in human cancer specimens and their respective surrounding normal tissue from breast and tongue, as an innovative approach for identifying particular features of malignancy. This methodology has been successfully used by the authors in human cells and was the first study of human tissues by neutron scattering techniques. A larger flexibility was observed for breast versus tongue tissues. Additionally, different dynamics were found for malignant and non-malignant specimens, depending on the tissue: higher plasticity for breast invasive cancer versus the normal, and an opposite effect for tongue. The data were interpreted in the light of two different water populations within the samples: one displaying bulk-like dynamics (extracellular and intracellular/cytoplasmic) and another with constrained flexibility (extracellular/interstitial and intracellular/hydration layers).

Influence of water models on water movement through AQP1

Water diffusion through membrane proteins is a key aspect of cellular function. Essential processes of cellular metabolism are driven by osmotic pressure, which depends on water channels. Membrane proteins such as aquaporins (AQPs) are responsible for enabling water permeation through the cell membrane. AQPs are highly selective, allowing only water and relatively small polar molecules to cross the membrane. Experimentally, estimation of water flux through membrane proteins is still a challenge, and hence, accurate simulations of water permeation are of particular importance. We present a numerical study of water diffusion through AQP1 comparing three water models: TIP3P, OPC, and TIP4P/2005. Bulk diffusion, diffusion permeability, and osmotic permeability are computed and compared among all models. The results show that there are significant differences between TIP3P (a particularly widespread model for simulations of biological systems) and the more recently developed TIP4P/2005 and OPC models. We demonstrate that OPC and TIP4P/2005 reproduce protein-water interactions and dynamics in very good agreement with experimental data. From this study, we find that the choice of the water model has a significant effect on the computed water dynamics as well as its molecular behavior within a biological nanopore.

Disentangling Self-Atomic Motions in Polyisobutylene by Molecular Dynamics Simulations

We present fully atomistic molecular dynamics simulations on polyisobutylene (PIB) in a wide temperature range above the glass transition. The cell is validated by direct comparison of magnitudes computed from the simulation and measured by neutron scattering on protonated samples reported in previous works. Once the reliability of the simulation is assured, we exploit the information in the atomic trajectories to characterize the dynamics of the different kinds of atoms in PIB. All of them, including main-chain carbons, show a crossover from Gaussian to non-Gaussian behavior in the intermediate scattering function that can be described in terms of the anomalous jump diffusion model. The full characterization of the methyl-group hydrogen motions requires accounting for rotational motions. We show that the usually assumed statistically independence of rotational and segmental motions fails in this case. We apply the rotational rate distribution model to correlation functions calculated for the relative positions of methyl-group hydrogens with respect to the carbon atom at which they are linked. The contributions to the vibrational density of states are also discussed. We conclude that methyl-group rotations are coupled with the main-chain dynamics. Finally, we revise in the light of the simulations the hypothesis and conclusions made in previously reported neutron scattering investigations on protonated samples trying to address the origin of the dielectric β-process.

Anesthetics significantly increase the amount of intramembrane water in lipid membranes

The potency of anesthesia was directly linked to the partitioning of the drug molecules in cell membranes by Meyer and Overton. Many molecules interact with lipid bilayers and lead to structural and functional changes. It remains an open question which change in membrane properties is responsible for a potential anesthetic effect or if anesthetics act by binding to direct targets. We studied the effect of ethanol, diethyl ether and isoflurane on the water distribution in lipid bilayers by combining all-atom molecular dynamics simulations and neutron diffraction experiments. The simulations show strong membrane-drug interactions with partitioning coefficients of 38%, 92% and 100% for ethanol, diethyl ether and isoflurane, respectively, and provide evidence for an increased water partitioning in the membrane core. The amount of intramembrane water molecules was experimentally determined by selectively deuterium labeling lipids, anesthetic drug and water molecules in neutron diffraction experiments. Four additional water molecules per lipid were observed in the presence of ethanol. Diethyl ether and isoflurane were found to significantly increase the amount of intramembrane water by 25% (8 water molecules). This increase in intramembrane water may contribute to the non-specific interactions between anesthetics and lipid membranes.

Quasi-Elastic Neutron Scattering Studies on Hydration Water in Phospholipid Membranes

The dynamic behavior of hydration water in phospholipid membranes has been investigated to understand the relationship between water and biological molecules using various experimental techniques. Quasi-elastic neutron scattering (QENS) is an effective method for this purpose because the dynamic behaviors of both water and lipid molecules could be identified by using selective deuteration. In addition, the measurable ranges from the 10⁻¹² to 10⁻⁹ s time scale and the 10⁻¹¹ to 10⁻⁸ m length scale are suitable to investigate the slowing down of water molecules due to their interaction with lipid membranes. In this mini-review, QENS experiments on the dynamic behavior of hydration water molecules in neighboring phospholipid membranes are summarized.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Exploring fast proton transfer events associated with lateral proton diffusion on the surface of membranes

Proton diffusion (PD) across biological membranes is a fundamental process in many biological systems, and much experimental and theoretical effort has been employed for deciphering it. Here, we report on a spectroscopic probe, which can be tightly tethered to the membrane, for following fast (nanosecond) proton transfer events on the surface of membranes. Our probe is composed of a photoacid that serves as our light-induced proton source for the initiation of the PD process. We use our probe to follow PD, and its pH dependence, on the surface of lipid vesicles composed of a zwitterionic headgroup, a negative headgroup, a headgroup that is composed only from the negative phosphate group, or a positive headgroup without the phosphate group. We reveal that the PD kinetic parameters are highly sensitive to the nature of the lipid headgroup, ranging from a fast lateral diffusion at some membranes to the escape of protons from surface to bulk (and vice versa) at others. By referring to existing theoretical models for membrane PD, we found that while some of our results confirm the quasi-equilibrium model, other results are in line with the nonequilibrium model.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Impact of Electric Fields on the Nanoscale Behavior of Lipid Monolayers at the Surface of Graphite in Solution

The nanoscale organization and dynamics of lipid molecules in self-assembled membranes is central to the biological function of cells and in the technological development of synthetic lipid structures as well as in devices such as biosensors. Here, we explore the nanoscale molecular arrangement and dynamics of lipids assembled in monolayers at the surface of highly ordered pyrolytic graphite (HOPG), in different ionic solutions, and under electrical potentials. Using a combination of atomic force microscopy and fluorescence recovery after photobleaching, we show that HOPG is able to support fully formed and fluid lipid membranes, but mesoscale order and corrugations can be observed depending on the type of the lipid considered (1,2-dioleoyl- sn-glycero-3-phosphocholine, 1,2-dioleoyl- sn-glycero-3-phospho-l-serine (DOPS), and 1,2-dioleoyl-3-trimethylammoniumpropane) and the ion present (Na⁺, Ca²⁺, Cl^-). Interfacial solvation forces and ion-specific effects dominate over the electrostatic changes induced by moderate electric fields (±1.0 V vs Ag/AgCl reference electrode) with particularly marked effects in the presence of calcium, and for DOPS. Our results provide insights into the interplay between the molecular, ionic, and electrostatic interactions and the formation of dynamical ordered structures in fluid lipid membranes.

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.

Influence of Cholesterol on the Dynamics of Hydration in Phospholipid Bilayers

We investigate the dynamics of interfacial waters in dipalmitoylphosphatidylcholine (DPPC) bilayers upon the addition of cholesterol, by molecular dynamics simulations. Our data reveal that the inclusion of cholesterol modifies the membrane aqueous interfacial dynamics: waters diffuse faster, their rotational decay time is shorter, and the DPPC/water hydrogen bond dynamics relaxes faster than in the pure DPPC membrane. The observed acceleration of the translational water dynamics agrees with recent experimental results, in which, by means of NMR techniques, an increment of the surface water diffusivity is measured upon the addition of cholesterol. A microscopic analysis of the lipid/water hydrogen bond network at the interfacial region suggests that the mechanism underlying the observed water mobility enhancement is given by the rupture of a fraction of interlipid water bridge hydrogen bonds connecting two different DPPC molecules, concomitant to the formation of new lipid/solvent bonds, whose dynamics is faster than that of the former. The consideration of a simple two-state model for the decay of the hydrogen bond correlation function yielded excellent results, obtaining two well-separated characteristic time scales: a slow one (∼250 ps) associated with bonds linking two DPPC molecules, and a fast one (∼15 ps), related to DPPC/solvent bonds.

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.

Description of Hydration Water in Protein (Green Fluorescent Protein) Solution

The structurally and dynamically perturbed hydration shells that surround proteins and biomolecules have a substantial influence upon their function and stability. This makes the extent and degree of water perturbation of practical interest for general biological study and industrial formulation. We present an experimental description of the dynamical perturbation of hydration water around green fluorescent protein in solution. Less than two shells (∼5.5 Å) were perturbed, with dynamics a factor of 2-10 times slower than bulk water, depending on their distance from the protein surface and the probe length of the measurement. This dependence on probe length demonstrates that hydration water undergoes subdiffusive motions (τ ∝ q^-2.5 for the first hydration shell, τ ∝ q^-2.3 for perturbed water in the second shell), an important difference with neat water, which demonstrates diffusive behavior (τ ∝ q^-2). These results help clarify the seemingly conflicting range of values reported for hydration water retardation as a logical consequence of the different length scales probed by the analytical techniques used.

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.

Amyloid-β(25-35) peptides aggregate into cross-β sheets in unsaturated anionic lipid membranes at high peptide concentrations

One of the hallmarks of Alzheimer's disease is the formation of protein plaques in the brain, which mainly consist of amyloid-β peptides of different lengths. While the role of these plaques in the pathology of the disease is not clear, the mechanism behind peptide aggregation is a topic of intense research and discussion. Because of their simplicity, synthetic membranes are promising model systems to identify the elementary processes involved. We prepared unsaturated zwitterionic/anionic lipid membranes made of 1-palmitoyl-2-oleoyl-sn-glycero-phosphocholine (POPC) and 1,2-dimyristoyl-sn-glycero-3-phospho-l-serine (DMPS) at concentrations of POPC/3 mol% DMPS containing 0 mol%, 3 mol%, 10 mol%, and 20 mol% amyloid-β25-35 peptides. Membrane-embedded peptide clusters were observed at peptide concentrations of 10 and 20 mol% with a typical cluster size of ∼11 μm. Cluster density increased with peptide concentration from 59 (±3) clusters per mm(2) to 920 (±64) clusters per mm(2), respectively. While monomeric peptides take an α-helical state when embedded in lipid bilayers at low peptide concentrations, the peptides in peptide clusters were found to form cross-β sheets and showed the characteristic pattern in X-ray experiments. The presence of the peptides was accompanied by an elastic distortion of the bilayers, which can induce a long range interaction between the peptides. The experimentally observed cluster patterns agree well with Monte Carlo simulations of long-range interacting peptides. This interaction may be the fundamental process behind cross-β sheet formation in membranes and these sheets may serve as seeds for further growth into amyloid fibrils.

Non-Brownian diffusion in lipid membranes: Experiments and simulations

The dynamics of constituents and the surface response of cellular membranes-also in connection to the binding of various particles and macromolecules to the membrane-are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder-through the addition of cholesterol or proteins-and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments-the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.