Thermal Phase Behavior of DMPG Bilayers in Aqueous Dispersions as Revealed by (2)H- and (31)P-NMR

The synthetic lipid 1,2-dimyristoyl-sn-3-phosphoglycerol (DMPG), when dispersed in water/NaCl exhibits a complex phase behavior caused by its almost unlimited swelling in excess water. Using deuterium ((2)H)- and phosphorus ((31)P)-NMR we have studied the molecular properties of DMPG/water/NaCl dispersions as a function of lipid and NaCl concentration. We have measured the order profile of the hydrophobic part of the lipid bilayer with deuterated DMPG while the orientation of the phosphoglycerol headgroup was deduced from the (31)P NMR chemical shielding anisotropy. At temperatures > 30 degrees C we observe well-resolved (2)H- and (31)P NMR spectra not much different from other liquid crystalline bilayers. From the order profiles it is possible to deduce the average length of the flexible fatty acyl chain. Unusual spectra are obtained in the temperature interval of 20-25 degrees C, indicating one or several phase transitions. The most dramatic changes are seen at low lipid concentration and low ionic strength. Under these conditions and at 25 degrees C, the phosphoglycerol headgroup rotates into the hydrocarbon layer and the hydrocarbon chains show larger flexing motions than at higher temperatures. The orientation of the phosphoglycerol headgroup depends on the bilayer surface charge and correlates with the degree of dissociation of DMPG-Na(+). The larger the negative surface charge, the more the headgroup rotates toward the nonpolar region.

Lipid bilayer pre-transition as the beginning of the melting process

We investigate the bilayer pre-transition exhibited by some lipids at temperatures below their main phase transition, and which is generally associated to the formation of periodic ripples in the membrane. Experimentally we focus on the anionic lipid dipalmytoylphosphatidylglycerol (DPPG) at different ionic strengths, and on the neutral lipid dipalmytoylphosphatidylcholine (DPPC). From the analysis of differential scanning calorimetry traces of the two lipids we find that both pre- and main transitions are part of the same melting process. Electron spin resonance of spin labels and excitation generalized polarization of Laurdan reveal the coexistence of gel and fluid domains at temperatures between the pre- and main transitions of both lipids, reinforcing the first finding. Also, the melting process of DPPG at low ionic strength is found to be less cooperative than that of DPPC. From the theoretical side, we introduce a statistical model in which a next-nearest-neighbor competing interaction is added to the usual two-state model. For the first time, modulated phases (ordered and disordered lipids periodically aligned) emerge between the gel and fluid phases as a natural consequence of the competition between lipid-lipid interactions. (C) 2009 Elsevier B.V. All rights reserved.

Measurement of repulsive forces between charged phospholipid bilayers

