Membrane-Binding Mechanism of a Peripheral Membrane Protein through Microsecond Molecular Dynamics Simulations

Model Mechanism for Lipid Uptake by the Human STARD2/PC-TP Phosphatidylcholine Transfer Protein

The human StAR-related lipid transfer domain protein 2 (STARD2), also known as phosphatidylcholine (PC) transfer protein, is a single-domain lipid transfer protein thought to transfer PC lipids between intracellular membranes. We performed extensive μs-long molecular dynamics simulations of STARD2 of its apo and holo forms in the presence or absence of complex lipid bilayers. The simulations in water reveal ligand-dependent conformational changes. In the 2 μs-long simulations of apo STARD2 in the presence of a lipid bilayer, we observed spontaneous reproducible PC lipid uptake into the protein hydrophobic cavity. We propose that the lipid extraction mechanism involves one to two metastable states stabilized by choline-tyrosine or choline-tryptophane cation-π interactions. Using free energy perturbation, we evaluate that PC-tyrosine cation-π interactions contribute 1.8 and 2.5 kcal/mol to the affinity of a PC-STARD2 metastable state, thus potentially providing a significant decrease of the energy barrier required for lipid desorption.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated (Myr) N-terminal matrix (MA) domain of Gag, which eventually multimerize on the membrane to form trimers and higher order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex transient dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between the MA monomers and MA-membrane remain elusive in the context of viral assembly and release. Our present study focuses on the membrane binding dynamics of a higher order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of an MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest complex protein dynamics during the formation of the immature HIV-1 lattice; while MA trimerization facilitates Myr insertion, MA trimer-trimer interactions in the immature lattice can hinder the same.

Molecular dynamics simulations of HIV-1 matrix-membrane interactions at different stages of viral maturation

Although the structural rearrangement of the membrane-bound matrix (MA) protein trimers upon HIV-1 maturation has been reported, the consequences of MA maturation on the MA-lipid interactions are not well understood. Long-timescale molecular dynamics simulations of the MA multimeric assemblies of immature and mature virus particles with our realistic asymmetric membrane model have explored MA-lipid interactions and lateral organization of lipids around MA complexes. The number of stable MA-phosphatidylserine and MA-phosphatidylinositol 4,5-bisphosphate (PIP2) interactions at the trimeric interface of the mature MA complex is observed to be greater compared to that of the immature MA complex. Our simulations identified an alternative PIP2-binding site in the immature MA complex where the multivalent headgroup of a PIP2 lipid with a greater negative charge binds to multiple basic amino acid residues such as ARG3 residues of both the MA monomers at the trimeric interface and highly basic region (HBR) residues (LYS29, LYS31) of one of the MA monomers. Our enhanced sampling simulations have explored the conformational space of phospholipids at different binding sites of the trimer-trimer interface of MA complexes that are not accessible by conventional unbiased molecular dynamics. Unlike the immature MA complex, the 2' acyl tail of two PIP2 lipids at the trimeric interface of the mature MA complex is observed to sample stable binding pockets of MA consisting of helix-4 residues. Together, our results provide molecular-level insights into the interactions of MA trimeric complexes with membrane and different lipid conformations at the specific binding sites of MA protein before and after viral maturation.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated(Myr) N-terminal matrix (MA) domain of Gag which eventually multimerize on the membrane to form trimers and higher-order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between MA monomers and MA-membrane still remain elusive. Our present study focuses on the membrane binding dynamics of a higher-order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of a MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest that MA trimerization facilitates Myr insertion, but MA trimer-trimer interactions in the lattice of immature HIV-1 particles can hinder the same. Additionally, local lipid density patterns of different lipid species provide a signature of the initial stage of lipid-domain formation upon membrane binding of the protein complex.

TOC

Measuring pico-Newton Forces with Lipid Anchors as Force Sensors in Molecular Dynamics Simulations

Binding forces between biomolecules are ubiquitous in nature but sometimes as weak as a few pico-Newtons (pN). In many cases, the binding partners are attached to biomembranes with the help of a lipid anchor. One important example are glycolipids that promote membrane adhesion through weak carbohydrate-carbohydrate binding between adjacent membranes. Here, we use molecular dynamics (MD) simulations to quantify the forces generated by bonds involving membrane-anchored molecules. We introduce a method in which the protrusion of the lipid anchors from the membrane acts as the force sensor. Our results with two different glycolipids reveal binding forces of up to 20 pN and corroborate the recent notion that carbohydrate-carbohydrate interactions are generic rather than specific.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Modeling the Membrane Binding Mechanism of a Lipid Transport Protein Osh4 to Single Membranes

