Molecular dynamics simulations of SOPS and sphingomyelin bilayers containing cholesterol

We performed a molecular dynamics simulation of an asymmetric bilayer that contained different lipid mixtures in its outer and inner leaflets. The outer leaflet contained a mixture of sphingomyelin (SM) with cholesterol and the inner leaflet a mixture of stearoyl-oleoyl-phosphatidylserine (SOPS) with cholesterol. For comparison purposes, we also performed two simulations on symmetric bilayers: the first simulation was performed on a bilayer containing a binary mixture of SOPS with cholesterol; the second contained a mixture of SM with cholesterol. We studied the hydrogen-bonding network of the bilayers in our simulations and the difference in the network properties in the monolayers either with SM or SOPS. We observed that in the asymmetric bilayer the properties of monolayers were the same as in the corresponding monolayers in the symmetric bilayers.     

Interactions of the Auxilin-1 PTEN-like Domain with Model Membranes Result in Nanoclustering of Phosphatidyl Inositol Phosphates

Auxilin-1 is a neuron-specific membrane-binding protein involved in a late stage of clathrin-mediated endocytosis. It recruits Hsc70, thus initiating uncoating of the clathrin-coated vesicles. Interactions of auxilin-1 with the vesicle membrane are crucial for this function and are mediated via an N-terminal PTEN-like domain. We have used multiscale molecular dynamics simulations to probe the interactions of the auxilin-1 PTEN-like domain with lipid bilayers containing differing phospholipid composition, including bilayers containing phosphatidyl inositol phosphates. Our results suggest a novel, to our knowledge, model for the auxilin/membrane encounter and subsequent interactions. Negatively charged lipids (especially PIP₂) enhance binding of auxilin to lipid bilayers and facilitate its correct orientation relative to the membrane. Mutations in three basic residues (R301E/R307E/K311E) of the C2 subdomain of the PTEN-like domain perturbed its interaction with the bilayer, changing its orientation. The interaction of membrane-bound auxilin-1 PTEN-like domain with negatively charged lipid headgroups results in nanoclustering of PIP₂ molecules in the adjacent bilayer leaflet.     

Effect of gibberellic acid and salicylic acid on the activity and electrophoretic pattern of IAA-oxidase during floral induction inImpatiens balsamina

IAA-oxidase activity increased in the stem as well as in the leaves of plants treated with GA₃, SA and GA₃ + SA during the early stages under inductive and non-inductive photoperiods, the activity being the highest in GA₃ + SA-treated plants. An isoenzyme of IAA-oxidase with Rm 0.15 developed in the stem as well as in the leaves subsequent to 1 or 2 inductive treatments. As this band persisted till the end of the experiment, it may be associated with the initiation as well as development of floral buds. Another band (Rm 0.30) appears to be associated with the phenol (SA) as it developed in the stem as well as in the leaves of SA- and GA₃ + SA-treated plants under both photoperiods. A band with Rm 0.60 developed in the leaves but not in the stem of GA₃-, SA- and GA₃ + SA-treated plants under both photoperiods.     

Methyl jasmonate and salicylic acid elicitation induces ginsenosides accumulation, enzymatic and non-enzymatic antioxidant in suspension culture Panax ginseng roots in bioreactors

The effects of methyl jasmonate (MJ) and salicylic acid (SA) on changes of the activities of major antioxidant enzymes, superoxide anion accumulation (O₂ ⁻), ascorbate, total glutathione (TG), malondialdehyde (MDA) content and ginsenoside accumulation were investigated in ginseng roots (Panax ginseng L.) in 4 l (working volume) air lift bioreactors. Single treatment of 200 μM MJ and SA to P. ginseng roots enhanced ginsenoside accumulation compared to the control and harvested 3, 5, 7 and 9 days after treatment. MJ and SA treatment induced an oxidative stress in P. ginseng roots, as shown by an increase in lipid peroxidation due to rise in O₂ ⁻ accumulation. Activity of superoxide dismutase (SOD) was inhibited in MJ-treated roots, while the activities of monodehydroascorbate reductase (MDHAR), dehydroascorbate reductase (DHAR), SOD, guaiacol peroxidase (G-POD), glutathione peroxidase (GPx) and glutathione reductase (GR) were induced in SA-treated roots. A strong decrease in the activity of catalase (CAT) was obtained in both MJ- and SA-treated roots. Activities of ascorbate peroxidase (APX) and glutathione S transferase (GST) were higher in MJ than SA while the contents of reduced ascorbate (ASC), redox state (ASC/(ASC+DHA)) and TG were higher in SA- than MJ-treated roots while oxidized ascorbate (DHA) decreased in both cases. The result of these analyses suggests that roots are better protected against the O₂ ⁻ stress, thus mitigating MJ and SA stress. The information obtained in this work is useful for efficient large-scale production of ginsenoside by plant-root cultures.     

