PubMed Health⌕ Search

Biomedical subjects

A Ben-Shaul

Publications and source records attributed to A Ben-Shaul.

At least 19 recordsLinked to original sources
PubMed2004

Tilt modulus of a lipid monolayer.

In addition to the familiar bending and stretching deformations, lipid monolayers and bilayers in their disordered state are often subjected to tilt deformations, occurring for instance in structural rearrangements accompanying membrane fusion, or upon insertion of "oblique" hydrophobic proteins into lipid bilayers. We study the elastic response of a flat lipid monolayer to a tilt deformation, using the spatial and conformational average of the chain end-to-end vector from the membrane normal to define a macroscopic membrane tilt. The physical origin and magnitude of the corresponding tilt modulus kappat is analyzed using two complementary theoretical approaches. The first is a phenomenological model showing that the tilt and bending deformations are decoupled and the effects of inter-chain correlations on the tilt modulus is small. The second is based on a molecular-level mean-field theory of chain packing, enabling numerical evaluation of the tilt modulus for realistic, multi-conformation, chain models. Both approaches reveal that the tilt modulus involves two major contributions. The first is elastic in origin, arising from the stretching of the hydrocarbon chains upon a tilt deformation and reflecting the loss of chain conformational freedom associated with chain stretching. The second, purely entropic, contribution results from the constraints imposed by a tilt deformation on the fluctuations of chain director orientations. Using the chain-packing theory we compute the two contributions numerically as a function of the cross-sectional area per chain. The elastic and entropic terms are shown to dominate the value of kappat for small and large areas per chain, respectively. For typical cross-sectional areas of lipid chains in biological membranes they are of comparable magnitude, yielding kappat approximately 0.2kBT/A2.

Computer Simulation↗
PubMed2004

Modeling of cationic lipid-DNA complexes.

Cationic lipid-DNA complexes, often referred to as lipoplexes, are formed spontaneously in aqueous solutions upon mixing DNA and liposomes composed of cationic and nonionic lipids. Understanding the mechanisms underlying lipoplex formation, structure and phase behavior is crucial for their further development and design as non-viral transfection vectors in gene therapy. From a physical point of view, lipoplexes are ordered, self-assembled, composite aggregates. Their preferred spatial geometry and phase behavior are governed by a delicate coupling between the electrostatic interactions which drive lipoplex formation and the elastic properties of the constituent lipid layers, both depending on the molecular nature and composition of the lipid mixture. In this review we outline some recent efforts to model the microscopic structure, energetic and phase behavior of cationic lipid-DNA mixtures, focusing on the two principal aggregation geometries: the lamellar (L(alpha)(C)), or "sandwich" complexes, and the hexagonal (H(II)(C)), or "honeycomb" complexes. We relate the structural and thermodynamic properties of these two "canonical" lipoplex morphologies to their appearance in phase diagrams of DNA-lipid mixtures, emphasizing the crucial role fulfilled by the molecular packing characteristics of the cationic and neutral lipids, as reflected in the curvature elastic properties of the mixed lipid layer.

Cations↗
PubMed2001

DNA packaging and ejection forces in bacteriophage.

We calculate the forces required to package (or, equivalently, acting to eject) DNA into (from) a bacteriophage capsid, as a function of the loaded (ejected) length, under conditions for which the DNA is either self-repelling or self-attracting. Through computer simulation and analytical theory, we find the loading force to increase more than 10-fold (to tens of piconewtons) during the final third of the loading process; correspondingly, the internal pressure drops 10-fold to a few atmospheres (matching the osmotic pressure in the cell) upon ejection of just a small fraction of the phage genome. We also determine an evolution of the arrangement of packaged DNA from toroidal to spool-like structures.

Bacteriophages↗
PubMed2000

Association entropy in adsorption processes.

