Protein-induced membrane curvature in coarse-grained simulations

Delineating the shape of COat Protein complex-II coated membrane bud

Curvature-generating proteins that direct membrane trafficking assemble on the surface of lipid bilayers to bud transport intermediates, which move protein and lipid cargoes from one cellular compartment to another. However, it remains unclear what controls the overall shape of the membrane bud once curvature induction has begun. In vitro experiments showed that excessive concentrations of the COPII protein Sar1 promoted the formation of membrane tubules from synthetic vesicles, while COPII-coated transport intermediates in cells are generally more spherical or lobed in shape. To understand the origin of these morphological differences, we employ atomistic, coarse-grained (CG), and continuum mesoscopic simulations of membranes in the presence of multiple curvature-generating proteins. We first characterize the membrane-bending ability of amphipathic peptides derived from the amino terminus of Sar1, as a function of interpeptide angle and concentration using an atomistic bicelle simulation protocol. Then, we employ CG simulations to reveal that Sec23 and Sec24 control the relative spacing between Sar1 protomers and form the inner-coat unit through an attachment with Sar1. Finally, using dynamical triangulated surface simulations based on the Helfrich Hamiltonian, we demonstrate that the uniform distribution of spacer molecules among curvature-generating proteins is crucial to the spherical budding of the membrane. Overall, our analyses suggest a new role for Sec23, Sec24, and cargo proteins in COPII-mediated membrane budding process in which they act as spacers to preserve a dispersed arrangement of Sar1 protomers and help determine the overall shape of the membrane bud.

Making the cut: Multiscale simulation of membrane remodeling

Biological membranes are dynamic heterogeneous materials, and their shape and organization are tightly coupled to the properties of the proteins in and around them. However, the length scales of lipid and protein dynamics are far below the size of membrane-bound organelles, much less an entire cell. Therefore, multiscale modeling approaches are often necessary to build a comprehensive picture of the interplay of these factors, and have provided critical insights into our understanding of membrane dynamics. Here, we review computational methods for studying membrane remodeling, as well as passive and active examples of protein-driven membrane remodeling. As the field advances towards the modeling of key aspects of organelles and whole cells - an increasingly accessible regime of study - we summarize here recent successes and offer comments on future trends.

Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates

Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.

Mechanism of Membrane Curvature Induced by SNX1: Insights from Molecular Dynamics Simulations

SNX proteins have been found to induce membrane remodeling to facilitate the generation of transport carriers in endosomal pathways. However, the molecular mechanism of membrane bending and the role of lipids in the bending process remain elusive. Here, we conducted coarse-grained molecular dynamics simulations to investigate the role of the three structural modules (PX, BAR, and AH) of SNX1 and the PI3P lipids in membrane deformation. We observed that the presence of all three domains is essential for SNX1 to achieve a stable membrane deformation. BAR is capable of remodeling the membrane through the charged residues on its concave surface, but it requires PX and AH to establish stable membrane binding. AH penetrates into the lipid membrane, thereby promoting the induction of membrane curvature; however, it is inadequate on its own to maintain membrane bending. PI3P lipids are also indispensable for membrane remodeling, as they play a dominant role in the interactions of lipids with the BAR domain. Our results enhance the comprehension of the molecular mechanism underlying SNX1-induced membrane curvature and help future studies of curvature-inducing proteins.

Membranes get in shape: Biophysics of curving bilayers

Membrane curvature is ubiquitous and essential in cell biology. Curved membranes have several distinct features, including specific protein and lipid sorting, distinct lipid ordering, and changes in transbilayer stress. Curvature also interplays with membrane tension to generate forces that change membrane shape. This research highlight summarizes recent contributions to this topic published in Biophysical Journal.

Amphipathic Helices Can Sense Both Positive and Negative Curvatures of Lipid Membranes

