1
|
Keller F, Alavizargar A, Wedlich-Söldner R, Heuer A. The impact of bilayer composition on the dimerization properties of the Slg1 stress sensor TMD from a multiscale analysis. Phys Chem Chem Phys 2023; 25:1299-1309. [PMID: 36533706 DOI: 10.1039/d2cp03497b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The impact of mutual interactions between the transmembrane domains of membrane proteins and lipids on bilayer properties has gained major attraction. Most simulation studies of membranes rely on the Martini force field, which has proven extremely helpful in providing molecular insights into realistic systems. Accordingly, an evaluation of the accuracy of the Martini force field is crucial to be able to correctly interpret the reported data. In this study, we combine atomistic and coarse-grained Martini simulations to investigate the properties of transmembrane domains (TMDs) in a model yeast membrane. The results show that the TMD binding state (monomeric and dimeric with positive or negative crossing angle) and the membrane composition significantly influence the properties around the TMDs and change TMD-TMD and TMD-lipid affinities. Furthermore, ergosterol (ERG) exhibits a strong affinity to TMD dimers. Importantly, the right-handed TMD dimer configuration is stabilized via TMD-TMD contacts by the addition of asymmetric anionic phosphatidylserine (PS). The coarse-grained simulations corroborate many of these findings, with two notable exceptions: a systematic overestimation of TMD-ERG interaction and lack of stabilization of the right-handed TMD dimers with the addition of PS.
Collapse
Affiliation(s)
- Fabian Keller
- Institut für Physikalische Chemie, Corrensstraße 28, Münster, Germany.
| | | | | | - Andreas Heuer
- Institut für Physikalische Chemie, Corrensstraße 28, Münster, Germany.
| |
Collapse
|
2
|
Joodaki F, Martin LM, Greenfield ML. Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:9481-9499. [PMID: 35901279 DOI: 10.1021/acs.langmuir.2c00483] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Å2. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.
Collapse
Affiliation(s)
- Faramarz Joodaki
- Department of Chemical Engineering, University of Rhode Island, 360 Fascitelli Center for Advanced Engineering, Kingston, Rhode Island 02881, United States
| | - Lenore M Martin
- Department of Cell and Molecular Biology, University of Rhode Island, 120 Flagg Road, Kingston, Rhode Island 02881, United States
| | - Michael L Greenfield
- Department of Chemical Engineering, University of Rhode Island, 360 Fascitelli Center for Advanced Engineering, Kingston, Rhode Island 02881, United States
| |
Collapse
|
3
|
Alavizargar A, Elting A, Wedlich-Söldner R, Heuer A. Lipid-Mediated Association of the Slg1 Transmembrane Domains in Yeast Plasma Membranes. J Phys Chem B 2022; 126:3240-3256. [PMID: 35446028 DOI: 10.1021/acs.jpcb.2c00192] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Clustering of transmembrane proteins underlies a multitude of fundamental biological processes at the plasma membrane (PM) such as receptor activation, lateral domain formation, and mechanotransduction. The self-association of the respective transmembrane domains (TMDs) has also been suggested to be responsible for the micron-scaled patterns seen for integral membrane proteins in the budding yeast PM. However, the underlying interplay between the local lipid composition and the TMD identity is still not mechanistically understood. In this work, we combined coarse-grained molecular dynamics simulations of simplified bilayer systems with high-resolution live-cell microscopy to analyze the distribution of a representative helical yeast TMD from the PM sensor Slg1 within different lipid environments. In our simulations, we specifically evaluated the effects of acyl chain saturation and anionic lipid head groups on the association of two TMDs. We found that weak lipid-protein interactions significantly affect the configuration of TMD dimers and the free energy of association. Increased amounts of unsaturated phospholipids (PLs) strongly reduced the helix-helix interaction, while the presence of anionic phosphatidylserine (PS) hardly affected the dimer formation. We could experimentally confirm this surprising lack of effect of PS using the network factor, a mesoscopic measure of PM pattern formation in yeast cells. Simulations also showed that the formation of TMD dimers in turn increased the order parameter of the surrounding lipids and induced long-range perturbations in lipid organization. In summary, our results shed new light on the mechanisms of lipid-mediated dimerization of TMDs in complex lipid mixtures.
Collapse
Affiliation(s)
- Azadeh Alavizargar
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| | - Annegret Elting
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Roland Wedlich-Söldner
- Institute of Cell Dynamics and Imaging, University of Muenster, Von-Esmarch-Str. 56, 48149 Muenster, Germany
| | - Andreas Heuer
- Institute of Physical Chemistry, University of Muenster, Corrensstr. 28/30, 48149 Muenster, Germany
| |
Collapse
|
4
|
Weigle AT, Carr M, Shukla D. Impact of Increased Membrane Realism on Conformational Sampling of Proteins. J Chem Theory Comput 2021; 17:5342-5357. [PMID: 34339605 DOI: 10.1021/acs.jctc.1c00276] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.
Collapse
Affiliation(s)
- Austin T Weigle
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Matthew Carr
- Independent Software Development Provider310 East Marlette Avenue, Phoenix, Arizona 85012, United States
| | - Diwakar Shukla
- Department of Chemical & Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Cancer Center at Illinois, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Center for Digital Agriculture, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,National Center for Supercomputing Applications, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,NIH Center for Macromolecular Modeling and Bioinformatics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
5
|
Wilson KA, Fairweather SJ, MacDermott-Opeskin HI, Wang L, Morris RA, O'Mara ML. The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane. J Chem Phys 2021; 154:095101. [PMID: 33685172 DOI: 10.1063/5.0040887] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.
Collapse
Affiliation(s)
- Katie A Wilson
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| | - Stephen J Fairweather
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| | - Hugo I MacDermott-Opeskin
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| | - Lily Wang
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| | - Richard A Morris
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| | - Megan L O'Mara
- Research School of Chemistry, College of Science, The Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
6
|
Kim M, Mun W, Jung WH, Lee J, Cho G, Kwon J, Ahn DJ, Mitchell RJ, Kim BS. Antimicrobial PEGtides: A Modular Poly(ethylene glycol)-Based Peptidomimetic Approach to Combat Bacteria. ACS NANO 2021; 15:9143-9153. [PMID: 33988968 DOI: 10.1021/acsnano.1c02644] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Despite their high potency, the widespread implementation of natural antimicrobial peptides is still challenging due to their low scalability and high hemolytic activities. Herein, we address these issues by employing a modular approach to mimic the key amino acid residues present in antimicrobial peptides, such as lysine, leucine, and serine, but on the highly biocompatible poly(ethylene glycol) (PEG) backbone. A series of these PEG-based peptides (PEGtides) were developed using functional epoxide monomers, corresponding to each key amino acid, with several possessing highly potent bactericidal activities and controlled selectivities, with respect to their hemolytic behavior. The critical role of the composition and the structure of the PEGtides in their selectivities was further supported by coarse-grained molecular dynamic simulations. This modular approach is anticipated to provide the design principles necessary for the future development of antimicrobial polymers.
Collapse
Affiliation(s)
- Minseong Kim
- Department of Chemistry, Yonsei University, Seoul 03722, Republic of Korea
| | | | | | - Joonhee Lee
- Department of Chemistry, Yonsei University, Seoul 03722, Republic of Korea
| | | | | | | | | | - Byeong-Su Kim
- Department of Chemistry, Yonsei University, Seoul 03722, Republic of Korea
| |
Collapse
|
7
|
Crowet JM, Buchoux S, Belloy N, Sarazin C, Lins L, Dauchez M. LIMONADA: A database dedicated to the simulation of biological membranes. J Comput Chem 2021; 42:1028-1033. [PMID: 33709443 DOI: 10.1002/jcc.26511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 02/16/2021] [Accepted: 02/19/2021] [Indexed: 11/06/2022]
Abstract
Cellular membranes are composed of a wide diversity of lipid species in varying proportions and these compositions are representative of the organism, cellular type and organelle to which they belong. Because models of these molecular systems simulated by MD steadily gain in size and complexity, they are increasingly representative of specific compositions and behaviors of biological membranes. Due to the number of lipid species involved, of force fields and topologies and because of the complexity of membrane objects that have been simulated, LIMONADA has been developed as an open database allowing to handle the various aspects of lipid membrane simulation. LIMONADA presents published membrane patches with their simulation files and the cellular membrane it models. Their compositions are then detailed based on the lipid identification from LIPID MAPS database plus the lipid topologies and the force field used. LIMONADA is freely accessible on the web at https://limonada.univ-reims.fr/.
Collapse
Affiliation(s)
- Jean-Marc Crowet
- Matrice Extracellulaire et Dynamique Cellulaire (UMR CNRS 7369), Chaire MAgICS, Université de Reims Champagne-Ardenne, Reims, France
| | - Sébastien Buchoux
- Unité de Génie Enzymatique et Cellulaire (GEC-UMR7025 CNRS/UPJV/UTC), Université de Picardie Jules Verne, Amiens, France
| | - Nicolas Belloy
- Matrice Extracellulaire et Dynamique Cellulaire (UMR CNRS 7369), Chaire MAgICS, Université de Reims Champagne-Ardenne, Reims, France
| | - Catherine Sarazin
- Unité de Génie Enzymatique et Cellulaire (GEC-UMR7025 CNRS/UPJV/UTC), Université de Picardie Jules Verne, Amiens, France
| | - Laurence Lins
- Laboratoire de Biophysique Moléculaire aux Interfaces, Université de Liège, Gembloux Agro-Bio Tech, Gembloux, Belgium
| | - Manuel Dauchez
- Matrice Extracellulaire et Dynamique Cellulaire (UMR CNRS 7369), Chaire MAgICS, Université de Reims Champagne-Ardenne, Reims, France
| |
Collapse
|
8
|
Jin Y, Liang Q, Tieleman DP. Interactions between Band 3 Anion Exchanger and Lipid Nanodomains in Ternary Lipid Bilayers: Atomistic Simulations. J Phys Chem B 2020; 124:3054-3064. [DOI: 10.1021/acs.jpcb.0c01055] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Yapan Jin
- Center for Statistical and Theoretical Condensed Matter Physics and Department of Physics, Zhejiang Normal University, Jinhua 321004, P. R. China
| | - Qing Liang
- Center for Statistical and Theoretical Condensed Matter Physics and Department of Physics, Zhejiang Normal University, Jinhua 321004, P. R. China
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Dr. NW, Calgary, AB T2N 1N4, Canada
| |
Collapse
|
9
|
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: 435] [Impact Index Per Article: 87.0] [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
|
10
|
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: 184] [Impact Index Per Article: 36.8] [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
|
11
|
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
|
12
|
Klauda JB. Perspective: Computational modeling of accurate cellular membranes with molecular resolution. J Chem Phys 2018; 149:220901. [DOI: 10.1063/1.5055007] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
13
|
Goossens K, De Winter H. Molecular Dynamics Simulations of Membrane Proteins: An Overview. J Chem Inf Model 2018; 58:2193-2202. [PMID: 30336018 DOI: 10.1021/acs.jcim.8b00639] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.
Collapse
Affiliation(s)
- Kenneth Goossens
- Department of Pharmaceutical Sciences, Laboratory of Medicinal Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Wilrijk , Belgium
| | - Hans De Winter
- Department of Pharmaceutical Sciences, Laboratory of Medicinal Chemistry , University of Antwerp , Universiteitsplein 1 , 2610 Wilrijk , Belgium
| |
Collapse
|
14
|
Kalli AC, Reithmeier RAF. Interaction of the human erythrocyte Band 3 anion exchanger 1 (AE1, SLC4A1) with lipids and glycophorin A: Molecular organization of the Wright (Wr) blood group antigen. PLoS Comput Biol 2018; 14:e1006284. [PMID: 30011272 PMCID: PMC6080803 DOI: 10.1371/journal.pcbi.1006284] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Revised: 08/07/2018] [Accepted: 06/08/2018] [Indexed: 11/29/2022] Open
Abstract
The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5–bisphosphate (PIP2) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP2 molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers. Human Band 3 (AE1, SLC4A1), an abundant 911 amino acid glycoprotein, catalyzes the exchange of bicarbonate and chloride across the red blood cell membrane, a process necessary for efficient respiration. Malfunction of Band 3 leads to inherited diseases such as Southeast Asian Ovalocytosis, hereditary spherocytosis and distal renal tubular acidosis. Despite much available structural and functional data about Band 3, key questions about the conformational changes associated with transport and the molecular details of its interaction with lipids and other proteins remain unanswered. In this study, we have used computer simulations to investigate the dynamics of Band 3 in lipid bilayers that resemble the red blood cell plasma membrane. Our results suggest that negatively charged phospholipids and cholesterol interact strongly with Band 3 forming an annulus around the protein. Glycophorin A (GPA) interacts with Band 3 to form the Wright (Wr) blood group antigen. We were able to recreate this complex and show that GPA promotes the clustering of Band 3 in red blood cell membranes. Understanding the molecular details of the interaction of Band 3 with GPA has provided new insights into the nature of the Wright blood group antigen.
Collapse
Affiliation(s)
- Antreas C. Kalli
- Leeds Institute of Cancer and Pathology, University of Leeds, Leeds, United Kingdom
- Astbury Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom
- * E-mail:
| | | |
Collapse
|
15
|
Corradi V, Mendez-Villuendas E, Ingólfsson HI, Gu RX, Siuda I, Melo MN, Moussatova A, DeGagné LJ, Sejdiu BI, Singh G, Wassenaar TA, Delgado Magnero K, Marrink SJ, Tieleman DP. Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins. ACS CENTRAL SCIENCE 2018; 4:709-717. [PMID: 29974066 PMCID: PMC6028153 DOI: 10.1021/acscentsci.8b00143] [Citation(s) in RCA: 212] [Impact Index Per Article: 35.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Indexed: 05/08/2023]
Abstract
Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.
Collapse
Affiliation(s)
- Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Eduardo Mendez-Villuendas
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Helgi I. Ingólfsson
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Ruo-Xu Gu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Iwona Siuda
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Manuel N. Melo
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Anastassiia Moussatova
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Lucien J. DeGagné
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Besian I. Sejdiu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Gurpreet Singh
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Tsjerk A. Wassenaar
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Karelia Delgado Magnero
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
- E-mail:
| |
Collapse
|
16
|
Ermakova E, Zuev Y. Interaction of Scots Pine Defensin with Model Membrane by Coarse-Grained Molecular Dynamics. J Membr Biol 2017; 250:205-216. [DOI: 10.1007/s00232-017-9950-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2016] [Accepted: 02/06/2017] [Indexed: 01/23/2023]
|
17
|
Bocharov EV, Mineev KS, Pavlov KV, Akimov SA, Kuznetsov AS, Efremov RG, Arseniev AS. Helix-helix interactions in membrane domains of bitopic proteins: Specificity and role of lipid environment. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1859:561-576. [PMID: 27884807 DOI: 10.1016/j.bbamem.2016.10.024] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Revised: 09/18/2016] [Accepted: 10/20/2016] [Indexed: 12/23/2022]
Abstract
Interaction between transmembrane helices often determines biological activity of membrane proteins. Bitopic proteins, a broad subclass of membrane proteins, form dimers containing two membrane-spanning helices. Some aspects of their structure-function relationship cannot be fully understood without considering the protein-lipid interaction, which can determine the protein conformational ensemble. Experimental and computer modeling data concerning transmembrane parts of bitopic proteins are reviewed in the present paper. They highlight the importance of lipid-protein interactions and resolve certain paradoxes in the behavior of such proteins. Besides, some properties of membrane organization provided a clue to understanding of allosteric interactions between distant parts of proteins. Interactions of these kinds appear to underlie a signaling mechanism, which could be widely employed in the functioning of many membrane proteins. Treatment of membrane proteins as parts of integrated fine-tuned proteolipid system promises new insights into biological function mechanisms and approaches to drug design. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.
Collapse
Affiliation(s)
- Eduard V Bocharov
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, Miklukho-Maklaya ul. 16/10, Moscow, 117997, Russian Federation; National Research Centre "Kurchatov Institute", Akad. Kurchatova pl. 1, Moscow, 123182, Russian Federation.
| | - Konstantin S Mineev
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, Miklukho-Maklaya ul. 16/10, Moscow, 117997, Russian Federation
| | - Konstantin V Pavlov
- Frumkin Institute of Physical Chemistry and Electrochemistry RAS, Leninskiy prospect 31/5, Moscow, 119071, Russian Federation
| | - Sergey A Akimov
- Frumkin Institute of Physical Chemistry and Electrochemistry RAS, Leninskiy prospect 31/5, Moscow, 119071, Russian Federation; National University of Science and Technology "MISiS", Leninskiy prospect 4, Moscow, 119049, Russian Federation
| | - Andrey S Kuznetsov
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, Miklukho-Maklaya ul. 16/10, Moscow, 117997, Russian Federation
| | - Roman G Efremov
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, Miklukho-Maklaya ul. 16/10, Moscow, 117997, Russian Federation; Higher School of Economics, Myasnitskaya ul. 20, Moscow, 101000, Russian Federation
| | - Alexander S Arseniev
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry RAS, Miklukho-Maklaya ul. 16/10, Moscow, 117997, Russian Federation.
| |
Collapse
|
18
|
Huang JY, Lin CY. Exploring the stochastic dynamics of correlated movement of receptor proteins in plasma membranes in vivo. J Chem Phys 2015; 143:225101. [PMID: 26671403 DOI: 10.1063/1.4936963] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ligand-induced receptor dimerization plays a crucial role in the signaling process of living cells. In this study, we developed a theoretical model and performed single-molecule tracking to explore the correlated diffusion processes of liganded epidermal growth factor receptors prior to dimer formation. We disclosed that both an attractive potential between liganded receptor proteins in proximity and correlated fluctuations in the local environments of the proteins play an important role to produce the observed correlated movement of the receptors. This result can serve as the foundation to shed light on the way in which receptor functions are regulated in plasma membranes in vivo.
Collapse
Affiliation(s)
- Jung Y Huang
- The T.K.B. Research Center of Photonics, Chiao Tung University, Hsinchu 300, Taiwan
| | - Chien Y Lin
- Department of Photonics, Chiao Tung University, Hsinchu 300, Taiwan
| |
Collapse
|
19
|
Sharma S, Kim BN, Stansfeld PJ, Sansom MSP, Lindau M. A Coarse Grained Model for a Lipid Membrane with Physiological Composition and Leaflet Asymmetry. PLoS One 2015; 10:e0144814. [PMID: 26659855 PMCID: PMC4681583 DOI: 10.1371/journal.pone.0144814] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2015] [Accepted: 11/24/2015] [Indexed: 11/24/2022] Open
Abstract
The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.
Collapse
Affiliation(s)
- Satyan Sharma
- Laboratory for Nanoscale Cell Biology, Max-Planck-Institute for Biophysical Chemistry, Göttingen, Germany
- * E-mail:
| | - Brian N. Kim
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York, United States of America
| | - Phillip J. Stansfeld
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, England, United Kingdom
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, England, United Kingdom
| | - Manfred Lindau
- Laboratory for Nanoscale Cell Biology, Max-Planck-Institute for Biophysical Chemistry, Göttingen, Germany
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York, United States of America
| |
Collapse
|