The association of two species to form a bound complex, e.g., the binding of a ligand to a protein or the adsorption of a peptide on a lipid membrane, involves an entropy loss, reflecting the conversion of free translational and rotational degrees of freedom into bound motions. Previous theoretical estimates of the standard entropy change in bimolecular binding processes, DeltaS(o), have been derived from the root-mean-square fluctuations in protein crystals, suggesting DeltaS(o) approximately -50 e.u., i.e., TDeltaS degrees approximately -25 kT = -15 kcal/mol. In this work we focus on adsorption, rather than binding processes. We first present a simple statistical-thermodynamic scheme for calculating the adsorption entropy, including its resolution into translational and rotational contributions, using the known distance-orientation dependent binding (adsorption) potential. We then utilize this scheme to calculate the free energy of interaction and entropy of pentalysine adsorption onto a lipid membrane, obtaining TDeltaS(o) approximately -1.7 kT approximately -1.3 kcal/mol. Most of this entropy change is due to the conversion of one free translation into a bound motion, the rest arising from the confinement of two rotational degrees of freedom. The smaller entropy loss in adsorption compared to binding processes arises partly because a smaller number of degrees of freedom become restricted, but mainly due to the fact that the binding potential is much "softer."

Adsorption↗
PubMed2000

Lipid demixing and protein-protein interactions in the adsorption of charged proteins on mixed membranes.

The adsorption free energy of charged proteins on mixed membranes, containing varying amounts of (oppositely) charged lipids, is calculated based on a mean-field free energy expression that accounts explicitly for the ability of the lipids to demix locally, and for lateral interactions between the adsorbed proteins. Minimization of this free energy functional yields the familiar nonlinear Poisson-Boltzmann equation and the boundary condition at the membrane surface that allows for lipid charge rearrangement. These two self-consistent equations are solved simultaneously. The proteins are modeled as uniformly charged spheres and the (bare) membrane as an ideal two-dimensional binary mixture of charged and neutral lipids. Substantial variations in the lipid charge density profiles are found when highly charged proteins adsorb on weakly charged membranes; the lipids, at a certain demixing entropy penalty, adjust their concentration in the vicinity of the adsorbed protein to achieve optimal charge matching. Lateral repulsive interactions between the adsorbed proteins affect the lipid modulation profile and, at high densities, result in substantial lowering of the binding energy. Adsorption isotherms demonstrating the importance of lipid mobility and protein-protein interactions are calculated using an adsorption equation with a coverage-dependent binding constant. Typically, at bulk-surface equilibrium (i.e., when the membrane surface is "saturated" by adsorbed proteins), the membrane charges are "overcompensated" by the protein charges, because only about half of the protein charges (those on the hemispheres facing the membrane) are involved in charge neutralization. Finally, it is argued that the formation of lipid-protein domains may be enhanced by electrostatic adsorption of proteins, but its origin (e.g., elastic deformations associated with lipid demixing) is not purely electrostatic.

Adsorption↗
PubMed2000

The phase behavior of cationic lipid-DNA complexes.

We present a theoretical analysis of the phase behavior of solutions containing DNA, cationic lipids, and nonionic (helper) lipids. Our model allows for five possible structures, treated as incompressible macroscopic phases: two lipid-DNA composite (lipoplex) phases, namely, the lamellar (L(alpha)(C)) and hexagonal (H(II)(C)) complexes; two binary (cationic/neutral) lipid phases, that is, the bilayer (L(alpha)) and inverse-hexagonal (H(II)) structures, and uncomplexed DNA. The free energy of the four lipid-containing phases is expressed as a sum of composition-dependent electrostatic, elastic, and mixing terms. The electrostatic free energies of all phases are calculated based on Poisson-Boltzmann theory. The phase diagram of the system is evaluated by minimizing the total free energy of the three-component mixture with respect to all the compositional degrees of freedom. We show that the phase behavior, in particular the preferred lipid-DNA complex geometry, is governed by a subtle interplay between the electrostatic, elastic, and mixing terms, which depend, in turn, on the lipid composition and lipid/DNA ratio. Detailed calculations are presented for three prototypical systems, exhibiting markedly different phase behaviors. The simplest mixture corresponds to a rigid planar membrane as the lipid source, in which case, only lamellar complexes appear in solution. When the membranes are "soft" (i.e., low bending modulus) the system exhibits the formation of both lamellar and hexagonal complexes, sometimes coexisting with each other, and with pure lipid or DNA phases. The last system corresponds to a lipid mixture involving helper lipids with strong propensity toward the inverse-hexagonal phase. Here, again, the phase diagram is rather complex, revealing a multitude of phase transitions and coexistences. Lamellar and hexagonal complexes appear, sometimes together, in different regions of the phase diagram.

