Theory of lipid polymorphism: application to phosphatidylethanolamine and phosphatidylserine

Pathological transitions in myelin membranes driven by environmental and multiple sclerosis conditions

In demyelination diseases, such as multiple sclerosis, the structure of the axons’ protective sheaths is disrupted. Due to the proximity of cytoplasmic myelin membrane to structural phase transition, minor alterations in the local environmental conditions can have devastating results. Using small-angle X-ray scattering and cryogenic transmission electron microscopy, we show that drastic structural reorganization and instabilities of myelin membrane are linked to specific environmental conditions and molecular composition in healthy and diseased states. These instabilities involve phase transition from the healthy lamellar membranes to pathological inverted hexagonal phase. These results highlight that local environmental conditions are critical for myelin function and should be considered as alternative routes for early pathology and as a means to avoid the initiation of demyelination.

Multiple sclerosis (MS) is an autoimmune disease, leading to the destruction of the myelin sheaths, the protective layers surrounding the axons. The etiology of the disease is unknown, although there are several postulated environmental factors that may contribute to it. Recently, myelin damage was correlated to structural phase transition from a healthy stack of lamellas to a diseased inverted hexagonal phase as a result of the altered lipid stoichiometry and low myelin basic protein (MBP) content. In this work, we show that environmental conditions, such as buffer salinity and temperature, induce the same pathological phase transition as in the case of the lipid composition in the absence of MBP. These phase transitions have different transition points, which depend on the lipid’s compositions, and are ion specific. In extreme environmental conditions, we find an additional dense lamellar phase and that the native lipid composition results in similar pathology as the diseased composition. These findings demonstrate that several local environmental changes can trigger pathological structural changes. We postulate that these structural modifications result in myelin membrane vulnerability to the immune system attacks and thus can help explain MS etiology.

How cholesterol could be drawn to the cytoplasmic leaf of the plasma membrane by phosphatidylethanolamine

In the mammalian plasma membrane, cholesterol can translocate rapidly between the exoplasmic and cytoplasmic leaves, so that its distribution between them should be given by the equality of its chemical potential in the leaves. Due to its favorable interaction with sphingomyelin, which is almost entirely in the outer leaf, one expects the great majority of cholesterol to be there also. Experimental results do not support this, implying that there is some mechanism attracting cholesterol to the inner leaf. We hypothesize that it is drawn there to reduce the bending free energy of the membrane caused by the presence of PE (phosphatidylethanolamine). It does this in two ways: first by simply diluting the amount of PE in the inner leaf, and second by ordering the tails of the PE to reduce its spontaneous curvature. Incorporating this mechanism into a model free energy for the bilayer, we find that between 50 and 60% of the total cholesterol should be in the inner leaf of human erythrocytes.

Unit cell structure of water-filled monoolein into inverted hexagonal (H(II)) mesophase modeled by molecular dynamics

The study investigates the unit cell structure of inverted hexagonal (H(II)) mesophase composed of monoolein (1-monoolein, GMO) and water using atomistic molecular dynamics methods without imposing any restraints on lipid and water molecules. Statistically meaningful and very contrast images of the radial mass density distribution, scrutinizing also the separate components water, monoolein, the polar headgroups of the lipids, the double bond, and the termini of the hydrocarbon chain (the tail), are obtained. The lipid/water interface structure is analyzed based on the obtained water density distribution, on the estimated number of hydrogen bonds per monoolein headgroup, and on the headgroup-water radial distribution functions. The headgroup mass density distribution demonstrates hexagonal shape of the monoolein/water interface that is well-defined at higher water/monoolein ratios. Water interacts with the headgroups by forming a three-layer diffusive mass density distribution, and each layer's shape is close to hexagonal, which is an indication of long-range structural interactions. It is found that the monoolein headgroups form a constant number of hydrogen bonds leaving an excessive amount of water molecules outside the first lipid coordination sphere. Furthermore, the quantity of water at the monoolein/water interface increases steadily upon extension of the unit cell, so the interface should have a very dynamic structure. Investigation of the hydrocarbon residues reveals high compression and well-expressed structuring of the tails. The tails form a very compressed and constrained structure of defined layers across the unit cell with properties corresponding to a more densely packed nonpolar liquid (oil). Due to the hexagonal shape the 2D packing frustration is constant and does not depend on the water content. All reported structural features are based on averaging of the atomic coordinates over the time-length of the simulation trajectories. That kind of processing allows the observation of the water/GMO interface shape and its stability and mobility at a time scale close to the ones of the intermolecular interactions.

