Membrane Bending Moduli of Coexisting Liquid Phases Containing Transmembrane Peptide

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

Binary bilayer simulations for partitioning within membranes

Membrane proteins (MPs) often show preference for one phase over the other, which is characterized by the partition coefficient, Kp. The physical mechanisms underlying Kp have been only inferred indirectly from experiments due to the unavailability of detailed structures and compositions of ordered phases. Molecular dynamics (MD) simulations can complement these details and thus, in principle, provide further insights into the partitioning of MPs between two phases. However, the application of MD has remained difficult due to long time scales required for equilibration and large system size for the phase stability, which have not been fully resolved even in free energy simulations. This chapter describes the recently developed binary bilayer simulation method, where the membrane is composed of two laterally attached membrane patches. The binary bilayer system (BBS) is designed to preserve the lateral packing of both phases in a significantly smaller size compared to that required for macroscopic phase separation. These characteristics are advantageous in partitioning simulations, as the length scale for diffusion across the system can be significantly smaller. Hence the BBS can be efficiently employed in both conventional MD and free energy simulations, though sampling in ordered phases remains difficult due to slow diffusion. Development of efficient lipid swapping methods and its combination with the BBS would be a useful approach for partitioning in coexisting phases.

Unsaturated Lipids Facilitate Partitioning of Transmembrane Peptides into the Liquid Ordered Phase

The affinity of single-pass transmembrane (TM) proteins for ordered membrane phases has been reported to depend on their lipidation, TM length, and lipid accessible surface area. In this work, the raft affinities of the TM domain of the linker for activation of T cells and its depalmitoylated variant are assessed using free energy simulations in a binary bilayer system composed of two laterally patched bilayers of ternary liquid ordered (Lo) and liquid disordered (Ld) phases. These phases are modeled by distinct compositions of distearoylphosphatidylcholine, palmitoyloleoylphosphatidylcholine (POPC), and cholesterol, and the simulations were carried out for 4.5 μs/window. Both peptides are shown to preferentially partition into the Ld phase in agreement with model membrane experiments and previous simulations on ternary lipid mixtures but not with measurements on giant plasma membrane vesicles where the Lo is slightly preferred. However, the 500 ns average relaxation time of lipid rearrangement around the peptide precluded a quantitative analysis of free energy differences arising from peptide palmitoylation and two different lipid compositions. When in the Lo phase, peptides reside in regions rich in POPC and interact preferentially with its unsaturated tail. Hence, the detailed substructure of the Lo phase is an important modulator of peptide partitioning, in addition to the inherent properties of the peptide.

Impact of transmembrane peptides on individual lipid motions and collective dynamics of lipid bilayers

The fluid nature of lipid bilayers is indispensable for the dynamic regulation of protein function and membrane morphology in biological membranes. Membrane-spanning domains of proteins interact with surrounding lipids and alter the physical properties of lipid bilayers. However, there is no comprehensive view of the effects of transmembrane proteins on the membrane's physical properties. Here, we investigated the effects of transmembrane peptides with different flip-flop-promoting abilities on the dynamics of a lipid bilayer employing complemental fluorescence and neutron scattering techniques. The quasi-elastic neutron scattering and fluorescence experiments revealed that lateral diffusion of the lipid molecules and the acyl chain motions were inhibited by the inclusion of transmembrane peptides. The neutron spin-echo spectroscopy measurements indicated that the lipid bilayer became more rigid but more compressible and the membrane viscosity increased when the transmembrane peptides were incorporated into the membrane. These results suggest that the inclusion of rigid transmembrane structures hinders individual and collective lipid motions by slowing down lipid diffusion and increasing interleaflet coupling. The present study provides a clue for understanding how the local interactions between lipids and proteins change the collective dynamics of the lipid bilayers, and therefore, the function of biological membranes.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

The Effects of Cholesterol Oxidation on Erythrocyte Plasma Membranes: A Monolayer Study

