Effect of membrane tension on the physical properties of DOPC lipid bilayer membrane

Assessing Membrane Fluidity and Visualizing Fluid Membrane Domains in Bacteria Using Fluorescent Membrane Dyes

Membrane fluidity is a key parameter of bacterial membranes that undergoes quick adaptation in response to environmental challenges and has recently emerged as an important factor in the antibacterial mechanism of membrane-targeting antibiotics. The specific level of membrane fluidity is not uniform across the bacterial cell membrane. Rather, specialized microdomains associated with different cellular functions can exhibit fluidity values that significantly deviate from the average. Assessing changes in the overall membrane fluidity and formation of membrane microdomains is therefore pivotal to understand both the functional organization of the bacterial cell membrane as well as antibiotic mechanisms. Here we describe how two fluorescent membrane dyes, laurdan and DiIC12, can be employed to assess membrane fluidity in living bacteria. We focus on Bacillus subtilis, since this organism has been relatively well-studied with respect to membrane domains. However, we also describe how these assays can be adapted for other bacteria such as Staphylococcus aureus and Streptococcus pneumoniae.

Distribution of mechanical stress in the Escherichia coli cell envelope

The cell envelope in Gram-negative bacteria comprises two distinct membranes with a cell wall between them. There has been a growing interest in understanding the mechanical adaptation of this cell envelope to the osmotic pressure (or turgor pressure), which is generated by the difference in the concentration of solutes between the cytoplasm and the external environment. However, it remains unexplored how the cell wall, the inner membrane (IM), and the outer membrane (OM) effectively protect the cell from this pressure by bearing the resulting surface tension, thus preventing the formation of inner membrane bulges, abnormal cell morphology, spheroplasts and cell lysis. In this study, we have used molecular dynamics (MD) simulations combined with experiments to resolve how and to what extent models of the IM, OM, and cell wall respond to changes in surface tension. We calculated the area compressibility modulus of all three components in simulations from tension-area isotherms. Experiments on monolayers mimicking individual leaflets of the IM and OM were also used to characterize their compressibility. While the membranes become softer as they expand, the cell wall exhibits significant strain stiffening at moderate to high tensions. We integrate these results into a model of the cell envelope in which the OM and cell wall share the tension at low turgor pressure (0.3 atm) but the tension in the cell wall dominates at high values (>1 atm).

Accurate In Silico Modeling of Asymmetric Bilayers Based on Biophysical Principles

Technological advances in the last decade have enabled the study of ever more complex and physiologically relevant model membranes to help dispel the mystery surrounding the role of plasma membrane asymmetry in various cellular processes. The slowly accumulating body of experimental data is fueling renewed interest in and the need for computational methods to support interpretations and address a wide range of problems that are still not amenable to direct experimental study. The specific appeal of molecular dynamics simulations for this purpose lies in their ability to access information at atomic resolution, which is useful for the formulation of testable mechanistic hypotheses. But, the range of questions that can be addressed reliably with such simulations is determined by the appropriate construction and simulation of asymmetric bilayer models. One essential way to achieve this goal is to follow rigorous biophysical criteria and principles. In this context, we show that the requirement for a robust comparison between the properties of simulated asymmetric and symmetric model membranes is for the tension in each bilayer leaflet to be zero. Commonly used methods for constructing asymmetric bilayers, including matching the average areas of the leaflets from the corresponding symmetric systems, do not ensure zero leaflet tension, thus precluding physically realistic changes in the areas of the two leaflets. We present, to our knowledge, a new method for identifying the ideal lipid packing in bilayers with different leaflet compositions that achieves the zero-tension goal, and discuss the basic principles underlying the biophysically correct computational study of asymmetric membranes.

Fate of Liposomes in the Presence of Phospholipase C and D: From Atomic to Supramolecular Lipid Arrangement

