1
|
In Silico Prediction of the Phosphorylation of NS3 as an Essential Mechanism for Dengue Virus Replication and the Antiviral Activity of Quercetin. BIOLOGY 2021; 10:biology10101067. [PMID: 34681164 PMCID: PMC8570334 DOI: 10.3390/biology10101067] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 10/09/2021] [Accepted: 10/11/2021] [Indexed: 11/25/2022]
Abstract
Simple Summary Dengue is a mosquito-borne virus that infects up to 400 million people worldwide annually. Dengue infection triggers high fever, severe body aches, rash, low platelet count, and could lead to Dengue hemorrhagic fever (DHF) in some cases. There is currently no cure, nor a broadly effective vaccine. The interaction of two viral proteins, nonstructural Proteins 3 and 5 (NS3 and NS5), is required for viral replication in the infected host’s cells. Our computational modeling of NS3 suggested that phosphorylation of a serine residue at position 137 of NS3 by a specific c-Jun N-terminal kinase (JNK) enhances viral replication by increasing the interaction of NS3 and NS5 through structural changes in amino acid residues 49–95. Experimental studies have shown that inhibition of JNK prevents viral replication and have suggested that the plants’ flavonoid Quercetin, Agathis flavone, and Myricetin inhibit Dengue infection. Our molecular simulations revealed that Quercetin binds NS3 and obstructs serine 137 phosphorylation, which may decrease viral replication. This work offers a molecular mechanism that can be used for anti-Dengue drug development. Abstract Dengue virus infection is a global health problem for which there have been challenges to obtaining a cure. Current vaccines and anti-viral drugs can only be narrowly applied in ongoing clinical trials. We employed computational methods based on structure-function relationships between human host kinases and viral nonstructural protein 3 (NS3) to understand viral replication inhibitors’ therapeutic effect. Phosphorylation at each of the two most evolutionarily conserved sites of NS3, serine 137 and threonine 189, compared to the unphosphorylated state were studied with molecular dynamics and docking simulations. The simulations suggested that phosphorylation at serine 137 caused a more remarkable structural change than phosphorylation at threonine 189, specifically located at amino acid residues 49–95. Docking studies supported the idea that phosphorylation at serine 137 increased the binding affinity between NS3 and nonstructural Protein 5 (NS5), whereas phosphorylation at threonine 189 decreased it. The interaction between NS3 and NS5 is essential for viral replication. Docking studies with the antiviral plant flavonoid Quercetin with NS3 indicated that Quercetin physically occluded the serine 137 phosphorylation site. Taken together, these findings suggested a specific site and mechanism by which Quercetin inhibits dengue and possible other flaviviruses.
Collapse
|
2
|
Ivanov A, Ramanathan P, Parry C, Ilinykh PA, Lin X, Petukhov M, Obukhov Y, Ammosova T, Amarasinghe GK, Bukreyev A, Nekhai S. Global phosphoproteomic analysis of Ebola virions reveals a novel role for VP35 phosphorylation-dependent regulation of genome transcription. Cell Mol Life Sci 2020; 77:2579-2603. [PMID: 31562565 PMCID: PMC7101265 DOI: 10.1007/s00018-019-03303-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2019] [Revised: 08/26/2019] [Accepted: 09/16/2019] [Indexed: 12/27/2022]
Abstract
Ebola virus (EBOV) causes severe human disease with a high case fatality rate. The balance of evidence implies that the virus circulates in bats. The molecular basis for host-viral interactions, including the role for phosphorylation during infections, is largely undescribed. To address this, and to better understand the biology of EBOV, the phosphorylation of EBOV proteins was analyzed in virions purified from infected monkey Vero-E6 cells and bat EpoNi/22.1 cells using high-resolution mass spectrometry. All EBOV structural proteins were detected with high coverage, along with phosphopeptides. Phosphorylation sites were identified in all viral structural proteins. Comparison of EBOV protein phosphorylation in monkey and bat cells showed only partial overlap of phosphorylation sites, with shared sites found in NP, VP35, and VP24 proteins, and no common sites in the other proteins. Three-dimensional structural models were built for NP, VP35, VP40, GP, VP30 and VP24 proteins using available crystal structures or by de novo structure prediction to elucidate the potential role of the phosphorylation sites. Phosphorylation of one of the identified sites in VP35, Thr-210, was demonstrated to govern the transcriptional activity of the EBOV polymerase complex. Thr-210 phosphorylation was also shown to be important for VP35 interaction with NP. This is the first study to compare phosphorylation of all EBOV virion proteins produced in primate versus bat cells, and to demonstrate the role of VP35 phosphorylation in the viral life cycle. The results uncover a novel mechanism of EBOV transcription and identify novel targets for antiviral drug development.
Collapse
Affiliation(s)
- Andrey Ivanov
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA
| | - Palaniappan Ramanathan
- Department of Pathology, University of Texas, Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77574-0609, USA
| | - Christian Parry
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA
- Department of Microbiology, Howard University, Washington, D.C., 20059, USA
| | - Philipp A Ilinykh
- Department of Pathology, University of Texas, Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77574-0609, USA
| | - Xionghao Lin
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA
- College of Dentistry, Howard University, Washington, D.C., 20059, USA
| | - Michael Petukhov
- Division of Molecular and Radiation Biophysics, Russian Nuclear Physics Institute Named After B. P. Konstantinov, National Research Center "Kurchatov Institute", Gatchina, 188300, Russia
- Russian Scientific Center of Radiology and Surgical Technologies Named After A. M. Granov, St. Petersburg, 197758, Russia
| | - Yuri Obukhov
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA
| | - Tatiana Ammosova
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA
- Department of Medicine, Howard University, Washington, D.C., 20059, USA
| | - Gaya K Amarasinghe
- Department of Pathology and Immunology, Washington University School of Medicine, St Louis, MO, 63110, USA
| | - Alexander Bukreyev
- Department of Pathology, University of Texas, Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77574-0609, USA.
- Department of Microbiology and Immunology, University of Texas, Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77574-0609, USA.
- Galveston National Laboratory, University of Texas, Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77574-0609, USA.
| | - Sergei Nekhai
- Center for Sickle Cell Disease, Howard University, 2201 Georgia Ave., N.W., Suite 321D, Washington, D.C., 20059, USA.
- Department of Microbiology, Howard University, Washington, D.C., 20059, USA.
- Department of Medicine, Howard University, Washington, D.C., 20059, USA.
| |
Collapse
|
3
|
Astl L, Verkhivker GM. Dynamic View of Allosteric Regulation in the Hsp70 Chaperones by J-Domain Cochaperone and Post-Translational Modifications: Computational Analysis of Hsp70 Mechanisms by Exploring Conformational Landscapes and Residue Interaction Networks. J Chem Inf Model 2020; 60:1614-1631. [DOI: 10.1021/acs.jcim.9b01045] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Lindy Astl
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States
| | - Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, California 92866, United States
- Depatment of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
4
|
Phosphorylated VP30 of Marburg Virus Is a Repressor of Transcription. J Virol 2018; 92:JVI.00426-18. [PMID: 30135121 DOI: 10.1128/jvi.00426-18] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Accepted: 08/06/2018] [Indexed: 12/29/2022] Open
Abstract
The filoviruses Marburg virus (MARV) and Ebola virus (EBOV) cause hemorrhagic fever in humans and nonhuman primates, with high case fatality rates. MARV VP30 is known to be phosphorylated and to interact with nucleoprotein (NP), but its role in regulation of viral transcription is disputed. Here, we analyzed phosphorylation of VP30 by mass spectrometry, which resulted in identification of multiple phosphorylated amino acids. Modeling the full-length three-dimensional structure of VP30 and mapping the identified phosphorylation sites showed that all sites lie in disordered regions, mostly in the N-terminal domain of the protein. Minigenome analysis of the identified phosphorylation sites demonstrated that phosphorylation of a cluster of amino acids at positions 46 through 53 inhibits transcription. To test the effect of VP30 phosphorylation on its interaction with other MARV proteins, coimmunoprecipitation analyses were performed. They demonstrated the involvement of VP30 phosphorylation in interaction with two other proteins of the MARV ribonucleoprotein complex, NP and VP35. To identify the role of protein phosphatase 1 (PP1) in the identified effects, a small molecule, 1E7-03, targeting a noncatalytic site of the enzyme that previously was shown to increase EBOV VP30 phosphorylation was used. Treatment of cells with 1E7-03 increased phosphorylation of VP30 at a cluster of phosphorylated amino acids from Ser-46 to Thr-53, reduced transcription of MARV minigenome, enhanced binding to NP and VP35, and dramatically reduced replication of infectious MARV particles. Thus, MARV VP30 phosphorylation can be targeted for development of future antivirals such as PP1-targeting compounds. IMPORTANCE The largest outbreak of MARV occurred in Angola in 2004 to 2005 and had a 90% case fatality rate. There are no approved treatments available for MARV. Development of antivirals as therapeutics requires a fundamental understanding of the viral life cycle. Because of the close similarity of MARV to another member of Filoviridae family, EBOV, it was assumed that the two viruses have similar mechanisms of regulation of transcription and replication. Here, characterization of the role of VP30 and its phosphorylation sites in transcription of the MARV genome demonstrated differences from those of EBOV. The identified phosphorylation sites appeared to inhibit transcription and appeared to be involved in interaction with both NP and VP35 ribonucleoproteins. A small molecule targeting PP1 inhibited transcription of the MARV genome, effectively suppressing replication of the viral particles. These data demonstrate the possibility developing antivirals based on compounds targeting PP1.
Collapse
|
5
|
Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics. Sci Rep 2018; 8:12976. [PMID: 30154518 PMCID: PMC6113251 DOI: 10.1038/s41598-018-31234-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Accepted: 08/13/2018] [Indexed: 12/14/2022] Open
Abstract
The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.
Collapse
|
6
|
Stetz G, Tse A, Verkhivker GM. Dissecting Structure-Encoded Determinants of Allosteric Cross-Talk between Post-Translational Modification Sites in the Hsp90 Chaperones. Sci Rep 2018; 8:6899. [PMID: 29720613 PMCID: PMC5932063 DOI: 10.1038/s41598-018-25329-4] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Accepted: 04/19/2018] [Indexed: 01/19/2023] Open
Abstract
Post-translational modifications (PTMs) represent an important regulatory instrument that modulates structure, dynamics and function of proteins. The large number of PTM sites in the Hsp90 proteins that are scattered throughout different domains indicated that synchronization of multiple PTMs through a combinatorial code can be invoked as an important mechanism to orchestrate diverse chaperone functions and recognize multiple client proteins. In this study, we have combined structural and coevolutionary analysis with molecular simulations and perturbation response scanning analysis of the Hsp90 structures to characterize functional role of PTM sites in allosteric regulation. The results reveal a small group of conserved PTMs that act as global mediators of collective dynamics and allosteric communications in the Hsp90 structures, while the majority of flexible PTM sites serve as sensors and carriers of the allosteric structural changes. This study provides a comprehensive structural, dynamic and network analysis of PTM sites across Hsp90 proteins, identifying specific role of regulatory PTM hotspots in the allosteric mechanism of the Hsp90 cycle. We argue that plasticity of a combinatorial PTM code in the Hsp90 may be enacted through allosteric coupling between effector and sensor PTM residues, which would allow for timely response to structural requirements of multiple modified enzymes.
Collapse
Affiliation(s)
- Gabrielle Stetz
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
| | - Amanda Tse
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California, United States of America.
- Chapman University School of Pharmacy, Irvine, California, United States of America.
| |
Collapse
|
7
|
Karasev DA, Veselova DA, Veselovsky AV, Sobolev BN, Zgoda VG, Archakov AI. Spatial features of proteins related to their phosphorylation and associated structural changes. Proteins 2017; 86:13-20. [DOI: 10.1002/prot.25397] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2017] [Revised: 09/13/2017] [Accepted: 10/04/2017] [Indexed: 11/06/2022]
Affiliation(s)
- Dmitry A. Karasev
- Department of Bioinformatics; Institute of Biomedical Chemistry (IBMC); Moscow Russia
- Department of Biochemistry; Pirogov Russian National Research Medical University (RNRMU); Moscow Russia
| | - Darya A. Veselova
- Department of Bioinformatics; Institute of Biomedical Chemistry (IBMC); Moscow Russia
| | | | - Boris N. Sobolev
- Department of Bioinformatics; Institute of Biomedical Chemistry (IBMC); Moscow Russia
| | - Victor G. Zgoda
- Department of Bioinformatics; Institute of Biomedical Chemistry (IBMC); Moscow Russia
| | - Alexander I. Archakov
- Department of Bioinformatics; Institute of Biomedical Chemistry (IBMC); Moscow Russia
| |
Collapse
|
8
|
Johnson JR, Santos SD, Johnson T, Pieper U, Strumillo M, Wagih O, Sali A, Krogan NJ, Beltrao P. Prediction of Functionally Important Phospho-Regulatory Events in Xenopus laevis Oocytes. PLoS Comput Biol 2015; 11:e1004362. [PMID: 26312481 PMCID: PMC4552029 DOI: 10.1371/journal.pcbi.1004362] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2014] [Accepted: 05/27/2015] [Indexed: 01/10/2023] Open
Abstract
The African clawed frog Xenopus laevis is an important model organism for studies in developmental and cell biology, including cell-signaling. However, our knowledge of X. laevis protein post-translational modifications remains scarce. Here, we used a mass spectrometry-based approach to survey the phosphoproteome of this species, compiling a list of 2636 phosphosites. We used structural information and phosphoproteomic data for 13 other species in order to predict functionally important phospho-regulatory events. We found that the degree of conservation of phosphosites across species is predictive of sites with known molecular function. In addition, we predicted kinase-protein interactions for a set of cell-cycle kinases across all species. The degree of conservation of kinase-protein interactions was found to be predictive of functionally relevant regulatory interactions. Finally, using comparative protein structure models, we find that phosphosites within structured domains tend to be located at positions with high conformational flexibility. Our analysis suggests that a small class of phosphosites occurs in positions that have the potential to regulate protein conformation.
Collapse
Affiliation(s)
- Jeffrey R Johnson
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, California, United States of America
| | - Silvia D Santos
- Quantitative Cell Biology group, MRC Clinical Sciences Centre, Imperial College, London, United Kingdom
| | - Tasha Johnson
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, California, United States of America
| | - Ursula Pieper
- Department of Bioengineering and Therapeutic Sciences, California Institute for Quantitative Biosciences, Byers Hall at Mission Bay, University of California, San Francisco, San Francisco, California, United States of America; Department of Pharmaceutical Chemistry, California Institute for Quantitative Biosciences, Byers Hall at Mission Bay, University of California, San Francisco, San Francisco, California, United States of America
| | - Marta Strumillo
- European Molecular Biology Laboratory, Genome Biology Unit, Heidelberg, Germany and European Bioinformatics Institute (EMBL-EBI), Cambridge, United Kingdom
| | - Omar Wagih
- European Molecular Biology Laboratory, Genome Biology Unit, Heidelberg, Germany and European Bioinformatics Institute (EMBL-EBI), Cambridge, United Kingdom
| | - Andrej Sali
- Department of Bioengineering and Therapeutic Sciences, California Institute for Quantitative Biosciences, Byers Hall at Mission Bay, University of California, San Francisco, San Francisco, California, United States of America; Department of Pharmaceutical Chemistry, California Institute for Quantitative Biosciences, Byers Hall at Mission Bay, University of California, San Francisco, San Francisco, California, United States of America
| | - Nevan J Krogan
- Department of Cellular and Molecular Pharmacology, University of California, San Francisco, San Francisco, California, United States of America; Gladstone Institutes, San Francisco, California, United States of America
| | - Pedro Beltrao
- European Molecular Biology Laboratory, Genome Biology Unit, Heidelberg, Germany and European Bioinformatics Institute (EMBL-EBI), Cambridge, United Kingdom; iBiMED and Department of Health Sciences, University of Aveiro, Aveiro, Portugal
| |
Collapse
|