Effect of electrostatic interactions on phase stability of cubic phases of biomembranes

We investigated effect of electrostatic interactions due to surfacecharges on structures and stability of cubic phases of monoolein (MO)membrane using the small-angle X-ray scattering method. Firstly, wechanged the surface charge density of the membrane by usingdioleoylphosphatidic acid (DOPA). As increasing DOPA concentration in themembrane at 30 wt % lipid concentration, a Q(224) to Q(229) phasetransition occurred at 0.6 mol % DOPA, and at and above 25 mol %, DOPA/MOmembranes were in the L(α) phase. NaCl in the bulk phase reduced theeffect of DOPA. These results indicate that as the electrostaticinteractions increase, the most stable phase changes as follows: Q(224)⇒ Q(229) ⇒ L(α). The increase in DOPAconcentration reduced the absolute value of spontaneous curvature of themembrane, | H(0) |. Secondly, we changed the surface charge of themembrane by adding a de novo designed peptide, which has netpositive charges and a binding site on the electrically neutral membraneinterface. The peptide-1 (WLFLLKKK) induced a Q(224) to Q(229)phase transition in the MO membrane at low peptide concentration. As NaClconcentration increases, the MO/peptide-1 membrane changed from Q(229)to Q(224) phase. The increase in peptide-1 concentration reduced |H(0) |. Based on these results, the stability of the cubic phases and themechanism of phase transition between cubic phase and L(α) phase arediscussed.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: theory versus experiment

Aqueous dispersions of phosphatidylglycerol (PG) lipids may present an anomalous chain-melting transition at low ionic strengths, as seen by different experimental techniques such as calorimetry or light scattering. The anomaly disappears at high ionic strengths or for longer acyl-chain lengths. In this article, we use a statistical model for the bilayer that distinguishes both lipid chain and headgroup states in order to compare model and experimental thermotropic and electrical properties. The effective van der Waals interactions among hydrophobic chains compete with the electrostatic repulsions between polar headgroups, which may be ionized (counterion dissociated) or electrically neutral (associated with counterions). Electric degrees of freedom introduce new thermotropic charge-ordered phases in which headgroup charges may be spatially ordered, depending on the electrolyte ionic strength, introducing a new rationale for experimental data on PGs. The thermal phases presented by the model for different chain lengths, at fixed ionic strength, compare well with an experimental phase diagram constructed on the basis of differential scanning calorimetry profiles. In the case of dispersions of DMPG (dimyristoyl phosphatidylglycerol) with added monovalent salt, the model properties reproduce the main features displayed by data from differential scanning calorimetry as well as the characteristic profile for the degree of ionization of the bilayer surface across the anomalous transition region, obtained from the theoretical interpretation of electrokinetic (conductivity and electrophoretic mobility) measurements.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: a statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types--each one associated, respectively, with the polar-headgroup and the acyl-chain states--which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Fourth-order curvature energy model for the stability of bicontinuous inverted cubic phases in amphiphile-water systems

The bicontinuous inverted cubic (Q(II)) phases of amphiphiles in water have many practical applications. It is necessary to understand the stability of these phases as a function of composition and ambient conditions in order to make the best use of them. Moreover, many biomembrane lipids and some biomembrane lipid extracts form Q(II) phases. The stability of Q(II) phases in a given lipid composition is closely related to the susceptibility of that composition to membrane fusion: changes in composition that stabilize Q(II) phases usually increase the rate of membrane fusion. However, the factors determining Q(II) phase stability are not fully understood. Previously, an expression was derived for the curvature free energy of Q(II) phases with respect to that of the lamellar (L(alpha)) phase using a model for the curvature energy with terms up to fourth order in curvature as formulated by Mitov. Here this model is extended to account for the effects of water content on Q(II) phase stability. It is shown that the observed L(alpha)/Q(II) phase-transition temperature, transition enthalpy, and transition kinetics are all sensitive to water content. The same observables also become sensitive to small noncurvature energy contributions to the total free-energy difference between the Q(II) and L(alpha) phases, especially the unbinding energy in the L(alpha) phase. These predictions rationalize earlier observations of Q(II) phase formation in N-monomethylated dioleoylphosphatidylethanolamine that otherwise appear to be inconsistent. The model also provides a fundamental explanation of the hysteresis typically observed in transitions between the L(alpha) and Q(II) phases. It is an accurate model of Q(II) phase stability when the ratio of the volume fraction of the lipid in the Q(II) phase unit cell is < or = 0.5.