Understanding the origins of lipid membrane bilayer rearrangement in response to external stimuli is an essential component of cell biology and the bottom-up design of liposomes for biomedical applications. The enzymes phospholipase C and D (PLC and PLD) both cleave the phosphorus-oxygen bonds of phosphate esters in phosphatidylcholine (PC) lipids. The atomic position of this hydrolysis reaction has huge implications for the stability of PC-containing self-assembled structures, such as the cell wall and lipid-based vesicle drug delivery vectors. While PLC converts PC to diacylglycerol (DAG), the interaction of PC with PLD produces phosphatidic acid (PA). Here we present a combination of small-angle scattering data and all-atom molecular dynamics simulations, providing insights into the effects of atomic-scale reorganization on the supramolecular assembly of PC membrane bilayers upon enzyme-mediated incorporation of DAG or PA. We observed that PC liposomes completely disintegrate in the presence of PLC, as conversion of PC to DAG progresses. At lower concentrations, DAG molecules within fluid PC bilayers form hydrogen bonds with backbone carbonyl oxygens in neighboring PC molecules and burrow into the hydrophobic region. This leads initially to membrane thinning followed by a swelling of the lamellar phase with increased DAG. At higher DAG concentrations, localized membrane tension causes a change in lipid phase from lamellar to the hexagonal and micellar cubic phases. Molecular dynamics simulations show that this destabilization is also caused in part by the decreased ability of DAG-containing PC membranes to coordinate sodium ions. Conversely, PLD-treated PC liposomes remain stable up to extremely high conversions to PA. Here, the negatively charged PA headgroup attracts significant amounts of sodium ions from the bulk solution to the membrane surface, leading to a swelling of the coordinated water layer. These findings are a vital step toward a fundamental understanding of the degradation behavior of PC lipid membranes in the presence of these clinically relevant enzymes, and toward the rational design of diagnostic and drug delivery technologies for phospholipase-dysregulation-based diseases.

Droplet microfluidics for the construction of compartmentalised model membranes

The design of membrane-based constructs with multiple compartments is of increasing importance given their potential applications as microreactors, as artificial cells in synthetic-biology, as simplified cell models, and as drug delivery vehicles. The emergence of droplet microfluidics as a tool for their construction has allowed rapid scale-up in generation throughput, scale-down of size, and control over gross membrane architecture. This is true on several levels: size, level of compartmentalisation and connectivity of compartments can all be programmed to various degrees. This tutorial review explains and explores the reasons behind this. We discuss microfluidic strategies for the generation of a family of compartmentalised systems that have lipid membranes as the basic structural motifs, where droplets are either the fundamental building blocks, or are precursors to the membrane-bound compartments. We examine the key properties associated with these systems (including stability, yield, encapsulation efficiency), discuss relevant device fabrication technologies, and outline the technical challenges. In doing so, we critically review the state-of-play in this rapidly advancing field.

High-Speed Microscopy of Diffusion in Pore-Spanning Lipid Membranes

Pore-spanning membranes (PSMs) provide a highly attractive model system for investigating fundamental processes in lipid bilayers. We measure and compare lipid diffusion in the supported and suspended regions of PSMs prepared on a microfabricated porous substrate. Although some properties of the suspended regions in PSMs have been characterized using fluorescence studies, it has not been possible to examine the mobility of membrane components on the supported membrane parts. Here, we resolve this issue by employing interferometric scattering microscopy (iSCAT). We study the location-dependent diffusion of DOPE 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine) lipids (DOPE) labeled with gold nanoparticles in (1,2-dioleoyl-sn-glycero-3-phosphocholine) (DOPC) bilayers prepared on holey silicon nitride substrates that were either (i) oxygen-plasma-treated or (ii) functionalized with gold and 6-mercapto-1-hexanol. For both substrate treatments, diffusion in regions suspended on pores with diameters of 5 μm is found to be free. In the case of functionalization with gold and 6-mercapto-1-hexanol, similar diffusion coefficients are obtained for both the suspended and the supported regions, whereas for oxygen-plasma-treated surfaces, diffusion is almost 4 times slower in the supported parts of the membranes. We attribute this reduced diffusion on the supported parts in the case of oxygen-plasma-treated surfaces to larger membrane-substrate interactions, which lead to a higher membrane tension in the freestanding membrane parts. Furthermore, we find clear indications for a decrease of the diffusion constant in the freestanding regions away from the pore center. We provide a detailed characterization of the diffusion behavior in these membrane systems and discuss future directions.