All atom molecular dynamics simulations are utilized to unravel the binding mechanism of yeast oxysterol binding protein (Osh4) to model membranes with varying anionic lipid concentration using all-atom (AA) and the highly mobile membrane mimetic (HMMM) representations. For certain protein-lipid interactions, an improved force field description is used (CUFIX) to accurately describe lipid-protein electrostatic interactions. Our detailed computational studies have identified a single, β-crease orientated, membrane bound conformation of Osh4 for all anionic membranes. The penetration of the PHE-239 residue below the membrane phosphate plane is the characteristic signature of the membrane-bound state of Osh4. As the phenylalanine loop anchors itself deeply in the membrane; the other regions of the Osh4, namely, ALPS motif, ß6- ß7 loop, ß14- ß15 loop and ß16- ß17 loop, maximize their contact with the membrane. Furthermore, loose lipid packing and higher mobility of HMMM enables stronger association of ALPS motif with the membrane lipids through its hydrophobic surface and after the HMMM is converted to AA and equilibrated the binding is 2-3 times stronger compared to simulations started with the AA representation yielding the major importance of the ALPS motif to binding. Quantitative estimation of binding energy revealed that the phenylalanine loop plays a crucial role in stable membrane attachment of Osh4 and contributes significantly toward overall binding process. The CUFIX parameters provide a more balanced picture of hydrophobic and electrostatic interactions between the protein and the membrane which differs from our past work that showed salt bridges alone stabilized Osh4-membrane contact. Our study provides a comprehensive picture of the binding mechanism of Osh4 with model single membranes and thus, understanding of the initial interactions is important for elucidating the biological function of this protein to shuttle lipids between organelles.

Coordination of inter-organelle communication and lipid fluxes by OSBP-related proteins

Oxysterol-binding protein (OSBP) and OSBP-related proteins (ORPs) constitute one of the largest families of lipid-binding/transfer proteins (LTPs) in eukaryotes. The current view is that many of them mediate inter-organelle lipid transfer over membrane contact sites (MCS). The transfer occurs in several cases in a 'counter-current' fashion: A lipid such as cholesterol or phosphatidylserine (PS) is transferred against its concentration gradient driven by transport of a phosphoinositide in the opposite direction. In this way ORPs are envisioned to maintain the distinct organelle lipid compositions, with impacts on multiple organelle functions. However, the functions of ORPs extend beyond lipid homeostasis to regulation of processes such as cell survival, proliferation and migration. Important expanding areas of mammalian ORP research include their roles in viral and bacterial infections, cancers, and neuronal function. The yeast OSBP homologue (Osh) proteins execute multifaceted functions in sterol and glycerophospholipid homeostasis, post-Golgi vesicle transport, phosphatidylinositol-4-phosphate, sphingolipid and target of rapamycin (TOR) signalling, and cell cycle control. These observations identify ORPs as lipid transporters and coordinators of signals with an unforeseen variety of cellular processes. Understanding their activities not only enlightens the biology of the living cell but also allows their employment as targets of new therapeutic approaches for disease.

Computational Approaches to Investigate and Design Lipid-binding Domains for Membrane Biosensing

Association of proteins with cellular membranes is critical for signaling and membrane trafficking processes. Many peripheral lipid-binding domains have been identified in the last few decades and have been investigated for their specific lipid-sensing properties using traditional in vivo and in vitro studies. However, several knowledge gaps remain owing to intrinsic limitations of these methodologies. Thus, novel approaches are necessary to further our understanding in lipid-protein biology. This review briefly discusses lipid-binding domains that act as specific lipid biosensors and provides a broad perspective on the computational approaches such as molecular dynamics (MD) simulations and machine learning (ML)-based techniques that can be used to study protein-membrane interactions. We also highlight the need for de novo design of proteins that elicit specific lipid-binding properties.

Membrane hydrophobicity determines the activation free energy of passive lipid transport

