Quantitative determination of conformational disorder in the acyl chains of phospholipid bilayers by infrared spectroscopy

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.

Conformational order of n-alkyl modified silica gels as evaluated by Fourier transform infrared spectroscopy

The conformational behaviour of non-deuterated and selectively deuterated alkyl modified silica gels in the dry state is examined by variable temperature FT IR spectroscopy. In the present study, three systems are considered, which are distinguished by the length of the tethered alkyl chains (C9Hl9-, C18H37-, C22H45-). The desired information is obtained by the analysis of various conformational-sensitive IR bands, including CH2 wagging, CD2 stretching and CD2 rocking bands. The analysis of the CH2 wagging bands provides the relative amounts (i.e., integral numbers over the whole chain) of the kink/gauche-trans-gauche, double-gauche and end-gauche conformers in the tethered alkyl chains. From the analysis of the CD2 stretching and CD2 rocking bands information about the conformational order at a specific deuterated methylene segment is available. Here, the CD2 rocking band data are used to determine the amount of gauche conformers at the deuterated carbon positions C-4 and C-6, and C-12. It is found that the conformational order critically depends on the actual alkyl chain length, chain position and sample temperature. Particular emphasis is given to the impact of the external pressure during sample preparation on the alkyl chain conformations, about which so far no information is available. It is observed that the samples prepared as KBr pellets, which experienced a pressure of about 10 kbar, are characterised by a lower amount of gauche conformers. This substantial increase of conformational order is attributed to better alkyl chain packing along with a gain of intermolecular chain interactions.

Reorganizational dynamics of multilamellar lipid bilayer assemblies using continuously scanning Fourier transform infrared spectroscopic imaging

We employ an implementation of rapid-scan Fourier transform infrared (FT-IR) microspectroscopic imaging to acquire time-resolved images for assessing the non-repetitive reorganizational dynamics of aqueous dispersions of multilamellar lipid vesicles (MLVs) derived from distearoylphosphatidylcholine (DSPC). The spatially and temporally resolved images allow direct and simultaneous determinations of various physical and chemical properties of the MLVs, including the main thermal gel to liquid crystalline phase transition, comparisons of vesicle diffusion rates in both phases and the variation in lipid bilayer packing properties between the inner and outer lamellae defining the vesicle. Specifically, in the lipid liquid crystalline phase, the inner bilayers of the MLVs are more intermolecularly ordered than the outer regions, while the intramolecular acyl chain order/disorder parameters, reflecting the overall characteristics of the fluid phase, remain uniform across the vesicle diameter. In contrast, the lipid vesicle gel phase displays no intermolecular or intramolecular dependence as a function of distance from the MLV center.

Effect of selected environmental and physico-chemical factors on bacterial cytoplasmic membranes

Membranes lipids are one of the most adaptable molecules in response to perturbations. Even subtle changes of the composition of acyl chains or head groups can alter the packing arrangements of lipids within the bilayer. This changes the balance between bilayer and nonbilayer lipids, serving to affect bilayer stability and fluidity, as well as altering lipid-protein interactions. External factors can also change membrane fluidity and lipid composition; including temperature, chemicals, ions, pressure, nutrients and the growth phase of the microbial culture. Various biophysical techniques have been used to monitor fluidity changes within the bacterial membrane. In this review, bacterial cytoplasmic membrane changes and related functional effects will be examined as well as the use of fluorescence polarization methods and examples of data obtained from research with bacteria.

The structure of DNA within cationic lipid/DNA complexes

The structure of DNA within CLDCs used for gene delivery is controversial. Previous studies using CD have been interpreted to indicate that the DNA is converted from normal B to C form in complexes. This investigation reexamines this interpretation using CD of model complexes, FTIR as well as Raman spectroscopy and molecular dynamics simulations to address this issue. CD spectra of supercoiled plasmid DNA undergo a significant loss of rotational strength in the signal near 275 nm upon interaction with either the cationic lipid dimethyldioctadecylammonium bromide or 1,2-dioleoyltrimethylammonium propane. This loss of rotational strength is shown, however, by both FTIR and Raman spectroscopy to occur within the parameters of the B-type conformation. Contributions of absorption flattening and differential scattering to the CD spectra of complexes are unable to account for the observed spectra. Model studies of the CD of complexes prepared from synthetic oligonucleotides of varying length suggest that significant reductions in rotational strength can occur within short stretches of DNA. Furthermore, some alteration in the hydrogen bonding of bases within CLDCs is indicated in the FTIR and Raman spectroscopy results. In addition, alterations in base stacking interactions as well as hydrogen bonding are suggested by molecular dynamics simulations. A global interpretation of all of the data suggests the DNA component of CLDCs remains in a variant B form in which base/base interactions are perturbed.