Single-lipid tracking on nanoscale membrane buds: The effects of curvature on lipid diffusion and sorting

Nanoscale membrane curvature in cells is critical for endocytosis/exocytosis and membrane trafficking. However, the biophysical ramifications of nanoscale membrane curvature on the behavior of lipids remain poorly understood. Here, we created an experimental model system of membrane curvature at a physiologically-relevant scale and obtained nanoscopic information on single-lipid distributions and dynamics. Supported lipid bilayers were created over 50 and 70 nm radius nanoparticles to create membrane buds. Single-molecule localization microscopy was performed with diverse mixtures of fluorescent and non-fluorescent lipids. Variations in lipid acyl tales length, saturation, head-group, and fluorescent labeling strategy were tested while maintaining a single fluid lipid phase throughout the membrane. Monte Carlo simulations were used to fit our experimental results and quantify the effects of curvature on the lipid diffusion and sorting. Whereas varying the composition of the non-fluorescent lipids yielded minimal changes to the curvature effects, the labeling strategy of the fluorescent lipids yielded highly varying effects of curvature. Most conditions yield single-population Brownian diffusion throughout the membrane; however, curvature-induced lipid sorting, slowing, and aggregation were observed in some conditions. Head-group labeled lipids such as DPPE-Texas Red and POPE-Rhodamine diffused >2.4× slower on the curved vs. the planar membranes; tail-labeled lipids such as NBD-PPC, TopFluor-PPC, and TopFluor-PIP2, as well as DiIC₁₂ and DiIC₁₈ displayed no significant changes in diffusion due to the membrane curvature. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Lipid composition dictates serum stability of reconstituted high-density lipoproteins: implications for in vivo applications

Nanolipoprotein particles (NLPs) are reconstituted high-density lipoproteins, consisting of a phospholipid bilayer stabilized by an apolipoprotein scaffold protein. This class of nanoparticle has been a vital tool in the study of membrane proteins, and in recent years has been increasingly used for in vivo applications. Previous work demonstrated that the composition of the lipid bilayer component affects the stability of these particles in serum solutions. In the current study, NLPs assembled with phosphatidylcholine lipids featuring different acyl chain structures were systematically tested to understand the effect that lipid composition has on NLP stability in both neat serum and cell culture media supplemented with 10% serum by volume. The time at which 50% of the particles dissociate, as well as the fraction of the initial population that remains resistant to dissociation, were correlated to key parameters obtained from all-atom simulations of the corresponding lipid bilayers. A significant correlation was observed between the compressibility modulus of the lipid bilayer and particle stability in these complex biological milieu. These results can be used as a reference to tune the stability of these versatile biological nanoparticles for in vitro and in vivo applications.

Effect of 5-trans Isomer of Arachidonic Acid on Model Liposomal Membranes Studied by a Combined Simulation and Experimental Approach

