Effect of trehalose on a phospholipid membrane under mechanical stress

Explicit solvent molecular dynamics simulations were used to investigate at atomic resolution the effect of trehalose on a hydrated phospholipid bilayer under mechanical stress (stretching forces imposed in the form of negative lateral pressure). Simulations were performed in the absence or presence of trehalose at 325 K, and with different values for negative lateral pressure. In the concentration regime (2 molal) and range of lateral pressures (1 to −250 bar) investigated, trehalose was found to interact directly with the membrane, partially replacing water molecules in the formation of hydrogen bonds with the lipid headgroups. Similar to previous findings in the context of thermal stress, the number, degree of bridging, and reaching depth of these hydrogen bonds increased with the magnitude of perturbation. However, at the concentration considered, trehalose was not sufficient to preserve the integrity of the membrane structure and to prevent its extreme elongation (and possible disruption) under the effect of stretching forces.     

Molecular simulation study of phospholipid bilayers and insights of the interactions with disaccharides

Molecular simulations of hydrated dipalmitoylphosphatidylcholine lipid bilayers have been performed for temperatures in the range of 250-450 K. The area per headgroup increases with temperature from 58 to 77 A(2). Other properties such as hydration number, alkyl tail order parameter, diffusion coefficients, and radial distribution functions exhibit a clear dependence on temperature. Simulations of bilayers have also been performed in the presence of two disaccharides, namely trehalose and sucrose, at concentrations of up to 18 wt % (lipid-free basis). The simulated area per headgroup of the bilayer is not affected by the presence of the disaccharides, suggesting that the overall structure of the bilayer remains undisturbed. The results of simulations reveal that the interaction of disaccharide molecules with the bilayer occurs at the surface of the bilayer, and it is governed by the formation of multiple hydrogen bonds to specific groups of the lipid. Disaccharide molecules are observed to adopt specific conformations to fit onto the surface topology of the bilayer, often interacting with up to three different lipids simultaneously. At high disaccharide concentrations, the results of simulations indicate that disaccharides can serve as an effective replacement for water under anhydrous conditions, which helps explain their effectiveness as lyophilization agents for liposomes and cells.     

Under the influence of alcohol: the effect of ethanol and methanol on lipid bilayers

Extensive microscopic molecular dynamics simulations have been performed to study the effects of short-chain alcohols, methanol and ethanol, on two different fully hydrated lipid bilayer systems (POPC and DPPC) in the fluid phase at 323 K. It is found that ethanol has a stronger effect on the structural properties of the membranes. In particular, the bilayers become more fluid and permeable: ethanol molecules are able to penetrate through the membrane in typical timescales of approximately 200 ns, whereas for methanol that timescale is considerably longer, at least of the order of microseconds. A closer examination exposes a number of effects due to ethanol. Hydrogen-bonding analysis reveals that a large fraction of ethanols is involved in hydrogen bonds with lipids. This in turn is intimately coupled to the ordering of hydrocarbon chains: we find that binding to an ethanol decreases the order of the chains. We have also determined the dependence of lipid-chain ordering on ethanol concentration and found that to be nonmonotonous. Overall, we find good agreement with NMR and micropipette studies.     

Protective mechanisms of α,α-trehalose revealed by molecular dynamics simulations

Molecular dynamics simulations were carried out to model aqueous solution with different concentration of α,α-trehalose, one kind of non-reducing sugars possessing outstanding freeze-drying protective effect on biological system. The dihedral angles of the intraglycosidic linkage in trehalose were measured to estimate its structure rigidity. The dynamics and hydrogen bonding properties were studied by calculating the self-diffusion coefficient of trehalose and the distributions and lifetimes of various types of H-bonds in the solution. Through analysing the results as well as comparison with another common sugar sucrose, the freeze-drying protective mechanism of trehalose was explained at molecule level. First, trehalose is able to maintain the local structure around it as a frame due to its relatively rigid conformation. Second, the addition of trehalose restrains the water molecules from rearrangement as a result of low mobility, thus reduces the probability of freezing; trehalose has lower diffusion coefficient than water and bigger thermal diffusivity, which are favourable for vitrification. Third, the formation of H-bonds between trehalose and water and between trehalose molecules is the essence of the protective effect. Trehalose does not work via strengthening the H-bonds formed between water molecules (W–W H-bonds), instead of which it breaks the potential tetrahedral pattern of W–W H-bonds, thus suppresses the tendency of ice formation. It was also found that trehalose realises protective action better at higher concentration as far as this study is concerned.     

