Salt Bridge Formation between the I-BAR Domain and Lipids Increases Lipid Density and Membrane Curvature

Molecular Dynamics Reveal Key Steps in BAR-Related Membrane Remodeling

Endocytosis plays a complex role in pathogen-host interactions. It serves as a pathway for pathogens to enter the host cell and acts as a part of the immune defense mechanism. Endocytosis involves the formation of lipid membrane vesicles and the reshaping of the cell membrane, a task predominantly managed by proteins containing BAR (Bin1/Amphiphysin/yeast RVS167) domains. Insights into how BAR domains can remodel and reshape cell membranes provide crucial information on infections and can aid the development of treatment. Aiming at deciphering the roles of the BAR dimers in lipid membrane bending and remodeling, we conducted extensive all-atom molecular dynamics simulations and discovered that the presence of helix kinks divides the BAR monomer into two segments-the "arm segment" and the "core segment"-which exhibit distinct movement patterns. Contrary to the prior hypothesis of BAR domains working as a rigid scaffold, we found that it functions in an "Arms-Hands" mode. These findings enhance the understanding of endocytosis, potentially advancing research on pathogen-host interactions and aiding in the identification of new treatment strategies targeting BAR domains.

Binding of anisotropic curvature-inducing proteins onto membrane tubes

Bin/Amphiphysin/Rvs superfamily proteins and other curvature-inducing proteins have anisotropic shapes and anisotropically bend biomembranes. Here, we report how the anisotropic proteins bind the membrane tube and are orientationally ordered using mean-field theory including an orientation-dependent excluded volume. The proteins exhibit a second-order or first-order nematic transition with increasing protein density depending on the radius of the membrane tube. The tube curvatures for the maximum protein binding and orientational order are different and varied by the protein density and rigidity. As the external force along the tube axis increases, a first-order transition from a large tube radius with low protein density to a small radius with high density occurs once, and subsequently, the protein orientation tilts to the tube-axis direction. When an isotropic bending energy is used for the proteins with an elliptic shape, the force-dependence curves become symmetric and the first-order transition occurs twice. This theory quantitatively reproduces the results of meshless membrane simulation for short proteins, whereas deviations are seen for long proteins owing to the formation of protein clusters.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Mechanism of negative membrane curvature generation by I-BAR domains

The membrane sculpting ability of BAR domains has been attributed to the intrinsic curvature of their banana-shaped dimeric structure. However, there is often a mismatch between this intrinsic curvature and the diameter of the membrane tubules generated. I-BAR domains are especially mysterious since they are almost flat but generate high negative membrane curvature. Here, we use atomistic implicit-solvent computer modeling to show that the membrane bending of the IRSp53 I-BAR domain is dictated by its higher oligomeric structure, whose curvature is completely unrelated to the intrinsic curvature of the dimer. Two other I-BARs give similar results, whereas a flat F-BAR sheet develops a concave membrane-binding interface, consistent with its observed positive membrane curvature generation. Laterally interacting helical spirals of I-BAR dimers on tube interiors are stable and have an enhanced binding energy that is sufficient for membrane bending to experimentally observed tubule diameters at a reasonable surface density.

Full assembly of HIV-1 particles requires assistance of the membrane curvature factor IRSp53

During HIV-1 particle formation, the requisite plasma membrane curvature is thought to be solely driven by the retroviral Gag protein. Here, we reveal that the cellular I-BAR protein IRSp53 is required for the progression of HIV-1 membrane curvature to complete particle assembly. siRNA-mediated knockdown of IRSp53 gene expression induces a decrease in viral particle production and a viral bud arrest at half completion. Single-molecule localization microscopy at the cell plasma membrane shows a preferential localization of IRSp53 around HIV-1 Gag assembly sites. In addition, we observe the presence of IRSp53 in purified HIV-1 particles. Finally, HIV-1 Gag protein preferentially localizes to curved membranes induced by IRSp53 I-BAR domain on giant unilamellar vesicles. Overall, our data reveal a strong interplay between IRSp53 I-BAR and Gag at membranes during virus assembly. This highlights IRSp53 as a crucial host factor in HIV-1 membrane curvature and its requirement for full HIV-1 particle assembly.

Binding, unbinding and aggregation of crescent-shaped nanoparticles on nanoscale tubular membranes

Using molecular dynamics simulations of a coarse-grained implicit solvent model, we investigate the binding of crescent-shaped nanoparticles (NPs) on tubular lipid membranes. The NPs adhere to the membrane through their concave side. We found that the binding/unbinding transition is first-order, with the threshold binding energy being higher than the unbinding threshold, and the energy barrier between the bound and unbound states at the transition that increases with increasing the NP's arclength L_np or curvature mismatch μ = R_c/R_np, where R_c and R_np are the radii of curvature of the tubular membrane and the NP, respectively. Furthermore, we found that the threshold binding energy increases with increasing either L_np or μ. NPs with curvature larger than that of the tubule (μ > 1) lie perpendicularly to the tubule's axis. However, for μ smaller than a specific arclength-dependent mismatch μ*, the NPs are tilted with respect to the tubule's axis, with the tilt angle that increases with decreasing μ. We also investigated the self-assembly of the NPs on the tubule at relatively weak adhesion strength and found that for μ > 1 and high values of L_np, the NPs self-assemble into linear chains, and lie side-by-side. For μ < μ* and high L_np, the NPs also self-assemble into chains, while being tilted with respect to the tubule's axis.