Unsaturated fatty acids are found in humans predominantly in the cis configuration. Fatty acids in the trans configuration are primarily the result of human processing (trans fats), but can also be formed endogenously by radical stress. The cis-trans isomerization of fatty acids by free radicals could be connected to several pathologies. Trans fats have been linked to an increased risk of coronary artery disease; however, the reasons for the resulting pathogenesis remain unclear. Here, we investigate the effect of a mono-trans isomer of arachidonic acid (C20:4-5trans, 8cis, 11cis, 14cis) produced by free radicals in physiological concentration on a model erythrocyte membrane using a combined experimental and theoretical approach. Molecular Dynamics (MD) simulations of two model lipid bilayers containing arachidonic acid and its 5-trans isomer in 3 mol% were carried out for this purpose. The 5-trans isomer formation in the phospholipids was catalyzed by HOCH₂CH₂S· radicals, generated from the corresponding thiol by γ-irradiation, in multilamellar vesicles of SAPC. Large unilamellar vesicles were made by the extrusion method (LUVET) as a biomimetic model for cis-trans isomerization. Atomic Force Microscopy and Dynamic Light Scattering were used to measure the average size, morphology, and the z-potential of the liposomes. Both results from MD simulations and experiments are in agreement and indicate that the two model membranes display different physicochemical properties in that the bilayers containing the trans fatty acids were more ordered and more rigid than those containing solely the cis arachidonic acid. Correspondingly, the average size of the liposomes containing trans isomers was smaller than the ones without.

Enhancing the Stability and Immunomodulatory Activity of Liposomal Spherical Nucleic Acids through Lipid-Tail DNA Modifications

Liposomal spherical nucleic acids (LSNAs) are an attractive therapeutic platform for gene regulation and immunomodulation due to their biocompatibility, chemically tunable structures, and ability to enter cells rapidly without the need for ancillary transfection agents. Such structures consist of small (<100 nm) liposomal cores functionalized with a dense, highly oriented nucleic acid shell, both of which are key components in facilitating their biological activity. Here, the properties of LSNAs synthesized using conventional methods, anchoring cholesterol terminated oligonucleotides into a liposomal core, are compared to LSNAs made by directly modifying the surface of a liposomal core containing azide-functionalized lipids with dibenzocyclooctyl-terminated oligonucleotides. The surface densities of the oligonucleotides are measured for both types of LSNAs, with the lipid-modified structures having approximately twice the oligonucleotide surface coverage. The stabilities and cellular uptake properties of these structures are also evaluated. The higher density, lipid-functionalized structures are markedly more stable than conventional cholesterol-based structures in the presence of other unmodified liposomes and serum proteins as evidenced by fluorescence assays. Significantly, this new form of LSNA exhibits more rapid cellular uptake and increased sequence-specific toll-like receptor activation in immune reporter cell lines, making it a promising candidate for immunotherapy.

Tension Independence of Lipid Diffusion and Membrane Viscosity

The diffusion of biomolecules at lipid membranes is governed by the viscosity of the underlying two-dimensionally fluid lipid bilayer. For common three-dimensional fluids, viscosity can be modulated by hydrostatic pressure, and pressure-viscosity data have been measured for decades. Remarkably, the two-dimensional analogue of this relationship, the dependence of molecular mobility on tension, has to the best of our knowledge never been measured for lipid bilayers, limiting our understanding of cellular mechanotransduction as well as the fundamental fluid mechanics of membranes. Here we report both molecular-scale and mesoscopic measures of fluidity in giant lipid vesicles as a function of mechanical tension applied using micropipette aspiration. Both molecular-scale data, from fluorescence recovery after photobleaching, and micron-scale data, from tracking the diffusion of phase-separated domains, show a surprisingly weak dependence of viscosity on tension, in contrast to predictions of recent molecular dynamics simulations, highlighting fundamental gaps in our understanding of membrane fluidity.

Experimental Estimation of Membrane Tension Induced by Osmotic Pressure