Bilayer undulation dynamics in unilamellar phospholipid vesicles: Effect of temperature,cholesterol and trehalose

We report a combined dynamic light scattering (DLS) and neutron spin-echo (NSE) study on the local bilayer undulation dynamics of phospholipid vesicles composed of 1,2-dimyristoyl-glycero-3-phosphatidylcholine (DMPC) under the influence of temperature and the additives cholesterol and trehalose. The additives affect vesicle size and self-diffusion. Mechanical properties of the membrane and corresponding bilayer undulations are tuned by changing lipid headgroup or acyl chain properties through temperature or composition. On the local length scale, changes at the lipid headgroup influence the bilayer bending rigidity κ less than changes at the lipid acyl chain: We observe a bilayer softening around the main phase transition temperature T_m of the single lipid system, and stiffening when more cholesterol is added, in concordance with literature. Surprisingly, no effect on the mechanical properties of the vesicles is observed upon the addition of trehalose.     

Dependence of trehalose protective action on the initial phase state of dipalmitoylphosphatidylcholine bilayers

Dipalmitoylphosphatidylcholine (DPPC) bilayers hydrated in the presence of trehalose were equilibrated at various temperatures (4, 20, and 60 degrees C) corresponding to the crystalline Lc, gel L beta', and liquid-crystalline L alpha phases, respectively, and then desiccated at these temperatures or freeze-dried at -80 degrees C to ca. DPPC dihydrate. The thermotropic behavior of the resulting DPPC/trehalose mixtures was investigated by differential scanning calorimetry and found to be dependent not only on the trehalose concentration but also on the phase state of the hydrated bilayers prior to their drying. Trehalose was most effective when the desiccation was carried out from the L alpha phase at 60 degrees C. In this case, one trehalose molecule per two DPPC molecules was sufficient to depress the melting temperature from values typical of DPPC dihydrate to 45 degrees C. Trehalose's influence decreased when dried from the L beta' phase and was significantly less pronounced when dried from the Lc phase. These data show that trehalose's protective influence depends on the initial phase state of the lipid bilayer and reaches its maximum in the liquid-crystalline state. The possible role of this effect in anhydrobiosis is pointed out.     

Solution effects on the thermotropic phase transition of unilamellar liposomes

We have investigated the effect of two monosaccharides, glucose and fructose, and two disaccharides, sucrose and trehalose, on the thermotropic phase transition of unilamellar extruded vesicles of DPPC. All the sugars investigated raise the main transition temperature (Tm) of some fraction of the lipid, but there are differences between the effect of glucose and the other three sugars. At low concentrations of glucose, Tm is lowered. At high concentrations of glucose there are two transitions, one with a low Tm and one with a high Tm. The data suggest that at low concentrations, all of the glucose present may bind to the bilayer and increase headgroup spacing by physical intercalation or increased hydration. The appearance of a Tm above that of pure hydrated DPPC suggests the possibility of the dehydration of some other population of phospholipid molecules. The other three sugars increase Tm, but at high concentrations of trehalose, sucrose, and fructose a second peak occurs at a low Tm. The other sugars appear to dehydrate the bilayer at low concentrations, but may show some binding or increased hydration of some portion of the lipid at very high concentrations. The sugar effects on unilamellar vesicles are strikingly different from the effects of these sugars on multilamellar vesicles.     

Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study

