1
|
Rodríguez-Salazar CA, van Tol S, Mailhot O, Gonzalez-Orozco M, Galdino GT, Warren AN, Teruel N, Behera P, Afreen KS, Zhang L, Juelich TL, Smith JK, Zylber MI, Freiberg AN, Najmanovich RJ, Giraldo MI, Rajsbaum R. Ebola virus VP35 interacts non-covalently with ubiquitin chains to promote viral replication. PLoS Biol 2024; 22:e3002544. [PMID: 38422166 PMCID: PMC10942258 DOI: 10.1371/journal.pbio.3002544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 03/15/2024] [Accepted: 02/09/2024] [Indexed: 03/02/2024] Open
Abstract
Ebolavirus (EBOV) belongs to a family of highly pathogenic viruses that cause severe hemorrhagic fever in humans. EBOV replication requires the activity of the viral polymerase complex, which includes the cofactor and Interferon antagonist VP35. We previously showed that the covalent ubiquitination of VP35 promotes virus replication by regulating interactions with the polymerase complex. In addition, VP35 can also interact non-covalently with ubiquitin (Ub); however, the function of this interaction is unknown. Here, we report that VP35 interacts with free (unanchored) K63-linked polyUb chains. Ectopic expression of Isopeptidase T (USP5), which is known to degrade unanchored polyUb chains, reduced VP35 association with Ub and correlated with diminished polymerase activity in a minigenome assay. Using computational methods, we modeled the VP35-Ub non-covalent interacting complex, identified the VP35-Ub interacting surface, and tested mutations to validate the interface. Docking simulations identified chemical compounds that can block VP35-Ub interactions leading to reduced viral polymerase activity. Treatment with the compounds reduced replication of infectious EBOV in cells and in vivo in a mouse model. In conclusion, we identified a novel role of unanchored polyUb in regulating Ebola virus polymerase function and discovered compounds that have promising anti-Ebola virus activity.
Collapse
Affiliation(s)
- Carlos A. Rodríguez-Salazar
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
- Molecular Biology and Virology Laboratory, Faculty of Medicine and Health Sciences, Corporación Universitaria Empresarial Alexander von Humboldt, Armenia, Colombia
| | - Sarah van Tol
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Olivier Mailhot
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Maria Gonzalez-Orozco
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Gabriel T. Galdino
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Abbey N. Warren
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey, United States of America
| | - Natalia Teruel
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Padmanava Behera
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey, United States of America
| | - Kazi Sabrina Afreen
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey, United States of America
| | - Lihong Zhang
- Department of Pathology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Terry L. Juelich
- Department of Pathology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Jennifer K. Smith
- Department of Pathology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - María Inés Zylber
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Alexander N. Freiberg
- Department of Pathology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Rafael J. Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Maria I. Giraldo
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
| | - Ricardo Rajsbaum
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston, Texas, United States of America
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey, United States of America
| |
Collapse
|
2
|
Rodríguez-Salazar CA, van Tol S, Mailhot O, Galdino G, Teruel N, Zhang L, Warren AN, González-Orozco M, Freiberg AN, Najmanovich RJ, Giraldo MI, Rajsbaum R. Ebola Virus VP35 Interacts Non-Covalently with Ubiquitin Chains to Promote Viral Replication Creating New Therapeutic Opportunities. bioRxiv 2023:2023.07.14.549057. [PMID: 37503276 PMCID: PMC10369991 DOI: 10.1101/2023.07.14.549057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Ebolavirus (EBOV) belongs to a family of highly pathogenic viruses that cause severe hemorrhagic fever in humans. EBOV replication requires the activity of the viral polymerase complex, which includes the co-factor and Interferon antagonist VP35. We previously showed that the covalent ubiquitination of VP35 promotes virus replication by regulating interactions with the polymerase complex. In addition, VP35 can also interact non-covalently with ubiquitin (Ub); however, the function of this interaction is unknown. Here, we report that VP35 interacts with free (unanchored) K63-linked polyUb chains. Ectopic expression of Isopeptidase T (USP5), which is known to degrade unanchored polyUb chains, reduced VP35 association with Ub and correlated with diminished polymerase activity in a minigenome assay. Using computational methods, we modeled the VP35-Ub non-covalent interacting complex, identified the VP35-Ub interacting surface and tested mutations to validate the interface. Docking simulations identified chemical compounds that can block VP35-Ub interactions leading to reduced viral polymerase activity that correlated with reduced replication of infectious EBOV. In conclusion, we identified a novel role of unanchored polyUb in regulating Ebola virus polymerase function and discovered compounds that have promising anti-Ebola virus activity.
Collapse
Affiliation(s)
- Carlos A. Rodríguez-Salazar
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
- Molecular Biology and Virology Laboratory, Faculty of Medicine and Health Sciences, Corporación Universitaria Empresarial Alexander von Humboldt, Armenia 630003, Colombia
| | - Sarah van Tol
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
| | - Olivier Mailhot
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Gabriel Galdino
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Natalia Teruel
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - Lihong Zhang
- Department of Pathology, University of Texas Medical Branch, Galveston 77555, Texas, USA
| | - Abbey N. Warren
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey 07103
| | - María González-Orozco
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
| | - Alexander N. Freiberg
- Department of Pathology, University of Texas Medical Branch, Galveston 77555, Texas, USA
| | - Rafael J. Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, Canada
| | - María I. Giraldo
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
| | - Ricardo Rajsbaum
- Department of Microbiology and Immunology, University of Texas Medical Branch, Galveston 77555, Texas, USA
- Center for Virus-Host-Innate Immunity and Department of Medicine; Rutgers Biomedical and Health Sciences, Institute for Infectious and Inflammatory Diseases, Rutgers University, Newark, New Jersey 07103
| |
Collapse
|
3
|
Mailhot O, Frappier V, Major F, Najmanovich RJ. Sequence-sensitive elastic network captures dynamical features necessary for miR-125a maturation. PLoS Comput Biol 2022; 18:e1010777. [PMID: 36516216 PMCID: PMC9797095 DOI: 10.1371/journal.pcbi.1010777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 12/28/2022] [Accepted: 11/29/2022] [Indexed: 12/15/2022] Open
Abstract
The Elastic Network Contact Model (ENCoM) is a coarse-grained normal mode analysis (NMA) model unique in its all-atom sensitivity to the sequence of the studied macromolecule and thus to the effect of mutations. We adapted ENCoM to simulate the dynamics of ribonucleic acid (RNA) molecules, benchmarked its performance against other popular NMA models and used it to study the 3D structural dynamics of human microRNA miR-125a, leveraging high-throughput experimental maturation efficiency data of over 26 000 sequence variants. We also introduce a novel way of using dynamical information from NMA to train multivariate linear regression models, with the purpose of highlighting the most salient contributions of dynamics to function. ENCoM has a similar performance profile on RNA than on proteins when compared to the Anisotropic Network Model (ANM), the most widely used coarse-grained NMA model; it has the advantage on predicting large-scale motions while ANM performs better on B-factors prediction. A stringent benchmark from the miR-125a maturation dataset, in which the training set contains no sequence information in common with the testing set, reveals that ENCoM is the only tested model able to capture signal beyond the sequence. This ability translates to better predictive power on a second benchmark in which sequence features are shared between the train and test sets. When training the linear regression model using all available data, the dynamical features identified as necessary for miR-125a maturation point to known patterns but also offer new insights into the biogenesis of microRNAs. Our novel approach combining NMA with multivariate linear regression is generalizable to any macromolecule for which relatively high-throughput mutational data is available.
Collapse
Affiliation(s)
- Olivier Mailhot
- Department of Biochemistry and Molecular Medicine, Université de Montréal, Montreal, QC, Canada
- Department of Computer Science and Operations Research, Université de Montréal, Montreal, QC, Canada
- Institute for Research in Immunology and Cancer, Université de Montréal, Montreal, QC, Canada
- Department of Pharmacology and Physiology, Université de Montréal, Montreal, QC, Canada
| | - Vincent Frappier
- Generate Biomedicines, Cambridge, Massachusetts, United States of America
| | - François Major
- Department of Computer Science and Operations Research, Université de Montréal, Montreal, QC, Canada
- Institute for Research in Immunology and Cancer, Université de Montréal, Montreal, QC, Canada
- * E-mail: (FM); (RJN)
| | - Rafael J. Najmanovich
- Department of Pharmacology and Physiology, Université de Montréal, Montreal, QC, Canada
- * E-mail: (FM); (RJN)
| |
Collapse
|
4
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
5
|
Yan LH, Le Roux A, Boyapelly K, Lamontagne AM, Archambault MA, Picard-Jean F, Lalonde-Seguin D, St-Pierre E, Najmanovich RJ, Fortier LC, Lafontaine D, Marsault É. Purine analogs targeting the guanine riboswitch as potential antibiotics against Clostridioides difficile. Eur J Med Chem 2017; 143:755-768. [PMID: 29220796 DOI: 10.1016/j.ejmech.2017.11.079] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2017] [Revised: 11/22/2017] [Accepted: 11/27/2017] [Indexed: 12/14/2022]
Abstract
Riboswitches recently emerged as possible targets for the development of alternative antimicrobial approaches. Guanine-sensing riboswitches in the bacterial pathogen Clostridioides difficile (formerly known as Clostridium difficile) constitute potential targets based on their involvement in the regulation of basal metabolic control of purine compounds. In this study, we deciphered the structure-activity relationship of several guanine derivatives on the guanine riboswitch and determined their antimicrobial activity. We describe the synthesis of purine analogs modified in ring B as well as positions 2 and 6. Their biological activity was determined by measuring their affinity for the C. difficile guanine riboswitch and their inhibitory effect on bacterial growth, including a counter-screen to discriminate against riboswitch-independent antibacterial effects. Altogether, our results suggest that improvements in riboswitch binding affinity in vitro do not necessarily translate into improved antibacterial activity in bacteria, despite the fact that some structure-activity relationship was observed at least with respect to binding affinity.
Collapse
Affiliation(s)
- Lok-Hang Yan
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Pharmacology-Physiology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Antoine Le Roux
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Pharmacology-Physiology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Kumaraswamy Boyapelly
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Pharmacology-Physiology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Anne-Marie Lamontagne
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Biology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Marie-Ann Archambault
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Biology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Frédéric Picard-Jean
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Biology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - David Lalonde-Seguin
- Department of Microbiology and Infectious Diseases, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Emilie St-Pierre
- Department of Microbiology and Infectious Diseases, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Rafael J Najmanovich
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Biochemistry, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada
| | - Louis-Charles Fortier
- Department of Microbiology and Infectious Diseases, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada.
| | - Daniel Lafontaine
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Biology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada.
| | - Éric Marsault
- Institut de Pharmacologie de Sherbrooke, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada; Department of Pharmacology-Physiology, Université de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec, J1H 5N4, Canada.
| |
Collapse
|
6
|
Chénard T, Guénard F, Vohl MC, Carpentier A, Tchernof A, Najmanovich RJ. Remodeling adipose tissue through in silico modulation of fat storage for the prevention of type 2 diabetes. BMC Syst Biol 2017; 11:60. [PMID: 28606124 PMCID: PMC5468946 DOI: 10.1186/s12918-017-0438-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/22/2016] [Accepted: 06/05/2017] [Indexed: 11/10/2022]
Abstract
BACKGROUND Type 2 diabetes is one of the leading non-infectious diseases worldwide and closely relates to excess adipose tissue accumulation as seen in obesity. Specifically, hypertrophic expansion of adipose tissues is related to increased cardiometabolic risk leading to type 2 diabetes. Studying mechanisms underlying adipocyte hypertrophy could lead to the identification of potential targets for the treatment of these conditions. RESULTS We present iTC1390adip, a highly curated metabolic network of the human adipocyte presenting various improvements over the previously published iAdipocytes1809. iTC1390adip contains 1390 genes, 4519 reactions and 3664 metabolites. We validated the network obtaining 92.6% accuracy by comparing experimental gene essentiality in various cell lines to our predictions of biomass production. Using flux balance analysis under various test conditions, we predict the effect of gene deletion on both lipid droplet and biomass production, resulting in the identification of 27 genes that could reduce adipocyte hypertrophy. We also used expression data from visceral and subcutaneous adipose tissues to compare the effect of single gene deletions between adipocytes from each compartment. CONCLUSIONS We generated a highly curated metabolic network of the human adipose tissue and used it to identify potential targets for adipose tissue metabolic dysfunction leading to the development of type 2 diabetes.
Collapse
Affiliation(s)
- Thierry Chénard
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, Canada
| | - Frédéric Guénard
- Institute of Nutrition and Functional Foods, Université Laval, Quebec City, Canada
| | - Marie-Claude Vohl
- Institute of Nutrition and Functional Foods, Université Laval, Quebec City, Canada.,School of Nutrition, Université Laval, Quebec City, Canada
| | - André Carpentier
- Division of Endocrinology, Department of Medicine, Centre de recherche du CHUS, Université de Sherbrooke, Sherbrooke, Canada
| | - André Tchernof
- School of Nutrition, Université Laval, Quebec City, Canada.,Centre de Recherche de l'Institut universitaire de cardiologie et de pneumologie de Québec, Quebec City, QC, Canada
| | - Rafael J Najmanovich
- Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Montreal, QC, Canada.
| |
Collapse
|
7
|
Chartier M, Morency LP, Zylber MI, Najmanovich RJ. Large-scale detection of drug off-targets: hypotheses for drug repurposing and understanding side-effects. BMC Pharmacol Toxicol 2017; 18:18. [PMID: 28449705 PMCID: PMC5408384 DOI: 10.1186/s40360-017-0128-7] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2016] [Accepted: 02/28/2017] [Indexed: 01/21/2023] Open
Abstract
Background Promiscuity in molecular interactions between small-molecules, including drugs, and proteins is widespread. Such unintended interactions can be exploited to suggest drug repurposing possibilities as well as to identify potential molecular mechanisms responsible for observed side-effects. Methods We perform a large-scale analysis to detect binding-site molecular interaction field similarities between the binding-sites of the primary target of 400 drugs against a dataset of 14082 cavities within 7895 different proteins representing a non-redundant dataset of all proteins with known structure. Statistically-significant cases with high levels of similarities represent potential cases where the drugs that bind the original target may in principle bind the suggested off-target. Such cases are further analysed with docking simulations to verify if indeed the drug could, in principle, bind the off-target. Diverse sources of data are integrated to associated potential cross-reactivity targets with side-effects. Results We observe that promiscuous binding-sites tend to display higher levels of hydrophobic and aromatic similarities. Focusing on the most statistically significant similarities (Z-score ≥ 3.0) and corroborating docking results (RMSD < 2.0 Å), we find 2923 cases involving 140 unique drugs and 1216 unique potential cross-reactivity protein targets. We highlight a few cases with a potential for drug repurposing (acetazolamide as a chorismate pyruvate lyase inhibitor, raloxifene as a bacterial quorum sensing inhibitor) as well as to explain the side-effects of zanamivir and captopril. A web-interface permits to explore the detected similarities for each of the 400 binding-sites of the primary drug targets and visualise them for the most statistically significant cases. Conclusions The detection of molecular interaction field similarities provide the opportunity to suggest drug repurposing opportunities as well as to identify potential molecular mechanisms responsible for side-effects. All methods utilized are freely available and can be readily applied to new query binding-sites. All data is freely available and represents an invaluable source to identify further candidates for repurposing and suggest potential mechanisms responsible for side-effects. Electronic supplementary material The online version of this article (doi:10.1186/s40360-017-0128-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Matthieu Chartier
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Québec, Canada
| | - Louis-Philippe Morency
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Québec, Canada
| | - María Inés Zylber
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Québec, Canada.,Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Québec, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Québec, Canada. .,Department of Pharmacology and Physiology, Faculty of Medicine, Université de Montréal, Québec, Canada.
| |
Collapse
|
8
|
Gaudreault F, Morency LP, Najmanovich RJ. NRGsuite: a PyMOL plugin to perform docking simulations in real time using FlexAID. Bioinformatics 2015; 31:3856-8. [PMID: 26249810 PMCID: PMC4653388 DOI: 10.1093/bioinformatics/btv458] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Accepted: 08/02/2015] [Indexed: 11/13/2022] Open
Abstract
Ligand protein docking simulations play a fundamental role in understanding molecular recognition. Herein we introduce the NRGsuite, a PyMOL plugin that permits the detection of surface cavities in proteins, their refinements, calculation of volume and use, individually or jointly, as target binding-sites for docking simulations with FlexAID. The NRGsuite offers the users control over a large number of important parameters in docking simulations including the assignment of flexible side-chains and definition of geometric constraints. Furthermore, the NRGsuite permits the visualization of the docking simulation in real time. The NRGsuite give access to powerful docking simulations that can be used in structure-guided drug design as well as an educational tool. The NRGsuite is implemented in Python and C/C++ with an easy to use package installer. The NRGsuite is available for Windows, Linux and MacOS. Availability and implementation: http://bcb.med.usherbrooke.ca/flexaid. Contact:rafael.najmanovich@usherbroke.ca Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Francis Gaudreault
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Canada
| | - Louis-Philippe Morency
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Canada
| |
Collapse
|
9
|
Frappier V, Chartier M, Najmanovich RJ. ENCoM server: exploring protein conformational space and the effect of mutations on protein function and stability. Nucleic Acids Res 2015; 43:W395-400. [PMID: 25883149 PMCID: PMC4489264 DOI: 10.1093/nar/gkv343] [Citation(s) in RCA: 118] [Impact Index Per Article: 13.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2015] [Revised: 03/27/2015] [Accepted: 04/06/2015] [Indexed: 11/24/2022] Open
Abstract
ENCoM is a coarse-grained normal mode analysis method recently introduced that unlike previous such methods is unique in that it accounts for the nature of amino acids. The inclusion of this layer of information was shown to improve conformational space sampling and apply for the first time a coarse-grained normal mode analysis method to predict the effect of single point mutations on protein dynamics and thermostability resulting from vibrational entropy changes. Here we present a web server that allows non-technical users to have access to ENCoM calculations to predict the effect of mutations on thermostability and dynamics as well as to generate geometrically realistic conformational ensembles. The server is accessible at: http://bcb.med.usherbrooke.ca/encom.
Collapse
Affiliation(s)
- Vincent Frappier
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Quebec, J1H 5N4, Canada
| | - Matthieu Chartier
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Quebec, J1H 5N4, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Quebec, J1H 5N4, Canada
| |
Collapse
|
10
|
Abstract
Small-molecule protein docking is an essential tool in drug design and to understand molecular recognition. In the present work we introduce FlexAID, a small-molecule docking algorithm that accounts for target side-chain flexibility and utilizes a soft scoring function, i.e. one that is not highly dependent on specific geometric criteria, based on surface complementarity. The pairwise energy parameters were derived from a large dataset of true positive poses and negative decoys from the PDBbind database through an iterative process using Monte Carlo simulations. The prediction of binding poses is tested using the widely used Astex dataset as well as the HAP2 dataset, while performance in virtual screening is evaluated using a subset of the DUD dataset. We compare FlexAID to AutoDock Vina, FlexX, and rDock in an extensive number of scenarios to understand the strengths and limitations of the different programs as well as to reported results for Glide, GOLD, and DOCK6 where applicable. The most relevant among these scenarios is that of docking on flexible non-native-complex structures where as is the case in reality, the target conformation in the bound form is not known a priori. We demonstrate that FlexAID, unlike other programs, is robust against increasing structural variability. FlexAID obtains equivalent sampling success as GOLD and performs better than AutoDock Vina or FlexX in all scenarios against non-native-complex structures. FlexAID is better than rDock when there is at least one critical side-chain movement required upon ligand binding. In virtual screening, FlexAID results are lower on average than those of AutoDock Vina and rDock. The higher accuracy in flexible targets where critical movements are required, intuitive PyMOL-integrated graphical user interface and free source code as well as precompiled executables for Windows, Linux, and Mac OS make FlexAID a welcome addition to the arsenal of existing small-molecule protein docking methods.
Collapse
Affiliation(s)
- Francis Gaudreault
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Quebec J1H5N4, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, Quebec J1H5N4, Canada
| |
Collapse
|
11
|
Abstract
Motivation: The field of structural bioinformatics and computational biophysics has undergone a revolution in the last 10 years. Developments that are captured annually through the 3DSIG meeting, upon which this article reflects. Results: An increase in the accessible data, computational resources and methodology has resulted in an increase in the size and resolution of studied systems and the complexity of the questions amenable to research. Concomitantly, the parameterization and efficiency of the methods have markedly improved along with their cross-validation with other computational and experimental results. Conclusion: The field exhibits an ever-increasing integration with biochemistry, biophysics and other disciplines. In this article, we discuss recent achievements along with current challenges within the field. Contact:Rafael.Najmanovich@USherbrooke.ca
Collapse
Affiliation(s)
- Ilan Samish
- Department of Plant Sciences, Weizmann Institute of Science, Rehovot, 76100, Israel, Ort Braude College, Karmiel, 2161002, Israel, Office of the Director, National Institutes of Health, Bethesda, MD 20814, USA and Department of Biochemistry, University of Sherbrooke, Sherbrooke, J1H 5N4, Canada Department of Plant Sciences, Weizmann Institute of Science, Rehovot, 76100, Israel, Ort Braude College, Karmiel, 2161002, Israel, Office of the Director, National Institutes of Health, Bethesda, MD 20814, USA and Department of Biochemistry, University of Sherbrooke, Sherbrooke, J1H 5N4, Canada
| | - Philip E Bourne
- Department of Plant Sciences, Weizmann Institute of Science, Rehovot, 76100, Israel, Ort Braude College, Karmiel, 2161002, Israel, Office of the Director, National Institutes of Health, Bethesda, MD 20814, USA and Department of Biochemistry, University of Sherbrooke, Sherbrooke, J1H 5N4, Canada
| | - Rafael J Najmanovich
- Department of Plant Sciences, Weizmann Institute of Science, Rehovot, 76100, Israel, Ort Braude College, Karmiel, 2161002, Israel, Office of the Director, National Institutes of Health, Bethesda, MD 20814, USA and Department of Biochemistry, University of Sherbrooke, Sherbrooke, J1H 5N4, Canada
| |
Collapse
|
12
|
Frappier V, Najmanovich RJ. A coarse-grained elastic network atom contact model and its use in the simulation of protein dynamics and the prediction of the effect of mutations. PLoS Comput Biol 2014; 10:e1003569. [PMID: 24762569 PMCID: PMC3998880 DOI: 10.1371/journal.pcbi.1003569] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Accepted: 02/25/2014] [Indexed: 11/18/2022] Open
Abstract
Normal mode analysis (NMA) methods are widely used to study dynamic aspects of protein structures. Two critical components of NMA methods are coarse-graining in the level of simplification used to represent protein structures and the choice of potential energy functional form. There is a trade-off between speed and accuracy in different choices. In one extreme one finds accurate but slow molecular-dynamics based methods with all-atom representations and detailed atom potentials. On the other extreme, fast elastic network model (ENM) methods with Cα-only representations and simplified potentials that based on geometry alone, thus oblivious to protein sequence. Here we present ENCoM, an Elastic Network Contact Model that employs a potential energy function that includes a pairwise atom-type non-bonded interaction term and thus makes it possible to consider the effect of the specific nature of amino-acids on dynamics within the context of NMA. ENCoM is as fast as existing ENM methods and outperforms such methods in the generation of conformational ensembles. Here we introduce a new application for NMA methods with the use of ENCoM in the prediction of the effect of mutations on protein stability. While existing methods are based on machine learning or enthalpic considerations, the use of ENCoM, based on vibrational normal modes, is based on entropic considerations. This represents a novel area of application for NMA methods and a novel approach for the prediction of the effect of mutations. We compare ENCoM to a large number of methods in terms of accuracy and self-consistency. We show that the accuracy of ENCoM is comparable to that of the best existing methods. We show that existing methods are biased towards the prediction of destabilizing mutations and that ENCoM is less biased at predicting stabilizing mutations.
Collapse
Affiliation(s)
- Vincent Frappier
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| | - Rafael J Najmanovich
- Department of Biochemistry, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| |
Collapse
|
13
|
Glouzon JPS, Bolduc F, Wang S, Najmanovich RJ, Perreault JP. Deep-sequencing of the peach latent mosaic viroid reveals new aspects of population heterogeneity. PLoS One 2014; 9:e87297. [PMID: 24498066 PMCID: PMC3907566 DOI: 10.1371/journal.pone.0087297] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Accepted: 12/24/2013] [Indexed: 01/04/2023] Open
Abstract
Viroids are small circular single-stranded infectious RNAs characterized by a relatively high mutation level. Knowledge of their sequence heterogeneity remains largely elusive and previous studies, using Sanger sequencing, were based on a limited number of sequences. In an attempt to address sequence heterogeneity from a population dynamics perspective, a GF305-indicator peach tree was infected with a single variant of the Avsunviroidae family member Peach latent mosaic viroid (PLMVd). Six months post-inoculation, full-length circular conformers of PLMVd were isolated and deep-sequenced. We devised an original approach to the bioinformatics refinement of our sequence libraries involving important phenotypic data, based on the systematic analysis of hammerhead self-cleavage activity. Two distinct libraries yielded a total of 3,939 different PLMVd variants. Sequence variants exhibiting up to ∼17% of mutations relative to the inoculated viroid were retrieved, clearly illustrating the high level of divergence dynamics within a unique population. While we initially assumed that most positions of the viroid sequence would mutate, we were surprised to discover that ∼50% of positions remained perfectly conserved, including several small stretches as well as a small motif reminiscent of a GNRA tetraloop which are the result of various selective pressures. Using a hierarchical clustering algorithm, the different variants harvested were subdivided into 7 clusters. We found that most sequences contained an average of 4.6 to 6.4 mutations compared to the variant used to initially inoculate the plant. Interestingly, it was possible to reconstitute and compare the sequence evolution of each of these clusters. In doing so, we identified several key mutations. This study provides a reliable pipeline for the treatment of viroid deep-sequencing. It also sheds new light on the extent of sequence variation that a viroid population can sustain, and which may give rise to a quasispecies.
Collapse
Affiliation(s)
- Jean-Pierre Sehi Glouzon
- Département d’informatique, Faculté des Sciences, Université de Sherbrooke, Sherbrooke, Québec, Canada
- Département de biochimie, Faculté de médecine et des sciences de la santé, Pavillon de Recherche Appliquée au Cancer, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - François Bolduc
- Département de biochimie, Faculté de médecine et des sciences de la santé, Pavillon de Recherche Appliquée au Cancer, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Shengrui Wang
- Département d’informatique, Faculté des Sciences, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Rafael J. Najmanovich
- Département de biochimie, Faculté de médecine et des sciences de la santé, Pavillon de Recherche Appliquée au Cancer, Université de Sherbrooke, Sherbrooke, Québec, Canada
- * E-mail: (RJN); (JPP)
| | - Jean-Pierre Perreault
- Département de biochimie, Faculté de médecine et des sciences de la santé, Pavillon de Recherche Appliquée au Cancer, Université de Sherbrooke, Sherbrooke, Québec, Canada
- * E-mail: (RJN); (JPP)
| |
Collapse
|
14
|
Han GW, Bakolitsa C, Miller MD, Kumar A, Carlton D, Najmanovich RJ, Abdubek P, Astakhova T, Axelrod HL, Chen C, Chiu HJ, Clayton T, Das D, Deller MC, Duan L, Ernst D, Feuerhelm J, Grant JC, Grzechnik A, Jaroszewski L, Jin KK, Johnson HA, Klock HE, Knuth MW, Kozbial P, Krishna SS, Marciano D, McMullan D, Morse AT, Nigoghossian E, Okach L, Reyes R, Rife CL, Sefcovic N, Tien HJ, Trame CB, van den Bedem H, Weekes D, Xu Q, Hodgson KO, Wooley J, Elsliger MA, Deacon AM, Godzik A, Lesley SA, Wilson IA. Structures of the first representatives of Pfam family PF06938 (DUF1285) reveal a new fold with repeated structural motifs and possible involvement in signal transduction. Acta Crystallogr Sect F Struct Biol Cryst Commun 2010; 66:1218-25. [PMID: 20944214 PMCID: PMC2954208 DOI: 10.1107/s1744309109050416] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2009] [Accepted: 11/23/2009] [Indexed: 11/14/2022]
Abstract
The crystal structures of SPO0140 and Sbal_2486 were determined using the semiautomated high-throughput pipeline of the Joint Center for Structural Genomics (JCSG) as part of the NIGMS Protein Structure Initiative (PSI). The structures revealed a conserved core with domain duplication and a superficial similarity of the C-terminal domain to pleckstrin homology-like folds. The conservation of the domain interface indicates a potential binding site that is likely to involve a nucleotide-based ligand, with genome-context and gene-fusion analyses additionally supporting a role for this family in signal transduction, possibly during oxidative stress.
Collapse
Affiliation(s)
- Gye Won Han
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Constantina Bakolitsa
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
| | - Mitchell D. Miller
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Abhinav Kumar
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Dennis Carlton
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Rafael J. Najmanovich
- Département de Biochimie, Faculté de Médecine, Université de Sherbrooke, Sherbrooke (Quebec), Canada
| | - Polat Abdubek
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Tamara Astakhova
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Herbert L. Axelrod
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Connie Chen
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Hsiu-Ju Chiu
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Thomas Clayton
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Debanu Das
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Marc C. Deller
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Lian Duan
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Dustin Ernst
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Julie Feuerhelm
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Joanna C. Grant
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Anna Grzechnik
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Lukasz Jaroszewski
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Kevin K. Jin
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Hope A. Johnson
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Heath E. Klock
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Mark W. Knuth
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Piotr Kozbial
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
| | - S. Sri Krishna
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - David Marciano
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Daniel McMullan
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Andrew T. Morse
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Edward Nigoghossian
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Linda Okach
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Ron Reyes
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Christopher L. Rife
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Natasha Sefcovic
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
| | - Henry J. Tien
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Christine B. Trame
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Henry van den Bedem
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Dana Weekes
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
| | - Qingping Xu
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Keith O. Hodgson
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Photon Science, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - John Wooley
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Marc-André Elsliger
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| | - Ashley M. Deacon
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Adam Godzik
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Program on Bioinformatics and Systems Biology, Burnham Institute for Medical Research, La Jolla, CA, USA
- Center for Research in Biological Systems, University of California, San Diego, La Jolla, CA, USA
| | - Scott A. Lesley
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
- Protein Sciences Department, Genomics Institute of the Novartis Research Foundation, San Diego, CA, USA
| | - Ian A. Wilson
- Joint Center for Structural Genomics, http://www.jcsg.org, USA
- Department of Molecular Biology, The Scripps Research Institute, La Jolla, CA, USA
| |
Collapse
|
15
|
Andreini C, Bertini I, Cavallaro G, Najmanovich RJ, Thornton JM. Structural analysis of metal sites in proteins: non-heme iron sites as a case study. J Mol Biol 2009; 388:356-80. [PMID: 19265704 DOI: 10.1016/j.jmb.2009.02.052] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2008] [Revised: 02/19/2009] [Accepted: 02/19/2009] [Indexed: 11/24/2022]
Abstract
In metalloproteins, the protein environment modulates metal properties to achieve the required goal, which can be protein stabilization or function. The analysis of metal sites at the atomic level of detail provided by protein structures can thus be of benefit in functional and evolutionary studies of proteins. In this work, we propose a structural bioinformatics approach to the study of metalloproteins based on structural templates of metal sites that include the PDB coordinates of protein residues forming the first and the second coordination sphere of the metal. We have applied this approach to non-heme iron sites, which have been analyzed at various levels. Templates of sites located in different protein domains have been compared, showing that similar sites can be found in unrelated proteins as the result of convergent evolution. Templates of sites located in proteins of a large superfamily have been compared, showing possible mechanisms of divergent evolution of proteins to achieve different functions. Furthermore, template comparisons have been used to predict the function of uncharacterized proteins, showing that similarity searches focused on metal sites can be advantageously combined with typical whole-domain comparisons. Structural templates of metal sites, finally, may constitute the basis for a systematic classification of metalloproteins in databases.
Collapse
Affiliation(s)
- Claudia Andreini
- Magnetic Resonance Center (CERM)-University of Florence, Via L. Sacconi 6, 50019 Sesto Fiorentino, Italy
| | | | | | | | | |
Collapse
|
16
|
Najmanovich RJ, Allali-Hassani A, Morris RJ, Dombrovsky L, Pan PW, Vedadi M, Plotnikov AN, Edwards A, Arrowsmith C, Thornton JM. Analysis of binding site similarity, small-molecule similarity and experimental binding profiles in the human cytosolic sulfotransferase family. Bioinformatics 2007; 23:e104-9. [PMID: 17237076 DOI: 10.1093/bioinformatics/btl292] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
MOTIVATION In the present work we combine computational analysis and experimental data to explore the extent to which binding site similarities between members of the human cytosolic sulfotransferase family correlate with small-molecule binding profiles. Conversely, from a small-molecule point of view, we explore the extent to which structural similarities between small molecules correlate to protein binding profiles. RESULTS The comparison of binding site structural similarities and small-molecule binding profiles shows that proteins with similar small-molecule binding profiles tend to have a higher degree of binding site similarity but the latter is not sufficient to predict small-molecule binding patterns, highlighting the difficulty of predicting small-molecule binding patterns from sequence or structure. Likewise, from a small-molecule perspective, small molecules with similar protein binding profiles tend to be topologically similar but topological similarity is not sufficient to predict their protein binding patterns. These observations have important consequences for function prediction and drug design.
Collapse
Affiliation(s)
- Rafael J Najmanovich
- European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus Cambridge CB10 1SD, UK.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
17
|
Abstract
The accurate identification of ligand binding sites in protein structures can be valuable in determining protein function. Once the binding site is known, it becomes easier to perform in silico and experimental procedures that may allow the ligand type and the protein function to be determined. For example, binding pocket shape analysis relies heavily on the correct localization of the ligand binding site. We have developed SURFNET-ConSurf, a modular, two-stage method for identifying the location and shape of potential ligand binding pockets in protein structures. In the first stage, the SURFNET program identifies clefts in the protein surface that are potential binding sites. In the second stage, these clefts are trimmed in size by cutting away regions distant from highly conserved residues, as defined by the ConSurf-HSSP database. The largest clefts that remain tend to be those where ligands bind. To test the approach, we analyzed a nonredundant set of 244 protein structures from the PDB and found that SURFNET-ConSurf identifies a ligand binding pocket in 75% of them. The trimming procedure reduces the original cleft volumes by 30% on average, while still encompassing an average 87% of the ligand volume. From the analysis of the results we conclude that for those cases in which the ligands are found in large, highly conserved clefts, the combined SURFNET-ConSurf method gives pockets that are a better match to the ligand shape and location. We also show that this approach works better for enzymes than for nonenzyme proteins.
Collapse
Affiliation(s)
- Fabian Glaser
- European Bioinformatics Institute, European Molecular Biology Laboratory, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom.
| | | | | | | | | |
Collapse
|
18
|
Najmanovich RJ, Torrance JW, Thornton JM. Prediction of protein function from structure: insights from methods for the detection of local structural similarities. Biotechniques 2005; 38:847, 849, 851. [PMID: 16018542 DOI: 10.2144/05386te01] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
19
|
Morris RJ, Najmanovich RJ, Kahraman A, Thornton JM. Real spherical harmonic expansion coefficients as 3D shape descriptors for protein binding pocket and ligand comparisons. Bioinformatics 2005; 21:2347-55. [PMID: 15728116 DOI: 10.1093/bioinformatics/bti337] [Citation(s) in RCA: 145] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION An increasing number of protein structures are being determined for which no biochemical characterization is available. The analysis of protein structure and function assignment is becoming an unexpected challenge and a major bottleneck towards the goal of well-annotated genomes. As shape plays a crucial role in biomolecular recognition and function, the examination and development of shape description and comparison techniques is likely to be of prime importance for understanding protein structure-function relationships. RESULTS A novel technique is presented for the comparison of protein binding pockets. The method uses the coefficients of a real spherical harmonics expansion to describe the shape of a protein's binding pocket. Shape similarity is computed as the L2 distance in coefficient space. Such comparisons in several thousands per second can be carried out on a standard linux PC. Other properties such as the electrostatic potential fit seamlessly into the same framework. The method can also be used directly for describing the shape of proteins and other molecules. AVAILABILITY A limited version of the software for the real spherical harmonics expansion of a set of points in PDB format is freely available upon request from the authors. Binding pocket comparisons and ligand prediction will be made available through the protein structure annotation pipeline Profunc (written by Roman Laskowski) which will be accessible from the EBI website shortly.
Collapse
Affiliation(s)
- Richard J Morris
- EMBL-EBI, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | | | | | | |
Collapse
|