Osmotic pressure (Π) induces the stretching of plasma membranes of cells or lipid membranes of vesicles, which plays various roles in physiological functions. However, there have been no experimental estimations of the membrane tension of vesicles upon exposure to Π. In this report, we estimated experimentally the lateral tension of the membranes of giant unilamellar vesicles (GUVs) when they were transferred into a hypotonic solution. First, we investigated the effect of Π on the rate constant, k_p, of constant-tension (σ_ex)-induced rupture of dioleoylphosphatidylcholine (DOPC)-GUVs using the method developed by us recently. We obtained the σ_ex dependence of k_p in GUVs under Π and by comparing this result with that in the absence of Π, we estimated the tension of the membrane due to Π at the swelling equilibrium, σ_osm^eq. Next, we measured the volume change of DOPC-GUVs under small Π. The experimentally obtained values of σ_osm^eq and the volume change agreed with their theoretical values within the limits of the experimental errors. Finally, we investigated the characteristics of the Π-induced pore formation in GUVs. The σ_osm^eq corresponding to the threshold Π at which pore formation is induced is similar to the threshold tension of the σ_ex-induced rupture. The time course of the radius change of GUVs in the Π-induced pore formation depends on the total membrane tension, σ_t; for small σ_t, the radius increased with time to an equilibrium one, which remained constant for a long time until pore formation, but for large σ_t, the radius increased with time and pore formation occurred before the swelling equilibrium was reached. Based on these results, we discussed the σ_osm^eq and the Π-induced pore formation in lipid membranes.

Biomimetic Phospholipid Membrane Organization on Graphene and Graphene Oxide Surfaces: A Molecular Dynamics Simulation Study

Supported phospholipid membrane patches stabilized on graphene surfaces have shown potential in sensor device functionalization, including biosensors and biocatalysis. Lipid dip-pen nanolithography (L-DPN) is a method useful in generating supported membrane structures that maintain lipid functionality, such as exhibiting specific interactions with protein molecules. Here, we have integrated L-DPN, atomic force microscopy, and coarse-grained molecular dynamics simulation methods to characterize the molecular properties of supported lipid membranes (SLMs) on graphene and graphene oxide supports. We observed substantial differences in the topologies of the stabilized lipid structures depending on the nature of the surface (polar graphene oxide vs nonpolar graphene). Furthermore, the addition of water to SLM systems resulted in large-scale reorganization of the lipid structures, with measurable effects on lipid lateral mobility within the supported membranes. We also observed reduced lipid ordering within the supported structures relative to free-standing lipid bilayers, attributed to the strong hydrophobic interactions between the lipids and support. Together, our results provide insight into the molecular effects of graphene and graphene oxide surfaces on lipid bilayer membranes. This will be important in the design of these surfaces for applications such as biosensor devices.

Whole-GUV patch-clamping

Studying how the membrane modulates ion channel and transporter activity is challenging because cells actively regulate membrane properties, whereas existing in vitro systems have limitations, such as residual solvent and unphysiologically high membrane tension. Cell-sized giant unilamellar vesicles (GUVs) would be ideal for in vitro electrophysiology, but efforts to measure the membrane current of intact GUVs have been unsuccessful. In this work, two challenges for obtaining the "whole-GUV" patch-clamp configuration were identified and resolved. First, unless the patch pipette and GUV pressures are precisely matched in the GUV-attached configuration, breaking the patch membrane also ruptures the GUV. Second, GUVs shrink irreversibly because the membrane/glass adhesion creating the high-resistance seal (>1 GΩ) continuously pulls membrane into the pipette. In contrast, for cell-derived giant plasma membrane vesicles (GPMVs), breaking the patch membrane allows the GPMV contents to passivate the pipette surface, thereby dynamically blocking membrane spreading in the whole-GMPV mode. To mimic this dynamic passivation mechanism, beta-casein was encapsulated into GUVs, yielding a stable, high-resistance, whole-GUV configuration for a range of membrane compositions. Specific membrane capacitance measurements confirmed that the membranes were truly solvent-free and that membrane tension could be controlled over a physiological range. Finally, the potential for ion transport studies was tested using the model ion channel, gramicidin, and voltage-clamp fluorometry measurements were performed with a voltage-dependent fluorophore/quencher pair. Whole-GUV patch-clamping allows ion transport and other voltage-dependent processes to be studied while controlling membrane composition, tension, and shape.

