Dynamics and state of lipid bilayer-internal water unraveled with solution state 1H dynamic nuclear polarization

The dynamics and state of lipid bilayer-internal hydration water of unilamellar lipid vesicles dispersed in solutions is characterized. This study was enabled by a recently developed technique based on Overhauser dynamic nuclear polarization (DNP)-driven amplification of (1)H nuclear magnetic resonance (NMR) signal of hydration water. This technique can, in the full presence of bulk water, selectively quantify the translational dynamics of hydration water within ～10 ? around spin labels that are specifically introduced to the local volume of interest within the lipid bilayer. With this approach, the local apparent diffusion coefficients of internal water at different depths of the lipid bilayer were determined. The modulation of these values as a response to external stimuli, such as the addition of sodium chloride or ethanol and the lipid phase transitions, that alter the fluctuations of bilayer interfaces together with the activation energy values of water diffusivity shows that water is not individually and homogeneously solvating lipid's hydrocarbon tails in the lipid bilayer. We provide experimental evidence that instead, water and the lipid membrane comprise a heterogeneous system whose constituents include transient hydrophobic water pores or water structures traversing the lipid bilayer. We show how these transient pore structures, as key vehicles for passive water transport can better reconcile our experimental data with existing literature data on lipid bilayer hydration and dynamics.     

Novel two-band ratiometric fluorescence probes with different location and orientation in phospholipid membranes

3-hydroxyflavone (3-HF) derivatives are very attractive fluorescence sensors due to their ability to respond to small changes in their microenvironment via a dramatic alteration of the relative intensities of their two well-separated emission bands. We developed fluorescence probes with locations at different depths and orientations of 3-HF moiety in the phospholipid bilayer, which determine their fluorescence behavior. While the spectral shifts of the probes correlate with their binding site polarity, the intensity ratio is a complex parameter that is also sensitive to the local hydration. We demonstrate that even the deeply located probes sense this hydration effect, which can be modulated by the charge of the lipid heads and is anisotropic with respect to the bilayer plane. Thus the two-band ratiometric fluorescence probes can provide multiparametric information on the properties of lipid membranes at different depths.     

Heating and temperature gradients of lipid bilayer samples induced by RF irradiation in MAS solid‐state NMR experiments

The MAS solid‐state NMR has been a powerful technique for studying membrane proteins within the native‐like lipid bilayer environment. In general, RF irradiation in MAS NMR experiments can heat and potentially destroy expensive membrane protein samples. However, under practical MAS NMR experimental conditions, detailed characterization of RF heating effect of lipid bilayer samples is still lacking. Herein, using ¹H chemical shift of water for temperature calibration, we systematically study the dependence of RF heating on hydration levels and salt concentrations of three lipids in MAS NMR experiments. Under practical ¹H decoupling conditions used in biological MAS NMR experiments, three lipids show different dependence of RF heating on hydration levels as well as salt concentrations, which are closely associated with the properties of lipids. The maximum temperature elevation of about 10 °C is similar for the three lipids containing 200% hydration, which is much lower than that in static solid‐state NMR experiments. The RF heating due to salt is observed to be less than that due to hydration, with a maximum temperature elevation of less than 4 °C in the hydrated samples containing 120 mmol l^?1 of salt. Upon RF irradiation, the temperature gradient across the sample is observed to be greatly increased up to 20 °C, as demonstrated by the remarkable broadening of ¹H signal of water. Based on detailed characterization of RF heating effect, we demonstrate that RF heating and temperature gradient can be significantly reduced by decreasing the hydration levels of lipid bilayer samples from 200% to 30%. Copyright © 2016 John Wiley & Sons, Ltd.     

Ab initio electron correlated studies on the intracluster reaction of NO+ (H2O)(n) → H3O+ (H2O)(n-2) (HONO) (n = 4 and 5)

