Simulations of simple linoleic acid-containing lipid membranes and models for the soybean plasma membranes

Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer

Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Ų. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.

Structural and functional evaluation mammalian and plant lipoxygenases upon association with nanodics as membrane mimetics

Lipoxygenases (LOX) are a family lipid oxygenating enzymes that can generate bioactive lipids of clinical relevance from polyunsaturated fatty acids. Most LOXs display a Ca²⁺-dependent association with membranes for their activity. Nanodiscs (ND) are stable self-assembled discoidal fragments of lipid bilayers that can mimic the plasma membrane. In this study, we evaluated the association of mammalian 15-LOXs (ALOX15 and ALOX15B) and soybean LOX-1 with NDs (LOX-ND), their enzymatic activities and inhibition. Mammalian LOXs associated with NDs showed better retention of enzymatic function compared to soybean LOX-1. Treatment of both LOX-NDs and free enzymes with the pan-LOX inhibitor nordihydroguaiaretic acid (NDGA) showed an approximately 5-fold more effective inhibition of the enzymes associated with NDs compared to the free form. NDs are easy to generate membrane mimics that can be used as an effective tool to determine enzymatic function and inhibition of membrane associated proteins.

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

All-Atom Modeling of Complex Cellular Membranes

Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.

Symmetric and Asymmetric Models for the Arabidopsis thaliana Plasma Membrane: A Simulation Study

Arabidopsis thaliana is an important model organism, which has attracted many biologists. While most research efforts have been on studying the genetics and proteins of this organism, a systematic study of its lipidomics is lacking. Here, we present a novel, asymmetric model of its cell membrane with its lipid composition consisting of five glycerophospholipids, two glycolipids, and sitosterol determined from multiple independent experiments. A typical lipid type in plant membranes is glycosyl inositol phosphoryl ceramide (GIPC), which accounts for about 10% of the total lipids in the outer leaflet in our model. Two symmetric models representing the inner and outer leaflets of the membrane were built and simulated until equilibrium was reached and then combined to form the asymmetric model. Our results indicate that the outer leaflet is more rigid and tightly packed compared to the inner leaflet. Pressure profiles for the two leaflets are overall similar though the outer leaflet exhibits larger oscillations. A special focus on lipid organization is discussed and the interplay between glycolipids and sitosterols is found to be important. The current model provides a baseline for future modeling of similar membranes and can be used to study partitioning of small molecules in the membrane or further developed to study the interaction between plant membrane proteins and lipids.

The structure and flexibility analysis of the Arabidopsis synaptotagmin 1 reveal the basis of its regulation at membrane contact sites

Non-vesicular lipid transfer at ER and plasma membrane (PM) contact sites (CS) is crucial for the maintenance of membrane lipid homeostasis. Extended synaptotagmins (E-Syts) play a central role in this process as they act as molecular tethers of ER and PM and as lipid transfer proteins between these organelles. E-Syts are proteins constitutively anchored to the ER through an N-terminal hydrophobic segment and bind the PM via a variable number of C-terminal C2 domains. Synaptotagmins (SYTs) are the plant orthologous of E-Syts and regulate the ER-PM communication in response to abiotic stress. Combining different structural and biochemical techniques, we demonstrate that the binding of SYT1 to lipids occurs through a Ca2+-dependent lipid-binding site and by a site for phosphorylated forms of phosphatidylinositol, thus integrating two different molecular signals in response to stress. In addition, we show that SYT1 displays three highly flexible hinge points that provide conformational freedom to facilitate lipid extraction, protein loading, and subsequent transfer between PM and ER.

Plant Sterol Clustering Correlates with Membrane Microdomains as Revealed by Optical and Computational Microscopy

Local inhomogeneities in lipid composition play a crucial role in the regulation of signal transduction and membrane traffic. This is particularly the case for plant plasma membrane, which is enriched in specific lipids, such as free and conjugated forms of phytosterols and typical phytosphingolipids. Nevertheless, most evidence for microdomains in cells remains indirect, and the nature of membrane inhomogeneities has been difficult to characterize. We used a new push–pull pyrene probe and fluorescence lifetime imaging microscopy (FLIM) combined with all-atom multiscale molecular dynamics simulations to provide a detailed view on the interaction between phospholipids and phytosterol and the effect of modulating cellular phytosterols on membrane-associated microdomains and phase separation formation. Our understanding of the organization principles of biomembranes is limited mainly by the challenge to measure distributions and interactions of lipids and proteins within the complex environment of living cells. Comparing phospholipids/phytosterol compositions typical of liquid-disordered (Ld) and liquid-ordered (Lo) domains, we furthermore show that phytosterols play crucial roles in membrane homeostasis. The simulation work highlights how state-of-the-art modeling alleviates some of the prior concerns and how unrefuted discoveries can be made through a computational microscope. Altogether, our results support the role of phytosterols in the lateral structuring of the PM of plant cells and suggest that they are key compounds for the formation of plant PM microdomains and the lipid-ordered phase.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

Simulations of Diabetic and Non-Diabetic Peripheral Nerve Myelin Lipid Bilayers

