Temperature change of the ripple structure in fully hydrated dimyristoylphosphatidylcholine/cholesterol multibilayers.

The ripple structure was studied as a function of temperature in fully hydrated dimyristoylphosphatidylcholine (DMPC)/cholesterol multibilayers using synchrotron x-ray small-angle diffraction and freeze-fracture electron microscopy. In the presence of cholesterol, the ripple structure appears below the pretransition temperature of pure DMPC multibilayers. In this temperature range the ripple periodicity is relatively large (25-30 nm) and rapidly decreases with increasing temperature. In this region, defined as region I, we observed coexistence of the P beta' phase and the L beta' phase. The large ripple periodicity is caused by the formation of the P beta' phase region in which cholesterol is concentrated and the L beta' phase region from which cholesterol is excluded. An increase in ripple periodicity also takes place in the narrow temperature range just below the main transition temperature. We define this temperature region as region III, where the ripple periodicity increases dramatically toward the main transition temperature. In region II, between regions I and III, the ripple periodicity decreases gradually with temperature. This behavior is quite similar to that of pure DMPC. Temperature-versus-ripple periodicity curves are parallel among pure DMPC and DMPCs with various cholesterol contents. We explain this behavior in terms of a model proposed by other workers.     

Thermotropic changes in the conductivity of black bilayers from egg phosphatidylethanolamine

Current fluctuations in black bilayers from phosphatidyl ethanolamine obtained from egg lecithin were registered in the temperature region of the main phase transition of this phospholipid and the bilayer--hexagonal phase transition about 35 degrees; they correspond to the conductivity changes of hundreds of pSm. This transition takes place in the same temperature region as shown by 31P-NMR and depolarization of light-scattering method in phosphatidyl ethanolamine multilamellar liposomes. The scheme of bilayer transformation into hexagonal phase in the temperature region of lipid polymorphic transition is discussed.     

Indirect evidence for lipid-domain formation in the transition region of phospholipid bilayers by two-probe fluorescence energy transfer.

The fluorescence energy transfer between two lipid probes, N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)-1, 2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine (donor) and N-(Lissamine rhodamine B sulfonyl)-1, 2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine (acceptor), incorporated into 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine unilamellar and multilamellar lipid bilayers, is studied in the temperature region of the main phase transition. The two probes display different relative solubilities in the gel and fluid lipid-bilayer phases. A distinct maximum in the fluorescence intensity of the donor is observed in the transition region, indicating that the two probes are demixing and hence increasing their average separation. The observation is interpreted in terms of dynamic segregation of the two probes into coexisting gel and fluid lipid domains that are formed dynamically in the transition region due to strong density fluctuations. The interpretation of the experimental observations is supported by a detailed theoretical calculation using computer simulation of a microscopic model that takes full account of diffusion of the two probes and the fluctuations of gel and fluid lipid domains characteristic of the main phase transition.     

Redetermination of the pressure dependence of the lipid bilayer phase transition.

The effect of pressure on the phase transition temperature for the dipalmitoyllecithin bilayer was redetermined by following the volume change accompanying the transition. These measurements were carried out isothermally with the transition from the ordered to the disordered phase induced by decreasing the pressure. This contrasts with our previous measurements which were carried out at constant pressure and increasing temperature. The transition at every temperature was sharp and confirmed our previous observation that the volume change associated with the transition (0.033 mL g-1) is invariant with pressure. However, our present measurements, in contrast to our previous results, indicate that dP m/dTm at all pressures is in agreement with the 1 atm value of delta H/Tm delta V within experimental error where Tm and Pm are the temperature and pressure of the phase transition, respectively. These results, which are now in agreement with all other known pressure data, indicate that the entropy change associated with the transition is invariant with pressure.     

Properties of bilayer membranes in the phase transition or phase separation region

The increase in passive permeability of bilayer membranes near the phase transition temperature is usually explained as caused by either the increase in the amount of ‘boundary lipid’ present in the membrane, or by the increase in lateral compressibility of the membrane. Since both the amount of ‘boundary lipid’ and the lateral compressibility show a similar anomaly near the transition temperature, it is difficult to distinguish experimentally between the two proposed mechanisms.We have examined some details of both of the proposed pictures. The fluid-solid boundary energy, neglected in previous work, has been computed as a function of the domain size. For a single component uncharged lipid bilayer, the results rule out the existence of even loosely defined solid domains in a fluid phase, or vice versa. Thermodynamic fluctuations, which are responsible for anomalous behaviour near the phase transition temperature, are not intense enough to approximate the formation of a domain of the opposite phase.Turning next to lateral compressibility of bilayer membranes we have considered two-component mixtures in the phase separation region. We present the first calculation of lateral compressibility for such systems. The behaviour shows interesting anomalies, which should correlate with existing and future data on transport across membranes.     