New insights into shear-sensitivity in dinoflagellate microalgae

A modification of a flow contraction device was used to subject shear-sensitive microalgae to well-defined hydrodynamic forces. The aim of the study was to elucidate if the inhibition of shear-induced growth commonly observed in dinoflagellate microalgae is in effect due to cell fragility that results in cell breakage even at low levels of turbulence. The microalgae assayed did not show any cell breakage even at energy dissipation rates (EDR) around 10(12)Wm(-3), implausible in culture devices. Conversely, animal cells, tested for comparison purposes, showed high physical cell damage at average EDR levels of 10(7)Wm(-3). Besides, very short exposures to high levels of EDR promoted variations in the membrane fluidity of the microalgae assayed, which might trigger mechanosensory cellular mechanisms. Average EDR values of only about 4·10(5)Wm(-3) increased cell membrane fluidity in microalgae whereas, in animal cells, they did not.

Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

LipidBuilder: A Framework To Build Realistic Models for Biological Membranes

The physical and chemical characterization of biological membranes is of fundamental importance for understanding the functional role of lipid bilayers in shaping cells and organelles, steering vesicle trafficking and promoting membrane-protein signaling. Molecular dynamics simulations stand as a powerful tool to probe the properties of membranes at atomistic level. However, the biological membrane is highly complex, and closely mimicking its physiological constitution in silico is not a straightforward task. Here, we present LipidBuilder, a framework for creating and storing models of biologically relevant phospholipid species with acyl tails of heterogeneous composition. LipidBuilder also enables the assembly of these database-stored lipids into realistic bilayers featuring asymmetric distribution on layer leaflets and concentration of given membrane constituents as defined, for example, by lipidomics experiments. The ability of LipidBuilder to assemble robust membrane models was validated by simulating membranes of homogeneous lipid composition for which experimental data are available. Furthermore, taking advantage of the extensive lipid headgroup repertoire, we assembled models of membranes of heterogeneous nature as naturally found in viral (phage PRD1), bacterial (Salmonella enterica, Laurinavicius , S. ; Kakela , R. ; Somerharju , P. ; Bamford , D. H. ; Virology 2004 , 322 , 328 - 336 ) and plant (Chlorella kessleri, Rezanka , T. ; Podojil , M. ; J. Chromatogr. 1989 , 463 , 397 - 408 ) organisms. These realistic membrane models were built using a near-exact lipid composition revealed from analytical chemistry experiments. We suggest LipidBuilder as a useful tool to model biological membranes of near-biological complexity, and as a robust complement to the current efforts to characterize the biophysical properties of biological membrane using molecular simulation.

Synthesis and proton NMR spectroscopy of intra-vesicular gamma-aminobutyric acid (GABA)

We report the synthesis of vesicles containing gamma-aminobutyric acid (GABA), and their proton nuclear magnetic resonance ((1)H NMR) spectra. These vesicles were constructed to more closely mimic the intracellular environment wherein GABA exists. For this study, these GABA-containing vesicles were examined under (1)H NMR as a potential platform for future studies on the differences between aqueous phantoms, ex vivo brain extracts, and in vivo magnetic resonance spectroscopy results. We found that intra-vesicular GABA faithfully yielded the chemical shifts and J-coupling constants of free aqueous GABA, alongside the chemical shift signals of the vesicle wall.

Elastic deformation and area per lipid of membranes: atomistic view from solid-state deuterium NMR spectroscopy