Biophysical Phenomena↗
PubMed1999

Molecular theory of lipid-protein interaction and the Lalpha-HII transition.

We present a molecular-level theory for lipid-protein interaction and apply it to the study of lipid-mediated interactions between proteins and the protein-induced transition from the planar bilayer (Lalpha) to the inverse-hexagonal (HII) phase. The proteins are treated as rigid, membrane-spanning, hydrophobic inclusions of different size and shape, e.g., "cylinder-like," "barrel-like," or "vase-like." We assume strong hydrophobic coupling between the protein and its neighbor lipids. This means that, if necessary, the flexible lipid chains surrounding the protein will stretch, compress, and/or tilt to bridge the hydrophobic thickness mismatch between the protein and the unperturbed bilayer. The system free energy is expressed as an integral over local molecular contributions, the latter accounting for interheadgroup repulsion, hydrocarbon-water surface energy, and chain stretching-tilting effects. We show that the molecular interaction constants are intimately related to familiar elastic (continuum) characteristics of the membrane, such as the bending rigidity and spontaneous curvature, as well as to the less familiar tilt modulus. The equilibrium configuration of the membrane is determined by minimizing the free energy functional, subject to boundary conditions dictated by the size, shape, and spatial distribution of inclusions. A similar procedure is used to calculate the free energy and structure of peptide-free and peptide-rich hexagonal phases. Two degrees of freedom are involved in the variational minimization procedure: the local length and local tilt angle of the lipid chains. The inclusion of chain tilt is particularly important for studying noncylindrical (for instance, barrel-like) inclusions and analyzing the structure of the HII lipid phase; e.g., we find that chain tilt relaxation implies strong faceting of the lipid monolayers in the hexagonal phase. Consistent with experiment, we find that only short peptides (large negative mismatch) can induce the Lalpha --> HII transition. At the transition, a peptide-poor Lalpha phase coexists with a peptide-rich HII phase.

Elasticity↗
PubMed1998

Structure, stability, and thermodynamics of lamellar DNA-lipid complexes.

We develop a statistical thermodynamic model for the phase evolution of DNA-cationic lipid complexes in aqueous solution, as a function of the ratios of charged to neutral lipid and charged lipid to DNA. The complexes consist of parallel strands of DNA intercalated in the water layers of lamellar stacks of mixed lipid bilayers, as determined by recent synchrotron x-ray measurements. Elastic deformations of the DNA and the lipid bilayers are neglected, but DNA-induced spatial inhomogeneities in the bilayer charge densities are included. The relevant nonlinear Poisson-Boltzmann equation is solved numerically, including self-consistent treatment of the boundary conditions at the polarized membrane surfaces. For a wide range of lipid compositions, the phase evolution is characterized by three regions of lipid to DNA charge ratio, rho: 1) for low rho, the complexes coexist with excess DNA, and the DNA-DNA spacing in the complex, d, is constant; 2) for intermediate rho, including the isoelectric point rho = 1, all of the lipid and DNA in solution is incorporated into the complex, whose inter-DNA distance d increases linearly with rho; and 3) for high rho, the complexes coexist with excess liposomes (whose lipid composition is different from that in the complex), and their spacing d is nearly, but not completely, independent of rho. These results can be understood in terms of a simple charging model that reflects the competition between counterion entropy and inter-DNA (rho < 1) and interbilayer (rho > 1) repulsions. Finally, our approach and conclusions are compared with theoretical work by others, and with relevant experiments.

Biophysical Phenomena↗
PubMed1997

DNA-lipid complexes: stability of honeycomb-like and spaghetti-like structures.