Phase diagrams for impure lipid systems. Application to lipid/anaesthetic mixtures

A physical model is presented to describe theoretically the temperature-dependent interactions of lipid bilayers with small molecules such as anaesthetics. Based on an earlier model, a triangular lattice in which each site is occupied by a single lipid chain is constructed and the small (anaesthetic) molecules are assumed to occupy interstitial sites in the centre of each lattice triangle. The phase characteristics of such lipid/anaesthetic mixtures are described in terms of the interaction parameters between lipid-lipid, lipid-anaesthetic and anaesthetic-anaesthetic molecules. Depending on the chemical nature of the interacting species the following three models are formulated: Model I. An interstitial model in which the only perturbation is in the head-group region of the bilayer and direct interactions between neighbouring anaesthetic molecules are taken into account. Model II. Here, only hydrophobic interactions between anaesthetics and lipids are considered. Model III. Both van der Waals' and coulombic interactions are taken into account. Phase diagrams for the three models are obtained by numerical calculation over a wide range of interaction parameters. It is shown that in all three models, lateral phase separation takes place due to the presence of anaesthetics. The heat of transition, however, is found to be virtually independent of the anaesthetic concentration.     

Factors determining pressure perturbation calorimetry measurements: evidence for the formation of metastable states at lipid phase transitions

The factors that influence the application of pressure perturbation calorimetry in studying the volume change of the phase transition of lipids are discussed. These factors include a correction for the temperature-shift induced by perturbation, the kinetic irreversibility of the phase transition and the magnitude of the pressure perturbation. We take into account the fact that the dependence of the phase transition temperature on pressure will affect the temperature-shift induced by pressure. As a result, there is a discrepancy between the compression part of the cycle and the expansion. In addition, sequential cycles lead to a gradual loss in magnitude of the heat effect upon pressure perturbation. We suggest that these phenomena can be explained by the formation of a metastable glass-like state that converts to a stable phase at temperatures removed from the region of the phase transition.     

Intrinsic molecules in lipid membranes change the lipid-domain interfacial area: cholesterol at domain interfaces

A theoretical analysis of the effects of intrinsic molecules on the lateral density fluctuations in lipid bilayer membranes is carried out by means of computer simulations on a microscopic interaction model of the gel-to-fluid chain-melting phase transition. The inhomogeneous equilibrium structures of gel and fluid domains, which in previous work (Cruzeiro-Hansson, L. and Mouritsen, O.G. (1988) Biochim. Biophys. Acta 944, 63-72) were shown to characterize the transition region of pure lipid membranes, are here shown to be enhanced by intrinsic molecules such as cholesterol. Cholesterol is found to increase the interfacial area and to accumulate in the interfaces. The interfacial area, the average cluster size, the lateral compressibility, and the membrane area are calculated as functions of temperature and cholesterol concentration. It is shown that the enhancement by cholesterol of the lateral density fluctuations and the lipid-domain interfacial area is most pronounced away from the transition temperature. The implications of the results are discussed in relation to passive ion permeability and function of interfacially active enzymes such as phospholipase.     

Quantification of the effects of chain length and concentration on the thermal behavior of elastin-like polypeptides

At a specific temperature, elastin-like polypeptides (ELPs) undergo a sharp solubility transition that can be exploited in a variety of applications in biotechnology and medicine. The temperature of the transition varies with ELP sequence, molecular weight, and concentration. We present a single equation of three parameters that quantitatively predicts the transition temperature as a function of ELP length and concentration for an ELP of a fixed composition. This model should be useful both for the design of new ELP sequences that exhibit a desired transition temperature and for the selection of variables to trigger the phase transition of an ELP for a given application.     

Simultaneous modeling of phase and calorimetric behavior in an amphiphilic peptide/phospholipid model membrane