Hydrated nitrosonium ion clusters NO(+)(H(2)O)(n) (n = 4 and 5) were investigated by using MP2/aug-cc-pVTZ level of theory to clarify isomeric reaction pathways for formation of HONO and fully hydrated hydride ions. We found some new isomers and transition state structures in each hydration number, whose lowest activation energies of the intracluster reactions were found to be 4.1 and 3.4 kcal mol(-1) for n = 4 and n = 5, respectively. These thermodynamic properties and full quantum mechanical molecular dynamics simulation suggest that product isomers with HONO and fully hydrated hydride ions can be obtained at n = 4 and n = 5 in terms of excess hydration binding energies which can overcome these activation barriers.     

Raman Spectroscopic Signature of Ectoine Conformations in Bulk Solution and Crystalline State

Recent crystallographic results revealed conformational changes of zwitterionic ectoine upon hydration. By means of confocal Raman spectroscopy and density functional theory calculations, we present a detailed study of this transformation process as part of a Fermi resonance analysis. The corresponding findings highlight that all resonant couplings are lifted upon exposure to water vapor as a consequence of molecular binding processes. The importance of the involved molecular groups for water binding and conformational changes upon hydration is discussed. Our approach further shows that the underlying rapid process can be reversed by carbon dioxide saturated atmospheres. For the first time, we also confirm that the conformational state of ectoine in aqueous bulk solution coincides with crystalline ectoine in its dihydrate state, thereby highlighting the important role of a few bound water molecules.     

Efficient formation of giant liposomes through the gentle hydration of phosphatidylcholine films doped with sugar

Giant liposomes, or giant vesicles, are cell-size (approximately 5-100 microm) compartments enclosed with phospholipid bilayers, and have often been used in biological research. They are usually generated using hydration methods, "electroformation" and "gentle hydration (or natural swelling)", in which dry lamellar films of phospholipids are hydrated with aqueous solutions. In gentle hydration, however, giant liposomes are difficult to prepare from an electrostatically neutral phospholipid because lipid lamellae cannot repel each other. In this study, we demonstrate the efficient formation of giant liposomes using the gentle hydration of neutral phospholipid (dioleoyl phosphatidylcholine, DOPC) dry films doped with nonelectrolytic monosaccharides (glucose, mannose, and fructose). A mixture of DOPC and such a sugar in an organic solvent (chloroform/methanol) was evaporated to form the films, which were then hydrated with distilled water or Tris buffers containing sodium chloride. Under these conditions, giant liposomes spontaneously formed rapidly and assumed a swollen cell-sized spherical shape with low lamellarity, whereas giant liposomes from pure DOPC films had multilamellar lipid layers, miscellaneous shapes and smaller sizes. This observation indicates that giant unilamellar vesicles (GUVs) of DOPC can be obtained efficiently through the gentle hydration of sugar-containing lipid dry films because repulsion between lipid lamellae is enhanced by the osmosis induced by dissolved sugar.     

Accurate predictions of nonpolar solvation free energies require explicit consideration of binding-site hydration

Continuum solvation methods are frequently used to increase the efficiency of computational methods to estimate free energies. In this paper, we have evaluated how well such methods estimate the nonpolar solvation free-energy change when a ligand binds to a protein. Three different continuum methods at various levels of approximation were considered, viz., the polarized continuum model (PCM), a method based on cavity and dispersion terms (CD), and a method based on a linear relation to the solvent-accessible surface area (SASA). Formally rigorous double-decoupling thermodynamic integration was used as a benchmark for the continuum methods. We have studied four protein-ligand complexes with binding sites of varying solvent exposure, namely the binding of phenol to ferritin, a biotin analogue to avidin, 2-aminobenzimidazole to trypsin, and a substituted galactoside to galectin-3. For ferritin and avidin, which have relatively hidden binding sites, rather accurate nonpolar solvation free energies could be obtained with the continuum methods if the binding site is prohibited to be filled by continuum water in the unbound state, even though the simulations and experiments show that the ligand replaces several water molecules upon binding. For the more solvent exposed binding sites of trypsin and galectin-3, no accurate continuum estimates could be obtained, even if the binding site was allowed or prohibited to be filled by continuum water. This shows that continuum methods fail to give accurate free energies on a wide range of systems with varying solvent exposure because they lack a microscopic picture of binding-site hydration as well as information about the entropy of water molecules that are in the binding site before the ligand binds. Consequently, binding affinity estimates based upon continuum solvation methods will give absolute binding energies that may differ by up to 200 kJ/mol depending on the method used. Moreover, even relative energies between ligands with the same scaffold may differ by up to 75 kJ/mol. We have tried to improve the continuum solvation methods by adding information about the solvent exposure of the binding site or the hydration of the binding site, and the results are promising at least for this small set of complexes.     

