Molecular dynamics simulation study of the interaction of trehalose with lipid membranes

The interactions of the cryoprotective agent trehalose with a lipid membrane made of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine at 323 K were studied by means of molecular dynamics simulations. It was observed that trehalose binds to the phospholipid headgroups with its main axis parallel to the membrane normal. Trehalose establishes hydrogen bonds with the carbonyl and phosphate groups and replaces water molecules from the lipid headgroup. Notably, the number of hydrogen bonds (HBs) that the membrane made with its environment was conserved after trehalose binding. The HBs between lipid and trehalose have a longer lifetime than those established between lipid and water. The binding of the sugar does not produce changes either in the lipid area or in the lipid order parameter. The effect of trehalose on the dipole potential is in agreement with experimental results. The contribution of the different components to the membrane dipole potential was analyzed. It was observed that the binding of trehalose produces changes in the different components and the sugar itself contributes to the surface potential due to the polarization of its hydroxyl in the interface.     

A molecular dynamics study of ethanol-water hydrogen bonding in binary structure I clathrate hydrate with CO2

Guest-host hydrogen bonding in clathrate hydrates occurs when in addition to the hydrophilic moiety which causes the molecule to form hydrates under high pressure-low temperature conditions, the guests contain a hydrophilic, hydrogen bonding functional group. In the presence of carbon dioxide, ethanol clathrate hydrate has been synthesized with 10% of large structure I (sI) cages occupied by ethanol. In this work, we use molecular dynamics simulations to study hydrogen bonding structure and dynamics in this binary sI clathrate hydrate in the temperature range of 100-250 K. We observe that ethanol forms long-lived (>500 ps) proton-donating and accepting hydrogen bonds with cage water molecules from both hexagonal and pentagonal faces of the large cages while maintaining the general cage integrity of the sI clathrate hydrate. The presence of the nondipolar CO(2) molecules stabilizes the hydrate phase, despite the strong and prevalent alcohol-water hydrogen bonding. The distortions of the large cages from the ideal form, the radial distribution functions of the guest-host interactions, and the ethanol guest dynamics are characterized in this study. In previous work through dielectric and NMR relaxation time studies, single crystal x-ray diffraction, and molecular dynamics simulations we have observed guest-water hydrogen bonding in structure II and structure H clathrate hydrates. The present work extends the observation of hydrogen bonding to structure I hydrates.     

Stabilizing role of lithium in structures of complex oxide compounds as an instrument for crystal chemical design

Various hydrogen bond lifetime distribution functions, used to describe the breaking and formation dynamics of these bonds in a computer experiment, are examined and relationships between them are found. The procedures for calculating these functions by the molecular dynamics method are described and the results for water models of 3456 molecules at 310 K are reported. The peak of short-lived spurious H-bonds, which results from short-time violations of hydrogen bonding criteria induced by dynamic intermolecular vibrations of molecules, prevails in the types of distributions most often referred to in the literature. A special distribution that appears to have not been used before is proposed. Along with short-lived bonds, it manifests long-lived ones whose lifetime is determined by the genuine, or random, hydrogen bond breaking rather than by dynamic. A technique to exclude dynamic effects and reveal the genuine H-bond breaking is proposed. This allows the evaluation of the average lifetime of “true” H-bonds that turns out to exceed 3 ps.     

Anomalous dielectric relaxation of water molecules at the surface of an aqueous micelle