Differential scanning calorimetry (DSC) experiments have been performed on the amphiphilic peptide/1,2-bis(perdeuteriopalmitoyl)-sn-glycero-3-phosphocholine system for which partial phase diagrams have been measured by deuterium nuclear magnetic resonance. The solute concentration dependence of the transition enthalpy in such systems is often interpreted in terms of an annulus of lipid withdrawn, by the solvent, from participation in the transition while the bulk lipid melts with its fully enthalpy. This idea is equivalent to postulating ideal mixing between the lipid and the peptide/lipid complex, and there is little justification for such an assumption. Adaptation of regular solution theory to this system demonstrates that the peptide concentration dependence of the transition enthalpies can be incorporated into a thermodynamic model which reproduces the observed phase behavior fairly well without postulating that a complexing annulus of lipid around the peptide be withdrawn from participating in the chain-melting transition. The model parameters determined by simultaneous fitting of the phase behavior and transition enthalpies are used to simulate the DSC scan shapes. The asymmetry of the calorimetric scans for chi 2 less than or equal to 0.02 is reproduced by the model, but a broad component observed for higher concentration is not. In light of the results presented here, previous analyses of the calorimetric behavior of two-component systems in terms of symmetric transitions which do not account for the possible extent of a region of two-phase equilibrium must be questioned.     

Biophysical perturbations induced by ethylazinphos in lipid membranes

Perturbations induced by ethylazinphos on the physical organization of dipalmitoylphosphatidylcholine (DPPC) and DPPC/cholesterol membranes were studied by differential scanning calorimetry (DSC) and fluorescence polarization of 2-, 6-, 12-(9-anthroyloxy) stearic acids and 16-(9-anthroyloxy) palmitic acid. Ethylazinphos (50 and 100 microM) increases the fluorescence polarization of the probes, either in the gel or in the fluid phase of DPPC bilayers, and this concentration dependent effect decreases from the surface to the bilayer core. Additionally, the insecticide displaces the phase transition to a lower temperature range and broadens the transition profile of DPPC. A shifting and broadening of the phase transition is also observed by DSC. Furthermore at insecticide/lipid molar ratios higher than 1/7, DSC thermograms, in addition to the normal transition centered at 41 degrees C, also display a new phase transition centered at 45.5 degrees C. The enthalpy of this new transition increases with insecticide concentration, with a corresponding decrease of the main transition enthalpy. Ethylazinphos in DPPC bilayers with low cholesterol (< or = 20 mol%) perturbs the membrane organization as described above for pure DPPC. However, cholesterol concentrations higher than 20 mol% prevent insecticide interaction, as revealed by fluorescence polarization and DSC data. Apparently, cholesterol significantly modulates insecticide interaction by competition for similar distribution domains in the membrane. The present results strongly support our previous hypothesis that ethylazinphos locates in the cooperativity region, i.e. the region of C1-C9 atoms of the acyl chains, and extends to the lipid-water interface, where it increases lipid packing order sensed across all the thickness of the bilayer. Additionally, and, on the basis of DSC data, a lateral regionalization of ethylazinphos is here tentatively suggested.     

Consideration of the Possibility that the slow step in protein denaturation reactions is due to cis-trans isomerism of proline residues.

A model is proposed to account for the observation that the denaturation of small proteins apparently occurs in two kinetic phases. It is suggested that only one of these phases--the fast one--is actually an unfolding process. The slow phase is assumed to arise from the cis-trans isomerism of proline residues in the denaturated protein. From model compound data, it is shown that the expected rate for isomerism is in satisfactory agreement with the rates actually observed for protein folding. It is also shown that a simple model of protein unfolding based on the isomerism concept is very successful in accounting for many known experimental characteristics of the kinetics and thermodynamic of protein denaturation. Thus, the model is able to predict that two kinetic phases will be seen in the transition region while none are seen in the base-line regions, that both the fast and slow refolding phases lead to the native protein as the product, that the fast phase becomes the only observable phase for jumps ending far in the denatured base-line region, that most or all small proteins show a limiting low-temperature activation energy of ca. 20,000 cal, and that the relaxtion time for the slow phase seen in cytochrome c denaturation is much shorter than for all other small proteins. By utilizing "double-jump" experiments, it is shown directly that the slow phase is not part of the unfolding process but that it corresponds to a transition among two or more denatured forms which have identical spectroscopic (286.5 nm) properties. Thus, the slow relaxation is "invisible" except in the transition region where it couples to the fast unfolding equilibrium. Finally, since the present model assumes that only one of the major kinetic phases seen in denaturation reactions is concerned with the denaturation process per se, it is in agreement with numerous thermodynamic studies which show consistency with the two-state model for unfolding.     