Enzyme-catalyzed hydrolysis of the supported phospholipid bilayers studied by atomic force microscopy

Atomic force microscopy (AFM) is employed to reveal the morphological changes of the supported phospholipid bilayers hydrolyzed by a phospholipase A₂ (PLA₂) enzyme in a buffer solution at room temperature. Based on the high catalytic selectivity of PLA₂ toward l-enantiomer phospholipids, five kinds of supported bilayers made of l- and d-dipalmitoylphosphatidylcholines (DPPC), including l-DPPC (upper leaflet adjacent to solution)/l-DPPC (bottom leaflet) (or l/l in short), l/d, d/l, d/d, and racemic ld/ld, were prepared on a mica surface in gel-phase, to explicate the kinetics and mechanism of the enzyme-induced hydrolysis reaction in detail. AFM observations for the l/l bilayer show that the hydrolysis rate for l-DPPC is significantly increased by PLA₂ and most of the hydrolysis products desorb from substrate surface in 40 min. As d-enantiomers are included in the bilayer, the hydrolysis rate is largely decreased in comparison with the l/l bilayer. The time used to hydrolyze the as-prepared bilayers by PLA₂ increases in the sequence of l/l, l/d, ld/ld, and d/l (d/d is inert to the enzyme action). d-enantiomers in the enantiomer hybrid bilayers remain on the mica surface at the end of the hydrolysis reaction. It was confirmed that the hydrolysis reaction catalyzed by PLA₂ preferentially occurs at the edges of pits or defects on the bilayer surface. The bilayer structures are preserved during the hydrolysis process. Based on these observations, a novel kinetics model is proposed to quantitatively account for the PLA₂-catalyzed hydrolysis of the supported phospholipid bilayers. The model simulation demonstrates that PLA₂ mainly binds with lipids at the perimeter of defects in the upper leaflet and leads to a hydrolysis reaction, yielding species soluble to the solution phase. The lipid molecules underneath subsequently flip up to the upper leaflet to maintain the hydrophilicity of the bilayer structure. Our analysis shows that d-enantiomers in the hybrid bilayers considerably reduce the hydrolysis rate by its ineffective binding with PLA₂.     

E. coli Outer Membrane and Interactions with OmpLA

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.     

E.?coli Outer Membrane and Interactions with OmpLA

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.     

Simulations of membranes and other interfacial systems using P2(1) and Pc periodic boundary conditions

We demonstrate the ease and utility of simulating heterogeneous interfacial systems with P2(1) and Pc periodic boundary conditions which allow, for example, lipids in a membrane to switch leaflets. In preliminary tests, P2(1) was shown to yield equivalent results to P1 in simulations of bulk water, a water/vacuum interface, and pure DPPC bilayers with an equal number of lipids per leaflet; equivalence of Pc and P1 was also demonstrated for the former two systems. P2(1) was further tested in simulations involving the spreading of an octane film on water, and equilibration of a DPPC bilayer from an initial condition containing different numbers of lipids in the two leaflets. Lastly, a simulation in P2(1) of a DOPC/melittin membrane showed significant passage of lipids to the melittin-containing leaflet from the initial distribution, and lends insight into the condensation of lipids by melittin.     

Acyl-chain methyl distributions of liquid-ordered and -disordered membranes