The collective behavior of lipids with diverse chemical and physical features determines a membrane's thermodynamic properties. Yet, the influence of lipid physicochemical properties on lipid dynamics, in particular interbilayer transport, remains underexplored. Here, we systematically investigate how the activation free energy of passive lipid transport depends on lipid chemistry and membrane phase. Through all-atom molecular dynamics simulations of 11 chemically distinct glycerophospholipids, we determine how lipid acyl chain length, unsaturation, and headgroup influence the free energy barriers for two elementary steps of lipid transport: lipid desorption, which is rate limiting, and lipid insertion into a membrane. Consistent with previous experimental measurements, we find that lipids with longer, saturated acyl chains have increased activation free energies compared to lipids with shorter, unsaturated chains. Lipids with different headgroups exhibit a range of activation free energies; however, no clear trend based solely on chemical structure can be identified, mirroring difficulties in the interpretation of previous experimental results. Compared to liquid-crystalline phase membranes, gel phase membranes exhibit substantially increased free energy barriers. Overall, we find that the activation free energy depends on a lipid's local hydrophobic environment in a membrane and that the free energy barrier for lipid insertion depends on a membrane's interfacial hydrophobicity. Both of these properties can be altered through changes in lipid acyl chain length, lipid headgroup, and membrane phase. Thus, the rate of lipid transport can be tuned through subtle changes in local membrane composition and order, suggesting an unappreciated role for nanoscale membrane domains in regulating cellular lipid dynamics.

Molecular Bases of the Membrane Association Mechanism Potentiating Antibiotic Resistance by New Delhi Metallo-β-lactamase 1

Resistance to last-resort carbapenem antibiotics is an increasing threat to human health, as it critically limits therapeutic options. Metallo-β-lactamases (MBLs) are the largest family of carbapenemases, enzymes that inactivate these drugs. Among MBLs, New Delhi metallo-β-lactamase 1 (NDM-1) has experienced the fastest and largest worldwide dissemination. This success has been attributed to the fact that NDM-1 is a lipidated protein anchored to the outer membrane of bacteria, while all other MBLs are soluble periplasmic enzymes. By means of a combined experimental and computational approach, we show that NDM-1 interacts with the surface of bacterial membranes in a stable, defined conformation, in which the active site is not occluded by the bilayer. Although the lipidation is required for a long-lasting interaction, the globular domain of NDM-1 is tuned to interact specifically with the outer bacterial membrane. In contrast, this affinity is not observed for VIM-2, a natively soluble MBL. Finally, we identify key residues involved in the membrane interaction with NDM-1, which constitute potential targets for developing therapeutic strategies able to combat resistance granted by this enzyme.

Binding mechanism of the matrix domain of HIV-1 gag on lipid membranes

Specific protein-lipid interactions are critical for viral assembly. We present a molecular dynamics simulation study on the binding mechanism of the membrane targeting domain of HIV-1 Gag protein. The matrix (MA) domain drives Gag onto the plasma membrane through electrostatic interactions at its highly-basic-region (HBR), located near the myristoylated (Myr) N-terminus of the protein. Our study suggests Myr insertion is involved in the sorting of membrane lipids around the protein-binding site to prepare it for viral assembly. Our realistic membrane models confirm interactions with PIP₂ and PS lipids are highly favored around the HBR and are strong enough to keep the protein bound even without Myr insertion. We characterized Myr insertion events from microsecond trajectories and examined the membrane response upon initial membrane targeting by MA. Insertion events only occur with one of the membrane models, showing a combination of surface charge and internal membrane structure modulate this process.

Structural and Functional Specialization of OSBP-Related Proteins

Lipids are precisely distributed in the eukaryotic cell where they help to define organelle identity and function, in addition to their structural role. Once synthesized, many lipids must be delivered to other compartments by non-vesicular routes, a process that is undertaken by proteins called Lipid Transfer Proteins (LTPs). OSBP and the closely-related ORP and Osh proteins constitute a major, evolutionarily conserved family of LTPs in eukaryotes. Most of these target one or more subcellular regions, and membrane contact sites in particular, where two organelle membranes are in close proximity. It was initially thought that such proteins were strictly dedicated to sterol sensing or transport. However, over the last decade, numerous studies have revealed that these proteins have many more functions, and we have expanded our understanding of their mechanisms. In particular, many of them are lipid exchangers that exploit PI(4)P or possibly other phosphoinositide gradients to directionally transfer sterol or PS between two compartments. Importantly, these transfer activities are tightly coupled to processes such as lipid metabolism, cellular signalling and vesicular trafficking. This review describes the molecular architecture of OSBP/ORP/Osh proteins, showing how their specific structural features and internal configurations impart unique cellular functions.