Cholesterol plays a key role in the molecular and mesoscopic organisation of lipid membranes and it is expected that changes in its molecular structure (e.g., through environmental factors such as oxidative stress) may affect adversely membrane properties and function. In this study, we present evidence that oxidation of cholesterol has significant effects on the mechanical properties, molecular and mesoscopic organisation and lipid-sterol interactions in condensed monolayers composed of the main species found in the inner leaflet of the erythrocyte membrane. Using a combination of experimental methods (static area compressibility, surface dilatational rheology, fluorescence microscopy, and surface sensitive X-ray techniques) and atomistic molecular dynamics simulations, we show that oxidation of cholesterol to 7-ketocholesterol leads to stiffening of the monolayer (under both static and dynamic conditions), significant changes in the monolayer microdomain organisation, disruption in the van der Waals, electrostatic and hydrophobic interactions between the sterol and the other lipid species, and the lipid membrane hydration. Surface sensitive X-ray techniques reveal that, whilst the molecular packing mode is not significantly affected by cholesterol oxidation in these condensed phases, there are subtle changes in membrane thickness and a significant decrease in the coherence length in monolayers containing 7-ketocholesterol.

Curved membrane structures induced by native lipids in giant vesicles

Native lipids in cell-membrane support crucial functions like intercell communication via their ability to deform into curved membrane structures. Cell membrane mimicking Giant unilamellar vesicles (GUV) is imperative in understanding native lipid's role in membrane transformation however remains challenging to assemble. We construct two giant vesicle models mimicking bacterial inner-membrane (IM) and outer-membrane (OM) under physiological conditions using single-step gel-assisted lipid swelling. IM vesicles composed of native bacterial lipids undergo small-scale membrane remodeling into bud and short-nanotube structures. In contrast, OM vesicles asymmetrically assembled from Lipopolysaccharide (LPS) and bacterial lipids underwent global membrane deformation under controlled osmotic stress. Remarkably, highly-curved structures mimicking cell-membrane architectures, including daughter vesicle networks interconnected by necks and nano-tubes ranging from micro to nanoscale, are generated in OM vesicles at osmotic stress comparable to that applied in IM vesicles. Further, we provide a quantitative description of the membrane structures by experimentally determining membrane elastic parameters, i.e., neck curvature and bending rigidity. We can conclude that a larger spontaneous curvature estimated from the neck curvature and softer membranes in OM vesicles is responsible for large-scale deformation compared to IM vesicles. Our findings will help comprehend the shape dynamics of complex native bacterial lipid membranes.

The influence of phosphatidylserine localisation and lipid phase on membrane remodelling by the ESCRT-II/ESCRT-III complex

The endosomal sorting complex required for transport (ESCRT) organises in supramolecular structures on the surface of lipid bilayers to drive membrane invagination and scission of intraluminal vesicles (ILVs), a process also controlled by membrane mechanics. However, ESCRT association with the membrane is also mediated by electrostatic interactions with anionic phospholipids. Phospholipid distribution within natural biomembranes is inhomogeneous due to, for example, the formation of lipid rafts and curvature-driven lipid sorting. Here, we have used phase-separated giant unilamellar vesicles (GUVs) to investigate the link between phosphatidylserine (PS)-rich lipid domains and ESCRT activity. We employ GUVs composed of phase separating lipid mixtures, where unsaturated DOPS and saturated DPPS lipids are incorporated individually or simultaneously to enhance PS localisation in liquid disordered (L_d) and/or liquid ordered (L_o) domains, respectively. PS partitioning between the coexisting phases is confirmed by a fluorescent Annexin V probe. Ultimately, we find that ILV generation promoted by ESCRTs is significantly enhanced when PS lipids localise within L_d domains. However, the ILVs that form are rich in L_o lipids. We interpret this surprising observation as preferential recruitment of the L_o phase beneath the ESCRT complex due to its increased rigidity, where the L_d phase is favoured in the neck of the resultant buds to facilitate the high membrane curvature in these regions of the membrane during the ILV formation process. L_d domains offer lower resistance to membrane bending, demonstrating a mechanism by which the composition and mechanics of membranes can be coupled to regulate the location and efficiency of ESCRT activity.

The membrane transporter lactose permease increases lipid bilayer bending rigidity