Lipid-Composition-Mediated Forces Can Stabilize Tubular Assemblies of I-BAR Proteins

Collective action by inverse-Bin/Amphiphysin/Rvs (I-BAR) domains drive micron-scale membrane remodeling. The macroscopic curvature sensing and generation behavior of I-BAR domains is well characterized, and computational models have suggested various mechanisms on simplified membrane systems, but there remain missing connections between the complex environment of the cell and the models proposed thus far. Here, we show a connection between the role of protein curvature and lipid clustering in the relaxation of large membrane deformations. When we include phosphatidylinositol 4,5-bisphosphate-like lipids that preferentially interact with the charged ends of an I-BAR domain, we find clustering of phosphatidylinositol 4,5-bisphosphate-like lipids that induce a directional membrane-mediated interaction between membrane-bound I-BAR domains. Lipid clusters mediate I-BAR domain interactions and cause I-BAR domain aggregates that would not arise through membrane fluctuation-based or curvature-based interactions. Inside of membrane protrusions, lipid cluster-mediated interaction draws long side-by-side aggregates together, resulting in more cylindrical protrusions as opposed to bulbous, irregularly shaped protrusions.

Mechanisms of negative membrane curvature sensing and generation by ESCRT III subunit Snf7

Certain proteins have the propensity to bind to negatively curved membranes and generate negative membrane curvature. The mechanism of action of these proteins is much less studied and understood than those that sense and generate positive curvature. In this work, we use implicit membrane modeling to explore the mechanism of an important negative curvature sensing and generating protein: the main ESCRT III subunit Snf7. We find that Snf7 monomers alone can sense negative curvature and that curvature sensitivity increases for dimers and trimers. We have observed spontaneous bending of Snf7 oligomers into circular structures with preferred radius of ~20 nm. The preferred curvature of Snf7 filaments is further confirmed by the simulations of preformed spirals on a cylindrical membrane surface. Snf7 filaments cannot bind with the same interface to flat and curved membranes. We find that even when a filament has the preferred radius, it is always less stable on the flat membrane surface than on the interior cylindrical membrane surface. This provides an additional energy for membrane bending which has not been considered in the spiral spring model. Furthermore, the rings on the cylindrical spirals are bridged together by helix 4 and hence are extra stabilized compared to the spirals on the flat membrane surface.

Curvature induction and sensing of the F-BAR protein Pacsin1 on lipid membranes via molecular dynamics simulations

F-Bin/Amphiphysin/Rvs (F-BAR) domain proteins play essential roles in biological processes that involve membrane remodelling, such as endocytosis and exocytosis. It has been shown that such proteins transform the lipid membrane into tubes. Notably, Pacsin1 from the Pacsin/Syndapin subfamily has the ability to transform the membrane into various morphologies: striated tubes, featureless wide and thin tubes, and pearling vesicles. The molecular mechanism of this interesting ability remains elusive. In this study, we performed all-atom (AA) and coarse-grained (CG) molecular dynamics simulations to investigate the curvature induction and sensing mechanisms of Pacsin1 on a membrane. From AA simulations, we show that Pacsin1 has internal structural flexibility. In CG simulations with parameters tuned from the AA simulations, spontaneous assembly of two Pacsin1 dimers through lateral interaction is observed. Based on the complex structure, we show that the regularly assembled Pacsin1 dimers bend a tensionless membrane. We also show that a single Pacsin1 dimer senses the membrane curvature, binding to a buckled membrane with a preferred curvature. These results provide molecular insights into polymorphic membrane remodelling.

Shape transition from elliptical to cylindrical membrane tubes induced by chiral crescent-shaped protein rods

Proteins often form chiral assembly structures on a biomembrane. However, the role of the chirality in the interaction with an achiral membrane is poorly understood. Here, we report how chirality of crescent-shaped protein rods changes their assembly and tubulation using meshless membrane simulations. The achiral rods deformed the membrane tube into an elliptical shape by stabilizing the edges of the ellipse. In contrast, the chiral rods formed a helical assembly that generated a cylindrical membrane tube with a constant radius in addition to the elliptical tube. This helical assembly could be further stabilized by the direct side-to-side attraction between the protein rods. The chirality also promotes the tubulation from a flat membrane. These results agree with experimental findings of the constant radius of membrane tubules induced by the Bin/Amphiphysin/Rvs (BAR) superfamily proteins.

Unusual Organization of I-BAR Proteins on Tubular and Vesicular Membranes