Sphingomyelin is enriched within lipid microdomains of the cell membrane termed lipid rafts. These microdomains play a part in regulating a variety of cellular events. Computer simulations of the hydrogen-bonding properties of sphingolipids, believed to be central to the organization of these domains, can delineate the possible molecular interactions that underlie this lipid structure. We have therefore used molecular dynamics simulations to unravel the hydrogen-bonding behavior of palmitoylsphingomyelin (PSM). A series of eight simulations of 3 ns each of a single PSM molecule in water showed that the sphingosine OH and NH groups can form hydrogen bonds with the phosphate oxygens of their own polar head, in agreement with NMR data. Simulations of PSM in a bilayer assembly were carried out for 8 ns with three different force field parameterizations. The major physico-chemical parameters of the simulated bilayer agree with those established experimentally. The sphingosine OH group was mainly involved in intramolecular hydrogen bonds, in contrast to the almost exclusive intermolecular hydrogen bonds formed by the amide NH moiety. During the bilayer simulations the intermolecular hydrogen bonds among lipids formed a dynamic network characterized by the presence of hydrogen-bonded lipid clusters of up to nine PSM molecules.     

Molecular dynamics simulation of dipalmitoylphosphatidylcholine membrane with cholesterol sulfate

Using the molecular dynamics simulation technique, we studied the changes occurring in a dipalmitoylphosphatidylcholine (DPPC):cholesterol (CH) membrane at 50 mol% sterol when cholesterol is replaced with cholesterol sulfate (CS). Our simulations were performed at constant pressure and temperature on a nanosecond time scale. We found that 1) the area per DPPC:CS heterodimer is greater than the area of the DPPC:CH heterodimer; 2) CS increases ordering of DPPC acyl chains, but to a lesser extent than CH; 3) the number of hydrogen bonds between DPPC and water is decreased in a CS-containing membrane, but CS forms more water hydrogen bonds than CH; and 4) the membrane dipole potential reverses its sign for a DPPC-CS membrane compared to a DPPC-CH bilayer. We also studied the changes occurring in lipid headgroup conformations and determined the location of CS molecules in the membrane. Our results are in good agreement with the data available from experiments.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

The cryoprotectant trehalose destabilises the bilayer organisation of Escherichia coli-derived membrane systems at elevated temperatures as determined by H and P-NMR

In this study, ²H and ³¹P-NMR techniques were used to study the effects of trehalose and glycerol on phase transitions and lipid acyl chain order of membrane systems derived from cells of E. coli unsaturated fatty acid auxotroph strain K1059, which was grown in the presence of [11,11-²H₂]-oleic acid or [11,11-²H₂]-elaidic acid. From an analysis of the temperature dependence of the quadrupolar splitting it could be concluded that neither 1 M trehalose or glycerol generally had any significant effect on the temperature of the lamellar gel to liquid-crystalline phase transition. In the case of the oleate-containing hydrated total lipid extract, glycerol but not trehalose caused a 5°C increase of this transition temperature. In general, both cryoprotectants induced an ordering of the acyl chains in the liquid-crystalline state. Trehalose and glycerol both decrease the bilayer to non-bilayer transition temperature of the hydrated lipid extract of oleate-grown cells by about 5°C, but only trehalose in addition induces an isotropic to hexagonal (H_II) phase transition. In the biological membranes, trehalose and not glycerol destabilised the lipid bilayer, and in the case of the E. coli spheroplasts, part of the induced non-bilayer structures is ascribed to a hexagonal (H_II) phase in analogy with the total lipids. Interestingly, 1 mM Mg²⁺ was a prerequisite for the destabilisation of the lipid bilayer. In the hydrated total lipid extract of E. coli grown on the more ordered elaidic acid, both transition temperatures were shifted about 20°C upwards compared with the oleate-containing lipid, but the effect of trehalose on the lipid phase behaviour was similar. The bilayer destabilising ability of trehalose might have implications for the possible protection of biological systems by (cryo-)protectants during dehydration, in that protection is unlikely to be caused by preventing the occurrence of polymorphic phase transitions.     

How alcohol chain-length and concentration modulate hydrogen bond formation in a lipid bilayer