Curvature sensing is an essential ability of biomolecules to preferentially localize to membrane regions of a specific curvature. It has been shown that amphipathic helices (AHs), helical peptides with both hydrophilic and hydrophobic regions, could sense a positive membrane curvature. The origin of this AH sensing has been attributed to their ability to exploit lipid-packing defects that are enhanced in regions of positive curvature. In this study, we revisit an alternative framework where AHs act as sensors of local internal stress within the membrane, suggesting the possibility of an AH sensing a negative membrane curvature. Using molecular dynamics simulations, we gradually tuned the hydrophobicity of AHs, thereby adjusting their insertion depth so that the curvature preference of AHs is switched from positive to negative. This study suggests that highly hydrophobic AHs could preferentially localize proteins to regions of a negative membrane curvature.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Simulation study of membrane bending by protein crowding: a case study with the epsin N-terminal homology domain

The mechanisms by which peripheral membrane proteins generate curvature is currently an active area of research. One of the proposed mechanisms is amphipathic insertion or the 'wedge' mechanism in which the protein shallowly inserts an amphipathic helix inside the membrane to drive the curvature. However, recent experimental studies have challenged the efficiency of the 'wedge' mechanism as it requires unusual protein densities. These studies proposed an alternative mechanism, namely 'protein-crowding', in which the lateral pressure generated by the random collisions among the membrane bound proteins drives the bending. In this study, we employ atomistic and coarse-grained molecular dynamics simulations to investigate the effects of amphipathic insertion and protein crowding on the membrane surface. Considering epsin N-terminal homology (ENTH) domain as a model protein, we show that amphipathic insertion is not essential for membrane bending. Our results suggest that ENTH domains can aggregate on the membrane surface by employing another structured region (H₃ helix). And this protein crowding decreases the cohesive energy of the lipid tails which causes a significant decrease in the membrane bending rigidity. The ENTH domain can generate a similar degree of membrane curvature irrespective of the activity of its H₀ helix. Our results are consistent with the recent experimental results.

Molecular mechanism of GTP binding- and dimerization-induced enhancement of Sar1-mediated membrane remodeling

The Sar1 GTPase initiates coat protein II (COPII)-mediated protein transport by generating membrane curvature at subdomains on the endoplasmic reticulum, where it is activated by the guanine nucleotide exchange factor (GEF) Sec12. Crystal structures of GDP- and GTP-bound forms of Sar1 suggest that it undergoes a conformational switch in which GTP binding enhances the exposure of an amino-terminal amphipathic helix necessary for efficient membrane penetration. However, key residues in the amino terminus were not resolved in crystal structures, and experimental studies have suggested that the amino terminus of Sar1 is solvent-exposed in the absence of a membrane, even in the GDP-bound state. Therefore, the molecular mechanism by which GTP binding activates the membrane-remodeling activity of Sar1 remains unclear. Using atomistic molecular dynamics simulations, we compare the membrane-binding and curvature generation activities of Sar1 in its GDP- and GTP-bound states. We show that in the GTP-bound state, Sar1 inserts into the membrane with its complete (residues 1 to 23) amphipathic amino-terminal helix, while Sar1-GDP binds to the membrane only through its first 12 residues. Such differential membrane-binding modes translate into significant differences in the protein volume inserted into the membrane. As a result, Sar1-GTP generates positive membrane curvature 10 to 20 times higher than Sar1-GDP. Dimerization of the GTP-bound form of Sar1 further amplifies curvature generation. Taken together, our results present a detailed molecular mechanism for how the nucleotide-bound state of Sar1 regulates its membrane-binding and remodeling activities in a concentration-dependent manner, paving the way toward a better understanding COPII-mediated membrane transport.

Magainin 2 and PGLa in bacterial membrane mimics IV: Membrane curvature and partitioning