Hydration changes accompanying nucleic acid intercalation reactions:volumetric characterizations

We use high precision ultrasonic velocimetric and densimetric techniques to determine at 25 degrees C the changes in volume, deltaV, and adiabatic compressibility, deltaK(S), that accompany the binding of ethidium to the poly(rA)poly(rU), poly(dAdT)poly(dAdT), poly(dGdC)poly(dGdC), and poly(dIdC)poly(dIdC) duplexes, as well as to the poly(rU)poly(rA)poly(rU) triplex. The binding of ethidium to each of the duplexes and the triplex is accompanied by negative changes in volume, deltaV, and adiabatic compressibility, deltaK(S). We discuss the basis for relating macroscopic and microscopic properties, particularly, emphasizing how measured changes in volume and compressibility can be quantitatively interpreted in terms of the differential hydration properties of DNA and RNA structures in their ligand-free and ligand-bound states. We also estimate the entropic cost of intercalation-induced changes in hydration of each of the nucleic acid structures and the drug. In general, our results emphasize the vital role of hydration in modulating the energetics of drug-DNA binding, while also underscoring the fact that hydration must be carefully taken into account in analysis and prediction of the energetics of nucleic acid recognition.     

Dehydration of solute-lipid systems: hydration forces analysis

Sorption isotherms were obtained for a range of lipid/sugar/water mixtures. These were analysed using a simple hydration forces formalism. The results demonstrate that this simple analysis can be used to estimate dehydration parameters for these relatively complex systems. This in turn provides some insight into the location and role of sugars in the hydration behaviour of lipid systems. The relevance of these results to the phase behaviour of lipid/sugar mixtures during dehydration are discussed.     

Dipolar Relaxation Dynamics at the Active Site of an ATPase Regulated by Membrane Lateral Pressure

The active transport of ions across biological membranes requires their hydration shell to interact with the interior of membrane proteins. However, the influence of the external lipid phase on internal dielectric dynamics is hard to access by experiment. Using the octahelical transmembrane architecture of the copper‐transporting P_1B‐type ATPase from Legionella pneumophila as a model structure, we have established the site‐specific labeling of internal cysteines with a polarity‐sensitive fluorophore. This enabled dipolar relaxation studies in a solubilized form of the protein and in its lipid‐embedded state in nanodiscs. Time‐dependent fluorescence shifts revealed the site‐specific hydration and dipole mobility around the conserved ion‐binding motif. The spatial distribution of both features is shaped significantly and independently of each other by membrane lateral pressure.     

A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer

A series of molecular dynamics simulations of dimyristoylphosphatidylcholine bilayers, with different levels of hydration and temperature, were performed to examine the influence of hydration on properties of lipid membranes. Structural and dynamical properties such as area per lipid, electron densities, order parameters for all CH bonds and water, diffusion, and reorientation autocorrelation functions were determined and were all found to be affected by changes in the hydration level. The simulations give an overall picture of the bilayer going to a more ordered state when the hydration level is reduced. Lipid headgroups were found to adopt an orientation more parallel to the membrane plane when the water content was decreased. Dynamical properties such as lipid diffusion and relaxation of reorientation time correlation functions were found to become slower with the removal of water. Our simulation results generally agree with experimental data in cases where such data are available. One important conclusion drawn is that while structural properties are affected only moderately dynamical properties are affected very strongly by a decrease of water content.     

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. II. Two-dimensional infrared and peak shift observables within different theoretical approximations