This article reviews the application of solid-state ²H nuclear magnetic resonance (NMR) spectroscopy for investigating the deformation of lipid bilayers at the atomistic level. For liquid-crystalline membranes, the average structure is manifested by the segmental order parameters (SCD) of the lipids. Solid-state ²H NMR yields observables directly related to the stress field of the lipid bilayer. The extent to which lipid bilayers are deformed by osmotic pressure is integral to how lipid-protein interactions affect membrane functions. Calculations of the average area per lipid and related structural properties are pertinent to bilayer remodeling and molecular dynamics (MD) simulations of membranes. To establish structural quantities, such as area per lipid and volumetric bilayer thickness, a mean-torque analysis of ²H NMR order parameters is applied. Osmotic stress is introduced by adding polymer solutions or by gravimetric dehydration, which are thermodynamically equivalent. Solid-state NMR studies of lipids under osmotic stress probe membrane interactions involving collective bilayer undulations, order-director fluctuations, and lipid molecular protrusions. Removal of water yields a reduction of the mean area per lipid, with a corresponding increase in volumetric bilayer thickness, by up to 20% in the liquid-crystalline state. Hydrophobic mismatch can shift protein states involving mechanosensation, transport, and molecular recognition by G-protein-coupled receptors. Measurements of the order parameters versus osmotic pressure yield the elastic area compressibility modulus and the corresponding bilayer thickness at an atomistic level. Solid-state ²H NMR thus reveals how membrane deformation can affect protein conformational changes within the stress field of the lipid bilayer.

Membrane Remodeling by Surface-Bound Protein Aggregates: Insights from Coarse-Grained Molecular Dynamics Simulation

The mechanism of curvature generation in membranes has been studied for decades due to its important role in many cellular functions. However, it is not clear if, or how, aggregates of lipid-anchored proteins might affect the geometry and elastic property of membranes. As an initial step toward addressing this issue, we performed structural, geometrical, and stress field analyses of coarse-grained molecular dynamics trajectories of a domain-forming bilayer in which an aggregate of lipidated proteins was asymmetrically bound. The results suggest a general mechanism whereby asymmetric incorporation of lipid-modified protein aggregates curve multidomain membranes primarily by expanding the surface area of the monolayer in which the lipid anchor is inserted.

Surface-Tension Replica-Exchange Molecular Dynamics Method for Enhanced Sampling of Biological Membrane Systems

Conformational sampling is fundamentally important for simulating complex biomolecular systems. The generalized-ensemble algorithm, especially the temperature replica-exchange molecular dynamics method (T-REMD), is one of the most powerful methods to explore structures of biomolecules such as proteins, nucleic acids, carbohydrates, and also of lipid membranes. T-REMD simulations have focused on soluble proteins rather than membrane proteins or lipid bilayers, because explicit membranes do not keep their structural integrity at high temperature. Here, we propose a new generalized-ensemble algorithm for membrane systems, which we call the surface-tension REMD method. Each replica is simulated in the NPγT ensemble, and surface tensions in a pair of replicas are exchanged at certain intervals to enhance conformational sampling of the target membrane system. We test the method on two biological membrane systems: a fully hydrated DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine) lipid bilayer and a WALP23-POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane system. During these simulations, a random walk in surface tension space is realized. Large-scale lateral deformation (shrinking and stretching) of the membranes takes place in all of the replicas without collapse of the lipid bilayer structure. There is accelerated lateral diffusion of DPPC lipid molecules compared with conventional MD simulation, and a much wider range of tilt angle of the WALP23 peptide is sampled due to large deformation of the POPC lipid bilayer and through peptide-lipid interactions. Our method could be applicable to a wide variety of biological membrane systems.

The stretching elasticity of biomembranes determines their line tension and bending rigidity

In this work, some implications of a recent model for the mechanical behavior of biological membranes (Deseri et al. in Continuum Mech Thermodyn 20(5):255-273, 2008) are exploited by means of a prototypical one-dimensional problem. We show that the knowledge of the membrane stretching elasticity permits to establish a precise connection among surface tension, bending rigidities and line tension during phase transition phenomena. For a specific choice of the stretching energy density, we evaluate these quantities in a membrane with coexistent fluid phases, showing a satisfactory comparison with the available experimental measurements. Finally, we determine the thickness profile inside the boundary layer where the order-disorder transition is observed.