A lipid matrix model of membrane raft structure

Domains in cell membranes are created by lipid-lipid interactions and are referred to as membrane rafts. Reliable isolation methods have been developed which have shown that rafts from the same membranes have different proteins and can be sub-fractionated by immunoaffinity methods. Analysis of these raft subfractions shows that they are also comprised of different molecular species of lipids. The major lipid classes present are phospholipids, glycosphingolipids and cholesterol. Model studies show that mixtures of phospholipids, particularly sphingomyelin, and cholesterol form liquid-ordered phase with properties intermediate between a gel and fluid phase. This type of liquid-ordered phase dominates theories of domain formation and raft structure in biological membranes. Recently it has been shown that sphingolipids with long (22-26C) N-acyl fatty acids form quasi-crystalline bilayer structures with diacylphospholipids that have well-defined stoichiometries. A two tier heuristic model of membrane raft structure is proposed in which liquid-ordered phase created by a molecular complex between sphingolipids with hydrocarbon chains of approximately equal length and cholesterol acts as a primary staging area for selecting raft proteins. Tailoring of the lipid anchors of raft proteins takes place at this site. Assembly of lipid-anchored proteins on a scaffold of sphingolipids with asymmetric hydrocarbon chains and phospholipids arranged in a quasi-crystalline bilayer structure serves to concentrate and orient the proteins in a manner that couples them functionally within the membrane. Specificity is inherent in the quasi-crystalline lipid structure of liquid-ordered matrices formed by both types of complex into which protein lipid anchors are interpolated. An interaction between the sugar residues of the glycolipids and the raft proteins provides an additional level of specificity that distinguishes one raft from another.

Ionic control of the metastable inner leaflet of the plasma membrane: Fusions natural and artefactual

The phospholipids of the inner and outer leaflets of the plasma membrane face chemically very different environments, and are specialized to serve different needs. While lipids of the outer leaflet are inherently stable in a lamellar (bilayer) phase, the main lipid of the inner layer, phosphatidylethanolamine (PE), does not form a lamellar phase unless evenly mixed with phosphatidylserine (PS(-)). This mixture can be readily perturbed by factors that include an influx of Ca(2+) that chelates the negatively charged PS(-), thereby destabilizing PE. The implications of this metastability of the inner leaflet for vesicular trafficking, and experimentally for the isolation of detergent-resistant membrane domains (DRMs) at physiological temperature, are considered.

Structure and organization of nanosized-inclusion-containing bilayer membranes

Based on a considerable amount of experimental evidence for lateral organization of lipid membranes which share astonishingly similar features in the presence of different inclusions, we use a hybrid self-consistent field theory (SCFT)/density-functional theory (DFT) approach to deal with bilayer membranes embedded by nanosized inclusions and explain experimental findings. Here, the hydrophobic inclusions are simple models of hydrophobic drugs or other nanoparticles for biomedical applications. It is found that lipid/inclusion-rich domains are formed at moderate inclusion concentrations and disappear with the increase in the concentration of inclusions. At high inclusion content, chaining of inclusions occurs due to the effective depletion attraction between inclusions mediated by lipids. Meanwhile, the increase in the concentration of inclusions can also cause thickening of the membrane and the distribution of inclusions undergoes a layering transition from one-layer structure located in the bilayer midplane to two-layer structure arranged into the two leaflets of a bilayer. Our theoretical predictions address the complex interactions between membranes and inclusions suggesting a unifying mechanism which reflects the competition between the conformational entropy of lipids favoring the formation of lipid- and inclusion-rich domains in lipids and the steric repulsion of inclusions leading to the uniform dispersion.

HIV-1 matrix organizes as a hexamer of trimers on membranes containing phosphatidylinositol-(4,5)-bisphosphate

