1
|
Poorinmohammad N, Salavati R. Prioritization of Trypanosoma brucei editosome protein interactions interfaces at residue resolution through proteome-scale network analysis. BMC Mol Cell Biol 2024; 25:3. [PMID: 38279116 PMCID: PMC10811811 DOI: 10.1186/s12860-024-00499-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 01/19/2024] [Indexed: 01/28/2024] Open
Abstract
BACKGROUND Trypanosoma brucei is the causative agent for trypanosomiasis in humans and livestock, which presents a growing challenge due to drug resistance. While identifying novel drug targets is vital, the process is delayed due to a lack of functional information on many of the pathogen's proteins. Accordingly, this paper presents a computational framework for prioritizing drug targets within the editosome, a vital molecular machinery responsible for mitochondrial RNA processing in T. brucei. Importantly, this framework may eliminate the need for prior gene or protein characterization, potentially accelerating drug discovery efforts. RESULTS By integrating protein-protein interaction (PPI) network analysis, PPI structural modeling, and residue interaction network (RIN) analysis, we quantitatively ranked and identified top hub editosome proteins, their key interaction interfaces, and hotspot residues. Our findings were cross-validated and further prioritized by incorporating them into gene set analysis and differential expression analysis of existing quantitative proteomics data across various life stages of T. brucei. In doing so, we highlighted PPIs such as KREL2-KREPA1, RESC2-RESC1, RESC12A-RESC13, and RESC10-RESC6 as top candidates for further investigation. This includes examining their interfaces and hotspot residues, which could guide drug candidate selection and functional studies. CONCLUSION RNA editing offers promise for target-based drug discovery, particularly with proteins and interfaces that play central roles in the pathogen's life cycle. This study introduces an integrative drug target identification workflow combining information from the PPI network, PPI 3D structure, and reside-level information of their interface which can be applicable to diverse pathogens. In the case of T. brucei, via this pipeline, the present study suggested potential drug targets with residue-resolution from RNA editing machinery. However, experimental validation is needed to fully realize its potential in advancing urgently needed antiparasitic drug development.
Collapse
Affiliation(s)
- Naghmeh Poorinmohammad
- Institute of Parasitology, McGill University, Ste. Anne de Bellevue, Montreal, Quebec, H9X 3V9, Canada
| | - Reza Salavati
- Institute of Parasitology, McGill University, Ste. Anne de Bellevue, Montreal, Quebec, H9X 3V9, Canada.
- Department of Biochemistry, McGill University, Montreal, Quebec, H3G 1Y6, Canada.
| |
Collapse
|
2
|
Pandey M, Shah SK, Gromiha MM. Computational approaches for identifying disease-causing mutations in proteins. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2023; 139:141-171. [PMID: 38448134 DOI: 10.1016/bs.apcsb.2023.11.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
Advancements in genome sequencing have expanded the scope of investigating mutations in proteins across different diseases. Amino acid mutations in a protein alter its structure, stability and function and some of them lead to diseases. Identification of disease-causing mutations is a challenging task and it will be helpful for designing therapeutic strategies. Hence, mutation data available in the literature have been curated and stored in several databases, which have been effectively utilized for developing computational methods to identify deleterious mutations (drivers), using sequence and structure-based properties of proteins. In this chapter, we describe the contents of specific databases that have information on disease-causing and neutral mutations followed by sequence and structure-based properties. Further, characteristic features of disease-causing mutations will be discussed along with computational methods for identifying cancer hotspot residues and disease-causing mutations in proteins.
Collapse
Affiliation(s)
- Medha Pandey
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India
| | - Suraj Kumar Shah
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India
| | - M Michael Gromiha
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, India; International Research Frontiers Initiative, School of Computing, Tokyo Institute of Technology, Yokohama, Japan.
| |
Collapse
|
3
|
Nussinov R, Zhang M, Maloney R, Liu Y, Tsai CJ, Jang H. Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs. J Mol Biol 2022; 434:167569. [PMID: 35378118 PMCID: PMC9398924 DOI: 10.1016/j.jmb.2022.167569] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2022] [Revised: 03/11/2022] [Accepted: 03/25/2022] [Indexed: 01/12/2023]
Abstract
Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.
Collapse
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA; Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel.
| | - Mingzhen Zhang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| | - Ryan Maloney
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| | - Yonglan Liu
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| | - Chung-Jung Tsai
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| | - Hyunbum Jang
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Laboratory of Cancer Immunometabolism, National Cancer Institute, Frederick, MD 21702, USA
| |
Collapse
|
4
|
Xiao F, Zhou Z, Song X, Gan M, Long J, Verkhivker G, Hu G. Dissecting mutational allosteric effects in alkaline phosphatases associated with different Hypophosphatasia phenotypes: An integrative computational investigation. PLoS Comput Biol 2022; 18:e1010009. [PMID: 35320273 PMCID: PMC8979438 DOI: 10.1371/journal.pcbi.1010009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Revised: 04/04/2022] [Accepted: 03/10/2022] [Indexed: 11/18/2022] Open
Abstract
Hypophosphatasia (HPP) is a rare inherited disorder characterized by defective bone mineralization and is highly variable in its clinical phenotype. The disease occurs due to various loss-of-function mutations in ALPL, the gene encoding tissue-nonspecific alkaline phosphatase (TNSALP). In this work, a data-driven and biophysics-based approach is proposed for the large-scale analysis of ALPL mutations-from nonpathogenic to severe HPPs. By using a pipeline of synergistic approaches including sequence-structure analysis, network modeling, elastic network models and atomistic simulations, we characterized allosteric signatures and effects of the ALPL mutations on protein dynamics and function. Statistical analysis of molecular features computed for the ALPL mutations showed a significant difference between the control, mild and severe HPP phenotypes. Molecular dynamics simulations coupled with protein structure network analysis were employed to analyze the effect of single-residue variation on conformational dynamics of TNSALP dimers, and the developed machine learning model suggested that the topological network parameters could serve as a robust indicator of severe mutations. The results indicated that the severity of disease-associated mutations is often linked with mutation-induced modulation of allosteric communications in the protein. This study suggested that ALPL mutations associated with mild and more severe HPPs can exert markedly distinct effects on the protein stability and long-range network communications. By linking the disease phenotypes with dynamic and allosteric molecular signatures, the proposed integrative computational approach enabled to characterize and quantify the allosteric effects of ALPL mutations and role of allostery in the pathogenesis of HPPs.
Collapse
Affiliation(s)
- Fei Xiao
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Ziyun Zhou
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Xingyu Song
- Department of Chemistry, Multiscale Research Institute of Complex Systems and Institute of Biomedical Sciences, Fudan University, Shanghai, China
| | - Mi Gan
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Jie Long
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
| | - Gennady Verkhivker
- Department of Computational and Data Sciences, Chapman University, One University Drive, Orange, California, United States of America
- Department of Biomedical and Pharmaceutical Sciences, Chapman University Pharmacy School 9401 Jeronimo Rd, Irvine, California, United States of America
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou, China
- * E-mail:
| |
Collapse
|
5
|
Zhu Y, Ye F, Zhou Z, Liu W, Liang Z, Hu G. Insights into Conformational Dynamics and Allostery in DNMT1-H3Ub/USP7 Interactions. Molecules 2021; 26:molecules26175153. [PMID: 34500587 PMCID: PMC8434485 DOI: 10.3390/molecules26175153] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Revised: 08/23/2021] [Accepted: 08/23/2021] [Indexed: 01/03/2023] Open
Abstract
DNA methyltransferases (DNMTs) including DNMT1 are a conserved family of cytosine methylases that play crucial roles in epigenetic regulation. The versatile functions of DNMT1 rely on allosteric networks between its different interacting partners, emerging as novel therapeutic targets. In this work, based on the modeling structures of DNMT1-ubiquitylated H3 (H3Ub)/ubiquitin specific peptidase 7 (USP7) complexes, we have used a combination of elastic network models, molecular dynamics simulations, structural residue perturbation, network modeling, and pocket pathway analysis to examine their molecular mechanisms of allosteric regulation. The comparative intrinsic and conformational dynamics analysis of three DNMT1 systems has highlighted the pivotal role of the RFTS domain as the dynamics hub in both intra- and inter-molecular interactions. The site perturbation and network modeling approaches have revealed the different and more complex allosteric interaction landscape in both DNMT1 complexes, involving the events caused by mutational hotspots and post-translation modification sites through protein-protein interactions (PPIs). Furthermore, communication pathway analysis and pocket detection have provided new mechanistic insights into molecular mechanisms underlying quaternary structures of DNMT1 complexes, suggesting potential targeting pockets for PPI-based allosteric drug design.
Collapse
Affiliation(s)
- Yu Zhu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China; (Y.Z.); (Z.Z.); (W.L.)
| | - Fei Ye
- College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China;
| | - Ziyun Zhou
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China; (Y.Z.); (Z.Z.); (W.L.)
| | - Wanlin Liu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China; (Y.Z.); (Z.Z.); (W.L.)
| | - Zhongjie Liang
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China; (Y.Z.); (Z.Z.); (W.L.)
- Correspondence: (Z.L.); (G.H.)
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China; (Y.Z.); (Z.Z.); (W.L.)
- Correspondence: (Z.L.); (G.H.)
| |
Collapse
|
6
|
Echave J. Fast computational mutation-response scanning of proteins. PeerJ 2021; 9:e11330. [PMID: 33976988 PMCID: PMC8067912 DOI: 10.7717/peerj.11330] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 03/31/2021] [Indexed: 12/21/2022] Open
Abstract
Studying the effect of perturbations on protein structure is a basic approach in protein research. Important problems, such as predicting pathological mutations and understanding patterns of structural evolution, have been addressed by computational simulations that model mutations using forces and predict the resulting deformations. In single mutation-response scanning simulations, a sensitivity matrix is obtained by averaging deformations over point mutations. In double mutation-response scanning simulations, a compensation matrix is obtained by minimizing deformations over pairs of mutations. These very useful simulation-based methods may be too slow to deal with large proteins, protein complexes, or large protein databases. To address this issue, I derived analytical closed formulas to calculate the sensitivity and compensation matrices directly, without simulations. Here, I present these derivations and show that the resulting analytical methods are much faster than their simulation counterparts.
Collapse
Affiliation(s)
- Julian Echave
- Instituto de Ciencias Físicas, Escuela de Ciencia y Tecnología, Universidad Nacional de San Martín, San Martín, Buenos Aires, Argentina
| |
Collapse
|
7
|
Sora V, Sanchez D, Papaleo E. Bcl-xL Dynamics under the Lens of Protein Structure Networks. J Phys Chem B 2021; 125:4308-4320. [PMID: 33848145 DOI: 10.1021/acs.jpcb.0c11562] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Understanding the finely orchestrated interactions leading to or preventing programmed cell death (apoptosis) is of utmost importance in cancer research because the failure of these systems could eventually lead to the onset of the disease. In this regard, the maintenance of a delicate balance between the promoters and inhibitors of mitochondrial apoptosis is crucial, as demonstrated by the interplay among the Bcl-2 family members. In particular, B-cell lymphoma extra-large (Bcl-xL) is a target of interest due to the forefront role of its dysfunctions in cancer development. Bcl-xL prevents apoptosis by binding both the pro-apoptotic BH3-only proteins, like PUMA, and the noncanonical partners, such as p53, at different sites. An allosteric communication between the BH3-only protein binding pocket and the p53 binding site, mediating the release of p53 from Bcl-xL upon PUMA binding, has been postulated and supported by nuclear magnetic resonance and other biophysical data. The molecular details of this mechanism, especially at the residue level, remain unclear. In this work, we investigated the distal communication between these two sites in Bcl-xL in its free state and when bound to PUMA. We also evaluated how missense mutations of Bcl-xL found in cancer samples might impair this communication and therefore the allosteric mechanism. We employed all-atom explicit solvent microsecond molecular dynamics simulations, analyzed through a Protein Structure Network approach and integrated with calculations of changes in free energies upon cancer-related mutations identified by genomics studies. We found a subset of candidate residues responsible for both maintaining protein stability and for conveying structural information between the two binding sites and hypothesized possible communication routes between specific residues at both sites.
Collapse
Affiliation(s)
- Valentina Sora
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark.,Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Dionisio Sanchez
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, 2100 Copenhagen, Denmark.,Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| |
Collapse
|
8
|
Abstract
Network theory provides one of the most potent analysis tools for the study of complex systems. In this paper, we illustrate the network-based perspective in drug research and how it is coherent with the new paradigm of drug discovery. We first present data sources from which networks are built, then show some examples of how the networks can be used to investigate drug-related systems. A section is devoted to network-based inference applications, i.e., prediction methods based on interactomes, that can be used to identify putative drug-target interactions without resorting to 3D modeling. Finally, we present some aspects of Boolean networks dynamics, anticipating that it might become a very potent modeling framework to develop in silico screening protocols able to simulate phenotypic screening experiments. We conclude that network applications integrated with machine learning and 3D modeling methods will become an indispensable tool for computational drug discovery in the next years.
Collapse
Affiliation(s)
- Maurizio Recanatini
- Department of Pharmacy and
Biotechnology, Alma Mater Studiorum—University of Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Chiara Cabrelle
- Department of Pharmacy and
Biotechnology, Alma Mater Studiorum—University of Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| |
Collapse
|
9
|
Xiao F, Song X, Tian P, Gan M, Verkhivker GM, Hu G. Comparative Dynamics and Functional Mechanisms of the CYP17A1 Tunnels Regulated by Ligand Binding. J Chem Inf Model 2020; 60:3632-3647. [PMID: 32530640 DOI: 10.1021/acs.jcim.0c00447] [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/16/2022]
Abstract
As an important member of cytochrome P450 (CYP) enzymes, CYP17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones, making it an attractive therapeutic target. The emerging structural information about CYP17A1 and the growing number of inhibitors for these enzymes call for a systematic strategy to delineate and classify mechanisms of ligand transport through tunnels that control catalytic activity. In this work, we applied an integrated computational strategy to different CYP17A1 systems with a panel of ligands to systematically study at the atomic level the mechanism of ligand-binding and tunneling dynamics. Atomistic simulations and binding free energy computations identify the dynamics of dominant tunnels and characterize energetic properties of critical residues responsible for ligand binding. The common transporting pathways including S, 3, and 2c tunnels were identified in CYP17A1 binding systems, while the 2c tunnel is a newly formed pathway upon ligand binding. We employed and integrated several computational approaches including the analysis of functional motions and sequence conservation, atomistic modeling of dynamic residue interaction networks, and perturbation response scanning analysis to dissect ligand tunneling mechanisms. The results revealed the hinge-binding and sliding motions as main functional modes of the tunnel dynamic, and a group of mediating residues as key regulators of tunnel conformational dynamics and allosteric communications. We have also examined and quantified the mutational effects on the tunnel composition, conformational dynamics, and long-range allosteric behavior. The results of this investigation are fully consistent with the experimental data, providing novel rationale to the experiments and offering valuable insights into the relationships between the structure and function of the channel networks and a robust atomistic model of activation mechanisms and allosteric interactions in CYP enzymes.
Collapse
Affiliation(s)
- Fei Xiao
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Xingyu Song
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Peiyi Tian
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Mi Gan
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Chapman University, One University Drive, Orange, California 92866, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University Pharmacy School, 9401 Jeronimo Rd, Irvine, California 92618, United States
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| |
Collapse
|
10
|
Jalalypour F, Sensoy O, Atilgan C. Perturb-Scan-Pull: A Novel Method Facilitating Conformational Transitions in Proteins. J Chem Theory Comput 2020; 16:3825-3841. [PMID: 32324386 DOI: 10.1021/acs.jctc.9b01222] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Conformational transitions in proteins facilitate precise physiological functions. Therefore, it is crucial to understand the mechanisms underlying these processes to modulate protein function. Yet, studying structural and dynamical properties of proteins is notoriously challenging due to the complexity of the underlying potential energy surfaces (PES). We have previously developed the perturbation-response scanning (PRS) method to identify key residues that participate in the communication network responsible for specific conformational transitions. PRS is based on a residue-by-residue scan of the protein to determine the subset of residues/forces which provide the closest conformational change leading to a target conformational state, inasmuch as linear response theory applies to these motions. Here, we develop a novel method to further evaluate if conformational transitions may be triggered on the PES. We aim to study functionally relevant conformational transitions in proteins by using results obtained from PRS and feeding them as inputs to steered molecular dynamics simulations. The success and the transferability of the method are evaluated on three protein systems having different complexities of motion on the PES: calmodulin, adenylate kinase, and bacterial ferric binding protein. We find that the method captures the target conformation, while providing key residues and the optimum paths with relatively low free energy profiles.
Collapse
Affiliation(s)
- Farzaneh Jalalypour
- Faculty of Engineering and Natural Sciences, Sabanci University, 34956, Istanbul, Turkey
| | - Ozge Sensoy
- School of Engineering and Natural Sciences, Istanbul Medipol University, 34810, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, 34956, Istanbul, Turkey.,Sabanci University Nanotechnology Research and Application Center, SUNUM, 34956, Istanbul, Turkey
| |
Collapse
|
11
|
Liang Z, Zhu Y, Long J, Ye F, Hu G. Both intra and inter-domain interactions define the intrinsic dynamics and allosteric mechanism in DNMT1s. Comput Struct Biotechnol J 2020; 18:749-764. [PMID: 32280430 PMCID: PMC7132064 DOI: 10.1016/j.csbj.2020.03.016] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 03/17/2020] [Accepted: 03/17/2020] [Indexed: 01/03/2023] Open
Abstract
Dynamics and allosteric potentials of the RFTS domain are proposed. Hinge sites located at the RFTS-CD interface are key regulators for inter-domain interactions. Network analysis reveals local allosteric networks and inter-domain communication pathways in DNMT1. A potential allosteric site at the TRD interface for DNMT1 is identified.
DNA methyltransferase 1 (DNMT1), a large multidomain enzyme, is believed to be involved in the passive transmission of genomic methylation patterns via methylation maintenance. Yet, the molecular mechanism of interaction networks underlying DNMT1 structures, dynamics, and its biological significance has yet to be fully characterized. In this work, we used an integrated computational strategy that combined coarse-grained and atomistic simulations with coevolution information and network modeling of the residue interactions for the systematic investigation of allosteric dynamics in DNMT1. The elastic network modeling has proposed that the high plasticity of RFTS has strengthened the correlated behaviors of DNMT1 structures through the hinge sites located at the RFTS-CD interface, which mediate the collective motions between domains. The perturbation response scanning (PRS) analysis combined with the enrichment analysis of disease mutations have further highlighted the allosteric potential of the RFTS domain. Furthermore, the long-range paths connect the intra-domain interactions through the TRD interface and catalytic interface, emphasizing some key inter-domain interactions as the bridges in the global allosteric regulation of DNMT1. The observed interplay between conserved intra-domain networks and dynamical plasticity encoded by inter-domain interactions provides insights into the intrinsic dynamics and functional evolution, as well as the design of allosteric modulators of DNMT1 based on the TRD interface.
Collapse
Affiliation(s)
- Zhongjie Liang
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Yu Zhu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Jie Long
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Fei Ye
- College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| |
Collapse
|
12
|
Understanding the Pyrimethamine Drug Resistance Mechanism via Combined Molecular Dynamics and Dynamic Residue Network Analysis. Molecules 2020; 25:molecules25040904. [PMID: 32085470 PMCID: PMC7070769 DOI: 10.3390/molecules25040904] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2019] [Revised: 02/11/2020] [Accepted: 02/14/2020] [Indexed: 02/08/2023] Open
Abstract
In this era of precision medicine, insights into the resistance mechanism of drugs are integral for the development of potent therapeutics. Here, we sought to understand the contribution of four point mutations (N51I, C59R, S108N, and I164L) within the active site of the malaria parasite enzyme dihydrofolate reductase (DHFR) towards the resistance of the antimalarial drug pyrimethamine. Homology modeling was used to obtain full-length models of wild type (WT) and mutant DHFR. Molecular docking was employed to dock pyrimethamine onto the generated structures. Subsequent all-atom molecular dynamics (MD) simulations and binding free-energy computations highlighted that pyrimethamine's stability and affinity inversely relates to the number of mutations within its binding site and, hence, resistance severity. Generally, mutations led to reduced binding affinity to pyrimethamine and increased conformational plasticity of DHFR. Next, dynamic residue network analysis (DRN) was applied to determine the impact of mutations and pyrimethamine binding on communication dispositions of DHFR residues. DRN revealed residues with distinctive communication profiles, distinguishing WT from drug-resistant mutants as well as pyrimethamine-bound from pyrimethamine-free models. Our results provide a new perspective on the understanding of mutation-induced drug resistance.
Collapse
|
13
|
Nussinov R, Tsai C, Jang H. Autoinhibition can identify rare driver mutations and advise pharmacology. FASEB J 2019; 34:16-29. [DOI: 10.1096/fj.201901341r] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2019] [Revised: 09/18/2019] [Accepted: 10/09/2019] [Indexed: 12/16/2022]
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section Basic Science Program Frederick National Laboratory for Cancer Research Frederick MD USA
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine Tel Aviv University Tel Aviv Israel
| | - Chung‐Jung Tsai
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine Tel Aviv University Tel Aviv Israel
| | - Hyunbum Jang
- Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine Tel Aviv University Tel Aviv Israel
| |
Collapse
|
14
|
Nussinov R, Tsai CJ, Jang H. Does Ras Activate Raf and PI3K Allosterically? Front Oncol 2019; 9:1231. [PMID: 31799192 PMCID: PMC6874141 DOI: 10.3389/fonc.2019.01231] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Accepted: 10/28/2019] [Indexed: 12/11/2022] Open
Abstract
The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP2 binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.
Collapse
Affiliation(s)
- Ruth Nussinov
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States.,Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Chung-Jung Tsai
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States
| | - Hyunbum Jang
- Cancer and Inflammation Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, National Cancer Institute at Frederick, Frederick, MD, United States
| |
Collapse
|
15
|
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
|
16
|
Novel and Improved Crystal Structures of H. influenzae, E. coli and P. aeruginosa Penicillin-Binding Protein 3 (PBP3) and N. gonorrhoeae PBP2: Toward a Better Understanding of β-Lactam Target-Mediated Resistance. J Mol Biol 2019; 431:3501-3519. [PMID: 31301409 DOI: 10.1016/j.jmb.2019.07.010] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2019] [Revised: 07/02/2019] [Accepted: 07/02/2019] [Indexed: 01/26/2023]
Abstract
Even with the emergence of antibiotic resistance, penicillin and the wider family of β-lactams have remained the single most important family of antibiotics. The periplasmic/extra-cytoplasmic targets of penicillin are a family of enzymes with a highly conserved catalytic activity involved in the final stage of bacterial cell wall (peptidoglycan) biosynthesis. Named after their ability to bind penicillin, rather than their catalytic activity, these key targets are called penicillin-binding proteins (PBPs). Resistance is predominantly mediated by reducing the target drug concentration via β-lactamases; however, naturally transformable bacteria have also acquired target-mediated resistance by inter-species recombination. Here we focus on structural based interpretations of amino acid alterations associated with the emergence of resistance within clinical isolates and include new PBP3 structures along with new, and improved, PBP-β-lactam co-structures.
Collapse
|
17
|
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
|
18
|
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
|
19
|
Nussinov R, Tsai CJ, Shehu A, Jang H. Computational Structural Biology: Successes, Future Directions, and Challenges. Molecules 2019; 24:molecules24030637. [PMID: 30759724 PMCID: PMC6384756 DOI: 10.3390/molecules24030637] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2019] [Revised: 02/05/2019] [Accepted: 02/10/2019] [Indexed: 02/06/2023] Open
Abstract
Computational biology has made powerful advances. Among these, trends in human health have been uncovered through heterogeneous 'big data' integration, and disease-associated genes were identified and classified. Along a different front, the dynamic organization of chromatin is being elucidated to gain insight into the fundamental question of genome regulation. Powerful conformational sampling methods have also been developed to yield a detailed molecular view of cellular processes. when combining these methods with the advancements in the modeling of supramolecular assemblies, including those at the membrane, we are finally able to get a glimpse into how cells' actions are regulated. Perhaps most intriguingly, a major thrust is on to decipher the mystery of how the brain is coded. Here, we aim to provide a broad, yet concise, sketch of modern aspects of computational biology, with a special focus on computational structural biology. We attempt to forecast the areas that computational structural biology will embrace in the future and the challenges that it may face. We skirt details, highlight successes, note failures, and map directions.
Collapse
Affiliation(s)
- Ruth Nussinov
- Computational Structural Biology Section, Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA.
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel.
| | - Chung-Jung Tsai
- Computational Structural Biology Section, Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA.
| | - Amarda Shehu
- Departments of Computer Science, Department of Bioengineering, and School of Systems Biology, George Mason University, Fairfax, VA 22030, USA.
| | - Hyunbum Jang
- Computational Structural Biology Section, Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA.
| |
Collapse
|
20
|
Precision medicine review: rare driver mutations and their biophysical classification. Biophys Rev 2019; 11:5-19. [PMID: 30610579 PMCID: PMC6381362 DOI: 10.1007/s12551-018-0496-2] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2018] [Accepted: 12/18/2018] [Indexed: 02/07/2023] Open
Abstract
How can biophysical principles help precision medicine identify rare driver mutations? A major tenet of pragmatic approaches to precision oncology and pharmacology is that driver mutations are very frequent. However, frequency is a statistical attribute, not a mechanistic one. Rare mutations can also act through the same mechanism, and as we discuss below, “latent driver” mutations may also follow the same route, with “helper” mutations. Here, we review how biophysics provides mechanistic guidelines that extend precision medicine. We outline principles and strategies, especially focusing on mutations that drive cancer. Biophysics has contributed profoundly to deciphering biological processes. However, driven by data science, precision medicine has skirted some of its major tenets. Data science embodies genomics, tissue- and cell-specific expression levels, making it capable of defining genome- and systems-wide molecular disease signatures. It classifies cancer driver genes/mutations and affected pathways, and its associated protein structural data guide drug discovery. Biophysics complements data science. It considers structures and their heterogeneous ensembles, explains how mutational variants can signal through distinct pathways, and how allo-network drugs can be harnessed. Biophysics clarifies how one mutation—frequent or rare—can affect multiple phenotypic traits by populating conformations that favor interactions with other network modules. It also suggests how to identify such mutations and their signaling consequences. Biophysics offers principles and strategies that can help precision medicine push the boundaries to transform our insight into biological processes and the practice of personalized medicine. By contrast, “phenotypic drug discovery,” which capitalizes on physiological cellular conditions and first-in-class drug discovery, may not capture the proper molecular variant. This is because variants of the same protein can express more than one phenotype, and a phenotype can be encoded by several variants.
Collapse
|