51
|
Verkhivker GM, Di Paola L. Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies. J Phys Chem B 2021; 125:4596-4619. [PMID: 33929853 PMCID: PMC8098774 DOI: 10.1021/acs.jpcb.1c00395] [Citation(s) in RCA: 47] [Impact Index Per Article: 15.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Revised: 04/16/2021] [Indexed: 02/07/2023]
Abstract
Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck Center for Science and Engineering, Schmid
College of Science and Technology, Chapman University, One
University Drive, Orange, California 92866, United States
- Department of Biomedical and Pharmaceutical Sciences,
Chapman University School of Pharmacy, Irvine, California
92618, United States
| | - Luisa Di Paola
- Unit of Chemical-Physics Fundamentals in Chemical
Engineering, Department of Engineering, Università Campus Bio-Medico
di Roma, via Álvaro del Portillo 21, 00128 Rome,
Italy
| |
Collapse
|
52
|
Wang S, Ma C, Zeng A. Dynamic energy correlation analysis of E. coli aspartokinase III and alteration of allosteric regulation by manipulating energy transduction pathways. Eng Life Sci 2021; 21:314-323. [PMID: 33976604 PMCID: PMC8092979 DOI: 10.1002/elsc.202000065] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2020] [Revised: 11/20/2020] [Accepted: 01/13/2021] [Indexed: 01/14/2023] Open
Abstract
Conformational change associated with allosteric regulation in a protein is ultimately driven by energy transformation. However, little is known about the latter process. In this work, we combined steered molecular dynamics simulations and sequence conservation analysis to investigate the conformational changes and energy transformation in the allosteric enzyme aspartokinase III (AK III) from Escherichia coli. Correlation analysis of energy change at residue level indicated significant transformation between electrostatic energy and dihedral angle energy during the allosteric regulation. Key amino acid residues located in the corresponding energy transduction pathways were identified by dynamic energy correlation analysis. To verify their functions, residues with a high energy correlation in the pathways were altered and their effects on allosteric regulation of AKIII were determined. This study sheds new insights into energy transformation during allosteric regulation of AK III and proposes a strategy to identify key residues that are involved in intramolecular energy transduction and thus in driving the allosteric process.
Collapse
Affiliation(s)
- Shizhen Wang
- Department of Chemical and Biochemical EngineeringCollege of Chemistry and Chemical EngineeringXiamen UniversityXiamenP. R. China
- Institute of Bioprocess and Biosystems EngineeringHamburg University of TechnologyHamburgGermany
| | - Chengwei Ma
- Institute of Bioprocess and Biosystems EngineeringHamburg University of TechnologyHamburgGermany
| | - An‐Ping Zeng
- Institute of Bioprocess and Biosystems EngineeringHamburg University of TechnologyHamburgGermany
| |
Collapse
|
53
|
Prabantu VM, Naveenkumar N, Srinivasan N. Influence of Disease-Causing Mutations on Protein Structural Networks. Front Mol Biosci 2021; 7:620554. [PMID: 33778000 PMCID: PMC7987782 DOI: 10.3389/fmolb.2020.620554] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Accepted: 12/17/2020] [Indexed: 01/18/2023] Open
Abstract
The interactions between residues in a protein tertiary structure can be studied effectively using the approach of protein structure network (PSN). A PSN is a node-edge representation of the structure with nodes representing residues and interactions between residues represented by edges. In this study, we have employed weighted PSNs to understand the influence of disease-causing mutations on proteins of known 3D structures. We have used manually curated information on disease mutations from UniProtKB/Swiss-Prot and their corresponding protein structures of wildtype and disease variant from the protein data bank. The PSNs of the wildtype and disease-causing mutant are compared to analyse variation of global and local dissimilarity in the overall network and at specific sites. We study how a mutation at a given site can affect the structural network at a distant site which may be involved in the function of the protein. We have discussed specific examples of the disease cases where the protein structure undergoes limited structural divergence in their backbone but have large dissimilarity in their all atom networks and vice versa, wherein large conformational alterations are observed while retaining overall network. We analyse the effect of variation of network parameters that characterize alteration of function or stability.
Collapse
Affiliation(s)
| | - Nagarajan Naveenkumar
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India.,National Centre for Biological Sciences, TIFR, Bangalore, India.,Bharathidasan University, Tiruchirappalli, India
| | | |
Collapse
|
54
|
Dodd T, Yao XQ, Hamelberg D, Ivanov I. Subsets of adjacent nodes (SOAN): a fast method for computing suboptimal paths in protein dynamic networks. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1893847] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Thomas Dodd
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Ivaylo Ivanov
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| |
Collapse
|
55
|
Abstract
Allosteric regulation in proteins is fundamental to many important biological processes. Allostery has been employed to control protein functions by regulating protein activity. Engineered allosteric regulation allows controlling protein activity in subsecond time scale and has a broad range of applications, from dissecting spatiotemporal dynamics in biochemical cascades to applications in biotechnology and medicine. Here, we review the concept of allostery in proteins and various approaches to identify allosteric sites and pathways. We then provide an overview of strategies and tools used in allosteric protein regulation and their utility in biological applications. We highlight various classes of proteins, where regulation is achieved through allostery. Finally, we analyze the current problems, critical challenges, and future prospective in achieving allosteric regulation in proteins.
Collapse
Affiliation(s)
| | - Jiaxing Chen
- Department of Pharmacology, Penn State College of Medicine, Hershey, Pennsylvania 17033-0850, United States
| | - Nikolay V Dokholyan
- Department of Pharmacology, Penn State College of Medicine, Hershey, Pennsylvania 17033-0850, United States
- Departments of Biochemistry & Molecular Biology, Penn State College of Medicine, Hershey, Pennsylvania 17033-0850, United States
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States
- Department of Biomedical Engineering, Penn State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
56
|
Bhattacharjee S, Bhattacharyya R, Sengupta J. Dynamics and electrostatics define an allosteric druggable site within the receptor-binding domain of SARS-CoV-2 spike protein. FEBS Lett 2021; 595:442-451. [PMID: 33449359 PMCID: PMC8014131 DOI: 10.1002/1873-3468.14038] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Revised: 12/11/2020] [Accepted: 01/07/2021] [Indexed: 02/07/2023]
Abstract
The pathogenesis of the SARS‐CoV‐2 virus initiates through recognition of the angiotensin‐converting enzyme 2 (ACE2) receptor of the host cells by the receptor‐binding domain (RBD) located at the spikes of the virus. Here, using molecular dynamics simulations, we have demonstrated the allosteric crosstalk within the RBD in the apo‐ and the ACE2 receptor‐bound states, revealing the contribution of the dynamics‐based correlated motions and the electrostatic energy perturbations to this crosstalk. While allostery, based on correlated motions, dominates inherent distal communication in the apo‐RBD, the electrostatic energy perturbations determine favorable pairwise crosstalk within the RBD residues upon binding to ACE2. Interestingly, the allosteric path is composed of residues which are evolutionarily conserved within closely related coronaviruses, pointing toward the biological relevance of the communication and its potential as a target for drug development.
Collapse
Affiliation(s)
- Sayan Bhattacharjee
- Division of Structural Biology and Bioinformatics, CSIR-Indian Institute of Chemical Biology, Kolkata, India
| | - Rajanya Bhattacharyya
- Division of Structural Biology and Bioinformatics, CSIR-Indian Institute of Chemical Biology, Kolkata, India
| | - Jayati Sengupta
- Division of Structural Biology and Bioinformatics, CSIR-Indian Institute of Chemical Biology, Kolkata, India.,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
| |
Collapse
|
57
|
Verkhivker GM, Di Paola L. Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches. J Phys Chem B 2021; 125:850-873. [PMID: 33448856 PMCID: PMC7839160 DOI: 10.1021/acs.jpcb.0c10637] [Citation(s) in RCA: 44] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Revised: 01/08/2021] [Indexed: 12/13/2022]
Abstract
The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Schmid College of Science and
Technology, Chapman University, One University Drive, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Luisa Di Paola
- Unit
of Chemical-Physics Fundamentals in Chemical Engineering, Department
of Engineering, Università Campus
Bio-Medico di Roma, via
Álvaro del Portillo 21, 00128 Rome, Italy
| |
Collapse
|
58
|
Sathanur AV, Baker NA. A clustering-based biased Monte Carlo approach to protein titration curve prediction. PROCEEDINGS OF THE ... INTERNATIONAL CONFERENCE ON MACHINE LEARNING AND APPLICATIONS. INTERNATIONAL CONFERENCE ON MACHINE LEARNING AND APPLICATIONS 2020; 2020:10.1109/icmla51294.2020.00037. [PMID: 34661203 PMCID: PMC8513769 DOI: 10.1109/icmla51294.2020.00037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In this work, we developed an efficient approach to compute ensemble averages in systems with pairwise-additive energetic interactions between the entities. Methods involving full enumeration of the configuration space result in exponential complexity. Sampling methods such as Markov Chain Monte Carlo (MCMC) algorithms have been proposed to tackle the exponential complexity of these problems; however, in certain scenarios where significant energetic coupling exists between the entities, the efficiency of the such algorithms can be diminished. We used a strategy to improve the efficiency of MCMC by taking advantage of the cluster structure in the interaction energy matrix to bias the sampling. We pursued two different schemes for the biased MCMC runs and show that they are valid MCMC schemes. We used both synthesized and real-world systems to show the improved performance of our biased MCMC methods when compared to the regular MCMC method. In particular, we applied these algorithms to the problem of estimating protonation ensemble averages and titration curves of residues in a protein.
Collapse
Affiliation(s)
| | - Nathan A Baker
- Pacific Northwest National Laboratory, Richland, WA, USA
| |
Collapse
|
59
|
Verkhivker GM. Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins. J Proteome Res 2020; 19:4587-4608. [PMID: 33006900 PMCID: PMC7640983 DOI: 10.1021/acs.jproteome.0c00654] [Citation(s) in RCA: 46] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2020] [Indexed: 12/13/2022]
Abstract
The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.
Collapse
Affiliation(s)
- 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
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
60
|
Villani G. A Time-Dependent Quantum Approach to Allostery and a Comparison With Light-Harvesting in Photosynthetic Phenomenon. Front Mol Biosci 2020; 7:156. [PMID: 33005625 PMCID: PMC7483663 DOI: 10.3389/fmolb.2020.00156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Accepted: 06/19/2020] [Indexed: 11/26/2022] Open
Abstract
The allosteric effect is one of the most important processes in regulating the function of proteins, and the elucidation of this phenomenon plays a significant role in understanding emergent behaviors in biological regulation. In this process, a perturbation, generated by a ligand in a part of the macromolecule (the allosteric site), moves along this system and reaches a specific (active) site, dozens of Ångströms away, with a great efficiency. The dynamics of this perturbation in the macromolecule can model precisely the allosteric process. In this article, we will be studying the general characteristics of allostery, using a time-dependent quantum approach to obtain rules that apply to this kind of process. Considering the perturbation as a wave that moves within the molecular system, we will characterize the allosteric process with three of the properties of this wave in the active site: (1) ta, the characteristic time for reaching that site, (2) Aa, the amplitude of the wave in this site, and (3) Ba, its corresponding spectral broadening. These three parameters, together with the process mechanism and the perturbation efficiency in the process, can describe the phenomenon. One of the main purposes of this paper is to link the parameters ta, Aa, and Ba and the perturbation efficiency to the characteristics of the system. There is another fundamental process for life that has some characteristics similar to allostery: the light-harvesting (LH) process in photosynthesis. Here, as in allostery, two distant macromolecular sites are involved—two sites dozens of Ångströms away. In both processes, it is particularly important that the perturbation is distributed efficiently without dissipating in the infinite degrees of freedom within the macromolecule. The importance of considering quantum effects in the LH process is well documented in literature, and the quantum coherences are experimentally proven by time-dependent spectroscopic techniques. Given the existing similarities between these two processes in macromolecules, in this work, we suggest using Quantum Mechanics (QM) to study allostery.
Collapse
Affiliation(s)
- Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici (UOS Pisa) - CNR, Area della Ricerca di Pisa, Pisa, Italy
| |
Collapse
|
61
|
Felline A, Seeber M, Fanelli F. webPSN v2.0: a webserver to infer fingerprints of structural communication in biomacromolecules. Nucleic Acids Res 2020; 48:W94-W103. [PMID: 32427333 PMCID: PMC7319592 DOI: 10.1093/nar/gkaa397] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2020] [Revised: 04/21/2020] [Accepted: 05/17/2020] [Indexed: 12/25/2022] Open
Abstract
A mixed Protein Structure Network (PSN) and Elastic Network Model-Normal Mode Analysis (ENM-NMA)-based strategy (i.e. PSN-ENM) was developed to investigate structural communication in bio-macromolecules. Protein Structure Graphs (PSGs) are computed on a single structure, whereas information on system dynamics is supplied by ENM-NMA. The approach was implemented in a webserver (webPSN), which was significantly updated herein. The webserver now handles both proteins and nucleic acids and relies on an internal upgradable database of network parameters for ions and small molecules in all PDB structures. Apart from the radical restyle of the server and some changes in the calculation setup, other major novelties concern the possibility to: a) compute the differences in nodes, links, and communication pathways between two structures (i.e. network difference) and b) infer links, hubs, communities, and metapaths from consensus networks computed on a number of structures. These new features are useful to identify commonalties and differences between two different functional states of the same system or structural-communication signatures in homologous or analogous systems. The output analysis relies on 3D-representations, interactive tables and graphs, also available for download. Speed and accuracy make this server suitable to comparatively investigate structural communication in large sets of bio-macromolecular systems. URL: http://webpsn.hpc.unimore.it.
Collapse
Affiliation(s)
- Angelo Felline
- Department of Life Sciences, University of Modena and Reggio Emilia, Modena 41125, Italy
| | - Michele Seeber
- Department of Life Sciences, University of Modena and Reggio Emilia, Modena 41125, Italy
| | - Francesca Fanelli
- Department of Life Sciences, University of Modena and Reggio Emilia, Modena 41125, Italy.,Center for Neuroscience and Neurotechnology, University of Modena and Reggio Emilia, Modena 41125, Italy
| |
Collapse
|
62
|
Astl L, Stetz G, Verkhivker GM. Dissecting Molecular Principles of the Hsp90 Chaperone Regulation by Allosteric Modulators Using a Hierarchical Simulation Approach and Network Modeling of Allosteric Interactions: Conformational Selection Dictates the Diversity of Protein Responses and Ligand-Specific Functional Mechanisms. J Chem Theory Comput 2020; 16:6656-6677. [PMID: 32941034 DOI: 10.1021/acs.jctc.0c00503] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Conformational plasticity of the Hsp90 molecular chaperones underlies the diversity of functional mechanisms that these versatile molecular machines employ to coordinate their vast protein clientele in the cellular environment. Despite a steady progress in studies of the Hsp90 machinery, a great deal remains unknown about molecular principles and ligand-specific functional mechanisms of the Hsp90 regulation by allosteric modulators that attracted significant attention because of their therapeutic potential. Due to structural complexity and dynamic nature of the Hsp90 responses to allosteric modulators, the atomistic details about the mode of action of these small molecules continue to be fairly scarce and controversial. In this work, we employ an integrative strategy that encompassed atomistic simulations of the Hsp90 proteins and hierarchical modeling of Hsp90-ligand binding with network analysis to explore functional mechanisms of the Hsp90 regulation by a panel of allosteric modulators (novobiocin, KU-135, KU-174, and KU-32) with different models of action. The results show that functional mechanisms of allosteric modulation in the Hsp90 proteins may be driven by conformational selection principles in which ligands elicit pre-existing states of the unbound chaperone to drive ligand-specific protein responses and distinct scenarios of Hsp90 regulation. We found that novobiocin can selectively sequester an ensemble of open chaperone conformations and inhibit the progression of the functional cycle through a cascade of cumulative dynamic changes. In contrast, KU-32 displayed unique preferences toward partially closed dynamic states, inducing robust allosteric signaling and stimulation of the ATPase cycle. The proposed model of the Hsp90 regulation by allosteric modulators reconciled diverse experimental data and showed that allosteric modulators may operate via targeted exploitation of dynamic landscapes eliciting vastly different protein responses and diverse mechanisms of action.
Collapse
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
| | - Gabrielle Stetz
- 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.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
63
|
Bhardwaj VK, Purohit R. Structural changes induced by substitution of amino acid 129 in the coat protein of Cucumber mosaic virus. Genomics 2020; 112:3729-3738. [DOI: 10.1016/j.ygeno.2020.04.023] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Revised: 04/01/2020] [Accepted: 04/24/2020] [Indexed: 01/06/2023]
|
64
|
Wang J, Jain A, McDonald LR, Gambogi C, Lee AL, Dokholyan NV. Mapping allosteric communications within individual proteins. Nat Commun 2020; 11:3862. [PMID: 32737291 PMCID: PMC7395124 DOI: 10.1038/s41467-020-17618-2] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2019] [Accepted: 06/30/2020] [Indexed: 02/05/2023] Open
Abstract
Allostery in proteins influences various biological processes such as regulation of gene transcription and activities of enzymes and cell signaling. Computational approaches for analysis of allosteric coupling provide inexpensive opportunities to predict mutations and to design small-molecule agents to control protein function and cellular activity. We develop a computationally efficient network-based method, Ohm, to identify and characterize allosteric communication networks within proteins. Unlike previously developed simulation-based approaches, Ohm relies solely on the structure of the protein of interest. We use Ohm to map allosteric networks in a dataset composed of 20 proteins experimentally identified to be allosterically regulated. Further, the Ohm allostery prediction for the protein CheY correlates well with NMR CHESCA studies. Our webserver, Ohm.dokhlab.org, automatically determines allosteric network architecture and identifies critical coupled residues within this network.
Collapse
Affiliation(s)
- Jian Wang
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA, 17033-0850, USA
| | - Abha Jain
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599-7363, USA
| | - Leanna R McDonald
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599-7363, USA
| | - Craig Gambogi
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599-7363, USA
| | - Andrew L Lee
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599-7363, USA
| | - Nikolay V Dokholyan
- Department of Pharmacology, Penn State College of Medicine, Hershey, PA, 17033-0850, USA.
- Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC, 27599, USA.
- Departments of Biochemistry & Molecular Biology, Penn State College of Medicine, Hershey, PA, 17033-0850, USA.
| |
Collapse
|
65
|
Verkhivker GM, Agajanian S, Hu G, Tao P. Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning. Front Mol Biosci 2020; 7:136. [PMID: 32733918 PMCID: PMC7363947 DOI: 10.3389/fmolb.2020.00136] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Accepted: 06/08/2020] [Indexed: 12/12/2022] Open
Abstract
Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.
Collapse
Affiliation(s)
- Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, United States
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, United States
| | - Steve Agajanian
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, United States
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Peng Tao
- Department of Chemistry, Center for Drug Discovery, Design, and Delivery (CD4), Center for Scientific Computation, Southern Methodist University, Dallas, TX, United States
| |
Collapse
|
66
|
Miotto M, Olimpieri PP, Di Rienzo L, Ambrosetti F, Corsi P, Lepore R, Tartaglia GG, Milanetti E. Insights on protein thermal stability: a graph representation of molecular interactions. Bioinformatics 2020; 35:2569-2577. [PMID: 30535291 PMCID: PMC6662296 DOI: 10.1093/bioinformatics/bty1011] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Revised: 10/29/2018] [Accepted: 12/07/2018] [Indexed: 11/14/2022] Open
Abstract
Motivation Understanding the molecular mechanisms of thermal stability is a challenge in protein biology. Indeed, knowing the temperature at which proteins are stable has important theoretical implications, which are intimately linked with properties of the native fold, and a wide range of potential applications from drug design to the optimization of enzyme activity. Results Here, we present a novel graph-theoretical framework to assess thermal stability based on the structure without any a priori information. In this approach we describe proteins as energy-weighted graphs and compare them using ensembles of interaction networks. Investigating the position of specific interactions within the 3D native structure, we developed a parameter-free network descriptor that permits to distinguish thermostable and mesostable proteins with an accuracy of 76% and area under the receiver operating characteristic curve of 78%. Availability and implementation Code is available upon request to edoardo.milanetti@uniroma1.it Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mattia Miotto
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Center for Life Nano Science@Sapienza, Instituto Italiano di Tecnologia, Viale Regina Elena, 291 Roma (RM), Italy.,Soft and Living Matter Laboratory, Institute of Nanotechnology, Consiglio Nazionale delle Ricerche, Rome, Italy
| | | | - Lorenzo Di Rienzo
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy
| | - Francesco Ambrosetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Padualaan 8, Utrecht, the Netherlands
| | - Pietro Corsi
- Department of Science, Università degli Studi "Roma Tre", via della Vasca Navale 84, Rome, Italy
| | - Rosalba Lepore
- Biozentrum, University of Basel, Klingelbergstrasse 50-70, CH-4056 Basel, Switzerland.,SIB Swiss Institute of Bioinformatics, Biozentrum, University of Basel, Klingelbergstrasse 50-70, CH-4056 Basel, Switzerland
| | - Gian Gaetano Tartaglia
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader St. 88, Barcelona, Spain.,Institucio' Catalana de Recerca i Estudis Avancats (ICREA), 23 Passeig Lluìs Companys, Barcelona, Spain.,Department of Biology and Biotechnology, Sapienza University of Rome, Piazzale Aldo Moro 5, Rome, Italy
| | - Edoardo Milanetti
- Department of Physics, Sapienza University, Piazzale Aldo Moro 5, Rome, Italy.,Center for Life Nano Science@Sapienza, Instituto Italiano di Tecnologia, Viale Regina Elena, 291 Roma (RM), Italy
| |
Collapse
|
67
|
Stetz G, Astl L, Verkhivker GM. Exploring Mechanisms of Communication Switching in the Hsp90-Cdc37 Regulatory Complexes with Client Kinases through Allosteric Coupling of Phosphorylation Sites: Perturbation-Based Modeling and Hierarchical Community Analysis of Residue Interaction Networks. J Chem Theory Comput 2020; 16:4706-4725. [PMID: 32492340 DOI: 10.1021/acs.jctc.0c00280] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Understanding molecular principles underlying chaperone-based modulation of kinase client activity is critically important to dissect functions and activation mechanisms of many oncogenic proteins. The recent experimental studies have suggested that phosphorylation sites in the Hsp90 and Cdc37 proteins can serve as conformational communication switches of chaperone regulation and kinase interactions. However, a mechanism of allosteric coupling between phosphorylation sites in the Hsp90 and Cdc37 during client binding is poorly understood, and the molecular signatures underpinning specific roles of phosphorylation sites in the Hsp90 regulation remain unknown. In this work, we employed a combination of evolutionary analysis, coarse-grained molecular simulations together with perturbation-based network modeling and scanning of the unbound and bound Hsp90 and Cdc37 structures to quantify allosteric effects of phosphorylation sites and identify unique signatures that are characteristic for communication switches of kinase-specific client binding. By using network-based metrics of the dynamic intercommunity bridgeness and community centrality, we characterize specific signatures of phosphorylation switches involved in allosteric regulation. Through perturbation-based analysis of the dynamic residue interaction networks, we show that mutations of kinase-specific phosphorylation switches can induce long-range effects and lead to a global rewiring of the allosteric network and signal transmission in the Hsp90-Cdc37-kinase complex. We determine a specific group of phosphorylation sites in the Hsp90 where mutations may have a strong detrimental effect on allosteric interaction network, providing insight into the mechanism of phosphorylation-induced communication switching. The results demonstrate that kinase-specific phosphorylation switches of communications in the Hsp90 may be partly predisposed for their regulatory role based on preexisting allosteric propensities.
Collapse
Affiliation(s)
- Gabrielle Stetz
- 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
| | - 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.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
68
|
Astl L, Stetz G, Verkhivker GM. Allosteric Mechanism of the Hsp90 Chaperone Interactions with Cochaperones and Client Proteins by Modulating Communication Spines of Coupled Regulatory Switches: Integrative Atomistic Modeling of Hsp90 Signaling in Dynamic Interaction Networks. J Chem Inf Model 2020; 60:3616-3631. [DOI: 10.1021/acs.jcim.0c00380] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [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
| | - Gabrielle Stetz
- 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
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California92618, United States
| |
Collapse
|
69
|
Chakrabarty B, Naganathan V, Garg K, Agarwal Y, Parekh N. NAPS update: network analysis of molecular dynamics data and protein-nucleic acid complexes. Nucleic Acids Res 2020; 47:W462-W470. [PMID: 31106363 PMCID: PMC6602509 DOI: 10.1093/nar/gkz399] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Revised: 04/30/2019] [Accepted: 05/07/2019] [Indexed: 02/04/2023] Open
Abstract
Network theory is now a method of choice to gain insights in understanding protein structure, folding and function. In combination with molecular dynamics (MD) simulations, it is an invaluable tool with widespread applications such as analyzing subtle conformational changes and flexibility regions in proteins, dynamic correlation analysis across distant regions for allosteric communications, in drug design to reveal alternative binding pockets for drugs, etc. Updated version of NAPS now facilitates network analysis of the complete repertoire of these biomolecules, i.e., proteins, protein–protein/nucleic acid complexes, MD trajectories, and RNA. Various options provided for analysis of MD trajectories include individual network construction and analysis of intermediate time-steps, comparative analysis of these networks, construction and analysis of average network of the ensemble of trajectories and dynamic cross-correlations. For protein–nucleic acid complexes, networks of the whole complex as well as that of the interface can be constructed and analyzed. For analysis of proteins, protein–protein complexes and MD trajectories, network construction based on inter-residue interaction energies with realistic edge-weights obtained from standard force fields is provided to capture the atomistic details. Updated version of NAPS also provides improved visualization features, interactive plots and bulk execution. URL: http://bioinf.iiit.ac.in/NAPS/
Collapse
Affiliation(s)
- Broto Chakrabarty
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology - Hyderabad 500032, India
| | - Varun Naganathan
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology - Hyderabad 500032, India
| | - Kanak Garg
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology - Hyderabad 500032, India
| | - Yash Agarwal
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology - Hyderabad 500032, India
| | - Nita Parekh
- Centre for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology - Hyderabad 500032, India
| |
Collapse
|
70
|
Lake PT, Davidson RB, Klem H, Hocky GM, McCullagh M. Residue-Level Allostery Propagates through the Effective Coarse-Grained Hessian. J Chem Theory Comput 2020; 16:3385-3395. [PMID: 32251581 DOI: 10.1021/acs.jctc.9b01149] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The long-ranged coupling between residues that gives rise to allostery in a protein is built up from short-ranged physical interactions. Computational tools used to predict this coupling and its functional relevance have relied on the application of graph theoretical metrics to residue-level correlations measured from all-atom molecular dynamics simulations. The short-ranged interactions that yield these long-ranged residue-level correlations are quantified by the effective coarse-grained Hessian. Here we compute an effective harmonic coarse-grained Hessian from simulations of a benchmark allosteric protein, IGPS, and demonstrate the improved locality of this graph Laplacian over two other connectivity matrices. Additionally, two centrality metrics are developed that indicate the direct and indirect importance of each residue at producing the covariance between the effector binding pocket and the active site. The residue importance indicated by these two metrics is corroborated by previous mutagenesis experiments and leads to unique functional insights; in contrast to previous computational analyses, our results suggest that fP76-hK181 is the most important contact for conveying direct allosteric paths across the HisF-HisH interface. The connectivity around fD98 is found to be important at affecting allostery through indirect means.
Collapse
Affiliation(s)
- Peter T Lake
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Russell B Davidson
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Heidi Klem
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| | - Glen M Hocky
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Martin McCullagh
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523, United States
| |
Collapse
|
71
|
Dawson W, Mohr S, Ratcliff LE, Nakajima T, Genovese L. Complexity Reduction in Density Functional Theory Calculations of Large Systems: System Partitioning and Fragment Embedding. J Chem Theory Comput 2020; 16:2952-2964. [PMID: 32216343 DOI: 10.1021/acs.jctc.9b01152] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
With the development of low order scaling methods for performing Kohn-Sham density functional theory, it is now possible to perform fully quantum mechanical calculations of systems containing tens of thousands of atoms. However, with an increase in the size of the system treated comes an increase in complexity, making it challenging to analyze such large systems and determine the cause of emergent properties. To address this issue, in this paper, we present a systematic complexity reduction methodology which can break down large systems into their constituent fragments and quantify interfragment interactions. The methodology proposed here requires no a priori information or user interaction, allowing a single workflow to be automatically applied to any system of interest. We apply this approach to a variety of different systems and show how it allows for the derivation of new system descriptors, the design of QM/MM partitioning schemes, and the novel application of graph metrics to molecules and materials.
Collapse
Affiliation(s)
- William Dawson
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
| | - Stephan Mohr
- Barcelona Supercomputing Center (BSC), 08034 Barcelona, Spain
| | - Laura E Ratcliff
- Department of Materials, Imperial College London, London SW7 2AZ, United Kingdom
| | | | - Luigi Genovese
- Université Grenoble Alpes, INAC-MEM, L_Sim, Grenoble F-38000, France.,CEA, INAC-MEM, L_Sim, Grenoble F-38000, France
| |
Collapse
|
72
|
Miotto M, Di Rienzo L, Corsi P, Ruocco G, Raimondo D, Milanetti E. Simulated Epidemics in 3D Protein Structures to Detect Functional Properties. J Chem Inf Model 2020; 60:1884-1891. [PMID: 32011881 DOI: 10.1021/acs.jcim.9b01027] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The outcome of an epidemic is closely related to the network of interactions between individuals. Likewise, protein functions depend on the 3D arrangement of their residues and the underlying energetic interaction network. Borrowing ideas from the theoretical framework that has been developed to address the spreading of real diseases, we study for the first time the diffusion of a fictitious epidemic inside the protein nonbonded interaction network, aiming to study network features and properties. Our approach allows us to probe the overall stability and the capability of propagating information in complex 3D structures, proving to be very efficient in addressing different problems, from the assessment of thermal stability to the identification of functional sites.
Collapse
Affiliation(s)
- Mattia Miotto
- Department of Physics, Sapienza University, Rome 00185, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Rome 00161, Italy
| | | | - Pietro Corsi
- Department of Science, Roma Tre University, Rome 00154, Italy
| | - Giancarlo Ruocco
- Department of Physics, Sapienza University, Rome 00185, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Rome 00161, Italy
| | - Domenico Raimondo
- Department of Molecular Medicine, Sapienza University, Rome 00161, Italy
| | - Edoardo Milanetti
- Department of Physics, Sapienza University, Rome 00185, Italy.,Center for Life Nanoscience, Istituto Italiano di Tecnologia, Rome 00161, Italy
| |
Collapse
|
73
|
Lakhani B, Thayer KM, Black E, Beveridge DL. Spectral analysis of molecular dynamics simulations on PDZ: MD sectors. J Biomol Struct Dyn 2020; 38:781-790. [PMID: 31262238 PMCID: PMC7307555 DOI: 10.1080/07391102.2019.1588169] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 02/23/2019] [Indexed: 02/06/2023]
Abstract
The idea of protein "sectors" posits that sparse subsets of amino acid residues form cooperative networks that are key elements of protein stability, ligand binding, and allosterism. To date, protein sectors have been calculated by the statistical coupling analysis (SCA) method of Ranganathan and co-workers via the spectral analysis of conservation-weighted evolutionary covariance matrices obtained from a multiple sequence alignments of homologous families of proteins. SCA sectors, a knowledge-based protocol, have been indentified with functional properties and allosterism for a number of systems. In this study, we investigate the utility of the sector idea for the analysis of physics-based molecular dynamics (MD) trajectories of proteins. Our test case for this procedure is PSD95- PDZ3, one of the smallest proteins for which allosterism has been observed. It has served previously as a model system for a number of prediction algorithms, and is well characterized by X-ray crystallography, NMR spectroscopy and site specific mutagenisis. All-atom MD simulations were performed for a total of 500 nanoseconds using AMBER, and MD-calculated covariance matrices for the fluctuations of residue displacements and non-bonded interaction energies were subjected to spectral analysis in a manner analogous to that of SCA. The composition of MD sectors was compared with results from SCA, site specific mutagenesis, and allosterism. The concordance indicates that MD sectors are a viable protocol for analyzing MD trajectories and provide insight into the physical origin of the phenomenon.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Bharat Lakhani
- Program in Molecular Biophysics, Wesleyan University, Middletown CT 06459, USA
- Department of Molecular Biology & Biochemistry, Wesleyan University, Middletown CT 06459, USA
| | - Kelly M. Thayer
- Program in Molecular Biophysics, Wesleyan University, Middletown CT 06459, USA
- Chemistry Department, Wesleyan University, Middletown CT 06459, USA
- Department of Mathematics and Computer Science, Wesleyan University, Middletown CT 06459, USA
| | - Emily Black
- Program in Molecular Biophysics, Wesleyan University, Middletown CT 06459, USA
| | - David L. Beveridge
- Program in Molecular Biophysics, Wesleyan University, Middletown CT 06459, USA
- Chemistry Department, Wesleyan University, Middletown CT 06459, USA
| |
Collapse
|
74
|
Sheik Amamuddy O, Veldman W, Manyumwa C, Khairallah A, Agajanian S, Oluyemi O, Verkhivker GM, Tastan Bishop Ö. Integrated Computational Approaches and Tools forAllosteric Drug Discovery. Int J Mol Sci 2020; 21:E847. [PMID: 32013012 PMCID: PMC7036869 DOI: 10.3390/ijms21030847] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2019] [Revised: 01/20/2020] [Accepted: 01/21/2020] [Indexed: 12/16/2022] Open
Abstract
Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.
Collapse
Affiliation(s)
- Olivier Sheik Amamuddy
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Wayde Veldman
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Colleen Manyumwa
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Afrah Khairallah
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Steve Agajanian
- 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, CA 92866, USA; (S.A.); (O.O.)
| | - Odeyemi Oluyemi
- 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, CA 92866, USA; (S.A.); (O.O.)
| | - 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, CA 92866, USA; (S.A.); (O.O.)
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| |
Collapse
|
75
|
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
|
76
|
Heifetz A, Sladek V, Townsend-Nicholson A, Fedorov DG. Characterizing Protein-Protein Interactions with the Fragment Molecular Orbital Method. Methods Mol Biol 2020; 2114:187-205. [PMID: 32016895 DOI: 10.1007/978-1-0716-0282-9_13] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Proteins are vital components of living systems, serving as building blocks, molecular machines, enzymes, receptors, ion channels, sensors, and transporters. Protein-protein interactions (PPIs) are a key part of their function. There are more than 645,000 reported disease-relevant PPIs in the human interactome, but drugs have been developed for only 2% of these targets. The advances in PPI-focused drug discovery are highly dependent on the availability of structural data and accurate computational tools for analysis of this data. Quantum mechanical approaches are often too expensive computationally, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. FMO provides essential information for PPI drug discovery, namely, identification of key interactions formed between residues of two proteins, including their strength (in kcal/mol) and their chemical nature (electrostatic or hydrophobic). In this chapter, we have demonstrated how three different FMO-based approaches (pair interaction energy analysis (PIE analysis), subsystem analysis (SA) and analysis of protein residue networks (PRNs)) have been applied to study PPI in three protein-protein complexes.
Collapse
Affiliation(s)
| | - Vladimir Sladek
- Institute of Chemistry, Centre for Glycomics, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Andrea Townsend-Nicholson
- Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London, UK
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan.
| |
Collapse
|
77
|
Summers TJ, Daniel BP, Cheng Q, DeYonker NJ. Quantifying Inter-Residue Contacts through Interaction Energies. J Chem Inf Model 2019; 59:5034-5044. [DOI: 10.1021/acs.jcim.9b00804] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Thomas J. Summers
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Baty P. Daniel
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Qianyi Cheng
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Nathan J. DeYonker
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| |
Collapse
|
78
|
Establishing Computational Approaches Towards Identifying Malarial Allosteric Modulators: A Case Study of Plasmodium falciparum Hsp70s. Int J Mol Sci 2019; 20:ijms20225574. [PMID: 31717270 PMCID: PMC6887781 DOI: 10.3390/ijms20225574] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2019] [Revised: 10/24/2019] [Accepted: 10/27/2019] [Indexed: 02/07/2023] Open
Abstract
Combating malaria is almost a never-ending battle, as Plasmodium parasites develop resistance to the drugs used against them, as observed recently in artemisinin-based combination therapies. The main concern now is if the resistant parasite strains spread from Southeast Asia to Africa, the continent hosting most malaria cases. To prevent catastrophic results, we need to find non-conventional approaches. Allosteric drug targeting sites and modulators might be a new hope for malarial treatments. Heat shock proteins (HSPs) are potential malarial drug targets and have complex allosteric control mechanisms. Yet, studies on designing allosteric modulators against them are limited. Here, we identified allosteric modulators (SANC190 and SANC651) against P. falciparum Hsp70-1 and Hsp70-x, affecting the conformational dynamics of the proteins, delicately balanced by the endogenous ligands. Previously, we established a pipeline to identify allosteric sites and modulators. This study also further investigated alternative approaches to speed up the process by comparing all atom molecular dynamics simulations and dynamic residue network analysis with the coarse-grained (CG) versions of the calculations. Betweenness centrality (BC) profiles for PfHsp70-1 and PfHsp70-x derived from CG simulations not only revealed similar trends but also pointed to the same functional regions and specific residues corresponding to BC profile peaks.
Collapse
|
79
|
Sinha N, Chowdhury S, Sarkar RR. Molecular basis of drug resistance in smoothened receptor: An
in silico
study of protein resistivity and specificity. Proteins 2019; 88:514-526. [DOI: 10.1002/prot.25830] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 08/28/2019] [Accepted: 09/17/2019] [Indexed: 12/15/2022]
Affiliation(s)
- Noopur Sinha
- Chemical Engineering and Process Development DivisionCSIR‐National Chemical Laboratory Pune Maharashtra India
- Academy of Scientific and Innovative Research (AcSIR) Ghaziabad India
| | - Saikat Chowdhury
- Chemical Engineering and Process Development DivisionCSIR‐National Chemical Laboratory Pune Maharashtra India
- Academy of Scientific and Innovative Research (AcSIR) Ghaziabad India
| | - Ram Rup Sarkar
- Chemical Engineering and Process Development DivisionCSIR‐National Chemical Laboratory Pune Maharashtra India
- Academy of Scientific and Innovative Research (AcSIR) Ghaziabad India
| |
Collapse
|
80
|
Serçinoglu O, Ozbek P. gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations. Nucleic Acids Res 2019; 46:W554-W562. [PMID: 29800260 PMCID: PMC6030995 DOI: 10.1093/nar/gky381] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 05/22/2018] [Indexed: 11/12/2022] Open
Abstract
Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.
Collapse
Affiliation(s)
- Onur Serçinoglu
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| |
Collapse
|
81
|
Astl L, Verkhivker GM. Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks. Biochim Biophys Acta Gen Subj 2019:S0304-4165(19)30179-5. [PMID: 31330173 DOI: 10.1016/j.bbagen.2019.07.008] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2019] [Revised: 07/15/2019] [Accepted: 07/17/2019] [Indexed: 02/07/2023]
Abstract
BACKGROUND Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks. METHODS In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems. RESULTS The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators. CONCLUSIONS We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases. GENERAL SIGNIFICANCE This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.
Collapse
Affiliation(s)
- Lindy Astl
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, United States of America
| | - Gennady M Verkhivker
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, United States of America; Department of Pharmacology, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, United States of America.
| |
Collapse
|
82
|
Aydınkal RM, Serçinoğlu O, Ozbek P. ProSNEx: a web-based application for exploration and analysis of protein structures using network formalism. Nucleic Acids Res 2019; 47:W471-W476. [PMID: 31114881 PMCID: PMC6602423 DOI: 10.1093/nar/gkz390] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2019] [Revised: 04/17/2019] [Accepted: 05/09/2019] [Indexed: 01/14/2023] Open
Abstract
ProSNEx (Protein Structure Network Explorer) is a web service for construction and analysis of Protein Structure Networks (PSNs) alongside amino acid flexibility, sequence conservation and annotation features. ProSNEx constructs a PSN by adding nodes to represent residues and edges between these nodes using user-specified interaction distance cutoffs for either carbon-alpha, carbon-beta or atom-pair contact networks. Different types of weighted networks can also be constructed by using either (i) the residue-residue interaction energies in the format returned by gRINN, resulting in a Protein Energy Network (PEN); (ii) the dynamical cross correlations from a coarse-grained Normal Mode Analysis (NMA) of the protein structure; (iii) interaction strength. Upon construction of the network, common network metrics (such as node centralities) as well as shortest paths between nodes and k-cliques are calculated. Moreover, additional features of each residue in the form of conservation scores and mutation/natural variant information are included in the analysis. By this way, tool offers an enhanced and direct comparison of network-based residue metrics with other types of biological information. ProSNEx is free and open to all users without login requirement at http://prosnex-tool.com.
Collapse
Affiliation(s)
- Rasim Murat Aydınkal
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Ali Nihat Gokyigit Foundation, Etiler, Istanbul 34340, Turkey
| | - Onur Serçinoğlu
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
- Department of Bioengineering, Faculty of Engineering, Recep Tayyip Erdoğan University, Rize 53100, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Faculty of Engineering, Marmara University, Kadikoy, Istanbul 34722, Turkey
| |
Collapse
|
83
|
Yao XQ, Momin M, Hamelberg D. Establishing a Framework of Using Residue–Residue Interactions in Protein Difference Network Analysis. J Chem Inf Model 2019; 59:3222-3228. [DOI: 10.1021/acs.jcim.9b00320] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Mohamed Momin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
84
|
Xia Q, Ding Y. Thermostability of Lipase A and Dynamic Communication Based on Residue Interaction Network. Protein Pept Lett 2019; 26:702-716. [PMID: 31215367 DOI: 10.2174/0929866526666190617091812] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Revised: 04/10/2019] [Accepted: 04/25/2019] [Indexed: 11/22/2022]
Abstract
OBJECTIVE Dynamic communication caused by mutation affects protein stability. The main objective of this study is to explore how mutations affect communication and to provide further insight into the relationship between heat resistance and signal propagation of Bacillus subtilis lipase (Lip A). METHODS The relationship between dynamic communication and Lip A thermostability is studied by long-time MD simulation and residue interaction network. The Dijkstra algorithm is used to get the shortest path of each residue pair. Subsequently, time-series frequent paths and spatio-temporal frequent paths are mined through an Apriori-like algorithm. RESULTS Time-series frequent paths show that the communication between residue pairs, both in wild-type lipase (WTL) and mutant 6B, becomes chaotic with an increase in temperature; however, more residues in 6B can maintain stable communication at high temperature, which may be associated with the structural rigidity. Furthermore, spatio-temporal frequent paths reflect the interactions among secondary structures. For WTL at 300K, β7, αC, αB, the longest loop, αA and αF contact frequently. The 310-helix between β3 and αA is penetrated by spatio-temporal frequent paths. At 400K, only αC can be frequently transmitted. For 6B, when at 300K, αA and αF are in more tight contact by spatio-temporal frequent paths though I157M and N166Y. Moreover, the rigidity of the active site His156 and the C-terminal of Lip A are increased, as reflected by the spatio-temporal frequent paths. At 400K, αA and αF, 310-helix between β3 and αA, the longest loop, and the loop where the active site Asp133 is located can still maintain stable communication. CONCLUSION From the perspective of residue dynamic communication, it is obviously found that mutations cause changes in interactions between secondary structures and enhance the rigidity of the structure, contributing to the thermal stability and functional activity of 6B.
Collapse
Affiliation(s)
- Qian Xia
- Laboratory of Media Design and Software Technology, Jiangnan University, Wuxi, Jiangsu, 214122, China
| | - Yanrui Ding
- Laboratory of Media Design and Software Technology, Jiangnan University, Wuxi, Jiangsu, 214122, China.,Key Laboratory of Industrial Biotechnology, Jiangnan University, Wuxi, Jiangsu, 214122, China
| |
Collapse
|
85
|
Agajanian S, Oluyemi O, Verkhivker GM. Integration of Random Forest Classifiers and Deep Convolutional Neural Networks for Classification and Biomolecular Modeling of Cancer Driver Mutations. Front Mol Biosci 2019; 6:44. [PMID: 31245384 PMCID: PMC6579812 DOI: 10.3389/fmolb.2019.00044] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 05/23/2019] [Indexed: 12/21/2022] Open
Abstract
Development of machine learning solutions for prediction of functional and clinical significance of cancer driver genes and mutations are paramount in modern biomedical research and have gained a significant momentum in a recent decade. In this work, we integrate different machine learning approaches, including tree based methods, random forest and gradient boosted tree (GBT) classifiers along with deep convolutional neural networks (CNN) for prediction of cancer driver mutations in the genomic datasets. The feasibility of CNN in using raw nucleotide sequences for classification of cancer driver mutations was initially explored by employing label encoding, one hot encoding, and embedding to preprocess the DNA information. These classifiers were benchmarked against their tree-based alternatives in order to evaluate the performance on a relative scale. We then integrated DNA-based scores generated by CNN with various categories of conservational, evolutionary and functional features into a generalized random forest classifier. The results of this study have demonstrated that CNN can learn high level features from genomic information that are complementary to the ensemble-based predictors often employed for classification of cancer mutations. By combining deep learning-generated score with only two main ensemble-based functional features, we can achieve a superior performance of various machine learning classifiers. Our findings have also suggested that synergy of nucleotide-based deep learning scores and integrated metrics derived from protein sequence conservation scores can allow for robust classification of cancer driver mutations with a limited number of highly informative features. Machine learning predictions are leveraged in molecular simulations, protein stability, and network-based analysis of cancer mutations in the protein kinase genes to obtain insights about molecular signatures of driver mutations and enhance the interpretability of cancer-specific classification models.
Collapse
Affiliation(s)
- Steve Agajanian
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, United States
| | - Odeyemi Oluyemi
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, United States
| | - Gennady M Verkhivker
- Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, United States
| |
Collapse
|
86
|
Astl L, Verkhivker GM. Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications. J Chem Theory Comput 2019; 15:3362-3380. [PMID: 31017783 DOI: 10.1021/acs.jctc.9b00119] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.
Collapse
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.,Department of Biomedical and Pharmaceutical Sciences , Chapman University School of Pharmacy , Irvine , California 92618 , United States
| |
Collapse
|
87
|
Smith IN, Thacker S, Seyfi M, Cheng F, Eng C. Conformational Dynamics and Allosteric Regulation Landscapes of Germline PTEN Mutations Associated with Autism Compared to Those Associated with Cancer. Am J Hum Genet 2019; 104:861-878. [PMID: 31006514 PMCID: PMC6506791 DOI: 10.1016/j.ajhg.2019.03.009] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2018] [Accepted: 03/08/2019] [Indexed: 01/07/2023] Open
Abstract
Individuals with germline PTEN tumor-suppressor variants have PTEN hamartoma tumor syndrome (PHTS). Clinically, PHTS has variable presentations; there are distinct subsets of PHTS-affected individuals, such as those diagnosed with autism spectrum disorder (ASD) or cancer. It remains unclear why mutations in one gene can lead to such seemingly disparate phenotypes. Therefore, we sought to determine whether it is possible to predict a given PHTS-affected individual's a priori risk of ASD, cancer, or the co-occurrence of both phenotypes. By integrating network proximity analysis performed on the human interactome, molecular simulations, and residue-interaction networks, we demonstrate the role of conformational dynamics in the structural communication and long-range allosteric regulation of germline PTEN variants associated with ASD or cancer. We show that the PTEN interactome shares significant overlap with the ASD and cancer interactomes, providing network-based evidence that PTEN is a crucial player in the biology of both disorders. Importantly, this finding suggests that a germline PTEN variant might perturb the ASD or cancer networks differently, thus favoring one disease outcome at any one time. Furthermore, protein-dynamic structural-network analysis reveals small-world structural communication mediated by highly conserved functional residues and potential allosteric regulation of PTEN. We identified a salient structural-communication pathway that extends across the inter-domain interface for cancer-only mutations. In contrast, the structural-communication pathway is predominantly restricted to the phosphatase domain for ASD-only mutations. Our integrative approach supports the prediction and potential modulation of the relevant conformational states that influence structural communication and long-range perturbations associated with mutational effects that lead to PTEN-ASD or PTEN-cancer phenotypes.
Collapse
Affiliation(s)
- Iris Nira Smith
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Stetson Thacker
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA
| | - Marilyn Seyfi
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA
| | - Feixiong Cheng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA; Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Charis Eng
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Cleveland Clinic Lerner College of Medicine, Case Western Reserve University, Cleveland, OH 44195, USA; Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA; Taussig Cancer Institute, Cleveland Clinic, Cleveland, OH 44195, USA; Department of Genetics and Genome Sciences, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| |
Collapse
|
88
|
Johnson LE, Ginovska B, Fenton AW, Raugei S. Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example. Biophys J 2019; 116:1598-1608. [PMID: 31010662 DOI: 10.1016/j.bpj.2019.03.026] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Revised: 03/14/2019] [Accepted: 03/21/2019] [Indexed: 12/14/2022] Open
Abstract
Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.
Collapse
Affiliation(s)
- Lewis E Johnson
- Department of Chemistry, University of Washington, Seattle, Washington; Physical and Computational Sciences Directorate, Pacific Northwestern National Laboratory, Richland, Washington
| | - Bojana Ginovska
- Physical and Computational Sciences Directorate, Pacific Northwestern National Laboratory, Richland, Washington
| | - Aron W Fenton
- Department of Biochemistry and Molecular Biology, The University of Kansas Medical Center, Kansas City, Kansas
| | - Simone Raugei
- Physical and Computational Sciences Directorate, Pacific Northwestern National Laboratory, Richland, Washington.
| |
Collapse
|
89
|
Liang Z, Verkhivker GM, Hu G. Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications. Brief Bioinform 2019; 21:815-835. [DOI: 10.1093/bib/bbz029] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Revised: 02/04/2019] [Accepted: 02/21/2019] [Indexed: 12/24/2022] Open
Abstract
Abstract
Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.
Collapse
Affiliation(s)
- Zhongjie Liang
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA, USA
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA, USA
| | - Guang Hu
- School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| |
Collapse
|
90
|
Dissecting a novel allosteric mechanism of cruzain: A computer-aided approach. PLoS One 2019; 14:e0211227. [PMID: 30682119 PMCID: PMC6347273 DOI: 10.1371/journal.pone.0211227] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Accepted: 01/09/2019] [Indexed: 02/08/2023] Open
Abstract
Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases.
Collapse
|
91
|
Astl L, Tse A, Verkhivker GM. Interrogating Regulatory Mechanisms in Signaling Proteins by Allosteric Inhibitors and Activators: A Dynamic View Through the Lens of Residue Interaction Networks. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2019; 1163:187-223. [DOI: 10.1007/978-981-13-8719-7_9] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
92
|
Sladek V, Tokiwa H, Shimano H, Shigeta Y. Protein Residue Networks from Energetic and Geometric Data: Are They Identical? J Chem Theory Comput 2018; 14:6623-6631. [PMID: 30500196 DOI: 10.1021/acs.jctc.8b00733] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Protein residue networks (PRN) from energetic and geometric data are probably not identical. PRNs constructed from ab initio pair interaction energies are analyzed for the first time and compared to PRN based on center of mass separation. We use modern, previously unused algorithms such as global and local efficiencies to quantitatively confirm that both types of PRNs do exhibit small-world character. The main novelty finding is that interaction energy-based PRNs preserve small-world character even when clustered. A node hierarchy independent of the cutoff energy used for the edge creation is characteristic for them. Efficiency centrality identifies hubs responsible for such behavior. The interaction energy-based PRNs seem to comply with the scale-free network model with respect to efficiency centrality distribution as opposed to distance based PRNs. Community detection is introduced into protein network research as an extension beyond cluster analysis to study tertiary and quaternary structures.
Collapse
Affiliation(s)
- Vladimir Sladek
- Institute of Chemistry - Centre for Glycomics , Dubravska cesta 9 , 84538 Bratislava , Slovakia.,Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan
| | - Hiroaki Tokiwa
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Chemistry , Rikkyo University , Nishi-Ikebukuro , Toshima, Tokyo 171-8501 , Japan
| | - Hitoshi Shimano
- Agency for Medical Research and Development (AMED) , Chiyoda-ku , Japan.,Department of Internal Medicine, Faculty of Medicine , University of Tsukuba , 1-1-1 Tennodai , Tsukuba, Ibaraki 305-8575 , Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences , University of Tsukuba , Tennodai 1-1-1 , Tsukuba, Ibaraki 305-8577 , Japan
| |
Collapse
|
93
|
Sahu TK, Pradhan D, Rao AR, Jena L. In silico site-directed mutagenesis of neutralizing mAb 4C4 and analysis of its interaction with G-H loop of VP1 to explore its therapeutic applications against FMD. J Biomol Struct Dyn 2018; 37:2641-2651. [PMID: 30051760 DOI: 10.1080/07391102.2018.1494631] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Investigating the behaviour of bio-molecules through computational mutagenesis is gaining interest to facilitate the development of new therapeutic solutions for infectious diseases. The antigenetically variant genotypes of foot and mouth disease virus (FMDV) and their subsequent infections are challenging to tackle with traditional vaccination. In such scenario, neutralizing antibodies might provide an alternate solution to manage the FMDV infection. Thus, we have analysed the interaction of the mAb 4C4 with a synthetic G-H loop of FMDV-VP1 through in silico mutagenesis and molecular modelling. Initially, a set of 25,434 mutants were designed and the mutants having better energetic stability than 4C4 were clustered based on sequence identity. The best mutant representing each cluster was selected and evaluated for its binding affinity with the antigen in terms of docking scores, interaction energy and binding energy. Six mutants have confirmed better binding affinities towards the antigen than 4C4. Further, interaction of these mutants with the natural G-H loop that is bound to mAb SD6 was also evaluated. One 4C4 variant having mutations at the positions 2034(N→L), 2096(N→C), 2098(D→Y), 2532(T→K) and 2599(A→G) has revealed better binding affinities towards both the synthetic and natural G-H loops than 4C4 and SD6, respectively. A molecular dynamic simulation for 50 ns was conducted for mutant and wild-type antibody structures which supported the pre-simulation results. Therefore, these mutations on mAb 4C4 are believed to provide a better antibody-based therapeutic option for FMD. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Tanmaya Kumar Sahu
- a Centre for Agricultural Bioinformatics , ICAR-Indian Agricultural Statistics Research Institute , New Delhi , Delhi , India
| | - Dibyabhaba Pradhan
- b Biomedical Informatics Centre , ICMR-National Institute of Pathology , New Delhi , Delhi , India.,c ICMR-Computational Genomics Centre , Indian Council of Medical Research , New Delhi , Delhi , India
| | - Atmakuri Ramakrishna Rao
- a Centre for Agricultural Bioinformatics , ICAR-Indian Agricultural Statistics Research Institute , New Delhi , Delhi , India
| | - Lingaraj Jena
- d Bioinformatics Centre , Mahatma Gandhi Institute of Medical Sciences , Sevagram , Maharashtra , India
| |
Collapse
|
94
|
Verkhivker GM. Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication. Biochim Biophys Acta Gen Subj 2018; 1863:210-225. [PMID: 30339916 DOI: 10.1016/j.bbagen.2018.10.009] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2018] [Revised: 10/06/2018] [Accepted: 10/13/2018] [Indexed: 12/19/2022]
Abstract
In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, United States; Department of Pharmacology, Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, 9500 Gilman Drive, La Jolla, CA 92093, USA.
| |
Collapse
|
95
|
Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function. Curr Opin Struct Biol 2018; 54:1-9. [PMID: 30268910 PMCID: PMC6420056 DOI: 10.1016/j.sbi.2018.09.004] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Revised: 08/13/2018] [Accepted: 09/10/2018] [Indexed: 01/12/2023]
Abstract
A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.
Collapse
|
96
|
Agajanian S, Odeyemi O, Bischoff N, Ratra S, Verkhivker GM. Machine Learning Classification and Structure–Functional Analysis of Cancer Mutations Reveal Unique Dynamic and Network Signatures of Driver Sites in Oncogenes and Tumor Suppressor Genes. J Chem Inf Model 2018; 58:2131-2150. [DOI: 10.1021/acs.jcim.8b00414] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Steve Agajanian
- Graduate Program in Computational and Data Sciences, Department of Computational Sciences, Schmid College of Science and Technology, Chapman University, One University
Drive, Orange, California 92866, United States
| | - Oluyemi Odeyemi
- Graduate Program in Computational and Data Sciences, Department of Computational Sciences, Schmid College of Science and Technology, Chapman University, One University
Drive, Orange, California 92866, United States
| | - Nathaniel Bischoff
- Graduate Program in Computational and Data Sciences, Department of Computational Sciences, Schmid College of Science and Technology, Chapman University, One University
Drive, Orange, California 92866, United States
| | - Simrath Ratra
- Graduate Program in Computational and Data Sciences, Department of Computational Sciences, 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, Department of Computational Sciences, Schmid College of Science and Technology, Chapman University, One University
Drive, Orange, California 92866, United States
- Chapman University, School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
97
|
Yao XQ, Momin M, Hamelberg D. Elucidating Allosteric Communications in Proteins with Difference Contact Network Analysis. J Chem Inf Model 2018; 58:1325-1330. [DOI: 10.1021/acs.jcim.8b00250] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xin-Qiu Yao
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Mohamed Momin
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| | - Donald Hamelberg
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-3965, United States
| |
Collapse
|
98
|
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
|
99
|
Verkhivker GM. Dynamics-based community analysis and perturbation response scanning of allosteric interaction networks in the TRAP1 chaperone structures dissect molecular linkage between conformational asymmetry and sequential ATP hydrolysis. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2018; 1866:899-912. [PMID: 29684503 DOI: 10.1016/j.bbapap.2018.04.008] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 02/10/2018] [Revised: 04/12/2018] [Accepted: 04/13/2018] [Indexed: 01/29/2023]
Abstract
Allosteric interactions of the Hsp90 chaperones with cochaperones and diverse protein clients can often exhibit distinct asymmetric features that determine regulatory mechanisms and cellular functions in many signaling networks. The recent crystal structures of the mitochondrial Hsp90 isoform TRAP1 in complexes with ATP analogs have provided first evidence of significant asymmetry in the closed dimerized state that triggers independent activity of the chaperone protomers, whereby preferential hydrolysis of the buckled protomer is followed by conformational flipping between protomers and hydrolysis of the second protomer. Despite significant insights in structural characterizations of the TRAP1 chaperone, the atomistic details and mechanics of allosteric interactions that couple sequential ATP hydrolysis with asymmetric conformational switching in the TRAP1 protomers remain largely unknown. In this work, we explored atomistic and coarse-grained simulations of the TRAP1 dimer structures in combination with the ensemble-based network modeling and perturbation response scanning of residue interaction networks to probe salient features underlying allosteric signaling mechanism. This study has revealed that key effector sites that orchestrate allosteric interactions occupy the ATP binding region and N-terminal interface of the buckled protomer, whereas the main sensors of allosteric signals that drive functional conformational changes during ATPase cycle are consolidated near the client binding region of the straight protomer, channeling the energy of ATP hydrolysis for client remodeling. The community decomposition analysis of the interaction networks and reconstruction of allosteric communication pathways in the TRAP1 structures have quantified mechanism of allosteric regulation, revealing control points and interactions that coordinate asymmetric switching during ATP hydrolysis.
Collapse
Affiliation(s)
- Gennady M Verkhivker
- Department of Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, CA 92866, United States; Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, United States.
| |
Collapse
|
100
|
Site-Mutation of Hydrophobic Core Residues Synchronically Poise Super Interleukin 2 for Signaling: Identifying Distant Structural Effects through Affordable Computations. Int J Mol Sci 2018; 19:ijms19030916. [PMID: 29558421 PMCID: PMC5877777 DOI: 10.3390/ijms19030916] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Revised: 03/15/2018] [Accepted: 03/17/2018] [Indexed: 12/17/2022] Open
Abstract
A superkine variant of interleukin-2 with six site mutations away from the binding interface developed from the yeast display technique has been previously characterized as undergoing a distal structure alteration which is responsible for its super-potency and provides an elegant case study with which to get insight about how to utilize allosteric effect to achieve desirable protein functions. By examining the dynamic network and the allosteric pathways related to those mutated residues using various computational approaches, we found that nanosecond time scale all-atom molecular dynamics simulations can identify the dynamic network as efficient as an ensemble algorithm. The differentiated pathways for the six core residues form a dynamic network that outlines the area of structure alteration. The results offer potentials of using affordable computing power to predict allosteric structure of mutants in knowledge-based mutagenesis.
Collapse
|