Dielectric relaxation of aqueous solutions of micelles, proteins, and many complex systems shows an anomalous dispersion at frequencies intermediate between those corresponding to the rotational motion of bulk water and that of the organized assembly or macromolecule. The precise origin of this anomalous dispersion is not well-understood. In this work we employ large scale atomistic molecular dynamics simulations to investigate the dielectric relaxation (DR) of water molecules in an aqueous micellar solution of cesium pentadecafluorooctanoate. The simulations clearly show the presence of a slow component in the moment-moment time correlation function [PhiMW(t)] of water molecules, with a time constant of about 40 ps, in contrast to only 9 ps for bulk water. Interestingly, the orientational time correlation function [Cmu(t)] of individual water molecules at the surface exhibits a component with a time constant of about 19 ps. We show that these two time constants can be related by the well-known micro-macrorelations of statistical mechanics. In addition, the reorientation of surface water molecules exhibits a very slow component that decays with a time constant of about 500 ps. An analysis of hydrogen bond lifetime and of the rotational relaxation in the coordinate frame fixed on the micellar body seems to suggest that the 500 ps component owes its origin to the existence of an extended hydrogen bond network of water molecules at the surface. However, this ultraslow component is not found in the total moment-moment time correlation function of water molecules in the solution. The slow DR of hydration water is found to be well correlated with the slow solvation dynamics of cesium ions at the water-micelle interface.     

OPTICAL SPECTROSCOPY OF BILAYER MEMBRANES

Abstract— The absorption and fluorescence spectroscopy of natural and model bilayer lipid membranes is reviewed. Basic structural features of biological membranes and the relative advantages of black lipid membranes (BLM) and of liposomes are discussed. Theoretical considerations show that the wavelengths of absorption maxima are affected by the refractive index and dielectric constant of the medium surrounding the chromophore. Techniques of obtaining photoelectric action spectra, direct absorption spectra, and reflection spectra of BLM are described. Polarized spectra can give information about the orientation of membrane constituents and show, for example, that the porphyrin ring of chlorophyll in BLM is tilted at 45 ± 5° to the membrane surface. Absorption maxima of chlorophyll in BLM are compared with solution spectra of various chlorophyll adducts and aggregates. It is concluded that chlorophyll in BLM exists largely as solvated monomer and dimer (or oligomer), depending on concentration, and is not coordinated with water. From the theory of fluorescence spectroscopy it follows that aggregation and the polarity of the environment affect the fluorescence yield and lifetime of a membrane component, and also the wavelength of its emission maximum. The microviscosity of the membrane matrix affects the anisotropy of fluorescence. Techniques of steady-state fluorescence spectroscopy and of fluorescence lifetime measurements are reviewed. Examples of the use of fluorescent probes in membrane studies are given. Certain probes such as anilinonaphthalene sulfonate (ANS) preferentially bind to membrane proteins. The location of a probe in a particular membrane region can be pinpointed from its fluorescence yield and emission maximum. The orientation of the hydrocarbon chains of membrane lipids has been found, from fluorescence polarization of certain probes, to be normal to the membrane surface as postulated a priori on the basis of the lipid bilayer model. Anisotropy of fluorescence shows that elongated probe molecules rotate rapidly about their long axes when surrounded by phospholipids but become immobilized when bound to proteins. Changes in intensity and anisotropy of fluorescence as function of temperature have demonstrated the existence of phase transitions and phase equilibria of membrane lipids. Excimer fluorescence has been used as a measure of the available lipid core volume of membranes. Mechanisms of energy transfer between membrane components are reviewed. The theoretical dependence of energy transfer on distance and orientation for several rigid and fluid membrane models is discussed in terms of the structural information it can provide. Fluorescence sensitization resulting from energy transfer within and across bilayer membranes has been demonstrated in various systems. Quantitative measurement of energy transfer efficiency in BLM has shown that such transfer is about five times more efficient than in solutions at comparable donor-acceptor distances. Lipid membranes can be viewed as structures which maintain their components at high concentrations, in a reactive state, and at favourable orientations.     

Interaction of zervamicin IIB with lipid bilayers. Molecular dynamics study

In this work we have studied the interaction of zervamicin IIB (ZrvIIB) with the model membranes of eukaryotes and prokaryotes using all-atom molecular dynamics. In all our simulations zervamicin molecule interacted only with lipid headgroups but did not penetrate the hydrophobic core of the bilayers. During the interaction with the prokaryotic membrane zervamicin placed by its N-termini towards the lipids and rotated at an angle of 40° relatively to the bilayer surface. In the case of eukaryotic membrane zervamicin stayed in the water and located parallel to the membrane surface. We compared hydrogen bonds between peptide and lipids and concluded that interactions of ZrvIIB with prokaryotic membrane are stronger than those with eukaryotic one. Also it was shown that two zervamicin molecules formed dimer and penetrated deeper in the area of lipid headgroups.     