Constant normal pressure, constant surface tension, and constant temperature molecular dynamics simulation of hydrated 1,2-dilignoceroylphosphatidylcholine monolayer

A constant normal pressure, constant surface tension, and constant temperature (NP(N)gammaT) molecular dynamics (MD) simulation of the liquid condensed phase of a 1,2-dilignoceroylphosphatidylcholine (DLGPC) monolayer has been performed at 293.15 K. A DLGPC molecule has two saturated 24-carbon acyl chains, giving the hydrocarbon core thickness of the monolayer approximately 28 A, which is close to the hydrocarbon core thickness of a membrane of a living system. NP(N)gammaT ensemble was used to reproduce the experimental observations, such as area/lipid, because surface tension is an essential factor in determining the monolayer structure. Data analysis on DLGPC/water monolayer shows that various liquid condensed-phase properties of the monolayer have been well reproduced from the simulation, indicating that surface tension 22.9 mN/M used in the simulation is an appropriate condition for the condensed-phase NP(N)gammaT simulation. The simulation results suggest that this long-chain phospholipid monolayer shares many structural characteristics with typical short-chain 1,2-diacylphosphatidylcholine systems, such as DPPC/water monolayer in the condensed phase and DPPC/water bilayer in the gel phase. Furthermore, it was found that DLGPC/water monolayer has almost completely rotationally disordered acyl chains, which have not been observed so far in short-chain 1,2-diacylphosphatidylcholine/water bilayers. This study indicates the good biological relevance of the DLGPC/water monolayer which might be useful in protein/lipid studies to reveal protein structure and protein/lipid interactions in a membrane environment.

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.

Infrared spectroscopic characterization of the interaction of cationic lipids with plasmid DNA

Fourier transform infrared spectroscopy was used to characterize the interaction of the cationic lipids 1,2-dioleoyl-3-trimethylammonium-propane and dioctadecyldimethylammonium bromide with plasmid DNA. The effect of incorporating the neutral colipids cholesterol and dioleoylphosphatidylethanolamine on this interaction was also examined. Additionally, dynamic and phase analysis light scattering were used to monitor the size and zeta potential of the resulting complexes under conditions similar to the Fourier transform infrared measurements. Results suggest that upon interaction of cationic lipids with DNA, the DNA remains in the B form. Distinct changes in the frequency of several infrared bands arising from the DNA bases, however, suggest perturbation of their hydration upon interaction with cationic lipids. A direct interaction of the lipid ammonium headgroup with and dehydration of the DNA phosphate is observed when DNA is complexed with these lipids. Changes in the apolar regions of the lipid bilayer are minimal, whereas the interfacial regions of the membrane show changes in hydration or molecular packing. Incorporation of helper lipids into the cationic membranes results in increased conformational disorder of the apolar region and further dehydration of the interfacial region. Changes in the hydration of the DNA bases were also observed as the molar ratio of helper lipid in the membranes was increased.

Molecular dynamics simulations of stratum corneum lipid models: fatty acids and cholesterol

We report the results of an investigation on stratum corneum lipids, which present the main barrier of the skin. Molecular dynamics simulations, thermal analysis and FTIR measurements were applied. The primary objective of this work was to study the effect of cholesterol on skin structure and dynamics. Two molecular models were constructed, a free fatty acid bilayer (stearic acid, palmitic acid) and a fatty acid/cholesterol mixture at a 1:1 molar ratio. Our simulations were performed at constant pressure and temperature on a nanosecond time scale. The resulting model structures were characterized by calculating surface areas per headgroup, conformational properties, atom densities and order parameters of the fatty acids. Analysis of the simulations indicates that the free fatty acid fraction of stratum corneum lipids stays in a highly ordered crystalline state at skin temperatures. The phase behavior is strongly influenced when cholesterol is added. Cholesterol smoothes the rigid phases of the fatty acids: the order of the hydrocarbon tails (mainly of the last eight bonds) is reduced, the area per molecule becomes larger, the fraction of trans dihedrals is lower and the hydrophobic thickness is reduced. The simulation results are in good agreement with our experimental data from FTIR analysis and NIR-FT Raman spectroscopy.