Cellular life relies on membranes, which provide a resilient and adaptive cell boundary. Many essential processes depend upon the ease with which the membrane is able to deform and bend, features that can be characterized by the bending rigidity. Quantitative investigations of such mechanical properties of biological membranes have primarily been undertaken in solely lipid bilayers and frequently in the absence of buffers. In contrast, much less is known about the influence of integral membrane proteins on bending rigidity under physiological conditions. We focus on an exemplar member of the ubiquitous major facilitator superfamily of transporters and assess the influence of lactose permease on the bending rigidity of lipid bilayers. Fluctuation analysis of giant unilamellar vesicles (GUVs) is a useful means to measure bending rigidity. We find that using a hydrogel substrate produces GUVs that are well suited to fluctuation analysis. Moreover, the hydrogel method is amenable to both physiological salt concentrations and anionic lipids, which are important to mimic key aspects of the native lactose permease membrane. Varying the fraction of the anionic lipid in the lipid mixture DOPC/DOPE/DOPG allows us to assess the dependence of membrane bending rigidity on the topology and concentration of an integral membrane protein in the lipid bilayer of GUVs. The bending rigidity gradually increases with the incorporation of lactose permease, but there is no further increase with greater amounts of the protein in the membrane.

Single-molecule fluorescence vistas of how lipids regulate membrane proteins

The study of membrane proteins is undergoing a golden era, and we are gaining unprecedented knowledge on how this key group of proteins works. However, we still have only a basic understanding of how the chemical composition and the physical properties of lipid bilayers control the activity of membrane proteins. Single-molecule (SM) fluorescence methods can resolve sample heterogeneity, allowing to discriminate between the different molecular populations that biological systems often adopt. This short review highlights relevant examples of how SM fluorescence methodologies can illuminate the different ways in which lipids regulate the activity of membrane proteins. These studies are not limited to lipid molecules acting as ligands, but also consider how the physical properties of the bilayer can be determining factors on how membrane proteins function.

Improved growth of Escherichia coli in aminoglycoside antibiotics by the zor-orz toxin-antitoxin system

Type I toxin-antitoxin systems consist of a small protein (under 60 amino acids) whose overproduction can result in cell growth stasis or death, and a small RNA that represses translation of the toxin mRNA. Despite their potential toxicity, type I toxin proteins are increasingly linked to improved survival of bacteria in stressful environments and antibiotic persistence. While the interaction of toxin mRNAs with their cognate antitoxin sRNAs in some systems are well characterized, additional translational control of many toxins and their biological roles are not well understood. Using an ectopic overexpression system, we show that the efficient translation of a chromosomally encoded type I toxin, ZorO, requires mRNA processing of its long 5' untranslated region (UTR; Δ28 UTR). The severity of ZorO induced toxicity on growth inhibition, membrane depolarization, and ATP depletion were significantly increased if expressed from the Δ28 UTR versus the full-length UTR. ZorO did not form large pores as evident via a liposomal leakage assay, in vivo morphological analyses, and measurement of ATP loss. Further, increasing the copy number of the entire zor-orz locus significantly improved growth of bacterial cells in the presence of kanamycin and increased the minimum inhibitory concentration against kanamycin and gentamycin; however, no such benefit was observed against other antibiotics. This supports a role for the zor-orz locus as a protective measure against specific stress agents and is likely not part of a general stress response mechanism. Combined, these data shed more insights into the possible native functions for type I toxin proteins. IMPORTANCE Bacterial species can harbor gene pairs known as type I toxin-antitoxin systems where one gene encodes a small protein that is toxic to the bacteria producing it and a second gene that encodes a small RNA antitoxin to prevent toxicity. While artificial overproduction of type I toxin proteins can lead to cell growth inhibition and cell lysis, the endogenous translation of type I toxins appears to be tightly regulated. Here, we show translational regulation controls production of the ZorO type I toxin and prevents subsequent negative effects on the cell. Further, we demonstrate a role for zorO and its cognate antitoxin in improved growth of E. coli in the presence of aminoglycoside antibiotics.

The Effect of Phosphatidylserine on a pH-Responsive Peptide Is Defined by Its Noninserting End