The Effects of Various Alcohols on the Stability of Melittin: A Molecular Dynamics Study

In this study, 200 ps molecular dynamics simulations were conducted to investigate the effects of various alcohols on the structural stability of melittin. The averaged helicity of melittin remained 80% in pure butanol, whereas it was below 60% both in pure water and in pure methanol. The α‐helix propensity of melittin increased with the aliphatic chain length of the alcohol. Charge‐charge interaction between Lys21 and Arg24 and polar‐nonpolar interaction between Trp19 and Arg22 are probably responsible for the higher structural integrity of the C‐terminal α‐helix over the N‐terminal one. The weaker dielectric constant of longer aliphatic chain length of alcohol possibly reduces the hydrogen bonding between amide protons and surrounding solvent molecules and simultaneously promotes the intramolecular hydrogen bonding in melittin and therefore stabilizes the secondary structure of melittin. The effect of various alcohols on stabilizing melittin is most likely due to their ability to form clusters on the surface of melittin effectively, favoring the formation of intramolecular hydrogen bonds instead of intermolecular hydrogen bonds and promoting the formation of stable α‐helices.     

Hydrophilic carbon nanotube membrane enhanced interfacial evaporation for desalination

Carbon nanotube-based(CNT-based) interfacial evaporation material is one of the most potential materials for solar desalination. Here, we studied the evaporation rate of the CNT-based membranes with different hydrophilic and hydrophobic chemical modified surfaces using molecular dynamic simulations.We found that the hydrogen bonding density among water molecules at the interface is a key factor in enhancing the evaporation rate. For a hydrophilic CNT-based membrane, the strong interactions betwe...     

Description of collective effects in computer models of water

The dynamics of a system containing 3456 water molecules in a cubic cell with periodic boundary conditions at 297 K was simulated. The time dependence of distances between oxygen atoms was examined for many pairs of molecules. These distances often oscillate around a certain average value over long periods of time (10 ps and longer). These average values can be about 2.8 Å (hydrogen bond) or much larger, up to 12–13 Å and more. This suggests that big groups of molecules are involved in a concerted motion. Lists of hydrogen bonds in 50 configurations divided by an interval of about 1 ps are compared. The average lifetime of a hydrogen bond is about 7 ps. The network of hydrogen bonds is colored according to their lifetimes for one of the configurations. The bonds that live longer than 7 ps form an infinite cluster. The bonds that live longer than 8 ps join to form a great number of finite clusters including several hundreds of nodes (molecules). These clusters contain few closed cycles. Even the bonds that live longer than 20 ps are united into clusters each containing two or three nodes (molecules). The self-diffusion coefficient for molecules involved in long-lived bonds is likely to be slightly smaller than that for molecules which do not participate in these bonds.     

Long-lived interfacial vibrations of water

We have observed long-lived OH-stretch (nu(OH)) excitations (v = 1) in water during ultrafast laser ablation by a mid-infrared pulse tuned to the nu(OH) absorption maximum. The spectrum of excitations is measured using incoherent anti-Stokes Raman spectroscopy. Relative to the equilibrium water spectrum, these excitations evidence a narrowed (100 cm(-1) fwhm) and blue-shifted (3600 cm(-1) peak) transition. The excited-state lifetime is T1 > 200 ps, compared to 0.2 ps in bulk water. In the early stages of the ablation process, the water mean density decreases rapidly, which breaks up the hydrogen bonding. The long-lived species is attributed to nu(OH) excitations on water molecules associated with interfaces, having broken hydrogen bonds which cannot be rapidly reformed as in the liquid state.     

Hydration of DNA-binding biological active compounds: EHF dielectrometry and molecular modeling results