A central feature of the lipid raft concept is the formation of cholesterol-rich lipid domains. The introduction of relatively rigid cholesterol molecules into fluid liquid-disordered (L_d) phospholipid bilayers can produce liquid-ordered (L_o) mixtures in which the rigidity of cholesterol causes partial ordering of the flexible hydrocarbon acyl chains of the phospholipids. Several lines of evidence support this concept, but direct structural information about L_o membranes is lacking. Here we present the structure of L_o membranes formed from cholesterol and dioleoylphosphatidylcholine (DOPC). Specific deuteration of the DOPC acyl-chain methyl groups and neutron diffraction measurements reveal an extraordinary disorder of the acyl chains of neat L_d DOPC bilayers. The disorder is so great that >20% of the methyl groups are in intimate contact with water in the bilayer interface. The ordering of the DOPC acyl chains by cholesterol leads to retraction of the methyl groups away from the interface. Molecular dynamics simulations based on experimental systems reveal asymmetric transbilayer distributions of the methyl groups associated with each bilayer leaflet.     

Regulation of TRPV1 Ion Channel by Phosphoinositide (4,5)-Bisphosphate: THE ROLE OF MEMBRANE ASYMMETRY*

Membrane asymmetry is essential for generating second messengers that act in the cytosol and for trafficking of membrane proteins and membrane lipids, but the role of asymmetry in regulating membrane protein function remains unclear. Here we show that the signaling lipid phosphoinositide 4,5-bisphosphate (PI(4,5)P₂) has opposite effects on the function of TRPV1 ion channels depending on which leaflet of the cell membrane it resides in. We observed potentiation of capsaicin-activated TRPV1 currents by PI(4,5)P₂ in the intracellular leaflet of the plasma membrane but inhibition of capsaicin-activated currents when PI(4,5)P₂ was in both leaflets of the membrane, although much higher concentrations of PI(4,5)P₂ in the extracellular leaflet were required for inhibition compared with the concentrations of PI(4,5)P₂ in the intracellular leaflet that produced activation. Patch clamp fluorometry using a synthetic PI(4,5)P₂ whose fluorescence reports its concentration in the membrane indicates that PI(4,5)P₂ must incorporate into the extracellular leaflet for its inhibitory effects to be observed. The asymmetry-dependent effect of PI(4,5)P₂ may resolve the long standing controversy about whether PI(4,5)P₂ is an activator or inhibitor of TRPV1. Our results also underscore the importance of membrane asymmetry and the need to consider its influence when studying membrane proteins reconstituted into synthetic bilayers.     

A molecular dynamics study of the response of lipid bilayers and monolayers to trehalose

Surface tensions evaluated from molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine bilayers and monolayers at surface areas/lipid of 54, 64, and 80 A2 are uniformly lowered 4-8 dyn/cm upon addition of trehalose in a 1:2 trehalose/lipid ratio. Constant surface tension simulations of bilayers yield the complementary result: an increase in surface area consistent with the surface pressure-surface area (pi-A) isotherms. Hydrogen bonding by trehalose, replacement of waters in the headgroup region, and modulation of the dipole potential are all similar in bilayers and monolayers at the same surface area. These results strongly support the assumption that experimental measurements on the interactions of surface active components such as trehalose with monolayers can yield quantitative insight to their effects on bilayers. The simulations also indicate that the 20-30 dyn/cm difference in surface tension of the bilayer leaflet and monolayer arises from differences in the chain regions, not the headgroup/water interfaces.     

Development of a CP P NMR Broadline Simulation Methodology for Studying the Interactions of Antihypertensive AT1 Antagonist Losartan with Phospholipid Bilayers

A cross-polarization (CP) ³¹P NMR broadline simulation methodology was developed for studying the effects of drugs in phospholipids bilayers. Based on seven-parameter fittings, this methodology provided information concerning the conformational changes and dynamics effects of losartan in the polar region of the dipalmitoylphosphatidylcholine bilayers. The test molecule for this study was losartan, an antihypertensive drug known to exert its effect on AT₁ transmembrane receptors. The results were complemented and compared with those of differential scanning calorimetry, solid-state ¹³C NMR spectroscopy, Raman spectroscopy, and electron spin resonance. More specifically, these physical chemical methodologies indicated that the amphipathic losartan molecule interacts with the hydrophilic-head zone of the lipid bilayers. The CP ³¹P NMR broadline simulations showed that the lipid molecules in the bilayers containing losartan displayed greater collective tilt compared to the tilt displayed by the load-free bilayers, indicating improved packing. The Raman results displayed a decrease in the trans/gauche ratio and increased intermolecular interactions of the acyl chains in the liquid crystalline phase. Additional evidence, suggesting that losartan possibly anchors in the realm of the headgroup, was derived from upfield shift of the average chemical shift σ^iso of the ³¹P signal in the presence of losartan and from shift of the observed peak at 715 cm⁻¹ attributed to C-N stretching in the Raman spectra.     

