Multiscale simulations of protein-facilitated membrane remodeling

SPICA Force Field for Proteins and Peptides

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

On the Role of Curved Membrane Nanodomains, and Passive and Active Skeleton Forces in the Determination of Cell Shape and Membrane Budding

Biological membranes are composed of isotropic and anisotropic curved nanodomains. Anisotropic membrane components, such as Bin/Amphiphysin/Rvs (BAR) superfamily protein domains, could trigger/facilitate the growth of membrane tubular protrusions, while isotropic curved nanodomains may induce undulated (necklace-like) membrane protrusions. We review the role of isotropic and anisotropic membrane nanodomains in stability of tubular and undulated membrane structures generated or stabilized by cyto- or membrane-skeleton. We also describe the theory of spontaneous self-assembly of isotropic curved membrane nanodomains and derive the critical concentration above which the spontaneous necklace-like membrane protrusion growth is favorable. We show that the actin cytoskeleton growth inside the vesicle or cell can change its equilibrium shape, induce higher degree of segregation of membrane nanodomains or even alter the average orientation angle of anisotropic nanodomains such as BAR domains. These effects may indicate whether the actin cytoskeleton role is only to stabilize membrane protrusions or to generate them by stretching the vesicle membrane. Furthermore, we demonstrate that by taking into account the in-plane orientational ordering of anisotropic membrane nanodomains, direct interactions between them and the extrinsic (deviatoric) curvature elasticity, it is possible to explain the experimentally observed stability of oblate (discocyte) shapes of red blood cells in a broad interval of cell reduced volume. Finally, we present results of numerical calculations and Monte-Carlo simulations which indicate that the active forces of membrane skeleton and cytoskeleton applied to plasma membrane may considerably influence cell shape and membrane budding.

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.

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.

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.

Scaffolding the cup-shaped double membrane in autophagy

Autophagy is a physiological process for the recycling and degradation of cellular materials. Forming the autophagosome from the phagophore, a cup-shaped double-membrane vesicle, is a critical step in autophagy. The origin of the cup shape of the phagophore is poorly understood. In yeast, fusion of a small number of Atg9-containing vesicles is considered a key step in autophagosome biogenesis, aided by Atg1 complexes (ULK1 in mammals) localized at the preautophagosomal structure (PAS). In particular, the S-shaped Atg17-Atg31-Atg29 subcomplex of Atg1 is critical for phagophore nucleation at the PAS. To study this process, we simulated membrane remodeling processes in the presence and absence of membrane associated Atg17. We show that at least three vesicles need to fuse to induce the phagophore shape, consistent with experimental observations. However, fusion alone is not sufficient. Interactions with 34-nm long, S-shaped Atg17 complexes are required to overcome a substantial kinetic barrier in the transition to the cup-shaped phagophore. Our finding rationalizes the recruitment of Atg17 complexes to the yeast PAS, and their unusual shape. In control simulations without Atg17, with weakly binding Atg17, or with straight instead of S-shaped Atg17, the membrane shape transition did not occur. We confirm the critical role of Atg17-membrane interactions experimentally by showing that mutations of putative membrane interaction sites result in reduction or loss of autophagic activity in yeast. Fusion of a small number of vesicles followed by Atg17-guided membrane shape-remodeling thus emerges as a viable route to phagophore formation.

Systems biology of cellular membranes: a convergence with biophysics

Systems biology and systems medicine have played an important role in the last two decades in shaping our understanding of biological processes. While systems biology is synonymous with network maps and '-omics' approaches, it is not often associated with mechanical processes. Here, we make the case for considering the mechanical and geometrical aspects of biological membranes as a key step in pushing the frontiers of systems biology of cellular membranes forward. We begin by introducing the basic components of cellular membranes, and highlight their dynamical aspects. We then survey the functions of the plasma membrane and the endomembrane system in signaling, and discuss the role and origin of membrane curvature in these diverse cellular processes. We further give an overview of the experimental and modeling approaches to study membrane phenomena. We close with a perspective on the converging futures of systems biology and membrane biophysics, invoking the need to include physical variables such as location and geometry in the study of cellular membranes. WIREs Syst Biol Med 2017, 9:e1386. doi: 10.1002/wsbm.1386 For further resources related to this article, please visit the WIREs website.

On the role of external force of actin filaments in the formation of tubular protrusions of closed membrane shapes with anisotropic membrane components

Biological membranes are composed of different components and there is no a priori reason to assume that all components are isotropic. It was previously shown that the anisotropic properties of membrane components may explain the stability of membrane tubular protrusions even without the application of external force. Our theoretical study focuses on the role of anisotropic membrane components in the stability of membrane tubular structures generated or stabilized by actin filaments. We show that the growth of the actin cytoskeleton inside the vesicle can induce the partial lateral segregation of different membrane components. The entropy of mixing of membrane components hinders the total lateral segregation of the anisotropic and isotropic membrane components. Self-assembled aggregates formed by anisotropic membrane components facilitate the growth of long membrane tubular protrusions. Protrusive force generated by actin filaments favors strong segregation of membrane components by diminishing the opposing effect of mixing entropy.