We report hydration properties of several DNA-binding ligands (pharmaceutical drug caffeine; mutagens proflavine, ethidium bromide, propidium iodide; polyamines putrescine, spermine), which are able to interact with DNA in different modes via external binding, intercalation or minor groove binding. We show that the detection of the bound water molecules and the estimation of their amount in aqueous solutions of ligands can be efficiently carried out using the measurements of complex dielectric permittivity of the solutions in the millimeter range of radio waves. Our dielectrometric data are combined with the results of the molecular modeling including quantum chemical calculations and Monte Carlo simulations. We show that number of water molecules able to form hydrogen bonds with donor–acceptor groups of ligands correlates with hydration numbers taken from the literature or obtained in EHF dielectrometry experiment. The latter indicates that the EHF dielectrometry method is sensitive for the tightly bound water molecules.     

Design and synthesis of sphingomyelin-cholesterol conjugates and their formation of ordered membranes

A lipid raft is a cholesterol (Chol)-rich microdomain floating in a sea of lipid bilayers. Although Chol is thought to interact preferentially with sphingolipids such as sphingomyelin (SM), rather than with glycerophospholipids, the origin of the specific interaction has remained unresolved, primarily because of the high mobility of lipid molecules and weak intermolecular interactions. In this study, we synthesized SM-Chol conjugates with functionally designed linker portions to restrain Chol mobility and examined their formation of ordered membranes by a detergent insolubility assay, fluorescence anisotropy experiments, and fluorescence-quenching assay. In all of the tests, membranes prepared from the conjugates showed properties of ordered domains comparable to a SM-Chol (1:1) membrane. To gain insight into the structure of bilayers composed from the conjugates, we performed molecular dynamics simulations with 64 molecules of the conjugates, which suggested that the conjugates form a stable bilayer structure by bending at the linker portion and, mostly, reproduce the hydrogen bonds between the SM and Chol portions. These results imply that the molecular recognition between SM and Chol in an ordered domain is essentially reproduced by the conjugated molecules and, thus, demonstrates that these conjugate molecules could potentially serve as molecular probes for understanding molecular recognition in lipid rafts.     

Tautomerism in 7-hydroxyquinoline: a combined experimental and theoretical study in water

Tautomerism in the ground and excited states of 7-hydroxyquinoline (7HQ) was studied in different solvents using steady-state and lifetime spectroscopic measurements, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. Equilibrium between the enol and the keto/zwitterion tautomers exists in 7HQ, which is solvent-dependent. Of the solvents used in this study, only in water does the absorbance spectrum of 7HQ show absorption from both the enol and zwitterion tautomers. In addition, in aqueous media, fluorescence is observed from the zwitterion tautomer only, which is attributed to self-quenching of the enol fluorescence by energy transfer to the ground-state zwitterion tautomer and energetically favorable excited-state proton transfer. Solvation of the hydrogen bonding sites of 7HQ was studied in binary mixtures of 1,4-dioxane and water, and three water molecules were estimated to connect the polar sites and induce intermolecular proton transfer. The results are confirmed by DFT calculations showing that three water molecules are the minimum number required to form a stable solvent wire. Mapping the water density around the polar sites using MD simulations shows well-defined hydrogen bonds around the amino and hydroxyl groups of the enol tautomer and slightly less well-defined hydrogen bonds for the zwitterion tautomer. The presence of three-member water wires connecting the polar centers in 7HQ is evident in the MD simulations. The results point to the unique spectral signatures of 7HQ in water, which make this molecule a potential probe to detect the presence of water in nanocavities of macromolecules.     

Properties of a water layer on hydrophilic and hydrophobic self-assembled monolayer surfaces: A molecular dynamics study