The human immunodeficiency virus type 1 (HIV-1) matrix (MA) protein represents the N-terminal domain of the HIV-1 precursor Gag (PrGag) protein and carries an N-terminal myristate (Myr) group. HIV-1 MA fosters PrGag membrane binding, as well as assembly of envelope (Env) proteins into virus particles, and recent studies have shown that HIV-1 MA preferentially directs virus assembly at plasma membrane sites enriched in cholesterol and phosphatidylinositol-(4,5)-bisphosphate (PI[4,5]P(2)). To characterize the membrane binding of MA and PrGag proteins, we have examined how Myr-MA proteins, and proteins composed of Myr-MA and its neighbor Gag capsid (CA) protein associate on membranes containing cholesterol and PI[4,5]P(2). Our results indicate that Myr-MA assembles as a hexamer of trimers on such membranes, and imply that MA trimers interconnect CA hexamer rings in immature virus particles. Our observations suggest a model for the organization of PrGag proteins, and for MA-Env protein interactions.

Isolation at physiological temperature of detergent-resistant membranes with properties expected of lipid rafts: the influence of buffer composition

The failure of most non-ionic detergents to release patches of DRM (detergent-resistant membrane) at 37 °C undermines the claim that DRMs consist of lipid nanodomains that exist in an Lo (liquid ordered) phase on the living cell surface. In the present study, we have shown that inclusion of cations (Mg2+, K+) to mimic the intracellular environment stabilizes membranes during solubilization sufficiently to allow the isolation of DRMs at 37 °C, using either Triton X-100 or Brij 96. These DRMs are sensitive to chelation of cholesterol, maintain outside-out orientation of membrane glycoproteins, have prolonged (18 h) stability at 37 °C, and are vesicles or sheets up to 150–200 nm diameter. DRMs containing GPI (glycosylphosphatidylinositol)-anchored proteins PrP (prion protein) and Thy-1 can be separated by immunoaffinity isolation, in keeping with their separate organization and trafficking on the neuronal surface. Thy-1, but not PrP, DRMs are associated with actin. EM (electron microscopy) immunohistochemistry shows most PrP, and some Thy-1, to be clustered on DRMs, again maintaining their organization on the neuronal surface. For DRMs labelled for either protein, the bulk of the surface of the DRM is not labelled, indicating that the GPI-anchored protein is a minor component of its lipid domain. These 37 °C DRMs thus have properties expected of raft membrane, yet pose more questions about how proteins are organized within these nanodomains.

Food structure and functionality: a soft matter perspective

The structure and functionality of foods are described from the perspective of recent advances in soft condensed matter physics. An overview is given of the structure and properties of food materials in terms of the physically relevant length scales. Recent developments in the understanding of the physics of gels, micelles, liquid crystals, biopolymer complexes and amorphous carbohydrates are presented.

Anomalous phase sequences in lyotropic liquid crystals

We present a coarse-grained model in order to describe the unusual sequence of mesophases observed in aqueous solutions of nonionic lipids, such as monoolein. The lipid molecules are modeled as a rigid head and a flexible Gaussian tail, and water is treated explicitly. A key component of the model is thermally reversible hydrogen bonding between the lipid head and water resulting in changes in both head volume and the interactions of the hydrated head with its surroundings. Phase diagrams obtained from unit-cell self-consistent field simulations capture the qualitative thermotropic and lyotropic phase behavior of the monoolein-water system. The unusual phase sequences result from a competition between hydrogen bond formation, changes in head volume and interactions, lipid tail entropy, and the hydrophobic effect.

Improved dissipative particle dynamics simulations of lipid bilayers

The authors introduce a new parameterization for the dissipative particle dynamics simulations of lipid bilayers. In this parameterization, the conservative pairwise forces between beads of the same type in two different hydrophobic chains are chosen to be less repulsive than the water-water interaction, but the intrachain bead interactions are the same as the water-water interaction. For a certain range of parameters, the new bilayer can only be stretched up to 30% before it ruptures. Membrane tension, density profiles, and the in-plane lipid diffusion coefficient of the new bilayer are discussed in detail. They find two kinds of finite size effects that influence the membrane tension: lateral finite size effects, for which larger membranes rupture at a smaller stretch, and transverse finite size effects, for which tensionless bilayers are more compact in larger systems. These finite size effects become rather small when the simulation box is sufficiently large.

Alcohols reduce lateral membrane pressures: predictions from molecular theory