Clustering and dynamics of crowded proteins near membranes and their influence on membrane bending

Atomistic molecular dynamics simulations of concentrated protein solutions in the presence of a phospholipid bilayer are presented to gain insights into the dynamics and interactions at the cytosol-membrane interface. The main finding is that proteins that are not known to specifically interact with membranes are preferentially excluded from the membrane, leaving a depletion zone near the membrane surface. As a consequence, effective protein concentrations increase, leading to increased protein contacts and clustering, whereas protein diffusion becomes faster near the membrane for proteins that do occasionally enter the depletion zone. Since protein-membrane contacts are infrequent and short-lived in this study, the structure of the lipid bilayer remains largely unaffected by the crowded protein solution, but when proteins do contact lipid head groups, small but statistically significant local membrane curvature is induced, on average.

An electrostatic switching mechanism to control the lipid transfer activity of Osh6p

A central assumption is that lipid transfer proteins (LTPs) bind transiently to organelle membranes to distribute lipids in the eukaryotic cell. Osh6p and Osh7p are yeast LTPs that transfer phosphatidylserine (PS) from the endoplasmic reticulum (ER) to the plasma membrane (PM) via PS/phosphatidylinositol-4-phosphate (PI4P) exchange cycles. It is unknown how, at each cycle, they escape from the electrostatic attraction of the PM, highly anionic, to return to the ER. Using cellular and in vitro approaches, we show that Osh6p reduces its avidity for anionic membranes once it captures PS or PI4P, due to a molecular lid closing its lipid-binding pocket. Thus, Osh6p maintains its transport activity between ER- and PM-like membranes. Further investigations reveal that the lid governs the membrane docking and activity of Osh6p because it is anionic. Our study unveils how an LTP self-limits its residency time on membranes, via an electrostatic switching mechanism, to transfer lipids efficiently.

Osh6p and Osh7p are yeast lipid transfer proteins (LTPs) that must transiently interact with membranes but how they escape from the electrostatic attraction of the plasma membrane is unclear. Here authors show that Osh6p reduces its avidity for anionic membranes once it captures PS or PI4P, due to a molecular lid closing its lipid-binding pocket.

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.

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.

Structural determinants of the interaction between influenza A virus matrix protein M1 and lipid membranes

Influenza A virus is a pathogen responsible for severe seasonal epidemics threatening human and animal populations every year. One of the ten major proteins encoded by the viral genome, the matrix protein M1, is abundantly produced in infected cells and plays a structural role in determining the morphology of the virus. During assembly of new viral particles, M1 is recruited to the host cell membrane where it associates with lipids and other viral proteins. The structure of M1 is only partially known. In particular, structural details of M1 interactions with the cellular plasma membrane as well as M1-protein interactions and multimerization have not been clarified, yet. In this work, we employed a set of complementary experimental and theoretical tools to tackle these issues. Using raster image correlation, surface plasmon resonance and circular dichroism spectroscopies, we quantified membrane association and oligomerization of full-length M1 and of different genetically engineered M1 constructs (i.e., N- and C-terminally truncated constructs and a mutant of the polybasic region, residues 95-105). Furthermore, we report novel information on structural changes in M1 occurring upon binding to membranes. Our experimental results are corroborated by an all-atom model of the full-length M1 protein bound to a negatively charged lipid bilayer.

Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins

Within eukaryotic cells, biochemical reactions need to be organized on the surface of membrane compartments that use distinct lipid constituents to dynamically modulate the functions of integral proteins or influence the selective recruitment of peripheral membrane effectors. As a result of these complex interactions, a variety of human pathologies can be traced back to improper communication between proteins and membrane surfaces; either due to mutations that directly alter protein structure or as a result of changes in membrane lipid composition. Among the known structural lipids found in cellular membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the membrane-anchored precursor of low-abundance regulatory lipids, the polyphosphoinositides (PPIn), which have restricted distributions within specific subcellular compartments. The ability of PPIn lipids to function as signaling platforms relies on both non-specific electrostatic interactions and the selective stereospecific recognition of PPIn headgroups by specialized protein folds. In this chapter, we will attempt to summarize the structural diversity of modular PPIn-interacting domains that facilitate the reversible recruitment and conformational regulation of peripheral membrane proteins. Outside of protein folds capable of capturing PPIn headgroups at the membrane interface, recent studies detailing the selective binding and bilayer extraction of PPIn species by unique functional domains within specific families of lipid-transfer proteins will also be highlighted. Overall, this overview will help to outline the fundamental physiochemical mechanisms that facilitate localized interactions between PPIn lipids and the wide-variety of PPIn-binding proteins that are essential for the coordinate regulation of cellular metabolism and membrane dynamics.