The microscopic behaviors of a water layer on different hydrophilic and hydrophobic surfaces of well ordered self-assembled monolayers (SAMs) are studied by molecular dynamics simulations. The SAMs consist of 18-carbon alkyl chains bound to a silicon(111) substrate, and the characteristic of its surface is tuned from hydrophobic to hydrophilic by using different terminal functional groups ( CH 3 , COOH). In the simulation, the properties of water membranes adjacent to the surfaces of SAMs were reported by comparing pure water in mobility, structure, and orientational ordering of water molecules. The results suggest that the mobility of water molecules adjacent to hydrophilic surface becomes weaker and the molecules have a better ordering. The distribution of hydrogen bonds indicates that the number of water-water hydrogen bonds per water molecule tends to be lower. However, the mobility of water molecules and distribution of hydrogen bonds of a water membrane in hydropho- bic system are nearly the same as those in pure water system. In addition, hydrogen bonds are mainly formed between the hydroxyl of the COOH group and water molecules in a hydrophilic system, which is helpful in understanding the structure of interfacial water.     

Strikingly Different Effects of Hydrogen Bonding on the Photodynamics of Individual Nucleobases in DNA: Comparison of Guanine and Cytosine

Ab initio surface hopping dynamics calculations were performed to study the photophysical behavior of cytosine and guanine embedded in DNA using a quantum mechanical/molecular mechanics (QM/MM) approach. It was found that the decay rates of photo excited cytosine and guanine were affected in a completely different way by the hydrogen bonding to the DNA environment. In case of cytosine, the geometrical restrictions exerted by the hydrogen bonds did not influence the relaxation time of cytosine significantly due to the generally small cytosine ring puckering required to access the crossing region between excited and ground state. On the contrary, the presence of hydrogen bonds significantly altered the photodynamics of guanine. The analysis of the dynamics indicates that the major contribution to the lifetime changes comes from the interstrand hydrogen bonds. These bonds considerably restricted the out-of-plane motions of the NH(2) group of guanine which are necessary for the ultrafast decay to the ground state. As a result, only a negligible amount of trajectories decayed into the ground state for guanine embedded in DNA within the simulation time of 0.5 ps, while for comparison, the isolated guanine relaxed to the ground state with a lifetime of about 0.22 ps. These examples show that, in addition to phenomena related to electronic interactions between nucleobases, there also exist relatively simple mechanisms in DNA by which the lifetime of a nucleobase is significantly enhanced as compared to the gas phase.     

Role of bound water in biological membrane structure: fluorescence and infrared studies.

Bound water is a major component of biological membranes and is required for the structural stability of the lipid bilayer. It has also been postulated that it is involved in water transport, membrane fusion, and mobility of membrane proteins and lipids. We have measured the fluorescence emission of membrane-bound 1-anilino-8-naphthalenesulfonate (ANS) and the infrared spectra of membranes, both as a function of hydration. ANS fluorescence is sensitive to polarity and fluidity of the membrane-aqueous interface, while infrared absorption is sensitive to the hydrogen bonding and vibrational motion of water and membrane proteins and lipids. The fluorescence results provide evidence of increasing rigidity and/or decreasing polarity of the membrane-aqueous interface with removal of water. The membrane infrared spectra show prominent hydration-dependent changes in a number of bands with possible assignments to cholesterol (vinyl CH bend, OH stretch), protein (amide A, II, V), and bound water (OH stretch). Further characterization of the bound water should allow its incorporation into current models of membrane structure and give insight into the role of membrane hydration in cell surface function.     

Sensitivity of hydrogen bond lifetime dynamics to the presence of ethanol at the interface of a phospholipid bilayer

Atomistic molecular dynamics simulations of a fully hydrated liquid crystalline lamellar phase of a dimyrystoylphosphatidylcholine lipid bilayer containing ethanol at 1:1 composition as well as of the pure lamellar phase of the bilayer have been performed. Detailed analyses have been carried out to investigate the effects of ethanol, if any, on the lifetime dynamics of lipid-water and water-water hydrogen bonds in the hydration layer of the lipid headgroups. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed the identification of the bound states at the bilayer interface and the quantification of the dynamic equilibrium between the bound and the free water molecules, in terms of time-dependent relaxation rates. The calculations show that the overall relaxation of phosphate-water hydrogen bonds is faster in the presence of ethanol. Studies of the residence time and the number fluctuation of the hydration layer water molecules reveal that the presence of ethanol molecules decreases the rigidity of the lipid hydration layer.     

