1
|
Talandashti R, Moqadam M, Reuter N. Model Mechanism for Lipid Uptake by the Human STARD2/PC-TP Phosphatidylcholine Transfer Protein. J Phys Chem Lett 2024; 15:8287-8295. [PMID: 39143857 PMCID: PMC11331517 DOI: 10.1021/acs.jpclett.4c01743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2024] [Revised: 07/11/2024] [Accepted: 07/16/2024] [Indexed: 08/16/2024]
Abstract
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.
Collapse
Affiliation(s)
- Reza Talandashti
- Department
of Chemistry, University of Bergen, Bergen 5020, Norway
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5020, Norway
| | - Mahmoud Moqadam
- Department
of Chemistry, University of Bergen, Bergen 5020, Norway
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5020, Norway
| | - Nathalie Reuter
- Department
of Chemistry, University of Bergen, Bergen 5020, Norway
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5020, Norway
| |
Collapse
|
2
|
Banerjee P, Monje-Galvan V, Voth GA. Cooperative Membrane Binding of HIV-1 Matrix Proteins. J Phys Chem B 2024; 128:2595-2606. [PMID: 38477117 PMCID: PMC10962350 DOI: 10.1021/acs.jpcb.3c06222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 02/24/2024] [Accepted: 02/27/2024] [Indexed: 03/14/2024]
Abstract
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.
Collapse
Affiliation(s)
- Puja Banerjee
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| | | | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Banerjee P, Qu K, Briggs JAG, Voth GA. Molecular dynamics simulations of HIV-1 matrix-membrane interactions at different stages of viral maturation. Biophys J 2024; 123:389-406. [PMID: 38196190 PMCID: PMC10870173 DOI: 10.1016/j.bpj.2024.01.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 12/05/2023] [Accepted: 01/04/2024] [Indexed: 01/11/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Puja Banerjee
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Kun Qu
- Infectious Diseases Translational Research Programme, Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore
| | - John A G Briggs
- Department of Cell and Virus Structure, Max Planck Institute of Biochemistry, Planegg, Germany
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
4
|
Banerjee P, Monje-Galvan V, Voth GA. Cooperative Membrane Binding of HIV-1 Matrix Proteins. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.22.559012. [PMID: 37790356 PMCID: PMC10542177 DOI: 10.1101/2023.09.22.559012] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
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
Collapse
|
5
|
Kav B, Weikl TR, Schneck E. Measuring pico-Newton Forces with Lipid Anchors as Force Sensors in Molecular Dynamics Simulations. J Phys Chem B 2023; 127:4081-4089. [PMID: 37127845 PMCID: PMC10184124 DOI: 10.1021/acs.jpcb.3c00063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Batuhan Kav
- Max Planck Institute of Colloids and Interfaces, 14467, Potsdam, Germany
- Institute of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, 52428 Jülich, Germany
| | - Thomas R Weikl
- Max Planck Institute of Colloids and Interfaces, 14467, Potsdam, Germany
| | - Emanuel Schneck
- Max Planck Institute of Colloids and Interfaces, 14467, Potsdam, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| |
Collapse
|
6
|
Larsen A, John L, Sansom M, Corey R. Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us? Biosci Rep 2022; 42:BSR20211406. [PMID: 35297484 PMCID: PMC9008707 DOI: 10.1042/bsr20211406] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 03/14/2022] [Accepted: 03/16/2022] [Indexed: 12/04/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Laura H. John
- Department of Biochemistry, University of Oxford, Oxford, U.K
| | | | - Robin A. Corey
- Department of Biochemistry, University of Oxford, Oxford, U.K
| |
Collapse
|
7
|
Karmakar S, Klauda JB. Modeling the Membrane Binding Mechanism of a Lipid Transport Protein Osh4 to Single Membranes. Biophys J 2022; 121:1560-1575. [PMID: 35247338 PMCID: PMC9072576 DOI: 10.1016/j.bpj.2022.03.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Revised: 09/13/2021] [Accepted: 03/01/2022] [Indexed: 11/26/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering; Biophysics Graduate Program University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
8
|
Arora A, Taskinen JH, Olkkonen VM. Coordination of inter-organelle communication and lipid fluxes by OSBP-related proteins. Prog Lipid Res 2022; 86:101146. [PMID: 34999137 DOI: 10.1016/j.plipres.2022.101146] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 12/10/2021] [Accepted: 01/03/2022] [Indexed: 12/31/2022]
Abstract
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.
Collapse
Affiliation(s)
- Amita Arora
- Minerva Foundation Institute for Medical Research, and Department of Anatomy, Faculty of Medicine, University of Helsinki, Finland
| | - Juuso H Taskinen
- Minerva Foundation Institute for Medical Research, and Department of Anatomy, Faculty of Medicine, University of Helsinki, Finland
| | - Vesa M Olkkonen
- Minerva Foundation Institute for Medical Research, and Department of Anatomy, Faculty of Medicine, University of Helsinki, Finland.
| |
Collapse
|
9
|
Srinivasan S, Vanni S. Computational Approaches to Investigate and Design Lipid-binding Domains for Membrane Biosensing. Chimia (Aarau) 2021; 75:1031-1036. [PMID: 34920773 DOI: 10.2533/chimia.2021.1031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Switzerland;,
| |
Collapse
|
10
|
Rogers JR, Espinoza Garcia G, Geissler PL. Membrane hydrophobicity determines the activation free energy of passive lipid transport. Biophys J 2021; 120:3718-3731. [PMID: 34302793 PMCID: PMC8456290 DOI: 10.1016/j.bpj.2021.07.016] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 07/07/2021] [Accepted: 07/19/2021] [Indexed: 12/16/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Julia R Rogers
- Department of Chemistry, University of California Berkeley, Berkeley, California.
| | | | - Phillip L Geissler
- Department of Chemistry, University of California Berkeley, Berkeley, California; Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California.
| |
Collapse
|
11
|
Moqadam M, Tubiana T, Moutoussamy EE, Reuter N. Membrane models for molecular simulations of peripheral membrane proteins. ADVANCES IN PHYSICS: X 2021. [DOI: 10.1080/23746149.2021.1932589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Affiliation(s)
- Mahmoud Moqadam
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Thibault Tubiana
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Emmanuel E. Moutoussamy
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
- Department of Biological Sciences, University of Bergen, Bergen, Norway
| | - Nathalie Reuter
- Department of Chemistry, University of Bergen, Bergen, Norway
- Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| |
Collapse
|
12
|
Prunotto A, Bahr G, González LJ, Vila AJ, Dal Peraro M. Molecular Bases of the Membrane Association Mechanism Potentiating Antibiotic Resistance by New Delhi Metallo-β-lactamase 1. ACS Infect Dis 2020; 6:2719-2731. [PMID: 32865963 DOI: 10.1021/acsinfecdis.0c00341] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Alessio Prunotto
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Guillermo Bahr
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Lisandro J. González
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Alejandro J. Vila
- Instituto de Biología Molecular y Celular de Rosario (IBR, CONICET-UNR), S2000EXF Rosario, Argentina
- Area Biofísica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
| | - Matteo Dal Peraro
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
13
|
Monje-Galvan V, Voth GA. Binding mechanism of the matrix domain of HIV-1 gag on lipid membranes. eLife 2020; 9:58621. [PMID: 32808928 PMCID: PMC7476761 DOI: 10.7554/elife.58621] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Accepted: 08/14/2020] [Indexed: 01/10/2023] Open
Abstract
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 PIP2 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.
Collapse
Affiliation(s)
- Viviana Monje-Galvan
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and The James Franck Institute, The University of Chicago, Chicago, United States
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and The James Franck Institute, The University of Chicago, Chicago, United States
| |
Collapse
|
14
|
Delfosse V, Bourguet W, Drin G. Structural and Functional Specialization of OSBP-Related Proteins. ACTA ACUST UNITED AC 2020. [DOI: 10.1177/2515256420946627] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
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.
Collapse
Affiliation(s)
- Vanessa Delfosse
- Centre de Biochimie Structurale, Inserm, CNRS, Univ Montpellier, Montpellier, France
| | - William Bourguet
- Centre de Biochimie Structurale, Inserm, CNRS, Univ Montpellier, Montpellier, France
| | - Guillaume Drin
- CNRS, Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d’Azur, Valbonne, France
| |
Collapse
|
15
|
Clustering and dynamics of crowded proteins near membranes and their influence on membrane bending. Proc Natl Acad Sci U S A 2019; 116:24562-24567. [PMID: 31740611 DOI: 10.1073/pnas.1910771116] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
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.
Collapse
|
16
|
An electrostatic switching mechanism to control the lipid transfer activity of Osh6p. Nat Commun 2019; 10:3926. [PMID: 31477717 PMCID: PMC6718676 DOI: 10.1038/s41467-019-11780-y] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Accepted: 08/01/2019] [Indexed: 02/01/2023] Open
Abstract
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.
Collapse
|
17
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 422] [Impact Index Per Article: 84.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
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.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
18
|
Enkavi G, Javanainen M, Kulig W, Róg T, Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem Rev 2019; 119:5607-5774. [PMID: 30859819 PMCID: PMC6727218 DOI: 10.1021/acs.chemrev.8b00538] [Citation(s) in RCA: 175] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Indexed: 12/23/2022]
Abstract
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.
Collapse
Affiliation(s)
- Giray Enkavi
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Matti Javanainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Institute
of Organic Chemistry and Biochemistry of the Czech Academy
of Sciences, Flemingovo naḿesti 542/2, 16610 Prague, Czech Republic
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Waldemar Kulig
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Tomasz Róg
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Ilpo Vattulainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
- MEMPHYS-Center
for Biomembrane Physics
| |
Collapse
|
19
|
Muller MP, Jiang T, Sun C, Lihan M, Pant S, Mahinthichaichan P, Trifan A, Tajkhorshid E. Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation. Chem Rev 2019; 119:6086-6161. [PMID: 30978005 PMCID: PMC6506392 DOI: 10.1021/acs.chemrev.8b00608] [Citation(s) in RCA: 153] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- Melanie P. Muller
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Tao Jiang
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Chang Sun
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Muyun Lihan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shashank Pant
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Paween Mahinthichaichan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Anda Trifan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Emad Tajkhorshid
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
20
|
Höfer CT, Di Lella S, Dahmani I, Jungnick N, Bordag N, Bobone S, Huang Q, Keller S, Herrmann A, Chiantia S. Structural determinants of the interaction between influenza A virus matrix protein M1 and lipid membranes. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2019; 1861:1123-1134. [PMID: 30902626 DOI: 10.1016/j.bbamem.2019.03.013] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 03/16/2019] [Indexed: 11/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- C T Höfer
- Institute for Biology, IRI Life Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115, Berlin, Germany
| | - S Di Lella
- Institute for Biology, IRI Life Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115, Berlin, Germany
| | - I Dahmani
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - N Jungnick
- Institute for Biology, IRI Life Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115, Berlin, Germany
| | - N Bordag
- Leibniz-Institute for Molecular Pharmacology (FMP), Biophysics of Membrane Proteins, Robert-Roessle-Str. 10, 13125 Berlin, Germany
| | - S Bobone
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Q Huang
- School of Life Sciences, Fudan University, 220 Handan Rd, WuJiaoChang, Yangpu Qu, Shanghai Shi 200433, China
| | - S Keller
- Molecular Biophysics, Technische Universität Kaiserslautern (TUK), Erwin-Schrödinger-Str. 13, 67663 Kaiserslautern, Germany
| | - A Herrmann
- Institute for Biology, IRI Life Sciences, Humboldt-Universität zu Berlin, Invalidenstraße 42, 10115, Berlin, Germany.
| | - S Chiantia
- University of Potsdam, Institute of Biochemistry and Biology, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany.
| |
Collapse
|
21
|
Pemberton JG, Balla T. Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2019; 1111:77-137. [PMID: 30483964 DOI: 10.1007/5584_2018_288] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
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.
Collapse
Affiliation(s)
- Joshua G Pemberton
- Section on Molecular Signal Transduction, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| | - Tamas Balla
- Section on Molecular Signal Transduction, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA.
| |
Collapse
|
22
|
Monje-Galvan V, Klauda JB. Preferred Binding Mechanism of Osh4’s Amphipathic Lipid-Packing Sensor Motif, Insights from Molecular Dynamics. J Phys Chem B 2018; 122:9713-9723. [DOI: 10.1021/acs.jpcb.8b07067] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
23
|
Roberts MF, Khan HM, Goldstein R, Reuter N, Gershenson A. Search and Subvert: Minimalist Bacterial Phosphatidylinositol-Specific Phospholipase C Enzymes. Chem Rev 2018; 118:8435-8473. [DOI: 10.1021/acs.chemrev.8b00208] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Mary F. Roberts
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | | | - Rebecca Goldstein
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | | | - Anne Gershenson
- Department of Biochemistry and Molecular Biology, University of Massachusetts Amherst, Amherst, Massachusetts 01003, United States
| |
Collapse
|
24
|
Kerr D, Tietjen GT, Gong Z, Tajkhorshid E, Adams EJ, Lee KYC. Sensitivity of peripheral membrane proteins to the membrane context: A case study of phosphatidylserine and the TIM proteins. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:2126-2133. [PMID: 29920237 DOI: 10.1016/j.bbamem.2018.06.010] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2018] [Revised: 06/11/2018] [Accepted: 06/12/2018] [Indexed: 12/17/2022]
Abstract
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.
Collapse
Affiliation(s)
- Daniel Kerr
- Program in Biophysical Sciences, Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, United States of America
| | - Gregory T Tietjen
- Program in Biophysical Sciences, Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, United States of America
| | - Zhiliang Gong
- Department of Chemistry, The University of Chicago, Chicago, IL, United States of America
| | - Emad Tajkhorshid
- Department of Biochemistry, Center for Biophysics and Quantitative Biology, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
| | - Erin J Adams
- Department of Biochemistry and Molecular Biology and Committee on Immunology, The University of Chicago, Chicago, IL, United States of America
| | - Ka Yee C Lee
- Program in Biophysical Sciences, Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, United States of America; Department of Chemistry, The University of Chicago, Chicago, IL, United States of America; James Franck Institute, The University of Chicago, Chicago, IL, United States of America.
| |
Collapse
|
25
|
Khan HM, He T, Fuglebakk E, Grauffel C, Yang B, Roberts MF, Gershenson A, Reuter N. A Role for Weak Electrostatic Interactions in Peripheral Membrane Protein Binding. Biophys J 2016; 110:1367-78. [PMID: 27028646 PMCID: PMC4816757 DOI: 10.1016/j.bpj.2016.02.020] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2015] [Revised: 01/08/2016] [Accepted: 02/02/2016] [Indexed: 12/01/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Hanif M Khan
- Department of Molecular Biology, University of Bergen, Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Tao He
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts
| | - Edvin Fuglebakk
- Department of Molecular Biology, University of Bergen, Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Cédric Grauffel
- Department of Molecular Biology, University of Bergen, Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway
| | - Boqian Yang
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts; Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts
| | - Mary F Roberts
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts
| | - Anne Gershenson
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, Massachusetts
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen, Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Bergen, Norway.
| |
Collapse
|
26
|
An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:3093-3104. [PMID: 27664502 DOI: 10.1016/j.bbamem.2016.09.016] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2016] [Revised: 09/14/2016] [Accepted: 09/19/2016] [Indexed: 12/21/2022]
Abstract
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 (SCD), and spin-lattice relaxation time (T1)) 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.
Collapse
|
27
|
Qi Y, Cheng X, Lee J, Vermaas JV, Pogorelov TV, Tajkhorshid E, Park S, Klauda JB, Im W. CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model. Biophys J 2016; 109:2012-22. [PMID: 26588561 DOI: 10.1016/j.bpj.2015.10.008] [Citation(s) in RCA: 78] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2015] [Revised: 10/06/2015] [Accepted: 10/09/2015] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Yifei Qi
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Xi Cheng
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Jumin Lee
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Josh V Vermaas
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois; Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois
| | - Taras V Pogorelov
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois; School of Chemical Sciences, University of Illinois at Urbana-Champaign, Urbana, Illinois; National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois
| | - Emad Tajkhorshid
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois; Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois
| | - Soohyung Park
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering and the Biophysics Program, The University of Maryland, College Park, Maryland
| | - Wonpil Im
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas, Lawrence, Kansas.
| |
Collapse
|
28
|
Pöyry S, Vattulainen I. Role of charged lipids in membrane structures - Insight given by simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2322-2333. [PMID: 27003126 DOI: 10.1016/j.bbamem.2016.03.016] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2015] [Revised: 03/09/2016] [Accepted: 03/11/2016] [Indexed: 01/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sanja Pöyry
- Department of Physics, Tampere University of Technology, POB 692, FI-33101 Tampere, Finland
| | - Ilpo Vattulainen
- Department of Physics, Tampere University of Technology, POB 692, FI-33101 Tampere, Finland; MEMPHYS - Center for Biomembrane Physics, University of Southern Denmark, Odense, Denmark; Department of Physics, University of Helsinki, POB 64, FI-00014 Helsinki, Finland.
| |
Collapse
|
29
|
Monje-Galvan V, Klauda JB. Peripheral membrane proteins: Tying the knot between experiment and computation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:1584-93. [PMID: 26903211 DOI: 10.1016/j.bbamem.2016.02.018] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2015] [Revised: 02/05/2016] [Accepted: 02/12/2016] [Indexed: 01/31/2023]
Abstract
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.
Collapse
Affiliation(s)
- Viviana Monje-Galvan
- Department of Chemical and Biomolecular Engineering, College Park, MD 20742, USA
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, College Park, MD 20742, USA; Biophysics Program, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
30
|
Lee J, Cheng X, Swails JM, Yeom MS, Eastman PK, Lemkul JA, Wei S, Buckner J, Jeong JC, Qi Y, Jo S, Pande VS, Case DA, Brooks CL, MacKerell AD, Klauda JB, Im W. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J Chem Theory Comput 2016. [PMID: 26631602 DOI: 10.1021/acs.jctc.5b00935/asset/images/large/ct-2015-00935e0005.jpeg] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/31/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jumin Lee
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Xi Cheng
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Jason M Swails
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854, United States
| | - Min Sun Yeom
- Korean Institute of Science and Technology Information , Yuseong-gu, Daejeon 305-806, Korea
| | - Peter K Eastman
- Department of Bioengineering, Stanford University , Stanford, California 94035, United States
| | - Justin A Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Shuai Wei
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Joshua Buckner
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Jong Cheol Jeong
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Harvard Medical School , Boston, Massachusetts 02215, United States
| | - Yifei Qi
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory , 9700 Cass Avenue, Building 240, Argonne, Illinois 60439, United States
| | - Vijay S Pande
- Department of Bioengineering, Stanford University , Stanford, California 94035, United States
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854, United States
| | - Charles L Brooks
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering and the Biophysics Program, University of Maryland , College Park, Maryland 20742, United States
| | - Wonpil Im
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| |
Collapse
|
31
|
Lee J, Cheng X, Swails JM, Yeom MS, Eastman PK, Lemkul JA, Wei S, Buckner J, Jeong JC, Qi Y, Jo S, Pande VS, Case DA, Brooks CL, MacKerell AD, Klauda JB, Im W. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J Chem Theory Comput 2015; 12:405-13. [PMID: 26631602 PMCID: PMC4712441 DOI: 10.1021/acs.jctc.5b00935] [Citation(s) in RCA: 2237] [Impact Index Per Article: 248.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
![]()
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.
Collapse
Affiliation(s)
- Jumin Lee
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Xi Cheng
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Jason M Swails
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854, United States
| | - Min Sun Yeom
- Korean Institute of Science and Technology Information , Yuseong-gu, Daejeon 305-806, Korea
| | - Peter K Eastman
- Department of Bioengineering, Stanford University , Stanford, California 94035, United States
| | - Justin A Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Shuai Wei
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Joshua Buckner
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Jong Cheol Jeong
- Cancer Research Institute, Beth Israel Deaconess Cancer Center, Harvard Medical School , Boston, Massachusetts 02215, United States
| | - Yifei Qi
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| | - Sunhwan Jo
- Leadership Computing Facility, Argonne National Laboratory , 9700 Cass Avenue, Building 240, Argonne, Illinois 60439, United States
| | - Vijay S Pande
- Department of Bioengineering, Stanford University , Stanford, California 94035, United States
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers University , Piscataway, New Jersey 08854, United States
| | - Charles L Brooks
- Department of Chemistry and the Biophysics Program, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering and the Biophysics Program, University of Maryland , College Park, Maryland 20742, United States
| | - Wonpil Im
- Department of Molecular Biosciences and Center for Computational Biology, The University of Kansas , Lawrence, Kansas 66047, United States
| |
Collapse
|
32
|
Two homologous neutrophil serine proteases bind to POPC vesicles with different affinities: When aromatic amino acids matter. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:3191-202. [PMID: 25218402 DOI: 10.1016/j.bbamem.2014.09.003] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2014] [Revised: 08/29/2014] [Accepted: 09/03/2014] [Indexed: 11/22/2022]
Abstract
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.
Collapse
|
33
|
A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:2520-9. [PMID: 24953542 DOI: 10.1016/j.bbamem.2014.06.010] [Citation(s) in RCA: 74] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/19/2014] [Revised: 06/07/2014] [Accepted: 06/11/2014] [Indexed: 11/21/2022]
Abstract
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.
Collapse
|
34
|
Galassi VV, Villarreal MA, Posada V, Montich GG. Interactions of the fatty acid-binding protein ReP1-NCXSQ with lipid membranes. Influence of the membrane electric field on binding and orientation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:910-20. [DOI: 10.1016/j.bbamem.2013.11.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2013] [Revised: 11/12/2013] [Accepted: 11/13/2013] [Indexed: 11/17/2022]
|
35
|
Grauffel C, Yang B, He T, Roberts MF, Gershenson A, Reuter N. Cation-π interactions as lipid-specific anchors for phosphatidylinositol-specific phospholipase C. J Am Chem Soc 2013; 135:5740-50. [PMID: 23506313 PMCID: PMC3797534 DOI: 10.1021/ja312656v] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
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.
Collapse
Affiliation(s)
- Cédric Grauffel
- Department of Molecular Biology, University of Bergen, Norway
- Computational Biology Unit, Uni Research, Bergen, Norway
| | - Boqian Yang
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, U.S.A
- Department of Chemistry, Boston College, Chestnut Hill, U.S.A
| | - Tao He
- Department of Chemistry, Boston College, Chestnut Hill, U.S.A
| | - Mary F. Roberts
- Department of Chemistry, Boston College, Chestnut Hill, U.S.A
| | - Anne Gershenson
- Department of Biochemistry and Molecular Biology, University of Massachusetts, Amherst, U.S.A
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen, Norway
- Computational Biology Unit, Uni Research, Bergen, Norway
| |
Collapse
|