The acidity-triggered rational membrane (ATRAM) peptide was designed to target acidic diseases such as cancer. An acidic extracellular medium, such as that found in aggressive tumors, drives the protonation of the glutamic acids in ATRAM, leading to the membrane translocation of its C-terminus and the formation of a transmembrane helix. Compared to healthy cells, cancerous cells often increase exposure of the negatively charged phosphatidylserine (PS) on the outer leaflet of the plasma membrane. Here we use a reconstituted vesicle system to explore how PS influences the interaction of ATRAM with membranes. To explore this, we used two new variants of ATRAM, termed K2-ATRAM and Y-ATRAM, with small modifications at the noninserting N-terminus. We observed that the effect of PS on the membrane insertion pK and lipid partitioning hinged on the sequence of the noninserting end. Our data additionally indicate that the effect of PS on the insertion pK does not merely depend on electrostatics, but it is multifactorial. Here we show how small sequence changes can impact the interaction of a peptide with membranes of mixed lipid composition. These data illustrate how model studies using neutral bilayers, which do not mimic the negative charge found in the plasma membrane of cancer cells, may fail to capture important aspects of the interaction of anticancer peptides with tumor cells. This information can guide the design of therapeutic peptides that target the acidic environments of different diseased tissues.

Curvature-driven adsorption of cationic nanoparticles to phase boundaries in multicomponent lipid bilayers

Understanding the interactions between surface-functionalized gold nanoparticles (NPs) and lipid bilayers is necessary to guide the design of NPs for biomedical applications. Recent experiments found that cationic NPs adsorb more strongly to phase-separated multicomponent lipid bilayers than single-component liquid-disordered bilayers, suggesting that phase separation affects NP-bilayer interactions. In this work, we use coarse-grained molecular dynamics simulations to investigate the effect of lipid phase behavior on the adsorption of small cationic NPs. We first determined the free energy change for adsorbing a NP to one-phase liquid-disordered and one-phase liquid-ordered bilayers. The simulations indicate that NP adsorption depends on a competition between favorable NP-lipid interactions and the unfavorable curvature deformation of the bilayer, resulting in stronger interactions with the liquid-disordered bilayer due to its lower bending modulus. We then measured the free energy change associated with moving a NP across the surface of a phase-separated bilayer and identified a free energy minimum at the phase boundary. The free energy minimum is attributed to the thickness gradient between the two phases that enables favorable NP-lipid interactions without necessitating large curvature deformations. The simulation results thus indicate that the intrinsic curvature present at phase boundaries drives preferential interactions with surface-adsorbed NPs, providing new insight into the forces that drive NP behavior at multicomponent, phase-separated lipid bilayers.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Transmembrane domain dimerization induces cholesterol rafts in curved lipid bilayers

Are the dimerization of transmembrane (TM) domains and the reorganization of the lipid bilayer two independent events?

A Rationale for Mesoscopic Domain Formation in Biomembranes

Cell plasma membranes display a dramatically rich structural complexity characterized by functional sub-wavelength domains with specific lipid and protein composition. Under favorable experimental conditions, patterned morphologies can also be observed in vitro on model systems such as supported membranes or lipid vesicles. Lipid mixtures separating in liquid-ordered and liquid-disordered phases below a demixing temperature play a pivotal role in this context. Protein-protein and protein-lipid interactions also contribute to membrane shaping by promoting small domains or clusters. Such phase separations displaying characteristic length-scales falling in-between the nanoscopic, molecular scale on the one hand and the macroscopic scale on the other hand, are named mesophases in soft condensed matter physics. In this review, we propose a classification of the diverse mechanisms leading to mesophase separation in biomembranes. We distinguish between mechanisms relying upon equilibrium thermodynamics and those involving out-of-equilibrium mechanisms, notably active membrane recycling. In equilibrium, we especially focus on the many mechanisms that dwell on an up-down symmetry breaking between the upper and lower bilayer leaflets. Symmetry breaking is an ubiquitous mechanism in condensed matter physics at the heart of several important phenomena. In the present case, it can be either spontaneous (domain buckling) or explicit, i.e., due to an external cause (global or local vesicle bending properties). Whenever possible, theoretical predictions and simulation results are confronted to experiments on model systems or living cells, which enables us to identify the most realistic mechanisms from a biological perspective.