Solvent‐Modulated Influence of Intramolecular Hydrogen‐Bonding on the Conformational Properties of the Hydroxymethyl Group in Glucose and Galactose: A Molecular Dynamics Simulation Study

Intramolecular hydrogen‐bonding (H‐bonding) is commonly regarded as a major determinant of the conformation of (bio)molecules. However, in an aqueous environment, solvent‐exposed H‐bonds are likely to represent only a marginal (possibly adverse) conformational driving as well as steering force. For example, the hydroxymethyl rotamers of glucose and galactose permitting the formation of an intramolecular H‐bond with the adjacent hydroxyl group are not favored in water but, in the opposite, least populated. This is because the solvent‐exposed H‐bond is dielectrically screened as well as subject to intense H‐bonding competition by the water molecules. In the present study, the effect of a decrease in the solvent polarity on this rotameric equilibrium is probed using molecular dynamics simulation. This is done by considering six physical solvents (H₂O, DMSO , MeOH , CHC l₃, CC l₄, and vacuum), along with 19 artificial water‐like solvent models for which the dielectric permittivity and H‐bonding capacity can be modulated independently via a scaling of the O–H distance and of the atomic partial charges. In the high polarity solvents, the intramolecular H‐bond is observed, but arises as an opportunistic consequence of the proximity of the H‐bonding partners in a given rotameric state. Only when the polarity of the solvent is decreased does the intramolecular H‐bond start to induce a conformational pressure on the rotameric equilibrium. The artificial solvent series also reveals that the effects of the solvent permittivity and of its H‐bonding capacity mutually enhance each other, with a slightly larger influence of the permittivity. The hydroxymethyl conformation in hexopyranoses appears to be particularly sensitive to solvent‐polarity effects because the H‐bond involving the hydroxymethyl group is only one out of up to five H‐bonds capable of forming a network around the ring.     

Structure and dynamics of the hydration shells of the Al3+ ion

First principles simulations of the hydration shells surrounding Al3+ ions are reported for temperatures near 300 degrees C. The predicted six water molecules in the octahedral first hydration shell were found to be trigonally coordinated via hydrogen bonds to 12 s shell water molecules in agreement with the putative structure used to analyze the x-ray data, but in disagreement with the results reported from conventional molecular dynamics using two-and three-body potentials. Bond lengths and angles of the water molecules in the first and second hydration shells and the average radii of these shells also agreed very well with the results of the x-ray analysis. Water transfers into and out of the second solvation shell were observed to occur on a picosecond time scale via a dissociative mechanism. Beyond the second shell the bonding pattern substantially returned to the tetrahedral structure of bulk water. Most of the simulations were done with 64 solvating water molecules (20 ps). Limited simulations with 128 water molecules (7 ps) were also carried out. Results agreed as to the general structure of the solvation region and were essentially the same for the first and second shell. However, there were differences in hydrogen bonding and Al-O radial distribution function in the region just beyond the second shell. At the end of the second shell a nearly zero minimum in the Al-O radial distribution was found for the 128 water system. This minimum is less pronounced minimum found for the 64 water system, which may indicate that sizes larger than 64 may be required to reliably predict behavior in this region.     

Dielectric Properties of Solutions Isomorphous with Water

The dielectric properties of aqueous solutions of ammonia, hydrogen fluoride, hydrogen peroxide, and ammonium halides are determined by the number of hydrogen bonds formed by the solute species with the molecules forming the water framework. For such solutions, the product of the concentration coefficient of the dielectric permittivity by the ratio of the number of hydrogen bonds per solute molecule and per water molecule is constant. In aqueous solutions with a water-like structure, the mechanism of the dielectric polarization is the same as in pure water (proton-activation mechanism) and is determined by proton fluctuations on the line of hydrogen bonds.__________Translated from Zhurnal Obshchei Khimii, Vol. 75, No. 1, 2005, pp. 39–45.Original Russian Text Copyright © 2005 by Potapov, Parkhomenko.