The diversity of the liquid ordered (Lo) phase of phosphatidylcholine/cholesterol membranes: a variable temperature multinuclear solid-state NMR and x-ray diffraction study

To investigate the properties of a pure liquid ordered (Lo) phase in a model membrane system, a series of saturated phosphatidylcholines combined with cholesterol were examined by variable temperature multinuclear (1H, 2H, 13C, 31P) solid-state NMR spectroscopy and x-ray scattering. Compositions with cholesterol concentrations>or=40 mol %, well within the Lo phase region, are shown to exhibit changes in properties as a function of temperature and cholesterol content. The 2H-NMR data of both cholesterol and phospholipids were used to more accurately map the Lo phase boundary. It has been established that the gel-Lo phase coexistence extends to 60 mol % cholesterol and a modified phase diagram is presented. Combined 1H-, 2H-, 13C-NMR, and x-ray scattering data indicate that there are large changes within the Lo phase region, in particular, 1H-magic angle spinning NMR and wide-angle x-ray scattering were used to examine the in-plane intermolecular spacing, which approaches that of a fluid Lalpha phase at high temperature and high cholesterol concentrations. Although it is well known for cholesterol to broaden the gel-to-fluid transition temperature, we have observed, from the 13C magic angle spinning NMR data, that the glycerol region can still undergo a "melting", though this is broadened with increasing cholesterol content and changes with phospholipid chain length. Also from 2H-NMR order parameter data it was observed that the effect of temperature on chain length became smaller with increasing cholesterol content. Finally, from the cholesterol order parameter, it has been previously suggested that it is possible to determine the degree to which cholesterol associates with different phospholipids. However, we have found that by taking into account the relative temperature above the phase boundary this relationship may not be correct.     

Infrared study of phospholipid hydration. New thermodynamic data about the main phase transition of saturated phosphatidylcholine water multidispersions

Infrared spectra measurements are performed with fully hydrated phospholipid multilamellar dispersions. The temperature profiles of the delta (OH) bending vibration at 1645 cm-1, corresponding to the water molecules, show three marked transitions in such temperature ranges that they can be respectively associated with the fusion of free water molecules, with pretransition phenomena and with the main phase transition. From absorbance variations and calorimetric data, main phase transition enthalpies and entropies were calculated for both the acyl chains and the polar head-water regions. Negative contributions of the polar head-water region are consistent with previous results, indicating a more ordered water layer in the fluid phase. On the other hand, a decrease of the amount of bound water molecules with increasing chain length is pointed out.     

ESR studies of the dynamics of dipalmitoylglycerophosphocholine-sulfatide model membranes

The effects of sulfatide on the fluidity and surface dynamics of bilayered and micellar model membranes of dipalmitoylglycerophosphocholine containing sulfatide were studied as a function of gel-to-liquid-crystalline state of the lipids by electron spin resonance. 5- and 15-nitroxystearic acid were employed as spinlabel probes for the region close to the surface and that close to the nonpolar core of lipid structures. The sulfatide effect is completely different above and below the gel-to-liquid-crystalline phase transition point, the glycolipid promoting a more disordered state below it and having a condensing effect above the phase transition temperature.     

A field test of a model for the stability of androdioecy in the freshwater shrimp,Eulimnadia texana

The evolution of hermaphroditism from dioecy is a poorly studied transition. Androdioecy (the coexistence of males and hermaphrodites) has been suggested as an intermediate step in this evolutionary transition or could be a stable reproductive mode. Freshwater crustaceans in the genus Eulimnadia have reproduced via androdioecy for 24+ million years and thus are excellent organisms to test models of the stability of androdioecy. Two related models that allow for the stable maintenance of males and hermaphrodites rely on the counterbalancing of three life history parameters. We tested these models in the field over three field seasons and compared the results to previous laboratory estimates of these three parameters. Male and hermaphroditic ratios within years were not well predicted using either the simpler original model or a version of this model updated to account for differences between hermaphroditic types (‘monogenic’ and ‘amphigenic’ hermaphrodites). Using parameter estimates of the previous year to predict the next year's sex ratios revealed a much better fit to the original relative to the updated version of the model. Therefore, counter to expectations, accounting for differences between the two hermaphroditic types did not improve the fit of these models. At the moment, we lack strong evidence that the long‐term maintenance of androdioecy in these crustaceans is the result of a balancing of life history parameters; other factors, such as metapopulation dynamics or evolutionary constraints, may better explain the 24+ million year maintenance of androdioecy in clam shrimp.     