We previously reported that the synergistically enhanced antimicrobial activity of magainin 2 (MG2a) and PGLa is related to membrane adhesion and fusion. Here, we demonstrate that equimolar mixtures of MG2a and L18W-PGLa induce positive monolayer curvature stress and sense, at the same time, positive mean and Gaussian bilayer curvatures already at low amounts of bound peptide. The combination of both abilities-membrane curvature sensing and inducing-is most likely the base for the synergistically enhanced peptide activity. In addition, our coarse-grained simulations suggest that fusion stalks are promoted by decreasing the free-energy barrier for their formation rather than by stabilizing their shape. We also interrogated peptide partitioning as a function of lipid and peptide concentration using tryptophan fluorescence spectroscopy and peptide-induced leakage of dyes from lipid vesicles. In agreement with a previous report, we find increased membrane partitioning of L18W-PGLa in the presence of MG2a. However, this effect does not prevail to lipid concentrations higher than 1 mM, above which all peptides associate with the lipid bilayers. This implies that synergistic effects of MG2a and L18W-PGLa in previously reported experiments with lipid concentrations >1 mM are due to peptide-induced membrane remodeling and not their specific membrane partitioning.

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

Simulation study of domain formation in a model bacterial membrane

Recent experimental studies revealed that functional membrane microdomains (FMMs) are formed in prokaryotic cells which are structurally and functionally similar to the lipid rafts formed in eukaryotic cells. In this study, we employ coarse-grained molecular dynamics simulations to investigate the mechanism of domain formation and its physiochemical properties in a model methicillin-resistant staphylococcus aureus (MRSA) cell membrane. We find that domains are formed through lateral segregation of staphyloxanthin (STX), a carotenoid which shields the bacteria from the host's immune because of its antioxidant nature. Simulation results suggest that membrane integrity increases with the size of the domain, which is assessed by computing bond order parameter of the lipid tails, membrane expansion modulus and water permeability across the membrane. Various membrane domain proteins such as flotillin-like protein floA and penicillin binding protein (PBP2a) preferentially bind with the STX and accumulate in the membrane domain which is consistent with the recent experimental results.

Coacervation of poly-electrolytes in the presence of lipid bilayers: mutual alteration of structure and morphology

Many intrinsically disordered peptides have been shown to undergo liquid-liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E30) and poly-lysine (K30) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting.

Cholesterol-Mediated Clustering of the HIV Fusion Protein gp41 in Lipid Bilayers

The envelope glycoprotein (Env) of the human immunodeficient virus (HIV-1) is known to cluster on the viral membrane surface to attach to target cells and cause membrane fusion for HIV-1 infection. However, the molecular structural mechanisms that drive Env clustering remain opaque. Here, we use solid-state NMR spectroscopy and molecular dynamics (MD) simulations to investigate nanometer-scale clustering of the membrane-proximal external region (MPER) and transmembrane domain (TMD) of gp41, the fusion protein component of Env. Using ¹⁹F solid-state NMR experiments of mixed fluorinated peptides, we show that MPER-TMD trimers form clusters with interdigitated MPER helices in cholesterol-containing membranes. Inter-trimer ¹⁹F-¹⁹F cross peaks, which are indicative of spatial contacts within ∼2 nm, are observed in cholesterol-rich virus-mimetic membranes but are suppressed in cholesterol-free model membranes. Water-peptide and lipid-peptide cross peaks in 2D ¹H-¹⁹F correlation spectra indicate that the MPER is well embedded in model phosphocholine membranes but is more exposed to the surface of the virus-mimetic membrane. These experimental results are reproduced in coarse-grained and atomistic molecular dynamics simulations, which suggest that the effects of cholesterol on gp41 clustering is likely via indirect modulation of the MPER orientation. Cholesterol binding to the helix-turn-helix region of the MPER-TMD causes a parallel orientation of the MPER with the membrane surface, thus allowing MPERs of neighboring trimers to interact with each other to cause clustering. These solid-state NMR data and molecular dynamics simulations suggest that MPER and cholesterol cooperatively govern the clustering of gp41 trimers during virus-cell membrane fusion.