Protein-mediated membrane remodeling is a ubiquitous and critical process for proper cellular function. Inverse Bin/Amphiphysin/Rvs (I-BAR) domains drive local membrane deformation as a precursor to large-scale membrane remodeling. We employ a multiscale approach to provide the molecular mechanism of unusual I-BAR domain-driven membrane remodeling at a low protein surface concentration with near-atomistic detail. We generate a bottom-up coarse-grained model that demonstrates similar membrane-bound I-BAR domain aggregation behavior as our recent Mesoscopic Membrane with Explicit Proteins model. Together, these models bridge several length scales and reveal an aggregation behavior of I-BAR domains. We find that at low surface coverage (i.e., low bound protein density), I-BAR domains form transient, tip-to-tip strings on periodic flat membrane sheets. Inside of lipid bilayer tubules, we find linear aggregates parallel to the axis of the tubule. Finally, we find that I-BAR domains form tip-to-tip aggregates around the edges of membrane domes. These results are supported by in vitro experiments showing low curvature bulges surrounded by I-BAR domains on giant unilamellar vesicles. Overall, our models reveal new I-BAR domain aggregation behavior in membrane tubules and on the surface of vesicles at low surface concentration that add insight into how I-BAR domain proteins may contribute to certain aspects of membrane remodeling in cells.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

Organizing membrane-curving proteins: the emerging dynamical picture

Lipid membranes play key roles in cells, such as in trafficking, division, infection, remodeling of organelles, among others. The key step in all these processes is creating membrane curvature, typically under the control of many anchored, adhered or included proteins. However, it has become clear that the membrane itself can mediate the interactions among proteins to produce highly ordered assemblies. Computer simulations are ideally suited to investigate protein organization and the dynamics of membrane remodeling at near-micron scales, something that is extremely challenging to tackle experimentally. We review recent computational efforts in modeling protein-caused membrane deformation mechanisms, specifically focusing on coarse-grained simulations. We highlight work that exposed the membrane-mediated ordering of proteins into lines, meshwork, spirals and other assemblies, in what seems to be a very generic mechanism driven by a combination of short and long-ranged forces. Modulating the mechanical properties of membranes is an underexplored signaling mechanism in various processes deserving of more attention in the near future.

Biophysics of membrane curvature remodeling at molecular and mesoscopic lengthscales

At the micron scale, where cell organelles display an amazing complexity in their shape and organization, the physical properties of a biological membrane can be better-understood using continuum models subject to thermal (stochastic) undulations. Yet, the chief orchestrators of these complex and intriguing shapes are a specialized class of membrane associating often peripheral proteins called curvature remodeling proteins (CRPs) that operate at the molecular level through specific protein-lipid interactions. We review multiscale methodologies to model these systems at the molecular as well as at the mesoscopic and cellular scales, and also present a free energy perspective of membrane remodeling through the organization and assembly of CRPs. We discuss the morphological space of nearly planar to highly curved membranes, methods to include thermal fluctuations, and review studies that model such proteins as curvature fields to describe the emergent curved morphologies. We also discuss several mesoscale models applied to a variety of cellular processes, where the phenomenological parameters (such as curvature field strength) are often mapped to models of real systems based on molecular simulations. Much insight can be gained from the calculation of free energies of membranes states with protein fields, which enable accurate mapping of the state and parameter values at which the membrane undergoes morphological transformations such as vesiculation or tubulation. By tuning the strength, anisotropy, and spatial organization of the curvature-field, one can generate a rich array of membrane morphologies that are highly relevant to shapes of several cellular organelles. We review applications of these models to budding of vesicles commonly seen in cellular signaling and trafficking processes such as clathrin mediated endocytosis, sorting by the ESCRT protein complexes, and cellular exocytosis regulated by the exocyst complex. We discuss future prospects where such models can be combined with other models for cytoskeletal assembly, and discuss their role in understanding the effects of cell membrane tension and the mechanics of the extracellular microenvironment on cellular processes.

Membrane re-modelling by BAR domain superfamily proteins via molecular and non-molecular factors

Lipid membranes are structural components of cell surfaces and intracellular organelles. Alterations in lipid membrane shape are accompanied by numerous cellular functions, including endocytosis, intracellular transport, and cell migration. Proteins containing Bin-Amphiphysin-Rvs (BAR) domains (BAR proteins) are unique, because their structures correspond to the membrane curvature, that is, the shape of the lipid membrane. BAR proteins present at high concentration determine the shape of the membrane, because BAR domain oligomers function as scaffolds that mould the membrane. BAR proteins co-operate with various molecular and non-molecular factors. The molecular factors include cytoskeletal proteins such as the regulators of actin filaments and the membrane scission protein dynamin. Lipid composition, including saturated or unsaturated fatty acid tails of phospholipids, also affects the ability of BAR proteins to mould the membrane. Non-molecular factors include the external physical forces applied to the membrane, such as tension and friction. In this mini-review, we will discuss how the BAR proteins orchestrate membrane dynamics together with various molecular and non-molecular factors.

The Unsolved Problem of How Cells Sense Micron-Scale Curvature

Membrane curvature is a fundamental feature of cells and their organelles. Much of what we know about how cells sense curved surfaces comes from studies examining nanometer-sized molecules on nanometer-scale curvatures. We are only just beginning to understand how cells recognize curved topologies at the micron scale. In this review, we provide the reader with an overview of our current understanding of how cells sense and respond to micron-scale membrane curvature.