A molecular level theory is presented for the thermodynamic stability of two (similar) types of structural complexes formed by (either single strand or supercoiled) DNA and cationic liposomes, both involving a monolayer-coated DNA as the central structural unit. In the "spaghetti" complex the central unit is surrounded by another, oppositely curved, monolayer, thus forming a bilayer mantle. The "honeycomb" complex is a bundle of hexagonally packed DNA-monolayer units. The formation free energy of these complexes, starting from a planar cationic/neutral lipid bilayer and bare DNA, is expressed as a sum of electrostatic, bending, mixing, and (for the honeycomb) chain frustration contributions. The electrostatic free energy is calculated using the Poisson-Boltzmann equation. The bending energy of the mixed lipid layers is treated in the quadratic curvature approximation with composition-dependent bending rigidity and spontaneous curvature. Ideal lipid mixing is assumed within each lipid monolayer. We found that the most stable monolayer-coated DNA units are formed when the charged/neutral lipid composition corresponds (nearly) to charge neutralization; the optimal monolayer radius corresponds to close DNA-monolayer contact. These conclusions are also valid for the honeycomb complex, as the chain frustration energy is found to be negligible. Typically, the stabilization energies for these structures are on the order of 1 k(B)T/A of DNA length, reflecting mainly the balance between the electrostatic and bending energies. The spaghetti complexes are less stable due to the additional bending energy of the external monolayer. A thermodynamic analysis is presented for calculating the equilibrium lipid compositions when the complexes coexist with excess bilayer.

Computer Simulation↗
PubMed1996

Statistical thermodynamic analysis of peptide and protein insertion into lipid membranes.

A statistical thermodynamic approach is used to analyze the various contributions to the free energy change associated with the insertion of proteins and protein fragments into lipid bilayers. The partition coefficient that determines the equilibrium distribution of proteins between the membrane and the solution is expressed as the ratio between the partition functions of the protein in the two phases. It is shown that when all of the relevant degrees of freedom (i.e., those that change their character upon insertion into the membrane) can be treated classically, the partition coefficient is fully determined by the ratio of the configurational integrals and thus does not involve any mass-dependent factors, a conclusion that is also valid for related processes such as protein adsorption on a membrane surface or substrate binding to proteins. The partition coefficient, and hence the transfer free energy, depend only on the potential energy of the protein in the membrane. Expressing this potential as a sum of a "static" term, corresponding to the equilibrium (minimal free energy) configuration of the protein in the membrane, and a "dynamical" term representing fluctuations around the equilibrium configuration, we show that the static term contains the "solvation" and "lipid perturbation" contributions to the transfer free energy. The dynamical term is responsible for the "immobilization" free energy, reflecting the loss of translational and rotational entropy of the protein upon incorporation into the membrane. Based on a recent molecular theory of lipid-protein interactions, the lipid perturbation and immobilization contributions are then expressed in terms of the elastic deformation free energy resulting from the perturbation of the lipid environment by the foreign (protein) inclusion. The model is formulated for cylindrically shaped proteins, and numerical estimates are given for the insertion of an alpha-helical peptide into a lipid bilayer. The immobilization free energy is shown to be considerably smaller than in previous estimates of this quantity, and the origin of the difference is discussed in detail.

Biophysical Phenomena↗
PubMed1996

Free-energy determinants of alpha-helix insertion into lipid bilayers.

A detailed treatment is provided of the various free-energy terms that contribute to the transfer of a polyalanine alpha-helix from the aqueous phase into lipid bilayers. In agreement with previous work, the hydrophobic effect is found to provide the major driving force for helix insertion. However, an opposing effect of comparable magnitude is also identified and is attributed to the large free-energy penalty associated with the desolvation of peptide hydrogen bonds on transfer to the low dielectric environment of the bilayer. Lipid perturbation effects as well as the entropy loss associated with helix immobilization in the bilayer are also evaluated. Two configurations of a membrane-bound 25mer polyalanine helix were found to be lower in free energy than the isolated helix in the aqueous phase. The first corresponds to the case of vertical insertion, in which a helix terminus protrudes from each side of the bilayer. The second minimum is for the case of horizontal insertion, for which the helix is adsorbed upon the surface of the bilayer. The calculated free-energy minima are found to be in good agreement with recent measurements of related systems. Large free-energy barriers resulting from desolvation of unsatisfied hydrogen-bonding groups at the helix termini are obtained for both insertion processes. The barriers for insertion are significantly reduced if the helix termini are assumed to be "capped" through the formation of hydrogen bonds with polar sidechains. For uncapped helices, our results support recently proposed models in which helices are inserted by first adsorbing on the membrane surface and then having one terminus "swing around" so as to penetrate the bilayer.