A macroscopic description of lipid bilayer phase transitions of mixed-chain phosphatidylcholines: chain-length and chain-asymmetry dependence

A macroscopic model is presented to quantitatively describe lipid bilayer gel to fluid phase transitions. In this model, the Gibbs potential of the lipid bilayer is expressed in terms of a single order parameter q, the average chain orientational order parameter. The Gibbs potential is based on molecular mean-field and statistical mechanical calculations of inter and intrachain interactions. Chain-length and chain-asymmetry are incorporated into the Gibbs potential so that one equation provides an accurate description of mixed-chain phosphatidylcholines of a single class. Two general classes of lipids are studied in this work: lipid bilayers of partially or noninterdigitated gel phases, and bilayers of mixed interdigitated gel phases. The model parameters are obtained by fitting the transition temperature and enthalpy data of phosphatidylcholines to the model. The proposed model provides estimates for the transition temperature and enthalpy, van der Waals energy, number of gauche bonds, chain orientational order parameter, and bond rotational and excluded volume entropies, achieving excellent agreement with existing data obtained with various techniques.

Area per lipid and acyl length distributions in fluid phosphatidylcholines determined by (2)H NMR spectroscopy

Deuterium ((2)H) NMR spectroscopy provides detailed information regarding the structural fluctuations of lipid bilayers, including both the equilibrium properties and dynamics. Experimental (2)H NMR measurements for the homologous series of 1, 2-diacyl-sn-glycero-3-phosphocholines with perdeuterated saturated chains (from C12:0 to C18:0) have been performed on randomly oriented, fully hydrated multilamellar samples. For each lipid, the C-D bond order parameters have been calculated from de-Paked (2)H NMR spectra as a function of temperature. The experimental order parameters were analyzed using a mean-torque potential model for the acyl chain segment distributions, and comparison was made with the conventional diamond lattice approach. Statistical mechanical principles were used to relate the measured order parameters to the lipid bilayer structural parameters: the hydrocarbon thickness and the mean interfacial area per lipid. At fixed temperature, the area decreases with increasing acyl length, indicating increased van der Waals attraction for longer lipid chains. However, the main effect of increasing the acyl chain length is on the hydrocarbon thickness rather than on the area per lipid. Expansion coefficients of the structural parameters are reported and interpreted using an empirical free energy function that describes the force balance in fluid bilayers. At the same absolute temperature, the phosphatidylcholine (PC) series exhibits a universal chain packing profile that differs from that of phosphatidylethanolamines (PE). Hence, the lateral packing of phospholipids is more sensitive to the headgroup methylation than to the acyl chain length. A fit to the area per lipid for the PC series using the empirical free energy function shows that the PE area represents a limiting value for the packing of fluid acyl chains.

Molecular dynamics and (2)H-NMR study of the influence of an amphiphilic peptide on membrane order and dynamics

A molecular dynamics simulation of a fully hydrated model membrane consisting of 12 molecules of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine, one amphiphilic peptide with the sequence acetyl-Lys-Lys-Gly-Leu(16)-Lys-Lys-Ala-amide, and 593 water molecules was performed for 1.06 ns (Belohorcova, K., J. H. Davis, T. B. Woolf, and B. Roux. 1997. Biophys. J. 73:3039-3055). The analysis presented here is primarily focused on the phospholipid component and the results are compared with experimental (2)H-NMR studies of the lipid component of mixtures of the same peptide and lipid at a molar ratio of 1:32, and with earlier studies of closely related peptide/lipid mixtures. The phospholipid chain and headgroup isomer populations and isomerization rates compare favorably with previous simulations and experimental measurements. Of particular interest is the effect of the peptide on the phospholipid headgroup and hydrocarbon chain orientational order calculated from the simulation, which also agree well with experimental measurements performed on this and closely related systems. Comparison of the experimental results with the simulations not only shows that there is significant agreement between the two methods, but also provides new insight into the effect of the peptide on the lipid dynamics. In particular, these results confirm that a membrane spanning peptide has little effect on lipid chain order, and bilayer thickness if its hydrophobic length closely matches the lipid hydrocarbon thickness. In addition, we find that the peptide can have a strong ordering effect if it is longer than the lipid hydrophobic thickness.