Sensitivity of peripheral membrane proteins to the membrane context: A case study of phosphatidylserine and the TIM proteins

There is a diverse class of peripheral membrane-binding proteins that specifically bind phosphatidylserine (PS), a lipid that signals apoptosis or cell fusion depending on the membrane context of its presentation. PS-receptors are specialized for particular PS-presenting pathways, indicating that they might be sensitive to the membrane context. In this review, we describe a combination of thermodynamic, structural, and computational techniques that can be used to investigate the mechanisms underlying this sensitivity. As an example, we focus on three PS-receptors of the T-cell Immunoglobulin and Mucin containing (TIM) protein family, which we have previously shown to differ in their sensitivity to PS surface density.

A Role for Weak Electrostatic Interactions in Peripheral Membrane Protein Binding

Bacillus thuringiensis phosphatidylinositol-specific phospholipase C (BtPI-PLC) is a secreted virulence factor that binds specifically to phosphatidylcholine (PC) bilayers containing negatively charged phospholipids. BtPI-PLC carries a negative net charge and its interfacial binding site has no obvious cluster of basic residues. Continuum electrostatic calculations show that, as expected, nonspecific electrostatic interactions between BtPI-PLC and membranes vary as a function of the fraction of anionic lipids present in the bilayers. Yet they are strikingly weak, with a calculated ΔGel below 1 kcal/mol, largely due to a single lysine (K44). When K44 is mutated to alanine, the equilibrium dissociation constant for small unilamellar vesicles increases more than 50 times (∼2.4 kcal/mol), suggesting that interactions between K44 and lipids are not merely electrostatic. Comparisons of molecular-dynamics simulations performed using different lipid compositions reveal that the bilayer composition does not affect either hydrogen bonds or hydrophobic contacts between the protein interfacial binding site and bilayers. However, the occupancies of cation-π interactions between PC choline headgroups and protein tyrosines vary as a function of PC content. The overall contribution of basic residues to binding affinity is also context dependent and cannot be approximated by a rule-of-thumb value because these residues can contribute to both nonspecific electrostatic and short-range protein-lipid interactions. Additionally, statistics on the distribution of basic amino acids in a data set of membrane-binding domains reveal that weak electrostatics, as observed for BtPI-PLC, might be a less unusual mechanism for peripheral membrane binding than is generally thought.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (S_CD), and spin-lattice relaxation time (T₁)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipid's mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model