Water permeability of asymmetric planar lipid bilayers: leaflets of different composition offer independent and additive resistances to permeation

To understand how plasma membranes may limit water flux, we have modeled the apical membrane of MDCK type 1 cells. Previous experiments demonstrated that liposomes designed to mimic the inner and outer leaflet of this membrane exhibited 18-fold lower water permeation for outer leaflet lipids than inner leaflet lipids (Hill, W.G., and M.L. Zeidel. 2000. J. Biol. Chem. 275:30176-30185), confirming that the outer leaflet is the primary barrier to permeation. If leaflets in a bilayer resist permeation independently, the following equation estimates single leaflet permeabilities: 1/P(AB) = 1/P(A) + 1/P(B) (Eq. l), where P(AB) is the permeability of a bilayer composed of leaflets A and B, P(A) is the permeability of leaflet A, and P(B) is the permeability of leaflet B. Using for the MDCK leaflet-specific liposomes gives an estimated value for the osmotic water permeability (P(f)) of 4.6 x 10(-4) cm/s (at 25 degrees C) that correlated well with experimentally measured values in intact cells. We have now constructed both symmetric and asymmetric planar lipid bilayers that model the MDCK apical membrane. Water permeability across these bilayers was monitored in the immediate membrane vicinity using a Na+-sensitive scanning microelectrode and an osmotic gradient induced by addition of urea. The near-membrane concentration distribution of solute was used to calculate the velocity of water flow (Pohl, P., S.M. Saparov, and Y.N. Antonenko. 1997. Biophys. J. 72:1711-1718). At 36 degrees C, P(f) was 3.44 +/- 0.35 x 10(-3) cm/s for symmetrical inner leaflet membranes and 3.40 +/- 0.34 x 10(-4) cm/s for symmetrical exofacial membranes. From, the estimated permeability of an asymmetric membrane is 6.2 x 10(-4) cm/s. Water permeability measured for the asymmetric planar bilayer was 6.7 +/- 0.7 x 10(-4) cm/s, which is within 10% of the calculated value. Direct experimental measurement of P(f) for an asymmetric planar membrane confirms that leaflets in a bilayer offer independent and additive resistances to water permeation and validates the use of.     

Molecular studies of the gel to liquid-crystalline phase transition for fully hydrated DPPC and DPPE bilayers

Molecular dynamics simulations were used for a comprehensive study of the structural properties of saturated lipid bilayers, DPPC and DPPE, near the main phase transition. Though the chemical structure of DPPC and DPPE are largely similar (they only differ in the choline and ethanolamine groups), their transformation process from a gel to a liquid-crystalline state is contrasting. For DPPC, three distinct structures can be identified relative to the melting temperature (T_m): below T_m with “mixed” domains consisting of lipids that are tilted with partial overlap of the lipid tails between leaflet; near T_m with a slight increase in the average area per lipid, resulting in a rearrangement of the lipid tails and an increase in the bilayer thickness; and above T_m with unhindered lipid tails in random motion resulting in an increase in %gauche formed and increase in the level of interdigitation between lipid leaflets. For DPPE, the structures identified were below T_m with “ordered” domains consisting of slightly tilted lipid tails and non-overlapping lipid tails between leaflets, near T_m with minimal rearrangement of the lipids as the bilayer thickness reduces slightly with increasing temperature, and above T_m with unhindered lipid tails as that for DPPC. For DPPE, most of the lipid tails do not overlap as observed to DPPC, which is due to the tight packing of the DPPE molecules. The non-overlapping behavior of DPPE above T_m is confirmed from the density profile of the terminal carbon atoms in each leaflet, which shows a narrow distribution near the center of the bilayer core. This study also demonstrates that atomistic simulations are capable of capturing the phase transition behavior of lipid bilayers, providing a rich set of molecular and structural information at and near the transition state.     

Constant surface tension molecular dynamics simulations of lipid bilayers with trehalose

