1
|
Kliche J, Simonetti L, Krystkowiak I, Kuss H, Diallo M, Rask E, Nilsson J, Davey NE, Ivarsson Y. Proteome-scale characterisation of motif-based interactome rewiring by disease mutations. Mol Syst Biol 2024:10.1038/s44320-024-00055-4. [PMID: 39009827 DOI: 10.1038/s44320-024-00055-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 06/14/2024] [Accepted: 06/28/2024] [Indexed: 07/17/2024] Open
Abstract
Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.
Collapse
Affiliation(s)
- Johanna Kliche
- Department of Chemistry - BMC, Box 576, Husargatan 3, 751 23, Uppsala, Sweden
| | - Leandro Simonetti
- Department of Chemistry - BMC, Box 576, Husargatan 3, 751 23, Uppsala, Sweden
| | - Izabella Krystkowiak
- Division of Cancer Biology, Institute of Cancer Research, Chester Beatty Laboratories, 237 Fulham Road, SW3 6JB, Chelsea, London, UK
| | - Hanna Kuss
- Department of Chemistry - BMC, Box 576, Husargatan 3, 751 23, Uppsala, Sweden
- University of Münster, Institute of Pharmaceutical and Medicinal Chemistry, DE-48149, Münster, Germany
| | - Marcel Diallo
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3B, 2200, Copenhagen, Denmark
| | - Emma Rask
- Department of Chemistry - BMC, Box 576, Husargatan 3, 751 23, Uppsala, Sweden
| | - Jakob Nilsson
- Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3B, 2200, Copenhagen, Denmark
| | - Norman E Davey
- Division of Cancer Biology, Institute of Cancer Research, Chester Beatty Laboratories, 237 Fulham Road, SW3 6JB, Chelsea, London, UK.
| | - Ylva Ivarsson
- Department of Chemistry - BMC, Box 576, Husargatan 3, 751 23, Uppsala, Sweden.
| |
Collapse
|
2
|
González-Avendaño M, López J, Vergara-Jaque A, Cerda O. The power of computational proteomics platforms to decipher protein-protein interactions. Curr Opin Struct Biol 2024; 88:102882. [PMID: 39003917 DOI: 10.1016/j.sbi.2024.102882] [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: 03/28/2024] [Revised: 05/31/2024] [Accepted: 06/19/2024] [Indexed: 07/16/2024]
Abstract
Adopting computational tools for analyzing extensive biological datasets has profoundly transformed our understanding and interpretation of biological phenomena. Innovative platforms have emerged, providing automated analysis to unravel essential insights about proteins and the complexities of their interactions. These computational advancements align with traditional studies, which employ experimental techniques to discern and quantify physical and functional protein-protein interactions (PPIs). Among these techniques, tandem mass spectrometry is notably recognized for its precision and sensitivity in identifying PPIs. These approaches might serve as important information enabling the identification of PPIs with potential pharmacological significance. This review aims to convey our experience using computational tools for detecting PPI networks and offer an analysis of platforms that facilitate predictions derived from experimental data.
Collapse
Affiliation(s)
- Mariela González-Avendaño
- Center for Bioinformatics, Simulation and Modeling (CBSM), Faculty of Engineering, Universidad de Talca, Talca, Chile; Millennium Nucleus of Ion Channel-Associated Diseases (MiNICAD), Santiago, Chile
| | - Joaquín López
- Millennium Nucleus of Ion Channel-Associated Diseases (MiNICAD), Santiago, Chile; Program of Cellular and Molecular Biology, Institute of Biomedical Sciences (ICBM), Faculty of Medicine, Universidad de Chile, Santiago, Chile
| | - Ariela Vergara-Jaque
- Center for Bioinformatics, Simulation and Modeling (CBSM), Faculty of Engineering, Universidad de Talca, Talca, Chile; Millennium Nucleus of Ion Channel-Associated Diseases (MiNICAD), Santiago, Chile.
| | - Oscar Cerda
- Millennium Nucleus of Ion Channel-Associated Diseases (MiNICAD), Santiago, Chile; Program of Cellular and Molecular Biology, Institute of Biomedical Sciences (ICBM), Faculty of Medicine, Universidad de Chile, Santiago, Chile.
| |
Collapse
|
3
|
Zhou Y, Myung Y, Rodrigues CM, Ascher D. DDMut-PPI: predicting effects of mutations on protein-protein interactions using graph-based deep learning. Nucleic Acids Res 2024; 52:W207-W214. [PMID: 38783112 PMCID: PMC11223791 DOI: 10.1093/nar/gkae412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 04/30/2024] [Accepted: 05/02/2024] [Indexed: 05/25/2024] Open
Abstract
Protein-protein interactions (PPIs) play a vital role in cellular functions and are essential for therapeutic development and understanding diseases. However, current predictive tools often struggle to balance efficiency and precision in predicting the effects of mutations on these complex interactions. To address this, we present DDMut-PPI, a deep learning model that efficiently and accurately predicts changes in PPI binding free energy upon single and multiple point mutations. Building on the robust Siamese network architecture with graph-based signatures from our prior work, DDMut, the DDMut-PPI model was enhanced with a graph convolutional network operated on the protein interaction interface. We used residue-specific embeddings from ProtT5 protein language model as node features, and a variety of molecular interactions as edge features. By integrating evolutionary context with spatial information, this framework enables DDMut-PPI to achieve a robust Pearson correlation of up to 0.75 (root mean squared error: 1.33 kcal/mol) in our evaluations, outperforming most existing methods. Importantly, the model demonstrated consistent performance across mutations that increase or decrease binding affinity. DDMut-PPI offers a significant advancement in the field and will serve as a valuable tool for researchers probing the complexities of protein interactions. DDMut-PPI is freely available as a web server and an application programming interface at https://biosig.lab.uq.edu.au/ddmut_ppi.
Collapse
Affiliation(s)
- Yunzhuo Zhou
- The Australian Centre for Ecogenomics, School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, Queensland 4072, Australia
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
| | - YooChan Myung
- The Australian Centre for Ecogenomics, School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, Queensland 4072, Australia
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
| | - Carlos H M Rodrigues
- The Australian Centre for Ecogenomics, School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, Queensland 4072, Australia
| | - David B Ascher
- The Australian Centre for Ecogenomics, School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, Queensland 4072, Australia
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
| |
Collapse
|
4
|
Passi G, Lieberman S, Zahdeh F, Murik O, Renbaum P, Beeri R, Linial M, May D, Levy-Lahad E, Schneidman-Duhovny D. Discovering predisposing genes for hereditary breast cancer using deep learning. Brief Bioinform 2024; 25:bbae346. [PMID: 39038933 PMCID: PMC11262808 DOI: 10.1093/bib/bbae346] [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: 01/16/2024] [Revised: 04/18/2024] [Accepted: 07/04/2024] [Indexed: 07/24/2024] Open
Abstract
Breast cancer (BC) is the most common malignancy affecting Western women today. It is estimated that as many as 10% of BC cases can be attributed to germline variants. However, the genetic basis of the majority of familial BC cases has yet to be identified. Discovering predisposing genes contributing to familial BC is challenging due to their presumed rarity, low penetrance, and complex biological mechanisms. Here, we focused on an analysis of rare missense variants in a cohort of 12 families of Middle Eastern origins characterized by a high incidence of BC cases. We devised a novel, high-throughput, variant analysis pipeline adapted for family studies, which aims to analyze variants at the protein level by employing state-of-the-art machine learning models and three-dimensional protein structural analysis. Using our pipeline, we analyzed 1218 rare missense variants that are shared between affected family members and classified 80 genes as candidate pathogenic. Among these genes, we found significant functional enrichment in peroxisomal and mitochondrial biological pathways which segregated across seven families in the study and covered diverse ethnic groups. We present multiple evidence that peroxisomal and mitochondrial pathways play an important, yet underappreciated, role in both germline BC predisposition and BC survival.
Collapse
Affiliation(s)
- Gal Passi
- The Rachel and Selim Benin School of Computer Science and Engineering, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Sari Lieberman
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
- Faculty of Medicine, The Hebrew University of Jerusalem, Ein Kerem PO Box 12271 Jerusalem 9112102, Israel
| | - Fouad Zahdeh
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
| | - Omer Murik
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
| | - Paul Renbaum
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
| | - Rachel Beeri
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
| | - Michal Linial
- Department of Biological Chemistry, Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Dalit May
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
- Clalit Health Services, Jerusalem, Israel
| | - Ephrat Levy-Lahad
- The Fuld Family Medical Genetics Institute, Shaare Zedek Medical Center 12 Bayit St., Jerusalem 9103101, Israel
- The Eisenberg R&D Authority, Shaare Zedek Medical Center, 12 Bayit St., Jerusalem 9103101, Israel
- Faculty of Medicine, The Hebrew University of Jerusalem, Ein Kerem PO Box 12271 Jerusalem 9112102, Israel
| | - Dina Schneidman-Duhovny
- The Rachel and Selim Benin School of Computer Science and Engineering, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| |
Collapse
|
5
|
Gurusinghe SNS, Shifman JM. Cold Spot SCANNER: Colab Notebook for predicting cold spots in protein-protein interfaces. BMC Bioinformatics 2024; 25:172. [PMID: 38689238 PMCID: PMC11061940 DOI: 10.1186/s12859-024-05796-5] [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: 05/31/2023] [Accepted: 04/22/2024] [Indexed: 05/02/2024] Open
Abstract
BACKGROUND Protein-protein interactions (PPIs) are conveyed through binding interfaces or surface patches on proteins that become buried upon binding. Structural and biophysical analysis of many protein-protein interfaces revealed certain unique features of these surfaces that determine the energetics of interactions and play a critical role in protein evolution. One of the significant aspects of binding interfaces is the presence of binding hot spots, where mutations are highly deleterious for binding. Conversely, binding cold spots are positions occupied by suboptimal amino acids and several mutations in such positions could lead to affinity enhancement. While there are many software programs for identification of hot spot positions, there is currently a lack of software for cold spot detection. RESULTS In this paper, we present Cold Spot SCANNER, a Colab Notebook, which scans a PPI binding interface and identifies cold spots resulting from cavities, unfavorable charge-charge, and unfavorable charge-hydrophobic interactions. The software offers a Py3DMOL-based interface that allows users to visualize cold spots in the context of the protein structure and generates a zip file containing the results for easy download. CONCLUSIONS Cold spot identification is of great importance to protein engineering studies and provides a useful insight into protein evolution. Cold Spot SCANNER is open to all users without login requirements and can be accessible at: https://colab. RESEARCH google.com/github/sagagugit/Cold-Spot-Scanner/blob/main/Cold_Spot_Scanner.ipynb .
Collapse
Affiliation(s)
- Sagara N S Gurusinghe
- Department of Biological Chemistry, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Julia M Shifman
- Department of Biological Chemistry, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel.
| |
Collapse
|
6
|
Schweke H, Pacesa M, Levin T, Goverde CA, Kumar P, Duhoo Y, Dornfeld LJ, Dubreuil B, Georgeon S, Ovchinnikov S, Woolfson DN, Correia BE, Dey S, Levy ED. An atlas of protein homo-oligomerization across domains of life. Cell 2024; 187:999-1010.e15. [PMID: 38325366 DOI: 10.1016/j.cell.2024.01.022] [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] [Received: 06/08/2023] [Revised: 11/03/2023] [Accepted: 01/15/2024] [Indexed: 02/09/2024]
Abstract
Protein structures are essential to understanding cellular processes in molecular detail. While advances in artificial intelligence revealed the tertiary structure of proteins at scale, their quaternary structure remains mostly unknown. We devise a scalable strategy based on AlphaFold2 to predict homo-oligomeric assemblies across four proteomes spanning the tree of life. Our results suggest that approximately 45% of an archaeal proteome and a bacterial proteome and 20% of two eukaryotic proteomes form homomers. Our predictions accurately capture protein homo-oligomerization, recapitulate megadalton complexes, and unveil hundreds of homo-oligomer types, including three confirmed experimentally by structure determination. Integrating these datasets with omics information suggests that a majority of known protein complexes are symmetric. Finally, these datasets provide a structural context for interpreting disease mutations and reveal coiled-coil regions as major enablers of quaternary structure evolution in human. Our strategy is applicable to any organism and provides a comprehensive view of homo-oligomerization in proteomes.
Collapse
Affiliation(s)
- Hugo Schweke
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Martin Pacesa
- Laboratory of Protein Design and Immunoengineering, École Polytechnique Fédérale de Lausanne and Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Tal Levin
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Casper A Goverde
- Laboratory of Protein Design and Immunoengineering, École Polytechnique Fédérale de Lausanne and Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Prasun Kumar
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK; School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK; Bristol BioDesign Institute, University of Bristol, Life Sciences Building, Bristol BS8 1TQ, UK; Max Planck-Bristol Centre for Minimal Biology, University of Bristol, Cantock's Close, Bristol BS8 1TS, UK
| | - Yoan Duhoo
- Protein Production and Structure Characterization Core Facility (PTPSP), School of Life Sciences, École polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Lars J Dornfeld
- Laboratory of Protein Design and Immunoengineering, École Polytechnique Fédérale de Lausanne and Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Benjamin Dubreuil
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Sandrine Georgeon
- Laboratory of Protein Design and Immunoengineering, École Polytechnique Fédérale de Lausanne and Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Sergey Ovchinnikov
- John Harvard Distinguished Science Fellowship Program, Harvard University, Cambridge, MA, USA
| | - Derek N Woolfson
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK; School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK; Bristol BioDesign Institute, University of Bristol, Life Sciences Building, Bristol BS8 1TQ, UK; Max Planck-Bristol Centre for Minimal Biology, University of Bristol, Cantock's Close, Bristol BS8 1TS, UK.
| | - Bruno E Correia
- Laboratory of Protein Design and Immunoengineering, École Polytechnique Fédérale de Lausanne and Swiss Institute of Bioinformatics, Lausanne, Switzerland.
| | - Sucharita Dey
- Department of Bioscience and Bioengineering, Indian Institute of Technology Jodhpur, Rajasthan, India.
| | - Emmanuel D Levy
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Rehovot, Israel.
| |
Collapse
|
7
|
Rana MM, Nguyen DD. Geometric Graph Learning to Predict Changes in Binding Free Energy and Protein Thermodynamic Stability upon Mutation. J Phys Chem Lett 2023; 14:10870-10879. [PMID: 38032742 DOI: 10.1021/acs.jpclett.3c02679] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
Accurate prediction of binding free energy changes upon mutations is vital for optimizing drugs, designing proteins, understanding genetic diseases, and cost-effective virtual screening. While machine learning methods show promise in this domain, achieving accuracy and generalization across diverse data sets remains a challenge. This study introduces Geometric Graph Learning for Protein-Protein Interactions (GGL-PPI), a novel approach integrating geometric graph representation and machine learning to forecast mutation-induced binding free energy changes. GGL-PPI leverages atom-level graph coloring and multiscale weighted colored geometric subgraphs to capture structural features of biomolecules, demonstrating superior performance on three standard data sets, namely, AB-Bind, SKEMPI 1.0, and SKEMPI 2.0 data sets. The model's efficacy extends to predicting protein thermodynamic stability in a blind test set, providing unbiased predictions for both direct and reverse mutations and showcasing notable generalization. GGL-PPI's precision in predicting changes in binding free energy and stability due to mutations enhances our comprehension of protein complexes, offering valuable insights for drug design endeavors.
Collapse
Affiliation(s)
- Md Masud Rana
- Department of Mathematics, University of Kentucky, Lexington, Kentucky 40506, United States
| | - Duc Duy Nguyen
- Department of Mathematics, University of Kentucky, Lexington, Kentucky 40506, United States
| |
Collapse
|
8
|
Sun H, Wang J, Wu H, Lin S, Chen J, Wei J, Lv S, Xiong Y, Wei DQ. A Multimodal Deep Learning Framework for Predicting PPI-Modulator Interactions. J Chem Inf Model 2023; 63:7363-7372. [PMID: 38037990 DOI: 10.1021/acs.jcim.3c01527] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2023]
Abstract
Protein-protein interactions (PPIs) are essential for various biological processes and diseases. However, most existing computational methods for identifying PPI modulators require either target structure or reference modulators, which restricts their applicability to novel PPI targets. To address this challenge, we propose MultiPPIMI, a sequence-based deep learning framework that predicts the interaction between any given PPI target and modulator. MultiPPIMI integrates multimodal representations of PPI targets and modulators and uses a bilinear attention network to capture intermolecular interactions. Experimental results on our curated benchmark data set show that MultiPPIMI achieves an average AUROC of 0.837 in three cold-start scenarios and an AUROC of 0.994 in the random-split scenario. Furthermore, the case study shows that MultiPPIMI can assist molecular docking simulations in screening inhibitors of Keap1/Nrf2 PPI interactions. We believe that the proposed method provides a promising way to screen PPI-targeted modulators.
Collapse
Affiliation(s)
- Heqi Sun
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jianmin Wang
- The Interdisciplinary Graduate Program in Integrative Biotechnology and Translational Medicine, Yonsei University, Incheon 21983, Republic of Korea
| | - Hongyan Wu
- Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
| | - Shenggeng Lin
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Junwei Chen
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Jinghua Wei
- Department of Chemistry, University of Toronto, Toronto M5R 0A3, Canada
| | - Shuai Lv
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Yi Xiong
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
- Shanghai Artificial Intelligence Laboratory, Shanghai 200232, China
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
- Peng Cheng National Laboratory, Shenzhen 518055, China
- Zhongjing Research and Industrialization Institute of Chinese Medicine, Nanyang 473006, China
| |
Collapse
|
9
|
Wen B, Hu S, Yin J, Wu J, Guo W. Molecular Evolution and Protein Structure Variation of Dkk Family. Genes (Basel) 2023; 14:1863. [PMID: 37895211 PMCID: PMC10606412 DOI: 10.3390/genes14101863] [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] [Received: 08/24/2023] [Revised: 09/17/2023] [Accepted: 09/22/2023] [Indexed: 10/29/2023] Open
Abstract
Dkks have inhibitory effects on the Wnt signaling pathway, which is involved in the development of skin and its appendages and the regulation of hair growth. The nucleotide sequences were compared and analyzed to further investigate the relationship between the structure and function of the Dkk gene family and vertebrate epidermal hair. The analysis of the molecular evolution of the Dkk family revealed that the evolution rate of the genes changed significantly after speciation, with the Aves and Reptilia branches showing accelerated evolution. Additionally, positive selection was observed at specific sites. The tertiary structure of the protein was also predicted. The analysis of the functional divergence of the Dkk family revealed that the functional divergence coefficient of each gene was greater than 0, with most of the functional divergence sites were located in the Cys-2 domain and a few in the Cys-1 domain. This suggests that the amino acid and functional divergence sites may play a role in regulating the binding of the Dkk family to LRP5/6, and thus affect the inhibition of Wnt signaling, leading to different functions of Dkk1, Dkk2, and Dkk4 in the development of skin hair follicles. In addition, the Dkk families of Aves and Reptilia may have undergone adaptive evolution and functional divergence.
Collapse
Affiliation(s)
- Binhong Wen
- College of Animal Science and Technology, Inner Mongolia Minzu University, Tongliao 028000, China;
| | - Sile Hu
- College of Life Science, Inner Mongolia Minzu University, Tongliao 028000, China;
| | - Jun Yin
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot 010018, China;
| | - Jianghong Wu
- College of Animal Science and Technology, Inner Mongolia Minzu University, Tongliao 028000, China;
| | - Wenrui Guo
- College of Veterinary Medicine, Inner Mongolia Agricultural University, Hohhot 010018, China
| |
Collapse
|
10
|
Kaur S, Vashistt J, Changotra H. Identification of molecular signatures and molecular dynamics simulation of highly deleterious missense variants of key autophagy regulator beclin 1: a computational based approach. J Biomol Struct Dyn 2023:1-14. [PMID: 37640005 DOI: 10.1080/07391102.2023.2252097] [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: 05/09/2023] [Accepted: 08/21/2023] [Indexed: 08/31/2023]
Abstract
Beclin 1 is a key autophagy regulator that also plays significant roles in other intracellular processes such as vacuolar protein sorting. Beclin 1 protein functions as a scaffold in the formation of a multiprotein assemblage during autophagy. Beclin 1 is involved in various diseases such as cancers, neurodegenerative and autophagy-related disorders. In this study, we have used various in silico tools to scan beclin 1 at the molecular level to find its molecular signatures. We have predicted and analysed deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) of beclin 1 causing alterations in its structure and also affecting its interactions with other proteins. In total, twelve coding region deleterious variants were predicted using sequence-based tools and nine were predicted using various structure-based tools. The molecular dynamics (MD) simulations revealed an altered stability of the native structure due to the introduction of mutations. Destabilization of beclin 1 ECD domain was observed due to nsSNPs W300R and E302K. Beclin 1 deleterious nsSNPs were predicted to show significant effects on beclin 1 interactions with ATG14L1, UVRAG and VPS34 proteins and were also predicted to alter the protein-protein interface of beclin 1 complexes. Additionally, beclin 1 was predicted to have thirty-one potential phosphorylation and three ubiquitination sites. In conclusion, the molecular details of beclin 1 could help in the better understanding of its functioning. The study of nsSNPs and their effect on beclin 1 and its interactions might aid in understanding the basis of anomalies caused due to beclin 1 dysfunction.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sargeet Kaur
- Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India
| | - Jitendraa Vashistt
- Department of Biotechnology and Bioinformatics, Jaypee University of Information Technology, Waknaghat, Solan, India
| | - Harish Changotra
- Department of Molecular Biology and Biochemistry, Guru Nanak Dev University, Amritsar, India
| |
Collapse
|
11
|
Mohseni Behbahani Y, Laine E, Carbone A. Deep Local Analysis deconstructs protein-protein interfaces and accurately estimates binding affinity changes upon mutation. Bioinformatics 2023; 39:i544-i552. [PMID: 37387162 DOI: 10.1093/bioinformatics/btad231] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023] Open
Abstract
MOTIVATION The spectacular recent advances in protein and protein complex structure prediction hold promise for reconstructing interactomes at large-scale and residue resolution. Beyond determining the 3D arrangement of interacting partners, modeling approaches should be able to unravel the impact of sequence variations on the strength of the association. RESULTS In this work, we report on Deep Local Analysis, a novel and efficient deep learning framework that relies on a strikingly simple deconstruction of protein interfaces into small locally oriented residue-centered cubes and on 3D convolutions recognizing patterns within cubes. Merely based on the two cubes associated with the wild-type and the mutant residues, DLA accurately estimates the binding affinity change for the associated complexes. It achieves a Pearson correlation coefficient of 0.735 on about 400 mutations on unseen complexes. Its generalization capability on blind datasets of complexes is higher than the state-of-the-art methods. We show that taking into account the evolutionary constraints on residues contributes to predictions. We also discuss the influence of conformational variability on performance. Beyond the predictive power on the effects of mutations, DLA is a general framework for transferring the knowledge gained from the available non-redundant set of complex protein structures to various tasks. For instance, given a single partially masked cube, it recovers the identity and physicochemical class of the central residue. Given an ensemble of cubes representing an interface, it predicts the function of the complex. AVAILABILITY AND IMPLEMENTATION Source code and models are available at http://gitlab.lcqb.upmc.fr/DLA/DLA.git.
Collapse
Affiliation(s)
- Yasser Mohseni Behbahani
- Laboratory of Computational and Quantitative Biology (LCQB), UMR 7238, Sorbonne Université, CNRS, IBPS, Paris 75005, France
| | - Elodie Laine
- Laboratory of Computational and Quantitative Biology (LCQB), UMR 7238, Sorbonne Université, CNRS, IBPS, Paris 75005, France
| | - Alessandra Carbone
- Laboratory of Computational and Quantitative Biology (LCQB), UMR 7238, Sorbonne Université, CNRS, IBPS, Paris 75005, France
| |
Collapse
|
12
|
David A, Sternberg MJE. Protein structure-based evaluation of missense variants: Resources, challenges and future directions. Curr Opin Struct Biol 2023; 80:102600. [PMID: 37126977 DOI: 10.1016/j.sbi.2023.102600] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Revised: 03/30/2023] [Accepted: 03/31/2023] [Indexed: 05/03/2023]
Abstract
We provide an overview of the methods that can be used for protein structure-based evaluation of missense variants. The algorithms can be broadly divided into those that calculate the difference in free energy (ΔΔG) between the wild type and variant structures and those that use structural features to predict the damaging effect of a variant without providing a ΔΔG. A wide range of machine learning approaches have been employed to develop those algorithms. We also discuss challenges and opportunities for variant interpretation in view of the recent breakthrough in three-dimensional structural modelling using deep learning.
Collapse
Affiliation(s)
- Alessia David
- Centre for Integrative Systems Biology and Bioinformatics, Department of Life Sciences, Imperial College London, London, SW7 2AZ, UK.
| | - Michael J E Sternberg
- Centre for Integrative Systems Biology and Bioinformatics, Department of Life Sciences, Imperial College London, London, SW7 2AZ, UK
| |
Collapse
|
13
|
Pennica C, Hanna G, Islam SA, JE Sternberg M, David A. Missense3D-PPI: a web resource to predict the impact of missense variants at protein interfaces using 3D structural data. J Mol Biol 2023. [DOI: 10.1016/j.jmb.2023.168060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
|
14
|
Koşaca M, Yılmazbilek İ, Karaca E. PROT-ON: A structure-based detection of designer PROTein interface MutatiONs. Front Mol Biosci 2023; 10:1063971. [PMID: 36936988 PMCID: PMC10018488 DOI: 10.3389/fmolb.2023.1063971] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Accepted: 01/31/2023] [Indexed: 03/06/2023] Open
Abstract
The mutation-induced changes across protein-protein interfaces have often been observed to lead to severe diseases. Therefore, several computational tools have been developed to predict the impact of such mutations. Among these tools, FoldX and EvoEF1 stand out as fast and accurate alternatives. Expanding on the capabilities of these tools, we have developed the PROT-ON (PROTein-protein interface mutatiONs) framework, which aims at delivering the most critical protein interface mutations that can be used to design new protein binders. To realize this aim, PROT-ON takes the 3D coordinates of a protein dimer as an input. Then, it probes all possible interface mutations on the selected protein partner with EvoEF1 or FoldX. The calculated mutational energy landscape is statistically analyzed to find the most enriching and depleting mutations. Afterward, these extreme mutations are filtered out according to stability and optionally according to evolutionary criteria. The final remaining mutation list is presented to the user as the designer mutation set. Together with this set, PROT-ON provides several residue- and energy-based plots, portraying the synthetic energy landscape of the probed mutations. The stand-alone version of PROT-ON is deposited at https://github.com/CSB-KaracaLab/prot-on. The users can also use PROT-ON through our user-friendly web service http://proton.tools.ibg.edu.tr:8001/ (runs with EvoEF1 only). Considering its speed and the range of analysis provided, we believe that PROT-ON presents a promising means to estimate designer mutations.
Collapse
Affiliation(s)
- Mehdi Koşaca
- Izmir Biomedicine and Genome Center, Dokuz Eylul Health Campus, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
| | - İrem Yılmazbilek
- Izmir Biomedicine and Genome Center, Dokuz Eylul Health Campus, Izmir, Türkiye
- Middle East Technical University, Ankara, Türkiye
| | - Ezgi Karaca
- Izmir Biomedicine and Genome Center, Dokuz Eylul Health Campus, Izmir, Türkiye
- Izmir International Biomedicine and Genome Institute, Dokuz Eylul University, Izmir, Türkiye
- *Correspondence: Ezgi Karaca,
| |
Collapse
|
15
|
Xue B, Li R, Ma H, Rahaman A, Kumar V. Comprehensive mapping of mutations in the C9ORF72 that affect folding and binding to SMCR8 protein. Process Biochem 2022. [DOI: 10.1016/j.procbio.2022.07.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
16
|
Gurusinghe SN, Oppenheimer B, Shifman JM. Cold spots are universal in protein–protein interactions. Protein Sci 2022; 31:e4435. [PMID: 36173158 PMCID: PMC9490803 DOI: 10.1002/pro.4435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Revised: 07/22/2022] [Accepted: 08/26/2022] [Indexed: 12/02/2022]
Abstract
Proteins interact with each other through binding interfaces that differ greatly in size and physico‐chemical properties. Within the binding interface, a few residues called hot spots contribute the majority of the binding free energy and are hence irreplaceable. In contrast, cold spots are occupied by suboptimal amino acids, providing possibility for affinity enhancement through mutations. In this study, we identify cold spots due to cavities and unfavorable charge interactions in multiple protein–protein interactions (PPIs). For our cold spot analysis, we first use a small affinity database of PPIs with known structures and affinities and then expand our search to nearly 4000 homo‐ and heterodimers in the Protein Data Bank (PDB). We observe that cold spots due to cavities are present in nearly all PPIs unrelated to their binding affinity, while unfavorable charge interactions are relatively rare. We also find that most cold spots are located in the periphery of the binding interface, with high‐affinity complexes showing fewer centrally located colds spots than low‐affinity complexes. A larger number of cold spots is also found in non‐cognate interactions compared to their cognate counterparts. Furthermore, our analysis reveals that cold spots are more frequent in homo‐dimeric complexes compared to hetero‐complexes, likely due to symmetry constraints imposed on sequences of homodimers. Finally, we find that glycines, glutamates, and arginines are the most frequent amino acids appearing at cold spot positions. Our analysis emphasizes the importance of cold spot positions to protein evolution and facilitates protein engineering studies directed at enhancing binding affinity and specificity in a wide range of applications.
Collapse
Affiliation(s)
- Sagara N.S. Gurusinghe
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| | - Ben Oppenheimer
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| | - Julia M. Shifman
- Department of Biological Chemistry The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem Jerusalem Israel
| |
Collapse
|
17
|
Liu Y, Yeung WSB, Chiu PCN, Cao D. Computational approaches for predicting variant impact: An overview from resources, principles to applications. Front Genet 2022; 13:981005. [PMID: 36246661 PMCID: PMC9559863 DOI: 10.3389/fgene.2022.981005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 08/08/2022] [Indexed: 11/13/2022] Open
Abstract
One objective of human genetics is to unveil the variants that contribute to human diseases. With the rapid development and wide use of next-generation sequencing (NGS), massive genomic sequence data have been created, making personal genetic information available. Conventional experimental evidence is critical in establishing the relationship between sequence variants and phenotype but with low efficiency. Due to the lack of comprehensive databases and resources which present clinical and experimental evidence on genotype-phenotype relationship, as well as accumulating variants found from NGS, different computational tools that can predict the impact of the variants on phenotype have been greatly developed to bridge the gap. In this review, we present a brief introduction and discussion about the computational approaches for variant impact prediction. Following an innovative manner, we mainly focus on approaches for non-synonymous variants (nsSNVs) impact prediction and categorize them into six classes. Their underlying rationale and constraints, together with the concerns and remedies raised from comparative studies are discussed. We also present how the predictive approaches employed in different research. Although diverse constraints exist, the computational predictive approaches are indispensable in exploring genotype-phenotype relationship.
Collapse
Affiliation(s)
- Ye Liu
- Shenzhen Key Laboratory of Fertility Regulation, Reproductive Medicine Center, The University of Hong Kong-Shenzhen Hospital, Shenzhen, China
| | - William S. B. Yeung
- Shenzhen Key Laboratory of Fertility Regulation, Reproductive Medicine Center, The University of Hong Kong-Shenzhen Hospital, Shenzhen, China
- Department of Obstetrics and Gynaecology, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
| | - Philip C. N. Chiu
- Shenzhen Key Laboratory of Fertility Regulation, Reproductive Medicine Center, The University of Hong Kong-Shenzhen Hospital, Shenzhen, China
- Department of Obstetrics and Gynaecology, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
- *Correspondence: Philip C. N. Chiu, ; Dandan Cao,
| | - Dandan Cao
- Shenzhen Key Laboratory of Fertility Regulation, Reproductive Medicine Center, The University of Hong Kong-Shenzhen Hospital, Shenzhen, China
- *Correspondence: Philip C. N. Chiu, ; Dandan Cao,
| |
Collapse
|
18
|
Identification of the most damaging nsSNPs in the human CFL1 gene and their functional and structural impacts on cofilin-1 protein. Gene 2022; 819:146206. [PMID: 35092861 DOI: 10.1016/j.gene.2022.146206] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2021] [Revised: 11/04/2021] [Accepted: 01/13/2022] [Indexed: 01/28/2023]
Abstract
The cofilin-1 protein, encoded by CFL1, is an actin-binding protein that regulates F-actin depolymerization and nucleation activity through phosphorylation and dephosphorylation. CFL1 has been implicated in the development of neurodegenerative diseases (Alzheimer's disease and Huntington's disease), neuronal migration disorders (lissencephaly, epilepsy, and schizophrenia), and neural tube closure defects. Mutations in CFL1 have been associated with impaired neural crest cell migration and neural tube closure defects. In our study, various computational approaches were utilized to explore single-nucleotide polymorphisms (SNPs) in CFL1. The Variation Viewer and gnomAD databases were used to retrieve CFL1 SNPs, including 46 nonsynonymous SNPs (nsSNPs). The functional and structural annotation of SNPs was performed using 12 sequence-based web applications, which identified 20 nsSNPs as being the most likely to be deleterious or disease-causing. The conservation of cofilin-1 protein structures was illustrated using the ConSurf and PROSITE web servers, which projected the 12 most deleterious nsSNPs onto conserved domains, with the potential to disrupt the protein's functionality. These 12 nsSNPs were selected for protein structure construction, and the DynaMut/DUET servers predicted that the protein variants V7G, L84P, and L99A were the most likely to be damaging to the cofilin-1 protein structure or function. The evaluation of molecular docking studies demonstrated that the L99A and L84P cofilin-1 variants reduce the binding affinity for actin compared with the native cofilin-1 structure, and molecular dynamic simulation studies confirmed that these variants might destabilize the protein structure. The consequences of putative mutations on protein-protein interactions and post-translational modification sites in the cofilin-1 protein structure were analyzed. This study represents the first complete approach to understanding the effects of nsSNPs within the actin-depolymerizing factor/cofilin family, which suggested that SNPs resulting in L84P (rs199716082) and L99A (rs267603119) variants represent significant CFL1 mutations associated with disease development.
Collapse
|
19
|
Chen YC, Chen YH, Wright JD, Lim C. PPI-Hotspot DB: Database of Protein-Protein Interaction Hot Spots. J Chem Inf Model 2022; 62:1052-1060. [PMID: 35147037 DOI: 10.1021/acs.jcim.2c00025] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Single-point mutations of certain residues (so-called hot spots) impair/disrupt protein-protein interactions (PPIs), leading to pathogenesis and drug resistance. Conventionally, a PPI-hot spot is identified when its replacement decreased the binding free energy significantly, generally by ≥2 kcal/mol. The relatively few mutations with such a significant binding free energy drop limited the number of distinct PPI-hot spots. By defining PPI-hot spots based on mutations that have been manually curated in UniProtKB to significantly impair/disrupt PPIs in addition to binding free energy changes, we have greatly expanded the number of distinct PPI-hot spots by an order of magnitude. These experimentally determined PPI-hot spots along with available structures have been collected in a database called PPI-HotspotDB. We have applied the PPI-HotspotDB to create a nonredundant benchmark, PPI-Hotspot+PDBBM, for assessing methods to predict PPI-hot spots using the free structure as input. PPI-HotspotDB will benefit the design of mutagenesis experiments and development of PPI-hot spot prediction methods. The database and benchmark are freely available at https://ppihotspot.limlab.dnsalias.org.
Collapse
Affiliation(s)
- Yao Chi Chen
- Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
| | - Yu-Hsien Chen
- Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
| | - Jon D Wright
- Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
| | - Carmay Lim
- Institute of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan.,Department of Chemistry, National Tsing Hua University, Hsinchu 300, Taiwan
| |
Collapse
|
20
|
Abrusán G, Ascher DB, Inouye M. Known allosteric proteins have central roles in genetic disease. PLoS Comput Biol 2022; 18:e1009806. [PMID: 35139069 PMCID: PMC10138267 DOI: 10.1371/journal.pcbi.1009806] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Revised: 04/27/2023] [Accepted: 01/05/2022] [Indexed: 12/15/2022] Open
Abstract
Allostery is a form of protein regulation, where ligands that bind sites located apart from the active site can modify the activity of the protein. The molecular mechanisms of allostery have been extensively studied, because allosteric sites are less conserved than active sites, and drugs targeting them are more specific than drugs binding the active sites. Here we quantify the importance of allostery in genetic disease. We show that 1) known allosteric proteins are central in disease networks, contribute to genetic disease and comorbidities much more than non-allosteric proteins, and there is an association between being allosteric and involvement in disease; 2) they are enriched in many major disease types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes, or diseases of the central nervous system; 3) variants from cancer genome-wide association studies are enriched near allosteric proteins, indicating their importance to polygenic traits; and 4) the importance of allosteric proteins in disease is due, at least partly, to their central positions in protein-protein interaction networks, and less due to their dynamical properties.
Collapse
Affiliation(s)
- György Abrusán
- Cambridge Baker Systems Genomics Initiative, Department of Public Health and Primary Care, School of Medicine, University of Cambridge, Cambridge, United Kingdom
- * E-mail:
| | - David B. Ascher
- Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom
- Structural Biology and Bioinformatics, Department of Biochemistry, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Australia
| | - Michael Inouye
- Cambridge Baker Systems Genomics Initiative, Department of Public Health and Primary Care, School of Medicine, University of Cambridge, Cambridge, United Kingdom
- Cambridge Baker Systems Genomics Initiative, Baker Heart and Diabetes Institute, Melbourne, Australia
- British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge, United Kingdom
- British Heart Foundation Centre of Research Excellence, University of Cambridge, Cambridge, United Kingdom
- Health Data Research UK Cambridge, Wellcome Genome Campus and University of Cambridge, Cambridge, United Kingdom
- The Alan Turing Institute, London, United Kingdom
| |
Collapse
|
21
|
Xiong D, Lee D, Li L, Zhao Q, Yu H. Implications of disease-related mutations at protein-protein interfaces. Curr Opin Struct Biol 2022; 72:219-225. [PMID: 34959033 PMCID: PMC8863207 DOI: 10.1016/j.sbi.2021.11.012] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 11/01/2021] [Accepted: 11/18/2021] [Indexed: 02/03/2023]
Abstract
Protein-protein interfaces have been attracting great attention owing to their critical roles in protein-protein interactions and the fact that human disease-related mutations are generally enriched in them. Recently, substantial research progress has been made in this field, which has significantly promoted the understanding and treatment of various human diseases. For example, many studies have discovered the properties of disease-related mutations. Besides, as more large-scale experimental data become available, various computational approaches have been proposed to advance our understanding of disease mutations from the data. Here, we overview recent advances in characteristics of disease-related mutations at protein-protein interfaces, mutation effects on protein interactions, and investigation of mutations on specific diseases.
Collapse
Affiliation(s)
- Dapeng Xiong
- Department of Computational Biology, Cornell University, Ithaca, NY 14853, USA,Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY 14853, USA
| | - Dongjin Lee
- Department of Computational Biology, Cornell University, Ithaca, NY 14853, USA,Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY 14853, USA
| | - Le Li
- Department of Computational Biology, Cornell University, Ithaca, NY 14853, USA,Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY 14853, USA
| | - Qiuye Zhao
- Department of Computational Biology, Cornell University, Ithaca, NY 14853, USA,Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY 14853, USA
| | - Haiyuan Yu
- Department of Computational Biology, Cornell University, Ithaca, NY 14853, USA,Weill Institute for Cell and Molecular Biology, Cornell University, Ithaca, NY 14853, USA
| |
Collapse
|
22
|
The properties of human disease mutations at protein interfaces. PLoS Comput Biol 2022; 18:e1009858. [PMID: 35120134 PMCID: PMC8849535 DOI: 10.1371/journal.pcbi.1009858] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Revised: 02/16/2022] [Accepted: 01/24/2022] [Indexed: 12/27/2022] Open
Abstract
The assembly of proteins into complexes and their interactions with other biomolecules are often vital for their biological function. While it is known that mutations at protein interfaces have a high potential to be damaging and cause human genetic disease, there has been relatively little consideration for how this varies between different types of interfaces. Here we investigate the properties of human pathogenic and putatively benign missense variants at homomeric (isologous and heterologous), heteromeric, DNA, RNA and other ligand interfaces, and at different regions in proteins with respect to those interfaces. We find that different types of interfaces vary greatly in their propensity to be associated with pathogenic mutations, with homomeric heterologous and DNA interfaces being particularly enriched in disease. We also find that residues that do not directly participate in an interface, but are close in three-dimensional space, show a significant disease enrichment. Finally, we observe that mutations at different types of interfaces tend to have distinct property changes when undergoing amino acid substitutions associated with disease, and that this is linked to substantial variability in their identification by computational variant effect predictors. Nearly all proteins interact with other molecules as part of their biological function. For example, proteins can interact with other copies of the same type of protein, with different proteins, with DNA, or with small ligand molecules. Many mutations at protein interfaces, the regions of proteins that interact with other molecules, are known to cause human genetic disease. In this study, we first investigate how different types of protein interfaces have different tendencies to be associated with disease. We also show that the closer a mutation is to an interface, the more likely it is to cause disease. Finally, we study how mutations at different types of interfaces tend to be associated with different changes in amino acid properties, which appears to influence our ability to computationally predict the effects of mutations. Ultimately, we hope that consideration of protein interface properties will eventually improve our ability to identify new disease-causing mutations.
Collapse
|
23
|
Probing altered enzyme activity in the biochemical characterization of cancer. Biosci Rep 2022; 42:230680. [PMID: 35048115 PMCID: PMC8819661 DOI: 10.1042/bsr20212002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 01/10/2022] [Accepted: 01/19/2022] [Indexed: 11/30/2022] Open
Abstract
Enzymes have evolved to catalyze their precise reactions at the necessary rates, locations, and time to facilitate our development, to respond to a variety of insults and challenges, and to maintain a healthy, balanced state. Enzymes achieve this extraordinary feat through their unique kinetic parameters, myriad regulatory strategies, and their sensitivity to their surroundings, including substrate concentration and pH. The Cancer Genome Atlas (TCGA) highlights the extraordinary number of ways in which the finely tuned activities of enzymes can be disrupted, contributing to cancer development and progression often due to somatic and/or inherited genetic alterations. Rather than being limited to the domain of enzymologists, kinetic constants such as kcat, Km, and kcat/Km are highly informative parameters that can impact a cancer patient in tangible ways—these parameters can be used to sort tumor driver mutations from passenger mutations, to establish the pathways that cancer cells rely on to drive patients’ tumors, to evaluate the selectivity and efficacy of anti-cancer drugs, to identify mechanisms of resistance to treatment, and more. In this review, we will discuss how changes in enzyme activity, primarily through somatic mutation, can lead to altered kinetic parameters, new activities, or changes in conformation and oligomerization. We will also address how changes in the tumor microenvironment can affect enzymatic activity, and briefly describe how enzymology, when combined with additional powerful tools, and can provide us with tremendous insight into the chemical and molecular mechanisms of cancer.
Collapse
|
24
|
Rajkumar MS, Garg R, Jain M. Genome-wide discovery of DNA polymorphisms via resequencing of chickpea cultivars with contrasting response to drought stress. PHYSIOLOGIA PLANTARUM 2022; 174:e13611. [PMID: 34957568 DOI: 10.1111/ppl.13611] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Revised: 11/29/2021] [Accepted: 12/08/2021] [Indexed: 06/14/2023]
Abstract
Drought stress limits plant growth, resulting in a significant yield loss in chickpea. The diversification in genome sequence and selective sweep of allele(s) in different genotypes of a crop plant may play an important role in the determination of agronomic traits, including drought stress response. We investigated, via whole genome resequencing, the DNA polymorphisms between two sets of chickpea genotypes with contrasting drought stress responses (3 drought-sensitive vs. 6 drought-tolerant). In total, 36,406 single nucleotide polymorphisms (SNPs) and 3407 insertions or deletions (InDels) differentiating drought-sensitive and drought-tolerant chickpea genotypes were identified. Interestingly, most (91%) of these DNA polymorphisms were located in chromosomes 1 and 4. The genes harboring DNA polymorphisms in their promoter and/or coding regions and exhibiting differential expression under control and/or drought stress conditions between/within the drought-sensitive and tolerant genotypes were found implicated in the stress response. Furthermore, we identified DNA polymorphisms within the cis-regulatory motifs in the promoter region of abiotic stress-related and QTL-associated genes, which might contribute to the differential expression of the candidate drought-responsive genes. In addition, we revealed the effect of nonsynonymous SNPs on mutational sensitivity and stability of the encoded proteins. Taken together, we identified DNA polymorphisms having relevance in drought stress response and revealed candidate genes to engineer drought tolerance in chickpea.
Collapse
Affiliation(s)
- Mohan Singh Rajkumar
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rohini Garg
- Department of Life Sciences, School of Natural Sciences, Shiv Nadar University, Gautam Buddha Nagar, India
| | - Mukesh Jain
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
25
|
Rodrigues CHM, Pires DEV, Ascher DB. pdCSM-PPI: Using Graph-Based Signatures to Identify Protein-Protein Interaction Inhibitors. J Chem Inf Model 2021; 61:5438-5445. [PMID: 34719929 DOI: 10.1021/acs.jcim.1c01135] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Protein-protein interactions are promising sites for development of selective drugs; however, they have generally been viewed as challenging targets. Molecules targeting protein-protein interactions tend to be larger and more lipophilic than other drug-like molecules, mimicking the properties of interacting interfaces. Here, we propose a machine learning approach that uses a graph-based representation of small molecules to guide identification of inhibitors modulating protein-protein interactions, pdCSM-PPI. This approach was applied to 21 different PPI targets. We developed interaction-specific models that were able to accurately identify active compounds achieving MCC and F1 scores up to 1, and Pearson's correlations up to 0.87, outperforming previous approaches. Using insights from these individual models, we developed a generic protein-protein interaction modulator predictive model, which accurately predicted IC50 with a Pearson's correlation of 0.64 on a low redundancy blind test. Importantly, we were able to accurately identify active from inactive compounds, achieving an AUC of 0.77 and sensitivity and specificity of 76% and 78%, respectively. We believe pdCSM-PPI will be an important tool to help guide more efficient screening of new PPI inhibitors; it is freely available as an easy-to-use web server and API at http://biosig.unimelb.edu.au/pdcsm_ppi.
Collapse
Affiliation(s)
- Carlos H M Rodrigues
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Parkville 3052, Victoria Australia.,Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne 3004, Victoria, Australia.,School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane 4072, Australia
| | - Douglas E V Pires
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Parkville 3052, Victoria Australia.,Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne 3004, Victoria, Australia.,School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane 4072, Australia.,School of Computing and Information Systems, University of Melbourne, Parkville 3052, Victoria, Australia
| | - David B Ascher
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Parkville 3052, Victoria Australia.,Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne 3004, Victoria, Australia.,School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane 4072, Australia
| |
Collapse
|
26
|
Ozturk K, Carter H. Predicting functional consequences of mutations using molecular interaction network features. Hum Genet 2021; 141:1195-1210. [PMID: 34432150 PMCID: PMC8873243 DOI: 10.1007/s00439-021-02329-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Accepted: 07/31/2021] [Indexed: 12/13/2022]
Abstract
Variant interpretation remains a central challenge for precision medicine. Missense variants are particularly difficult to understand as they change only a single amino acid in a protein sequence yet can have large and varied effects on protein activity. Numerous tools have been developed to identify missense variants with putative disease consequences from protein sequence and structure. However, biological function arises through higher order interactions among proteins and molecules within cells. We therefore sought to capture information about the potential of missense mutations to perturb protein interaction networks by integrating protein structure and interaction data. We developed 16 network-based annotations for missense mutations that provide orthogonal information to features classically used to prioritize variants. We then evaluated them in the context of a proven machine-learning framework for variant effect prediction across multiple benchmark datasets to demonstrate their potential to improve variant classification. Interestingly, network features resulted in larger performance gains for classifying somatic mutations than for germline variants, possibly due to different constraints on what mutations are tolerated at the cellular versus organismal level. Our results suggest that modeling variant potential to perturb context-specific interactome networks is a fruitful strategy to advance in silico variant effect prediction.
Collapse
Affiliation(s)
- Kivilcim Ozturk
- Division of Medical Genetics, Department of Medicine, University of California San Diego, La Jolla, CA, USA.,Bioinformatics and Systems Biology Program, University of California San Diego, La Jolla, CA, USA
| | - Hannah Carter
- Division of Medical Genetics, Department of Medicine, University of California San Diego, La Jolla, CA, USA. .,Bioinformatics and Systems Biology Program, University of California San Diego, La Jolla, CA, USA. .,Moores Cancer Center, University of California San Diego, La Jolla, CA, USA.
| |
Collapse
|
27
|
Rodrigues CHM, Pires DEV, Ascher DB. mmCSM-PPI: predicting the effects of multiple point mutations on protein-protein interactions. Nucleic Acids Res 2021; 49:W417-W424. [PMID: 33893812 PMCID: PMC8262703 DOI: 10.1093/nar/gkab273] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Revised: 03/18/2021] [Accepted: 04/15/2021] [Indexed: 11/16/2022] Open
Abstract
Protein-protein interactions play a crucial role in all cellular functions and biological processes and mutations leading to their disruption are enriched in many diseases. While a number of computational methods to assess the effects of variants on protein-protein binding affinity have been proposed, they are in general limited to the analysis of single point mutations and have been shown to perform poorly on independent test sets. Here, we present mmCSM-PPI, a scalable and effective machine learning model for accurately assessing changes in protein-protein binding affinity caused by single and multiple missense mutations. We expanded our well-established graph-based signatures in order to capture physicochemical and geometrical properties of multiple wild-type residue environments and integrated them with substitution scores and dynamics terms from normal mode analysis. mmCSM-PPI was able to achieve a Pearson's correlation of up to 0.75 (RMSE = 1.64 kcal/mol) under 10-fold cross-validation and 0.70 (RMSE = 2.06 kcal/mol) on a non-redundant blind test, outperforming existing methods. Our method is freely available as a user-friendly and easy-to-use web server and API at http://biosig.unimelb.edu.au/mmcsm_ppi.
Collapse
Affiliation(s)
- Carlos H M Rodrigues
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria, Australia
- Structural Biology and Bioinformatics, Department of Biochemistry and Pharmacology, University of Melbourne, Melbourne, Victoria, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Melbourne, Victoria, Australia
| | - Douglas E V Pires
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria, Australia
- Structural Biology and Bioinformatics, Department of Biochemistry and Pharmacology, University of Melbourne, Melbourne, Victoria, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Melbourne, Victoria, Australia
- School of Computing and Information Systems, University of Melbourne, Melbourne, Victoria, Australia
| | - David B Ascher
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria, Australia
- Structural Biology and Bioinformatics, Department of Biochemistry and Pharmacology, University of Melbourne, Melbourne, Victoria, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, Melbourne, Victoria, Australia
- Department of Biochemistry, University of Cambridge, Cambridge, UK
| |
Collapse
|
28
|
The structure-based cancer-related single amino acid variation prediction. Sci Rep 2021; 11:13599. [PMID: 34193921 PMCID: PMC8245468 DOI: 10.1038/s41598-021-92793-w] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 06/16/2021] [Indexed: 11/09/2022] Open
Abstract
Single amino acid variation (SAV) is an amino acid substitution of the protein sequence that can potentially influence the entire protein structure or function, as well as its binding affinity. Protein destabilization is related to diseases, including several cancers, although using traditional experiments to clarify the relationship between SAVs and cancer uses much time and resources. Some SAV prediction methods use computational approaches, with most predicting SAV-induced changes in protein stability. In this investigation, all SAV characteristics generated from protein sequences, structures and the microenvironment were converted into feature vectors and fed into an integrated predicting system using a support vector machine and genetic algorithm. Critical features were used to estimate the relationship between their properties and cancers caused by SAVs. We describe how we developed a prediction system based on protein sequences and structure that is capable of distinguishing if the SAV is related to cancer or not. The five-fold cross-validation performance of our system is 89.73% for the accuracy, 0.74 for the Matthews correlation coefficient, and 0.81 for the F1 score. We have built an online prediction server, CanSavPre ( http://bioinfo.cmu.edu.tw/CanSavPre/ ), which is expected to become a useful, practical tool for cancer research and precision medicine.
Collapse
|
29
|
Missense3D-DB web catalogue: an atom-based analysis and repository of 4M human protein-coding genetic variants. Hum Genet 2021; 140:805-812. [PMID: 33502607 PMCID: PMC8052235 DOI: 10.1007/s00439-020-02246-z] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Accepted: 12/07/2020] [Indexed: 12/22/2022]
Abstract
The interpretation of human genetic variation is one of the greatest challenges of modern genetics. New approaches are urgently needed to prioritize variants, especially those that are rare or lack a definitive clinical interpretation. We examined 10,136,597 human missense genetic variants from GnomAD, ClinVar and UniProt. We were able to perform large-scale atom-based mapping and phenotype interpretation of 3,960,015 of these variants onto 18,874 experimental and 84,818 in house predicted three-dimensional coordinates of the human proteome. We demonstrate that 14% of amino acid substitutions from the GnomAD database that could be structurally analysed are predicted to affect protein structure (n = 568,548, of which 566,439 rare or extremely rare) and may, therefore, have a yet unknown disease-causing effect. The same is true for 19.0% (n = 6266) of variants of unknown clinical significance or conflicting interpretation reported in the ClinVar database. The results of the structural analysis are available in the dedicated web catalogue Missense3D-DB ( http://missense3d.bc.ic.ac.uk/ ). For each of the 4 M variants, the results of the structural analysis are presented in a friendly concise format that can be included in clinical genetic reports. A detailed report of the structural analysis is also available for the non-experts in structural biology. Population frequency and predictions from SIFT and PolyPhen are included for a more comprehensive variant interpretation. This is the first large-scale atom-based structural interpretation of human genetic variation and offers geneticists and the biomedical community a new approach to genetic variant interpretation.
Collapse
|
30
|
Pathogenic missense protein variants affect different functional pathways and proteomic features than healthy population variants. PLoS Biol 2021; 19:e3001207. [PMID: 33909605 PMCID: PMC8110273 DOI: 10.1371/journal.pbio.3001207] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 05/10/2021] [Accepted: 03/26/2021] [Indexed: 12/27/2022] Open
Abstract
Missense variants are present amongst the healthy population, but some of them are causative of human diseases. A classification of variants associated with “healthy” or “diseased” states is therefore not always straightforward. A deeper understanding of the nature of missense variants in health and disease, the cellular processes they may affect, and the general molecular principles which underlie these differences is essential to offer mechanistic explanations of the true impact of pathogenic variants. Here, we have formalised a statistical framework which enables robust probabilistic quantification of variant enrichment across full-length proteins, their domains, and 3D structure-defined regions. Using this framework, we validate and extend previously reported trends of variant enrichment in different protein structural regions (surface/core/interface). By examining the association of variant enrichment with available functional pathways and transcriptomic and proteomic (protein half-life, thermal stability, abundance) data, we have mined a rich set of molecular features which distinguish between pathogenic and population variants: Pathogenic variants mainly affect proteins involved in cell proliferation and nucleotide processing and are enriched in more abundant proteins. Additionally, rare population variants display features closer to common than pathogenic variants. We validate the association between these molecular features and variant pathogenicity by comparing against existing in silico variant impact annotations. This study provides molecular details into how different proteins exhibit resilience and/or sensitivity towards missense variants and provides the rationale to prioritise variant-enriched proteins and protein domains for therapeutic targeting and development. The ZoomVar database, which we created for this study, is available at fraternalilab.kcl.ac.uk/ZoomVar. It allows users to programmatically annotate missense variants with protein structural information and to calculate variant enrichment in different protein structural regions. How do can one improve the classification of genetic variants as harmful or harmless? This study uses a robust statistical analysis to exploit the interplay between protein structure, proteomic measurements and functional pathways to enable better discrimination between missense variants in health and disease.
Collapse
|
31
|
Sitani D, Giorgetti A, Alfonso-Prieto M, Carloni P. Robust principal component analysis-based prediction of protein-protein interaction hot spots. Proteins 2021; 89:639-647. [PMID: 33458895 DOI: 10.1002/prot.26047] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Revised: 12/28/2020] [Accepted: 12/31/2020] [Indexed: 12/21/2022]
Abstract
Proteins often exert their function by binding to other cellular partners. The hot spots are key residues for protein-protein binding. Their identification may shed light on the impact of disease associated mutations on protein complexes and help design protein-protein interaction inhibitors for therapy. Unfortunately, current machine learning methods to predict hot spots, suffer from limitations caused by gross errors in the data matrices. Here, we present a novel data pre-processing pipeline that overcomes this problem by recovering a low rank matrix with reduced noise using Robust Principal Component Analysis. Application to existing databases shows the predictive power of the method.
Collapse
Affiliation(s)
- Divya Sitani
- JARA-Institute: Molecular Neuroscience and Neuroimaging, Institute for Neuroscience and Medicine INM-11/JARA-BRAIN Institute JBI-2, Forschungszentrum Jülich GmbH, Jülich, Germany.,Department of Biology, RWTH Aachen University, Aachen, Germany
| | - Alejandro Giorgetti
- Institute for Advanced Simulations IAS-5 / Institute for Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich GmbH, Jülich, Germany.,Department of Biotechnology, University of Verona, Verona, Italy
| | - Mercedes Alfonso-Prieto
- Institute for Advanced Simulations IAS-5 / Institute for Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich GmbH, Jülich, Germany.,Cécile and Oskar Vogt Institute for Brain Research, University Hospital Düsseldorf, Medical Faculty, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Paolo Carloni
- JARA-Institute: Molecular Neuroscience and Neuroimaging, Institute for Neuroscience and Medicine INM-11/JARA-BRAIN Institute JBI-2, Forschungszentrum Jülich GmbH, Jülich, Germany.,Institute for Advanced Simulations IAS-5 / Institute for Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich GmbH, Jülich, Germany.,Department of Physics, RWTH Aachen University, Aachen, Germany.,JARA-HPC, IAS-5/INM-9 Computational Biomedicine, Forschungszentrum Jülich GmbH, Jülich, Germany
| |
Collapse
|
32
|
Bermúdez-Guzmán L, Jimenez-Huezo G, Arguedas A, Leal A. Mutational survivorship bias: The case of PNKP. PLoS One 2020; 15:e0237682. [PMID: 33332469 PMCID: PMC7746193 DOI: 10.1371/journal.pone.0237682] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2020] [Accepted: 11/23/2020] [Indexed: 01/21/2023] Open
Abstract
The molecular function of a protein relies on its structure. Understanding how variants alter structure and function in multidomain proteins is key to elucidate the generation of a pathological phenotype. However, one may fall into the logical bias of assessing protein damage only based on the variants that are visible (survivorship bias), which can lead to partial conclusions. This is the case of PNKP, an important nuclear and mitochondrial DNA repair enzyme with both kinase and phosphatase function. Most variants in PNKP are confined to the kinase domain, leading to a pathological spectrum of three apparently distinct clinical entities. Since proteins and domains may have a different tolerability to variation, we evaluated whether variants in PNKP are under survivorship bias. Here, we provide the evidence that supports a higher tolerance in the kinase domain even when all variants reported are deleterious. Instead, the phosphatase domain is less tolerant due to its lower variant rates, a higher degree of sequence conservation, lower dN/dS ratios, and the presence of more disease-propensity hotspots. Together, our results support previous experimental evidence that demonstrated that the phosphatase domain is functionally more necessary and relevant for DNA repair, especially in the context of the development of the central nervous system. Finally, we propose the term "Wald’s domain" for future studies analyzing the possible survivorship bias in multidomain proteins.
Collapse
Affiliation(s)
- Luis Bermúdez-Guzmán
- Section of Genetics and Biotechnology, School of Biology, University de Costa Rica, San Pedro, San José, Costa Rica
| | - Gabriel Jimenez-Huezo
- Section of Genetics and Biotechnology, School of Biology, University de Costa Rica, San Pedro, San José, Costa Rica
| | - Andrés Arguedas
- School of Statistics, University de Costa Rica, San Pedro, San José, Costa Rica
| | - Alejandro Leal
- Section of Genetics and Biotechnology, School of Biology, University de Costa Rica, San Pedro, San José, Costa Rica
| |
Collapse
|
33
|
Fang Z, Chen S, Manchanda Y, Bitsi S, Pickford P, David A, Shchepinova MM, Corrêa Jr IR, Hodson DJ, Broichhagen J, Tate EW, Reimann F, Salem V, Rutter GA, Tan T, Bloom SR, Tomas A, Jones B. Ligand-Specific Factors Influencing GLP-1 Receptor Post-Endocytic Trafficking and Degradation in Pancreatic Beta Cells. Int J Mol Sci 2020; 21:E8404. [PMID: 33182425 PMCID: PMC7664906 DOI: 10.3390/ijms21218404] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 11/01/2020] [Accepted: 11/02/2020] [Indexed: 12/30/2022] Open
Abstract
The glucagon-like peptide-1 receptor (GLP-1R) is an important regulator of blood glucose homeostasis. Ligand-specific differences in membrane trafficking of the GLP-1R influence its signalling properties and therapeutic potential in type 2 diabetes. Here, we have evaluated how different factors combine to control the post-endocytic trafficking of GLP-1R to recycling versus degradative pathways. Experiments were performed in primary islet cells, INS-1 832/3 clonal beta cells and HEK293 cells, using biorthogonal labelling of GLP-1R to determine its localisation and degradation after treatment with GLP-1, exendin-4 and several further GLP-1R agonist peptides. We also characterised the effect of a rare GLP1R coding variant, T149M, and the role of endosomal peptidase endothelin-converting enzyme-1 (ECE-1), in GLP1R trafficking. Our data reveal how treatment with GLP-1 versus exendin-4 is associated with preferential GLP-1R targeting towards a recycling pathway. GLP-1, but not exendin-4, is a substrate for ECE-1, and the resultant propensity to intra-endosomal degradation, in conjunction with differences in binding affinity, contributes to alterations in GLP-1R trafficking behaviours and degradation. The T149M GLP-1R variant shows reduced signalling and internalisation responses, which is likely to be due to disruption of the cytoplasmic region that couples to intracellular effectors. These observations provide insights into how ligand- and genotype-specific factors can influence GLP-1R trafficking.
Collapse
Affiliation(s)
- Zijian Fang
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
- Wellcome Trust–Medical Research Council Cambridge Stem Cell Institute and Department of Haematology, University of Cambridge, Cambridge CB2 0AW, UK
| | - Shiqian Chen
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| | - Yusman Manchanda
- Section of Cell Biology and Functional Genomics, Imperial College London, London W12 0NN, UK; (Y.M.); (S.B.); (G.A.R.)
| | - Stavroula Bitsi
- Section of Cell Biology and Functional Genomics, Imperial College London, London W12 0NN, UK; (Y.M.); (S.B.); (G.A.R.)
| | - Philip Pickford
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| | - Alessia David
- Centre for Bioinformatics and System Biology, Department of Life Sciences, Imperial College London, London SW7 2BX, UK;
| | - Maria M. Shchepinova
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 80 Wood Lane, London W12 0BZ, UK; (M.M.S.); (E.W.T.)
| | | | - David J. Hodson
- Institute of Metabolism and Systems Research (IMSR), and Centre of Membrane Proteins and Receptors (COMPARE), University of Birmingham, Birmingham B15 2TT, UK;
- Centre for Endocrinology, Diabetes and Metabolism, Birmingham Health Partners, Birmingham B15 2TT, UK
| | - Johannes Broichhagen
- Department of Chemical Biology, Leibniz-Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany;
- Department of Chemical Biology, Max Planck Institute for Medical Research, 69120 Heidelberg, Germany
| | - Edward W. Tate
- Department of Chemistry, Imperial College London, Molecular Sciences Research Hub, 80 Wood Lane, London W12 0BZ, UK; (M.M.S.); (E.W.T.)
| | - Frank Reimann
- Institute of Metabolic Science & MRC Metabolic Diseases Unit, University of Cambridge, Addenbrooke’s Hospital, Hills Road, Cambridge CB2 0QQ, UK;
| | - Victoria Salem
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| | - Guy A. Rutter
- Section of Cell Biology and Functional Genomics, Imperial College London, London W12 0NN, UK; (Y.M.); (S.B.); (G.A.R.)
- Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore 639798, Singapore
| | - Tricia Tan
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| | - Stephen R. Bloom
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| | - Alejandra Tomas
- Section of Cell Biology and Functional Genomics, Imperial College London, London W12 0NN, UK; (Y.M.); (S.B.); (G.A.R.)
| | - Ben Jones
- Section of Endocrinology and Investigative Medicine, Imperial College London, London W12 0NN, UK; (Z.F.); (S.C.); (P.P.); (V.S.); (T.T.); (S.R.B.)
| |
Collapse
|
34
|
Khalighi S, Singh S, Varadan V. Untangling a complex web: Computational analyses of tumor molecular profiles to decode driver mechanisms. J Genet Genomics 2020; 47:595-609. [PMID: 33423960 PMCID: PMC7902422 DOI: 10.1016/j.jgg.2020.11.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2020] [Revised: 11/04/2020] [Accepted: 11/14/2020] [Indexed: 12/19/2022]
Abstract
Genome-scale studies focusing on molecular profiling of cancers across tissue types have revealed a plethora of aberrations across the genomic, transcriptomic, and epigenomic scales. The significant molecular heterogeneity across individual tumors even within the same tissue context complicates decoding the key etiologic mechanisms of this disease. Furthermore, it is increasingly likely that biologic mechanisms underlying the pathobiology of cancer involve multiple molecular entities interacting across functional scales. This has motivated the development of computational approaches that integrate molecular measurements with prior biological knowledge in increasingly intricate ways to enable the discovery of driver genomic aberrations across cancers. Here, we review diverse methodological approaches that have powered significant advances in our understanding of the genomic underpinnings of cancer at the cohort and at the individual tumor scales. We outline the key advances and challenges in the computational discovery of cancer mechanisms while motivating the development of systems biology approaches to comprehensively decode the biologic drivers of this complex disease.
Collapse
Affiliation(s)
- Sirvan Khalighi
- Division of General Medical Sciences-Oncology, Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Salendra Singh
- Division of General Medical Sciences-Oncology, Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA
| | - Vinay Varadan
- Division of General Medical Sciences-Oncology, Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH 44106, USA.
| |
Collapse
|
35
|
Gerasimavicius L, Liu X, Marsh JA. Identification of pathogenic missense mutations using protein stability predictors. Sci Rep 2020; 10:15387. [PMID: 32958805 PMCID: PMC7506547 DOI: 10.1038/s41598-020-72404-w] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Accepted: 08/31/2020] [Indexed: 12/17/2022] Open
Abstract
Attempts at using protein structures to identify disease-causing mutations have been dominated by the idea that most pathogenic mutations are disruptive at a structural level. Therefore, computational stability predictors, which assess whether a mutation is likely to be stabilising or destabilising to protein structure, have been commonly used when evaluating new candidate disease variants, despite not having been developed specifically for this purpose. We therefore tested 13 different stability predictors for their ability to discriminate between pathogenic and putatively benign missense variants. We find that one method, FoldX, significantly outperforms all other predictors in the identification of disease variants. Moreover, we demonstrate that employing predicted absolute energy change scores improves performance of nearly all predictors in distinguishing pathogenic from benign variants. Importantly, however, we observe that the utility of computational stability predictors is highly heterogeneous across different proteins, and that they are all inferior to the best performing variant effect predictors for identifying pathogenic mutations. We suggest that this is largely due to alternate molecular mechanisms other than protein destabilisation underlying many pathogenic mutations. Thus, better ways of incorporating protein structural information and molecular mechanisms into computational variant effect predictors will be required for improved disease variant prioritisation.
Collapse
Affiliation(s)
- Lukas Gerasimavicius
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Xin Liu
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK.
| |
Collapse
|
36
|
Insights into changes in binding affinity caused by disease mutations in protein-protein complexes. Comput Biol Med 2020; 123:103829. [DOI: 10.1016/j.compbiomed.2020.103829] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Revised: 05/20/2020] [Accepted: 05/20/2020] [Indexed: 01/11/2023]
|
37
|
Protein-Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations. Biomolecules 2020; 10:biom10081097. [PMID: 32722039 PMCID: PMC7463635 DOI: 10.3390/biom10081097] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2020] [Revised: 07/15/2020] [Accepted: 07/20/2020] [Indexed: 12/27/2022] Open
Abstract
Because proteins are fundamental to most biological processes, many genetic diseases can be traced back to single nucleotide variants (SNVs) that cause changes in protein sequences. However, not all SNVs that result in amino acid substitutions cause disease as each residue is under different structural and functional constraints. Influential studies have shown that protein–protein interaction interfaces are enriched in disease-associated SNVs and depleted in SNVs that are common in the general population. These studies focus primarily on folded (globular) protein domains and overlook the prevalent class of protein interactions mediated by intrinsically disordered regions (IDRs). Therefore, we investigated the enrichment patterns of missense mutation-causing SNVs that are associated with disease and cancer, as well as those present in the healthy population, in structures of IDR-mediated interactions with comparisons to classical globular interactions. When comparing the different categories of interaction interfaces, division of the interface regions into solvent-exposed rim residues and buried core residues reveal distinctive enrichment patterns for the various types of missense mutations. Most notably, we demonstrate a strong enrichment at the interface core of interacting IDRs in disease mutations and its depletion in neutral ones, which supports the view that the disruption of IDR interactions is a mechanism underlying many diseases. Intriguingly, we also found an asymmetry across the IDR interaction interface in the enrichment of certain missense mutation types, which may hint at an increased variant tolerance and urges further investigations of IDR interactions.
Collapse
|
38
|
Nisar H, Pasha U, Mirza MU, Abid R, Hanif K, Kadarmideen HN, Sadaf S. Impact of IL-17F 7488T/C Functional Polymorphism on Progressive Rheumatoid Arthritis: Novel Insight from the Molecular Dynamic Simulations. Immunol Invest 2020; 50:416-426. [PMID: 32543936 DOI: 10.1080/08820139.2020.1775642] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Resorption of bones and cartilage coupled with structural changes in the inflamed joints are the major hallmark of rheumatoid arthritis (RA). Genetic polymorphisms in pro-inflammatory interleukins (ILs) appear to play an important role in the susceptibility towards progressive RA. We therefore aimed to investigate the association of single nucleotide polymorphisms (SNP), present in the hotspot coding/promoter regions of IL-6, -17 and -18, with RA susceptibility or severity in a larger study cohort from Pakistan together with finding clues as to how a functional SNP impacts the predisposition towards RA. TaqMan SNP genotyping approach was first used to assess IL-6 (rs1800795), IL-17 F (rs763780), IL-17A (rs2275913), and IL-18 (rs1946518) polymorphisms in 310 subjects (150 RA and 160 control). Molecular dynamic simulations (MDS) of wild- and mutant-type IL-17A with corresponding receptor were thereafter performed using AMBER-16; Chimera 1.13 was used for analyses. Our results showed the association of two SNPs, namely IL-6 - 174 G/C [allelic (OR = 0.960, 95% CI = 0.929-0.992, p = .009)] and IL-17 F 7488 T/C [allelic (OR = 0.907, 95%CI = 0.861-0.954, p = .000)] with increased RA risk in Pakistani subjects. When mapped, IL-17 F 7488 T/C was found involved in His161→Arg161 change near the C-terminus of IL-17 F. Comparative MDS revealed enhanced stability of the mutant hence advocating a potential role of IL-17F functional SNP in RA susceptibility and/or severity. This study provides a novel structural insight for SNP-derived functional mutation and its overall impact on binding with heterotrimeric receptor complex of IL-17 receptor thereby opening new avenues for understanding the biochemical basis of the disease.
Collapse
Affiliation(s)
- Haseeb Nisar
- Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore, Pakistan.,Quantitative Genomics, Bioinformatics and Computational Biology Group, DTU Compute, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Usman Pasha
- Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore, Pakistan
| | - Muhammad Usman Mirza
- Department of Pharmaceutical and Pharmacological Sciences, Rega Institute for Medical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium
| | - Rizwan Abid
- Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore, Pakistan
| | - Kiran Hanif
- Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore, Pakistan
| | - Haja N Kadarmideen
- Quantitative Genomics, Bioinformatics and Computational Biology Group, DTU Compute, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Saima Sadaf
- Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore, Pakistan
| |
Collapse
|
39
|
Fleming JR, Rigden DJ, Mayans O. The importance of chain context in assessing small nucleotide variants in titin: in silico case study of the I10-I11 tandem and its arrhythmic right ventricular cardiomyopathy linked position T2580. J Biomol Struct Dyn 2020; 39:3480-3490. [PMID: 32396765 DOI: 10.1080/07391102.2020.1768148] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Non-synonymous small nucleotide variations (nsSNVs) in the giant muscle protein, titin, have key roles in the development of several myopathologies. Although there is considerable motive to screen at-risk individuals for nsSNVs, to identify patients in early disease stages while therapeutic intervention is still possible, the clinical significance of most titin variations remains unclear. Therefore, there is a growing need to establish methods to classify nsSNVs in a simple, economic and rapid manner. Due to its strong correlation to arrhythmogenic right ventricular cardiomyopathy (ARVC), one particular mutation in titin-T2580I, located in the I10 immunoglobulin domain-has received considerable attention. Here, we use the I10-I11 tandem as a case study to explore the possible benefits of considering the titin chain context-i.e. domain interfaces-in the assessment of titin nsSNVs. Specifically, we investigate which exchanges mimic the conformational molecular phenotype of the T2580I mutation at the I10-I11 domain interface. Then, we computed a residue stability landscape for domains alone and in tandem to define a Domain Interface Score (DIS) which identifies several hotspot residues. Our findings suggest that the T2580 position is highly sensitive to exchange and that any variant found in this position should be considered with care. Furthermore, we conclude that the consideration of the higher order structure of the titin chain is important to gain accurate insights into the vulnerability of positions in linker regions and that titin nsSNV prediction benefits from a contextual analysis. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- J R Fleming
- Department of Biology, University of Konstanz, Konstanz, Germany
| | - D J Rigden
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, UK
| | - O Mayans
- Department of Biology, University of Konstanz, Konstanz, Germany
| |
Collapse
|
40
|
Udhaya Kumar S, Thirumal Kumar D, Mandal PD, Sankar S, Haldar R, Kamaraj B, Walter CEJ, Siva R, George Priya Doss C, Zayed H. Comprehensive in silico screening and molecular dynamics studies of missense mutations in Sjogren-Larsson syndrome associated with the ALDH3A2 gene. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2020; 120:349-377. [PMID: 32085885 DOI: 10.1016/bs.apcsb.2019.11.004] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.
Collapse
Affiliation(s)
- S Udhaya Kumar
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - D Thirumal Kumar
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Pinky D Mandal
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Srivarshini Sankar
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Rishin Haldar
- School of Computer Science and Engineering, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Balu Kamaraj
- Department of Neuroscience Technology, College of Applied Medical Sciences, Imam Abdulrahman Bin Faisal University, Jubail, Saudi Arabia
| | - Charles Emmanuel Jebaraj Walter
- Department of Biotechnology, Sri Ramachandra Institute of Higher Education and Research (Deemed to be University), Chennai, India
| | - R Siva
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - C George Priya Doss
- School of BioSciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - Hatem Zayed
- Department of Biomedical Sciences, College of Health and Sciences, Qatar University, Doha, Qatar
| |
Collapse
|
41
|
Galano-Frutos JJ, García-Cebollada H, Sancho J. Molecular dynamics simulations for genetic interpretation in protein coding regions: where we are, where to go and when. Brief Bioinform 2019; 22:3-19. [PMID: 31813950 DOI: 10.1093/bib/bbz146] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Revised: 09/22/2019] [Accepted: 10/25/2019] [Indexed: 12/18/2022] Open
Abstract
The increasing ease with which massive genetic information can be obtained from patients or healthy individuals has stimulated the development of interpretive bioinformatics tools as aids in clinical practice. Most such tools analyze evolutionary information and simple physical-chemical properties to predict whether replacement of one amino acid residue with another will be tolerated or cause disease. Those approaches achieve up to 80-85% accuracy as binary classifiers (neutral/pathogenic). As such accuracy is insufficient for medical decision to be based on, and it does not appear to be increasing, more precise methods, such as full-atom molecular dynamics (MD) simulations in explicit solvent, are also discussed. Then, to describe the goal of interpreting human genetic variations at large scale through MD simulations, we restrictively refer to all possible protein variants carrying single-amino-acid substitutions arising from single-nucleotide variations as the human variome. We calculate its size and develop a simple model that allows calculating the simulation time needed to have a 0.99 probability of observing unfolding events of any unstable variant. The knowledge of that time enables performing a binary classification of the variants (stable-potentially neutral/unstable-pathogenic). Our model indicates that the human variome cannot be simulated with present computing capabilities. However, if they continue to increase as per Moore's law, it could be simulated (at 65°C) spending only 3 years in the task if we started in 2031. The simulation of individual protein variomes is achievable in short times starting at present. International coordination seems appropriate to embark upon massive MD simulations of protein variants.
Collapse
Affiliation(s)
- Juan J Galano-Frutos
- Protein Folding and Molecular Design (ProtMol)' group at BIFI, University of Zaragoza
| | | | - Javier Sancho
- Protein Folding and Molecular Design (ProtMol)' group at BIFI, University of Zaragoza
| |
Collapse
|
42
|
Michalak W, Tsiamis V, Schwämmle V, Rogowska-Wrzesińska A. ComplexBrowser: A Tool for Identification and Quantification of Protein Complexes in Large-scale Proteomics Datasets. Mol Cell Proteomics 2019; 18:2324-2334. [PMID: 31447428 PMCID: PMC6823858 DOI: 10.1074/mcp.tir119.001434] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Revised: 07/27/2019] [Indexed: 12/25/2022] Open
Abstract
We have developed ComplexBrowser, an open source, online platform for supervised analysis of quantitative proteomic data (label free and isobaric mass tag based) that focuses on protein complexes. The software uses manually curated information from CORUM and Complex Portal databases to identify protein complex components. For the first time, we provide a Complex Fold Change (CFC) factor that identifies up- and downregulated complexes based on the level of complex subunits coregulation. The software provides interactive visualization of protein complexes' composition and expression for exploratory analysis and incorporates a quality control step that includes normalization and statistical analysis based on the limma package. ComplexBrowser was tested on two published studies identifying changes in protein expression within either human adenocarcinoma tissue or activated mouse T-cells. The analysis revealed 1519 and 332 protein complexes, of which 233 and 41 were found coordinately regulated in the respective studies. The adopted approach provided evidence for a shift to glucose-based metabolism and high proliferation in adenocarcinoma tissues, and the identification of chromatin remodeling complexes involved in mouse T-cell activation. The results correlate with the original interpretation of the experiments and provide novel biological details about the protein complexes affected. ComplexBrowser is, to our knowledge, the first tool to automate quantitative protein complex analysis for high-throughput studies, providing insights into protein complex regulation within minutes of analysis.
Collapse
Affiliation(s)
- Wojciech Michalak
- Department of Biochemistry & Molecular Biology and VILLUM Center for Bioanalytical Sciences, University of Southern Denmark, Campusvej 55, DK-5230, Odense M, Denmark
| | - Vasileios Tsiamis
- Department of Biochemistry & Molecular Biology and VILLUM Center for Bioanalytical Sciences, University of Southern Denmark, Campusvej 55, DK-5230, Odense M, Denmark
| | - Veit Schwämmle
- Department of Biochemistry & Molecular Biology and VILLUM Center for Bioanalytical Sciences, University of Southern Denmark, Campusvej 55, DK-5230, Odense M, Denmark
| | - Adelina Rogowska-Wrzesińska
- Department of Biochemistry & Molecular Biology and VILLUM Center for Bioanalytical Sciences, University of Southern Denmark, Campusvej 55, DK-5230, Odense M, Denmark.
| |
Collapse
|
43
|
Sandoval JE, Reich NO. The R882H substitution in the human de novo DNA methyltransferase DNMT3A disrupts allosteric regulation by the tumor supressor p53. J Biol Chem 2019; 294:18207-18219. [PMID: 31640986 DOI: 10.1074/jbc.ra119.010827] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Revised: 10/18/2019] [Indexed: 12/14/2022] Open
Abstract
A myriad of protein partners modulate the activity of the human DNA methyltransferase 3A (DNMT3A), whose interactions with these other proteins are frequently altered during oncogenesis. We show here that the tumor suppressor p53 decreases DNMT3A activity by forming a heterotetramer complex with DNMT3A. Mutational and modeling experiments suggested that p53 interacts with the same region in DNMT3A as does the structurally characterized DNMT3L. We observed that the p53-mediated repression of DNMT3A activity is blocked by amino acid substitutions within this interface, but surprisingly, also by a distal DNMT3A residue, R882H. DNMT3A R882H occurs frequently in various cancers, including acute myeloid leukemia, and our results suggest that the effects of R882H and other DNMT3A mutations may go beyond changes in DNMT3A methylation activity. To further understand the dynamics of how protein-protein interactions modulate DNMT3A activity, we determined that p53 has a greater affinity for DNMT3A than for DNMT3L and that p53 readily displaces DNMT3L from the DNMT3A:DNMT3L heterotetramer. Interestingly, this occurred even when the preformed DNMT3A:DNMT3L complex was actively methylating DNA. The frequently identified p53 substitutions (R248W and R273H), whereas able to regulate DNMT3A function when forming the DNMT3A:p53 heterotetramer, no longer displaced DNMT3L from the DNMT3A:DNMT3L heterotetramer. The results of our work highlight the complex interplay between DNMT3A, p53, and DNMT3L and how these interactions are further modulated by clinically derived mutations in each of the interacting partners.
Collapse
Affiliation(s)
- Jonathan E Sandoval
- Department of Molecular, Cellular and Developmental Biology, University of California, Santa Barbara, California 93106-9510
| | - Norbert O Reich
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510.
| |
Collapse
|
44
|
Ozdemir ES, Gursoy A, Keskin O. Analysis of single amino acid variations in singlet hot spots of protein-protein interfaces. Bioinformatics 2019; 34:i795-i801. [PMID: 30423104 DOI: 10.1093/bioinformatics/bty569] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Motivation Single amino acid variations (SAVs) in protein-protein interaction (PPI) sites play critical roles in diseases. PPI sites (interfaces) have a small subset of residues called hot spots that contribute significantly to the binding energy, and they may form clusters called hot regions. Singlet hot spots are the single amino acid hot spots outside of the hot regions. The distribution of SAVs on the interface residues may be related to their disease association. Results We performed statistical and structural analyses of SAVs with literature curated experimental thermodynamics data, and demonstrated that SAVs which destabilize PPIs are more likely to be found in singlet hot spots rather than hot regions and energetically less important interface residues. In contrast, non-hot spot residues are significantly enriched in neutral SAVs, which do not affect PPI stability. Surprisingly, we observed that singlet hot spots tend to be enriched in disease-causing SAVs, while benign SAVs significantly occur in non-hot spot residues. Our work demonstrates that SAVs in singlet hot spot residues have significant effect on protein stability and function. Availability and implementation The dataset used in this paper is available as Supplementary Material. The data can be found at http://prism.ccbb.ku.edu.tr/data/sav/ as well. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- E Sila Ozdemir
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey
| | - Attila Gursoy
- Department of Computer Engineering, Koc University, Istanbul, Turkey.,Research Center for Translational Medicine (KUTTAM), Koc University, Istanbul, Turkey
| | - Ozlem Keskin
- Department of Chemical and Biological Engineering, Koc University, Istanbul, Turkey.,Research Center for Translational Medicine (KUTTAM), Koc University, Istanbul, Turkey
| |
Collapse
|
45
|
Rodrigues CHM, Myung Y, Pires DEV, Ascher DB. mCSM-PPI2: predicting the effects of mutations on protein-protein interactions. Nucleic Acids Res 2019; 47:W338-W344. [PMID: 31114883 PMCID: PMC6602427 DOI: 10.1093/nar/gkz383] [Citation(s) in RCA: 200] [Impact Index Per Article: 40.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Revised: 04/30/2019] [Accepted: 05/20/2019] [Indexed: 12/13/2022] Open
Abstract
Protein-protein Interactions are involved in most fundamental biological processes, with disease causing mutations enriched at their interfaces. Here we present mCSM-PPI2, a novel machine learning computational tool designed to more accurately predict the effects of missense mutations on protein-protein interaction binding affinity. mCSM-PPI2 uses graph-based structural signatures to model effects of variations on the inter-residue interaction network, evolutionary information, complex network metrics and energetic terms to generate an optimised predictor. We demonstrate that our method outperforms previous methods, ranking first among 26 others on CAPRI blind tests. mCSM-PPI2 is freely available as a user friendly webserver at http://biosig.unimelb.edu.au/mcsm_ppi2/.
Collapse
Affiliation(s)
- Carlos H M Rodrigues
- Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, Australia
- ACRF Facility for Innovative Cancer Drug Discovery, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Structural Biology and Bioinformatics, Baker Heart and Diabetes Institute, Melbourne, Australia
| | - Yoochan Myung
- Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, Australia
- ACRF Facility for Innovative Cancer Drug Discovery, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Structural Biology and Bioinformatics, Baker Heart and Diabetes Institute, Melbourne, Australia
| | - Douglas E V Pires
- Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, Australia
- ACRF Facility for Innovative Cancer Drug Discovery, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Structural Biology and Bioinformatics, Baker Heart and Diabetes Institute, Melbourne, Australia
| | - David B Ascher
- Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, Australia
- ACRF Facility for Innovative Cancer Drug Discovery, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Structural Biology and Bioinformatics, Baker Heart and Diabetes Institute, Melbourne, Australia
- Department of Biochemistry, University of Cambridge, Cambridge, UK
| |
Collapse
|
46
|
Bergendahl LT, Gerasimavicius L, Miles J, Macdonald L, Wells JN, Welburn JPI, Marsh JA. The role of protein complexes in human genetic disease. Protein Sci 2019; 28:1400-1411. [PMID: 31219644 DOI: 10.1002/pro.3667] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Accepted: 06/10/2019] [Indexed: 12/20/2022]
Abstract
Many human genetic disorders are caused by mutations in protein-coding regions of DNA. Taking protein structure into account has therefore provided key insight into the molecular mechanisms underlying human genetic disease. Although most studies have focused on the intramolecular effects of mutations, the critical role of the assembly of proteins into complexes is being increasingly recognized. Here, we review multiple ways in which consideration of protein complexes can help us to understand and explain the effects of pathogenic mutations. First, we discuss disorders caused by mutations that perturb intersubunit interactions in homomeric and heteromeric complexes. Second, we address how protein complex assembly can facilitate a dominant-negative mechanism, whereby mutated subunits can disrupt the activity of wild-type protein. Third, we show how mutations that change protein expression levels can lead to damaging stoichiometric imbalances. Finally, we review how mutations affecting different subunits of the same heteromeric complex often cause similar diseases, whereas mutations in different interfaces of the same subunit can cause distinct phenotypes.
Collapse
Affiliation(s)
- L Therese Bergendahl
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Lukas Gerasimavicius
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Jamilla Miles
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Lewis Macdonald
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| | - Jonathan N Wells
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York, 14850
| | - Julie P I Welburn
- Wellcome Trust Centre for Cell Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3BF, United Kingdom
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, United Kingdom
| |
Collapse
|
47
|
Konc J, Skrlj B, Erzen N, Kunej T, Janezic D. GenProBiS: web server for mapping of sequence variants to protein binding sites. Nucleic Acids Res 2019; 45:W253-W259. [PMID: 28498966 PMCID: PMC5570222 DOI: 10.1093/nar/gkx420] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Accepted: 05/02/2017] [Indexed: 02/02/2023] Open
Abstract
Discovery of potentially deleterious sequence variants is important and has wide implications for research and generation of new hypotheses in human and veterinary medicine, and drug discovery. The GenProBiS web server maps sequence variants to protein structures from the Protein Data Bank (PDB), and further to protein–protein, protein–nucleic acid, protein–compound, and protein–metal ion binding sites. The concept of a protein–compound binding site is understood in the broadest sense, which includes glycosylation and other post-translational modification sites. Binding sites were defined by local structural comparisons of whole protein structures using the Protein Binding Sites (ProBiS) algorithm and transposition of ligands from the similar binding sites found to the query protein using the ProBiS-ligands approach with new improvements introduced in GenProBiS. Binding site surfaces were generated as three-dimensional grids encompassing the space occupied by predicted ligands. The server allows intuitive visual exploration of comprehensively mapped variants, such as human somatic mis-sense mutations related to cancer and non-synonymous single nucleotide polymorphisms from 21 species, within the predicted binding sites regions for about 80 000 PDB protein structures using fast WebGL graphics. The GenProBiS web server is open and free to all users at http://genprobis.insilab.org.
Collapse
Affiliation(s)
- Janez Konc
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia.,University of Primorska, Faculty of Mathematics, Natural Sciences and Information Technologies, 6000 Koper, Slovenia
| | - Blaz Skrlj
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Nika Erzen
- National Institute of Chemistry, Hajdrihova 19, 1000 Ljubljana, Slovenia
| | - Tanja Kunej
- Biotechnical Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia
| | - Dusanka Janezic
- University of Primorska, Faculty of Mathematics, Natural Sciences and Information Technologies, 6000 Koper, Slovenia
| |
Collapse
|
48
|
Shim JE, Kim JH, Shin J, Lee JE, Lee I. Pathway-specific protein domains are predictive for human diseases. PLoS Comput Biol 2019; 15:e1007052. [PMID: 31075101 PMCID: PMC6530867 DOI: 10.1371/journal.pcbi.1007052] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Revised: 05/22/2019] [Accepted: 04/19/2019] [Indexed: 01/04/2023] Open
Abstract
Protein domains are basic functional units of proteins. Many protein domains are pervasive among diverse biological processes, yet some are associated with specific pathways. Human complex diseases are generally viewed as pathway-level disorders. Therefore, we hypothesized that pathway-specific domains could be highly informative for human diseases. To test the hypothesis, we developed a network-based scoring scheme to quantify specificity of domain-pathway associations. We first generated domain profiles for human proteins, then constructed a co-pathway protein network based on the associations between domain profiles. Based on the score, we classified human protein domains into pathway-specific domains (PSDs) and non-specific domains (NSDs). We found that PSDs contained more pathogenic variants than NSDs. PSDs were also enriched for disease-associated mutations that disrupt protein-protein interactions (PPIs) and tend to have a moderate number of domain interactions. These results suggest that mutations in PSDs are likely to disrupt within-pathway PPIs, resulting in functional failure of pathways. Finally, we demonstrated the prediction capacity of PSDs for disease-associated genes with experimental validations in zebrafish. Taken together, the network-based quantitative method of modeling domain-pathway associations presented herein suggested underlying mechanisms of how protein domains associated with specific pathways influence mutational impacts on diseases via perturbations in within-pathway PPIs, and provided a novel genomic feature for interpreting genetic variants to facilitate the discovery of human disease genes. Protein domains are basic functional units of proteins, yet domain-based pathway annotations for proteins are challenging tasks because many domains are pervasive among diverse pathways. Therefore, we developed a network-based scoring scheme to measure pathway specificity of domains, and then used it to identify pathway-specific domains. Surprisingly, we observed substantially more disease mutations in pathway-specific domains than non-specific domains. We found evidences that mutations of pathway-specific domains tend to perturb pathway integrity via disrupting within-pathway protein-protein interactions. We also demonstrated prediction capacity of pathway-specific domains for complex diseases with experimental validations. Our study demonstrated the usefulness of pathway information for protein domains in interpreting non-random distribution of disease mutations among domains and identification of disease genes and variants.
Collapse
Affiliation(s)
- Jung Eun Shim
- Department of Biotechnology, College of Life Science and Biotechnology, Yonsei University, Seoul, Korea
- Yonsei Biomedical Research Institute, Yonsei University College of Medicine, Seoul, Korea
| | - Ji Hyun Kim
- Department of Health Sciences and Technology, SAIHST, Sungkyunkwan University, Seoul, Korea
| | - Junha Shin
- Department of Biotechnology, College of Life Science and Biotechnology, Yonsei University, Seoul, Korea
| | - Ji Eun Lee
- Department of Health Sciences and Technology, SAIHST, Sungkyunkwan University, Seoul, Korea
- Samsung Biomedical Research Institute, Samsung Medical Center, Seoul, Korea
| | - Insuk Lee
- Department of Biotechnology, College of Life Science and Biotechnology, Yonsei University, Seoul, Korea
- Department of Biomedical Systems Informatics, Yonsei University College of Medicine, Seoul, Korea
- * E-mail:
| |
Collapse
|
49
|
Can Predicted Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated? J Mol Biol 2019; 431:2197-2212. [PMID: 30995449 PMCID: PMC6544567 DOI: 10.1016/j.jmb.2019.04.009] [Citation(s) in RCA: 283] [Impact Index Per Article: 56.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Revised: 03/18/2019] [Accepted: 04/07/2019] [Indexed: 01/29/2023]
Abstract
Knowledge of protein structure can be used to predict the phenotypic consequence of a missense variant. Since structural coverage of the human proteome can be roughly tripled to over 50% of the residues if homology-predicted structures are included in addition to experimentally determined coordinates, it is important to assess the reliability of using predicted models when analyzing missense variants. Accordingly, we assess whether a missense variant is structurally damaging by using experimental and predicted structures. We considered 606 experimental structures and show that 40% of the 1965 disease-associated missense variants analyzed have a structurally damaging change in the mutant structure. Only 11% of the 2134 neutral variants are structurally damaging. Importantly, similar results are obtained when 1052 structures predicted using Phyre2 algorithm were used, even when the model shares low (<40%) sequence identity to the template. Thus, structure-based analysis of the effects of missense variants can be effectively applied to homology models. Our in-house pipeline, Missense3D, for structurally assessing missense variants was made available at http://www.sbg.bio.ic.ac.uk/~missense3d.
Collapse
|
50
|
Structural and Computational Characterization of Disease-Related Mutations Involved in Protein-Protein Interfaces. Int J Mol Sci 2019; 20:ijms20071583. [PMID: 30934865 PMCID: PMC6479360 DOI: 10.3390/ijms20071583] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2019] [Revised: 03/26/2019] [Accepted: 03/27/2019] [Indexed: 11/24/2022] Open
Abstract
One of the known potential effects of disease-causing amino acid substitutions in proteins is to modulate protein-protein interactions (PPIs). To interpret such variants at the molecular level and to obtain useful information for prediction purposes, it is important to determine whether they are located at protein-protein interfaces, which are composed of two main regions, core and rim, with different evolutionary conservation and physicochemical properties. Here we have performed a structural, energetics and computational analysis of interactions between proteins hosting mutations related to diseases detected in newborn screening. Interface residues were classified as core or rim, showing that the core residues contribute the most to the binding free energy of the PPI. Disease-causing variants are more likely to occur at the interface core region rather than at the interface rim (p < 0.0001). In contrast, neutral variants are more often found at the interface rim or at the non-interacting surface rather than at the interface core region. We also found that arginine, tryptophan, and tyrosine are over-represented among mutated residues leading to disease. These results can enhance our understanding of disease at molecular level and thus contribute towards personalized medicine by helping clinicians to provide adequate diagnosis and treatments.
Collapse
|