Slow diffusion of the lipids in conventional all-atom simulations of membrane systems makes it difficult to sample large rearrangements of lipids and protein-lipid interactions. Recently, Tajkhorshid and co-workers developed the highly mobile membrane-mimetic (HMMM) model with accelerated lipid motion by replacing the lipid tails with small organic molecules. The HMMM model provides accelerated lipid diffusion by one to two orders of magnitude, and is particularly useful in studying membrane-protein associations. However, building an HMMM simulation system is not easy, as it requires sophisticated treatment of the lipid tails. In this study, we have developed CHARMM-GUI HMMM Builder (http://www.charmm-gui.org/input/hmmm) to provide users with ready-to-go input files for simulating HMMM membrane systems with/without proteins. Various lipid-only and protein-lipid systems are simulated to validate the qualities of the systems generated by HMMM Builder with focus on the basic properties and advantages of the HMMM model. HMMM Builder supports all lipid types available in CHARMM-GUI and also provides a module to convert back and forth between an HMMM membrane and a full-length membrane. We expect HMMM Builder to be a useful tool in studying membrane systems with enhanced lipid diffusion.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Peripheral membrane proteins: Tying the knot between experiment and computation

Experimental biology has contributed to answer questions about the morphology of a system and how molecules organize themselves to maintain a healthy functional cell. Single-molecule techniques, optical and magnetic experiments, and fluorescence microscopy have come a long way to probe structural and dynamical information at multiple scales. However, some details are simply too small or the processes are too short-lived to detect by experiments. Computational biology provides a bridge to understand experimental results at the molecular level, makes predictions that have not been seen in vivo, and motivates new fields of research. This review focuses on the advances on peripheral membrane proteins (PMPs) studies; what is known about their interaction with membranes, their role in cell biology, and some limitations that both experiment and computation still have to overcome to gain better structural and functional understanding of these PMPs. As many recent reviews have acknowledged, interdisciplinary efforts between experiment and computation are needed in order to have useful models that lead future directions in the study of PMPs. We present new results of a case study on a PMP that behaves as an intricate machine controlling lipid homeostasis between cellular organelles, Osh4 in yeast Saccharomyces cerevisiae. Molecular dynamics simulations were run to examine the interaction between the protein and membrane models that reflect the lipid diversity of the endoplasmic reticulum and trans-Golgi membranes. Our study is consistent with experimental data showing several residues that interact to smaller or larger extent with the bilayer upon stable binding (~200 ns into the trajectory). We identified PHE239 as a key residue stabilizing the protein-membrane interaction along with two other binding regions, the ALPS-like motif and the β6-β7 loops in the mouth region of the protein. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Two homologous neutrophil serine proteases bind to POPC vesicles with different affinities: When aromatic amino acids matter

Neutrophil serine proteases Proteinase 3 (PR3) and human neutrophil elastase (HNE) are homologous antibiotic serine proteases of the polymorphonuclear neutrophils. Despite sharing a 56% sequence identity they have been shown to have different functions and localizations in the neutrophils. In particular, and in contrast to HNE, PR3 has been detected at the outer leaflet of the plasma membrane and its membrane expression is a risk factor in a number of chronic inflammatory diseases. Although a plethora of studies performed in various cell-based assays have been reported, the mechanism by which PR3, and possibly HNE bind to simple membrane models remains unclear. We used surface plasmon resonance (SPR) experiments to measure and compare the affinity of PR3 and HNE for large unilamellar vesicles composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). We also conducted 500-nanosecond long molecular dynamics simulations of each enzyme at the surface of a POPC bilayer to map the interactions between proteins and lipids and rationalize the difference in affinity observed in the SPR experiment. We find that PR3 binds strongly to POPC large unilamellar vesicles (Kd=9.2×10(-7)M) thanks to the insertion of three phenylalanines, one tryptophan and one leucine beyond the phosphate groups of the POPC lipids. HNE binds in a significantly weaker manner (Kd>10(-5)M) making mostly electrostatic interactions via lysines and arginines and inserting only one leucine between the hydrophobic lipid tails. Our results support the early reports that PR3, unlike HNE, is able to directly and strongly anchor directly to the neutrophil membrane.

A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Kučerka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (DB), hydrophobic thickness (DC), and lipid volume (VL). The only minor discrepancy between simulation and experiment is a measure of (DB-DHH)/2 where DHH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.

Cation-π interactions as lipid-specific anchors for phosphatidylinositol-specific phospholipase C

Amphitropic proteins, such as the virulence factor phosphatidylinositol-specific phospholipase C (PI-PLC) from Bacillus thuringiensis , often depend on lipid-specific recognition of target membranes. However, the recognition mechanisms for zwitterionic lipids, such as phosphatidylcholine, which is enriched in the outer leaflet of eukaryotic cells, are not well understood. A 500 ns long molecular dynamics simulation of PI-PLC at the surface of a lipid bilayer revealed a strikingly high number of interactions between tyrosines at the interfacial binding site and lipid choline groups with structures characteristic of cation-π interactions. Membrane affinities of PI-PLC tyrosine variants mostly tracked the simulation results, falling into two classes: (i) those with minor losses in affinity, Kd(mutant)/Kd(wild-type) ≤ 5 and (ii) those where the apparent Kd was 50-200 times higher than wild-type. Estimating ΔΔG for these Tyr/PC interactions from the apparent Kd values reveals that the free energy associated with class I is ~1 kcal/mol, comparable to the value predicted by the Wimley-White hydrophobicity scale. In contrast, removal of class II tyrosines has a higher energy cost: ~2.5 kcal/mol toward pure PC vesicles. These higher energies correlate well with the occupancy of the cation-π adducts throughout the MD simulation. Together, these results strongly indicate that PI-PLC interacts with PC headgroups via cation-π interactions with tyrosine residues and suggest that cation-π interactions at the interface may be a mechanism for specific lipid recognition by amphitropic and membrane proteins.