Surface areas and fluctuations evaluated from 50 ns molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers in a 1:2 trehalose:lipid ratio carried out at surface tensions 10, 17 and 25 dyn/cm/leaflet are compared with those of pure bilayers under the same conditions. Trehalose increases the surface area, as consistent with the surface tension lowering observed in simulations at constant area. The system bulk elastic modulus K _b  = 1.5 ± 0.3 × 10¹⁰ dyn/cm². It is independent of bilayer surface area and trehalose content within statistical error. In contrast, the area elastic modulus K _a shows a strong area dependence. At 64 Å²/lipid (the experimental surface area), K _a  = 138 ± 26 dyn/cm for a pure DPPC bilayer and 82 ± 10 dyn/cm for one with trehalose; i.e. trehalose increases fluidity of the bilayer surface at this area per lipid.     

Transmembrane cholesterol migration in planar lipid membranes measured with Vibrio cholerae cytolysin as molecular tool

The rate of transbilayer movement (flip-flop) of cholesterol was estimated using planar bilayers with defined initial asymmetry, formed by the opposing monolayers technique. Vibrio cholerae cytolysin (VCC) was utilized as a molecular tool for measuring the cholesterol concentration in the cis leaflet of asymmetric bilayers. To quantify cholesterol flip-flop in planar lipid bilayers, a mathematical model was developed. It considers both the lateral diffusion rate of cholesterol within each monolayer and the flip-flop rate. The difference in initial and steady-state cholesterol contents in bilayer leaflets was used as a start point. Assuming the lateral diffusion coefficient to be of 1 × 10⁻⁸ cm² s⁻¹, the characteristic time of cholesterol flip-flop at 25 ± 2 °C was estimated as <10 s.     

Salicylic acid antagonism of EDS1‐driven cell death is important for immune and oxidative stress responses in Arabidopsis

Reactive oxygen species (ROS) have emerged as signals in the responses of plants to stress. Arabidopsis Enhanced Disease Susceptibility1 (EDS1) regulates defense and cell death against biotrophic pathogens and controls cell death propagation in response to chloroplast‐derived ROS. Arabidopsis Nudix hydrolase7 (nudt7) mutants are sensitized to photo‐oxidative stress and display EDS1‐dependent enhanced resistance, salicylic acid (SA) accumulation and initiation of cell death. Here we explored the relationship between EDS1, EDS1‐regulated SA and ROS by examining gene expression profiles, photo‐oxidative stress and resistance phenotypes of nudt7 mutants in combination with eds1 and the SA‐biosynthetic mutant, sid2. We establish that EDS1 controls steps downstream of chloroplast‐derived O₂^?? that lead to SA‐assisted H₂O₂ accumulation as part of a mechanism limiting cell death. A combination of EDS1‐regulated SA‐antagonized and SA‐promoted processes is necessary for resistance to host‐adapted pathogens and for a balanced response to photo‐oxidative stress. In contrast to SA, the apoplastic ROS‐producing enzyme NADPH oxidase RbohD promotes initiation of cell death during photo‐oxidative stress. Thus, chloroplastic O₂^?? signals are processed by EDS1 to produce counter‐balancing activities of SA and RbohD in the control of cell death. Our data strengthen the idea that EDS1 responds to the status of O₂^?? or O₂^??‐generated molecules to coordinate cell death and defense outputs. This activity may enable the plant to respond flexibly to different biotic and abiotic stresses in the environment.     

Coupled Diffusion of Peripherally Bound Peptides along the Outer and Inner Membrane Leaflets