In a previous report, we calculated the infrared absorption spectrum and both the isotropic and anisotropic pump-probe signals for the OD stretch of isotopically dilute water in dilauroylphosphatidylcholine (DLPC) multi-bilayers as a function of the lipid hydration level. These results were then compared to recent experimental measurements and are in generally good agreement. In this paper, we will further investigate the structure and dynamics of hydration water using molecular dynamics simulations and calculations of the two-dimensional infrared and vibrational echo peak shift observables for hydration water in DLPC membranes. These observables have not yet been measured experimentally, but future comparisons may provide insight into spectral diffusion processes and hydration water heterogeneity. We find that at low hydration levels the motion of water molecules inside the lipid membrane is significantly arrested, resulting in very slow spectral diffusion. At higher hydration levels, spectral diffusion is more rapid, but still slower than in bulk water. We also investigate the effects of several common approximations on the calculation of spectroscopic observables by computing these observables within multiple levels of theory. The impact of these approximations on the resulting spectra affects our interpretation of these measurements and reveals that, for example, the cumulant approximation, which may be valid for certain systems, is not a good approximation for a highly heterogeneous environment such as hydration water in lipid multi-bilayers.     

Interactions of alkali metal chlorides with phosphatidylcholine vesicles

We study the interaction of alkali metal chlorides with lipid vesicles made of palmitoyloleoylphosphatidylcholine (POPC). An elaborate set of techniques is used to investigate the binding process at physiological conditions. The alkali cation binding to POPC is characterized thermodynamically using isothermal titration calorimetry. The isotherms show that for all ions in the alkali group the binding process is endothermic, counterintuitively to what is expected for Coulomb interactions between the slightly negatively charged POPC liposomes and the cations. The process is entropy driven and presumably related to the liberation of water molecules from the hydration shells of the ions and the lipid headgroups. The measured molar enthalpies of the binding of the ions follows the Hofmeister series. The binding constants were also estimated, whereby lithium shows the strongest affinity to POPC membranes, followed by the rest of the ions according to the Hofmeister series. Cation adsorption increases the net surface potential of the vesicles as observed from electrophoretic mobility and zeta potential measurements. While lithium adsorption leads to slightly positive zeta potentials above a concentration of 100 mM, the adsorption of the rest of the ions mainly causes neutralization of the membrane. This is the first study characterizing the binding equilibrium of alkali metal chlorides to phosphatidylcholine membranes at physiological salt concentrations.     

INFLUENCE OF HYDRATION ON THE INTERNAL DYNAMICS OF HEN EGG WHITE LYSOZYME IN THE DRY STATE

Abstract— Proteins exist in a predominately aqueous solvent environment. Hydration of the protein surface significantly affects many aspects of the protein's structure and function; these effects may be related to the molecular dynamics of the protein. We have examined the influence of hydration on the internal dynamics of hen egg white lysozyme using room-temperature phosphorescence from the intrinsic tryptophan residues. Powders of lyophilized lysozyme were hydrated in a phosphorimeter using a flow system that allowed for continuous manipulation of relative humidity over the range 0–92%; this system allowed us to directly compare intensity differences that result from changes in hydration. Lysozyme phosphorescence intensity decreased as a function of hydration over the entire relative humidity range; the decrease was not linear but appeared to occur in distinct phases. The phosphorescence intensity decays were multiexponential over the hydration range studied, and hydration had the largest influence on the long lifetime component. These data suggest that the protein exists in multiple, static conformations in the dry state and that water binding to polar (as opposed to charged) sites on the protein surface induces local and/or global softening of the protein structure.     

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

We investigated molecular motions on a picosecond timescale of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) model membranes as a function of hydration by using elastic and quasielastic neutron scattering. Two different hydrations corresponding to approximately nine and twelve water molecules per lipid were studied, the latter being the fully hydrated state. In our study, we focused on head group motions by using chain deuterated lipids. Information on in-plane and out-of-plane motions could be extracted by using solid supported DMPC multilayers. Our studies confirm and complete former investigations by Ko?nig et al. [J. Phys. II (France) 2, 1589 (1992)] and Rheinsta?dter et al. [Phys. Rev. Lett. 101, 248106 (2008)] who described the dynamics of lipid membranes, but did not explore the influence of hydration on the head group dynamics as presented here. From the elastic data, a clear shift of the main phase transition from the P(β) ripple phase to the L(α) liquid phase was observed. Decreasing water content moves the transition temperature to higher temperatures. The quasielastic data permit a closer investigation of the different types of head group motion of the two samples. Two different models are needed to fit the elastic incoherent structure factor and corresponding radii were calculated. The presented data show the strong influence hydration has on the head group mobility of DMPC.     