We explore the effects of alcohols on fluid lipid bilayers using a molecular theory with a coarse-grained model. We show that the trends predicted from the theory in the changes in area per lipid, alcohol concentration in the bilayer, and area compressibility modulus, as a function of alcohol chain length and of the alcohol concentration in the solvent far from the bilayer, follow those found experimentally. We then use the theory to study the effect of added alcohol on the lateral pressure profile across the membrane, and find that added alcohol reduces the surface tensions at both the headgroup/solvent and headgroup/tailgroup interfaces, as well as the lateral pressures in the headgroup and tailgroup regions. These changes in lateral pressures could affect the conformations of membrane proteins, providing a nonspecific mechanism for the biological effects of alcohols on cells.

Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state

Protein-protein binding usually involves structural changes that may extend beyond the rearrangements on a local scale, and cannot be explained by a classical lock-and-key mechanism. Several models have been advanced to explain the flexible binding of proteins such as the induced fit mechanism where the ligand is postulated to induce a conformational change at the interaction site upon binding, or the preexisting equilibrium hypothesis that assumes that protein samples an ensemble of conformations at equilibrium conditions and that the ligand binds selectively to an active conformation. We explored the equilibrium motions of proteins that exhibit relatively large (nonlocal) conformational changes upon protein binding using the Gaussian network model and the anisotropic network model of protein dynamics. For four complexes, LIR-1/HLA-A2, Actin/DNase I, CDK2/cyclin, and CDK6/p16(INK4a), the motions calculated for the monomer exhibiting the largest conformational change, in its unbound (free) form, correlate with the experimentally observed structural changes upon binding. This study emphasizes the preexisting equilibrium/conformational selection as a mechanism for protein-protein interaction and lends support the concept that proteins, in their native conformation, are predisposed to undergo conformational fluctuations that are relevant to, or even required for, their biological functions.

Phase behavior of a blend of polymer-tethered nanoparticles with diblock copolymers

Using the self-consistent field theory (SCFT), we investigate the phase behavior of a mixture of diblock copolymers and nanoparticles with monodisperse polymer chains tethered to their surfaces. We assume the size of the nanoparticles to be much smaller than that of the attached polymer chains and therefore model the particles with their grafted polymer "shell" as star polymers. The polymer chains attached to the particles are of the same species as one of the blocks of the symmetric diblock copolymer. Of primary interest is how to tune the shell of the particle by changing both the length and number of tethered polymers in order to achieve higher loading of nanoparticles within an ordered structure without macrophase separation occurring. We find that the phase behavior of the system is very sensitive to the size of the particle including its tethered shell. The region of microphase separation is increased upon decreasing the star polymer size, which may be achieved by shortening and/or removing tethered polymer chains. To explore the possible structures in these systems we employ SCFT simulations that provide insight into the arrangement of the different species in these complex composites.

Density functional theory approach for coarse-grained lipid bilayers

Lipid bilayers are important inhomogeneous fluid systems that mediate the environment of cells and the interaction of cells with their environment. A variety of approaches have been taken to model the lipid molecules in bilayers, from all atom molecular dynamics to rigid body liquid crystals. In this paper we discuss the application of a density functional theory approach that treats the lipid molecules at the coarse-grained level of a freely jointed chain. This approach allows for compressibility effects, and can therefore be used to study not only the long range structure in lipid bilayers, but also the nanoscale structure induced in the bilayer when the lipids crystallize or when an inclusion (e.g., an embedded protein) is present. This paper presents a detailed analysis of fluid bilayers and lamellae predicted by the theory. In particular we locate solutions with zero surface tension. We calculate the phase diagram for all possible phases with planar symmetry, including uniform macrophases. Surprisingly, we find a first-order phase transition from the lamellar phase to an isolated bilayer phase on lowering the temperature. This transition appears to be driven by solvent packing effects. A further lowering of the temperature leads to a set of highly ordered bilayers.

Low pH Stabilizes the Inverted Hexagonal II Phase in Dipalmitoleoylphosphatidylethanolamine Membrane

Dipalmitoleoylphosphatidylethanolamine (DPOPE) membrane is in the L(α) phase in neutral pH at 20 °C. The results of small-angle X-ray scattering (SAXS) indicate that an L(α) to H(II) phase transition in DPOPE membranes occurred at pH 1.9 in the absence of salt, and at pH 2.8 in the presence of 0.5 M KCl, at fully hydrated condition at 20 °C. The spontaneous curvature of DPOPE monolayer membrane did not change with a decrease in pH values. To elucidate the mechanism, we have investigated the effect of the cationic dioctadecyldimethylammonium (DODMA) on the structure and phase behavior of DPOPE membrane. The result shows that DODMA stabilizes the H(II) phase rather than the L(α) phase in DPOPE membrane at its low concentrations. Based on these results, the H(II) phase stability of DPOPE membrane due to low pH is discussed in terms of the spontaneous curvature of the monolayer membrane and the packing energy of acyl chains in the membrane.