Biophysical Phenomena↗
PubMed1994

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1993

A molecular model for lipid-protein interaction in membranes: the role of hydrophobic mismatch.

The interaction free energy between a hydrophobic, transmembrane, protein and the surrounding lipid environment is calculated based on a microscopic model for lipid organization. The protein is treated as a rigid hydrophobic solute of thickness dP, embedded in a lipid bilayer of unperturbed thickness doL. The lipid chains in the immediate vicinity of the protein are assumed to adjust their length to that of the protein (e.g., they are stretched when dP > doL) in order to bridge over the lipid-protein hydrophobic mismatch (dP-doL). The bilayer's hydrophobic thickness is assumed to decay exponentially to its asymptotic, unperturbed, value. The lipid deformation free energy is represented as a sum of chain (hydrophobic core) and interfacial (head-group region) contributions. The chain contribution is calculated using a detailed molecular theory of chain packing statistics, which allows the calculation of conformational properties and thermodynamic functions (in a mean-field approximation) of the lipid tails. The tails are treated as single chain amphiphiles, modeled using the rotational isometric state scheme. The interfacial free energy is represented by a phenomenological expression, accounting for the opposing effects of head-group repulsions and hydrocarbon-water surface tension. The lipid deformation free energy delta F is calculated as a function of dP-doL. Most calculations are for C14 amphiphiles which, in the absence of a protein, pack at an average area per head-group ao approximately equal to 32 A2 (doL approximately 24.5 A), corresponding to the fluid state of the membrane. When dP = doL, delta F > 0 and is due entirely to the loss of conformational entropy experienced by the chains around the protein. When dP > doL, the interaction free energy is further increased due to the enhanced stretching of the tails. When dP < doL, chain flexibility (entropy) increases, but this contribution to delta F is overcounted by the increase in the interfacial free energy. Thus, delta F obtains a minimum at dP-doL approximately 0. These qualitative interpretations are supported by detailed numerical calculations of the various contributions to the interaction free energy, and of chain conformational properties. The range of the perturbation of lipid order extends typically over few molecular diameters. A rather detailed comparison of our approach to other models is provided in the discussion.

Biophysical Phenomena↗
PubMed1990

States of aggregation and phase transformations in mixtures of phosphatidylcholine and octyl glucoside.

The result of mixing varying concentrations of the nonionic detergent octyl glucoside (OG) with small unilamellar vesicles (SUV) of egg phosphatidylcholine (PC) made by sonication depends on the ratio between OG and PC in the mixed aggregates. When this molar ratio (Re) is lower than 1.4, the detergent partitions between the PC vesicles and the aqueous medium with a partition coefficient of K = 0.033 mM-1. As a consequence of introduction of OG into the bilayers, the vesicles grow in size. The resultant vesicles have a mean diameter that is an increasing function of Re and is independent of the total PC concentration. Experiments in which the vesicles were loaded with high molecular weight dextran prior to being exposed to OG suggest that the mechanism responsible for the size growth involves lipid transfer rather than fusion. Mixtures with Re values within the range of 1.4-3.2 separate into two macroscopic phases: The lower phase is clear but very viscous. It contains constant OG and PC concentrations and is characterized by an Re value of 3.2, independent of the composition of the whole dispersion. The upper phase contains vesicles of varying concentrations of OG and PC, but a constant Re of 1.4. When the saturating level of 1.4 OG molecules per PC molecule is approached, the concentration of OG monomers in the aqueous medium reaches the value of 16.6 +/- 0.3 mM, which is the apparent cmc of OG in the lipid-containing medium. OG-PC mixed micelles contain at least 3.2 OG molecules per PC molecule. The mixed micelles present at Re = 3.2 apparently have the shape of oblate ellipsoids with a minor axis of about 2 nm and two major axes of about 25 nm. The surface area of the mixed micelles at this point is just sufficient for them to undergo conversion into the smallest possible spherical vesicles of a radius of 12 nm. At Re values above 3.2, the major axis of the mixed micelles becomes smaller as Re increases, while at values of Re below 3.2 the micelles would have been expected to grow very rapidly with decreasing Re. This may explain the partial vesicle closure occurring below Re = 3.2.

Detergents↗