Energy relaxation dynamics of the hydration complex of hydroxide

We use polarization-resolved mid-infrared pump-probe spectroscopy to study the dynamics of the hydration shells of hydroxide ions (OH(-)). We excite the OH stretch vibrations of H(2)O molecules solvating the OH(-) ion and observe that this excitation decays with a relaxation time constant T(1) of 200 fs. This relaxation is followed by a thermalization process that becomes slower with increasing concentration of OH(-). The prethermalized state is observed to be anisotropic, showing that the energy of the excited OH stretch vibrations is dissipated within the hydration complex. The anisotropy of the prethermalized state decays both as a result of the reorientation of the OH(-) hydration complex and heat diffusion from the excited complexes to unexcited complexes. Modeling the anisotropy data at different concentrations allows for an accurate estimate of the number of water molecules in the hydration shell of OH(-), the reorientation dynamics of the OH(-) hydration complex, and the molecular-scale heat diffusivity.     

Controlled Release of a Sparingly Water‐Soluble Anticancer Drug through pH‐Responsive Functionalized Gold‐Nanoparticle‐Decorated Liposomes

The binding and detachment of carboxyl‐modified gold nanoparticles from liposomes is used for controlled drug delivery. This study reveals that the binding and detachment of nanoparticles from liposomes depends on the degree of hydration of the liposomes. Liposomes with a lower hydration level undergo stronger electrostatic interactions with negatively charged gold nanoparticles, thus leading to a slower detachment of the carboxyl‐modified gold nanoparticles under gastric conditions. Therefore, under gastric conditions, gold‐nanoparticle‐decorated dipalmitoylphosphatidylcholine (DPPC) liposomes exhibit an at least ten‐times‐slower drug release compared to gold‐nanoparticle‐decorated 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC) liposomes, although both liposomes in the bare state fail to pursue controlled release. Our study also reveals that one can modulate the drug‐release rate by simply varying the concentration of nanoparticles. This study highlights a novel strategy for the controlled release of drug molecules from liposomes.     

Systematic implicit solvent coarse graining of dimyristoylphosphatidylcholine lipids

We have used systematic structure‐based coarse graining to derive effective site–site potentials for a 10‐site coarse‐grained dimyristoylphosphatidylcholine (DMPC) lipid model and investigated their state point dependence. The potentials provide for the coarse‐grained model the same site–site radial distribution functions, bond and angle distributions as those computed in atomistic simulations carried out at four different lipid–water molar ratios. It was shown that there is a non‐negligible dependence of the effective potentials on the concentration at which they were generated, which is also manifested in the properties of the lipid bilayers simulated using these potentials. Thus, effective potentials computed at low lipid concentration favor to more condensed and ordered structure of the bilayer with lower average area per lipid, while potentials obtained at higher lipid concentrations provide more fluid‐like structure. The best agreement with the reference data and experiment was achieved using the set of potentials derived from atomistic simulations at 1:30 lipid:water molar ratio providing fully saturated hydration of DMPC lipids. Despite theoretical limitations of pairwise coarse‐grained potentials expressed in their state point dependence, all the resulting potentials provide a stable bilayer structure with correct partitioning of different lipid groups across the bilayer as well as acceptable values of the average lipid area, compressibility and orientational ordering. In addition to bilayer simulations, the model has proven its robustness in modeling of self‐aggregation of lipids from randomly dispersed solution to ordered bilayer structures, bicelles, and vesicles. © 2014 Wiley Periodicals, Inc.     

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. I. Infrared spectra and ultrafast pump-probe observables

The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.