Effect of electrostatic interactions on phase stability of cubic phases of membranes of monoolein/dioleoylphosphatidic acid mixtures

To elucidate effects of electrostatic interactions resulting from surface charges on structures and phase stability of cubic phases of lipid membranes, membranes of 1-monoolein (MO) and dioleoylphosphatidic acid (DOPA) (DOPA/MO membrane) mixtures have been investigated by small-angle x-ray scattering method. As increasing DOPA concentration in the DOPA/MO membrane at 30 wt% lipid concentration, a phase transition from Q(224) to Q(229) phase occurred at 0.6 mol% DOPA, and at and above 25 mol% DOPA, DOPA/MO membranes were in the L(alpha) phase. As NaCl concentration in the bulk phase increased, for 10% DOPA/90% MO membrane in excess water, a Q(229) to Q(224) phase transition occurred at 60 mM NaCl, and then a Q(224) to H(II) phase transition occurred at 1.2 M NaCl. Similarly, for 30% DOPA/70% MO membrane in excess water, at low NaCl concentrations it was in the L(alpha) phase, but at and above 0.50 M NaCl it was in the Q(224) phase, and then at 0.65 M NaCl a Q(224) to H(II) phase transition occurred. These results indicate that the electrostatic interactions in the membrane interface make the Q(229) phase more stable than the Q(224) phase, and that, at larger electrostatic interactions, the L(alpha) phase is more stable than the cubic phases (Q(224) and Q(229)). We have found that the addition of tetradecane to the MO membrane induced a Q(224)-to-H(II) phase transition and also that to the 30% DOPA/70% MO membrane induced an L(alpha)-to-H(II) phase transition. By using these membranes, the effect of the electrostatic interactions resulting from the membrane surface charge (DOPA) on the spontaneous curvature of the monolayer membrane has been investigated. The increase in DOPA concentration in the DOPA/MO membrane reduced the absolute value of spontaneous curvature of the membrane. In the 30% DOPA/70% MO membrane, the absolute value of spontaneous curvature of the membrane increased with an increase in NaCl concentration. On the basis of these new results, the phase stability of DOPA/MO membranes can be reasonably explained by the spontaneous curvature of the monolayer membrane and a curvature elastic energy of the membrane.

Theory of tunable pH-sensitive vesicles of anionic and cationic lipids or anionic and neutral lipids

The design of vesicles that become unstable at an easily tuned value of pH is of great interest for targeted drug delivery. We present a microscopic theory for two forms of such vesicles. A model of lipids introduced by us previously is applied to a system of ionizable anionic lipid and permanently charged cationic lipid. We calculate the pH at which the lamellar phase becomes unstable with respect to an inverted hexagonal one, a value that depends continuously on the system composition. Identifying this instability with that displayed by unilamellar vesicles undergoing fusion, we obtain very good agreement with the recent experimental data of Hafez, Ansell, and Cullis, (2000, Biophys. J. 79:1438-1446) on the pH at which fusion occurs versus vesicle composition. We explicate the mechanism in terms of the role of the counterions. This understanding suggests that a system of a neutral, nonlamellar-forming lipid stabilized by an anionic lipid would serve equally well for preparing tunable, pH-sensitive vesicles. Our calculations confirm this. Further, we show that both forms of vesicle have the desirable feature of exhibiting a regime in which the pH at instability is a rapidly varying function of the vesicle composition.

Aggregation models at high packing fraction

Dense phases of micellar aggregates have strong molecular correlation at two different levels: that of the molecules forming a micelle and that between micelles, leading to a possible phase transition from a micellar fluid to a micellar crystal. The global phase diagram may also include lamellar and other dense phases, which do not have a micellar structure. We present here a generic approach to deal with these systems through a two-level density-functional description, to first describe an isolated micellar aggregate and then the dense micellar system, obtaining the free energy in a self-consistent way from the molecular interactions. Nonmicellar dense phases are included with the same density-functional approach applied at the first level. The results are shown to be very accurate for a one-dimensional model with exact solution, and the method is then applied to a three-dimensional amphiphile model that had been successfully used to describe the properties of diluted amphiphile solutions.