Molecular dynamics simulations are used to measure the change in properties of a hydrated dipalmitoylphosphatidylcholine bilayer when solvated with ethanol, propanol, and butanol solutions. There are eight oxygen atoms in dipalmitoylphosphatidylcholine that serve as hydrogen bond acceptors, and two of the oxygen atoms participate in hydrogen bonds that exist for significantly longer time spans than the hydrogen bonds at the other six oxygen atoms for the ethanol and propanol simulations. We conclude that this is caused by the lipid head group conformation, where the two favored hydrogen-bonding sites are partially protected between the head group choline and the sn-2 carbonyl oxygen. We find that the concentration of the alcohol in the ethanol and propanol simulations does not have a significant influence on the locations of the alcohol/lipid hydrogen bonds, whereas the concentration does impact the locations of the butanol/lipid hydrogen bonds. The concentration is important for all three alcohol types when the lipid chain order is examined, where, with the exception of the high-concentration butanol simulation, the alcohol molecules having the longest hydrogen-bonding relaxation times at the favored carbonyl oxygen acceptor sites also have the largest order in the upper chain region. The lipid behavior in the high-concentration butanol simulation differs significantly from that of the other alcohol concentrations in the order parameter, head group rotational relaxation time, and alcohol/lipid hydrogen-bonding location and relaxation time. This appears to be the result of the system being very near to a phase transition, and one occurrence of lipid flip-flop is seen at this concentration.     

The modulation of membrane fluidity by hydrogenation processes. II. Homogeneous catalysis and model biomembranes

A homogeneous catalyst, chlorotris (triphenylphosphine) rhodium (I) has been incorporated into model biomembrane structures in the form of lipid bilayer dispersions in water. This enables the hydrogenation of the double bonds of the unsaturated lipids within the bilayers to be accomplished. To decide the optimum conditions for efficient hydrogenation the reaction conditions have been varied. The effect of catalyst concentration, hydrogen gas pressure and lipid composition (with and without cholesterol) have all been studied. The partition of the catalyst into the lipid medium was checked by rhodium analysis. The results show that an increase of catalyst concentration or an increase of hydrogen gas pressure leads to increasing rates of hydrogenation. Successful hydrogenation was accomplished with different types of lipid dispersions (mitochondrial, microsomal and erythrocyte lipids). A selectivity of the homogeneous hydrogenation process is indicated. The polyunsaturated fatty acyl residues are hydrogenated at an earlier stage and at a faster rate than the monoenoic acids. Furthermore, an increase in the proportion of cholesterol to lipid within the bilayer structures causes a progressive decrease in the rate of hydrogenation. The fluidity of the lipid bilayers can be altered to such an extent by the hydrogenation process that new sharp endotherms corresponding to the order-disorder transition of saturated lipids occur at temperatures as high as 319 K. Some potential uses of hydrogenation for the modulation of cell membrane fluidity are discussed as well as the design of new types of catalyst molecules.     

On the mechanism of non-electrolyte permeation through lipid bilayers and through biomembranes

The temperature effects on the permeation of polyhydroxy alcohols through the lipid bilayers of liposomes with a great variety in chemical composition were studied. Although important differences in the permeability of the various lipid bilayers were observed, Arrhenius plots demonstrated that the activation energy is independent of the degree of unsaturation or the presence of cholesterol in the paraffin barriers. The activation energies found for the penetration of a bilayer with a liquid paraffin core are 14.3 kcal for glycol, 19.4 kcal for glycerol, and 20.8 kcal for erythritol. These values are in agreement with the energies that can be expected for complete dehydration of the permeant molecules. The idea that the activation energy is determined by the number of hydrogen bonds with water is supported by the finding that a series of different diols did demonstrate practically identical activation energies. Studies on a number of biological membranes demonstrated the same activation energies for the penetration of glycerol and erythritol as found in the experiments with liposomes. These facts support the view that both the lipid bilayers and the biological membranes are penetrated by single fully dehydrated molecules.     

Molecular simulation study of the influence of small molecules on the dynamic and structural properties of phospholipid bilayers