The multilayered myelin sheath is a critical component of both central and peripheral nervous systems, forming a protective barrier against axonal damage and facilitating the movement of nervous impulses. It is primarily composed of cholesterol (CHL1), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), phosphatidylinositol (PI), sphingomyelin (SM), and galactosylceramide (GalCer) lipids. For rat sciatic nerve myelin (part of the peripheral nervous system, PNS), it has been found that cholesterol and unsaturated fatty acid contents are significantly lower in diabetic than in non-diabetic conditions. In this study, lipid compositions from experimental data are used to create four model rat sciatic nerve myelin lipid bilayers: PI-containing (non-diabetic and diabetic) and PS-containing (non-diabetic and diabetic), which were then simulated using the all-atom CHARMM36 force field. Simulation results of diabetic membranes indicate less rigid, more laterally expansive, and thinner bilayers as well as potentially reduced interactions between GalCer on opposing myelin leaflets, supporting a direct role of the cholesterol content decrease in instigating myelin deterioration and diabetic peripheral neuropathy. Compared to PI-lipids, PS-lipids were found to cause higher inter-lipid spacing and decreased order within membranes as a result of their smaller headgroup size and higher inter-lipid hydrogen bonding potential, which allow them to more frequently reside deeper in the membrane plane and produce pushing effects on other lipids. GalCer deuterium order parameters and non-diabetic headgroup-to-headgroup bilayer thicknesses were compared to experimental data, exhibiting close alignment, supporting the future usage of these models to study the PNS myelin sheath.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Characteristics of a Folate Receptor-α Anchored into a Multilipid Bilayer Obtained from Atomistic Molecular Dynamics Simulations

Thorough computational description of the properties of membrane-anchored protein receptors, which are important for example in the context of active targeting drug delivery, may be achieved by models representing as close as possible the immediate environment of these macromolecules. An all-atom bilayer, including 35 different lipid types asymmetrically distributed among the two monolayers, is suggested as a model neoplastic cell membrane. One molecule of folate receptor-α (FRα) is anchored into its outer leaflet, and the behavior of the system is explored by atomistic molecular dynamics simulations. The total number of atoms in the model is ∼185 000. Three 1-μs-long simulations are carried out, where physiological conditions (310 K and 1 bar) are maintained with three different pressure scaling schemes. To evaluate the structure and the phase state of the membrane, the density profiles of the system, the average area per lipid, and the deuterium order parameter of the lipid tails are calculated. The bilayer is in liquid ordered state, and the specific arrangement varies between the three trajectories. The changes in the structure of FRα are investigated and are found time- and ensemble-dependent. The volume of the ligand binding pocket fluctuates with time, but this variation remains independent of the more global structural alterations. The latter are mostly "waving" motions of the protein, which periodically approaches and retreats from the membrane. The semi-isotropic pressure scaling perturbs the receptor most significantly, while the isotropic algorithm induces rather slow changes. Maintaining constant nonzero surface tension leads to behavior closest to the experimentally observed one.

Styrene maleic acid copolymer induces pores in biomembranes

We investigated the interactions between styrene-maleic acid (SMA) copolymers and phospholipid bilayers, using confocal microscopy and surface acoustic wave resonance (SAR) sensing. For the first time we experimentally observed and followed pore formation by SMA copolymers in intact supported bilayers and unilamellar vesicles, showing that fluorescein, a water-soluble organic compound with a mean diameter of 6.9 Å, can traverse the membrane. Our findings are in agreement with recent theoretical predictions, which suggested that SMA copolymers may insert along the plane of the bilayer to form stable toroidal pores.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Plant lipids: Key players of plasma membrane organization and function

The plasma membrane (PM) is the biological membrane that separates the interior of all cells from the outside. The PM is constituted of a huge diversity of proteins and lipids. In this review, we will update the diversity of molecular species of lipids found in plant PM. We will further discuss how lipids govern global properties of the plant PM, explaining that plant lipids are unevenly distributed and are able to organize PM in domains. From that observation, it emerges a complex picture showing a spatial and multiscale segregation of PM components. Finally, we will discuss how lipids are key players in the function of PM in plants, with a particular focus on plant-microbe interaction, transport and hormone signaling, abiotic stress responses, plasmodesmata function. The last chapter is dedicated to the methods that the plant membrane biology community needs to develop to get a comprehensive understanding of membrane organization in plants.

Divide and Rule: Plant Plasma Membrane Organization

Since the publication of the fluid mosaic as a relevant model for biological membranes, accumulating evidence has revealed the outstanding complexity of the composition and organization of the plant plasma membrane (PM). Powerful new methodologies have uncovered the remarkable multiscale and multicomponent heterogeneity of PM subcompartmentalization, and this is emerging as a general trait with different features and properties. It is now evident that the dynamics of such a complex organization are intrinsically related to signaling pathways that regulate key physiological processes. Listing and linking recent progress in precisely qualifying these heterogeneities will help to draw an integrated picture of the plant PM. Understanding the key principles governing such a complex dynamic organization will contribute to deciphering the crucial role of the PM in cell physiology.