Electrostatic forces between membranes containing charged lipids were assumed to play an important role in influencing interactions between membranes long before quantitative measurements of such forces were available. ~ur measurements were designed to measure electrostatic forces between layers of lecithin charged with lipi~s carrying ionizable head groups. These experiments have shown that the interactions between charged lipid bila.yere are dominated by electrostatic forces only at separations greater than 30 A. At smaller separations the repulsion between charged bilayers is dominated by strong hydration forces. The net repulsive force between egg lecithin bilayers containing various amounts of cherged lipids (phosphatidylglycerol (PG) 5,10 ano 50 mole%, phosphatidyli. nosi tol (PI) 10 mole% and sodium oleate (Na-Ol) 3,5 and 10 mole%, where mole% gives the ratio of the number of moles' of .charged lipid to the total number of moles of all lipids present in the sample) was stuoied with the help ('If the osmotic streas technique described by LeNeveu et aI, (1977). Also, the forces between pure PG were j_nvestigated in the same manner. The results have been plotted showing variation of force as a function of bilay- _ er separation dw• All curVes 90 obtained called force curves, were found to be similar in sha.pe, showing two distinct regions, one when dw<.30 A is a region cf very rapid iiivariation of force with separation ( it is the region dominated by hydre,tion force) and second when dw> 40 A is a region of very slow variation of force with separB.tion ( it is the region dominated by the electrostatic force). Between these two regions there exists a transition area in which, in most systems studied, a phase separation of lipids into fractions containing different amounts of charged groups, was observed. A qualitative analysis showed that our results were v/ell described by the simple electrostatic double -le.yer theory. For quantitative agreement between measured and calculated force curves however, the charge density for the calculations had to be taken as half of that given by the number density of charged lipids present in the lecithin bilayers. It is not clear at the moment what causes such low apparent degree of ionization among the charged head groups, and further study is needed in this area.

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Surface instability of confined elastic bilayers: Theory and simulations

The surface of a soft elastic film becomes unstable and forms a self-organized undulating pattern because of adhesive interactions when it comes in contact proximity with a rigid surface. For a single film, the pattern length scale lambda, which is governed by the minimization of the elastic stored energy, gives lambda similar to 3h, where h is the film thickness. Based on a linear stability analysis and simulations of adhesion and debonding, we consider the contact instability of an elastic bilayer, which provides greater flexibility in the morphological control of interfacial instability. Unlike the case of a single film, the morphology of the contact instability patterns, debonding distance, and debonding force in a bilayer can be controlled in a nonlinear way by varying the thicknesses and shear moduli of the films. Interestingly, the pattern wavelength in a bilayer can be greatly increased or decreased compared to a single film when the adhesive contact is formed by the stiffer or the softer of the two films, respectively. In particular, lambda as small as 0.5h can be obtained. This indicates a new strategy for pattern miniaturization in elastic contact lithography.

Observation of reduced activation energy and the possible existence of decoupled pancake vortices in superconductor/ferromagnet bilayers

Investigations of different superconducting (S)/ferromagnetic (F) heterostructures grown by pulsed laser deposition reveal that the activation energy (U) for the vortex motion in a high T-c superconductor is reduced remarkably by the presence of F layers. The U exhibits a logarithmic dependence on the applied magnetic field in the S/F bilayers suggesting the existence of decoupled two-dimensional (2D) pancake vortices. This result is discussed in terms of the reduction in the effective S layer thickness and the weakening of the S coherence length due to the presence of F layers. In addition, the U and the superconducting T-c in YBa2Cu3O7-delta/La0.5Sr0.5CoO3 bilayers are observed to be much lower than in the YBa2Cu3O7-delta/La0.7Sr0.3MnO3 ones. This in turn suggests that the degree of spin polarization of the F layer might not play a crucial role for the suppression of superconductivity due to a spin polarized induced pair-breaking effect in S/F bilayers.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L-beta phase to the L-alpha phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT4 bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in kappa from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results. (C) 2010 American Institute of Physics.

Microscopic origin of the chirality driven morphologies of the amphiphilic monolayers and bilayers

It is widely known that the compressed monolayers and bilayers of chiral lipids or fatty acids form helical morphologies, while the corresponding racemic modification gives only flat platelets without twist. No molecular explanation of this phenomenon is yet available, although subtle interactions at the chiral centers have often been proposed as the driving force behind the morphology of the aggregate to form a particular shape. In the present study, the morphologies of the chiral amphiphilic assemblies have been predicted on the basis of an effective pair potential between the molecules, which depends on the relative sizes of the groups attached to the chiral centers, the orientation of the amphiphilic molecules and also on the distance between them. It is shown that fur a pair of same kind of enantiomers, the minimum energy conformation favours a twist angle between them. This twist between the neighbouring molecules gives rise to the helicity of the aggregate. The present theory also shows from the molecular considerations that for a pair of mirror-image isomers (i.e. the racemic modification) the minimum energy conformation corresponds to the zero angle between the molecules, thus giving rise to flat platelets as observed in experiments. Another fascinating aspect of such chirality driven helical structures is that the sense (or the handedness) of the helix is highly specific about the chirality of the monomer concerned. The molecular theory shows, for the first time, that the sense of the helical structures in many cases is determined by the sizes of the groups attached to the chiral centers and the effective potential between them. The predicted senses of the helical structures are in complete agreement with the experimental results.

Prediction of the senses of helical amphiphilic assemblies from effective intermolecular pair potential: Studies on chiral monolayers and bilayers

It is well-known that the senses (or the handedness) of the helical assemblies formed from compressed monolayers and bilayers of chiral amphiphiles are highly specific about the chirality of the monomers concerned. We present here a molecular approach that can successfully predict the senses of such helical morphologies. The present approach is based on a reduced tractable description in terms of an effective pair potential (EPP) which depends on the distance of separation and the relative orientations of the two amphiphiles. This approach explicitly considers the pairwise intermolecular interactions between the groups attached to the chiral centers of the two neighboring amphiphiles. It is found that for a pair of the same kind of enantiomers the minimum energy configuration favors a twist angle between molecules and that this twist from neighbor to neighbor gives rise to the helicity of the aggregate. From the known twist angles at the minimum energy configuration the successive arrangement of an array of molecules can be predicted. Therefore, the sense of the helicity can be predicted from the molecular interactions. The predicted senses of the helical structures are in complete agreement with all known experimental results.

Phase transitions in the anchored organic bilayers of long-chain alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4; n=12, 16, 18

The single perovskite slab alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4, n = 12, 16, 18, display two phase transitions, just above room temperature, associated with changes in the alkylammonium chains. We have followed these two phase transitions using scanning calorimetry, X-ray powder diffraction, and IR and Raman spectroscopies. We find the first phase transition to be associated with symmetry changes arising from a dynamic rotational disordering of the ammonium headgroup of the chain whereas the second transition, the melting of the chains in two dimensions, is characterized by an increased conformational disorder of the methylene units of the alkyl chains. We examine these phase transitions in light of the interesting optical properties of these materials, as well as the relevance of these systems as models for phase transitions in lipid bilayers.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.