Molecular simulations were performed for phospholipid bilayers as model cell membranes to understand their dynamic and structural properties, in particular how those are affected in the presence of the cell-preservation agents trehalose and DMSO. Simulations covered a large range of temperatures and concentrations of trehalose and DMSO. A number of properties were examined, including the area per head group, lipid-tail order parameter, density profile, conformation, and binding of trehalose. Our results indicate that trehalose specifically interacts with the phospholipid bilayer by binding to the polar head groups (H-bonds). This binding serves as a bridge between the lipids and prevents their aggregation. The hydrophilic nature of trehalose also makes it a water substitute when binding to the phospholipids. At high trehalose concentrations, trehalose forms a glassy matrix and preserves the bilayer structure by 'freezing' it in its natured state, as the core of the bilayer is found to remain fluid. Simulations with DMSO revealed that DMSO can penetrate and freely diffuse through the bilayer, which is attributed to its unique characteristic of interacting favorably with both hydrophilic and hydrophobic domains.     

How lipids and proteins interact in a membrane: a molecular approach

Membrane proteins in a biological membrane are surrounded by a shell or annulus of 'solvent' lipid molecules. These lipid molecules in general interact rather non-specifically with the protein molecules, although a few 'hot-spots' may be present on the protein where anionic lipids bind with high affinity. Because of the low structural specificity of most of the annular sites, the composition of the lipid annulus will be rather similar to the bulk lipid composition of the membrane. The structures of the solvent lipid molecules are important in determining the conformational state of a membrane protein, and hence its activity, through charge and hydrogen bonding interactions between the lipid headgroups and residues in the protein, and through hydrophobic matching between the protein and the surrounding lipid bilayer. Evidence is also accumulating for the presence of 'co-factor' lipid molecules binding with high specificity to membrane proteins, often between transmembrane alpha-helices, and often being essential for activity.     

Monotopic enzymes and lipid bilayers: a comparative study

Monotopic proteins make up a class of membrane proteins that bind tightly to, but do not span, cell membranes. We examine and compare how two monotopic proteins, monoamine oxidase B (MAO-B) and cyclooxygenase-2 (COX-2), interact with a phospholipid bilayer using molecular dynamics simulations. Both enzymes form between three and seven hydrogen bonds with the bilayer in our simulations with basic side chains accounting for the majority of these interactions. By analyzing lipid order parameters, we show that, to a first approximation, COX-2 disrupts only the upper leaflet of the bilayer. In contrast, the top and bottom halves of the lipid tails surrounding MAO-B are more and less ordered, respectively, than in the absence of the protein. Finally, we identify which residues are important in binding individual phospholipids by counting the number and type of lipid atoms that come close to each amino acid residue. The existing models that explain how these proteins bind to bilayers were proposed following inspection of the X-ray crystallographic structures. Our results support these models and suggest that basic residues contribute significantly to the binding of these monotopic proteins to bilayers through the formation of hydrogen bonds with phospholipids.     

Molecular dynamics simulation of a palmitoyl-oleoyl phosphatidylserine bilayer with Na+ counterions and NaCl

Two 40 ns molecular dynamics simulations of a palmitoyl-oleoyl phosphatidylserine (POPS) lipid bilayer in the liquid crystalline phase with Na(+) counterions and NaCl were carried out to investigate the structure of the negatively charged lipid bilayer and the effect of salt (NaCl) on the lipid bilayer structure. Na(+) ions were found to penetrate deep into the ester region of the water/lipid interface of the bilayer. Interaction of the Na(+) ions with the lipid bilayer is accompanied by a loss of water molecules around the ion and a simultaneous increase in the number of ester carbonyl oxygen atoms binding the ion, which define an octahedral and square pyramidal geometry. The amine group of the lipid molecule is involved in the formation of inter- and intramolecular hydrogen bonds with the carboxylate and the phosphodiester groups of the lipid molecule. The area per lipid of the POPS bilayer is unaffected by the presence of 0.15M NaCl. There is a small increase in the order parameter of carbon atoms in the beginning of the alkyl chain in the presence of NaCl. This is due to a greater number of Na(+) ions being coordinated by the ester carbonyl oxygen atoms in the water/lipid interface region of the POPS bilayer.