Transmembrane signaling implies that peripheral protein binding to one leaflet be detected by the opposite leaflet. Therefore, protein recruitment into preexisting cholesterol and sphingolipid rich platforms may be required. However, no clear molecular picture has evolved about how these rafts in both leaflets are connected. By using planar lipid bilayers, we show that the peripheral binding of a charged molecule (poly-lysine, PLL) is detected at the other side of the bilayer without involvement of raft lipids. The diffusion coefficient, D_P, of PLL differed by a factor of √2 when PLL absorbed to one or to both leaflets of planar membranes. Fluorescence correlation spectroscopy showed that the changes of the lipid diffusion coefficient, D_M, were even more pronounced. Although D_M remained larger than D_P on PLL binding to the first membrane leaflet, D_M dropped to D_P on PLL binding to both leaflets, which indicated that the lipids sandwiched between two PLL molecules had formed a nanodomain. Due to its small area of ∼20 nm² membrane electrostriction or leaflet interaction at bilayer midplane can only make a small contribution to interleaflet coupling. The tendency of the system to maximize the area where the membrane is free to undulate seems to be more important. As a spot with increased bending stiffness, the PLL bound patch in one leaflet attracts a stiffening additive on the other leaflet. That is to say, instead of suppressing undulations in two spots, two opposing PLL molecules migrate along a membrane at matching positions and suppress these undulations in a single spot. The gain in undulation energy is larger than the energy required for the alignment of two small PLL domains in opposite leafs and their coordinated diffusion. We propose that this type of mechanical interaction between two membrane separated ligands generally contributes to transmembrane signaling.     

Spontaneous Curvature,Differential Stress,and Bending Modulus of Asymmetric Lipid Membranes

Lipid bilayers can exhibit asymmetric states, in which the physical characteristics of one leaflet differ from those of the other. This most visibly manifests in a different lipid composition, but it can also involve opposing lateral stresses in each leaflet that combine to an overall vanishing membrane tension. Here, we use theoretical modeling and coarse-grained simulation to explore the interplay between a compositional asymmetry and a nonvanishing differential stress. Minimizing the total elastic energy leads to a preferred spontaneous curvature that balances torques due to both bending moments and differential stress, with sometimes unexpected consequences. For instance, asymmetric flat bilayers, whose specific areas in each leaflet are matched to those of corresponding tensionless symmetric flat membranes, still exhibit a residual differential stress because the conditions of vanishing area strain and vanishing bending moment differ. We also measure the curvature rigidity of asymmetric bilayers and find that a sufficiently strong differential stress, but not compositional asymmetry alone, can increase the bending modulus. The likely cause is a stiffening of the compressed leaflet, which appears to be related to its gel transition but not identical with it. We finally show that the impact of cholesterol on differential stress depends on the relative strength of elastic and thermodynamic driving forces: if cholesterol solvates equally well in both leaflets, it will redistribute to cancel both leaflet tensions almost completely, but if its partitioning free energy prefers one leaflet over the other, the resulting distribution bias may even create differential stress. Because cells keep most of their lipid bilayers in an asymmetric nonequilibrium steady state, our findings suggest that biomembranes are elastically more complex than previously thought: besides a spontaneous curvature, they might also exhibit significant differential stress, which could strongly affect their curvature energetics.     

Molecular Dynamics Simulations of the Bacterial UraA H+-Uracil Symporter in Lipid Bilayers Reveal a Closed State and a Selective Interaction with Cardiolipin

The Escherichia coli UraA H⁺-uracil symporter is a member of the nucleobase/ascorbate transporter (NAT) family of proteins, and is responsible for the proton-driven uptake of uracil. Multiscale molecular dynamics simulations of the UraA symporter in phospholipid bilayers consisting of: 1) 1-palmitoyl 2-oleoyl-phosphatidylcholine (POPC); 2) 1-palmitoyl 2-oleoyl-phosphatidylethanolamine (POPE); and 3) a mixture of 75% POPE, 20% 1-palmitoyl 2-oleoyl-phosphatidylglycerol (POPG); and 5% 1-palmitoyl 2-oleoyl-diphosphatidylglycerol/cardiolipin (CL) to mimic the lipid composition of the bacterial inner membrane, were performed using the MARTINI coarse-grained force field to self-assemble lipids around the crystal structure of this membrane transport protein, followed by atomistic simulations. The overall fold of the protein in lipid bilayers remained similar to the crystal structure in detergent on the timescale of our simulations. Simulations were performed in the absence of uracil, and resulted in a closed state of the transporter, due to relative movement of the gate and core domains. Anionic lipids, including POPG and especially CL, were found to associate with UraA, involving interactions between specific basic residues in loop regions and phosphate oxygens of the CL head group. In particular, three CL binding sites were identified on UraA: two in the inner leaflet and a single site in the outer leaflet. Mutation of basic residues in the binding sites resulted in the loss of CL binding in the simulations. CL may play a role as a “proton trap” that channels protons to and from this transporter within CL-enriched areas of the inner bacterial membrane.