Anisotropic motion of cholesterol in oriented DPPC bilayers studied by quasielastic neutron scattering: the liquid-ordered phase

Quasielastic neutron scattering (QENS) at two energy resolutions (1 and 14 microeV) was employed to study high-frequency cholesterol motion in the liquid ordered phase (lo-phase) of oriented multilayers of dipalmitoylphosphatidylcholine at three temperatures: T = 20 degrees C, T = 36 degrees C, and T = 50 degrees C. We studied two orientations of the bilayer stack with respect to the incident neutron beam. This and the two energy resolutions for each orientation allowed us to determine the cholesterol dynamics parallel to the normal of the membrane stack and in the plane of the membrane separately at two different time scales in the GHz range. We find a surprisingly high, model-independent motional anisotropy of cholesterol within the bilayer. The data analysis using explicit models of molecular motion suggests a superposition of two motions of cholesterol: an out-of-plane diffusion of the molecule parallel to the bilayer normal combined with a locally confined motion within the bilayer plane. The rather high amplitude of the out-of-plane diffusion observed at higher temperatures (T >/= 36 degrees C) strongly suggests that cholesterol can move between the opposite leaflets of the bilayer while it remains predominantly confined within its host monolayer at lower temperatures (T = 20 degrees C). The locally confined in-plane cholesterol motion is dominated by discrete, large-angle rotational jumps of the steroid body rather than a quasicontinous rotational diffusion by small angle jumps. We observe a significant increase of the rotational jump rate between T = 20 degrees C and T = 36 degrees C, whereas a further temperature increase to T = 50 degrees C leaves this rate essentially unchanged.

Two British families with variants of the 'cryohydrocytosis' form of hereditary stomatocytosis

We describe two British families with similar, dominantly-inherited, temperature-related variants of hereditary stomatocytosis, consistent with the original description of 'cryohydrocytosis'. The cells show a 5-6-fold increase in passive permeability at 37 degrees C with abnormal intracellular Na and K levels at 15-20 and 60-65 mmol/(l cells) respectively. Marked temperature effects were evident: lysis of red cells on storage in the cold was blatant and when whole heparinized blood was stored at room temperature, K accumulated in the plasma, producing 'pseudohyperkalaemia'. Studies of the temperature dependence of passive permeability showed that the minimum in the passive permeability, which is seen in normal cells at 8-10 degrees C, was shifted up to 23 degrees C in these abnormal cells, such that the permeability at 0 degrees C exceeded that at 37 degrees C. The abnormal temperature dependence in these genetically abnormal red cells strongly resembles that seen in normal cells when suspended in media in which either Na or Cl has been replaced by an organic cation or anion: it could be said these cells had a genetic mutation that somehow rendered the cell resistant to the stabilizing action of NaCl at low temperatures.

Infrared studies of protein-induced perturbation of lipids in lipoproteins and membranes

The paper reviews the main recent publications concerning infrared (IR) spectroscopy as applied to the study of lipid-protein interactions in model and cell membranes, lipoproteins, and related systems (e.g. lung surfactant). The review focuses mainly on transmission IR. Based on the available data, a number of general conclusions are presented on the perturbations caused by proteins on either the hydrocarbon chains, the polar headgroups or the interface region. Lipid-protein interactions in native cell membranes do not reveal significant differences from what is observed in semisynthetic model systems.

Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations

Lipid adhesion forces can be measured using several experimental techniques, but none of these techniques provide insight on the atomic level. Therefore, we performed extensive nonequilibrium molecular dynamics simulations of a phospholipid membrane in the liquid-crystalline phase out of which individual lipid molecules were pulled. In our method, as an idealization of the experimental setups, we have simply attached a harmonic spring to one of the lipid headgroup atoms. Upon retraction of the spring, the force needed to drag the lipid out of the membrane is recorded. By simulating different retraction rates, we were able to investigate the high pull rate part of the dynamical spectrum of lipid adhesion forces. We find that the adhesion force increases along the unbinding path, until the point of rupture is reached. The maximum value of the adhesion force, the rupture force, decreases as the pull rate becomes slower, and eventually enters a friction-dominated regime. The computed bond lengths depend on the rate of rupture, and show some scatter due to the nonequilibrium nature of the experiment. On average, the bond length increases from approximately 1.7 nm to 2.3 nm as the rates go down. Conformational analyses elucidate the detailed mechanism of lipid-membrane bond rupture. We present results of over 15 ns of membrane simulations. Implications for the interpretation and understanding of experimental rupture data are discussed.

Structure and dynamic properties of diunsaturated 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine lipid bilayer from molecular dynamics simulation

Unsaturated fatty acid chains are known to be an essential structural part of biomembranes, but only monounsaturated chains have been included in the molecular dynamics (MD) simulations of membrane systems. Here we present a 1-ns MD simulation for a diunsaturated 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC; 16:0/18:2[delta9,12]) bilayer. The structural behavior of the phosphatidylcholine headgroup, the glycerol backbone, and the hydrating water were assessed and found to be consistent with the existing information about similar systems from both experimental and computational studies. Further analysis was focused on the structure of the double bond region and the effects of the diunsaturation on the bilayer interior. The behavior of the diunsaturated sn-2 chains is affected by the tilted beginning of the chain and the four main conformations of the double bond region. The double bonds of the sn-2 chains also influenced the characteristics of the saturated chains in the sn-1 position. Furthermore, extreme conformations of the sn-2 chains existed that are likely to be related to the functional role of the double bonds. The results here point out the importance of polyunsaturation for the biological interpretations deduced from the membrane MD simulations.

Molecular dynamics simulations of six different fully hydrated monomeric conformers of Escherichia coli re-lipopolysaccharide in the presence and absence of Ca2+

Six previously published conformational models of Escherichia coli Re lipopolysaccharide (ReLPS) were subjected to molecular dynamics simulations using the CHARMM force field. The monomers of ReLPS were completely immersed in a water box. The dynamic behavior of the solvated models in the presence and absence of calcium cations was compared. The structure of the solvent shell was analyzed in terms of radial distribution functions. Diffusion coefficients and mean residence times were analyzed to characterize the dynamic behavior of the solvent. Order parameters and number of gauche defects were used for the description of the dynamics of the acyl chains. The cations are preferentially located between the carboxylate and phosphate groups of the headgroup. Their presence leads to a rigidification of the headgroup structure and alters the conformation of the backbone, thus influencing the structure and flexibility of the hydrophobic region as well. The effect of calcium on the backbone flexibility was measured in terms of glycosidic torsion angles. The six fatty acid chains of each ReLPS monomer adopt a highly ordered micromembrane structure. The packing parameter indicates that aggregation of these ReLPS monomers will lead to lamellar structures. Evaluation of all data enables us to present one conformation, C, which is thought to best represent the average structure of the ReLPS conformers.

Interdigitation of bilayers from ether lipid analogs: (R)-PAF, (R)-Lyso-PAF and the antineoplastic (R)-ET-18-OMe

The ether lipid analogs (R)-PAF, (R)-ET-18-OMe and (R)-Lyso-PAF can form bilayers. Using differential scanning calorimetry (DSC) and small angle X-ray diffraction, we have investigated bilayers formed entirely of these molecules. Small angle X-ray diffraction experiments strongly indicated that these bilayers are interdigitated at 98% relative humidity in the gel phase. The slight differences in the structures of the ether lipids affect their thermal behavior and packing properties. Of the three, (R)-Lyso-PAF shows the highest phase transition temperature, broadest phase transition peak and smallest enthalpy change. These characteristics may be attributed to hydrogen bonding of the -OH group at sn-2 position of its glycerol backbone to water or to other (R)-Lyso-PAF molecules. The interdigitated structure shows that these bilayers are similar to lipid bilayers and suggests that these molecules could easily insert themselves into membranes to reach their respective active sites or to perturb the membranes. The packing and thermal characteristics of (R)-Lyso-PAF may help explain its biological inactivity.

