1
|
Cheng P, Mao C, Tang J, Yang S, Cheng Y, Wang W, Gu Q, Han W, Chen H, Li S, Chen Y, Zhou J, Li W, Pan A, Zhao S, Huang X, Zhu S, Zhang J, Shu W, Wang S. Zero-shot prediction of mutation effects with multimodal deep representation learning guides protein engineering. Cell Res 2024; 34:630-647. [PMID: 38969803 PMCID: PMC11369238 DOI: 10.1038/s41422-024-00989-2] [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: 03/13/2024] [Accepted: 06/03/2024] [Indexed: 07/07/2024] Open
Abstract
Mutations in amino acid sequences can provoke changes in protein function. Accurate and unsupervised prediction of mutation effects is critical in biotechnology and biomedicine, but remains a fundamental challenge. To resolve this challenge, here we present Protein Mutational Effect Predictor (ProMEP), a general and multiple sequence alignment-free method that enables zero-shot prediction of mutation effects. A multimodal deep representation learning model embedded in ProMEP was developed to comprehensively learn both sequence and structure contexts from ~160 million proteins. ProMEP achieves state-of-the-art performance in mutational effect prediction and accomplishes a tremendous improvement in speed, enabling efficient and intelligent protein engineering. Specifically, ProMEP accurately forecasts mutational consequences on the gene-editing enzymes TnpB and TadA, and successfully guides the development of high-performance gene-editing tools with their engineered variants. The gene-editing efficiency of a 5-site mutant of TnpB reaches up to 74.04% (vs 24.66% for the wild type); and the base editing tool developed on the basis of a TadA 15-site mutant (in addition to the A106V/D108N double mutation that renders deoxyadenosine deaminase activity to TadA) exhibits an A-to-G conversion frequency of up to 77.27% (vs 69.80% for ABE8e, a previous TadA-based adenine base editor) with significantly reduced bystander and off-target effects compared to ABE8e. ProMEP not only showcases superior performance in predicting mutational effects on proteins but also demonstrates a great capability to guide protein engineering. Therefore, ProMEP enables efficient exploration of the gigantic protein space and facilitates practical design of proteins, thereby advancing studies in biomedicine and synthetic biology.
Collapse
Affiliation(s)
- Peng Cheng
- Bioinformatics Center of AMMS, Beijing, China
| | - Cong Mao
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China
| | - Jin Tang
- Zhejiang Lab, Hangzhou, Zhejiang, China
| | - Sen Yang
- Bioinformatics Center of AMMS, Beijing, China
| | - Yu Cheng
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China
| | - Wuke Wang
- Zhejiang Lab, Hangzhou, Zhejiang, China
| | - Qiuxi Gu
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China
| | - Wei Han
- Zhejiang Lab, Hangzhou, Zhejiang, China
| | - Hao Chen
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China
| | - Sihan Li
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China
| | | | | | - Wuju Li
- Bioinformatics Center of AMMS, Beijing, China
| | - Aimin Pan
- Zhejiang Lab, Hangzhou, Zhejiang, China
| | - Suwen Zhao
- iHuman Institute, ShanghaiTech University, Shanghai, China
- School of Life Science and Technology, ShanghaiTech University, Shanghai, China
| | - Xingxu Huang
- Zhejiang Lab, Hangzhou, Zhejiang, China
- School of Life Science and Technology, ShanghaiTech University, Shanghai, China
| | | | - Jun Zhang
- State Key Laboratory of Reproductive Medicine and Offspring Health, Women's Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing Medical University, Nanjing, Jiangsu, China.
| | - Wenjie Shu
- Bioinformatics Center of AMMS, Beijing, China.
| | | |
Collapse
|
2
|
Li SS, Liu ZM, Li J, Ma YB, Dong ZY, Hou JW, Shen FJ, Wang WB, Li QM, Su JG. Prediction of mutation-induced protein stability changes based on the geometric representations learned by a self-supervised method. BMC Bioinformatics 2024; 25:282. [PMID: 39198740 PMCID: PMC11360314 DOI: 10.1186/s12859-024-05876-6] [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: 02/27/2024] [Accepted: 07/19/2024] [Indexed: 09/01/2024] Open
Abstract
BACKGROUND Thermostability is a fundamental property of proteins to maintain their biological functions. Predicting protein stability changes upon mutation is important for our understanding protein structure-function relationship, and is also of great interest in protein engineering and pharmaceutical design. RESULTS Here we present mutDDG-SSM, a deep learning-based framework that uses the geometric representations encoded in protein structure to predict the mutation-induced protein stability changes. mutDDG-SSM consists of two parts: a graph attention network-based protein structural feature extractor that is trained with a self-supervised learning scheme using large-scale high-resolution protein structures, and an eXtreme Gradient Boosting model-based stability change predictor with an advantage of alleviating overfitting problem. The performance of mutDDG-SSM was tested on several widely-used independent datasets. Then, myoglobin and p53 were used as case studies to illustrate the effectiveness of the model in predicting protein stability changes upon mutations. Our results show that mutDDG-SSM achieved high performance in estimating the effects of mutations on protein stability. In addition, mutDDG-SSM exhibited good unbiasedness, where the prediction accuracy on the inverse mutations is as well as that on the direct mutations. CONCLUSION Meaningful features can be extracted from our pre-trained model to build downstream tasks and our model may serve as a valuable tool for protein engineering and drug design.
Collapse
Affiliation(s)
- Shan Shan Li
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China
- National Engineering Center for New Vaccine Research, Beijing, China
| | - Zhao Ming Liu
- National Engineering Center for New Vaccine Research, Beijing, China
- The Sixth Laboratory, National Vaccine and Serum Institute (NVSI), Beijing, China
| | - Jiao Li
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China
- National Engineering Center for New Vaccine Research, Beijing, China
| | - Yi Bo Ma
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China
- National Engineering Center for New Vaccine Research, Beijing, China
| | - Ze Yuan Dong
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China
- National Engineering Center for New Vaccine Research, Beijing, China
| | - Jun Wei Hou
- National Engineering Center for New Vaccine Research, Beijing, China
- The Sixth Laboratory, National Vaccine and Serum Institute (NVSI), Beijing, China
| | - Fu Jie Shen
- National Engineering Center for New Vaccine Research, Beijing, China
- The Sixth Laboratory, National Vaccine and Serum Institute (NVSI), Beijing, China
| | - Wei Bu Wang
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China
- National Engineering Center for New Vaccine Research, Beijing, China
| | - Qi Ming Li
- National Engineering Center for New Vaccine Research, Beijing, China.
- The Sixth Laboratory, National Vaccine and Serum Institute (NVSI), Beijing, China.
| | - Ji Guo Su
- High Performance Computing Center, National Vaccine and Serum Institute (NVSI), Beijing, China.
- National Engineering Center for New Vaccine Research, Beijing, China.
| |
Collapse
|
3
|
Badonyi M, Marsh JA. Proteome-scale prediction of molecular mechanisms underlying dominant genetic diseases. PLoS One 2024; 19:e0307312. [PMID: 39172982 PMCID: PMC11341024 DOI: 10.1371/journal.pone.0307312] [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/07/2024] [Accepted: 06/26/2024] [Indexed: 08/24/2024] Open
Abstract
Many dominant genetic disorders result from protein-altering mutations, acting primarily through dominant-negative (DN), gain-of-function (GOF), and loss-of-function (LOF) mechanisms. Deciphering the mechanisms by which dominant diseases exert their effects is often experimentally challenging and resource intensive, but is essential for developing appropriate therapeutic approaches. Diseases that arise via a LOF mechanism are more amenable to be treated by conventional gene therapy, whereas DN and GOF mechanisms may require gene editing or targeting by small molecules. Moreover, pathogenic missense mutations that act via DN and GOF mechanisms are more difficult to identify than those that act via LOF using nearly all currently available variant effect predictors. Here, we introduce a tripartite statistical model made up of support vector machine binary classifiers trained to predict whether human protein coding genes are likely to be associated with DN, GOF, or LOF molecular disease mechanisms. We test the utility of the predictions by examining biologically and clinically meaningful properties known to be associated with the mechanisms. Our results strongly support that the models are able to generalise on unseen data and offer insight into the functional attributes of proteins associated with different mechanisms. We hope that our predictions will serve as a springboard for researchers studying novel variants and those of uncertain clinical significance, guiding variant interpretation strategies and experimental characterisation. Predictions for the human UniProt reference proteome are available at https://osf.io/z4dcp/.
Collapse
Affiliation(s)
- Mihaly Badonyi
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, United Kingdom
| | - Joseph A. Marsh
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
4
|
Emanuelli G, Zhu J, Li W, Morrell NW, Marciniak SJ. Functional validation of EIF2AK4 (GCN2) missense variants associated with pulmonary arterial hypertension. Hum Mol Genet 2024; 33:1495-1505. [PMID: 38776952 PMCID: PMC11336063 DOI: 10.1093/hmg/ddae082] [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: 04/09/2024] [Accepted: 05/09/2024] [Indexed: 05/25/2024] Open
Abstract
Pulmonary arterial hypertension (PAH) is a disorder with a large genetic component. Biallelic mutations of EIF2AK4, which encodes the kinase GCN2, are causal in two ultra-rare subtypes of PAH, pulmonary veno-occlusive disease and pulmonary capillary haemangiomatosis. EIF2AK4 variants of unknown significance have also been identified in patients with classical PAH, though their relationship to disease remains unclear. To provide patients with diagnostic information and enable family testing, the functional consequences of such rare variants must be determined, but existing computational methods are imperfect. We applied a suite of bioinformatic and experimental approaches to sixteen EIF2AK4 variants that had been identified in patients. By experimentally testing the functional integrity of the integrated stress response (ISR) downstream of GCN2, we determined that existing computational tools have insufficient sensitivity to reliably predict impaired kinase function. We determined experimentally that several EIF2AK4 variants identified in patients with classical PAH had preserved function and are therefore likely to be non-pathogenic. The dysfunctional variants of GCN2 that we identified could be subclassified into three groups: misfolded, kinase-dead, and hypomorphic. Intriguingly, members of the hypomorphic group were amenable to paradoxical activation by a type-1½ GCN2 kinase inhibitor. This experiment approach may aid in the clinical stratification of EIF2AK4 variants and potentially identify hypomorophic alleles receptive to pharmacological activation.
Collapse
Affiliation(s)
- Giulia Emanuelli
- Cambridge Institute for Medical Research (CIMR), University of Cambridge, Keith Peters Building, Biomedical Campus, Hills Rd, Cambridge CB2 0XY, United Kingdom
| | - JiaYi Zhu
- Cambridge Institute for Medical Research (CIMR), University of Cambridge, Keith Peters Building, Biomedical Campus, Hills Rd, Cambridge CB2 0XY, United Kingdom
| | - Wei Li
- Victor Phillip Dahdaleh Heart and Lung Research Institute, University of Cambridge, Papworth Road, Trumpington, Cambridge CB2 0BB, United Kingdom
- Department of Medicine, University of Cambridge, Addenbrooke's Hospital (Box 157), Hills Road, Cambridge CB2 2QQ, United Kingdom
| | - Nicholas W Morrell
- Victor Phillip Dahdaleh Heart and Lung Research Institute, University of Cambridge, Papworth Road, Trumpington, Cambridge CB2 0BB, United Kingdom
- Department of Medicine, University of Cambridge, Addenbrooke's Hospital (Box 157), Hills Road, Cambridge CB2 2QQ, United Kingdom
- Royal Papworth Hospital NHS Foundation Trust, Papworth Rd, Trumpington, Cambridge CB2 0AY, United Kingdom
| | - Stefan J Marciniak
- Cambridge Institute for Medical Research (CIMR), University of Cambridge, Keith Peters Building, Biomedical Campus, Hills Rd, Cambridge CB2 0XY, United Kingdom
- Department of Medicine, University of Cambridge, Addenbrooke's Hospital (Box 157), Hills Road, Cambridge CB2 2QQ, United Kingdom
- Royal Papworth Hospital NHS Foundation Trust, Papworth Rd, Trumpington, Cambridge CB2 0AY, United Kingdom
| |
Collapse
|
5
|
Ozkan S, Padilla N, de la Cruz X. QAFI: a novel method for quantitative estimation of missense variant impact using protein-specific predictors and ensemble learning. Hum Genet 2024:10.1007/s00439-024-02692-z. [PMID: 39048855 DOI: 10.1007/s00439-024-02692-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Accepted: 07/14/2024] [Indexed: 07/27/2024]
Abstract
Next-generation sequencing (NGS) has revolutionized genetic diagnostics, yet its application in precision medicine remains incomplete, despite significant advances in computational tools for variant annotation. Many variants remain unannotated, and existing tools often fail to accurately predict the range of impacts that variants have on protein function. This limitation restricts their utility in relevant applications such as predicting disease severity and onset age. In response to these challenges, a new generation of computational models is emerging, aimed at producing quantitative predictions of genetic variant impacts. However, the field is still in its early stages, and several issues need to be addressed, including improved performance and better interpretability. This study introduces QAFI, a novel methodology that integrates protein-specific regression models within an ensemble learning framework, utilizing conservation-based and structure-related features derived from AlphaFold models. Our findings indicate that QAFI significantly enhances the accuracy of quantitative predictions across various proteins. The approach has been rigorously validated through its application in the CAGI6 contest, focusing on ARSA protein variants, and further tested on a comprehensive set of clinically labeled variants, demonstrating its generalizability and robust predictive power. The straightforward nature of our models may also contribute to better interpretability of the results.
Collapse
Affiliation(s)
- Selen Ozkan
- Research Unit in Clinical and Translational Bioinformatics, Vall d'Hebron Institute of Research (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Natàlia Padilla
- Research Unit in Clinical and Translational Bioinformatics, Vall d'Hebron Institute of Research (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Xavier de la Cruz
- Research Unit in Clinical and Translational Bioinformatics, Vall d'Hebron Institute of Research (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain.
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain.
| |
Collapse
|
6
|
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
|
7
|
Visani GM, Pun MN, Galvin W, Daniel E, Borisiak K, Wagura U, Nourmohammad A. HERMES: Holographic Equivariant neuRal network model for Mutational Effect and Stability prediction. ARXIV 2024:arXiv:2407.06703v1. [PMID: 39040640 PMCID: PMC11261993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
Predicting the stability and fitness effects of amino acid mutations in proteins is a cornerstone of biological discovery and engineering. Various experimental techniques have been developed to measure mutational effects, providing us with extensive datasets across a diverse range of proteins. By training on these data, traditional computational modeling and more recent machine learning approaches have advanced significantly in predicting mutational effects. Here, we introduce HERMES, a 3D rotationally equivariant structure-based neural network model for mutational effect and stability prediction. Pre-trained to predict amino acid propensity from its surrounding 3D structure, HERMES can be fine-tuned for mutational effects using our open-source code. We present a suite of HERMES models, pre-trained with different strategies, and fine-tuned to predict the stability effect of mutations. Benchmarking against other models shows that HERMES often outperforms or matches their performance in predicting mutational effect on stability, binding, and fitness. HERMES offers versatile tools for evaluating mutational effects and can be fine-tuned for specific predictive objectives.
Collapse
Affiliation(s)
- Gian Marco Visani
- Department of Computer Science and Engineering, University of Washington, Seattle, USA
| | - Michael N. Pun
- Department of Physics, University of Washington, 3910 15th Avenue Northeast, Seattle, WA 98195, USA
| | - William Galvin
- Department of Computer Science and Engineering, University of Washington, Seattle, USA
| | - Eric Daniel
- Department of Computer Science and Engineering, University of Washington, Seattle, USA
| | - Kevin Borisiak
- Department of Physics, University of Washington, 3910 15th Avenue Northeast, Seattle, WA 98195, USA
| | - Utheri Wagura
- Department of Physics, University of Washington, 3910 15th Avenue Northeast, Seattle, WA 98195, USA
- Department of Physics, Massachusetts Institute of Technology, 182 Memorial Dr, Cambridge, MA 02139
| | - Armita Nourmohammad
- Department of Computer Science and Engineering, University of Washington, Seattle, USA
- Department of Physics, University of Washington, 3910 15th Avenue Northeast, Seattle, WA 98195, USA
- Department of Applied Mathematics, University of Washington, Seattle, USA
- Fred Hutchinson cancer Research Center, 1100 Fairview ave N, Seattle, WA 98109, USA
| |
Collapse
|
8
|
McDonnell AF, Plech M, Livesey BJ, Gerasimavicius L, Owen LJ, Hall HN, FitzPatrick DR, Marsh JA, Kudla G. Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants. Mol Syst Biol 2024; 20:825-844. [PMID: 38849565 PMCID: PMC11219921 DOI: 10.1038/s44320-024-00043-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Revised: 04/05/2024] [Accepted: 05/14/2024] [Indexed: 06/09/2024] Open
Abstract
Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.
Collapse
Affiliation(s)
- Alexander F McDonnell
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Marcin Plech
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Benjamin J Livesey
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Lukas Gerasimavicius
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Liusaidh J Owen
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Hildegard Nikki Hall
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - David R FitzPatrick
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Grzegorz Kudla
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, EH4 2XU, UK.
| |
Collapse
|
9
|
Chu SKS, Narang K, Siegel JB. Protein stability prediction by fine-tuning a protein language model on a mega-scale dataset. PLoS Comput Biol 2024; 20:e1012248. [PMID: 39038042 PMCID: PMC11293664 DOI: 10.1371/journal.pcbi.1012248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 08/01/2024] [Accepted: 06/13/2024] [Indexed: 07/24/2024] Open
Abstract
Protein stability plays a crucial role in a variety of applications, such as food processing, therapeutics, and the identification of pathogenic mutations. Engineering campaigns commonly seek to improve protein stability, and there is a strong interest in streamlining these processes to enable rapid optimization of highly stabilized proteins with fewer iterations. In this work, we explore utilizing a mega-scale dataset to develop a protein language model optimized for stability prediction. ESMtherm is trained on the folding stability of 528k natural and de novo sequences derived from 461 protein domains and can accommodate deletions, insertions, and multiple-point mutations. We show that a protein language model can be fine-tuned to predict folding stability. ESMtherm performs reasonably on small protein domains and generalizes to sequences distal from the training set. Lastly, we discuss our model's limitations compared to other state-of-the-art methods in generalizing to larger protein scaffolds. Our results highlight the need for large-scale stability measurements on a diverse dataset that mirrors the distribution of sequence lengths commonly observed in nature.
Collapse
Affiliation(s)
- Simon K. S. Chu
- Biophysics Graduate Program, University of California Davis, Davis, California, United States of America
| | - Kush Narang
- College of Biological Sciences, University of California Davis, Davis, California, United States of America
| | - Justin B. Siegel
- Genome Center, University of California Davis, Davis, California, United States of America
- Department of Chemistry, University of California Davis, Davis, California, United States of America
- Department of Biochemistry and Molecular Medicine, University of California Davis, Davis, California, United States of America
| |
Collapse
|
10
|
Shorthouse D, Lister H, Freeman GS, Hall BA. Understanding large scale sequencing datasets through changes to protein folding. Brief Funct Genomics 2024:elae007. [PMID: 38521964 DOI: 10.1093/bfgp/elae007] [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: 10/08/2023] [Revised: 02/26/2024] [Accepted: 03/01/2024] [Indexed: 03/25/2024] Open
Abstract
The expansion of high-quality, low-cost sequencing has created an enormous opportunity to understand how genetic variants alter cellular behaviour in disease. The high diversity of mutations observed has however drawn a spotlight onto the need for predictive modelling of mutational effects on phenotype from variants of uncertain significance. This is particularly important in the clinic due to the potential value in guiding clinical diagnosis and patient treatment. Recent computational modelling has highlighted the importance of mutation induced protein misfolding as a common mechanism for loss of protein or domain function, aided by developments in methods that make large computational screens tractable. Here we review recent applications of this approach to different genes, and how they have enabled and supported subsequent studies. We further discuss developments in the approach and the role for the approach in light of increasingly high throughput experimental approaches.
Collapse
Affiliation(s)
- David Shorthouse
- School of Pharmacy, University College London, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Harris Lister
- Department of Medical Physics and Biomedical Engineering, Malet Place Engineering Building, University College London, Gower Street, London WC1E 6BT, UK
| | - Gemma S Freeman
- Department of Medical Physics and Biomedical Engineering, Malet Place Engineering Building, University College London, Gower Street, London WC1E 6BT, UK
| | - Benjamin A Hall
- Department of Medical Physics and Biomedical Engineering, Malet Place Engineering Building, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
11
|
Vincenzi M, Mercurio FA, Autiero I, Leone M. Cancer-Related Mutations in the Sam Domains of EphA2 Receptor and Ship2 Lipid Phosphatase: A Computational Study. Molecules 2024; 29:1024. [PMID: 38474536 DOI: 10.3390/molecules29051024] [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: 12/22/2023] [Revised: 02/09/2024] [Accepted: 02/22/2024] [Indexed: 03/14/2024] Open
Abstract
The lipid phosphatase Ship2 interacts with the EphA2 receptor by forming a heterotypic Sam (sterile alpha motif)-Sam complex. Ship2 works as a negative regulator of receptor endocytosis and consequent degradation, and anti-oncogenic effects in cancer cells should be induced by hindering its association with EphA2. Herein, a computational approach is presented to investigate the relationship between Ship2-Sam/EphA2-Sam interaction and cancer onset and further progression. A search was first conducted through the COSMIC (Catalogue of Somatic Mutations in Cancer) database to identify cancer-related missense mutations positioned inside or close to the EphA2-Sam and Ship2-Sam reciprocal binding interfaces. Next, potential differences in the chemical-physical properties of mutant and wild-type Sam domains were evaluated by bioinformatics tools based on analyses of primary sequences. Three-dimensional (3D) structural models of mutated EphA2-Sam and Ship2-Sam domains were built as well and deeply analysed with diverse computational instruments, including molecular dynamics, to classify potentially stabilizing and destabilizing mutations. In the end, the influence of mutations on the EphA2-Sam/Ship2-Sam interaction was studied through docking techniques. This in silico approach contributes to understanding, at the molecular level, the mutation/cancer relationship by predicting if amino acid substitutions could modulate EphA2 receptor endocytosis.
Collapse
Affiliation(s)
- Marian Vincenzi
- Institute of Biostructures and Bioimaging, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Flavia Anna Mercurio
- Institute of Biostructures and Bioimaging, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Ida Autiero
- Institute of Biostructures and Bioimaging, Via Pietro Castellino 111, 80131 Naples, Italy
| | - Marilisa Leone
- Institute of Biostructures and Bioimaging, Via Pietro Castellino 111, 80131 Naples, Italy
| |
Collapse
|
12
|
Clausen L, Voutsinos V, Cagiada M, Johansson KE, Grønbæk-Thygesen M, Nariya S, Powell RL, Have MKN, Oestergaard VH, Stein A, Fowler DM, Lindorff-Larsen K, Hartmann-Petersen R. A mutational atlas for Parkin proteostasis. Nat Commun 2024; 15:1541. [PMID: 38378758 PMCID: PMC10879094 DOI: 10.1038/s41467-024-45829-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Accepted: 02/01/2024] [Indexed: 02/22/2024] Open
Abstract
Proteostasis can be disturbed by mutations affecting folding and stability of the encoded protein. An example is the ubiquitin ligase Parkin, where gene variants result in autosomal recessive Parkinsonism. To uncover the pathological mechanism and provide comprehensive genotype-phenotype information, variant abundance by massively parallel sequencing (VAMP-seq) is leveraged to quantify the abundance of Parkin variants in cultured human cells. The resulting mutational map, covering 9219 out of the 9300 possible single-site amino acid substitutions and nonsense Parkin variants, shows that most low abundance variants are proteasome targets and are located within the structured domains of the protein. Half of the known disease-linked variants are found at low abundance. Systematic mapping of degradation signals (degrons) reveals an exposed degron region proximal to the so-called "activation element". This work provides examples of how missense variants may cause degradation either via destabilization of the native protein, or by introducing local signals for degradation.
Collapse
Affiliation(s)
- Lene Clausen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Vasileios Voutsinos
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Matteo Cagiada
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kristoffer E Johansson
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Martin Grønbæk-Thygesen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Snehal Nariya
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Rachel L Powell
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Magnus K N Have
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Amelie Stein
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Douglas M Fowler
- Department of Genome Sciences, University of Washington, Seattle, WA, USA.
- Department of Bioengineering, University of Washington, Seattle, WA, USA.
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Rasmus Hartmann-Petersen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
13
|
David C, Badonyi M, Kechiche R, Insalaco A, Zecca M, De Benedetti F, Orcesi S, Chiapparini L, Comoli P, Federici S, Gattorno M, Ginevrino M, Giorgio E, Matteo V, Moran-Alvarez P, Politano D, Prencipe G, Sirchia F, Volpi S, Masson C, Rice GI, Frémond ML, Lepelley A, Marsh JA, Crow YJ. Interface Gain-of-Function Mutations in TLR7 Cause Systemic and Neuro-inflammatory Disease. J Clin Immunol 2024; 44:60. [PMID: 38324161 PMCID: PMC10850255 DOI: 10.1007/s10875-024-01660-6] [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: 12/06/2023] [Accepted: 01/17/2024] [Indexed: 02/08/2024]
Abstract
TLR7 recognizes pathogen-derived single-stranded RNA (ssRNA), a function integral to the innate immune response to viral infection. Notably, TLR7 can also recognize self-derived ssRNA, with gain-of-function mutations in human TLR7 recently identified to cause both early-onset systemic lupus erythematosus (SLE) and neuromyelitis optica. Here, we describe two novel mutations in TLR7, F507S and L528I. While the L528I substitution arose de novo, the F507S mutation was present in three individuals from the same family, including a severely affected male, notably given that the TLR7 gene is situated on the X chromosome and that all other cases so far described have been female. The observation of mutations at residues 507 and 528 of TLR7 indicates the importance of the TLR7 dimerization interface in maintaining immune homeostasis, where we predict that altered homo-dimerization enhances TLR7 signaling. Finally, while mutations in TLR7 can result in SLE-like disease, our data suggest a broader phenotypic spectrum associated with TLR7 gain-of-function, including significant neurological involvement.
Collapse
Affiliation(s)
- Clémence David
- Laboratory of Neurogenetics and NeuroinflammationImagine Institute, INSERM UMR1163, Paris, France
| | - Mihaly Badonyi
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, UK
| | - Robin Kechiche
- Laboratory of Neurogenetics and NeuroinflammationImagine Institute, INSERM UMR1163, Paris, France
- Department of Paediatric Hematology-Immunology and Rheumatology, Necker-Enfants Malades Hospital, AP-HP, Paris, France
| | - Antonella Insalaco
- Division of Rheumatology, IRCCS Ospedale Pediatrico Bambino Gesù, Rome, Italy
| | - Marco Zecca
- Pediatric Haematology/Oncology, Fondazione IRCCS Policlinico San Matteo, Pavia, Italy
| | | | - Simona Orcesi
- Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
- Child Neurology and Psychiatry Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Luisa Chiapparini
- Neuroradiology Unit, Fondazione IRCCS Policlinico San Matteo, Pavia, Italy
| | - Patrizia Comoli
- Cell Factory, Fondazione IRCCS Policlinico San Matteo, Pavia, Italy
| | - Silvia Federici
- Division of Rheumatology, IRCCS Ospedale Pediatrico Bambino Gesù, Rome, Italy
| | - Marco Gattorno
- UOC Reumatologia E Malattie Autoinfiammatorie, IRCCS Istituto Giannina Gaslini, Genoa, Italy
| | - Monia Ginevrino
- Translational Cytogenomics Research Unit, Bambino Gesù Children's Hospital, IRCCS, Rome, Italy
| | - Elisa Giorgio
- Department of Molecular Medicine, University of Pavia, Pavia, Italy
- Medical Genetics Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Valentina Matteo
- Laboratory of Immuno-Rheumatology, Bambino Gesù Children's Hospital, IRCCS, Rome, Italy
| | | | - Davide Politano
- Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
- Child Neurology and Psychiatry Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Giusi Prencipe
- Laboratory of Immuno-Rheumatology, Bambino Gesù Children's Hospital, IRCCS, Rome, Italy
| | - Fabio Sirchia
- Department of Molecular Medicine, University of Pavia, Pavia, Italy
- Medical Genetics Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Stefano Volpi
- UOC Reumatologia E Malattie Autoinfiammatorie, IRCCS Istituto Giannina Gaslini, Genoa, Italy
- Dipartimento Di Neuroscienze, Riabilitazione, Oftalmologia, Genetica e Scienze Materno-Infantili (DINOGMI), Università Degli Studi Di Genova, Genoa, Italy
| | - Cécile Masson
- Bioinformatics Core Facility, Paris-Cité University-Structure Fédérative de Recherche Necker, INSERM US24/CNRS UMS3633, Paris, France
| | - Gillian I Rice
- Division of Evolution and Genomic Sciences, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester Academic Health Science Centre, Manchester, UK
| | - Marie-Louise Frémond
- Laboratory of Neurogenetics and NeuroinflammationImagine Institute, INSERM UMR1163, Paris, France
- Department of Paediatric Hematology-Immunology and Rheumatology, Necker-Enfants Malades Hospital, AP-HP, Paris, France
- Reference Center for Rheumatic, AutoImmune and Systemic Diseases in Children (RAISE), Paris, France
| | - Alice Lepelley
- Laboratory of Neurogenetics and NeuroinflammationImagine Institute, INSERM UMR1163, Paris, France
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, UK
| | - Yanick J Crow
- Laboratory of Neurogenetics and NeuroinflammationImagine Institute, INSERM UMR1163, Paris, France.
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, UK.
- University Paris Cité, Paris, France.
| |
Collapse
|
14
|
Pandey P, Alexov E. Most Monogenic Disorders Are Caused by Mutations Altering Protein Folding Free Energy. Int J Mol Sci 2024; 25:1963. [PMID: 38396641 PMCID: PMC10888012 DOI: 10.3390/ijms25041963] [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: 12/29/2023] [Revised: 01/31/2024] [Accepted: 02/02/2024] [Indexed: 02/25/2024] Open
Abstract
Revealing the molecular effect that pathogenic missense mutations have on the corresponding protein is crucial for developing therapeutic solutions. This is especially important for monogenic diseases since, for most of them, there is no treatment available, while typically, the treatment should be provided in the early development stages. This requires fast targeted drug development at a low cost. Here, we report an updated database of monogenic disorders (MOGEDO), which includes 768 proteins and the corresponding 2559 pathogenic and 1763 benign mutations, along with the functional classification of the corresponding proteins. Using the database and various computational tools that predict folding free energy change (ΔΔG), we demonstrate that, on average, 70% of pathogenic cases result in decreased protein stability. Such a large fraction indicates that one should aim at in silico screening for small molecules stabilizing the structure of the mutant protein. We emphasize that knowledge of ΔΔG is essential because one wants to develop stabilizers that compensate for ΔΔG, but do not make protein over-stable, since over-stable protein may be dysfunctional. We demonstrate that, by using ΔΔG and predicted solvent exposure of the mutation site, one can develop a predictive method that distinguishes pathogenic from benign mutations with a success rate even better than some of the leading pathogenicity predictors. Furthermore, hydrophobic-hydrophobic mutations have stronger correlations between folding free energy change and pathogenicity compared with others. Also, mutations involving Cys, Gly, Arg, Trp, and Tyr amino acids being replaced by any other amino acid are more likely to be pathogenic. To facilitate further detection of pathogenic mutations, the wild type of amino acids in the 768 proteins mentioned above was mutated to other 19 residues (14,847,817 mutations), the ΔΔG was calculated with SAAFEC-SEQ, and 5,506,051 mutations were predicted to be pathogenic.
Collapse
Affiliation(s)
| | - Emil Alexov
- Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA;
| |
Collapse
|
15
|
Utsuno Y, Hamada K, Hamanaka K, Miyoshi K, Tsuchimoto K, Sunada S, Itai T, Sakamoto M, Tsuchida N, Uchiyama Y, Koshimizu E, Fujita A, Miyatake S, Misawa K, Mizuguchi T, Kato Y, Saito K, Ogata K, Matsumoto N. Novel missense variants cause intermediate phenotypes in the phenotypic spectrum of SLC5A6-related disorders. J Hum Genet 2024; 69:69-77. [PMID: 38012394 DOI: 10.1038/s10038-023-01206-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2023] [Revised: 11/05/2023] [Accepted: 11/05/2023] [Indexed: 11/29/2023]
Abstract
SLC5A6 encodes the sodium-dependent multivitamin transporter, a transmembrane protein that uptakes biotin, pantothenic acid, and lipoic acid. Biallelic SLC5A6 variants cause sodium-dependent multivitamin transporter deficiency (SMVTD) and childhood-onset biotin-responsive peripheral motor neuropathy (COMNB), which both respond well to replacement therapy with the above three nutrients. SMVTD usually presents with various symptoms in multiple organs, such as gastrointestinal hemorrhage, brain atrophy, and global developmental delay, at birth or in infancy. Without nutrient replacement therapy, SMVTD can be lethal in early childhood. COMNB is clinically milder and has a later onset than SMVTD, at approximately 10 years of age. COMNB symptoms are mostly limited to peripheral motor neuropathy. Here we report three patients from one Japanese family harboring novel compound heterozygous missense variants in SLC5A6, namely NM_021095.4:c.[221C>T];[642G>C] p.[(Ser74Phe)];[(Gln214His)]. Both variants were predicted to be deleterious through multiple lines of evidence, including amino acid conservation, in silico predictions of pathogenicity, and protein structure considerations. Drosophila analysis also showed c.221C>T to be pathogenic. All three patients had congenital brain cysts on neonatal cranial imaging, but no other morphological abnormalities. They also had a mild motor developmental delay that almost completely resolved despite no treatment. In terms of severity, their phenotypes were intermediate between SMVTD and COMNB. From these findings we propose a new SLC5A6-related disorder, spontaneously remitting developmental delay with brain cysts (SRDDBC) whose phenotypic severity is between that of SMVTD and COMNB. Further clinical and genetic evidence is needed to support our suggestion.
Collapse
Affiliation(s)
- Yasuhiro Utsuno
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- Department of Obstetrics and Gynecology, Asahikawa Medical University, Hokkaido, Japan
| | - Keisuke Hamada
- Department of Biochemistry, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Kohei Hamanaka
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Keita Miyoshi
- Department of Chromosome Science, National Institute of Genetics, Research Organization of Information and Systems (ROIS), Shizuoka, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Shizuoka, Japan
| | - Keiji Tsuchimoto
- Department of Pediatrics, Kochi Health Sciences Center, Kochi, Japan
| | - Satoshi Sunada
- Department of Pediatrics, Kurashiki Central Hospital, Okayama, Japan
| | - Toshiyuki Itai
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Masamune Sakamoto
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- Department of Rare Disease Genomics, Yokohama City University Hospital, Yokohama, Japan
| | - Naomi Tsuchida
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- Department of Rare Disease Genomics, Yokohama City University Hospital, Yokohama, Japan
| | - Yuri Uchiyama
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- Department of Rare Disease Genomics, Yokohama City University Hospital, Yokohama, Japan
| | - Eriko Koshimizu
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Atsushi Fujita
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Satoko Miyatake
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- Department of Clinical Genetics, Yokohama City University Hospital, Yokohama, Japan
| | - Kazuharu Misawa
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan
| | - Takeshi Mizuguchi
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Yasuhito Kato
- Department of Obstetrics and Gynecology, Asahikawa Medical University, Hokkaido, Japan
| | - Kuniaki Saito
- Department of Chromosome Science, National Institute of Genetics, Research Organization of Information and Systems (ROIS), Shizuoka, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Shizuoka, Japan
| | - Kazuhiro Ogata
- Department of Biochemistry, Yokohama City University Graduate School of Medicine, Yokohama, Japan
| | - Naomichi Matsumoto
- Department of Human Genetics, Yokohama City University Graduate School of Medicine, Yokohama, Japan.
| |
Collapse
|
16
|
Sinclair M, Stein RA, Sheehan JH, Hawes EM, O’Brien RM, Tajkhorshid E, Claxton DP. Integrative analysis of pathogenic variants in glucose-6-phosphatase based on an AlphaFold2 model. PNAS NEXUS 2024; 3:pgae036. [PMID: 38328777 PMCID: PMC10849595 DOI: 10.1093/pnasnexus/pgae036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 01/09/2024] [Indexed: 02/09/2024]
Abstract
Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein glucose-6-phosphate catalytic subunit 1 (G6PC1) regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 causes glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. In this study, we determine the atomic interactions governing G6P binding as well as explore the perturbations imposed by disease-linked missense variants by subjecting an AlphaFold2 G6PC1 structural model to molecular dynamics simulations and in silico predictions of thermodynamic stability validated with robust in vitro and in situ biochemical assays. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. The introduction of GSD type 1a mutations modified the thermodynamic landscape, altered side chain packing and substrate-binding interactions, and induced trapping of catalytic intermediates. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm the active-site structural organization but also identify previously unobserved mechanistic contributions of catalytic and noncatalytic side chains.
Collapse
Affiliation(s)
- Matt Sinclair
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Richard A Stein
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Center for Applied Artificial Intelligence in Protein Dynamics, Vanderbilt University, Nashville, TN 37240, USA
| | - Jonathan H Sheehan
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37240, USA
- Division of Infectious Diseases, Department of Internal Medicine, Washington University School of Medicine, St Louis, MO 63110, USA
| | - Emily M Hawes
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
| | - Richard M O’Brien
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Derek P Claxton
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Center for Applied Artificial Intelligence in Protein Dynamics, Vanderbilt University, Nashville, TN 37240, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37240, USA
| |
Collapse
|
17
|
Pucci F, Zerihun MB, Rooman M, Schug A. pycofitness-Evaluating the fitness landscape of RNA and protein sequences. Bioinformatics 2024; 40:btae074. [PMID: 38335928 PMCID: PMC10881095 DOI: 10.1093/bioinformatics/btae074] [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: 05/29/2023] [Revised: 01/25/2024] [Accepted: 02/06/2024] [Indexed: 02/12/2024] Open
Abstract
MOTIVATION The accurate prediction of how mutations change biophysical properties of proteins or RNA is a major goal in computational biology with tremendous impacts on protein design and genetic variant interpretation. Evolutionary approaches such as coevolution can help solving this issue. RESULTS We present pycofitness, a standalone Python-based software package for the in silico mutagenesis of protein and RNA sequences. It is based on coevolution and, more specifically, on a popular inverse statistical approach, namely direct coupling analysis by pseudo-likelihood maximization. Its efficient implementation and user-friendly command line interface make it an easy-to-use tool even for researchers with no bioinformatics background. To illustrate its strengths, we present three applications in which pycofitness efficiently predicts the deleteriousness of genetic variants and the effect of mutations on protein fitness and thermodynamic stability. AVAILABILITY AND IMPLEMENTATION https://github.com/KIT-MBS/pycofitness.
Collapse
Affiliation(s)
- Fabrizio Pucci
- Computational Biology and Bioinformatics, Université Libre de Bruxelles, 1050 Brussels, Belgium
- Interuniversity Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
| | - Mehari B Zerihun
- John von Neumann Institute for Computing, Jülich Supercomputer Centre, 52428 Jülich, Germany
| | - Marianne Rooman
- Computational Biology and Bioinformatics, Université Libre de Bruxelles, 1050 Brussels, Belgium
- Interuniversity Institute of Bioinformatics in Brussels, 1050 Brussels, Belgium
| | - Alexander Schug
- John von Neumann Institute for Computing, Jülich Supercomputer Centre, 52428 Jülich, Germany
- Department of Biology, University of Duisburg-Essen, D-45141 Essen, Germany
| |
Collapse
|
18
|
Nourbakhsh M, Degn K, Saksager A, Tiberti M, Papaleo E. Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks. Brief Bioinform 2024; 25:bbad519. [PMID: 38261338 PMCID: PMC10805075 DOI: 10.1093/bib/bbad519] [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: 06/09/2023] [Revised: 11/27/2023] [Accepted: 12/11/2023] [Indexed: 01/24/2024] Open
Abstract
The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.
Collapse
Affiliation(s)
- Mona Nourbakhsh
- Cancer Systems Biology, Section for Bioinformatics, Department of Health Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Kristine Degn
- Cancer Systems Biology, Section for Bioinformatics, Department of Health Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Astrid Saksager
- Cancer Systems Biology, Section for Bioinformatics, Department of Health Technology, Technical University of Denmark, 2800 Lyngby, Denmark
| | - Matteo Tiberti
- Cancer Structural Biology, Danish Cancer Institute, 2100 Copenhagen, Denmark
| | - Elena Papaleo
- Cancer Systems Biology, Section for Bioinformatics, Department of Health Technology, Technical University of Denmark, 2800 Lyngby, Denmark
- Cancer Structural Biology, Danish Cancer Institute, 2100 Copenhagen, Denmark
| |
Collapse
|
19
|
Boujemaa M, Nouira F, Jandoubi N, Mejri N, Bouaziz H, Charfeddine C, Ben Nasr S, Labidi S, El Benna H, Berrazega Y, Rachdi H, Daoud N, Benna F, Haddaoui A, Abdelhak S, Samir Boubaker M, Boussen H, Hamdi Y. Uncovering the clinical relevance of unclassified variants in DNA repair genes: a focus on BRCA negative Tunisian cancer families. Front Genet 2024; 15:1327894. [PMID: 38313678 PMCID: PMC10834681 DOI: 10.3389/fgene.2024.1327894] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 01/02/2024] [Indexed: 02/06/2024] Open
Abstract
Introduction: Recent advances in sequencing technologies have significantly increased our capability to acquire large amounts of genetic data. However, the clinical relevance of the generated data continues to be challenging particularly with the identification of Variants of Uncertain Significance (VUSs) whose pathogenicity remains unclear. In the current report, we aim to evaluate the clinical relevance and the pathogenicity of VUSs in DNA repair genes among Tunisian breast cancer families. Methods: A total of 67 unsolved breast cancer cases have been investigated. The pathogenicity of VUSs identified within 26 DNA repair genes was assessed using different in silico prediction tools including SIFT, PolyPhen2, Align-GVGD and VarSEAK. Effects on the 3D structure were evaluated using the stability predictor DynaMut and molecular dynamics simulation with NAMD. Family segregation analysis was also performed. Results: Among a total of 37 VUSs identified, 11 variants are likely deleterious affecting ATM, BLM, CHEK2, ERCC3, FANCC, FANCG, MSH2, PMS2 and RAD50 genes. The BLM variant, c.3254dupT, is novel and seems to be associated with increased risk of breast, endometrial and colon cancer. Moreover, c.6115G>A in ATM and c.592+3A>T in CHEK2 were of keen interest identified in families with multiple breast cancer cases and their familial cosegregation with disease has been also confirmed. In addition, functional in silico analyses revealed that the ATM variant may lead to protein immobilization and rigidification thus decreasing its activity. We have also shown that FANCC and FANCG variants may lead to protein destabilization and alteration of the structure compactness which may affect FANCC and FANCG protein activity. Conclusion: Our findings revealed that VUSs in DNA repair genes might be associated with increased cancer risk and highlight the need for variant reclassification for better disease management. This will help to improve the genetic diagnosis and therapeutic strategies of cancer patients not only in Tunisia but also in neighboring countries.
Collapse
Affiliation(s)
- Maroua Boujemaa
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
| | - Fatma Nouira
- Laboratory of Bioactive Substances, Center of Biotechnology of Borj Cedria, University of Tunis El Manar, Hamam, Tunisia
| | - Nouha Jandoubi
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
| | - Nesrine Mejri
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Hanen Bouaziz
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Surgical Oncology Department, Salah Azaiez Institute of Cancer, Tunis, Tunisia
| | - Cherine Charfeddine
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- High Institute of Biotechnology of Sidi Thabet, Biotechpole of Sidi Thabet, University of Manouba, Ariana, Tunisia
| | - Sonia Ben Nasr
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Department of Medical Oncology, Military Hospital of Tunis, Tunis, Tunisia
| | - Soumaya Labidi
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Houda El Benna
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Yosra Berrazega
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Haifa Rachdi
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Nouha Daoud
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Farouk Benna
- Radiation Oncology Department, Salah Azaiez Institute, Tunis, Tunisia
| | | | - Sonia Abdelhak
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
| | - Mohamed Samir Boubaker
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Laboratory of Human and Experimental Pathology, Institut Pasteur de Tunis, Tunis, Tunisia
| | - Hamouda Boussen
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Medical Oncology Department, Abderrahman Mami Hospital, Faculty of Medicine Tunis, University Tunis El Manar, Tunis, Tunisia
| | - Yosr Hamdi
- Laboratory of Biomedical Genomics and Oncogenetics, LR20IPT05, Institut Pasteur de Tunis, University of Tunis El Manar, Tunis, Tunisia
- Laboratory of Human and Experimental Pathology, Institut Pasteur de Tunis, Tunis, Tunisia
| |
Collapse
|
20
|
Díaz-Rodríguez SM, Herrero-Turrión MJ, García-Peral C, Gómez-Nieto R. Delving into the significance of the His289Tyr single-nucleotide polymorphism in the glutamate ionotropic receptor kainate-1 ( Grik1) gene of a genetically audiogenic seizure model. Front Mol Neurosci 2024; 16:1322750. [PMID: 38249292 PMCID: PMC10797026 DOI: 10.3389/fnmol.2023.1322750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Accepted: 12/07/2023] [Indexed: 01/23/2024] Open
Abstract
Genetic abnormalities affecting glutamate receptors are central to excitatory overload-driven neuronal mechanisms that culminate in seizures, making them pivotal targets in epilepsy research. Increasingly used to advance this field, the genetically audiogenic seizure hamster from Salamanca (GASH/Sal) exhibits generalized seizures triggered by high-intensity acoustic stimulation and harbors significant genetic variants recently identified through whole-exome sequencing. Here, we addressed the influence of the missense single-nucleotide polymorphism (C9586732T, p.His289Tyr) in the glutamate receptor ionotropic kainate-1 (Grik1) gene and its implications for the GASH/Sal seizure susceptibility. Using a protein 3D structure prediction, we showed a potential effect of this sequence variation, located in the amino-terminal domain, on the stability and/or conformation of the kainate receptor subunit-1 protein (GluK1). We further employed a multi-technique approach, encompassing gene expression analysis (RT-qPCR), Western blotting, and immunohistochemistry in bright-field and confocal fluorescence microscopy, to investigate critical seizure-associated brain regions in GASH/Sal animals under seizure-free conditions compared to matched wild-type controls. We detected disruptions in the transcriptional profile of the Grik1 gene within the audiogenic seizure-associated neuronal network. Alterations in GluK1 protein levels were also observed in various brain structures, accompanied by an unexpected lower molecular weight band in the inferior and superior colliculi. This correlated with substantial disparities in GluK1-immunolabeling distribution across multiple brain regions, including the cerebellum, hippocampus, subdivisions of the inferior and superior colliculi, and the prefrontal cortex. Notably, the diffuse immunolabeling accumulated within perikarya, axonal fibers and terminals, exhibiting a prominent concentration in proximity to the cell nucleus. This suggests potential disturbances in the GluK1-trafficking mechanism, which could subsequently affect glutamate synaptic transmission. Overall, our study sheds light on the genetic underpinnings of seizures and underscores the importance of investigating the molecular mechanisms behind synaptic dysfunction in epileptic neural networks, laying a crucial foundation for future research and therapeutic strategies targeting GluK1-containing kainate receptors.
Collapse
Affiliation(s)
- Sandra M. Díaz-Rodríguez
- Institute of Neuroscience of Castilla y León (INCYL), University of Salamanca, Salamanca, Spain
- Institute for Biomedical Research of Salamanca (IBSAL), Salamanca, Spain
- Department of Cell Biology and Pathology, Faculty of Medicine, University of Salamanca, Salamanca, Spain
| | - M. Javier Herrero-Turrión
- Institute of Neuroscience of Castilla y León (INCYL), University of Salamanca, Salamanca, Spain
- Institute for Biomedical Research of Salamanca (IBSAL), Salamanca, Spain
- Neurological Tissue Bank INCYL (BTN-INCYL), Salamanca, Spain
| | - Carlos García-Peral
- Institute of Neuroscience of Castilla y León (INCYL), University of Salamanca, Salamanca, Spain
- Institute for Biomedical Research of Salamanca (IBSAL), Salamanca, Spain
| | - Ricardo Gómez-Nieto
- Institute of Neuroscience of Castilla y León (INCYL), University of Salamanca, Salamanca, Spain
- Institute for Biomedical Research of Salamanca (IBSAL), Salamanca, Spain
- Department of Cell Biology and Pathology, Faculty of Medicine, University of Salamanca, Salamanca, Spain
| |
Collapse
|
21
|
James JK, Norland K, Johar AS, Kullo IJ. Deep generative models of LDLR protein structure to predict variant pathogenicity. J Lipid Res 2023; 64:100455. [PMID: 37821076 PMCID: PMC10696256 DOI: 10.1016/j.jlr.2023.100455] [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: 03/22/2023] [Revised: 09/16/2023] [Accepted: 10/05/2023] [Indexed: 10/13/2023] Open
Abstract
The complex structure and function of low density lipoprotein receptor (LDLR) makes classification of protein-coding missense variants challenging. Deep generative models, including Evolutionary model of Variant Effect (EVE), Evolutionary Scale Modeling (ESM), and AlphaFold 2 (AF2), have enabled significant progress in the prediction of protein structure and function. ESM and EVE directly estimate the likelihood of a variant sequence but are purely data-driven and challenging to interpret. AF2 predicts LDLR structures, but variant effects are explicitly modeled by estimating changes in stability. We tested the effectiveness of these models for predicting variant pathogenicity compared to established methods. AF2 produced two distinct conformations based on a novel hinge mechanism. Within ESM's hidden space, benign and pathogenic variants had different distributions. In EVE, these distributions were similar. EVE and ESM were comparable to Polyphen-2, SIFT, REVEL, and Primate AI for predicting binary classifications in ClinVar. However, they were more strongly correlated with experimental measures of LDL uptake. AF2 poorly performed in these tasks. Using the UK Biobank to compare association with clinical phenotypes, ESM and EVE were more strongly associated with serum LDL-C than Polyphen-2. ESM was able to identify variants with more extreme LDL-C levels than EVE and had a significantly stronger association with atherosclerotic cardiovascular disease. In conclusion, AF2 predicted LDLR structures do not accurately model variant pathogenicity. ESM and EVE are competitive with prior scoring methods for prediction based on binary classifications in ClinVar but are superior based on correlations with experimental assays and clinical phenotypes.
Collapse
Affiliation(s)
- Jose K James
- Department of Cardiovascular Medicine, Mayo Clinic, Rochester, MN, USA
| | - Kristjan Norland
- Department of Cardiovascular Medicine, Mayo Clinic, Rochester, MN, USA
| | - Angad S Johar
- Department of Cardiovascular Medicine, Mayo Clinic, Rochester, MN, USA
| | - Iftikhar J Kullo
- Department of Cardiovascular Medicine, Mayo Clinic, Rochester, MN, USA; Gonda Vascular Center, Mayo Clinic, Rochester, MN, USA.
| |
Collapse
|
22
|
Pandey P, Alexov E. Most monogenic disorders are caused by mutations altering protein folding free energy. RESEARCH SQUARE 2023:rs.3.rs-3442589. [PMID: 37886551 PMCID: PMC10602188 DOI: 10.21203/rs.3.rs-3442589/v1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
Revealing the molecular effect that pathogenic missense mutations cause on the corresponding protein is crucial for developing therapeutic solutions. This is especially important for monogenic diseases since, for most of them, there is no treatment available, while typically, the treatment should be provided in the early development stages. This requires fast, targeted drug development at a low cost. Here, we report a database of monogenic disorders (MOGEDO), which includes 768 proteins, the corresponding 2559 pathogenic and 1763 benign mutations, along with the functional classification of the corresponding proteins. Using the database and various computational tools that predict folding free energy change (ΔΔG), we demonstrate that, on average, 70% of pathogenic cases result in decreased protein stability. Such a large fraction indicates that one should aim at in-silico screening for small molecules stabilizing the structure of the mutant protein. We emphasize that knowledge of ΔΔG is essential because one wants to develop stabilizers that compensate for ΔΔG but not to make protein over-stable since over-stable protein may be dysfunctional. We demonstrate that using ΔΔG and predicted solvent exposure of the mutation site; one can develop a predictive method that distinguishes pathogenic from benign mutation with a success rate even better than some of the leading pathogenicity predictors. Furthermore, hydrophobic-hydrophobic mutations have stronger correlations between folding free energy change and pathogenicity compared with others. Also, mutations involving Cys, Gly, Arg, Trp and Tyr amino acids being replaced by any other amino acid are more likely to be pathogenic. To facilitate further detection of pathogenic mutations, the wild type of amino acids in the 768 proteins mentioned above was mutated to other 19 residues (14,847,817 mutations), and the ΔΔG was calculated with SAAFEC-SEQ, and 5,506,051 mutations were predicted to be pathogenic.
Collapse
|
23
|
Avsar O. Identification of the effects of pathogenic genetic variations of human CYP2C9 and CYP2D6: an in silico approach. NUCLEOSIDES, NUCLEOTIDES & NUCLEIC ACIDS 2023; 43:356-376. [PMID: 37747773 DOI: 10.1080/15257770.2023.2262519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 09/19/2023] [Indexed: 09/26/2023]
Abstract
Genetic variations in the human cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily D member 6 (CYP2D6) genes may affect drug metabolism and lead to alterations in phenotypes. Genetic variations are associated with toxicity, adverse drug reactions, inefficient treatment. Various in silico tools were combined to investigate the deleterious effects of missense non-synonymous single nucleotide polymorphisms (nsSNPs) of the human CYP2C9 and CYP2D6. The structural and functional effects of the high-risk non-synonymous SNPs in the human CYP2C9 and CYP2D6 were predicted by numerous computational mutation analysis methods. Out of 24 pathogenic missense SNPs in the CYP2C9, 22 nsSNPs had a decreasing effect on protein stability and 13 SNPs were showed to be located at conserved positions. Out of 27 high-risk deleterious non-synonymous SNPs in the human CYP2D6, 21 SNPs decreased protein stability and 16 nsSNPs were predicted to be positioned at conserved regions. Our present study suggests that the identified functional SNPs may affect drug metabolism associated with CYP2C9 and CYP2D6 enzymes.
Collapse
Affiliation(s)
- Orcun Avsar
- Department of Molecular Biology and Genetics, Faculty of Arts and Sciences, Hitit University, Corum, Türkiye
| |
Collapse
|
24
|
Sinclair M, Stein RA, Sheehan JH, Hawes EM, O'Brien RM, Tajkhorshid E, Claxton DP. Molecular mechanisms of catalytic inhibition for active site mutations in glucose-6-phosphatase catalytic subunit 1 linked to glycogen storage disease. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.13.532485. [PMID: 36993754 PMCID: PMC10054992 DOI: 10.1101/2023.03.13.532485] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein G6PC1 regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 cause glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. Exploiting a computational model of G6PC1 derived from the groundbreaking structure prediction algorithm AlphaFold2 (AF2), we combine molecular dynamics (MD) simulations and computational predictions of thermodynamic stability with a robust in vitro screening platform to define the atomic interactions governing G6P binding as well as explore the energetic perturbations imposed by disease-linked variants. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. Introduction of GSD type 1a mutations into the G6PC1 sequence elicits changes in G6P binding energy, thermostability and structural properties, suggesting multiple pathways of catalytic impairment. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm active site structural organization but also suggest novel mechanistic contributions of catalytic and non-catalytic side chains.
Collapse
|
25
|
Ferreira-Martins AJ, Castaldoni R, Alencar BM, Ferreira MV, Nogueira T, Rios RA, Lopes TJS. Full-scale network analysis reveals properties of the FV protein structure organization. Sci Rep 2023; 13:9546. [PMID: 37308572 DOI: 10.1038/s41598-023-36528-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 06/05/2023] [Indexed: 06/14/2023] Open
Abstract
Blood coagulation is a vital process for humans and other species. Following an injury to a blood vessel, a cascade of molecular signals is transmitted, inhibiting and activating more than a dozen coagulation factors and resulting in the formation of a fibrin clot that ceases the bleeding. In this process, the Coagulation factor V (FV) is a master regulator, coordinating critical steps of this process. Mutations to this factor result in spontaneous bleeding episodes and prolonged hemorrhage after trauma or surgery. Although the role of FV is well characterized, it is unclear how single-point mutations affect its structure. In this study, to understand the effect of mutations, we created a detailed network map of this protein, where each node is a residue, and two residues are connected if they are in close proximity in the three-dimensional structure. Overall, we analyzed 63 point-mutations from patients and identified common patterns underlying FV deficient phenotypes. We used structural and evolutionary patterns as input to machine learning algorithms to anticipate the effects of mutations and anticipated FV-deficiency with fair accuracy. Together, our results demonstrate how clinical features, genetic data and in silico analysis are converging to enhance treatment and diagnosis of coagulation disorders.
Collapse
Affiliation(s)
| | | | - Brenno M Alencar
- Institute of Computing, Federal University of Bahia, Salvador, Brazil
| | - Marcos V Ferreira
- Institute of Computing, Federal University of Bahia, Salvador, Brazil
| | - Tatiane Nogueira
- Institute of Computing, Federal University of Bahia, Salvador, Brazil
| | - Ricardo A Rios
- Institute of Computing, Federal University of Bahia, Salvador, Brazil
| | - Tiago J S Lopes
- Center for Regenerative Medicine, National Centre for Child Health and Development Research Institute, 2-10-1 Okura, Setagaya, Tokyo, 157-8535, Japan.
| |
Collapse
|
26
|
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]
|
27
|
Abildgaard AB, Nielsen SV, Bernstein I, Stein A, Lindorff-Larsen K, Hartmann-Petersen R. Lynch syndrome, molecular mechanisms and variant classification. Br J Cancer 2023; 128:726-734. [PMID: 36434153 PMCID: PMC9978028 DOI: 10.1038/s41416-022-02059-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 10/31/2022] [Accepted: 11/02/2022] [Indexed: 11/27/2022] Open
Abstract
Patients with the heritable cancer disease, Lynch syndrome, carry germline variants in the MLH1, MSH2, MSH6 and PMS2 genes, encoding the central components of the DNA mismatch repair system. Loss-of-function variants disrupt the DNA mismatch repair system and give rise to a detrimental increase in the cellular mutational burden and cancer development. The treatment prospects for Lynch syndrome rely heavily on early diagnosis; however, accurate diagnosis is inextricably linked to correct clinical interpretation of individual variants. Protein variant classification traditionally relies on cumulative information from occurrence in patients, as well as experimental testing of the individual variants. The complexity of variant classification is due to (1) that variants of unknown significance are rare in the population and phenotypic information on the specific variants is missing, and (2) that individual variant testing is challenging, costly and slow. Here, we summarise recent developments in high-throughput technologies and computational prediction tools for the assessment of variants of unknown significance in Lynch syndrome. These approaches may vastly increase the number of interpretable variants and could also provide important mechanistic insights into the disease. These insights may in turn pave the road towards developing personalised treatment approaches for Lynch syndrome.
Collapse
Affiliation(s)
- Amanda B Abildgaard
- The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Sofie V Nielsen
- The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Inge Bernstein
- Department of Surgical Gastroenterology, Aalborg University Hospital, Aalborg, Denmark
- Institute of Clinical Medicine, Aalborg University Hospital, Aalborg University, Aalborg, Denmark
| | - Amelie Stein
- The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Rasmus Hartmann-Petersen
- The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
28
|
Mendelian inheritance revisited: dominance and recessiveness in medical genetics. Nat Rev Genet 2023:10.1038/s41576-023-00574-0. [PMID: 36806206 DOI: 10.1038/s41576-023-00574-0] [Citation(s) in RCA: 18] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/14/2022] [Indexed: 02/22/2023]
Abstract
Understanding the consequences of genotype for phenotype (which ranges from molecule-level effects to whole-organism traits) is at the core of genetic diagnostics in medicine. Many measures of the deleteriousness of individual alleles exist, but these have limitations for predicting the clinical consequences. Various mechanisms can protect the organism from the adverse effects of functional variants, especially when the variant is paired with a wild type allele. Understanding why some alleles are harmful in the heterozygous state - representing dominant inheritance - but others only with the biallelic presence of pathogenic variants - representing recessive inheritance - is particularly important when faced with the deluge of rare genetic alterations identified by high throughput DNA sequencing. Both awareness of the specific quantitative and/or qualitative effects of individual variants and the elucidation of allelic and non-allelic interactions are essential to optimize genetic diagnosis and counselling.
Collapse
|
29
|
V M DD, Sivaramakrishnan V, Arvind Kumar K. Structural systems biology approach delineate the functional implications of SNPs in exon junction complex interaction network. J Biomol Struct Dyn 2023; 41:11969-11986. [PMID: 36617892 DOI: 10.1080/07391102.2022.2164355] [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: 02/03/2022] [Accepted: 12/26/2022] [Indexed: 01/10/2023]
Abstract
In eukaryotes, transcripts that carry premature termination codons (PTC) leading to truncated proteins are degraded by the Nonsense Mediated Decay (NMD) machinery. Missense and nonsense Single Nucleotide Polymorphisms (SNPs) in proteins belonging to Exon junction complex (EJC) and up-frameshift protein (UPF) will compromise NMD leading to the accumulation of truncated proteins in various diseases. The EJC and UPF which are involved in NMD is a good model system to study the effect of SNPs at a system level. Despite the availability of crystal structures, computational tools, and data on mutational and deletion studies, with functional implications, an integrated effort to understand the impact of SNPs at the systems level is lacking. To study the functional consequences of missense SNPs, sequence-based techniques like SIFT and PolyPhen which classify SNPs as deleterious or non-deleterious and structure-based methods like FoldX which calculate the Delta Delta G, (ddGs, ∆∆G) are used. Using FoldX, the ddG for mutations with experimentally validated functional effects is calculated and compared with those calculated for SNPs in the same protein-protein interaction interface. Further, a model is conceived to explain the functional implications of SNPs based on the effects observed for known mutants. The results are visualized in a network format. The effects of nonsense mutations are discerned by comparing with deletion mutation studies and loss of interaction in the crystal structure. The present work not only integrates genomics, proteomics, and classical genetics with 'Structural Biology' but also helps to integrate it into a 'systems-level functional network'.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Datta Darshan V M
- Disease Biology Lab, Department of Biosciences, Sri Sathya Sai Institute of Higher Learning, Prasanthi Nilayam, Anantapur, Andhra Pradesh, India
| | - Venketesh Sivaramakrishnan
- Disease Biology Lab, Department of Biosciences, Sri Sathya Sai Institute of Higher Learning, Prasanthi Nilayam, Anantapur, Andhra Pradesh, India
| | - K Arvind Kumar
- Disease Biology Lab, Department of Biosciences, Sri Sathya Sai Institute of Higher Learning, Prasanthi Nilayam, Anantapur, Andhra Pradesh, India
- Department of Physiology and Biophysics, School of Medicine, Case Western Reserve University, Cleveland, OH, USA
| |
Collapse
|
30
|
Sora V, Laspiur AO, Degn K, Arnaudi M, Utichi M, Beltrame L, De Menezes D, Orlandi M, Stoltze UK, Rigina O, Sackett PW, Wadt K, Schmiegelow K, Tiberti M, Papaleo E. RosettaDDGPrediction for high-throughput mutational scans: From stability to binding. Protein Sci 2023; 32:e4527. [PMID: 36461907 PMCID: PMC9795540 DOI: 10.1002/pro.4527] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Revised: 11/25/2022] [Accepted: 11/25/2022] [Indexed: 12/05/2022]
Abstract
Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein-protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction.
Collapse
Affiliation(s)
- Valentina Sora
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Adrian Otamendi Laspiur
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Kristine Degn
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Matteo Arnaudi
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Mattia Utichi
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Ludovica Beltrame
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Dayana De Menezes
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Matteo Orlandi
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Ulrik Kristoffer Stoltze
- Department of Clinical GeneticsCopenhagen University Hospital RigshospitaletCopenhagenDenmark
- Department of Pediatrics and Adolescent MedicineUniversity Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Olga Rigina
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Peter Wad Sackett
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| | - Karin Wadt
- Department of Clinical GeneticsCopenhagen University Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Kjeld Schmiegelow
- Department of Pediatrics and Adolescent MedicineUniversity Hospital RigshospitaletCopenhagenDenmark
- Institute of Clinical Medicine, Faculty of MedicineUniversity of CopenhagenCopenhagenDenmark
| | - Matteo Tiberti
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
| | - Elena Papaleo
- Cancer Structural Biology, Danish Cancer Society Research CenterCopenhagenDenmark
- Cancer Systems Biology, Section for Bioinformatics, Department of Health and TechnologyTechnical University of DenmarkLyngbyDenmark
| |
Collapse
|
31
|
Keskin Karakoyun H, Yüksel ŞK, Amanoglu I, Naserikhojasteh L, Yeşilyurt A, Yakıcıer C, Timuçin E, Akyerli CB. Evaluation of AlphaFold structure-based protein stability prediction on missense variations in cancer. Front Genet 2023; 14:1052383. [PMID: 36896237 PMCID: PMC9988940 DOI: 10.3389/fgene.2023.1052383] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Accepted: 02/08/2023] [Indexed: 02/23/2023] Open
Abstract
Identifying pathogenic missense variants in hereditary cancer is critical to the efforts of patient surveillance and risk-reduction strategies. For this purpose, many different gene panels consisting of different number and/or set of genes are available and we are particularly interested in a panel of 26 genes with a varying degree of hereditary cancer risk consisting of ABRAXAS1, ATM, BARD1, BLM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, EPCAM, MEN1, MLH1, MRE11, MSH2, MSH6, MUTYH, NBN, PALB2, PMS2, PTEN, RAD50, RAD51C, RAD51D, STK11, TP53, and XRCC2. In this study, we have compiled a collection of the missense variations reported in any of these 26 genes. More than a thousand missense variants were collected from ClinVar and the targeted screen of a breast cancer cohort of 355 patients which contributed to this set with 160 novel missense variations. We analyzed the impact of the missense variations on protein stability by five different predictors including both sequence- (SAAF2EC and MUpro) and structure-based (Maestro, mCSM, CUPSAT) predictors. For the structure-based tools, we have utilized the AlphaFold (AF2) protein structures which comprise the first structural analysis of this hereditary cancer proteins. Our results agreed with the recent benchmarks that computed the power of stability predictors in discriminating the pathogenic variants. Overall, we reported a low-to-medium-level performance for the stability predictors in discriminating pathogenic variants, except MUpro which had an AUROC of 0.534 (95% CI [0.499-0.570]). The AUROC values ranged between 0.614-0.719 for the total set and 0.596-0.682 for the set with high AF2 confidence regions. Furthermore, our findings revealed that the confidence score for a given variant in the AF2 structure could alone predict pathogenicity more robustly than any of the tested stability predictors with an AUROC of 0.852. Altogether, this study represents the first structural analysis of the 26 hereditary cancer genes underscoring 1) the thermodynamic stability predicted from AF2 structures as a moderate and 2) the confidence score of AF2 as a strong descriptor for variant pathogenicity.
Collapse
Affiliation(s)
- Hilal Keskin Karakoyun
- Department of Biochemistry and Molecular Biology, Institute of Health Sciences, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| | - Şirin K Yüksel
- Department of Biochemistry and Molecular Biology, Institute of Health Sciences, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| | - Ilayda Amanoglu
- Department of Biostatistics and Bioinformatics, Institute of Health Sciences, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| | - Lara Naserikhojasteh
- Department of Biostatistics and Bioinformatics, Institute of Health Sciences, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| | - Ahmet Yeşilyurt
- Acibadem Labgen Genetic Diagnosis Centre, Acibadem Health Group, Istanbul, Türkiye
| | - Cengiz Yakıcıer
- Acibadem Pathology Laboratories, Acibadem Health Group, Istanbul, Türkiye
| | - Emel Timuçin
- Department of Biostatistics and Medical Informatics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| | - Cemaliye B Akyerli
- Department of Medical Biology, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Türkiye
| |
Collapse
|
32
|
Carpentier M, Chomilier J. Analyses of Mutation Displacements from Homology Models. Methods Mol Biol 2023; 2627:195-210. [PMID: 36959449 DOI: 10.1007/978-1-0716-2974-1_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/25/2023]
Abstract
Evaluation of the structural perturbations introduced by a single amino acid mutation is the main issue for protein structural biology. We propose here to present some recent advances in methods, allowing the splitting of distortion between the actual substitution effect and the contribution of the local flexibility of the position where the mutation occurs. Its main drawback is the need of many structures with a single mutation in each of them. To bypass this difficulty, we propose to use molecular modeling tools, with several software enabling us to build a model from a template, given the sequence. As a proof of concept, we rely on a gold standard, the human lysozyme. Both wild-type and three mutant structures are available in the PDB. Two of these mutations result in amyloid fibril formation, and the last one is neutral. As a conclusion, irrespective of the algorithm used for modeling, side chain conformations at the site of mutation are reliable, although long-range effects are out of reach of these tools.
Collapse
Affiliation(s)
- Mathilde Carpentier
- Institut Systématique Evolution Biodiversité (ISYEB), Sorbonne Université, MNHN, CNRS, EPHE, Paris, France.
| | - Jacques Chomilier
- Sorbonne Université, BiBiP, IMPMC, UMR 7590, CNRS, MNHN, Paris, France
| |
Collapse
|
33
|
Valanciute A, Nygaard L, Zschach H, Maglegaard Jepsen M, Lindorff-Larsen K, Stein A. Accurate protein stability predictions from homology models. Comput Struct Biotechnol J 2022; 21:66-73. [PMID: 36514339 PMCID: PMC9729920 DOI: 10.1016/j.csbj.2022.11.048] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 11/22/2022] [Accepted: 11/23/2022] [Indexed: 11/27/2022] Open
Abstract
Calculating changes in protein stability (ΔΔG) has been shown to be central for predicting the consequences of single amino acid substitutions in protein engineering as well as interpretation of genomic variants for disease risk. Structure-based calculations are considered most accurate, however the tools used to calculate ΔΔGs have been developed on experimentally resolved structures. Extending those calculations to homology models based on related proteins would greatly extend their applicability as large parts of e.g. the human proteome are not structurally resolved. In this study we aim to investigate the accuracy of ΔΔG values predicted on homology models compared to crystal structures. Specifically, we identified four proteins with a large number of experimentally tested ΔΔGs and templates for homology modeling across a broad range of sequence identities, and selected three methods for ΔΔG calculations to test. We find that ΔΔG-values predicted from homology models compare equally well to experimental ΔΔGs as those predicted on experimentally established crystal structures, as long as the sequence identity of the model template to the target protein is at least 40%. In particular, the Rosetta cartesian_ddg protocol is robust against the small perturbations in the structure which homology modeling introduces. In an independent assessment, we observe a similar trend when using ΔΔGs to categorize variants as low or wild-type-like abundance. Overall, our results show that stability calculations performed on homology models can substitute for those on crystal structures with acceptable accuracy as long as the model is built on a template with sequence identity of at least 40% to the target protein.
Collapse
Affiliation(s)
- Audrone Valanciute
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Lasse Nygaard
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Henrike Zschach
- Section for Computational and RNA Biology, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Michael Maglegaard Jepsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark,Corresponding authors.
| | - Amelie Stein
- Section for Computational and RNA Biology, Department of Biology, University of Copenhagen, Copenhagen, Denmark,Corresponding authors.
| |
Collapse
|
34
|
Hall HN, Bengani H, Hufnagel RB, Damante G, Ansari M, Marsh JA, Grimes GR, von Kriegsheim A, Moore D, McKie L, Rahmat J, Mio C, Blyth M, Keng WT, Islam L, McEntargart M, Mannens MM, Heyningen VV, Rainger J, Brooks BP, FitzPatrick DR. Monoallelic variants resulting in substitutions of MAB21L1 Arg51 Cause Aniridia and microphthalmia. PLoS One 2022; 17:e0268149. [PMID: 36413568 PMCID: PMC9681113 DOI: 10.1371/journal.pone.0268149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Accepted: 10/06/2022] [Indexed: 11/23/2022] Open
Abstract
Classical aniridia is a congenital and progressive panocular disorder almost exclusively caused by heterozygous loss-of-function variants at the PAX6 locus. We report nine individuals from five families with severe aniridia and/or microphthalmia (with no detectable PAX6 mutation) with ultrarare monoallelic missense variants altering the Arg51 codon of MAB21L1. These mutations occurred de novo in 3/5 families, with the remaining families being compatible with autosomal dominant inheritance. Mice engineered to carry the p.Arg51Leu change showed a highly-penetrant optic disc anomaly in heterozygous animals with severe microphthalmia in homozygotes. Substitutions of the same codon (Arg51) in MAB21L2, a close homolog of MAB21L1, cause severe ocular and skeletal malformations in humans and mice. The predicted nucleotidyltransferase function of MAB21L1 could not be demonstrated using purified protein with a variety of nucleotide substrates and oligonucleotide activators. Induced expression of GFP-tagged wildtype and mutant MAB21L1 in human cells caused only modest transcriptional changes. Mass spectrometry of immunoprecipitated protein revealed that both mutant and wildtype MAB21L1 associate with transcription factors that are known regulators of PAX6 (MEIS1, MEIS2 and PBX1) and with poly(A) RNA binding proteins. Arg51 substitutions reduce the association of wild-type MAB21L1 with TBL1XR1, a component of the NCoR complex. We found limited evidence for mutation-specific interactions with MSI2/Musashi-2, an RNA-binding proteins with effects on many different developmental pathways. Given that biallelic loss-of-function variants in MAB21L1 result in a milder eye phenotype we suggest that Arg51-altering monoallelic variants most plausibly perturb eye development via a gain-of-function mechanism.
Collapse
Affiliation(s)
- Hildegard Nikki Hall
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Hemant Bengani
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Robert B. Hufnagel
- National Eye Institute, National Institutes of Health, Bethesda, MD, United States of America
| | | | - Morad Ansari
- South East Scotland Genetic Service, Western General Hospital, Edinburgh, United Kingdom
| | - Joseph A. Marsh
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Graeme R. Grimes
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Alex von Kriegsheim
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - David Moore
- South East Scotland Genetic Service, Western General Hospital, Edinburgh, United Kingdom
| | - Lisa McKie
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Jamalia Rahmat
- Ophthalmology Department, Hospital Kuala Lumpur, Kuala Lumpur, Malaysia
| | - Catia Mio
- Department of Medicine, University of Udine, Udine, Italy
| | - Moira Blyth
- University of Leeds, St. James’s University Hospital, Leeds, United Kingdom
| | - Wee Teik Keng
- Department of Genetics, Kuala Lumpur Hospital, Kuala Lumpur, Malaysia
| | - Lily Islam
- West Midlands Regional Genetics Service, Birmingham Women’s and Children’s NHS Foundation Trust, Birmingham, England
| | - Meriel McEntargart
- Medical Genetics, St George’s University Hospitals NHS Foundation Trust, London, United Kingdom
| | - Marcel M. Mannens
- Genome Diagnostics laboratory, Department of Clinical Genetics, Amsterdam University Medical Center, Amsterdam, The Netherlands
| | - Veronica Van Heyningen
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| | - Joe Rainger
- Roslin Institute, University of Edinburgh, Edinburgh, United Kingdom
| | - Brian P. Brooks
- National Eye Institute, National Institutes of Health, Bethesda, MD, United States of America
| | - David R. FitzPatrick
- MRC Human Genetics Unit, MRC Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
35
|
Attard TJ, Welburn JPI, Marsh JA. Understanding molecular mechanisms and predicting phenotypic effects of pathogenic tubulin mutations. PLoS Comput Biol 2022; 18:e1010611. [PMID: 36206299 PMCID: PMC9581425 DOI: 10.1371/journal.pcbi.1010611] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Revised: 10/19/2022] [Accepted: 09/28/2022] [Indexed: 11/21/2022] Open
Abstract
Cells rely heavily on microtubules for several processes, including cell division and molecular trafficking. Mutations in the different tubulin-α and -β proteins that comprise microtubules have been associated with various diseases and are often dominant, sporadic and congenital. While the earliest reported tubulin mutations affect neurodevelopment, mutations are also associated with other disorders such as bleeding disorders and infertility. We performed a systematic survey of tubulin mutations across all isotypes in order to improve our understanding of how they cause disease, and increase our ability to predict their phenotypic effects. Both protein structural analyses and computational variant effect predictors were very limited in their utility for differentiating between pathogenic and benign mutations. This was even worse for those genes associated with non-neurodevelopmental disorders. We selected tubulin-α and -β disease mutations that were most poorly predicted for experimental characterisation. These mutants co-localise to the mitotic spindle in HeLa cells, suggesting they may exert dominant-negative effects by altering microtubule properties. Our results show that tubulin mutations represent a blind spot for current computational approaches, being much more poorly predicted than mutations in most human disease genes. We suggest that this is likely due to their strong association with dominant-negative and gain-of-function mechanisms.
Collapse
Affiliation(s)
- Thomas J. Attard
- Wellcome Trust Centre for Cell Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, Scotland, United Kingdom
| | - Julie P. I. Welburn
- Wellcome Trust Centre for Cell Biology, School of Biological Sciences, University of Edinburgh, Edinburgh, Scotland, United Kingdom
| | - Joseph A. Marsh
- MRC Human Genetics Unit, Institute of Genetics & Cancer, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
36
|
Woodard J, Iqbal S, Mashaghi A. Circuit topology predicts pathogenicity of missense mutations. Proteins 2022; 90:1634-1644. [PMID: 35394672 PMCID: PMC9543832 DOI: 10.1002/prot.26342] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Revised: 03/07/2022] [Accepted: 03/30/2022] [Indexed: 12/05/2022]
Abstract
The contact topology of a protein determines important aspects of the folding process. The topological measure of contact order has been shown to be predictive of the rate of folding. Circuit topology is emerging as another fundamental descriptor of biomolecular structure, with predicted effects on the folding rate. We analyze the residue‐based circuit topological environments of 21 K mutations labeled as pathogenic or benign. Multiple statistical lines of reasoning support the conclusion that the number of contacts in two specific circuit topological arrangements, namely inverse parallel and cross relations, with contacts involving the mutated residue have discriminatory value in determining the pathogenicity of human variants. We investigate how results vary with residue type and according to whether the gene is essential. We further explore the relationship to a number of structural features and find that circuit topology provides nonredundant information on protein structures and pathogenicity of mutations. Results may have implications for the polymer physics of protein folding and suggest that “local” topological information, including residue‐based circuit topology and residue contact order, could be useful in improving state‐of‐the‐art machine learning algorithms for pathogenicity prediction.
Collapse
Affiliation(s)
- Jaie Woodard
- Medical Systems Biophysics and Bioengineering, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University, Leiden, The Netherlands.,Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, USA
| | - Sumaiya Iqbal
- Center for the Development of Therapeutics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA.,Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA.,Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA.,Analytic and Translational Genetics Unit, Massachusetts General Hospital, Boston, Massachusetts, USA
| | - Alireza Mashaghi
- Medical Systems Biophysics and Bioengineering, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University, Leiden, The Netherlands.,Centre for Interdisciplinary Genome Research, Faculty of Science, Leiden University, Leiden, The Netherlands
| |
Collapse
|
37
|
Iqbal S, Ge F, Li F, Akutsu T, Zheng Y, Gasser RB, Yu DJ, Webb GI, Song J. PROST: AlphaFold2-aware Sequence-Based Predictor to Estimate Protein Stability Changes upon Missense Mutations. J Chem Inf Model 2022; 62:4270-4282. [PMID: 35973091 DOI: 10.1021/acs.jcim.2c00799] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
An essential step in engineering proteins and understanding disease-causing missense mutations is to accurately model protein stability changes when such mutations occur. Here, we developed a new sequence-based predictor for the protein stability (PROST) change (Gibb's free energy change, ΔΔG) upon a single-point missense mutation. PROST extracts multiple descriptors from the most promising sequence-based predictors, such as BoostDDG, SAAFEC-SEQ, and DDGun. RPOST also extracts descriptors from iFeature and AlphaFold2. The extracted descriptors include sequence-based features, physicochemical properties, evolutionary information, evolutionary-based physicochemical properties, and predicted structural features. The PROST predictor is a weighted average ensemble model based on extreme gradient boosting (XGBoost) decision trees and an extra-trees regressor; PROST is trained on both direct and hypothetical reverse mutations using the S5294 (S2647 direct mutations + S2647 inverse mutations). The parameters for the PROST model are optimized using grid searching with 5-fold cross-validation, and feature importance analysis unveils the most relevant features. The performance of PROST is evaluated in a blinded manner, employing nine distinct data sets and existing state-of-the-art sequence-based and structure-based predictors. This method consistently performs well on frataxin, S217, S349, Ssym, S669, Myoglobin, and CAGI5 data sets in blind tests and similarly to the state-of-the-art predictors for p53 and S276 data sets. When the performance of PROST is compared with the latest predictors such as BoostDDG, SAAFEC-SEQ, ACDC-NN-seq, and DDGun, PROST dominates these predictors. A case study of mutation scanning of the frataxin protein for nine wild-type residues demonstrates the utility of PROST. Taken together, these findings indicate that PROST is a well-suited predictor when no protein structural information is available. The source code of PROST, data sets, examples, and pretrained models along with how to use PROST are available at https://github.com/ShahidIqb/PROST and https://prost.erc.monash.edu/seq.
Collapse
Affiliation(s)
- Shahid Iqbal
- Department of Data Science and AI, Faculty of IT, Monash University, Clayton, Victoria 3800, Australia.,Biomedicine Discovery Institute and Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria 3800, Australia.,Monash Data Futures Institute, Monash University, Clayton, Victoria 3800, Australia
| | - Fang Ge
- School of Computer Science and Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
| | - Fuyi Li
- Biomedicine Discovery Institute and Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria 3800, Australia.,Monash Data Futures Institute, Monash University, Clayton, Victoria 3800, Australia
| | - Tatsuya Akutsu
- Bioinformatics Center, Institute for Chemical Research, Kyoto University, Kyoto 611-0011, Japan
| | - Yuanting Zheng
- Department of Veterinary Biosciences, Melbourne Veterinary School, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Robin B Gasser
- Department of Veterinary Biosciences, Melbourne Veterinary School, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Dong-Jun Yu
- School of Computer Science and Engineering, Nanjing University of Science and Technology, Nanjing 210094, China
| | - Geoffrey I Webb
- Department of Data Science and AI, Faculty of IT, Monash University, Clayton, Victoria 3800, Australia.,Monash Data Futures Institute, Monash University, Clayton, Victoria 3800, Australia
| | - Jiangning Song
- Department of Data Science and AI, Faculty of IT, Monash University, Clayton, Victoria 3800, Australia.,Biomedicine Discovery Institute and Department of Biochemistry and Molecular Biology, Monash University, Clayton, Victoria 3800, Australia.,Monash Data Futures Institute, Monash University, Clayton, Victoria 3800, Australia
| |
Collapse
|
38
|
Loss-of-function, gain-of-function and dominant-negative mutations have profoundly different effects on protein structure. Nat Commun 2022; 13:3895. [PMID: 35794153 PMCID: PMC9259657 DOI: 10.1038/s41467-022-31686-6] [Citation(s) in RCA: 74] [Impact Index Per Article: 37.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 06/29/2022] [Indexed: 12/12/2022] Open
Abstract
Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Taking protein structure into account has therefore provided great insight into the molecular mechanisms underlying human genetic disease. While there has been much focus on how mutations can disrupt protein structure and thus cause a loss of function (LOF), alternative mechanisms, specifically dominant-negative (DN) and gain-of-function (GOF) effects, are less understood. Here, we investigate the protein-level effects of pathogenic missense mutations associated with different molecular mechanisms. We observe striking differences between recessive vs dominant, and LOF vs non-LOF mutations, with dominant, non-LOF disease mutations having much milder effects on protein structure, and DN mutations being highly enriched at protein interfaces. We also find that nearly all computational variant effect predictors, even those based solely on sequence conservation, underperform on non-LOF mutations. However, we do show that non-LOF mutations could potentially be identified by their tendency to cluster in three-dimensional space. Overall, our work suggests that many pathogenic mutations that act via DN and GOF mechanisms are likely being missed by current variant prioritisation strategies, but that there is considerable scope to improve computational predictions through consideration of molecular disease mechanisms. Most known pathogenic mutations occur in protein-coding regions of DNA and change the way proteins are made. Here the authors analyse the locations of thousands of human disease mutations and their predicted effects on protein structure and show that,while loss-of-function mutations tend to be highly disruptive, non-loss-of-function mutations are in general much milder at a protein structural level.
Collapse
|
39
|
Panda G, Mishra N, Sharma D, Kutum R, Bhoyar RC, Jain A, Imran M, Senthilvel V, Divakar MK, Mishra A, Garg P, Banerjee P, Sivasubbu S, Scaria V, Ray A. Comprehensive Assessment of Indian Variations in the Druggable Kinome Landscape Highlights Distinct Insights at the Sequence, Structure and Pharmacogenomic Stratum. Front Pharmacol 2022; 13:858345. [PMID: 35865963 PMCID: PMC9294532 DOI: 10.3389/fphar.2022.858345] [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: 01/19/2022] [Accepted: 06/07/2022] [Indexed: 11/13/2022] Open
Abstract
India confines more than 17% of the world’s population and has a diverse genetic makeup with several clinically relevant rare mutations belonging to many sub-group which are undervalued in global sequencing datasets like the 1000 Genome data (1KG) containing limited samples for Indian ethnicity. Such databases are critical for the pharmaceutical and drug development industry where diversity plays a crucial role in identifying genetic disposition towards adverse drug reactions. A qualitative and comparative sequence and structural study utilizing variant information present in the recently published, largest curated Indian genome database (IndiGen) and the 1000 Genome data was performed for variants belonging to the kinase coding genes, the second most targeted group of drug targets. The sequence-level analysis identified similarities and differences among different populations based on the nsSNVs and amino acid exchange frequencies whereas a comparative structural analysis of IndiGen variants was performed with pathogenic variants reported in UniProtKB Humsavar data. The influence of these variations on structural features of the protein, such as structural stability, solvent accessibility, hydrophobicity, and the hydrogen-bond network was investigated. In-silico screening of the known drugs to these Indian variation-containing proteins reveals critical differences imparted in the strength of binding due to the variations present in the Indian population. In conclusion, this study constitutes a comprehensive investigation into the understanding of common variations present in the second largest population in the world and investigating its implications in the sequence, structural and pharmacogenomic landscape. The preliminary investigation reported in this paper, supporting the screening and detection of ADRs specific to the Indian population could aid in the development of techniques for pre-clinical and post-market screening of drug-related adverse events in the Indian population.
Collapse
Affiliation(s)
- Gayatri Panda
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla, India
| | - Neha Mishra
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla, India
| | - Disha Sharma
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Rintu Kutum
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
- Ashoka University, Sonipat, India
| | - Rahul C. Bhoyar
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Abhinav Jain
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Mohamed Imran
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Vigneshwar Senthilvel
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Mohit Kumar Divakar
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Anushree Mishra
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Parth Garg
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla, India
| | - Priyanka Banerjee
- Institute for Physiology, Charité-University Medicine Berlin, Berlin, Germany
| | - Sridhar Sivasubbu
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Vinod Scaria
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad, India
- CSIR-Institute of Genomics and Integrative Biology, Delhi, India
| | - Arjun Ray
- Department of Computational Biology, Indraprastha Institute of Information Technology, Okhla, India
- *Correspondence: Arjun Ray,
| |
Collapse
|
40
|
Ferla MP, Pagnamenta AT, Koukouflis L, Taylor JC, Marsden BD. Venus: Elucidating the Impact of Amino Acid Variants on Protein Function Beyond Structure Destabilisation. J Mol Biol 2022; 434:167567. [PMID: 35662467 PMCID: PMC9742853 DOI: 10.1016/j.jmb.2022.167567] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2021] [Revised: 03/11/2022] [Accepted: 03/22/2022] [Indexed: 12/15/2022]
Abstract
Exploring the functional effect of a non-synonymous coding variant at the protein level requires multiple pieces of information to be interpreted appropriately. This is particularly important when embarking on the study of a potentially pathogenic variant linked to a rare or monogenic disease. Whereas accurate protein stability predictions alone are generally informative, other effects, such as disruption of post-translational modifications or weakened ligand binding, may also contribute to the disease phenotype. Furthermore, consideration of nearby variants that are found in the healthy population may strengthen or refute a given mechanistic hypothesis. Whilst there are several bioinformatics tools available that score a genetic variant in terms of deleteriousness, there is no single tool that assembles multiple effects of a variant on the encoded protein, beyond structural stability, and presents them on the structure for inspection. Venus is a web application which, given a protein substitution, rapidly estimates the predicted effect on protein stability of the variant, flags if the variant affects a post-translational modification site, a predicted linear motif or known annotation, and determines the effect on protein stability of variants which affect nearby residues and have been identified in healthy populations. Venus is built upon Michelanglo and the results can be exported to it, allowing them to be annotated and shared with other researchers. Venus is freely accessible at https://venus.cmd.ox.ac.uk and its source code is openly available at https://github.com/CMD-Oxford/Michelanglo-and-Venus.
Collapse
Affiliation(s)
- Matteo P Ferla
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK; Oxford NIHR Biomedical Research Centre, Oxford, UK.
| | - Alistair T Pagnamenta
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK; Oxford NIHR Biomedical Research Centre, Oxford, UK. https://twitter.com/@alistairp2011
| | - Leonidas Koukouflis
- Centre for Medicines Discovery, University of Oxford, Old Road Campus Research Building, Oxford OX3 7DQ, UK
| | - Jenny C Taylor
- Wellcome Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK; Oxford NIHR Biomedical Research Centre, Oxford, UK
| | - Brian D Marsden
- Centre for Medicines Discovery, University of Oxford, Old Road Campus Research Building, Oxford OX3 7DQ, UK; Kennedy Institute of Rheumatology, University of Oxford, Oxford OX3 7FY, UK. https://twitter.com/@bmarsden19
| |
Collapse
|
41
|
Livesey BJ, Marsh JA. Interpreting protein variant effects with computational predictors and deep mutational scanning. Dis Model Mech 2022; 15:275742. [PMID: 35736673 PMCID: PMC9235876 DOI: 10.1242/dmm.049510] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Computational predictors of genetic variant effect have advanced rapidly in recent years. These programs provide clinical and research laboratories with a rapid and scalable method to assess the likely impacts of novel variants. However, it can be difficult to know to what extent we can trust their results. To benchmark their performance, predictors are often tested against large datasets of known pathogenic and benign variants. These benchmarking data may overlap with the data used to train some supervised predictors, which leads to data re-use or circularity, resulting in inflated performance estimates for those predictors. Furthermore, new predictors are usually found by their authors to be superior to all previous predictors, which suggests some degree of computational bias in their benchmarking. Large-scale functional assays known as deep mutational scans provide one possible solution to this problem, providing independent datasets of variant effect measurements. In this Review, we discuss some of the key advances in predictor methodology, current benchmarking strategies and how data derived from deep mutational scans can be used to overcome the issue of data circularity. We also discuss the ability of such functional assays to directly predict clinical impacts of mutations and how this might affect the future need for variant effect predictors.
Collapse
Affiliation(s)
- Benjamin J Livesey
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh EH4 2XU, UK
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh EH4 2XU, UK
| |
Collapse
|
42
|
Joshi A, Bhardwaj D, Kaushik A, Juneja VK, Taneja P, Thakur S, Kumra Taneja N. Advances in multi-omics based quantitative microbial risk assessment in the dairy sector: A semi-systematic review. Food Res Int 2022; 156:111323. [PMID: 35651076 DOI: 10.1016/j.foodres.2022.111323] [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: 11/09/2021] [Revised: 04/26/2022] [Accepted: 04/28/2022] [Indexed: 11/16/2022]
Abstract
With the increasing consumption of packaged and ready-to-eat food products, the risk of foodborne illness has drastically increased and so has the dire need for proper management. The conventional Microbial Risk Assessment (MRA) investigations require prior knowledge of process flow, exposure, and hazard assessment throughout the supply chain. These data are often generated using conventional microbiological approaches based either on shelf-life studies or specific spoilage organisms (SSOs), frequently overlooking crucial information such as antimicrobial resistance (AMR), biofilm formation, virulence factors and other physiological variations coupled with bio-chemical characteristics of food matrix. Additionally, the microbial risks in food are diverse and heterogenous, that might be an outcome of growth and activity of multiple microbial populations rather than a single species contamination. The uncertainty on the microbial source, time as well as point of entry into the food supply chain poses a constraint to the efficiency of preventive approaches and conventional MRA. In the last few decades, significant breakthroughs in molecular methods and continuously progressing bioinformatics tools have opened up a new horizon for risk analysis-based approaches in food safety. Real time polymerase chain reaction (qPCR) and kit-based assays provide better accuracy and precision with shorter processing time. Despite these improvements, the effect of complex food matrix on growth environment and recovery of pathogen is a persistent problem for risk assessors. The dairy industry is highly impacted by spoilage and pathogenic microorganisms. Therefore, this review discusses the evolution and recent advances in MRAmethodologies equipped with predictive interventions and "multi-omics" approach for robust MRA specifically targeting dairy products. It also highlights the limiting gap area and the opportunity for improvement in this field to ensure precision food safety.
Collapse
Affiliation(s)
- Akanksha Joshi
- Dept. of Basic and Applied Sciences, National Institute of Food Technology Entrepreneurship and Management, Haryana 131028, India
| | - Dinesh Bhardwaj
- Dept. of Basic and Applied Sciences, National Institute of Food Technology Entrepreneurship and Management, Haryana 131028, India
| | - Abhishek Kaushik
- Dept. of Basic and Applied Sciences, National Institute of Food Technology Entrepreneurship and Management, Haryana 131028, India
| | | | - Pankaj Taneja
- Department of Biotechnology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Sheetal Thakur
- Department of Food Science and Technology, MMICT & BM (HM), MMDU, Mullana, Ambala, Haryana, India
| | - Neetu Kumra Taneja
- Dept. of Basic and Applied Sciences, National Institute of Food Technology Entrepreneurship and Management, Haryana 131028, India; Center for Advance Translational Research in Food Nanobiotechnology (CATR-FNB), National Institute of Food Technology Entrepreneurship and Management, Haryana 131028, India.
| |
Collapse
|
43
|
Montanucci L, Capriotti E, Birolo G, Benevenuta S, Pancotti C, Lal D, Fariselli P. DDGun: an untrained predictor of protein stability changes upon amino acid variants. Nucleic Acids Res 2022; 50:W222-W227. [PMID: 35524565 PMCID: PMC9252764 DOI: 10.1093/nar/gkac325] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 04/15/2022] [Accepted: 05/04/2022] [Indexed: 01/22/2023] Open
Abstract
Estimating the functional effect of single amino acid variants in proteins is fundamental for predicting the change in the thermodynamic stability, measured as the difference in the Gibbs free energy of unfolding, between the wild-type and the variant protein (ΔΔG). Here, we present the web-server of the DDGun method, which was previously developed for the ΔΔG prediction upon amino acid variants. DDGun is an untrained method based on basic features derived from evolutionary information. It is antisymmetric, as it predicts opposite ΔΔG values for direct (A → B) and reverse (B → A) single and multiple site variants. DDGun is available in two versions, one based on only sequence information and the other one based on sequence and structure information. Despite being untrained, DDGun reaches prediction performances comparable to those of trained methods. Here we make DDGun available as a web server. For the web server version, we updated the protein sequence database used for the computation of the evolutionary features, and we compiled two new data sets of protein variants to do a blind test of its performances. On these blind data sets of single and multiple site variants, DDGun confirms its prediction performance, reaching an average correlation coefficient between experimental and predicted ΔΔG of 0.45 and 0.49 for the sequence-based and structure-based versions, respectively. Besides being used for the prediction of ΔΔG, we suggest that DDGun should be adopted as a benchmark method to assess the predictive capabilities of newly developed methods. Releasing DDGun as a web-server, stand-alone program and docker image will facilitate the necessary process of method comparison to improve ΔΔG prediction.
Collapse
Affiliation(s)
- Ludovica Montanucci
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, 9500 Euclid Avenue, Cleveland, OH 44195, USA
| | - Emidio Capriotti
- BioFolD Unit, Department of Pharmacy and Biotechnology (FaBiT), University of Bologna, Via F. Selmi 3, 40126 Bologna, Italy
| | - Giovanni Birolo
- Department of Medical Sciences, University of Torino, Via Santena 19, 10126, Torino, Italy
| | - Silvia Benevenuta
- Department of Medical Sciences, University of Torino, Via Santena 19, 10126, Torino, Italy
| | - Corrado Pancotti
- Department of Medical Sciences, University of Torino, Via Santena 19, 10126, Torino, Italy
| | - Dennis Lal
- Genomic Medicine Institute, Lerner Research Institute, Cleveland Clinic, 9500 Euclid Avenue, Cleveland, OH 44195, USA
| | - Piero Fariselli
- Department of Medical Sciences, University of Torino, Via Santena 19, 10126, Torino, Italy
| |
Collapse
|
44
|
Backwell L, Marsh JA. Diverse Molecular Mechanisms Underlying Pathogenic Protein Mutations: Beyond the Loss-of-Function Paradigm. Annu Rev Genomics Hum Genet 2022; 23:475-498. [PMID: 35395171 DOI: 10.1146/annurev-genom-111221-103208] [Citation(s) in RCA: 39] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Most known disease-causing mutations occur in protein-coding regions of DNA. While some of these involve a loss of protein function (e.g., through premature stop codons or missense changes that destabilize protein folding), many act via alternative molecular mechanisms and have dominant-negative or gain-of-function effects. In nearly all cases, these non-loss-of-function mutations can be understood by considering interactions of the wild-type and mutant protein with other molecules, such as proteins, nucleic acids, or small ligands and substrates. Here, we review the diverse molecular mechanisms by which pathogenic mutations can have non-loss-of-function effects, including by disrupting interactions, increasing binding affinity, changing binding specificity, causing assembly-mediated dominant-negative and dominant-positive effects, creating novel interactions, and promoting aggregation and phase separation. We believe that increased awareness of these diverse molecular disease mechanisms will lead to improved diagnosis (and ultimately treatment) of human genetic disorders. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Lisa Backwell
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, United Kingdom;
| | - Joseph A Marsh
- MRC Human Genetics Unit, Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, United Kingdom;
| |
Collapse
|
45
|
Ghadie MA, Xia Y. Are transient protein-protein interactions more dispensable? PLoS Comput Biol 2022; 18:e1010013. [PMID: 35404956 PMCID: PMC9000134 DOI: 10.1371/journal.pcbi.1010013] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2021] [Accepted: 03/11/2022] [Indexed: 12/12/2022] Open
Abstract
Protein-protein interactions (PPIs) are key drivers of cell function and evolution. While it is widely assumed that most permanent PPIs are important for cellular function, it remains unclear whether transient PPIs are equally important. Here, we estimate and compare dispensable content among transient PPIs and permanent PPIs in human. Starting with a human reference interactome mapped by experiments, we construct a human structural interactome by building three-dimensional structural models for PPIs, and then distinguish transient PPIs from permanent PPIs using several structural and biophysical properties. We map common mutations from healthy individuals and disease-causing mutations onto the structural interactome, and perform structure-based calculations of the probabilities for common mutations (assumed to be neutral) and disease mutations (assumed to be mildly deleterious) to disrupt transient PPIs and permanent PPIs. Using Bayes' theorem we estimate that a similarly small fraction (<~20%) of both transient and permanent PPIs are completely dispensable, i.e., effectively neutral upon disruption. Hence, transient and permanent interactions are subject to similarly strong selective constraints in the human interactome.
Collapse
Affiliation(s)
| | - Yu Xia
- Department of Bioengineering, McGill University, Montreal, Canada
| |
Collapse
|
46
|
Pan Q, Nguyen TB, Ascher DB, Pires DEV. Systematic evaluation of computational tools to predict the effects of mutations on protein stability in the absence of experimental structures. Brief Bioinform 2022; 23:bbac025. [PMID: 35189634 PMCID: PMC9155634 DOI: 10.1093/bib/bbac025] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2021] [Revised: 01/13/2022] [Accepted: 01/30/2022] [Indexed: 12/26/2022] Open
Abstract
Changes in protein sequence can have dramatic effects on how proteins fold, their stability and dynamics. Over the last 20 years, pioneering methods have been developed to try to estimate the effects of missense mutations on protein stability, leveraging growing availability of protein 3D structures. These, however, have been developed and validated using experimentally derived structures and biophysical measurements. A large proportion of protein structures remain to be experimentally elucidated and, while many studies have based their conclusions on predictions made using homology models, there has been no systematic evaluation of the reliability of these tools in the absence of experimental structural data. We have, therefore, systematically investigated the performance and robustness of ten widely used structural methods when presented with homology models built using templates at a range of sequence identity levels (from 15% to 95%) and contrasted performance with sequence-based tools, as a baseline. We found there is indeed performance deterioration on homology models built using templates with sequence identity below 40%, where sequence-based tools might become preferable. This was most marked for mutations in solvent exposed residues and stabilizing mutations. As structure prediction tools improve, the reliability of these predictors is expected to follow, however we strongly suggest that these factors should be taken into consideration when interpreting results from structure-based predictors of mutation effects on protein stability.
Collapse
Affiliation(s)
- Qisheng Pan
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane City, Queensland 4072, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, 30 Flemington Rd, Parkville, Victoria 3052, Australia
| | - Thanh Binh Nguyen
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane City, Queensland 4072, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, 30 Flemington Rd, Parkville, Victoria 3052, Australia
| | - David B Ascher
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane City, Queensland 4072, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, 30 Flemington Rd, Parkville, Victoria 3052, Australia
- Department of Biochemistry, University of Cambridge, 80 Tennis Ct Rd, Cambridge CB2 1GA, UK
| | - Douglas E V Pires
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Victoria 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane City, Queensland 4072, Australia
- Systems and Computational Biology, Bio21 Institute, University of Melbourne, 30 Flemington Rd, Parkville, Victoria 3052, Australia
- School of Computing and Information Systems, University of Melbourne, Melbourne, Victoria 3053, Australia
| |
Collapse
|
47
|
Population-level deficit of homozygosity unveils CPSF3 as an intellectual disability syndrome gene. Nat Commun 2022; 13:705. [PMID: 35121750 PMCID: PMC8817032 DOI: 10.1038/s41467-022-28330-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Accepted: 01/20/2022] [Indexed: 02/05/2023] Open
Abstract
AbstractPredicting the pathogenicity of biallelic missense variants can be challenging. Here, we use a deficit of observed homozygous carriers of missense variants, versus an expected number in a set of 153,054 chip-genotyped Icelanders, to identify potentially pathogenic genotypes. We follow three missense variants with a complete deficit of homozygosity and find that their pathogenic effect in homozygous state ranges from severe childhood disease to early embryonic lethality. One of these variants is in CPSF3, a gene not previously linked to disease. From a set of clinically sequenced Icelanders, and by sequencing archival samples targeted through the Icelandic genealogy, we find four homozygous carriers. Additionally, we find two homozygous carriers of Mexican descent of another missense variant in CPSF3. All six homozygous carriers of missense variants in CPSF3 show severe intellectual disability, seizures, microcephaly, and abnormal muscle tone. Here, we show how the absence of certain homozygous genotypes from a large population set can elucidate causes of previously unexplained recessive diseases and early miscarriage.
Collapse
|
48
|
Massively parallel phenotyping of coding variants in cancer with Perturb-seq. Nat Biotechnol 2022; 40:896-905. [DOI: 10.1038/s41587-021-01160-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2020] [Accepted: 11/11/2021] [Indexed: 02/08/2023]
|
49
|
Høie MH, Cagiada M, Beck Frederiksen AH, Stein A, Lindorff-Larsen K. Predicting and interpreting large-scale mutagenesis data using analyses of protein stability and conservation. Cell Rep 2022; 38:110207. [PMID: 35021073 DOI: 10.1016/j.celrep.2021.110207] [Citation(s) in RCA: 43] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2021] [Revised: 10/01/2021] [Accepted: 12/13/2021] [Indexed: 01/23/2023] Open
Abstract
Understanding and predicting the functional consequences of single amino acid changes is central in many areas of protein science. Here, we collect and analyze experimental measurements of effects of >150,000 variants in 29 proteins. We use biophysical calculations to predict changes in stability for each variant and assess them in light of sequence conservation. We find that the sequence analyses give more accurate prediction of variant effects than predictions of stability and that about half of the variants that show loss of function do so due to stability effects. We construct a machine learning model to predict variant effects from protein structure and sequence alignments and show how the two sources of information support one another and enable mechanistic interpretations. Together, our results show how one can leverage large-scale experimental assessments of variant effects to gain deeper and general insights into the mechanisms that cause loss of function.
Collapse
Affiliation(s)
- Magnus Haraldson Høie
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Matteo Cagiada
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Anders Haagen Beck Frederiksen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark
| | - Amelie Stein
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark.
| | - Kresten Lindorff-Larsen
- Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, DK-2200 Copenhagen N, Denmark.
| |
Collapse
|
50
|
Woodard J, Zheng W, Zhang Y. Protein structural features predict responsiveness to pharmacological chaperone treatment for three lysosomal storage disorders. PLoS Comput Biol 2021; 17:e1009370. [PMID: 34529671 PMCID: PMC8478239 DOI: 10.1371/journal.pcbi.1009370] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 09/28/2021] [Accepted: 08/21/2021] [Indexed: 12/15/2022] Open
Abstract
Three-dimensional structures of proteins can provide important clues into the efficacy of personalized treatment. We perform a structural analysis of variants within three inherited lysosomal storage disorders, comparing variants responsive to pharmacological chaperone treatment to those unresponsive to such treatment. We find that predicted ΔΔG of mutation is higher on average for variants unresponsive to treatment, in the case of datasets for both Fabry disease and Pompe disease, in line with previous findings. Using both a single decision tree and an advanced machine learning approach based on the larger Fabry dataset, we correctly predict responsiveness of three Gaucher disease variants, and we provide predictions for untested variants. Many variants are predicted to be responsive to treatment, suggesting that drug-based treatments may be effective for a number of variants in Gaucher disease. In our analysis, we observe dependence on a topological feature reporting on contact arrangements which is likely connected to the order of folding of protein residues, and we provide a potential justification for this observation based on steady-state cellular kinetics. Pharmacological chaperones are small molecule drugs that bind to proteins to help stabilize the folded state. One set of diseases for which this treatment has been effective is the lysosomal storage disorders, which are caused by defective lysosomal enzymes. However, not all genotypes are equally responsive to treatment. For instance, missense mutants that are particularly destabilized relative to WT are less likely to respond. The availability of datasets containing responsiveness data for large numbers of mutants, along with crystal structures of the protein involved in each disease, make machine learning methods incorporating sequence-based and structural data feasible. We hypothesize that data from two diseases, Fabry and Pompe disease, may be useful for predicting responsiveness of variants in the related Gaucher disease. Results suggest that many rare variants in Gaucher disease could be amenable to existing drugs. Results also suggest that drug responsiveness depends on protein topology in such a way that mutations in early-to-fold residues are more likely to be non-responsive to pharmacological chaperone treatment, which is consistent with a simple kinetic model of stability rescue. This study provides an example of how machine learning can be used to inform further studies towards personalized treatment in medicine.
Collapse
Affiliation(s)
- Jaie Woodard
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Wei Zheng
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Yang Zhang
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan, United States of America
- Department of Biological Chemistry, University of Michigan, Ann Arbor, Michigan, United States of America
- * E-mail:
| |
Collapse
|