1
|
Levintov L, Gorai B, Vashisth H. Spontaneous Dimerization and Distinct Packing Modes of Transmembrane Domains in Receptor Tyrosine Kinases. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.09.593448. [PMID: 38798363 PMCID: PMC11118388 DOI: 10.1101/2024.05.09.593448] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
The insulin receptor (IR) and the insulin-like growth factor-1 receptor (IGF1R) are homodimeric transmembrane glycoproteins that transduce signals across the membrane on binding of extracellular peptide ligands. The structures of IR/IGF1R fragments in apo and liganded states have revealed that the extracellular subunits of these receptors adopt Λ-shaped configurations to which are connected the intracellular tyrosine kinase (TK) domains. The binding of peptide ligands induces structural transitions in the extracellular subunits leading to potential dimerization of transmembrane domains (TMDs) and autophosphorylation in TKs. However, the activation mechanisms of IR/IGF1R, especially the role of TMDs in coordinating signal-inducing structural transitions, remain poorly understood, in part due to the lack of structures of full-length receptors in apo or liganded states. While atomistic simulations of IR/IGF1R TMDs showed that these domains can dimerize in single component membranes, spontaneous unbiased dimerization in a plasma membrane having physiologically representative lipid composition has not been observed. We address this limitation by employing coarse-grained (CG) molecular dynamics simulations to probe the dimerization propensity of IR/IGF1R TMDs. We observed that TMDs in both receptors spontaneously dimerized independent of their initial orientations in their dissociated states, signifying their natural propensity for dimerization. In the dimeric state, IR TMDs predominantly adopted X-shaped configurations with asymmetric helical packing and significant tilt relative to the membrane normal, while IGF1R TMDs adopted symmetric V-shaped or parallel configurations with either no tilt or a small tilt relative to the membrane normal. Our results suggest that IR/IGF1R TMDs spontaneously dimerize and adopt distinct dimerized configurations.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering and Bioengineering, University of New Hampshire, Durham 03824, New Hampshire, USA
| | - Biswajit Gorai
- Institute of Chemistry, Technical University of Berlin, Berlin 10623, Germany
| | - Harish Vashisth
- Department of Chemical Engineering and Bioengineering, University of New Hampshire, Durham 03824, New Hampshire, USA
- Department of Chemistry, University of New Hampshire, Durham 03824, New Hampshire, USA
- Integrated Applied Mathematics Program, University of New Hampshire, Durham 03824, New Hampshire, USA
- Molecular and Cellular Biotechnology Program, University of New Hampshire, Durham 03824, New Hampshire, USA
| |
Collapse
|
2
|
Ito S, Sugita Y. Free-energy landscapes of transmembrane homodimers by bias-exchange adaptively biased molecular dynamics. Biophys Chem 2024; 307:107190. [PMID: 38290241 DOI: 10.1016/j.bpc.2024.107190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2023] [Revised: 01/21/2024] [Accepted: 01/21/2024] [Indexed: 02/01/2024]
Abstract
Membrane proteins play essential roles in various biological functions within the cell. One of the most common functional regulations involves the dimerization of two single-pass transmembrane (TM) helices. Glycophorin A (GpA) and amyloid precursor protein (APP) form TM homodimers in the membrane, which have been investigated both experimentally and computationally. The homodimer structures are well characterized using only four collective variables (CVs) when each TM helix is stable. The CVs are the interhelical distance, the crossing angle, and the Crick angles for two TM helices. However, conformational sampling with multi-dimensional replica-exchange umbrella sampling (REUS) requires too many replicas to sample all the CVs for exploring the conformational landscapes. Here, we show that the bias-exchange adaptively biased molecular dynamics (BE-ABMD) with the four CVs effectively explores the free-energy landscapes of the TM helix dimers of GpA, wild-type APP and its mutants in the IMM1 implicit membrane. Compared to the original ABMD, the bias-exchange algorithm in BE-ABMD can provide a more rapidly converged conformational landscape. The BE-ABMD simulations could also reveal TM packing interfaces of the membrane proteins and the dependence of the free-energy landscapes on the membrane thickness. This approach is valuable for numerous other applications, including those involving explicit solvent and a lipid bilayer in all-atom force fields or Martini coarse-grained models, and enhances our understanding of protein-protein interactions in biological membranes.
Collapse
Affiliation(s)
- Shingo Ito
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| |
Collapse
|
3
|
Chen SY, Feilen LP, Chávez-Gutiérrez L, Steiner H, Zacharias M. Enzyme-substrate hybrid β-sheet controls geometry and water access to the γ-secretase active site. Commun Biol 2023; 6:670. [PMID: 37355752 PMCID: PMC10290658 DOI: 10.1038/s42003-023-05039-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 06/12/2023] [Indexed: 06/26/2023] Open
Abstract
γ-Secretase is an aspartyl intramembrane protease that cleaves the amyloid precursor protein (APP) involved in Alzheimer's disease pathology and other transmembrane proteins. Substrate-bound structures reveal a stable hybrid β-sheet immediately following the substrate scissile bond consisting of β1 and β2 from the enzyme and β3 from the substrate. Molecular dynamics simulations and enhanced sampling simulations demonstrate that the hybrid β-sheet stability is strongly correlated with the formation of a stable cleavage-compatible active geometry and it also controls water access to the active site. The hybrid β-sheet is only stable for substrates with 3 or more C-terminal residues beyond the scissile bond. The simulation model allowed us to predict the effect of Pro and Phe mutations that weaken the formation of the hybrid β-sheet which were confirmed by experimental testing. Our study provides a direct explanation why γ-secretase preferentially cleaves APP in steps of 3 residues and how the hybrid β-sheet facilitates γ-secretase proteolysis.
Collapse
Affiliation(s)
- Shu-Yu Chen
- Center of Functional Protein Assemblies, Technical University of Munich, Garching, Germany
| | - Lukas P Feilen
- German Center for Neurodegenerative Diseases (DZNE), Munich, Germany
| | - Lucía Chávez-Gutiérrez
- VIB-KU Leuven Center for Brain & Disease Research, Leuven, Belgium
- Department of Neurosciences, Leuven Research Institute for Neuroscience and Disease (LIND), KU Leuven, Leuven, Belgium
| | - Harald Steiner
- German Center for Neurodegenerative Diseases (DZNE), Munich, Germany
- Biomedical Center (BMC), Division of Metabolic Biochemistry, Faculty of Medicine, LMU Munich, Germany
| | - Martin Zacharias
- Center of Functional Protein Assemblies, Technical University of Munich, Garching, Germany.
| |
Collapse
|
4
|
Zang Y, Wang H, Hao D, Kang Y, Zhang J, Li X, Zhang L, Yang Z, Zhang S. p38α Kinase Auto-Activation through Its Conformational Transition Induced by Tyr323 Phosphorylation. J Chem Inf Model 2022; 62:6639-6648. [PMID: 36394912 DOI: 10.1021/acs.jcim.2c00236] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
p38α is a key serine/threonine kinase that can enable atypical auto-activation through Zap70 phosphorylation and initiate T cell receptor signaling. The auto-activation plays an important role in autoimmune diseases. Although the classical activation mechanism of p38α has been studied in-depth, the atypical activation mechanism of Y323 phosphorylation-induced p38α auto-activation remains largely unexplained, especially the regulatory effects of phosphorylation on different sites (Y323 vs T180). From the X-ray experimental data, we identified the inactive and active states of p38α using principal component analysis. To understand the auto-activation process and the internal driving mechanism, a computational paradigm that couples the targeted molecular dynamics simulations, the String Method, and the umbrella sampling strategy were employed to generate the conformational landscape of p38α, including p38α T180-Y323, p38α T180-pY323, and p38α pT180-pY323 systems (pT180/pY323: phosphorylated T180/Y323). We explored that pY323 could change the conformational distribution and promote the conformational transition of p38α from the inactive state to the active state. Auto-activation of p38α is regulated by pY323 through destabilization of the hydrophobic core structure and aided by R173. This study will further explain the conformational transition of p38α induced by Y323 phosphorylation and provide insights into the universal molecular auto-activation mechanism of the p38 subfamily at the atomic level.
Collapse
Affiliation(s)
- Yongjian Zang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - He Wang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Dongxiao Hao
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Ying Kang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Jianwen Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Lei Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| | - Shengli Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an710049, China
| |
Collapse
|
5
|
Ma Z, Shi S, Ren M, Pang C, Zhan Y, An H, Sun F. Molecular Mechanism of CD44 Homodimerization Modulated by Palmitoylation and Membrane Environments. Biophys J 2022; 121:2671-2683. [PMID: 35733341 DOI: 10.1016/j.bpj.2022.06.021] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 04/21/2022] [Accepted: 06/16/2022] [Indexed: 11/15/2022] Open
Abstract
The homodimerization of CD44 plays a key role in an intercellular-to-extracellular signal transduction and tumor progression. Acylated modification and specific membrane environments have been reported to mediate translocation and oligomerization of CD44, however, the underlying molecular mechanism remains elusive. In this study, extensive molecular dynamics simulations are performed to characterize the dimerization of palmitoylated CD44 variants in different bilayer environments. CD44 forms homodimer depending on the cysteines on the juxta-membrane domains, and the dimerization efficiency and packing configurations are defected by their palmitoylated modifications. In the phase-segregated (raft included) membrane, homodimerization of the palmitoylated CD44 is hardly observed, whereas PIP2 addition compensates to realize dimerization. However, the structure of CD44 homodimer formed in the phase-segregated bilayer turns susceptive and PIP2 addition allows for an extensive conformation of the cytoplasmic domain, a proposal prerequisite to access the cytoskeleton linker proteins. The results unravel a delicate competitive relationship between PIP2, lipid raft and palmitoylation in mediating protein homodimerization, which helps to clarify the dynamic dimer conformations and involved cellular signaling of the CD44 likewise proteins.
Collapse
Affiliation(s)
- Ziyi Ma
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China
| | - Sai Shi
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China; State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300401, China; Key Laboratory of Electromagnetic Field and Electrical Apparatus Reliability of Hebei Province, Hebei University of Technology, Tianjin 300401, China
| | - Meina Ren
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China
| | - Chunli Pang
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China
| | - Yong Zhan
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China; State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300401, China; Key Laboratory of Electromagnetic Field and Electrical Apparatus Reliability of Hebei Province, Hebei University of Technology, Tianjin 300401, China
| | - Hailong An
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China; State Key Laboratory of Reliability and Intelligence of Electrical Equipment, Hebei University of Technology, Tianjin 300401, China; Key Laboratory of Electromagnetic Field and Electrical Apparatus Reliability of Hebei Province, Hebei University of Technology, Tianjin 300401, China.
| | - Fude Sun
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin, 300401, China.
| |
Collapse
|
6
|
A Free-Energy Landscape Analysis of Calmodulin Obtained from an NMR Data-Utilized Multi-Scale Divide-and-Conquer Molecular Dynamics Simulation. Life (Basel) 2021; 11:life11111241. [PMID: 34833117 PMCID: PMC8617919 DOI: 10.3390/life11111241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 11/08/2021] [Accepted: 11/10/2021] [Indexed: 11/16/2022] Open
Abstract
Calmodulin (CaM) is a multifunctional calcium-binding protein, which regulates a variety of biochemical processes. CaM acts through its conformational changes and complex formation with its target enzymes. CaM consists of two globular domains (N-lobe and C-lobe) linked by an extended linker region. Upon calcium binding, the N-lobe and C-lobe undergo local conformational changes, followed by a major conformational change of the entire CaM to wrap the target enzyme. However, the regulation mechanisms, such as allosteric interactions, which regulate the large structural changes, are still unclear. In order to investigate the series of structural changes, the free-energy landscape of CaM was obtained by multi-scale divide-and-conquer molecular dynamics (MSDC-MD). The resultant free-energy landscape (FEL) shows that the Ca2+ bound CaM (holo-CaM) would take an experimentally famous elongated structure, which can be formed in the early stage of structural change, by breaking the inter-domain interactions. The FEL also shows that important interactions complete the structural change from the elongated structure to the ring-like structure. In addition, the FEL might give a guiding principle to predict mutational sites in CaM. In this study, it was demonstrated that the movement process of macroscopic variables on the FEL may be diffusive to some extent, and then, the MSDC-MD is suitable to the parallel computation.
Collapse
|
7
|
Ozkan A, Sitharam M, Flores-Canales JC, Prabhu R, Kurnikova M. Baseline Comparisons of Complementary Sampling Methods for Assembly Driven by Short-Ranged Pair Potentials toward Fast and Flexible Hybridization. J Chem Theory Comput 2021; 17:1967-1987. [PMID: 33576635 DOI: 10.1021/acs.jctc.0c00945] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This work measures baseline sampling characteristics that highlight fundamental differences between sampling methods for assembly driven by short-ranged pair potentials. Such granular comparison is essential for fast, flexible, and accurate hybridization of complementary methods. Besides sampling speed, efficiency, and accuracy of uniform grid coverage, other sampling characteristics measured are (i) accuracy of covering narrow low energy regions that have low effective dimension (ii) ability to localize sampling to specific basins, and (iii) flexibility in sampling distributions. As a proof of concept, we compare a recently developed geometric methodology EASAL (Efficient Atlasing and Search of Assembly Landscapes) and the traditional Monte Carlo (MC) method for sampling the energy landscape of two assembling trans-membrane helices, driven by short-range pair potentials. By measuring the above-mentioned sampling characteristics, we demonstrate that EASAL provides localized and accurate coverage of crucial regions of the energy landscape of low effective dimension, under flexible sampling distributions, with much fewer samples and computational resources than MC sampling. EASAL's empirically validated theoretical guarantees permit credible extrapolation of these measurements and comparisons to arbitrary number and size of assembling units. Promising avenues for hybridizing the complementary advantages of the two methods are discussed.
Collapse
Affiliation(s)
- Aysegul Ozkan
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | - Meera Sitharam
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | | | - Rahul Prabhu
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | - Maria Kurnikova
- Chemistry Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
8
|
Ramesh S, Park S, Call MJ, Im W, Call ME. Experimentally Guided Computational Methods Yield Highly Accurate Insights into Transmembrane Interactions within the T Cell Receptor Complex. J Phys Chem B 2020; 124:10303-10310. [PMID: 33030343 DOI: 10.1021/acs.jpcb.0c06403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding how molecular interactions within the plasma membrane govern assembly, clustering, and conformational changes in single-pass transmembrane (TM) receptors has long presented substantial experimental challenges. Our previous work on activating immune receptors has combined direct biochemical and biophysical characterizations with both independent and experimentally restrained computational methods to provide novel insights into the key TM interactions underpinning assembly and stability of complex, multisubunit receptor systems. The recently published cryo-EM structure of the intact T cell receptor (TCR)-CD3 complex provides a unique opportunity to test the models and predictions arising from these studies, and we find that they are accurate, which we attribute to robust simulation environments and careful consideration of limitations related to studying TM interactions in isolation from additional receptor domains. With this in mind, we revisit results in other immune receptors and look forward to how similar methods may be applied to understand receptors for which little or no structural information is currently available.
Collapse
Affiliation(s)
- Samyuktha Ramesh
- The Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria 3052, Australia.,Department of Medical Biology, University of Melbourne, Parkville, Victoria 3052, Australia
| | - Soohyung Park
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Melissa J Call
- The Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria 3052, Australia.,Department of Medical Biology, University of Melbourne, Parkville, Victoria 3052, Australia
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Matthew E Call
- The Walter and Eliza Hall Institute of Medical Research, Parkville, Victoria 3052, Australia.,Department of Medical Biology, University of Melbourne, Parkville, Victoria 3052, Australia
| |
Collapse
|
9
|
Martinotti C, Ruiz-Perez L, Deplazes E, Mancera RL. Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes. Chemphyschem 2020; 21:1486-1514. [PMID: 32452115 DOI: 10.1002/cphc.202000219] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Revised: 05/22/2020] [Indexed: 12/12/2022]
Abstract
Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.
Collapse
Affiliation(s)
- Carlo Martinotti
- School of Pharmacy and Biomedical Sciences, Curtin Health Innovation Research Institute and, Curtin Institute for Computation, Curtin University, Perth, WA 6845, Australia
| | - Lanie Ruiz-Perez
- School of Pharmacy and Biomedical Sciences, Curtin Health Innovation Research Institute and, Curtin Institute for Computation, Curtin University, Perth, WA 6845, Australia
| | - Evelyne Deplazes
- School of Life Sciences, University of Technology Sydney, Sydney, NSW 2007, Australia
| | - Ricardo L Mancera
- School of Pharmacy and Biomedical Sciences, Curtin Health Innovation Research Institute and, Curtin Institute for Computation, Curtin University, Perth, WA 6845, Australia
| |
Collapse
|
10
|
Shirts MR, Ferguson AL. Statistically Optimal Continuous Free Energy Surfaces from Biased Simulations and Multistate Reweighting. J Chem Theory Comput 2020; 16:4107-4125. [DOI: 10.1021/acs.jctc.0c00077] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Michael R. Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Andrew L. Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
11
|
Corey RA, Stansfeld PJ, Sansom MS. The energetics of protein-lipid interactions as viewed by molecular simulations. Biochem Soc Trans 2020; 48:25-37. [PMID: 31872229 PMCID: PMC7054751 DOI: 10.1042/bst20190149] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Revised: 12/06/2019] [Accepted: 12/10/2019] [Indexed: 12/20/2022]
Abstract
Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.
Collapse
Affiliation(s)
- Robin A. Corey
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| | - Phillip J. Stansfeld
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
- School of Life Sciences and Department of Chemistry, University of Warwick, Coventry CV4 7AL, U.K
| | - Mark S.P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K
| |
Collapse
|
12
|
Matsuoka D, Kamiya M, Sato T, Sugita Y. Role of the N-Terminal Transmembrane Helix Contacts in the Activation of FGFR3. J Comput Chem 2019; 41:561-572. [PMID: 31804721 DOI: 10.1002/jcc.26122] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2019] [Revised: 10/31/2019] [Accepted: 11/08/2019] [Indexed: 12/16/2022]
Abstract
Fibroblast growth factor receptor 3 (FGFR3) is a member of receptor tyrosine kinases, which is involved in skeletal cell growth, differentiation, and migration. FGFR3 transduces biochemical signals from the extracellular ligand-binding domain to the intracellular kinase domain through the conformational changes of the transmembrane (TM) helix dimer. Here, we apply generalized replica exchange with solute tempering method to wild type (WT) and G380R mutant (G380R) of FGFR3. The dimer interface in G380R is different from WT and the simulation results are in good agreement with the solid-state nuclear magnetic resonance (NMR) spectroscopy. TM helices in G380R are extended more than WT, and thereby, G375 in G380R contacts near the N-termini of the TM helix dimer. Considering that both G380R and G375C show the constitutive activation, the formation of the N-terminal contacts of the TM helices can be generally important for the activation mechanism. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Daisuke Matsuoka
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Saitama, 351-0198, Japan
| | - Motoshi Kamiya
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, 650-0047, Japan
| | - Takeshi Sato
- Division of Liberal Arts and Science, Kyoto Pharmaceutical University, Kyoto, 607-8414, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako-shi, Saitama, 351-0198, Japan.,Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, 650-0047, Japan.,Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystem Dynamics Research, Kobe, 650-0047, Japan
| |
Collapse
|
13
|
Mori T, Sugita Y. Implicit Micelle Model for Membrane Proteins Using Superellipsoid Approximation. J Chem Theory Comput 2019; 16:711-724. [DOI: 10.1021/acs.jctc.9b00783] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama, 351-0198, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako-shi, Saitama, 351-0198, Japan
- RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- RIKEN Center for Biosystems Dynamics Research, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
14
|
Balusek C, Hwang H, Lau CH, Lundquist K, Hazel A, Pavlova A, Lynch DL, Reggio PH, Wang Y, Gumbart JC. Accelerating Membrane Simulations with Hydrogen Mass Repartitioning. J Chem Theory Comput 2019; 15:4673-4686. [PMID: 31265271 PMCID: PMC7271963 DOI: 10.1021/acs.jctc.9b00160] [Citation(s) in RCA: 77] [Impact Index Per Article: 15.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.
Collapse
Affiliation(s)
| | | | - Chun Hon Lau
- Department of Physics , The Chinese University of Hong Kong , Shatin, NT, Hong Kong , People's Republic of China
| | | | | | | | - Diane L Lynch
- Department of Chemistry and Biochemistry , University of North Carolina , Greensboro , North Carolina 27402 , United States
| | - Patricia H Reggio
- Department of Chemistry and Biochemistry , University of North Carolina , Greensboro , North Carolina 27402 , United States
| | - Yi Wang
- Department of Physics , The Chinese University of Hong Kong , Shatin, NT, Hong Kong , People's Republic of China
| | | |
Collapse
|
15
|
Minh Hung H, Dieu Hang T, Nguyen MT. Structural Investigation of Human Prolactin Receptor Transmembrane Domain Homodimerization in a Membrane Environment through Multiscale Simulations. J Phys Chem B 2019; 123:4858-4866. [PMID: 31099581 DOI: 10.1021/acs.jpcb.9b01986] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
It is well established that prolactin (PRL) and its receptor (PRLR) are associated with hundreds of biological functions. They have been postulated to be linked to breast and prostate cancers, and PRLR signaling has attracted considerable medical and pharmaceutical interest in the development of compounds targeting PRLR. Dimerization of the receptor through its transmembrane (TM) domain is a key step for understanding its signaling and related issues. Our multiscale simulation results revealed that its TM domain can form dimers in a membrane environment with distinct states stabilized by different residue motifs. On the basis of the simulated data, an activation mechanism of PRL with the importance of two symmetrical tryptophan residues was proposed in detail to determine the conformational change of its receptor, which is essential for signal transduction. The better knowledge of PRLR structure and its protein-protein interaction can considerably contribute to a further understanding of PRLR signaling action and thereby help to develop some new PRLR signaling-based strategies for PRL-related diseases.
Collapse
Affiliation(s)
- Huynh Minh Hung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium.,Department of Chemistry , Quy Nhon University , Quy Nhon 590000 , Vietnam
| | - Tran Dieu Hang
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium.,Department of Chemistry , Quy Nhon University , Quy Nhon 590000 , Vietnam
| | - Minh Tho Nguyen
- Computational Chemistry Research Group , Ton Duc Thang University , Ho Chi Minh City 700000 Vietnam.,Faculty of Applied Sciences , Ton Duc Thang University , Ho Chi Minh City 700000 Vietnam
| |
Collapse
|
16
|
Domański J, Sansom MSP, Stansfeld PJ, Best RB. Balancing Force Field Protein-Lipid Interactions To Capture Transmembrane Helix-Helix Association. J Chem Theory Comput 2018; 14:1706-1715. [PMID: 29424543 PMCID: PMC5852462 DOI: 10.1021/acs.jctc.7b00983] [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: 09/21/2017] [Indexed: 01/21/2023]
Abstract
Atomistic simulations have recently been shown to be sufficiently accurate to reversibly fold globular proteins and have provided insights into folding mechanisms. Gaining similar understanding from simulations of membrane protein folding and association would be of great medical interest. All-atom simulations of the folding and assembly of transmembrane protein domains are much more challenging, not least due to very slow diffusion within the lipid bilayer membrane. Here, we focus on a simple and well-characterized prototype of membrane protein folding and assembly, namely the dimerization of glycophorin A, a homodimer of single transmembrane helices. We have determined the free energy landscape for association of the dimer using the CHARMM36 force field. We find that the native structure is a metastable state, but not stable as expected from experimental estimates of the dissociation constant and numerous experimental structures obtained under a variety of conditions. We explore two straightforward approaches to address this problem and demonstrate that they result in stable dimers with dissociation constants consistent with experimental data.
Collapse
Affiliation(s)
- Jan Domański
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
- Laboratory
of Chemical Physics, National Institute of Diabetes and Digestive
and Kidney Diseases, National Institutes
of Health, Bethesda, Maryland 20892-0520, United States
| | - Mark S. P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| | - Phillip J. Stansfeld
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| | - Robert B. Best
- Laboratory
of Chemical Physics, National Institute of Diabetes and Digestive
and Kidney Diseases, National Institutes
of Health, Bethesda, Maryland 20892-0520, United States
| |
Collapse
|
17
|
Pawar AB, Sengupta D. Effect of Membrane Composition on Receptor Association: Implications of Cancer Lipidomics on ErbB Receptors. J Membr Biol 2018; 251:359-368. [DOI: 10.1007/s00232-018-0015-1] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2017] [Accepted: 01/15/2018] [Indexed: 10/18/2022]
|
18
|
Hanke CA, Gohlke H. Tertiary Interactions in the Unbound Guanine-Sensing Riboswitch Focus Functional Conformational Variability on the Binding Site. J Chem Inf Model 2017; 57:2822-2832. [DOI: 10.1021/acs.jcim.7b00567] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Christian A. Hanke
- Mathematisch-Naturwissenschaftliche
Fakultät, Institut für Pharmazeutische und Medizinische
Chemie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
| | - Holger Gohlke
- Mathematisch-Naturwissenschaftliche
Fakultät, Institut für Pharmazeutische und Medizinische
Chemie, Heinrich-Heine-Universität Düsseldorf, 40225 Düsseldorf, Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC) & Institute for Complex Systems - Structural Biochemistry (ICS 6), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| |
Collapse
|
19
|
Kobayashi C, Jung J, Matsunaga Y, Mori T, Ando T, Tamura K, Kamiya M, Sugita Y. GENESIS 1.1: A hybrid-parallel molecular dynamics simulator with enhanced sampling algorithms on multiple computational platforms. J Comput Chem 2017; 38:2193-2206. [PMID: 28718930 DOI: 10.1002/jcc.24874] [Citation(s) in RCA: 111] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Revised: 06/08/2017] [Accepted: 06/09/2017] [Indexed: 01/09/2023]
Abstract
GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Yasuhiro Matsunaga
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,JST PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Tadashi Ando
- Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center Computational Biology Research Core, 1-6-5 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Department of Applied Electronics, Faculty of Industrial Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan.,Water Frontier Science and Technology Research Center, Research Institute for Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan.,Research Division of Multiscale Interfacial Thermofluid Dynamics, Research Institute for Science and Technology, Tokyo University of Science, 6-3-1 Niijuku, Katsushika-ku, Tokyo, 125-8585, Japan
| | - Koichi Tamura
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Motoshi Kamiya
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan.,Theoretical Molecular Science Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama, 351-0198, Japan.,Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center Computational Biology Research Core, 1-6-5 Minatojima-minamachi, Chuo-ku, Kobe, 650-0047, Japan
| |
Collapse
|
20
|
Conformational Heterogeneity of the HIV Envelope Glycan Shield. Sci Rep 2017; 7:4435. [PMID: 28667249 PMCID: PMC5493700 DOI: 10.1038/s41598-017-04532-9] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2017] [Accepted: 05/16/2017] [Indexed: 01/08/2023] Open
Abstract
To better understand the conformational properties of the glycan shield covering the surface of the HIV gp120/gp41 envelope (Env) trimer, and how the glycan shield impacts the accessibility of the underlying protein surface, we performed enhanced sampling molecular dynamics (MD) simulations of a model glycosylated HIV Env protein and related systems. Our simulation studies revealed a conformationally heterogeneous glycan shield with a network of glycan-glycan interactions more extensive than those observed to date. We found that partial preorganization of the glycans potentially favors binding by established broadly neutralizing antibodies; omission of several specific glycans could increase the accessibility of other glycans or regions of the protein surface to antibody or CD4 receptor binding; the number of glycans that can potentially interact with known antibodies is larger than that observed in experimental studies; and specific glycan conformations can maximize or minimize interactions with individual antibodies. More broadly, the enhanced sampling MD simulations described here provide a valuable tool to guide the engineering of specific Env glycoforms for HIV vaccine design.
Collapse
|
21
|
Computational Studies of Molecular Permeation through Connexin26 Channels. Biophys J 2017; 110:584-599. [PMID: 26840724 DOI: 10.1016/j.bpj.2015.11.3528] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2015] [Revised: 10/21/2015] [Accepted: 11/23/2015] [Indexed: 11/22/2022] Open
Abstract
A signal property of connexin channels is the ability to mediate selective diffusive movement of molecules through plasma membrane(s), but the energetics and determinants of molecular movement through these channels have yet to be understood. Different connexin channels have distinct molecular selectivities that cannot be explained simply on the basis of size or charge of the permeants. To gain insight into the forces and interactions that underlie selective molecular permeation, we investigated the energetics of two uncharged derivatized sugars, one permeable and one impermeable, through a validated connexin26 (Cx26) channel structural model, using molecular dynamics and associated analytic tools. The system is a Cx26 channel equilibrated in explicit membrane/solvent, shown by Brownian dynamics to reproduce key conductance characteristics of the native channel. The results are consistent with the known difference in permeability to each molecule. The energetic barriers extend through most of the pore length, rather than being highly localized as in ion-specific channels. There is little evidence for binding within the pore. Force decomposition reveals how, for each tested molecule, interactions with water and the Cx26 protein vary over the length of the pore and reveals a significant contribution from hydrogen bonding and interaction with K(+). The flexibility of the pore width varies along its length, and the tested molecules have differential effects on pore width as they pass through. Potential sites of interaction within the pore are defined for each molecule. The results suggest that for the tested molecules, differences in hydrogen bonding and entropic factors arising from permeant flexibility substantially contribute to the energetics of permeation. This work highlights factors involved in selective molecular permeation that differ from those that define selectivity among atomic ions.
Collapse
|
22
|
Periole X. Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations. Chem Rev 2016; 117:156-185. [PMID: 28073248 DOI: 10.1021/acs.chemrev.6b00344] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.
Collapse
Affiliation(s)
- Xavier Periole
- Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen , Nijenborgh 7, 9747AG Groningen, The Netherlands
| |
Collapse
|
23
|
Mori T, Miyashita N, Im W, Feig M, Sugita Y. Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:1635-51. [PMID: 26766517 DOI: 10.1016/j.bbamem.2015.12.032] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/01/2015] [Revised: 12/24/2015] [Accepted: 12/29/2015] [Indexed: 12/15/2022]
Abstract
This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
Affiliation(s)
- Takaharu Mori
- iTHES Research Group and Theoretical Molecular Science Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Naoyuki Miyashita
- Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Faculty of Biology-Oriented Science and Technology, KINDAI University, 930 Nishimitani, Kinokawa, Wakayama 649-6493, Japan
| | - Wonpil Im
- Department of Molecular Sciences and Center for Computational Biology, The University of Kansas, 2030 Becker Drive, Lawrence, KS 66047, United States
| | - Michael Feig
- Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, United States; Department of Chemistry, Michigan State University, East Lansing, MI 48824, United States
| | - Yuji Sugita
- iTHES Research Group and Theoretical Molecular Science Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan; Laboratory for Biomolecular Function Simulation, RIKEN Quantitative Biology Center, Integrated Innovation Building 7F, 6-7-1 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan; Department of Chemistry, Michigan State University, East Lansing, MI 48824, United States; Computational Biophysics Research Team, RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
| |
Collapse
|
24
|
Pothula KR, Solano CJF, Kleinekathöfer U. Simulations of outer membrane channels and their permeability. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1858:1760-71. [PMID: 26721326 DOI: 10.1016/j.bbamem.2015.12.020] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2015] [Revised: 12/15/2015] [Accepted: 12/17/2015] [Indexed: 12/25/2022]
Abstract
Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
Affiliation(s)
- Karunakar R Pothula
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Carlos J F Solano
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
25
|
Nishizawa M, Nishizawa K. Potential of mean force analysis of the self-association of leucine-rich transmembrane α-helices: difference between atomistic and coarse-grained simulations. J Chem Phys 2015; 141:075101. [PMID: 25149815 DOI: 10.1063/1.4891932] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Interaction of transmembrane (TM) proteins is important in many biological processes. Large-scale computational studies using coarse-grained (CG) simulations are becoming popular. However, most CG model parameters have not fully been calibrated with respect to lateral interactions of TM peptide segments. Here, we compare the potential of mean forces (PMFs) of dimerization of TM helices obtained using a MARTINI CG model and an atomistic (AT) Berger lipids-OPLS/AA model (AT(OPLS)). For helical, tryptophan-flanked, leucine-rich peptides (WL15 and WALP15) embedded in a parallel configuration in an octane slab, the AT(OPLS) PMF profiles showed a shallow minimum (with a depth of approximately 3 kJ/mol; i.e., a weak tendency to dimerize). A similar analysis using the CHARMM36 all-atom model (AT(CHARMM)) showed comparable results. In contrast, the CG analysis generally showed steep PMF curves with depths of approximately 16-22 kJ/mol, suggesting a stronger tendency to dimerize compared to the AT model. This CG > AT discrepancy in the propensity for dimerization was also seen for dilauroylphosphatidylcholine (DLPC)-embedded peptides. For a WL15 (and WALP15)/DLPC bilayer system, AT(OPLS) PMF showed a repulsive mean force for a wide range of interhelical distances, in contrast to the attractive forces observed in the octane system. The change from the octane slab to the DLPC bilayer also mitigated the dimerization propensity in the CG system. The dimerization energies of CG (AALALAA)3 peptides in DLPC and dioleoylphosphatidylcholine bilayers were in good agreement with previous experimental data. The lipid headgroup, but not the length of the lipid tails, was a key causative factor contributing to the differences between octane and DLPC. Furthermore, the CG model, but not the AT model, showed high sensitivity to changes in amino acid residues located near the lipid-water interface and hydrophobic mismatch between the peptides and membrane. These findings may help interpret CG and AT simulation results on membrane proteins.
Collapse
Affiliation(s)
- Manami Nishizawa
- Teikyo University School of Medical Technology, Itabashi, Tokyo, Japan
| | | |
Collapse
|
26
|
Modeling transmembrane domain dimers/trimers of plexin receptors: implications for mechanisms of signal transmission across the membrane. PLoS One 2015; 10:e0121513. [PMID: 25837709 PMCID: PMC4383379 DOI: 10.1371/journal.pone.0121513] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2015] [Accepted: 02/03/2015] [Indexed: 01/01/2023] Open
Abstract
Single-pass transmembrane (TM) receptors transmit signals across lipid bilayers by helix association or by configurational changes within preformed dimers. The structure determination for such TM regions is challenging and has mostly been accomplished by NMR spectroscopy. Recently, the computational prediction of TM dimer structures is becoming recognized for providing models, including alternate conformational states, which are important for receptor regulation. Here we pursued a strategy to predict helix oligomers that is based on packing considerations (using the PREDDIMER webserver) and is followed by a refinement of structures, utilizing microsecond all-atom molecular dynamics simulations. We applied this method to plexin TM receptors, a family of 9 human proteins, involved in the regulation of cell guidance and motility. The predicted models show that, overall, the preferences identified by PREDDIMER are preserved in the unrestrained simulations and that TM structures are likely to be diverse across the plexin family. Plexin-B1 and -B3 TM helices are regular and tend to associate, whereas plexin-A1, -A2, -A3, -A4, -C1 and -D1 contain sequence elements, such as poly-Glycine or aromatic residues that distort helix conformation and association. Plexin-B2 does not form stable dimers due to the presence of TM prolines. No experimental structural information on the TM region is available for these proteins, except for plexin-C1 dimeric and plexin-B1 - trimeric structures inferred from X-ray crystal structures of the intracellular regions. Plexin-B1 TM trimers utilize Ser and Thr sidechains for interhelical contacts. We also modeled the juxta-membrane (JM) region of plexin-C1 and plexin-B1 and show that it synergizes with the TM structures. The structure and dynamics of the JM region and TM-JM junction provide determinants for the distance and distribution of the intracellular domains, and for their binding partners relative to the membrane. The structures suggest experimental tests and will be useful for the interpretation of future studies.
Collapse
|
27
|
Pawar AB, Deshpande SA, Gopal SM, Wassenaar TA, Athale CA, Sengupta D. Thermodynamic and kinetic characterization of transmembrane helix association. Phys Chem Chem Phys 2015; 17:1390-8. [DOI: 10.1039/c4cp03732d] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The transient dimerization of transmembrane proteins is an important event in several cellular processes and here we use coarse-grain and meso-scale modeling methods to quantify their underlying dynamics.
Collapse
Affiliation(s)
| | | | | | - Tsjerk A. Wassenaar
- Department of Biology
- Computational Biology
- University of Erlangen-Nürnberg
- 91058 Erlangen
- Germany
| | | | | |
Collapse
|
28
|
Effect of Lipid Bilayer Composition on Membrane Protein Association. ADVANCES IN PLANAR LIPID BILAYERS AND LIPOSOMES 2015. [DOI: 10.1016/bs.adplan.2015.06.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
29
|
Prasanna X, Chattopadhyay A, Sengupta D. Cholesterol modulates the dimer interface of the β₂-adrenergic receptor via cholesterol occupancy sites. Biophys J 2014; 106:1290-300. [PMID: 24655504 DOI: 10.1016/j.bpj.2014.02.002] [Citation(s) in RCA: 122] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2013] [Revised: 02/01/2014] [Accepted: 02/03/2014] [Indexed: 11/26/2022] Open
Abstract
The β2-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β2-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β2-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.
Collapse
|
30
|
Chavent M, Chetwynd AP, Stansfeld PJ, Sansom MSP. Dimerization of the EphA1 receptor tyrosine kinase transmembrane domain: Insights into the mechanism of receptor activation. Biochemistry 2014; 53:6641-52. [PMID: 25286141 PMCID: PMC4298228 DOI: 10.1021/bi500800x] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
EphA1
is a receptor tyrosine kinase (RTK) that plays a key role
in developmental processes, including guidance of the migration of
axons and cells in the nervous system. EphA1, in common with other
RTKs, contains an N-terminal extracellular domain, a single transmembrane
(TM) α-helix, and a C-terminal intracellular kinase domain.
The TM helix forms a dimer, as seen in recent NMR studies. We have
modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain
(CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional
potential of mean force (PMF) for this system, based on interhelix
separation, has been calculated using CG MD simulations. This provides
a view of the free energy landscape for helix–helix interactions
of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest
two states, consistent with a rotation-coupled activation mechanism.
The more stable state corresponds to a right-handed helix dimer interacting
via an N-terminal glycine zipper motif, consistent with a recent NMR
structure (2K1K). A second metastable state corresponds to a structure in which
the glycine zipper motif is not involved. Analysis of unrestrained
CG MD simulations based on representative models from the PMF calculations
or on the NMR structure reveals possible pathways of interconversion
between these two states, involving helix rotations about their long
axes. This suggests that the interaction of TM helices in EphA1 dimers
may be intrinsically dynamic. This provides a potential mechanism
for signaling whereby extracellular events drive a shift in the repopulation
of the underlying TM helix dimer energy landscape.
Collapse
Affiliation(s)
- Matthieu Chavent
- Department of Biochemistry, University of Oxford , South Parks Road, Oxford OX1 3QU, United Kingdom
| | | | | | | |
Collapse
|
31
|
Park S, Im W. Theory of Adaptive Optimization for Umbrella Sampling. J Chem Theory Comput 2014; 10:2719-2728. [PMID: 25018672 PMCID: PMC4089916 DOI: 10.1021/ct500504g] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2014] [Indexed: 01/24/2023]
Abstract
We present a theory of adaptive optimization for umbrella sampling. With the analytical bias force constant obtained from the constrained thermodynamic length along the reaction coordinate, the windows are distributed to optimize the overlap between neighbors. Combining with the replica exchange method, we propose a method of adaptive window exchange umbrella sampling. The efficiency gain in sampling by the present method originates from the optimal window distribution, in which windows are concentrated to the region where the free energy is steep, as well as consequently improved random walk.
Collapse
Affiliation(s)
- Soohyung Park
- Department of Molecular Biosciences
and Center for Bioinformatics, The University
of Kansas, 2030 Becker
Drive, Lawrence, Kansas 66047, United States
| | - Wonpil Im
- Department of Molecular Biosciences
and Center for Bioinformatics, The University
of Kansas, 2030 Becker
Drive, Lawrence, Kansas 66047, United States
| |
Collapse
|
32
|
Reddy T, Manrique S, Buyan A, Hall BA, Chetwynd A, Sansom MSP. Primary and secondary dimer interfaces of the fibroblast growth factor receptor 3 transmembrane domain: characterization via multiscale molecular dynamics simulations. Biochemistry 2014; 53:323-32. [PMID: 24397339 PMCID: PMC4871223 DOI: 10.1021/bi401576k] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Receptor tyrosine kinases are single-pass membrane proteins that form dimers within the membrane. The interactions of their transmembrane domains (TMDs) play a key role in dimerization and signaling. Fibroblast growth factor receptor 3 (FGFR3) is of interest as a G380R mutation in its TMD is the underlying cause of ~99% of the cases of achondroplasia, the most common form of human dwarfism. The structural consequences of this mutation remain uncertain: the mutation shifts the position of the TMD relative to the lipid bilayer but does not alter the association free energy. We have combined coarse-grained and all-atom molecular dynamics simulations to study the dimerization of wild-type, heterodimer, and mutant FGFR3 TMDs. The simulations reveal that the helices pack together in the dimer to form a flexible interface. The primary packing mode is mediated by a Gx3G motif. There is also a secondary dimer interface that is more highly populated in heterodimer and mutant configurations that may feature in the molecular mechanism of pathology. Both coarse-grained and atomistic simulations reveal a significant shift of the G380R mutant dimer TMD relative to the bilayer to allow interactions of the arginine side chain with lipid headgroup phosphates.
Collapse
Affiliation(s)
- Tyler Reddy
- Department of Biochemistry, University of Oxford , South Parks Road, Oxford OX1 3QU, U.K
| | | | | | | | | | | |
Collapse
|