1
|
Versini R, Baaden M, Cavellini L, Cohen MM, Taly A, Fuchs PFJ. Lys716 in the transmembrane domain of yeast mitofusin Fzo1 modulates anchoring and fusion. Structure 2024:S0969-2126(24)00334-4. [PMID: 39299234 DOI: 10.1016/j.str.2024.08.017] [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: 12/22/2023] [Revised: 05/06/2024] [Accepted: 08/23/2024] [Indexed: 09/22/2024]
Abstract
Outer mitochondrial membrane fusion, a vital cellular process, is mediated by mitofusins. However, the underlying molecular mechanism remains elusive. We have performed extensive multiscale molecular dynamics simulations to predict a model of the transmembrane (TM) domain of the yeast mitofusin Fzo1. Coarse-grained simulations of the two TM domain helices, TM1 and TM2, reveal a stable interface, which is controlled by the charge status of residue Lys716. Atomistic replica-exchange simulations further tune our model, which is confirmed by a remarkable agreement with an independent AlphaFold2 (AF2) prediction of Fzo1 in complex with its fusion partner Ugo1. Furthermore, the presence of the TM domain destabilizes the membrane, even more if Lys716 is charged, which can be an asset for initiating fusion. The functional role of Lys716 was confirmed with yeast experiments, which show that mutating Lys716 to a hydrophobic residue prevents mitochondrial fusion.
Collapse
Affiliation(s)
- Raphaëlle Versini
- Laboratoire de Biochimie Théorique, CNRS, Université Paris Cité, 75005 Paris, France; Laboratoire des Biomolécules, LBM, Sorbonne Université, École normale supérieure, PSL University, CNRS, 75005 Paris, France
| | - Marc Baaden
- Laboratoire de Biochimie Théorique, CNRS, Université Paris Cité, 75005 Paris, France
| | - Laetitia Cavellini
- Laboratoire de Biologie Cellulaire et Moléculaire des Eucaryotes, Institut de Biologie Physico-Chimique, UMR 8226, CNRS, Sorbonne Université, Paris, France
| | - Mickaël M Cohen
- Laboratoire de Biologie Cellulaire et Moléculaire des Eucaryotes, Institut de Biologie Physico-Chimique, UMR 8226, CNRS, Sorbonne Université, Paris, France
| | - Antoine Taly
- Laboratoire de Biochimie Théorique, CNRS, Université Paris Cité, 75005 Paris, France.
| | - Patrick F J Fuchs
- Laboratoire des Biomolécules, LBM, Sorbonne Université, École normale supérieure, PSL University, CNRS, 75005 Paris, France; Université Paris Cité, 75006 Paris, France.
| |
Collapse
|
2
|
Liang J, Menon A, Tomco T, Bhattarai N, Smith IN, Khrestian M, Formica SV, Eng C, Buck M, Bekris LM. A Computational Approach in the Systematic Search of the Interaction Partners of Alternatively Spliced TREM2 Isoforms. Int J Mol Sci 2024; 25:9667. [PMID: 39273614 PMCID: PMC11395018 DOI: 10.3390/ijms25179667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2024] [Revised: 08/30/2024] [Accepted: 09/02/2024] [Indexed: 09/15/2024] Open
Abstract
Alzheimer's disease is the most common form of dementia, characterized by the pathological accumulation of amyloid-beta (Aβ) plaques and tau neurofibrillary tangles. Triggering receptor expressed on myeloid cells 2 (TREM2) is increasingly recognized as playing a central role in Aβ clearance and microglia activation in AD. The TREM2 gene transcriptional product is alternatively spliced to produce three different protein isoforms. The canonical TREM2 isoform binds to DAP12 to activate downstream pathways. However, little is known about the function or interaction partners of the alternative TREM2 isoforms. The present study utilized a computational approach in a systematic search for new interaction partners of the TREM2 isoforms by integrating several state-of-the-art structural bioinformatics tools from initial large-scale screening to one-on-one corroborative modeling and eventual all-atom visualization. CD9, a cell surface glycoprotein involved in cell-cell adhesion and migration, was identified as a new interaction partner for two TREM2 isoforms, and CALM, a calcium-binding protein involved in calcium signaling, was identified as an interaction partner for a third TREM2 isoform, highlighting the potential role of cell adhesion and calcium regulation in AD.
Collapse
Affiliation(s)
- Junyi Liang
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Aditya Menon
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Taylor Tomco
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Nisha Bhattarai
- Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Iris Nira Smith
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Maria Khrestian
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Shane V Formica
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Charis Eng
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
- Cleveland Clinic Lerner College of Medicine, Cleveland, OH 44195, USA
- Department of Neurosciences, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, Cleveland, OH 44106, USA
| | - Lynn M Bekris
- Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA
- Cleveland Clinic Lerner College of Medicine, Cleveland, OH 44195, USA
- Department of Neurosciences, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| |
Collapse
|
3
|
Majumder A, Straub JE. Machine Learning Derived Collective Variables for the Study of Protein Homodimerization in Membrane. J Chem Theory Comput 2024; 20:5774-5783. [PMID: 38918177 DOI: 10.1021/acs.jctc.4c00454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
The accurate calculation of equilibrium constants for protein-protein association is of fundamental importance to quantitative biology and remains an outstanding challenge for computational biophysics. Traditionally, equilibrium constants have been computed from one-dimensional free energy surfaces derived from sampling along a single collective variable. Importantly, recent advances in enhanced sampling methodology have facilitated the characterization of multidimensional free energy landscapes, often exposing multiple thermodynamically important minima missed by more restrictive sampling methods. A key to the effectiveness of this multidimensional sampling approach is the identification of collective variables that effectively define the configurational space of dissociated and associated states. Here we present the application of two machine learning methods for the unbiased determination of collective variables for enhanced sampling and analysis of protein-protein association. Our results both validate prior work, based on intuition derived collective variables, and demonstrate the effectiveness of the machine learning methods for the identification of collective variables for association reactions in complex biomolecular systems.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
4
|
Borges-Araújo L, Pereira GP, Valério M, Souza PCT. Assessing the Martini 3 protein model: A review of its path and potential. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2024; 1872:141014. [PMID: 38670324 DOI: 10.1016/j.bbapap.2024.141014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/13/2024] [Accepted: 04/09/2024] [Indexed: 04/28/2024]
Abstract
Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Gilberto P Pereira
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Mariana Valério
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France
| | - Paulo C T Souza
- Laboratoire de Biologie et Modélisation de la Cellule, CNRS, UMR 5239, Inserm, U1293, Université Claude Bernard Lyon 1, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France; Centre Blaise Pascal de Simulation et de Modélisation Numérique, Ecole Normale Supérieure de Lyon, 46 Allée d'Italie, 69364 Lyon, France.
| |
Collapse
|
5
|
Versini R, Sritharan S, Aykac Fas B, Tubiana T, Aimeur SZ, Henri J, Erard M, Nüsse O, Andreani J, Baaden M, Fuchs P, Galochkina T, Chatzigoulas A, Cournia Z, Santuz H, Sacquin-Mora S, Taly A. A Perspective on the Prospective Use of AI in Protein Structure Prediction. J Chem Inf Model 2024; 64:26-41. [PMID: 38124369 DOI: 10.1021/acs.jcim.3c01361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
AlphaFold2 (AF2) and RoseTTaFold (RF) have revolutionized structural biology, serving as highly reliable and effective methods for predicting protein structures. This article explores their impact and limitations, focusing on their integration into experimental pipelines and their application in diverse protein classes, including membrane proteins, intrinsically disordered proteins (IDPs), and oligomers. In experimental pipelines, AF2 models help X-ray crystallography in resolving the phase problem, while complementarity with mass spectrometry and NMR data enhances structure determination and protein flexibility prediction. Predicting the structure of membrane proteins remains challenging for both AF2 and RF due to difficulties in capturing conformational ensembles and interactions with the membrane. Improvements in incorporating membrane-specific features and predicting the structural effect of mutations are crucial. For intrinsically disordered proteins, AF2's confidence score (pLDDT) serves as a competitive disorder predictor, but integrative approaches including molecular dynamics (MD) simulations or hydrophobic cluster analyses are advocated for accurate dynamics representation. AF2 and RF show promising results for oligomeric models, outperforming traditional docking methods, with AlphaFold-Multimer showing improved performance. However, some caveats remain in particular for membrane proteins. Real-life examples demonstrate AF2's predictive capabilities in unknown protein structures, but models should be evaluated for their agreement with experimental data. Furthermore, AF2 models can be used complementarily with MD simulations. In this Perspective, we propose a "wish list" for improving deep-learning-based protein folding prediction models, including using experimental data as constraints and modifying models with binding partners or post-translational modifications. Additionally, a meta-tool for ranking and suggesting composite models is suggested, driving future advancements in this rapidly evolving field.
Collapse
Affiliation(s)
- Raphaelle Versini
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Sujith Sritharan
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Burcu Aykac Fas
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Thibault Tubiana
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
| | - Sana Zineb Aimeur
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Julien Henri
- Sorbonne Université, CNRS, Laboratoire de Biologie, Computationnelle et Quantitative UMR 7238, Institut de Biologie Paris-Seine, 4 Place Jussieu, F-75005 Paris, France
| | - Marie Erard
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Oliver Nüsse
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, 91405 Orsay, France
| | - Jessica Andreani
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
| | - Marc Baaden
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Patrick Fuchs
- Sorbonne Université, École Normale Supérieure, PSL University, CNRS, Laboratoire des Biomolécules, LBM, 75005 Paris, France
- Université de Paris, UFR Sciences du Vivant, 75013 Paris, France
| | - Tatiana Galochkina
- Université Paris Cité and Université des Antilles and Université de la Réunion, INSERM, BIGR, F-75014 Paris, France
| | - Alexios Chatzigoulas
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
- Department of Informatics and Telecommunications, National and Kapodistrian University of Athens, 15784 Athens, Greece
| | - Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
- Department of Informatics and Telecommunications, National and Kapodistrian University of Athens, 15784 Athens, Greece
| | - Hubert Santuz
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Sophie Sacquin-Mora
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| | - Antoine Taly
- Laboratoire de Biochimie Théorique, CNRS (UPR9080), Université Paris Cité, F-75005 Paris, France
| |
Collapse
|
6
|
Blazhynska M, Gumbart JC, Chen H, Tajkhorshid E, Roux B, Chipot C. A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes. J Chem Theory Comput 2023; 19:9077-9092. [PMID: 38091976 PMCID: PMC11145395 DOI: 10.1021/acs.jctc.3c00941] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.
Collapse
Affiliation(s)
- Marharyta Blazhynska
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, 837 State Street, Atlanta, Georgia 30332, United States
| | - Haochuan Chen
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 N. Mathews Avenue, Urbana, Illinois 61801, United States
- Department of Biochemistry, University of Illinois at Urbana-Champaign, 600 S. Mathews Avenue, Urbana, Illinois 61801, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 E. 57th Street W225, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n°7019, Université de Lorraine, B.P. 70239, Vandœuvre-lès-Nancy cedex 54506, France
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 N. Mathews Avenue, Urbana, Illinois 61801, United States
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 E. 57th Street W225, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Hawai'i at Ma̅noa, 2545 McCarthy Mall, Honolulu, Hawaii 96822, United States
| |
Collapse
|
7
|
Majumder A, Straub JE. The role of structural heterogeneity in the homodimerization of transmembrane proteins. J Chem Phys 2023; 159:134101. [PMID: 37782254 PMCID: PMC10547497 DOI: 10.1063/5.0159801] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Accepted: 07/27/2023] [Indexed: 10/03/2023] Open
Abstract
The equilibrium association of transmembrane proteins plays a fundamental role in membrane protein function and cellular signaling. While the study of the equilibrium binding of single pass transmembrane proteins has received significant attention in experiment and simulation, the accurate assessment of equilibrium association constants remains a challenge to experiment and simulation. In experiment, there remain wide variations in association constants derived from experimental studies of the most widely studied transmembrane proteins. In simulation, state-of-the art methods have failed to adequately sample the thermodynamically relevant structures of the dimer state ensembles using coarse-grained models. In addition, all-atom force fields often fail to accurately assess the relative free energies of the dimer and monomer states. Given the importance of this fundamental biophysical process, it is essential to address these shortcomings. In this work, we establish an effective computational protocol for the calculation of equilibrium association constants for transmembrane homodimer formation. A set of transmembrane protein homodimers, used in the parameterization of the MARTINI v3 force field, are simulated using metadynamics, based on three collective variables. The method is found to be accurate and computationally efficient, providing a standard to be used in the future simulation studies using coarse-grained or all-atom models.
Collapse
Affiliation(s)
- Ayan Majumder
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - John E. Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|