1
|
Kotev M, Diaz Gonzalez C. Molecular Dynamics and Other HPC Simulations for Drug Discovery. Methods Mol Biol 2024; 2716:265-291. [PMID: 37702944 DOI: 10.1007/978-1-0716-3449-3_12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
High performance computing (HPC) is taking an increasingly important place in drug discovery. It makes possible the simulation of complex biochemical systems with high precision in a short time, thanks to the use of sophisticated algorithms. It promotes the advancement of knowledge in fields that are inaccessible or difficult to access through experimentation and it contributes to accelerating the discovery of drugs for unmet medical needs while reducing costs. Herein, we report how computational performance has evolved over the past years, and then we detail three domains where HPC is essential. Molecular dynamics (MD) is commonly used to explore the flexibility of proteins, thus generating a better understanding of different possible approaches to modulate their activity. Modeling and simulation of biopolymer complexes enables the study of protein-protein interactions (PPI) in healthy and disease states, thus helping the identification of targets of pharmacological interest. Virtual screening (VS) also benefits from HPC to predict in a short time, among millions or billions of virtual chemical compounds, the best potential ligands that will be tested in relevant assays to start a rational drug design process.
Collapse
Affiliation(s)
- Martin Kotev
- Evotec SE, Integrated Drug Discovery, Molecular Architects, Campus Curie, Toulouse, France
| | | |
Collapse
|
2
|
Ozden B, Kryshtafovych A, Karaca E. The impact of AI-based modeling on the accuracy of protein assembly prediction: Insights from CASP15. Proteins 2023; 91:1636-1657. [PMID: 37861057 PMCID: PMC10873090 DOI: 10.1002/prot.26598] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 09/12/2023] [Accepted: 09/14/2023] [Indexed: 10/21/2023]
Abstract
In CASP15, 87 predictors submitted around 11 000 models on 41 assembly targets. The community demonstrated exceptional performance in overall fold and interface contact predictions, achieving an impressive success rate of 90% (compared to 31% in CASP14). This remarkable accomplishment is largely due to the incorporation of DeepMind's AF2-Multimer approach into custom-built prediction pipelines. To evaluate the added value of participating methods, we compared the community models to the baseline AF2-Multimer predictor. In over 1/3 of cases, the community models were superior to the baseline predictor. The main reasons for this improved performance were the use of custom-built multiple sequence alignments, optimized AF2-Multimer sampling, and the manual assembly of AF2-Multimer-built subcomplexes. The best three groups, in order, are Zheng, Venclovas, and Wallner. Zheng and Venclovas reached a 73.2% success rate over all (41) cases, while Wallner attained 69.4% success rate over 36 cases. Nonetheless, challenges remain in predicting structures with weak evolutionary signals, such as nanobody-antigen, antibody-antigen, and viral complexes. Expectedly, modeling large complexes also remains challenging due to their high memory compute demands. In addition to the assembly category, we assessed the accuracy of modeling interdomain interfaces in the tertiary structure prediction targets. Models on seven targets featuring 17 unique interfaces were analyzed. Best predictors achieved a 76.5% success rate, with the UM-TBM group being the leader. In the interdomain category, we observed that the predictors faced challenges, as in the case of the assembly category, when the evolutionary signal for a given domain pair was weak or the structure was large. Overall, CASP15 witnessed unprecedented improvement in interface modeling, reflecting the AI revolution seen in CASP14.
Collapse
Affiliation(s)
- Burcu Ozden
- Izmir Biomedicine and Genome Center, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
| | - Andriy Kryshtafovych
- Protein Structure Prediction Center, Genome and Biomedical Sciences Facilities, University of California, Davis, California, USA
| | - Ezgi Karaca
- Izmir Biomedicine and Genome Center, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
| |
Collapse
|
3
|
Ozden B, Kryshtafovych A, Karaca E. The Impact of AI-Based Modeling on the Accuracy of Protein Assembly Prediction: Insights from CASP15. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.10.548341. [PMID: 37503072 PMCID: PMC10369898 DOI: 10.1101/2023.07.10.548341] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
In CASP15, 87 predictors submitted around 11,000 models on 41 assembly targets. The community demonstrated exceptional performance in overall fold and interface contact prediction, achieving an impressive success rate of 90% (compared to 31% in CASP14). This remarkable accomplishment is largely due to the incorporation of DeepMind's AF2-Multimer approach into custom-built prediction pipelines. To evaluate the added value of participating methods, we compared the community models to the baseline AF2-Multimer predictor. In over 1/3 of cases the community models were superior to the baseline predictor. The main reasons for this improved performance were the use of custom-built multiple sequence alignments, optimized AF2-Multimer sampling, and the manual assembly of AF2-Multimer-built subcomplexes. The best three groups, in order, are Zheng, Venclovas and Wallner. Zheng and Venclovas reached a 73.2% success rate over all (41) cases, while Wallner attained 69.4% success rate over 36 cases. Nonetheless, challenges remain in predicting structures with weak evolutionary signals, such as nanobody-antigen, antibody-antigen, and viral complexes. Expectedly, modeling large complexes remains also challenging due to their high memory compute demands. In addition to the assembly category, we assessed the accuracy of modeling interdomain interfaces in the tertiary structure prediction targets. Models on seven targets featuring 17 unique interfaces were analyzed. Best predictors achieved the 76.5% success rate, with the UM-TBM group being the leader. In the interdomain category, we observed that the predictors faced challenges, as in the case of the assembly category, when the evolutionary signal for a given domain pair was weak or the structure was large. Overall, CASP15 witnessed unprecedented improvement in interface modeling, reflecting the AI revolution seen in CASP14.
Collapse
Affiliation(s)
- Burcu Ozden
- Izmir Biomedicine and Genome Center, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
| | - Andriy Kryshtafovych
- Protein Structure Prediction Center, Genome and Biomedical Sciences Facilities, University of California, Davis, California, USA
| | - Ezgi Karaca
- Izmir Biomedicine and Genome Center, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
| |
Collapse
|
4
|
Kurisaki I, Suzuki M. Simulation toolkits at the molecular scale for trans-scale thermal signaling. Comput Struct Biotechnol J 2023; 21:2547-2557. [PMID: 37102156 PMCID: PMC10123322 DOI: 10.1016/j.csbj.2023.03.040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 03/22/2023] [Accepted: 03/23/2023] [Indexed: 04/28/2023] Open
Abstract
Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.
Collapse
Affiliation(s)
- Ikuo Kurisaki
- Waseda Research Institute for Science and Engineering, Waseda University, Bldg. No.55, S Tower, 4th Floor, 3–4-1 Okubo Shinjuku-ku, Tokyo 169–8555, Japan
- Corresponding authors.
| | - Madoka Suzuki
- Institute for Protein Research, Osaka University, 3–2 Yamadaoka, Suita, Osaka 565–0871, Japan
- Corresponding authors.
| |
Collapse
|
5
|
Launay R, Teppa E, Esque J, André I. Modeling Protein Complexes and Molecular Assemblies Using Computational Methods. Methods Mol Biol 2023; 2553:57-77. [PMID: 36227539 DOI: 10.1007/978-1-0716-2617-7_4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Many biological molecules are assembled into supramolecular complexes that are necessary to perform functions in the cell. Better understanding and characterization of these molecular assemblies are thus essential to further elucidate molecular mechanisms and key protein-protein interactions that could be targeted to modulate the protein binding affinity or develop new binders. Experimental access to structural information on these supramolecular assemblies is often hampered by the size of these systems that make their recombinant production and characterization rather difficult. Computational methods combining both structural data, molecular modeling techniques, and sequence coevolution information can thus offer a good alternative to gain access to the structural organization of protein complexes and assemblies. Herein, we present some computational methods to predict structural models of the protein partners, to search for interacting regions using coevolution information, and to build molecular assemblies. The approach is exemplified using a case study to model the succinate-quinone oxidoreductase heterocomplex.
Collapse
Affiliation(s)
- Romain Launay
- Toulouse Biotechnology Institute, TBI, Université de Toulouse, CNRS, INRAE, INSA, Toulouse Cedex 04, France
| | - Elin Teppa
- Toulouse Biotechnology Institute, TBI, Université de Toulouse, CNRS, INRAE, INSA, Toulouse Cedex 04, France
| | - Jérémy Esque
- Toulouse Biotechnology Institute, TBI, Université de Toulouse, CNRS, INRAE, INSA, Toulouse Cedex 04, France.
| | - Isabelle André
- Toulouse Biotechnology Institute, TBI, Université de Toulouse, CNRS, INRAE, INSA, Toulouse Cedex 04, France.
| |
Collapse
|
6
|
Drake ZC, Seffernick JT, Lindert S. Protein complex prediction using Rosetta, AlphaFold, and mass spectrometry covalent labeling. Nat Commun 2022; 13:7846. [PMID: 36543826 PMCID: PMC9772387 DOI: 10.1038/s41467-022-35593-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Accepted: 12/09/2022] [Indexed: 12/24/2022] Open
Abstract
Covalent labeling (CL) in combination with mass spectrometry can be used as an analytical tool to study and determine structural properties of protein-protein complexes. However, data from these experiments is sparse and does not unambiguously elucidate protein structure. Thus, computational algorithms are needed to deduce structure from the CL data. In this work, we present a hybrid method that combines models of protein complex subunits generated with AlphaFold with differential CL data via a CL-guided protein-protein docking in Rosetta. In a benchmark set, the RMSD (root-mean-square deviation) of the best-scoring models was below 3.6 Å for 5/5 complexes with inclusion of CL data, whereas the same quality was only achieved for 1/5 complexes without CL data. This study suggests that our integrated approach can successfully use data obtained from CL experiments to distinguish between nativelike and non-nativelike models.
Collapse
Affiliation(s)
- Zachary C Drake
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, 43210, US
| | - Justin T Seffernick
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, 43210, US
| | - Steffen Lindert
- Department of Chemistry and Biochemistry, Ohio State University, Columbus, OH, 43210, US.
| |
Collapse
|
7
|
Bryant P, Pozzati G, Zhu W, Shenoy A, Kundrotas P, Elofsson A. Predicting the structure of large protein complexes using AlphaFold and Monte Carlo tree search. Nat Commun 2022; 13:6028. [PMID: 36224222 PMCID: PMC9556563 DOI: 10.1038/s41467-022-33729-4] [Citation(s) in RCA: 69] [Impact Index Per Article: 34.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 09/29/2022] [Indexed: 11/30/2022] Open
Abstract
AlphaFold can predict the structure of single- and multiple-chain proteins with very high accuracy. However, the accuracy decreases with the number of chains, and the available GPU memory limits the size of protein complexes which can be predicted. Here we show that one can predict the structure of large complexes starting from predictions of subcomponents. We assemble 91 out of 175 complexes with 10–30 chains from predicted subcomponents using Monte Carlo tree search, with a median TM-score of 0.51. There are 30 highly accurate complexes (TM-score ≥0.8, 33% of complete assemblies). We create a scoring function, mpDockQ, that can distinguish if assemblies are complete and predict their accuracy. We find that complexes containing symmetry are accurately assembled, while asymmetrical complexes remain challenging. The method is freely available and accesible as a Colab notebook https://colab.research.google.com/github/patrickbryant1/MoLPC/blob/master/MoLPC.ipynb. The accuracy of AlphaFold decreases with the number of protein chains and the available GPU memory limits the size of protein complexes that can be predicted. Here, the authors show that complexes with 10–30 chains can be assembled from predicted subcomponents using Monte Carlo tree search.
Collapse
Affiliation(s)
- Patrick Bryant
- Science for Life Laboratory, 172 21, Solna, Sweden. .,Department of Biochemistry and Biophysics, Stockholm University, 106 91, Stockholm, Sweden.
| | - Gabriele Pozzati
- Science for Life Laboratory, 172 21, Solna, Sweden.,Department of Biochemistry and Biophysics, Stockholm University, 106 91, Stockholm, Sweden
| | - Wensi Zhu
- Science for Life Laboratory, 172 21, Solna, Sweden.,Department of Biochemistry and Biophysics, Stockholm University, 106 91, Stockholm, Sweden
| | - Aditi Shenoy
- Science for Life Laboratory, 172 21, Solna, Sweden.,Department of Biochemistry and Biophysics, Stockholm University, 106 91, Stockholm, Sweden
| | - Petras Kundrotas
- Science for Life Laboratory, 172 21, Solna, Sweden.,Center for Computational Biology, The University of Kansas, Lawrence, KS, 66047, USA
| | - Arne Elofsson
- Science for Life Laboratory, 172 21, Solna, Sweden.,Department of Biochemistry and Biophysics, Stockholm University, 106 91, Stockholm, Sweden
| |
Collapse
|
8
|
Sen N, Madhusudhan MS. A structural database of chain–chain and domain–domain interfaces of proteins. Protein Sci 2022. [DOI: 10.1002/pro.4406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Neeladri Sen
- Indian Institute of Science Education and Research Pune India
- Institute of Structural and Molecular Biology University College London London UK
| | | |
Collapse
|
9
|
Begum R, Thota S, Abdulkadir A, Kaur G, Bagam P, Batra S. NADPH oxidase family proteins: signaling dynamics to disease management. Cell Mol Immunol 2022; 19:660-686. [PMID: 35585127 DOI: 10.1038/s41423-022-00858-1] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Accepted: 03/12/2022] [Indexed: 12/16/2022] Open
Abstract
Reactive oxygen species (ROS) are pervasive signaling molecules in biological systems. In humans, a lack of ROS causes chronic and extreme bacterial infections, while uncontrolled release of these factors causes pathologies due to excessive inflammation. Professional phagocytes such as neutrophils (PMNs), eosinophils, monocytes, and macrophages use superoxide-generating NADPH oxidase (NOX) as part of their arsenal of antimicrobial mechanisms to produce high levels of ROS. NOX is a multisubunit enzyme complex composed of five essential subunits, two of which are localized in the membrane, while three are localized in the cytosol. In resting phagocytes, the oxidase complex is unassembled and inactive; however, it becomes activated after cytosolic components translocate to the membrane and are assembled into a functional oxidase. The NOX isoforms play a variety of roles in cellular differentiation, development, proliferation, apoptosis, cytoskeletal control, migration, and contraction. Recent studies have identified NOX as a major contributor to disease pathologies, resulting in a shift in focus on inhibiting the formation of potentially harmful free radicals. Therefore, a better understanding of the molecular mechanisms and the transduction pathways involved in NOX-mediated signaling is essential for the development of new therapeutic agents that minimize the hyperproduction of ROS. The current review provides a thorough overview of the various NOX enzymes and their roles in disease pathophysiology, highlights pharmacological strategies, and discusses the importance of computational modeling for future NOX-related studies.
Collapse
Affiliation(s)
- Rizwana Begum
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA
| | - Shilpa Thota
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA
| | - Abubakar Abdulkadir
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA
| | - Gagandeep Kaur
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA.,Department of Environmental Medicine, University of Rochester, School of Medicine and Dentistry, Rochester, NY, 14642, USA
| | - Prathyusha Bagam
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA.,Division of Systems Biology, National Center for Toxicological Research, Jefferson, AR, 72079, USA
| | - Sanjay Batra
- Laboratory of Pulmonary Immunotoxicology, Department of Environmental Toxicology, Southern University and A&M College, Baton Rouge, LA, 70813, USA.
| |
Collapse
|
10
|
Kurisaki I, Tanaka S. Computational prediction of heteromeric protein complex disassembly order using hybrid Monte Carlo/molecular dynamics simulation. Phys Chem Chem Phys 2022; 24:10575-10587. [PMID: 35445673 DOI: 10.1039/d2cp00267a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The physicochemical entities comprising the biological phenomena in the cell form a network of biochemical reactions and the activity of such a network is regulated by multimeric protein complexes. Mass spectroscopy (MS) experiments and multimeric protein docking simulations based on structural bioinformatics techniques have revealed the molecular-level stoichiometry and static configuration of subcomplexes in their bound forms, thus revealing the subcomplex population and formation orders. Meanwhile, these methodologies are not designed to straightforwardly examine the temporal dynamics of multimeric protein assembly and disassembly, essential physicochemical properties to understand the functional expression mechanisms of proteins in the biological environment. To address this problem, we have developed an atomistic simulation in the framework of the hybrid Monte Carlo/molecular dynamics (hMC/MD) method and succeeded in observing the disassembly of a homomeric pentamer of the serum amyloid P component protein in an experimentally consistent order. In this study, we improved the hMC/MD method to examine the disassembly processes of the tryptophan synthase tetramer, a paradigmatic heteromeric protein complex in MS studies. We employed the likelihood-based selection scheme to determine a dissociation-prone subunit pair at every hMC/MD simulation cycle and achieved highly reliable predictions of the disassembly orders without a priori knowledge of the MS experiments and structural bioinformatics simulations. The success rate for the experimentally-observed disassembly order is over 0.9. We similarly succeeded in reliable predictions for three other tetrameric protein complexes. These achievements indicate the potential applicability of our hMC/MD approach as a general-purpose methodology to obtain microscopic and physicochemical insights into multimeric protein complex formation.
Collapse
Affiliation(s)
- Ikuo Kurisaki
- Department of Computational Science, Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan.
| | - Shigenori Tanaka
- Department of Computational Science, Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan.
| |
Collapse
|
11
|
Yang P, Ning K. How much metagenome data is needed for protein structure prediction: The advantages of targeted approach from the ecological and evolutionary perspectives. IMETA 2022; 1:e9. [PMID: 38867727 PMCID: PMC10989767 DOI: 10.1002/imt2.9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 12/23/2021] [Accepted: 01/04/2022] [Indexed: 06/14/2024]
Abstract
It has been proven that three-dimensional protein structures could be modeled by supplementing homologous sequences with metagenome sequences. Even though a large volume of metagenome data is utilized for such purposes, a significant proportion of proteins remain unsolved. In this review, we focus on identifying ecological and evolutionary patterns in metagenome data, decoding the complicated relationships of these patterns with protein structures, and investigating how these patterns can be effectively used to improve protein structure prediction. First, we proposed the metagenome utilization efficiency and marginal effect model to quantify the divergent distribution of homologous sequences for the protein family. Second, we proposed that the targeted approach effectively identifies homologous sequences from specified biomes compared with the untargeted approach's blind search. Finally, we determined the lower bound for metagenome data required for predicting all the protein structures in the Pfam database and showed that the present metagenome data is insufficient for this purpose. In summary, we discovered ecological and evolutionary patterns in the metagenome data that may be used to predict protein structures effectively. The targeted approach is promising in terms of effectively extracting homologous sequences and predicting protein structures using these patterns.
Collapse
Affiliation(s)
- Pengshuo Yang
- Key Laboratory of Molecular Biophysics of the Ministry of Education, Hubei Key Laboratory of Bioinformatics and Molecular‐Imaging, Department of Bioinformatics and Systems BiologyCenter of AI Biology, College of Life Science and Technology, Huazhong University of Science and TechnologyWuhanHubeiChina
| | - Kang Ning
- Key Laboratory of Molecular Biophysics of the Ministry of Education, Hubei Key Laboratory of Bioinformatics and Molecular‐Imaging, Department of Bioinformatics and Systems BiologyCenter of AI Biology, College of Life Science and Technology, Huazhong University of Science and TechnologyWuhanHubeiChina
| |
Collapse
|
12
|
Martino E, Chiarugi S, Margheriti F, Garau G. Mapping, Structure and Modulation of PPI. Front Chem 2021; 9:718405. [PMID: 34692637 PMCID: PMC8529325 DOI: 10.3389/fchem.2021.718405] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Accepted: 09/21/2021] [Indexed: 11/13/2022] Open
Abstract
Because of the key relevance of protein–protein interactions (PPI) in diseases, the modulation of protein-protein complexes is of relevant clinical significance. The successful design of binding compounds modulating PPI requires a detailed knowledge of the involved protein-protein system at molecular level, and investigation of the structural motifs that drive the association of the proteins at the recognition interface. These elements represent hot spots of the protein binding free energy, define the complex lifetime and possible modulation strategies. Here, we review the advanced technologies used to map the PPI involved in human diseases, to investigate the structure-function features of protein complexes, and to discover effective ligands that modulate the PPI for therapeutic intervention.
Collapse
Affiliation(s)
- Elisa Martino
- Laboratorio NEST, Scuola Normale Superiore, Pisa, Italy
| | - Sara Chiarugi
- Laboratorio NEST, Scuola Normale Superiore, Pisa, Italy.,BioStructures Lab, Istituto Italiano di Tecnologia (IIT@NEST), Pisa, Italy
| | | | - Gianpiero Garau
- BioStructures Lab, Istituto Italiano di Tecnologia (IIT@NEST), Pisa, Italy
| |
Collapse
|
13
|
Soto-Ospina A, Araque Marín P, Bedoya G, Sepulveda-Falla D, Villegas Lanau A. Protein Predictive Modeling and Simulation of Mutations of Presenilin-1 Familial Alzheimer's Disease on the Orthosteric Site. Front Mol Biosci 2021; 8:649990. [PMID: 34150846 PMCID: PMC8206637 DOI: 10.3389/fmolb.2021.649990] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Accepted: 04/22/2021] [Indexed: 11/13/2022] Open
Abstract
Alzheimer's disease pathology is characterized by β-amyloid plaques and neurofibrillary tangles. Amyloid precursor protein is processed by β and γ secretase, resulting in the production of β-amyloid peptides with a length ranging from 38 to 43 amino acids. Presenilin 1 (PS1) is the catalytic unit of γ-secretase, and more than 200 PS1 pathogenic mutations have been identified as causative for Alzheimer's disease. A complete monocrystal structure of PS1 has not been determined so far due to the presence of two flexible domains. We have developed a complete structural model of PS1 using a computational approach with structure prediction software. Missing fragments Met1-Glut72 and Ser290-Glu375 were modeled and validated by their energetic and stereochemical characteristics. Then, with the complete structure of PS1, we defined that these fragments do not have a direct effect in the structure of the pore. Next, we used our hypothetical model for the analysis of the functional effects of PS1 mutations Ala246GLu, Leu248Pro, Leu248Arg, Leu250Val, Tyr256Ser, Ala260Val, and Val261Phe, localized in the catalytic pore. For this, we used a quantum mechanics/molecular mechanics (QM/MM) hybrid method, evaluating modifications in the topology, potential surface density, and electrostatic potential map of mutated PS1 proteins. We found that each mutation exerts changes resulting in structural modifications of the active site and in the shape of the pore. We suggest this as a valid approach for functional studies of PS1 in view of the possible impact in substrate processing and for the design of targeted therapeutic strategies.
Collapse
Affiliation(s)
- Alejandro Soto-Ospina
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
| | - Pedronel Araque Marín
- School of Life Sciences, Research and Innovation in Chemistry Formulations Group, EIA University, Envigado, Colombia
| | - Gabriel Bedoya
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
| | - Diego Sepulveda-Falla
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
- Molecular Neuropathology of Alzheimer’s Disease, Institute of Neuropathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Andrés Villegas Lanau
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
| |
Collapse
|
14
|
Vedithi SC, Malhotra S, Acebrón-García-de-Eulate M, Matusevicius M, Torres PHM, Blundell TL. Structure-Guided Computational Approaches to Unravel Druggable Proteomic Landscape of Mycobacterium leprae. Front Mol Biosci 2021; 8:663301. [PMID: 34026836 PMCID: PMC8138464 DOI: 10.3389/fmolb.2021.663301] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Accepted: 04/12/2021] [Indexed: 02/02/2023] Open
Abstract
Leprosy, caused by Mycobacterium leprae (M. leprae), is treated with a multidrug regimen comprising Dapsone, Rifampicin, and Clofazimine. These drugs exhibit bacteriostatic, bactericidal and anti-inflammatory properties, respectively, and control the dissemination of infection in the host. However, the current treatment is not cost-effective, does not favor patient compliance due to its long duration (12 months) and does not protect against the incumbent nerve damage, which is a severe leprosy complication. The chronic infectious peripheral neuropathy associated with the disease is primarily due to the bacterial components infiltrating the Schwann cells that protect neuronal axons, thereby inducing a demyelinating phenotype. There is a need to discover novel/repurposed drugs that can act as short duration and effective alternatives to the existing treatment regimens, preventing nerve damage and consequent disability associated with the disease. Mycobacterium leprae is an obligate pathogen resulting in experimental intractability to cultivate the bacillus in vitro and limiting drug discovery efforts to repositioning screens in mouse footpad models. The dearth of knowledge related to structural proteomics of M. leprae, coupled with emerging antimicrobial resistance to all the three drugs in the multidrug therapy, poses a need for concerted novel drug discovery efforts. A comprehensive understanding of the proteomic landscape of M. leprae is indispensable to unravel druggable targets that are essential for bacterial survival and predilection of human neuronal Schwann cells. Of the 1,614 protein-coding genes in the genome of M. leprae, only 17 protein structures are available in the Protein Data Bank. In this review, we discussed efforts made to model the proteome of M. leprae using a suite of software for protein modeling that has been developed in the Blundell laboratory. Precise template selection by employing sequence-structure homology recognition software, multi-template modeling of the monomeric models and accurate quality assessment are the hallmarks of the modeling process. Tools that map interfaces and enable building of homo-oligomers are discussed in the context of interface stability. Other software is described to determine the druggable proteome by using information related to the chokepoint analysis of the metabolic pathways, gene essentiality, homology to human proteins, functional sites, druggable pockets and fragment hotspot maps.
Collapse
Affiliation(s)
- Sundeep Chaitanya Vedithi
- Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom,*Correspondence: Sundeep Chaitanya Vedithi,
| | - Sony Malhotra
- Rutherford Appleton Laboratory, Science and Technology Facilities Council, Oxon, United Kingdom
| | | | | | - Pedro Henrique Monteiro Torres
- Laboratório de Modelagem e Dinâmica Molecular, Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Tom L. Blundell
- Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom,Tom L. Blundell,
| |
Collapse
|
15
|
Kurisaki I, Tanaka S. Reaction Pathway Sampling and Free-Energy Analyses for Multimeric Protein Complex Disassembly by Employing Hybrid Configuration Bias Monte Carlo/Molecular Dynamics Simulation. ACS OMEGA 2021; 6:4749-4758. [PMID: 33644582 PMCID: PMC7905796 DOI: 10.1021/acsomega.0c05579] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Accepted: 01/27/2021] [Indexed: 05/08/2023]
Abstract
Physicochemical characterization of multimeric biomacromolecule assembly and disassembly processes is a milestone to understand the mechanisms for biological phenomena at the molecular level. Mass spectroscopy (MS) and structural bioinformatics (SB) approaches have become feasible to identify subcomplexes involved in assembly and disassembly, while they cannot provide atomic information sufficient for free-energy calculation to characterize transition mechanism between two different sets of subcomplexes. To combine observations derived from MS and SB approaches with conventional free-energy calculation protocols, we here designed a new reaction pathway sampling method by employing hybrid configuration bias Monte Carlo/molecular dynamics (hcbMC/MD) scheme and applied it to simulate the disassembly process of serum amyloid P component (SAP) pentamer. The results we obtained are consistent with those of the earlier MS and SB studies with respect to SAP subcomplex species and the initial stage of SAP disassembly processes. Furthermore, we observed a novel dissociation event, ring-opening reaction of SAP pentamer. Employing free-energy calculation combined with the hcbMC/MD reaction pathway trajectories, we moreover obtained experimentally testable observations on (1) reaction time of the ring-opening reaction and (2) importance of Asp42 and Lys117 for stable formation of SAP oligomer.
Collapse
|
16
|
Das S, Chakrabarti S. Classification and prediction of protein-protein interaction interface using machine learning algorithm. Sci Rep 2021; 11:1761. [PMID: 33469042 PMCID: PMC7815773 DOI: 10.1038/s41598-020-80900-2] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2020] [Accepted: 12/15/2020] [Indexed: 01/29/2023] Open
Abstract
Structural insight of the protein-protein interaction (PPI) interface can provide knowledge about the kinetics, thermodynamics and molecular functions of the complex while elucidating its role in diseases and further enabling it as a potential therapeutic target. However, owing to experimental lag in solving protein-protein complex structures, three-dimensional (3D) knowledge of the PPI interfaces can be gained via computational approaches like molecular docking and post-docking analyses. Despite development of numerous docking tools and techniques, success in identification of native like interfaces based on docking score functions is limited. Hence, we employed an in-depth investigation of the structural features of the interface that might successfully delineate native complexes from non-native ones. We identify interface properties, which show statistically significant difference between native and non-native interfaces belonging to homo and hetero, protein-protein complexes. Utilizing these properties, a support vector machine (SVM) based classification scheme has been implemented to differentiate native and non-native like complexes generated using docking decoys. Benchmarking and comparative analyses suggest very good performance of our SVM classifiers. Further, protein interactions, which are proven via experimental findings but not resolved structurally, were subjected to this approach where 3D-models of the complexes were generated and most likely interfaces were predicted. A web server called Protein Complex Prediction by Interface Properties (PCPIP) is developed to predict whether interface of a given protein-protein dimer complex resembles known protein interfaces. The server is freely available at http://www.hpppi.iicb.res.in/pcpip/ .
Collapse
Affiliation(s)
- Subhrangshu Das
- grid.417635.20000 0001 2216 5074Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, WB India
| | - Saikat Chakrabarti
- grid.417635.20000 0001 2216 5074Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, Kolkata, WB India
| |
Collapse
|
17
|
Abstract
Biological processes are often mediated by complexes formed between proteins and various biomolecules. The 3D structures of such protein-biomolecule complexes provide insights into the molecular mechanism of their action. The structure of these complexes can be predicted by various computational methods. Choosing an appropriate method for modelling depends on the category of biomolecule that a protein interacts with and the availability of structural information about the protein and its interacting partner. We intend for the contents of this chapter to serve as a guide as to what software would be the most appropriate for the type of data at hand and the kind of 3D complex structure required. Particularly, we have dealt with protein-small molecule ligand, protein-peptide, protein-protein, and protein-nucleic acid interactions.Most, if not all, model building protocols perform some sampling and scoring. Typically, several alternate conformations and configurations of the interactors are sampled. Each such sample is then scored for optimization. To boost the confidence in these predicted models, their assessment using other independent scoring schemes besides the inbuilt/default ones would prove to be helpful. This chapter also lists such software and serves as a guide to gauge the fidelity of modelled structures of biomolecular complexes.
Collapse
|
18
|
Dhawanjewar AS, Roy AA, Madhusudhan MS. A knowledge-based scoring function to assess quaternary associations of proteins. Bioinformatics 2020; 36:3739-3748. [PMID: 32246820 DOI: 10.1093/bioinformatics/btaa207] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 03/01/2020] [Accepted: 03/30/2020] [Indexed: 12/21/2022] Open
Abstract
MOTIVATION The elucidation of all inter-protein interactions would significantly enhance our knowledge of cellular processes at a molecular level. Given the enormity of the problem, the expenses and limitations of experimental methods, it is imperative that this problem is tackled computationally. In silico predictions of protein interactions entail sampling different conformations of the purported complex and then scoring these to assess for interaction viability. In this study, we have devised a new scheme for scoring protein-protein interactions. RESULTS Our method, PIZSA (Protein Interaction Z-Score Assessment), is a binary classification scheme for identification of native protein quaternary assemblies (binders/nonbinders) based on statistical potentials. The scoring scheme incorporates residue-residue contact preference on the interface with per residue-pair atomic contributions and accounts for clashes. PIZSA can accurately discriminate between native and non-native structural conformations from protein docking experiments and outperform other contact-based potential scoring functions. The method has been extensively benchmarked and is among the top 6 methods, outperforming 31 other statistical, physics based and machine learning scoring schemes. The PIZSA potentials can also distinguish crystallization artifacts from biological interactions. AVAILABILITY AND IMPLEMENTATION PIZSA is implemented as a web server at http://cospi.iiserpune.ac.in/pizsa and can be downloaded as a standalone package from http://cospi.iiserpune.ac.in/pizsa/Download/Download.html. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Abhilesh S Dhawanjewar
- Indian Institute of Science Education and Research, Pashan, Pune 411008, India.,School of Biological Sciences, University of Nebraska, Lincoln, NE 68588, USA
| | - Ankit A Roy
- Indian Institute of Science Education and Research, Pashan, Pune 411008, India
| | | |
Collapse
|
19
|
Roy AA, Dhawanjewar AS, Sharma P, Singh G, Madhusudhan MS. Protein Interaction Z Score Assessment (PIZSA): an empirical scoring scheme for evaluation of protein-protein interactions. Nucleic Acids Res 2020; 47:W331-W337. [PMID: 31114890 PMCID: PMC6602501 DOI: 10.1093/nar/gkz368] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Revised: 04/24/2019] [Accepted: 05/15/2019] [Indexed: 11/24/2022] Open
Abstract
Our web server, PIZSA (http://cospi.iiserpune.ac.in/pizsa), assesses the likelihood of protein–protein interactions by assigning a Z Score computed from interface residue contacts. Our score takes into account the optimal number of atoms that mediate the interaction between pairs of residues and whether these contacts emanate from the main chain or side chain. We tested the score on 174 native interactions for which 100 decoys each were constructed using ZDOCK. The native structure scored better than any of the decoys in 146 cases and was able to rank within the 95th percentile in 162 cases. This easily outperforms a competing method, CIPS. We also benchmarked our scoring scheme on 15 targets from the CAPRI dataset and found that our method had results comparable to that of CIPS. Further, our method is able to analyse higher order protein complexes without the need to explicitly identify chains as receptors or ligands. The PIZSA server is easy to use and could be used to score any input three-dimensional structure and provide a residue pair-wise break up of the results. Attractively, our server offers a platform for users to upload their own potentials and could serve as an ideal testing ground for this class of scoring schemes.
Collapse
Affiliation(s)
- Ankit A Roy
- Indian Institute of Science Education and Research, Pune, Dr Homi Bhabha Road, Pashan, Pune 411008, India
| | - Abhilesh S Dhawanjewar
- Indian Institute of Science Education and Research, Pune, Dr Homi Bhabha Road, Pashan, Pune 411008, India.,presently at School of Biological Sciences, University of Nebraska, Lincoln, NE 68588, USA
| | - Parichit Sharma
- Indian Institute of Science Education and Research, Pune, Dr Homi Bhabha Road, Pashan, Pune 411008, India.,presently at School of Informatics, Computing & Engineering, Department of Computer Science, Indiana University, Bloomington, IN 47408, USA
| | - Gulzar Singh
- Indian Institute of Science Education and Research, Pune, Dr Homi Bhabha Road, Pashan, Pune 411008, India
| | - M S Madhusudhan
- Indian Institute of Science Education and Research, Pune, Dr Homi Bhabha Road, Pashan, Pune 411008, India
| |
Collapse
|
20
|
Andreani J, Quignot C, Guerois R. Structural prediction of protein interactions and docking using conservation and coevolution. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1470] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Jessica Andreani
- Université Paris‐Saclay CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC) Gif‐sur‐Yvette France
| | - Chloé Quignot
- Université Paris‐Saclay CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC) Gif‐sur‐Yvette France
| | - Raphael Guerois
- Université Paris‐Saclay CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC) Gif‐sur‐Yvette France
| |
Collapse
|
21
|
Koukos P, Bonvin A. Integrative Modelling of Biomolecular Complexes. J Mol Biol 2020; 432:2861-2881. [DOI: 10.1016/j.jmb.2019.11.009] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2019] [Revised: 11/12/2019] [Accepted: 11/13/2019] [Indexed: 12/31/2022]
|
22
|
Kumar AP, Verma CS, Lukman S. Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design. Brief Bioinform 2020; 22:270-287. [PMID: 31950981 DOI: 10.1093/bib/bbz161] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Revised: 08/29/2019] [Accepted: 11/15/2019] [Indexed: 01/09/2023] Open
Abstract
Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.
Collapse
Affiliation(s)
- Ammu Prasanna Kumar
- Department of Chemistry, College of Arts and Sciences, Khalifa University, Abu Dhabi, United Arab Emirates.,Research Unit in Bioinformatics, Department of Biochemistry and Microbiology, Rhodes University, South Africa
| | - Chandra S Verma
- Bioinformatics Institute, Agency for Science, Technology and Research, Singapore.,Department of Biological Sciences, National University of Singapore, Singapore.,School of Biological Sciences, Nanyang Technological University, Singapore
| | - Suryani Lukman
- Department of Chemistry, College of Arts and Sciences, Khalifa University, Abu Dhabi, United Arab Emirates
| |
Collapse
|
23
|
da Fonseca-Pereira P, Souza PVL, Hou LY, Schwab S, Geigenberger P, Nunes-Nesi A, Timm S, Fernie AR, Thormählen I, Araújo WL, Daloso DM. Thioredoxin h2 contributes to the redox regulation of mitochondrial photorespiratory metabolism. PLANT, CELL & ENVIRONMENT 2020; 43:188-208. [PMID: 31378951 DOI: 10.1111/pce.13640] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Revised: 07/24/2019] [Accepted: 07/31/2019] [Indexed: 05/18/2023]
Abstract
Thioredoxins (TRXs) are important proteins involved in redox regulation of metabolism. In plants, it has been shown that the mitochondrial metabolism is regulated by the mitochondrial TRX system. However, the functional significance of TRX h2, which is found at both cytosol and mitochondria, remains unclear. Arabidopsis plants lacking TRX h2 showed delayed seed germination and reduced respiration alongside impaired stomatal and mesophyll conductance, without impacting photosynthesis under ambient O2 conditions. However, an increase in the stoichiometry of photorespiratory CO2 release was found during O2 -dependent gas exchange measurements in trxh2 mutants. Metabolite profiling of trxh2 leaves revealed alterations in key metabolites of photorespiration and in several metabolites involved in respiration and amino acid metabolism. Decreased abundance of serine hydroxymethyltransferase and glycine decarboxylase (GDC) H and L subunits as well as reduced NADH/NAD+ ratios were also observed in trxh2 mutants. We further demonstrated that the redox status of GDC-L is altered in trxh2 mutants in vivo and that recombinant TRX h2 can deactivate GDC-L in vitro, indicating that this protein is redox regulated by the TRX system. Collectively, our results demonstrate that TRX h2 plays an important role in the redox regulation of mitochondrial photorespiratory metabolism.
Collapse
Affiliation(s)
- Paula da Fonseca-Pereira
- Max-Planck Partner Group at the Departamento de Biologia Vegetal, Universidade Federal de Viçosa, Viçosa, 36570-900, Brazil
| | - Paulo V L Souza
- Departamento de Bioquímica e Biologia Molecular, Universidade Federal do Ceará, Fortaleza, 60451-970, Brazil
| | - Liang-Yu Hou
- Department Biology I, Ludwig-Maximilians-Universität München, Planegg-Martinsried, 82152, Germany
| | - Saskia Schwab
- Plant Physiology Department, University of Rostock, Rostock, D-18051, Germany
| | - Peter Geigenberger
- Department Biology I, Ludwig-Maximilians-Universität München, Planegg-Martinsried, 82152, Germany
| | - Adriano Nunes-Nesi
- Max-Planck Partner Group at the Departamento de Biologia Vegetal, Universidade Federal de Viçosa, Viçosa, 36570-900, Brazil
| | - Stefan Timm
- Plant Physiology Department, University of Rostock, Rostock, D-18051, Germany
| | - Alisdair R Fernie
- Department Willmitzer, Max-Planck-Institute of Molecular Plant Physiology, Potsdam-Golm, 14476, Germany
| | - Ina Thormählen
- Department Biology I, Ludwig-Maximilians-Universität München, Planegg-Martinsried, 82152, Germany
| | - Wagner L Araújo
- Max-Planck Partner Group at the Departamento de Biologia Vegetal, Universidade Federal de Viçosa, Viçosa, 36570-900, Brazil
| | - Danilo M Daloso
- Departamento de Bioquímica e Biologia Molecular, Universidade Federal do Ceará, Fortaleza, 60451-970, Brazil
| |
Collapse
|
24
|
Geng C, Jung Y, Renaud N, Honavar V, Bonvin AMJJ, Xue LC. iScore: a novel graph kernel-based function for scoring protein-protein docking models. Bioinformatics 2020; 36:112-121. [PMID: 31199455 PMCID: PMC6956772 DOI: 10.1093/bioinformatics/btz496] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Revised: 05/08/2019] [Accepted: 06/11/2019] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Protein complexes play critical roles in many aspects of biological functions. Three-dimensional (3D) structures of protein complexes are critical for gaining insights into structural bases of interactions and their roles in the biomolecular pathways that orchestrate key cellular processes. Because of the expense and effort associated with experimental determinations of 3D protein complex structures, computational docking has evolved as a valuable tool to predict 3D structures of biomolecular complexes. Despite recent progress, reliably distinguishing near-native docking conformations from a large number of candidate conformations, the so-called scoring problem, remains a major challenge. RESULTS Here we present iScore, a novel approach to scoring docked conformations that combines HADDOCK energy terms with a score obtained using a graph representation of the protein-protein interfaces and a measure of evolutionary conservation. It achieves a scoring performance competitive with, or superior to, that of state-of-the-art scoring functions on two independent datasets: (i) Docking software-specific models and (ii) the CAPRI score set generated by a wide variety of docking approaches (i.e. docking software-non-specific). iScore ranks among the top scoring approaches on the CAPRI score set (13 targets) when compared with the 37 scoring groups in CAPRI. The results demonstrate the utility of combining evolutionary, topological and energetic information for scoring docked conformations. This work represents the first successful demonstration of graph kernels to protein interfaces for effective discrimination of near-native and non-native conformations of protein complexes. AVAILABILITY AND IMPLEMENTATION The iScore code is freely available from Github: https://github.com/DeepRank/iScore (DOI: 10.5281/zenodo.2630567). And the docking models used are available from SBGrid: https://data.sbgrid.org/dataset/684). SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Cunliang Geng
- Bijvoet Center for Biomolecular Research, Faculty of Science – Chemistry, Utrecht University, Utrecht 3584 CH, The Netherlands
| | - Yong Jung
- Bioinformatics & Genomics Graduate Program, Pennsylvania State University, University Park, PA 16802, USA
- Artificial Intelligence Research Laboratory, Pennsylvania State University, University Park, PA 16823, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA
| | - Nicolas Renaud
- Netherlands eScience Center, Amsterdam 1098 XG, The Netherlands
| | - Vasant Honavar
- Bioinformatics & Genomics Graduate Program, Pennsylvania State University, University Park, PA 16802, USA
- Artificial Intelligence Research Laboratory, Pennsylvania State University, University Park, PA 16823, USA
- Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA
- Center for Big Data Analytics and Discovery Informatics, Pennsylvania State University, University Park, PA 16823, USA
- Institute for Cyberscience, University Park, PA 16802, USA
- Clinical and Translational Sciences Institute, University Park, PA 16802, USA
- College of Information Sciences & Technology, Pennsylvania State University, University Park, PA 16802, USA
| | - Alexandre M J J Bonvin
- Bijvoet Center for Biomolecular Research, Faculty of Science – Chemistry, Utrecht University, Utrecht 3584 CH, The Netherlands
| | - Li C Xue
- Bijvoet Center for Biomolecular Research, Faculty of Science – Chemistry, Utrecht University, Utrecht 3584 CH, The Netherlands
| |
Collapse
|
25
|
Guo F, Zou Q, Yang G, Wang D, Tang J, Xu J. Identifying protein-protein interface via a novel multi-scale local sequence and structural representation. BMC Bioinformatics 2019; 20:483. [PMID: 31874604 PMCID: PMC6929278 DOI: 10.1186/s12859-019-3048-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2019] [Accepted: 08/21/2019] [Indexed: 12/23/2022] Open
Abstract
Background Protein-protein interaction plays a key role in a multitude of biological processes, such as signal transduction, de novo drug design, immune responses, and enzymatic activities. Gaining insights of various binding abilities can deepen our understanding of the interaction. It is of great interest to understand how proteins in a complex interact with each other. Many efficient methods have been developed for identifying protein-protein interface. Results In this paper, we obtain the local information on protein-protein interface, through multi-scale local average block and hexagon structure construction. Given a pair of proteins, we use a trained support vector regression (SVR) model to select best configurations. On Benchmark v4.0, our method achieves average Irmsd value of 3.28Å and overall Fnat value of 63%, which improves upon Irmsd of 3.89Å and Fnat of 49% for ZRANK, and Irmsd of 3.99Å and Fnat of 46% for ClusPro. On CAPRI targets, our method achieves average Irmsd value of 3.45Å and overall Fnat value of 46%, which improves upon Irmsd of 4.18Å and Fnat of 40% for ZRANK, and Irmsd of 5.12Å and Fnat of 32% for ClusPro. The success rates by our method, FRODOCK 2.0, InterEvDock and SnapDock on Benchmark v4.0 are 41.5%, 29.0%, 29.4% and 37.0%, respectively. Conclusion Experiments show that our method performs better than some state-of-the-art methods, based on the prediction quality improved in terms of CAPRI evaluation criteria. All these results demonstrate that our method is a valuable technological tool for identifying protein-protein interface.
Collapse
Affiliation(s)
- Fei Guo
- College of Intelligence and Computing, Tianjin University, Tianjin, People's Republic of China.
| | - Quan Zou
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu, People's Republic of China
| | - Guang Yang
- School of Economics, Nankai University, Tianjin, People's Republic of China
| | - Dan Wang
- Department of Computer Science, City University of Hong Kong, Kowloon Tong, Hong Kong
| | - Jijun Tang
- College of Intelligence and Computing, Tianjin University, Tianjin, People's Republic of China.,Department of Computer Science and Engineering, University of South Carolina, Columbia, USA
| | - Junhai Xu
- College of Intelligence and Computing, Tianjin University, Tianjin, People's Republic of China
| |
Collapse
|
26
|
Siebenmorgen T, Zacharias M. Computational prediction of protein–protein binding affinities. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1448] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Till Siebenmorgen
- Physics Department T38 Technical University of Munich Garching Germany
| | - Martin Zacharias
- Physics Department T38 Technical University of Munich Garching Germany
| |
Collapse
|
27
|
DRAMP 2.0, an updated data repository of antimicrobial peptides. Sci Data 2019; 6:148. [PMID: 31409791 PMCID: PMC6692298 DOI: 10.1038/s41597-019-0154-y] [Citation(s) in RCA: 173] [Impact Index Per Article: 34.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Accepted: 07/17/2019] [Indexed: 12/20/2022] Open
Abstract
Data Repository of Antimicrobial Peptides (DRAMP, http://dramp.cpu-bioinfor.org/) is an open-access comprehensive database containing general, patent and clinical antimicrobial peptides (AMPs). Currently DRAMP has been updated to version 2.0, it contains a total of 19,899 entries (newly added 2,550 entries), including 5,084 general entries, 14,739 patent entries, and 76 clinical entries. The update covers new entries, structures, annotations, classifications and downloads. Compared with APD and CAMP, DRAMP contains 14,040 (70.56% in DRAMP) non-overlapping sequences. In order to facilitate users to trace original references, PubMed_ID of references have been contained in activity information. The data of DRAMP can be downloaded by dataset and activity, and the website source code is also available on dedicatedly designed download webpage. Although thousands of AMPs have been reported, only a few parts have entered clinical stage. In the paper, we described several AMPs in clinical trials, including their properties, indications and clinicaltrials.gov identifiers. Finally, we provide the applications of DRAMP in the development of AMPs.
Collapse
|
28
|
Developments in integrative modeling with dynamical interfaces. Curr Opin Struct Biol 2019; 56:11-17. [DOI: 10.1016/j.sbi.2018.10.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2018] [Revised: 10/26/2018] [Accepted: 10/27/2018] [Indexed: 11/19/2022]
|
29
|
Geng C, Xue LC, Roel‐Touris J, Bonvin AMJJ. Finding the ΔΔ
G
spot: Are predictors of binding affinity changes upon mutations in protein–protein interactions ready for it? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1410] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Cunliang Geng
- Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry Utrecht University Utrecht The Netherlands
| | - Li C. Xue
- Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry Utrecht University Utrecht The Netherlands
| | - Jorge Roel‐Touris
- Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry Utrecht University Utrecht The Netherlands
| | - Alexandre M. J. J. Bonvin
- Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry Utrecht University Utrecht The Netherlands
| |
Collapse
|
30
|
Hogues H, Gaudreault F, Corbeil CR, Deprez C, Sulea T, Purisima EO. ProPOSE: Direct Exhaustive Protein-Protein Docking with Side Chain Flexibility. J Chem Theory Comput 2018; 14:4938-4947. [PMID: 30107730 DOI: 10.1021/acs.jctc.8b00225] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Despite decades of development, protein-protein docking remains a largely unsolved problem. The main difficulties are the immense space spanned by the translational and rotational degrees of freedom and the prediction of the conformational changes of proteins upon binding. FFT is generally the preferred method to exhaustively explore the translation-rotation space at a fine grid resolution, albeit with the trade-off of approximating force fields with correlation functions. This work presents a direct search alternative that samples the states in Cartesian space at the same resolution and computational cost as standard FFT methods. Operating in real space allows the use of standard force field functional forms used in typical non-FFT methods as well as the implementation of strategies for focused exploration of conformational flexibility. Currently, a few misplaced side chains can cause docking programs to fail. This work specifically addresses the problem of side chain rearrangements upon complex formation. Based on the observation that most side chains retain their unbound conformation upon binding, each rigidly docked pose is initially scored ignoring up to a limited number of side chain overlaps which are resolved in subsequent repacking and minimization steps. On test systems where side chains are altered and backbones held in their bound state, this implementation provides significantly better native pose recovery and higher quality (lower RMSD) predictions when compared with five of the most popular docking programs. The method is implemented in the software program ProPOSE (Protein Pose Optimization by Systematic Enumeration).
Collapse
Affiliation(s)
- Hervé Hogues
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| | - Francis Gaudreault
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| | - Christopher R Corbeil
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| | - Christophe Deprez
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| | - Traian Sulea
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| | - Enrico O Purisima
- Human Health Therapeutics , National Research Council Canada , 6100 Royalmount Avenue , Montreal , Quebec H4P 2R2 , Canada
| |
Collapse
|
31
|
Guidolin D, Marcoli M, Tortorella C, Maura G, Agnati LF. G protein-coupled receptor-receptor interactions give integrative dynamics to intercellular communication. Rev Neurosci 2018; 29:703-726. [DOI: 10.1515/revneuro-2017-0087] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2017] [Accepted: 01/01/2018] [Indexed: 01/14/2023]
Abstract
Abstract
The proposal of receptor-receptor interactions (RRIs) in the early 1980s broadened the view on the role of G protein-coupled receptors (GPCR) in the dynamics of the intercellular communication. RRIs, indeed, allow GPCR to operate not only as monomers but also as receptor complexes, in which the integration of the incoming signals depends on the number, spatial arrangement, and order of activation of the protomers forming the complex. The main biochemical mechanisms controlling the functional interplay of GPCR in the receptor complexes are direct allosteric interactions between protomer domains. The formation of these macromolecular assemblies has several physiologic implications in terms of the modulation of the signaling pathways and interaction with other membrane proteins. It also impacts on the emerging field of connectomics, as it contributes to set and tune the synaptic strength. Furthermore, recent evidence suggests that the transfer of GPCR and GPCR complexes between cells via the exosome pathway could enable the target cells to recognize/decode transmitters and/or modulators for which they did not express the pertinent receptors. Thus, this process may also open the possibility of a new type of redeployment of neural circuits. The fundamental aspects of GPCR complex formation and function are the focus of the present review article.
Collapse
Affiliation(s)
- Diego Guidolin
- Department of Neuroscience , University of Padova, via Gabelli 65 , I-35121 Padova , Italy
| | - Manuela Marcoli
- Department of Pharmacy and Center of Excellence for Biomedical Research , University of Genova , I-16126 Genova , Italy
| | - Cinzia Tortorella
- Department of Neuroscience , University of Padova, via Gabelli 65 , I-35121 Padova , Italy
| | - Guido Maura
- Department of Pharmacy and Center of Excellence for Biomedical Research , University of Genova , I-16126 Genova , Italy
| | - Luigi F. Agnati
- Department of Biomedical Sciences , University of Modena and Reggio Emilia , I-41121 Modena , Italy
- Department of Neuroscience , Karolinska Institutet , S-17177 Stockholm , Sweden
| |
Collapse
|
32
|
Fraternali F. [Protein-protein interacting networks, their structures and disease-related mutations]. Biol Aujourdhui 2018; 211:223-228. [PMID: 29412132 DOI: 10.1051/jbio/2017031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2017] [Indexed: 11/14/2022]
Abstract
In recent years, the comparison of protein interactomes has identified conserved modules, that could represent functional nuclei with a common ancestry. Within this context, recent analyses of protein-protein interacting networks have led to a debate on the influence of the experimental method on the quality and biological pertinence of these data. It is crucial to understand the measure in which divergence between networks of different species reflect sampling biases in respective experimental methods, as opposed to topological features dictated by biological functionality. This aspect requires novel, precise and practical mathematical tools, to quantify and compare high resolution networks. To this end, we have studied the relationship between pools of aleatory graphs and real biological signalization networks, while stressing the number of graph cycles in the networks, which represent complexes in experimental protein interactomes. By combining methods for graph and algorithm dynamics to count the loops, we evaluate the relative importance of the loops in biological networks in comparison with network analyses.
Collapse
Affiliation(s)
- Franca Fraternali
- Randall Division of Cellular and Molecular Biology, King's College, London, UK
| |
Collapse
|
33
|
Abstract
ExoU is a type III-secreted cytotoxin expressing A2 phospholipase activity when injected into eukaryotic target cells by the bacterium Pseudomonas aeruginosa The enzymatic activity of ExoU is undetectable in vitro unless ubiquitin, a required cofactor, is added to the reaction. The role of ubiquitin in facilitating ExoU enzymatic activity is poorly understood but of significance for designing inhibitors to prevent tissue injury during infections with strains of P. aeruginosa producing this toxin. Most ubiquitin-binding proteins, including ExoU, demonstrate a low (micromolar) affinity for monoubiquitin (monoUb). Additionally, ExoU is a large and dynamic protein, limiting the applicability of traditional structural techniques such as NMR and X-ray crystallography to define this protein-protein interaction. Recent advancements in computational methods, however, have allowed high-resolution protein modeling using sparse data. In this study, we combine double electron-electron resonance (DEER) spectroscopy and Rosetta modeling to identify potential binding interfaces of ExoU and monoUb. The lowest-energy scoring model was tested using biochemical, biophysical, and biological techniques. To verify the binding interface, Rosetta was used to design a panel of mutations to modulate binding, including one variant with enhanced binding affinity. Our analyses show the utility of computational modeling when combined with sensitive biological assays and biophysical approaches that are exquisitely suited for large dynamic proteins.
Collapse
|
34
|
Ivanov SM, Cawley A, Huber RG, Bond PJ, Warwicker J. Protein-protein interactions in paralogues: Electrostatics modulates specificity on a conserved steric scaffold. PLoS One 2017; 12:e0185928. [PMID: 29016650 PMCID: PMC5634604 DOI: 10.1371/journal.pone.0185928] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2017] [Accepted: 09/21/2017] [Indexed: 12/05/2022] Open
Abstract
An improved knowledge of protein-protein interactions is essential for better understanding of metabolic and signaling networks, and cellular function. Progress tends to be based on structure determination and predictions using known structures, along with computational methods based on evolutionary information or detailed atomistic descriptions. We hypothesized that for the case of interactions across a common interface, between proteins from a pair of paralogue families or within a family of paralogues, a relatively simple interface description could distinguish between binding and non-binding pairs. Using binding data for several systems, and large-scale comparative modeling based on known template complex structures, it is found that charge-charge interactions (for groups bearing net charge) are generally a better discriminant than buried non-polar surface. This is particularly the case for paralogue families that are less divergent, with more reliable comparative modeling. We suggest that electrostatic interactions are major determinants of specificity in such systems, an observation that could be used to predict binding partners.
Collapse
Affiliation(s)
- Stefan M. Ivanov
- Manchester Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester, United Kingdom
| | - Andrew Cawley
- Manchester Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester, United Kingdom
| | - Roland G. Huber
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), Matrix, Singapore, Singapore
| | - Peter J. Bond
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), Matrix, Singapore, Singapore
- Department of Biological Sciences, National University of Singapore, Singapore, Singapore
| | - Jim Warwicker
- Manchester Institute of Biotechnology, School of Chemistry, The University of Manchester, Manchester, United Kingdom
- * E-mail:
| |
Collapse
|