The effect of high external pressure on DPPC-cholesterol multilamellar vesicles: a pressure-tuning Fourier transform infrared spectroscopy study

We have investigated the effect of incorporation of cholesterol on the barotropic phase behavior of aqueous dispersions of 1,2-dipalmitoylphosphatidylcholine (DPPC) using Fourier-transform infrared spectroscopy (FTIR) in combination with the diamond anvil technique. Infrared spectral parameters, such as the frequencies, intensities, bandshapes and band splittings have been used to detect structural and dynamical changes upon incorporation of cholesterol into the DPPC bilayer. Analysis of these spectral parameters yields information on conformer population, reorientational fluctuations, interchain interaction, hydrogen bonding, interdigitation packing, and phase transformations of the DPPC/cholesterol mixtures. We present FTIR data of aqueous DPPC dispersions at 0, 10, 20, 30, 40, and 50 mol% cholesterol in the pressure range from 0.001 to 20 kbar at two temperatures, 25 degrees C and 55 degrees C. In addition, comprehensive temperature dependent measurements in the range from 20 degrees C to 80 degrees C were performed at ambient pressure. Analysis of the CH2 symmetric and antisymmetric stretching modes yields information of the effect of cholesterol concentration on the phase transition phenomena occurring in the lipid bilayer. Observation of the correlation field splittings of the CH2 bending and rocking modes monitors structural changes and dynamical properties of the lipid mixtures. Cholesterol induces more orientational disorder of the lipid molecules in terms of an increase of the reorientational fluctuations of the molecules and twisting/torsion motions of the acyl chains in the gel phase even at elevated pressures. It therefore appears that one important role of cholesterol is to make the membrane insensitive to changes in external environment, such as high hydrostatic pressure. Increase of pressure leads to a decrease in half width of the C = O band contour of pure DPPC and of DPPC/cholesterol mixtures, especially for cholesterol concentrations equal and higher than 30 mol%, which might be due to a marked increase in free carbonyl groups. At high pressure part of the bound water from the interfacial zone of the membrane is withdrawn. Increase of the cholesterol concentration and increase in pressure have opposite effects on the population of free and hydrated carbonyl ester groups of DPPC in the gel phases.

Molecular dynamics simulations of adipocyte lipid-binding protein: effect of electrostatics and acyl chain unsaturation

Molecular dynamics (MD) simulations have been performed on adipocyte lipid-binding protein, using the apo and holo forms, bound with stearic and oleic acid. The contribution of electrostatics to protein dynamics and ligand stabilization was assayed by perturbing the electrostatic charge of Arg106 and Arg126 (positive-->neutral) and the fatty acid (132H) headgroup (negative-->neutral). MD simulations for charged holo forms demonstrated significantly greater electrostatic binding energy and a more stabilized hydrogen bond network than simulations performed using neutral forms. Electrostatics, however, appeared to have little effect on fatty acid behavior, e.g., fluctuation of the dihedral head group; number of dihedral transitions within the acyl chain; and change in the end-to-end distance for fatty acid. Instead, fatty acid behavior appeared to be dictated by the presence or absence of an unsaturated bond within the acyl chain. A significantly greater number of transitions were observed during MD simulations in oleic than stearic acid. In addition, significantly greater fluctuation was observed for oleic acid, within the C2 headgroup and C9 and C11 dihedrals (which lie adjacent to the olefin bond of oleic acid). The dynamic behavior of the acyl chain may thereby be more a property of van der Waals contact, and the degree of acyl chain unsaturation, than a function of electrostatics. In the absence of fatty acid, an increase in distance between guanidino carbon centered atoms of Arg126 and Arg106 was observed during MD simulations of the charged apo form. This effect not observed with the neutral apo form or in any of the holo complexes and, presumbably, was a result of repulsion between the negatively charged arginine sidechains. Conserved waters reflected substantially lower mean-square displacement (msd) in all simulations, except the neutral apo form. This suggests that the presence of either charged amino acids or lipid provides increased order for water within the binding pocket. These results provide a dynamic perspective of the interactive nature within the FABP binding pocket regulated in a complex manner by the electrostatics within the binding cavity, acyl chain structure and behavior, and water energetics.