1
|
Teruel N, Borges VM, Najmanovich R. Surfaces: a software to quantify and visualize interactions within and between proteins and ligands. Bioinformatics 2023; 39:btad608. [PMID: 37788107 PMCID: PMC10568369 DOI: 10.1093/bioinformatics/btad608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2023] [Revised: 08/23/2023] [Accepted: 09/29/2023] [Indexed: 10/05/2023] Open
Abstract
SUMMARY Computational methods for the quantification and visualization of the relative contribution of molecular interactions to the stability of biomolecular structures and complexes are fundamental to understand, modulate and engineer biological processes. Here, we present Surfaces, an easy to use, fast and customizable software for quantification and visualization of molecular interactions based on the calculation of surface areas in contact. Surfaces calculations shows equivalent or better correlations with experimental data as computationally expensive methods based on molecular dynamics. AVAILABILITY AND IMPLEMENTATION All scripts are available at https://github.com/NRGLab/Surfaces. Surface's documentation is available at https://surfaces-tutorial.readthedocs.io/en/latest/index.html.
Collapse
Affiliation(s)
- Natália Teruel
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal H3T 1J4, Canada
| | - Vinicius Magalhães Borges
- Department of Biomedical Sciences, Joan C. Edwards School of Medicine, Marshall University, Huntington, WV, USA
| | - Rafael Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal H3T 1J4, Canada
| |
Collapse
|
2
|
Teruel N, Mailhot O, Najmanovich RJ. Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants. PLoS Comput Biol 2021; 17:e1009286. [PMID: 34351895 PMCID: PMC8384204 DOI: 10.1371/journal.pcbi.1009286] [Citation(s) in RCA: 54] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 08/24/2021] [Accepted: 07/17/2021] [Indexed: 01/21/2023] Open
Abstract
The SARS-CoV-2 Spike protein needs to be in an open-state conformation to interact with ACE2 to initiate viral entry. We utilise coarse-grained normal mode analysis to model the dynamics of Spike and calculate transition probabilities between states for 17081 variants including experimentally observed variants. Our results correctly model an increase in open-state occupancy for the more infectious D614G via an increase in flexibility of the closed-state and decrease of flexibility of the open-state. We predict the same effect for several mutations on glycine residues (404, 416, 504, 252) as well as residues K417, D467 and N501, including the N501Y mutation recently observed within the B.1.1.7, 501.V2 and P1 strains. This is, to our knowledge, the first use of normal mode analysis to model conformational state transitions and the effect of mutations on such transitions. The specific mutations of Spike identified here may guide future studies to increase our understanding of SARS-CoV-2 infection mechanisms and guide public health in their surveillance efforts. The present work explores the molecular mechanisms underlying and potentially helping new strains of SARS-CoV-2 to gain an evolutionary advantage during the ongoing COVID-19 pandemics. We show how a computational method called normal mode analysis that treats protein dynamics in a simplified manner is capable to predict the higher propensity of the Spike protein to be in the open state in which it is capable to interact with the human ACE2 receptor and thus facilitate cell entry. Because the simulation of the simplified computational model is relatively less demanding on resources than alternative methods, we were able to simulate over 17000 mutations in the SARS-CoV-2 Spike protein to identify multiple mutations that if they were to appear as the virus continues to evolve, could confer an evolutionary advantage. As a matter of fact, our predictions foresaw the emergence of particular mutations such as N501Y that appeared in several variants of concern. Our results can inform public health regarding new variants and serves as a proof of concept for the application of normal mode analysis to study the effect of mutations on both, protein dynamics and conformational transitions in a high-throughput manner.
Collapse
Affiliation(s)
- Natália Teruel
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Olivier Mailhot
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
- Institute for Research in Immunology and Cancer (IRIC), Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Rafael J. Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
- * E-mail:
| |
Collapse
|
3
|
Karami Y, Bitard-Feildel T, Laine E, Carbone A. "Infostery" analysis of short molecular dynamics simulations identifies highly sensitive residues and predicts deleterious mutations. Sci Rep 2018; 8:16126. [PMID: 30382169 PMCID: PMC6208415 DOI: 10.1038/s41598-018-34508-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 10/17/2018] [Indexed: 11/09/2022] Open
Abstract
Characterizing a protein mutational landscape is a very challenging problem in Biology. Many disease-associated mutations do not seem to produce any effect on the global shape nor motions of the protein. Here, we use relatively short all-atom biomolecular simulations to predict mutational outcomes and we quantitatively assess the predictions on several hundreds of mutants. We perform simulations of the wild type and 175 mutants of PSD95’s third PDZ domain in complex with its cognate ligand. By recording residue displacements correlations and interactions, we identify “communication pathways” and quantify them to predict the severity of the mutations. Moreover, we show that by exploiting simulations of the wild type, one can detect 80% of the positions highly sensitive to mutations with a precision of 89%. Importantly, our analysis describes the role of these positions in the inter-residue communication and dynamical architecture of the complex. We assess our approach on three different systems using data from deep mutational scanning experiments and high-throughput exome sequencing. We refer to our analysis as “infostery”, from “info” - information - and “steric” - arrangement of residues in space. We provide a fully automated tool, COMMA2 (www.lcqb.upmc.fr/COMMA2), that can be used to guide medicinal research by selecting important positions/mutations.
Collapse
Affiliation(s)
- Yasaman Karami
- Sorbonne Université, CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), 75005, Paris, France
| | - Tristan Bitard-Feildel
- Sorbonne Université, CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), 75005, Paris, France.,Sorbonne Université, Institut des Sciences du Calcul et de des Données (ISCD), Paris, France
| | - Elodie Laine
- Sorbonne Université, CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), 75005, Paris, France.
| | - Alessandra Carbone
- Sorbonne Université, CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), 75005, Paris, France. .,Institut Universitaire de France (IUF), Paris, France.
| |
Collapse
|
4
|
Functional diversity of TMPRSS6 isoforms and variants expressed in hepatocellular carcinoma cell lines. Sci Rep 2018; 8:12562. [PMID: 30135444 PMCID: PMC6105633 DOI: 10.1038/s41598-018-30618-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2018] [Accepted: 08/02/2018] [Indexed: 02/06/2023] Open
Abstract
TMPRSS6, also known as matriptase-2, is a type II transmembrane serine protease that plays a major role in iron homeostasis by acting as a negative regulator of hepcidin production through cleavage of the BMP co-receptor haemojuvelin. Iron-refractory iron deficiency anaemia (IRIDA), an iron metabolism disorder, is associated with mutations in the TMPRSS6 gene. By analysing RNA-seq data encoding TMPRSS6 isoforms and other proteins involved in hepcidin production, we uncovered significant differences in expression levels between hepatocellular carcinoma (HCC) cell lines and normal human liver samples. Most notably, TMPRSS6 and HAMP expression was found to be much lower in HepG2 and Huh7 cells when compared to human liver samples. Furthermore, we characterized the common TMPRSS6 polymorphism V736A identified in Hep3B cells, the V795I mutation found in HepG2 cells, also associated with IRIDA, and the G603R substitution recently detected in two IRIDA patients. While variant V736A is as active as wild-type TMPRSS6, mutants V795I and G603R displayed significantly reduced proteolytic activity. Our results provide important information about commonly used liver cell models and shed light on the impact of two TMPRSS6 mutations associated with IRIDA.
Collapse
|
5
|
Chiang RZH, Gan SKE, Su CTT. A computational study for rational HIV-1 non-nucleoside reverse transcriptase inhibitor selection and the discovery of novel allosteric pockets for inhibitor design. Biosci Rep 2018; 38:BSR20171113. [PMID: 29437904 PMCID: PMC5835713 DOI: 10.1042/bsr20171113] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 01/29/2018] [Accepted: 01/29/2018] [Indexed: 12/15/2022] Open
Abstract
HIV drug resistant mutations that render the current Highly Active Anti-Retroviral Therapy (HAART) cocktail drugs ineffective are increasingly reported. To study the mechanisms of these mutations in conferring drug resistance, we computationally analyzed 14 reverse transcriptase (RT) structures of HIV-1 on the following parameters: drug-binding pocket volume, allosteric effects caused by the mutations, and structural thermal stability. We constructed structural correlation-based networks of the mutant RT-drug complexes and the analyses support the use of efavirenz (EFZ) as the first-line drug, given that cross-resistance is least likely to develop from EFZ-resistant mutations. On the other hand, rilpivirine (RPV)-resistant mutations showed the highest cross-resistance to the other non-nucleoside RT inhibitors. With significant drug cross-resistance associated with the known allosteric drug-binding site, there is a need to identify new allosteric druggable sites in the structure of RT. Through computational analyses, we found such a novel druggable pocket on the HIV-1 RT structure that is comparable with the original allosteric drug site, opening the possibility to the design of new inhibitors.
Collapse
Affiliation(s)
- Ron Zhi-Hui Chiang
- Bioinformatics Institute, Agency for Science, Technology, and Research (A*STAR), Singapore 138671
| | - Samuel Ken-En Gan
- Bioinformatics Institute, Agency for Science, Technology, and Research (A*STAR), Singapore 138671
- p53 Laboratory, Agency for Science, Technology, and Research (A*STAR), Singapore 138648
| | - Chinh Tran-To Su
- Bioinformatics Institute, Agency for Science, Technology, and Research (A*STAR), Singapore 138671
| |
Collapse
|
6
|
Vedithi SC, Malhotra S, Das M, Daniel S, Kishore N, George A, Arumugam S, Rajan L, Ebenezer M, Ascher DB, Arnold E, Blundell TL. Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae. Sci Rep 2018; 8:5016. [PMID: 29567948 PMCID: PMC5864748 DOI: 10.1038/s41598-018-23423-1] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2017] [Accepted: 03/09/2018] [Indexed: 11/09/2022] Open
Abstract
The rpoB gene encodes the β subunit of RNA polymerase holoenzyme in Mycobacterium leprae (M. leprae). Missense mutations in the rpoB gene were identified as etiological factors for rifampin resistance in leprosy. In the present study, we identified mutations corresponding to rifampin resistance in relapsed leprosy cases from three hospitals in southern India which treat leprosy patients. DNA was extracted from skin biopsies of 35 relapse/multidrug therapy non-respondent leprosy cases, and PCR was performed to amplify the 276 bp rifampin resistance-determining region of the rpoB gene. PCR products were sequenced, and mutations were identified in four out of the 35 cases at codon positions D441Y, D441V, S437L and H476R. The structural and functional effects of these mutations were assessed in the context of three-dimensional comparative models of wild-type and mutant M. leprae RNA polymerase holoenzyme (RNAP), based on the recently solved crystal structures of RNAP of Mycobacterium tuberculosis, containing a synthetic nucleic acid scaffold and rifampin. The resistance mutations were observed to alter the hydrogen-bonding and hydrophobic interactions of rifampin and the 5' ribonucleotide of the growing RNA transcript. This study demonstrates that rifampin-resistant strains of M. leprae among leprosy patients in southern India are likely to arise from mutations that affect the drug-binding site and stability of RNAP.
Collapse
Affiliation(s)
- Sundeep Chaitanya Vedithi
- Department of Biochemistry, University of Cambridge, Tennis Court Rd., Cambridge, CB2 1GA, UK. .,Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India.
| | - Sony Malhotra
- Department of Biochemistry, University of Cambridge, Tennis Court Rd., Cambridge, CB2 1GA, UK
| | - Madhusmita Das
- Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India
| | - Sheela Daniel
- Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India
| | - Nanda Kishore
- Department of Dermatology, Father Muller Medical College & Hospital, Mangalore, Karnataka, 575 002, India
| | - Anuja George
- Department of Dermatology, Trivandrum Medical College, Trivandrum, Kerala, 695011, India
| | - Shantha Arumugam
- Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India
| | - Lakshmi Rajan
- Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India
| | - Mannam Ebenezer
- Molecular Biology and Immunology Division, Schieffelin Institute of Health Research & Leprosy Center (SIH R & LC), Karigiri, Vellore, Tamil Nadu, 632106, India
| | - David B Ascher
- Department of Biochemistry and Molecular Biology, Bio21 Institute, University of Melbourne, Parkville, VIC 3052, Australia
| | - Eddy Arnold
- Center for Advanced Biotechnology and Medicine (CABM), and Rutgers University Department of Chemistry and Chemical Biology, 679 Hoes Lane, Piscataway, NJ, 08854, USA
| | - Tom L Blundell
- Department of Biochemistry, University of Cambridge, Tennis Court Rd., Cambridge, CB2 1GA, UK.
| |
Collapse
|