Applying the stochastic difference equation to DNA conformational transitions: a study of B-Z and B-A DNA transitions

Despite the existence of numerous models to account for the B-Z DNA transition, experimenters have not yet arrived at a conclusive answer to the structural and dynamical features of the B-Z transition. By applying the stochastic difference equation to simulate the B-Z DNA transition, we have shown that the stretched intermediate model of the B-Z transition is more probable than other B-Z transition models such as the Harvey model. This is accomplished by comparing potential energy profiles of various B-Z DNA transition models and calculating relative probabilities based on the stochastic difference equation with respect to length (SDEL) formalism. The results garnered in this article allow for new approaches in determining the structural transition of B-DNA to Z-DNA experimentally. We have also simulated the B-A DNA transition using the stochastic difference equation. Unlike the B-Z DNA transition, the mechanism for the B-A DNA transition is well established. The variation in the pseudorotation angle during the transition is in good agreement with experimental results. Qualitative features of the simulated B-A transition also agree well with experimental data. The SDEL approach is thus a suitable numerical technique to compute long-time molecular dynamics trajectory for DNA molecules.     

Polarizability rescaling and atom-based Thole scaling in the CHARMM Drude polarizable force field for ethers

Within the CHARMM polarizable force field based on the classical Drude oscillator, atomic polarizabilities are derived via fitting to ab initio calculated data on isolated gas phase molecules, with an empirical scaling factor applied to account for differences between the gas and condensed phases. In the development of polarizable models for the ethers, a polarizability scaling factor of 0.7 was previously applied [Vorobyov et al. J Comput Chem 3:1120–1133, 2007]. While the resulting force field models gave good agreement with a variety of experimental data, they systematically underestimated the liquid phase dielectric constants. Here, a new CHARMM polarizable model is developed for the ethers, employing a polarizability scaling factor of 0.85 and including atom-based Thole scale factors recently introduced into the CHARMM Drude polarizable force field [Harder et al. J Phys Chem B 112:3509-3521, 2008]. The new model offers a significant improvement in the reproduction of liquid phase dielectric constants, while maintaining the good agreement of the previous model with all other experimental and quantum mechanical data, highlighting the sensitivity of liquid phase properties to the choice of atomic polarizability parameters.     

Nonsymmetrical polymethine dyes as fluorescent probes for the study of microviscosity of membrane phospholipid bilayers

Nonsymmetrical polymethine dyes are shown to be applied as a new approach in the studies of phospholipid membrane microviscosity. The method requires determination of the intensity ratio for the long-wave length (Ig) and short-wave length (Ik) bands of fluorescence spectra in the region of 730-770 nm at exitation 600 nm. It allows determination of microviscosity by comparative measurements of the fluorescence parameter Ig/Ik in the model medium of the known viscosity (glycerol) and the object under study. Microviscosity in egg phosphatidylcholine vesicules (liposomes) is found to be 0.6-1.2 P. It depends on the surface curvature (size of vesicle), cholesterol content and temperature. It the studies of dimiristoylphosphatidyl choline liposomes the changes in microviscosity at the phase transition temperature are shown to be from 4.5 to 1.1 P. The transition temperature is 24.5 degrees C, the transition range being 2.2 degrees C. The results of this work demonstrate the advantages of the suggested approach to study mobility in phospholipid membranes and confirm it to be promising to study natural membranes and whole cells.     

Simulating the Entropic Collapse of Coarse-Grained Chromosomes

Depletion forces play a role in the compaction and decompaction of chromosomal material in simple cells, but it has remained debatable whether they are sufficient to account for chromosomal collapse. We present coarse-grained molecular dynamics simulations, which reveal that depletion-induced attraction is sufficient to cause the collapse of a flexible chain of large structural monomers immersed in a bath of smaller depletants. These simulations use an explicit coarse-grained computational model that treats both the supercoiled DNA structural monomers and the smaller protein crowding agents as combinatorial, truncated Lennard-Jones spheres. By presenting a simple theoretical model, we quantitatively cast the action of depletants on supercoiled bacterial DNA as an effective solvent quality. The rapid collapse of the simulated flexible chromosome at the predicted volume fraction of depletants is a continuous phase transition. Additional physical effects to such simple chromosome models, such as enthalpic interactions between structural monomers or chain rigidity, are required if the collapse is to be a first-order phase transition.