1
|
Maharaj AV, Ishida M, Rybak A, Elfeky R, Andrews A, Joshi A, Elmslie F, Joensuu A, Kantojärvi K, Jia RY, Perry JRB, O'Toole EA, McGuffin LJ, Hwa V, Storr HL. QSOX2 Deficiency-induced short stature, gastrointestinal dysmotility and immune dysfunction. Nat Commun 2024; 15:8420. [PMID: 39341815 PMCID: PMC11439042 DOI: 10.1038/s41467-024-52587-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 09/13/2024] [Indexed: 10/01/2024] Open
Abstract
Postnatal growth failure is often attributed to dysregulated somatotropin action, however marked genetic and phenotypic heterogeneity exist. We report five patients from three families who present with short stature, immune dysfunction, atopic eczema and gastrointestinal pathology associated with recessive variants in QSOX2. QSOX2 encodes a nuclear membrane protein linked to disulphide isomerase and oxidoreductase activity. Loss of QSOX2 disrupts Growth hormone-mediated STAT5B nuclear translocation despite enhanced Growth hormone-induced STAT5B phosphorylation. Moreover, patient-derived dermal fibroblasts demonstrate Growth hormone-induced mitochondriopathy and reduced mitochondrial membrane potential. Located at the nuclear membrane, QSOX2 acts as a gatekeeper for regulating stabilisation and import of phosphorylated-STAT5B. Altogether, QSOX2 deficiency modulates human growth by impairing Growth hormone-STAT5B downstream activities and mitochondrial dynamics, which contribute to multi-system dysfunction. Furthermore, our work suggests that therapeutic recombinant insulin-like growth factor-1 may circumvent the Growth hormone-STAT5B dysregulation induced by pathological QSOX2 variants and potentially alleviate organ specific disease.
Collapse
Affiliation(s)
- Avinaash V Maharaj
- Centre for Endocrinology, John Vane Science Centre, Queen Mary University of London, Charterhouse Square, London, UK.
| | - Miho Ishida
- Centre for Endocrinology, John Vane Science Centre, Queen Mary University of London, Charterhouse Square, London, UK
| | - Anna Rybak
- Gastroenterology Department, Great Ormond Street Hospital, London, UK
| | - Reem Elfeky
- Immunology Department, Great Ormond Street Hospital, London, UK
| | - Afiya Andrews
- Centre for Endocrinology, John Vane Science Centre, Queen Mary University of London, Charterhouse Square, London, UK
| | - Aakash Joshi
- St George's University Hospitals NHS Foundation Trust, London, UK
| | - Frances Elmslie
- Genomics Clinical Academic Group, St George's University Hospitals NHS Foundation Trust, London, UK
| | - Anni Joensuu
- THL Biobank, the Department of Knowledge Brokers, Finnish Institute for Health and Welfare, Helsinki, Finland
| | - Katri Kantojärvi
- THL Biobank, the Department of Knowledge Brokers, Finnish Institute for Health and Welfare, Helsinki, Finland
| | - Raina Y Jia
- MRC Epidemiology Unit, University of Cambridge, School of Clinical Medicine, Cambridge, UK
| | - John R B Perry
- MRC Epidemiology Unit, University of Cambridge, School of Clinical Medicine, Cambridge, UK
| | - Edel A O'Toole
- Centre for Cell Biology and Cutaneous Research, Blizard Institute, Queen Mary University of London, London, UK
| | - Liam J McGuffin
- School of Biological Sciences, University of Reading, Reading, UK
| | - Vivian Hwa
- Cincinnati Children's Hospital Medical Center, Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio, USA.
- Premium Research Institute for Human Metaverse Medicine (WPI-PRIMe), Osaka University, Suita, Osaka, Japan.
| | - Helen L Storr
- Centre for Endocrinology, John Vane Science Centre, Queen Mary University of London, Charterhouse Square, London, UK.
| |
Collapse
|
2
|
Khan MU, Sakhawat A, Rehman R, Wali AH, Ghani MU, Akram A, Javed MA, Ali Q, Yu-Ming Z, Ali D, Yu-Ming Z. Identification of novel natural compounds against CFTR p.Gly628Arg pathogenic variant. AMB Express 2024; 14:99. [PMID: 39249658 PMCID: PMC11383896 DOI: 10.1186/s13568-024-01762-9] [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: 06/20/2024] [Accepted: 08/27/2024] [Indexed: 09/10/2024] Open
Abstract
Cystic fibrosis transmembrane conductance regulator (CFTR) protein is an ion channel found in numerous epithelia and controls the flow of water and salt across the epithelium. The aim of our study to find natural compounds that can improve lung function for people with cystic fibrosis (CF) caused by the p.Gly628Arg (rs397508316) mutation of CFTR protein. The sequence of CFTR protein as a target structure was retrieved from UniProt and PDB database. The ligands that included Armepavine, Osthole, Curcumin, Plumbagine, Quercetin, and one Trikafta (R*) reference drug were screened out from PubChem database. Autodock vina software carried out docking, and binding energies between the drug and the target were included using docking-score. The following tools examined binding energy, interaction, stability, toxicity, and visualize protein-ligand complexes. The compounds having binding energies of -6.4, -5.1, -6.6, -5.1, and - 6.5 kcal/mol for Armepavine, Osthole, Curcumin, Plumbagine, Quercetin, and R*-drug, respectively with mutated CFTR (Gly628Arg) structure were chosen as the most promising ligands. The ligands bind to the mutated CFTR protein structure active sites in hydrophobic bonds, hydrogen bonds, and electrostatic interactions. According to ADMET analyses, the ligands Armepavine and Quercetin also displayed good pharmacokinetic and toxicity characteristics. An MD simulation for 200 ns was also established to ensure that Armepavine and Quercetin ligands attached to the target protein favorably and dynamically, and that protein-ligand complex stability was maintained. It is concluded that Armepavine and Quercetin have stronger capacity to inhibit the effect of mutated CFTR protein through improved trafficking and restoration of original function.
Collapse
Affiliation(s)
- Muhammad Umer Khan
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan.
| | - Azra Sakhawat
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Raima Rehman
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Abbas Haider Wali
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Muhammad Usman Ghani
- Precision Genomics Research Lab, Centre for Applied Molecular Biology, University of the Punjab, Lahore, Pakistan
| | - Areeba Akram
- Precision Genomics Research Lab, Centre for Applied Molecular Biology, University of the Punjab, Lahore, Pakistan
| | - Muhammad Arshad Javed
- Department of Plant Breeding and Genetics, Faculty of Agricultural Sciences, University of the Punjab, Lahore, Pakistan
| | - Qurban Ali
- Department of Plant Breeding and Genetics, Faculty of Agricultural Sciences, University of the Punjab, Lahore, Pakistan.
| | - Zhou Yu-Ming
- Department of Emergency, The Affiliated Ganzhou Hospital of Nanchang University, Ganzhou, 341000, Jiangxi Province, P.R. China
| | - Daoud Ali
- Department of Zoology, College of Science, King Saud University, PO Box 2455, Riyadh, 11451, Saudi Arabia
| | - Zhou Yu-Ming
- Department of Emergency, The Affiliated Ganzhou Hospital of Nanchang University, Ganzhou, 341000, Jiangxi Province, P.R. China
| |
Collapse
|
3
|
Velecký J, Berezný M, Musil M, Damborsky J, Bednar D, Mazurenko S. BenchStab: a tool for automated querying of web-based stability predictors. BIOINFORMATICS (OXFORD, ENGLAND) 2024; 40:btae553. [PMID: 39259175 PMCID: PMC11427696 DOI: 10.1093/bioinformatics/btae553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2024] [Revised: 08/02/2024] [Accepted: 09/10/2024] [Indexed: 09/12/2024]
Abstract
SUMMARY Protein design requires information about how mutations affect protein stability. Many web-based predictors are available for this purpose, yet comparing them or using them en masse is difficult. Here, we present BenchStab, a console tool/Python package for easy and quick execution of 19 predictors and result collection on a list of mutants. Moreover, the tool is easily extensible with additional predictors. We created an independent dataset derived from the FireProtDB and evaluated 24 different prediction methods. AVAILABILITY AND IMPLEMENTATION BenchStab is an open-source Python package available at https://github.com/loschmidt/BenchStab with a detailed README and example usage at https://loschmidt.chemi.muni.cz/benchstab. The BenchStab dataset is available on Zenodo: https://zenodo.org/records/10637728.
Collapse
Affiliation(s)
- Jan Velecký
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic
| | - Matej Berezný
- Department of Information Systems, Faculty of Information Technology, Brno University of Technology, 612 00 Brno, Czech Republic
| | - Milos Musil
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic
- Department of Information Systems, Faculty of Information Technology, Brno University of Technology, 612 00 Brno, Czech Republic
- International Clinical Research Centre, St. Anne's University Hospital, 656 91 Brno, Czech Republic
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic
- International Clinical Research Centre, St. Anne's University Hospital, 656 91 Brno, Czech Republic
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic
- International Clinical Research Centre, St. Anne's University Hospital, 656 91 Brno, Czech Republic
| | - Stanislav Mazurenko
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic
- International Clinical Research Centre, St. Anne's University Hospital, 656 91 Brno, Czech Republic
| |
Collapse
|
4
|
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
|
5
|
Kamal MM, Islam MN, Rabby MG, Zahid MA, Hasan MM. In Silico Functional and Structural Analysis of Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in Human Paired Box 4 Gene. Biochem Genet 2024; 62:2975-2998. [PMID: 38062275 DOI: 10.1007/s10528-023-10589-1] [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: 07/08/2023] [Accepted: 11/06/2023] [Indexed: 07/31/2024]
Abstract
In human genome, members of Paired box (PAX) transcription factor family are highly sequence-specific DNA-binding proteins. Among PAX gene family members, PAX4 gene has significant role in growth, proliferation, differentiation, and insulin secretion of pancreatic β-cells. Single nucleotide polymorphisms (SNPs) in PAX4 gene progress in the pathogenesis of various human diseases. Hence, the molecular mechanism of how these SNPs in PAX4 gene significantly progress diseases pathogenesis needs to be elucidated. For the reason, a series of bioinformatic analyzes were done to identify the SNPs of PAX4 gene that contribute in diseases pathogenesis. From the analyzes, 4145 SNPs (rsIDs) in PAX4 gene were obtained, where, 362 missense (8.73%), 169 synonymous (4.08%), and 2323 intron variants (56.04%). The rest SNPs were unspecified. Among the 362 missense variants, 118 nsSNPs were found as deleterious in SIFT analysis. Among those, 25 nsSNPs were most probably damaging and 23 were deleterious as observed in PolyPhen-2 and PROVEAN analyzes, respectively. Following all analyzes, 14 nsSNPs (rs149708455, rs115887120, rs147279315, rs35155575, rs370095957, rs373939873, rs145468905, rs121917718, rs2233580, rs3824004, rs372751660, rs369459316, rs375472849, rs372497946) were common and observed as deleterious, probably damaging, affective and diseases associated. Following structural analyzes, 11 nsSNPs guided proteins were found as most unstable and highly conserved. Among these, R20W, R39Q, R45Q, R60H, G65D, and A223D mutated proteins were highly harmful. Hence, the results from above-mentioned integrated comprehensive bioinformatic analyzes guide how different nsSNPs in PAX4 gene alter structural and functional characteristics of the protein that might progress diseases pathogenesis in human including type 2 diabetes.
Collapse
Affiliation(s)
- Md Mostafa Kamal
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Md Numan Islam
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh
- Department of Food Engineering, North Pacific International University of Bangladesh, Dhaka, Bangladesh
| | - Md Golam Rabby
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Md Ashrafuzzaman Zahid
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Md Mahmudul Hasan
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh.
| |
Collapse
|
6
|
Akram A, Sakhawat A, Ghani MU, Khan MU, Rehman R, Ali Q, Jin-Liang P, Ali D. Silibinins and curcumin as promising ligands against mutant cystic fibrosis transmembrane regulator protein. AMB Express 2024; 14:84. [PMID: 39043981 PMCID: PMC11266341 DOI: 10.1186/s13568-024-01742-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2024] [Accepted: 07/08/2024] [Indexed: 07/25/2024] Open
Abstract
Cystic Fibrosis Transmembrane Regulator (CFTR) is a significant protein that is responsible for the movement of ions across cell membranes. The cystic fibrosis (CF) occur due to the mutations in the CFTR gene as it produces the dysfunctional CFTR protein. The sequence of CFTR protein as a target structure was retrieved from UniProt and PDB database. The ligands selection was performed through virtual screening and top 3 ligands choose out of 65 ligands silibinins, curcumin, demethoxycurcumin were selected with a reference drug Trikafta (R*). According to docking, ADMET analyses, the natural ligands (Silibinins and Curcumin) displayed best binding energy, pharmacokinetic and free toxicity than other natural compounds and reference drug (R*). An MD simulation for 200 ns was also established to ensure that natural ligands (Silibinins and Curcumin) attached to the target protein favorably and dynamically, and that protein-ligand complex stability was maintained. It is concluded that silibinins and curcumins have a better capacity to decrease the effect of mutant CFTR protein through improved trafficking and the restoration of original function. In conclusion, in silico studies demonstrate the potential of silibinins and curcumin as therapeutic agents for cystic fibrosis, particularly for the D614G mutated protein. Their ability to increase CFTR function while reducing cellular stress and inflammation, together with their favorable safety profile and accessibility could make them valuable additions to cystic fibrosis treatment options. Further experimental and clinical validation will be required to fully realize their potential and include them into effective therapy regimens.
Collapse
Affiliation(s)
- Areeba Akram
- Precision Genomics Research Lab, Centre for Applied Molecular Biology, University of the Punjab, Lahore, Pakistan
| | - Azra Sakhawat
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan
| | - Muhammad Usman Ghani
- Precision Genomics Research Lab, Centre for Applied Molecular Biology, University of the Punjab, Lahore, Pakistan
| | - Muhammad Umer Khan
- Institute of Molecular Biology and Biotechnology, The University of Lahore, Lahore, Pakistan.
| | - Raima Rehman
- Centre of Excellence in Molecular Biology, University of the Punjab, Lahore, Pakistan
| | - Qurban Ali
- Department of Plant Breeding and Genetics, Faculty of Agricultural Sciences, University of the Punjab, Lahore, Pakistan.
| | - Peng Jin-Liang
- Department of Emergency, The Affiliated Ganzhou Hospital of Nanchang University, Ganzhou, 341000, Jiangxi, People's Republic of China
| | - Daoud Ali
- Department of Zoology, College of Science, King Saud University, PO Box 2455, Riyadh, 11451, Saudi Arabia
| |
Collapse
|
7
|
Diaz DJ, Gong C, Ouyang-Zhang J, Loy JM, Wells J, Yang D, Ellington AD, Dimakis AG, Klivans AR. Stability Oracle: a structure-based graph-transformer framework for identifying stabilizing mutations. Nat Commun 2024; 15:6170. [PMID: 39043654 PMCID: PMC11266546 DOI: 10.1038/s41467-024-49780-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: 10/04/2023] [Accepted: 06/14/2024] [Indexed: 07/25/2024] Open
Abstract
Engineering stabilized proteins is a fundamental challenge in the development of industrial and pharmaceutical biotechnologies. We present Stability Oracle: a structure-based graph-transformer framework that achieves SOTA performance on accurately identifying thermodynamically stabilizing mutations. Our framework introduces several innovations to overcome well-known challenges in data scarcity and bias, generalization, and computation time, such as: Thermodynamic Permutations for data augmentation, structural amino acid embeddings to model a mutation with a single structure, a protein structure-specific attention-bias mechanism that makes transformers a viable alternative to graph neural networks. We provide training/test splits that mitigate data leakage and ensure proper model evaluation. Furthermore, to examine our data engineering contributions, we fine-tune ESM2 representations (Prostata-IFML) and achieve SOTA for sequence-based models. Notably, Stability Oracle outperforms Prostata-IFML even though it was pretrained on 2000X less proteins and has 548X less parameters. Our framework establishes a path for fine-tuning structure-based transformers to virtually any phenotype, a necessary task for accelerating the development of protein-based biotechnologies.
Collapse
Affiliation(s)
- Daniel J Diaz
- UT Austin, Department of Computer Science, Austin, TX, 78712, USA.
- Intelligent Proteins, LLC, Austin, TX, 78712, USA.
- UT Austin, Department of Chemistry, Austin, TX, 78712, USA.
| | - Chengyue Gong
- UT Austin, Department of Computer Science, Austin, TX, 78712, USA
| | | | - James M Loy
- Intelligent Proteins, LLC, Austin, TX, 78712, USA
- UT Austin, Department of Molecular Biosciences, Austin, TX, 78712, USA
| | - Jordan Wells
- UT Austin, McKetta Department of Chemical Engineering, Austin, TX, 78712, USA
| | - David Yang
- UT Austin, Department of Molecular Biosciences, Austin, TX, 78712, USA
| | | | - Alexandros G Dimakis
- UT Austin, Chandra Family Department of Electrical and Computer Engineering, Austin, TX, 78712, USA
| | - Adam R Klivans
- UT Austin, Department of Computer Science, Austin, TX, 78712, USA
| |
Collapse
|
8
|
Mozibullah M, Khatun M, Sikder M, Islam M, Sharmin M. Identifying Oncogenic Missense Single Nucleotide Polymorphisms in Human SAT1 Gene Using Computational Algorithms and Molecular Dynamics Tools. Cancer Rep (Hoboken) 2024; 7:e2130. [PMID: 39041636 PMCID: PMC11264109 DOI: 10.1002/cnr2.2130] [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: 10/29/2023] [Revised: 05/04/2024] [Accepted: 06/30/2024] [Indexed: 07/24/2024] Open
Abstract
BACKGROUND The human SAT1 gene encodes spermidine/spermine N1-acetyltransferase 1 (SSAT1), a regulatory biological catalyst of polyamine catabolism. Numerous essential biological processes, such as cellular proliferation, differentiation, and survival, depend on polyamines like spermidine and spermine. Thus, SSAT1 is involved in key cellular activities such as proliferation and survival of cells and mediates various diseases including cancer. A plethora of studies established the involvement of missense single nucleotide polymorphisms (SNPs) in numerous pathological conditions due to their ability to adversely affect the structure and subsequent function of the protein. AIMS To date, an in silico study to identify the pathogenic missense SNPs of the human SAT1 gene has not been accomplished yet. This study aimed to filter the missense SNPs that were functionally detrimental and pathogenic. METHODS AND RESULTS The rs757435207 (I21N) was ascertained to be the most deleterious and pathogenic by all algorithmic tools. Stability and evolutionary conservation analysis tools also stated that I21N variant decreased the stability and was located in the highly conserved residue. Molecular dynamics simulation revealed that I21N caused substantial alterations in the conformational stability and dynamics of the SSAT1 protein. Consequently, the I21N variant could disrupt the native functional roles of the SSAT1 enzyme. CONCLUSION Therefore, the I21N variant was identified and concluded to be an oncogenic missense variant of the human SAT1 gene. Overall, the findings of this study would be a great directory of future experimental research to develop personalized medicine.
Collapse
Affiliation(s)
- Md. Mozibullah
- Department of Biochemistry and Molecular BiologyMawlana Bhashani Science and Technology UniversityTangailBangladesh
| | - Marina Khatun
- Department of Biochemistry and Molecular BiologyMawlana Bhashani Science and Technology UniversityTangailBangladesh
| | - Md. Asaduzzaman Sikder
- Department of Biochemistry and Molecular BiologyMawlana Bhashani Science and Technology UniversityTangailBangladesh
| | - Mohammod Johirul Islam
- Department of Biochemistry and Molecular BiologyMawlana Bhashani Science and Technology UniversityTangailBangladesh
| | - Mehbuba Sharmin
- Department of Biochemistry and Molecular BiologyMawlana Bhashani Science and Technology UniversityTangailBangladesh
| |
Collapse
|
9
|
Kamal MM, Teeya ST, Rahman MM, Talukder MEK, Sarmin S, Wani TA, Hasan MM. Prediction and assessment of deleterious and disease causing nonsynonymous single nucleotide polymorphisms (nsSNPs) in human FOXP4 gene: An in - silico study. Heliyon 2024; 10:e32791. [PMID: 38994097 PMCID: PMC11237951 DOI: 10.1016/j.heliyon.2024.e32791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Accepted: 06/10/2024] [Indexed: 07/13/2024] Open
Abstract
In humans, FOXP gene family is involved in embryonic development and cancer progression. The FOXP4 (Forkhead box protein P4) gene belongs to this FOXP gene family. FOXP4 gene plays a crucial role in oncogenesis. Single nucleotide polymorphisms are biological markers and common determinants of human diseases. Mutations can largely affect the function of the corresponding protein. Therefore, the molecular mechanism of nsSNPs in the FOXP4 gene needs to be elucidated. Initially, the SNPs of the FOXP4 gene were extracted from the dbSNP database and a total of 23124 SNPs was found, where 555 nonsynonymous, 20525 intronic, 1114 noncoding transcript, 334 synonymous were obtained and the rest were unspecified. Then, a series of bioinformatics tools (SIFT, PolyPhen2, SNAP2, PhD SNP, PANTHER, I-Mutant2.0, MUpro, GOR IV, ConSurf, NetSurfP 2.0, HOPE, DynaMut2, GeneMANIA, STRING and Schrodinger) were used to explore the effect of nsSNPs on FOXP4 protein function and structural stability. First, 555 nsSNPs were analyzed using SIFT, of which 57 were found as deleterious. Following, PolyPhen2, SNAP2, PhD SNP and PANTHER analyses, 10 nsSNPs (rs372762294, rs141899153, rs142575732, rs376938850, rs367607523, rs112517943, rs140387832, rs373949416, rs373949416 and rs376160648) were common and observed as deleterious, damaging and diseases associated. Following that, using I-Mutant2.0 and MUpro servers, 7 nsSNPs were found to be the most unstable. GOR IV predicted that these seven nsSNPs affect protein structure by altering the protein contents of alpha helixes, extended strands, and random coils. Following DynaMut2, 5 nsSNPs showed a decrease in the ΔΔG value compared with the wild-type and were found to be responsible for destabilizing the corresponding protein. GeneMANIA and STRING network revealed interaction of FOXP4 with other genes. Finally, molecular dynamics simulation analysis revealed consistent fluctuation in RMSD and RMSF values, Rg and hydrogen bonds in the mutant proteins compared with WT, which might alter the functional and structural stability of the corresponding protein. As a result, the aforementioned integrated comprehensive bioinformatic analyses provide insight into how various nsSNPs of the FOXP4 gene change the structural and functional properties of the corresponding protein, potentially proceeding with the pathophysiology of human diseases.
Collapse
Affiliation(s)
- Md Mostafa Kamal
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
- Laboratory of Computational Biology, Biological Solution Centre, Jashore, 7408, Bangladesh
| | - Shamiha Tabassum Teeya
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| | - Md Mahfuzur Rahman
- Department of Genetic Engineering & Biotechnology, Bangabandhu Sheikh Mujibur Rahman Maritime University, Dhaka, 1216, Bangladesh
| | - Md Enamul Kabir Talukder
- Department of Genetic Engineering & Biotechnology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
- Laboratory of Computational Biology, Biological Solution Centre, Jashore, 7408, Bangladesh
| | - Sonia Sarmin
- BIRTAN-Bangladesh Institute of Research and Training on Applied Nutrition, Jhenaidah, 7300, Bangladesh
| | - Tanveer A Wani
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, 11451, Saudi Arabia
| | - Md Mahmudul Hasan
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| |
Collapse
|
10
|
da Silva-Júnior AHP, de Oliveira Silva RC, Gurgel APAD, Barros-Júnior MR, Nascimento KCG, Santos DL, Pena LJ, Lima RDCP, Batista MVDA, Chagas BS, de Freitas AC. Identification and Functional Implications of the E5 Oncogene Polymorphisms of Human Papillomavirus Type 16. Trop Med Infect Dis 2024; 9:140. [PMID: 39058182 PMCID: PMC11281449 DOI: 10.3390/tropicalmed9070140] [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: 04/27/2024] [Revised: 06/12/2024] [Accepted: 06/21/2024] [Indexed: 07/28/2024] Open
Abstract
The persistence of the human papillomavirus type 16 (HPV16) infection on the cervical epithelium contributes to the progression of cervical cancer. Studies have demonstrated that HPV16 genetic variants may be associated with different risks of developing cervical cancer. However, the E5 oncoprotein of HPV16, which is related to several cellular mechanisms in the initial phases of the infection and thus contributes to carcinogenesis, is still little studied. Here we investigate the HPV16 E5 oncogene variants to assess the effects of different mutations on the biological function of the E5 protein. We detected and analyzed the HPV16 E5 oncogene polymorphisms and their phylogenetic relationships. After that, we proposed a tertiary structure analysis of the protein variants, preferential codon usage, and functional activity of the HPV16 E5 protein. Intra-type variants were grouped in the lineages A and D using in silico analysis. The mutations in E5 were located in the T-cell epitopes region. We therefore analyzed the interference of the HPV16 E5 protein in the NF-kB pathway. Our results showed that the variants HPV16E5_49PE and HPV16E5_85PE did not increase the potential of the pathway activation capacity. This study provides additional knowledge about the mechanisms of dispersion of the HPV16 E5 variants, providing evidence that these variants may be relevant to the modulation of the NF-κB signaling pathway.
Collapse
Affiliation(s)
- Antônio Humberto P. da Silva-Júnior
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Ruany Cristyne de Oliveira Silva
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Ana Pavla A. Diniz Gurgel
- Department of Engineering and Environment, Federal University of Paraiba, João Pessoa 58033-455, Paraíba, Brazil;
| | - Marconi Rêgo Barros-Júnior
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Kamylla Conceição Gomes Nascimento
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Daffany Luana Santos
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Lindomar J. Pena
- Laboratory of Virology and Experimental Therapy, Instituto Aggeu Magalhães (IAM), Oswaldo Cruz Foundation, Recife 50670-901, Pernambuco, Brazil;
| | - Rita de Cássia Pereira Lima
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Marcus Vinicius de Aragão Batista
- Laboratory of Molecular Genetics and Biotechnology (GMBio), Department of Biology, Federal University of Sergipe, São Cristóvão 49107-230, Sergipe, Brazil
| | - Bárbara Simas Chagas
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| | - Antonio Carlos de Freitas
- Laboratory of Molecular Studies and Experimental Therapy (LEMTE), Department of Genetics, Federal University of Pernambuco, Recife 50670-901, Pernambuco, Brazil; (A.H.P.d.S.-J.); (R.C.d.O.S.); (M.R.B.-J.); (K.C.G.N.); (D.L.S.); (R.d.C.P.L.); (B.S.C.)
| |
Collapse
|
11
|
Ruiz-Bayón A, Cara-Rodríguez C, Sarmiento-Mañús R, Muñoz-Viana R, Lozano FM, Ponce MR, Micol JL. Roles of the Arabidopsis KEULE Gene in Postembryonic Development. Int J Mol Sci 2024; 25:6667. [PMID: 38928373 PMCID: PMC11204279 DOI: 10.3390/ijms25126667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2024] [Revised: 06/14/2024] [Accepted: 06/16/2024] [Indexed: 06/28/2024] Open
Abstract
Cytokinesis in plant cells begins with the fusion of vesicles that transport cell wall materials to the center of the cell division plane, where the cell plate forms and expands radially until it fuses with the parental cell wall. Vesicle fusion is facilitated by trans-SNARE complexes, with assistance from Sec1/Munc18 (SM) proteins. The SNARE protein KNOLLE and the SM protein KEULE are required for membrane fusion at the cell plate. Due to the crucial function of KEULE, all Arabidopsis (Arabidopsis thaliana) keule mutants identified to date are seedling lethal. Here, we identified the Arabidopsis serrata4-1 (sea4-1) and sea4-2 mutants, which carry recessive, hypomorphic alleles of KEULE. Homozygous sea4-1 and sea4-2 plants are viable and fertile but have smaller rosettes and fewer leaves at bolting than the wild type. Their leaves are serrated, small, and wavy, with a complex venation pattern. The mutant leaves also develop necrotic patches and undergo premature senescence. RNA-seq revealed transcriptome changes likely leading to reduced cell wall integrity and an increase in the unfolded protein response. These findings shed light on the roles of KEULE in postembryonic development, particularly in the patterning of rosette leaves and leaf margins.
Collapse
Affiliation(s)
| | | | | | | | | | | | - José Luis Micol
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain; (A.R.-B.); (C.C.-R.); (R.S.-M.); (R.M.-V.); (F.M.L.); (M.R.P.)
| |
Collapse
|
12
|
Ho CH, Chen CW, Su PY. Membranome-based identification of amino acid substitution in Haemophilus influenzae multidrug efflux pump HmrM for reduced chloramphenicol susceptibility. Arch Microbiol 2024; 206:298. [PMID: 38860999 DOI: 10.1007/s00203-024-04025-0] [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/19/2024] [Revised: 05/21/2024] [Accepted: 05/30/2024] [Indexed: 06/12/2024]
Abstract
A decreased chloramphenicol susceptibility in Haemophilus influenzae is commonly caused by the activity of chloramphenicol acetyltransferases (CATs). However, the involvement of membrane proteins in chloramphenicol susceptibility in H. influenzae remains unclear. In this study, chloramphenicol susceptibility testing, whole-genome sequencing, and analyses of membrane-related genes were performed in 51 H. influenzae isolates. Functional complementation assays and structure-based protein analyses were conducted to assess the effect of proteins with sequence substitutions on the minimum inhibitory concentration (MIC) of chloramphenicol in CAT-negative H. influenzae isolates. Six isolates were resistant to chloramphenicol and positive for type A-2 CATs. Of these isolates, A3256 had a similar level of CAT activity but a higher chloramphenicol MIC relative to the other resistant isolates; it also had 163 specific variations in 58 membrane genes. Regarding the CAT-negative isolates, logistic regression and receiver operator characteristic curve analyses revealed that 48T > G (Asn16Lys), 85 C > T (Leu29Phe), and 88 C > A (Leu30Ile) in HI_0898 (emrA), and 86T > G (Phe29Cys) and 141T > A (Ser47Arg) in HI_1177 (artM) were associated with enhanced chloramphenicol susceptibility, whereas 997G > A (Val333Ile) in HI_1612 (hmrM) was associated with reduced chloramphenicol susceptibility. Furthermore, the chloramphenicol MIC was lower in the CAT-negative isolates with EmrA-Leu29Phe/Leu30Ile or ArtM-Ser47Arg substitution and higher in those with HmrM-Val333Ile substitution, relative to their counterparts. The Val333Ile substitution was associated with enhanced HmrM protein stability and flexibility and increased chloramphenicol MICs in CAT-negative H. influenzae isolates. In conclusion, the substitution in H. influenzae multidrug efflux pump HmrM associated with reduced chloramphenicol susceptibility was characterised.
Collapse
Affiliation(s)
- Cheng-Hsun Ho
- Department of Medical Laboratory Science, College of Medical Science and Technology, I-Shou University, No.8, Yida Road, Jiaosu Village, Yanchao District, Kaohsiung, 82445, Taiwan.
| | - Chi-Wei Chen
- Graduate Degree Program of Smart Healthcare & Bioinformatics, College of Medical Science and Technology, I-Shou University, Kaohsiung, Taiwan
- Department of Biomedical Engineering, College of Medical Science and Technology, I-Shou University, Kaohsiung, Taiwan
| | - Pei-Yi Su
- Department of Laboratory Medicine, E-DA Hospital, Kaohsiung, Taiwan
| |
Collapse
|
13
|
Mahur P, Sharma A, Jahan G, S G A, Kumar Singh A, Muthukumaran J, Jain M. Understanding Genetic Risks: Computational Exploration of Human β-Synuclein nsSNPs and their Potential Impact on Structural Alteration. Neurosci Lett 2024; 833:137826. [PMID: 38768940 DOI: 10.1016/j.neulet.2024.137826] [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: 04/12/2024] [Revised: 05/16/2024] [Accepted: 05/17/2024] [Indexed: 05/22/2024]
Abstract
Synucleins are pivotal in neurodegenerative conditions. Beta-synuclein (β-synuclein) is part of the synuclein protein family alongside alpha-synuclein (α-synuclein) and gamma-synuclein (γ-synuclein). These proteins, found mainly in brain tissue and cancers, are soluble and unstructured. β-synuclein shares significant similarity with α-synuclein, especially in their N-terminus, with a 90% match. However, their aggregation tendencies differ significantly. While α-synuclein aggregation is believed to be counteracted by β-synuclein, which occurs in conditions like Parkinson's disease, β-synuclein may counteract α-synuclein's toxic effects on the nervous system, offering potential treatment for neurodegenerative diseases. Under normal circumstances, β-synuclein may guard against disease by interacting with α-synuclein. Yet, in pathological environments with heightened levels or toxic substances, it might contribute to disease. Our research aims to explore potential harmful mutations in the β-synuclein using computational tools to predict their destabilizing impact on protein structure. Consensus analysis revealed rs1207608813 (A63P), rs1340051870 (S72F), and rs1581178262 (G36C) as deleterious. These findings highlight the intricate relationship between nsSNPs and protein function, shedding light on their potential implications in disease pathways. Understanding the structural consequences of nsSNPs is crucial for elucidating their role in pathogenesis and developing targeted therapeutic interventions. Our results offer a robust computational framework for identifying neurodegenerative disorder-related mutations from SNP datasets, potentially reducing the costs associated with experimental characterization.
Collapse
Affiliation(s)
- Pragati Mahur
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Abhishek Sharma
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Gulnaz Jahan
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Adithya S G
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Amit Kumar Singh
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India
| | - Jayaraman Muthukumaran
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India.
| | - Monika Jain
- Department of Biotechnology, Sharda School of Engineering and Technology, Sharda University, Greater Noida, Uttar Pradesh, India.
| |
Collapse
|
14
|
Velloso JPL, de Sá AGC, Pires DEV, Ascher DB. Engineering G protein-coupled receptors for stabilization. Protein Sci 2024; 33:e5000. [PMID: 38747401 PMCID: PMC11094779 DOI: 10.1002/pro.5000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2023] [Revised: 03/21/2024] [Accepted: 04/10/2024] [Indexed: 05/19/2024]
Abstract
G protein-coupled receptors (GPCRs) are one of the most important families of targets for drug discovery. One of the limiting steps in the study of GPCRs has been their stability, with significant and time-consuming protein engineering often used to stabilize GPCRs for structural characterization and drug screening. Unfortunately, computational methods developed using globular soluble proteins have translated poorly to the rational engineering of GPCRs. To fill this gap, we propose GPCR-tm, a novel and personalized structurally driven web-based machine learning tool to study the impacts of mutations on GPCR stability. We show that GPCR-tm performs as well as or better than alternative methods, and that it can accurately rank the stability changes of a wide range of mutations occurring in various types of class A GPCRs. GPCR-tm achieved Pearson's correlation coefficients of 0.74 and 0.46 on 10-fold cross-validation and blind test sets, respectively. We observed that the (structural) graph-based signatures were the most important set of features for predicting destabilizing mutations, which points out that these signatures properly describe the changes in the environment where the mutations occur. More specifically, GPCR-tm was able to accurately rank mutations based on their effect on protein stability, guiding their rational stabilization. GPCR-tm is available through a user-friendly web server at https://biosig.lab.uq.edu.au/gpcr_tm/.
Collapse
Affiliation(s)
- João Paulo L. Velloso
- School of Chemistry and Molecular Biosciences, The Australian Centre for EcogenomicsThe University of QueenslandBrisbaneQueenslandAustralia
- Computational Biology and Clinical InformaticsBaker Heart and Diabetes InstituteMelbourneVictoriaAustralia
- Baker Department of Cardiometabolic HealthThe University of MelbourneParkvilleVictoriaAustralia
| | - Alex G. C. de Sá
- School of Chemistry and Molecular Biosciences, The Australian Centre for EcogenomicsThe University of QueenslandBrisbaneQueenslandAustralia
- Computational Biology and Clinical InformaticsBaker Heart and Diabetes InstituteMelbourneVictoriaAustralia
- Baker Department of Cardiometabolic HealthThe University of MelbourneParkvilleVictoriaAustralia
| | - Douglas E. V. Pires
- School of Computing and Information SystemsThe University of MelbourneParkvilleVictoriaAustralia
| | - David B. Ascher
- School of Chemistry and Molecular Biosciences, The Australian Centre for EcogenomicsThe University of QueenslandBrisbaneQueenslandAustralia
- Computational Biology and Clinical InformaticsBaker Heart and Diabetes InstituteMelbourneVictoriaAustralia
- Baker Department of Cardiometabolic HealthThe University of MelbourneParkvilleVictoriaAustralia
| |
Collapse
|
15
|
Qiu Y, Huang T, Cai YD. Review of predicting protein stability changes upon variations. Proteomics 2024; 24:e2300371. [PMID: 38643379 DOI: 10.1002/pmic.202300371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Revised: 04/07/2024] [Accepted: 04/08/2024] [Indexed: 04/22/2024]
Abstract
Forecasting alterations in protein stability caused by variations holds immense importance. Improving the thermal stability of proteins is important for biomedical and industrial applications. This review discusses the latest methods for predicting the effects of mutations on protein stability, databases containing protein mutations and thermodynamic parameters, and experimental techniques for efficiently assessing protein stability in high-throughput settings. Various publicly available databases for protein stability prediction are introduced. Furthermore, state-of-the-art computational approaches for anticipating protein stability changes due to variants are reviewed. Each method's types of features, base algorithm, and prediction results are also detailed. Additionally, some experimental approaches for verifying the prediction results of computational methods are introduced. Finally, the review summarizes the progress and challenges of protein stability prediction and discusses potential models for future research directions.
Collapse
Affiliation(s)
- Yiling Qiu
- Bio-Med Big Data Center, CAS Key Laboratory of Computational Biology, Shanghai Institute of Nutrition and Health, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai, China
- School of Mathematics and Statistics, Guangdong University of Technology, Guangzhou, China
| | - Tao Huang
- Bio-Med Big Data Center, CAS Key Laboratory of Computational Biology, Shanghai Institute of Nutrition and Health, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai, China
| | - Yu-Dong Cai
- School of Life Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
16
|
Storozhuk O, Bruekner SR, Paul A, Lebbink JHG, Sixma TK, Friedhoff P. MutL Activates UvrD by Interaction Between the MutL C-terminal Domain and the UvrD 2B Domain. J Mol Biol 2024; 436:168589. [PMID: 38677494 DOI: 10.1016/j.jmb.2024.168589] [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/16/2024] [Revised: 04/14/2024] [Accepted: 04/19/2024] [Indexed: 04/29/2024]
Abstract
UvrD is a helicase vital for DNA replication and quality control processes. In its monomeric state, UvrD exhibits limited helicase activity, necessitating either dimerization or assistance from an accessory protein to efficiently unwind DNA. Within the DNA mismatch repair pathway, MutL plays a pivotal role in relaying the repair signal, enabling UvrD to unwind DNA from the strand incision site up to and beyond the mismatch. Although this interdependence is well-established, the precise mechanism of activation and the specific MutL-UvrD interactions that trigger helicase activity remain elusive. To address these questions, we employed site-specific crosslinking techniques using single-cysteine variants of MutL and UvrD followed by functional assays. Our investigation unveils that the C-terminal domain of MutL not only engages with UvrD but also acts as a self-sufficient activator of UvrD helicase activity on DNA substrates with 3'-single-stranded tails. Especially when MutL is covalently attached to the 2B or 1B domain the tail length can be reduced to a minimal substrate of 5 nucleotides without affecting unwinding efficiency.
Collapse
Affiliation(s)
- Olha Storozhuk
- Institute for Biochemistry, FB 08, Justus Liebig University, Heinrich-Buff-Ring 17, D-35392 Giessen, Germany
| | - Susanne R Bruekner
- Division of Biochemistry, Netherlands Cancer Institute and Oncode Institute, Amsterdam, the Netherlands
| | - Ankon Paul
- Department of Molecular Genetics, Oncode Institute, Erasmus MC Cancer Institute, Erasmus University Medical Center, Rotterdam, the Netherlands
| | - Joyce H G Lebbink
- Department of Molecular Genetics, Oncode Institute, Erasmus MC Cancer Institute, Erasmus University Medical Center, Rotterdam, the Netherlands; Department of Radiotherapy, Erasmus University Medical Center, Rotterdam, the Netherlands
| | - Titia K Sixma
- Division of Biochemistry, Netherlands Cancer Institute and Oncode Institute, Amsterdam, the Netherlands
| | - Peter Friedhoff
- Institute for Biochemistry, FB 08, Justus Liebig University, Heinrich-Buff-Ring 17, D-35392 Giessen, Germany.
| |
Collapse
|
17
|
Alwan IH, Aljubouri TRS, Al-Shuhaib MBS. A Novel Missense SNP in the Fatty Acid-Binding Protein 4 (FABP4) Gene is Associated with Growth Traits in Karakul and Awassi Sheep. Biochem Genet 2024; 62:1462-1484. [PMID: 37640973 DOI: 10.1007/s10528-023-10504-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Accepted: 08/14/2023] [Indexed: 08/31/2023]
Abstract
The fatty acid-binding protein 4 (FABP4) plays a crucial role in the transportation and metabolism of fatty acids. It binds to long-chain fatty acids and facilitates their transport within cells. FABP4 is involved in lipid metabolism, insulin sensitivity, inflammation, and energy homeostasis. This study was conducted to assess the association between the FABP4 gene and growth traits in Karakul and Awassi sheep. A PCR-single strand conformation polymorphism (SSCP) protocol was utilized to assess the polymorphism of FABP4 PCR products with growth traits measured at birth, three, six, nine, and twelve-month intervals. One non-synonymous SNP was identified in the second exon, in which threonine was converted to aspartate in the 61st position in FABP4 (p.61Thr > Asp). This novel SNP showed significant associations with all growth traits measured at all age intervals. The results showed that lambs with the TT genotype exhibited higher growth traits than those with the GT and GG genotypes, respectively. The conducted prediction showed a clearly deleterious effect of p.61Thr > Asp on FABP4 structure, which was accompanied by reduced fatty acid binding efficiency. Owing to the predicted damaging effects caused by p.61Thr > Asp on FABP, lower levels of lipid transport and its consequent increased weight gain and other growth trait indices are expected. Therefore, a putative mechanism through which lambs with these genotypes exhibit higher growth traits is proposed. The FABP4 gene is suggested as a promising marker to improve growth traits in Karakul and Awassi sheep. However, more research is required to validate this mechanism.
Collapse
Affiliation(s)
- Ibrahim H Alwan
- Department of Animal Production, College of Agriculture, Al-Qasim Green University, Al-Qasim, Babil, 51001, Iraq
| | - Thamer R S Aljubouri
- Department of Animal Production, College of Agriculture, Al-Qasim Green University, Al-Qasim, Babil, 51001, Iraq
| | - Mohammed Baqur S Al-Shuhaib
- Department of Animal Production, College of Agriculture, Al-Qasim Green University, Al-Qasim, Babil, 51001, Iraq.
| |
Collapse
|
18
|
Roy AS, Feroz T, Islam MK, Munim MA, Supti DA, Antora NJ, Al Reza H, Gosh S, Bahadur NM, Alam MR, Hossain MS. A computational approach for structural and functional analyses of disease-associated mutations in the human CYLD gene. Genomics Inform 2024; 22:4. [PMID: 38907316 PMCID: PMC11184958 DOI: 10.1186/s44342-024-00007-2] [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: 11/22/2023] [Accepted: 12/26/2023] [Indexed: 06/23/2024] Open
Abstract
Tumor suppressor cylindromatosis protein (CYLD) regulates NF-κB and JNK signaling pathways by cleaving K63-linked poly-ubiquitin chain from its substrate molecules and thus preventing the progression of tumorigenesis and metastasis of the cancer cells. Mutations in CYLD can cause aberrant structure and abnormal functionality leading to tumor formation. In this study, we utilized several computational tools such as PANTHER, PROVEAN, PredictSNP, PolyPhen-2, PhD-SNP, PON-P2, and SIFT to find out deleterious nsSNPs. We also highlighted the damaging impact of those deleterious nsSNPs on the structure and function of the CYLD utilizing ConSurf, I-Mutant, SDM, Phyre2, HOPE, Swiss-PdbViewer, and Mutation 3D. We shortlisted 18 high-risk nsSNPs from a total of 446 nsSNPs recorded in the NCBI database. Based on the conservation profile, stability status, and structural impact analysis, we finalized 13 nsSNPs. Molecular docking analysis and molecular dynamic simulation concluded the study with the findings of two significant nsSNPs (R830K, H827R) which have a remarkable impact on binding affinity, RMSD, RMSF, radius of gyration, and hydrogen bond formation during CYLD-ubiquitin interaction. The principal component analysis compared native and two mutants R830K and H827R of CYLD that signify structural and energy profile fluctuations during molecular dynamic (MD) simulation. Finally, the protein-protein interaction network showed CYLD interacts with 20 proteins involved in several biological pathways that mutations can impair. Considering all these in silico analyses, our study recommended conducting large-scale association studies of nsSNPs of CYLD with cancer as well as designing precise medications against diseases associated with these polymorphisms.
Collapse
Affiliation(s)
- Arpita Singha Roy
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Tasmiah Feroz
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Md Kobirul Islam
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Md Adnan Munim
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Dilara Akhter Supti
- Department of Food Technology & Nutrition Sciences, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Nusrat Jahan Antora
- Department of Genetic Engineering and Biotechnology, East West University, Dhaka, 1212, Bangladesh
| | - Hasan Al Reza
- Department of Genetic Engineering and Biotechnology, University of Dhaka, Dhaka, 1000, Bangladesh
| | - Supriya Gosh
- Department of Food Technology & Nutrition Sciences, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Newaz Mohammed Bahadur
- Department of Applied Chemistry and Chemical Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh
| | - Mohammad Rahanur Alam
- Department of Food Technology & Nutrition Sciences, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh.
| | - Md Shahadat Hossain
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, 3814, Bangladesh.
| |
Collapse
|
19
|
Kamal MM, Mia MS, Faruque MO, Rabby MG, Islam MN, Talukder MEK, Wani TA, Rahman MA, Hasan MM. In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene. Sci Rep 2024; 14:11607. [PMID: 38773180 PMCID: PMC11109216 DOI: 10.1038/s41598-024-62299-2] [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: 01/16/2024] [Accepted: 05/15/2024] [Indexed: 05/23/2024] Open
Abstract
Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.
Collapse
Affiliation(s)
- Md Mostafa Kamal
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| | - Md Sohel Mia
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| | - Md Omar Faruque
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| | - Md Golam Rabby
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh
| | - Md Numan Islam
- Department of Food Engineering, North Pacific International University of Bangladesh, Dhaka, Bangladesh
| | | | - Tanveer A Wani
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, 11451, Riyadh, Saudi Arabia
| | - M Atikur Rahman
- Department of Biological Sciences, Alabama State University, 915 S Jackson St, Montgomery, AL, 36104, USA.
| | - Md Mahmudul Hasan
- Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, 7408, Bangladesh.
| |
Collapse
|
20
|
Hoda A, Bixheku X, Lika Çekani M. Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR geneComputational analysis of bovine PKLR gene. J Biomol Struct Dyn 2024; 42:4155-4168. [PMID: 37278385 DOI: 10.1080/07391102.2023.2219315] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 05/23/2023] [Indexed: 06/07/2023]
Abstract
Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins' stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anila Hoda
- Agricultural University of Tirana, Tirana, Albania
| | | | | |
Collapse
|
21
|
Hoda A, Berisha B, Bixheku X, Zanchi FB. Comprehensive in silico analysis of prolactin receptor (PRLR) gene nonsynonymous single nucleotide polymorphisms (nsSNPs) reveals multifaceted impact on protein structure, function, and interactions. J Biomol Struct Dyn 2024:1-17. [PMID: 38656135 DOI: 10.1080/07391102.2024.2335295] [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: 01/27/2024] [Accepted: 03/20/2024] [Indexed: 04/26/2024]
Abstract
This study delves into the functional and structural implications of non-synonymous single nucleotide polymorphisms (nsSNPs) within the Prolactin Receptor (PRLR) gene. Thirteen deleterious nsSNPs were identified through bioinformatics tools, with SIFT predicting 168 out of 395 nsSNPs as detrimental, exhibiting tolerance index (TI) scores ranging from 0 to 0.05. Polyphen2 assigned likelihood scores >0.99 to all 13 nsSNPs, indicating high probability of harm, while Panther scores classified most nsSNPs as 'probably damaging', with specific mutations like W218R scoring 0.74, suggesting a higher impact. Stability analysis using DDG I-Mutant and DDG Mupro consistently predicted decreased stability for all mutations, with CUPSAT indicating mutations like V125G and W218R significantly decreasing stability. Structural analysis through DynaMut predicted destabilization for all mutations except L196I and L292H. MutPred2 highlighted structural alterations for all nsSNPs except L196I, L293V, R315W, and S353N. Domain analysis revealed key mutations within essential functional domains, with five nsSNPs located within Fibronectin type-III domains. Bayesian analysis through ConSurf identified 9 critical residues, with 11 nsSNPs exhibiting notably high conservation. STRING analysis unveiled a complex interaction network, indicating involvement in vital biological processes like lactation. Molecular dynamics (MD) simulations, spanning 100 nanoseconds, elucidated structural dynamics induced by detrimental missense SNPs. Post-translational modification (PTM) analysis identified specific mutations, such as R351, involved in methylation, while S353 was implicated in phosphorylation and glycosylation. These findings offer comprehensive insights into the molecular and phenotypic effects of deleterious nsSNPs in the PRLR gene, crucial for selective breeding.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anila Hoda
- Department Animal Sciences, Agricultural University of Tirana, Albania
| | - Bajram Berisha
- Department of Animal Physiology & Immunology, Life Science Center Weihenstephan, TU Munich (TUM), Germany
| | | | - Fernando Berton Zanchi
- Pesquisador - Fiocruz Rondônia. Laboratório de Bioinformática e Química Medicinal-LABIOQUIM, Centro de Estudos de Biomoléculas Aplicadas à Saúde - CEBio
| |
Collapse
|
22
|
Huang CK, Yu MC, Hung CS, Lin JC. Emerging insight of whole genome sequencing coupled with protein structure prediction into the pyrazinamide-resistance signature of Mycobacterium tuberculosis. Int J Antimicrob Agents 2024; 63:107053. [PMID: 38081550 DOI: 10.1016/j.ijantimicag.2023.107053] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2023] [Revised: 11/11/2023] [Accepted: 12/04/2023] [Indexed: 02/25/2024]
Abstract
Pyrazinamide (PZA) is considered to be a pivotal drug to shorten the treatment of both drug-susceptible and drug-resistant tuberculosis, but its use is challenged by the reliability of drug-susceptibility testing (DST). PZA resistance in Mycobacterium tuberculosis (MTB) is relevant to the amino acid substitution of pyrazinamidase that is responsible for the conversion of PZA to active pyrazinoic acid (POA). The single nucleotide variants (SNVs) within ribosomal protein S1 (rpsA) or aspartate decarboxylase (panD), the binding targets of POA, has been reported to drive the PZA-resistance signature of MTB. In this study, whole genome sequencing (WGS) was used to identify SNVs within the pncA, rpsA and panD genes in 100 clinical MTB isolates associated with DST results for PZA. The potential influence of high-confidence, interim-confidence or emerging variants on the interplay between target genes and PZA or POA was simulated computationally, and predicted with a protein structure modelling approach. The DST results showed weak agreement with the identification of high-confidence variants within the pncA gene (Cohen's kappa coefficient=0.58), the analytic results of WGS coupled with protein structure modelling on pncA mutants (Cohen's kappa coefficient=0.524) or related genes (Cohen's kappa coefficient=0.504). Taken together, these results suggest the practicable application of a genotypic-coupled bioinformatic approach to manage PZA-containing regimens for patients with MTB.
Collapse
Affiliation(s)
- Chun-Kai Huang
- Ph.D. Program in Medical Biotechnology, College of Medical Science and Technology, Taipei Medical University, Taipei, Taiwan; Department of Laboratory Medicine, Wan Fang Hospital, Taipei Medical University, Taipei, Taiwan
| | - Ming-Chih Yu
- Division of Pulmonary Medicine, Department of Internal Medicine, Wan Fang Hospital, Taipei Medical University, Taipei, Taiwan; School of Respiratory Therapy, College of Medicine, Taipei Medical University, Taipei, Taiwan; Pulmonary Research Centre, Wan Fang Hospital, Taipei Medical University, Taipei, Taiwan
| | - Ching-Sheng Hung
- Department of Laboratory Medicine, Wan Fang Hospital, Taipei Medical University, Taipei, Taiwan
| | - Jung-Chun Lin
- Ph.D. Program in Medical Biotechnology, College of Medical Science and Technology, Taipei Medical University, Taipei, Taiwan; Pulmonary Research Centre, Wan Fang Hospital, Taipei Medical University, Taipei, Taiwan; School of Medical Laboratory Science and Biotechnology, College of Medical Science and Technology, Taipei Medical University, Taipei, Taiwan.
| |
Collapse
|
23
|
Panchal NK, Samdani P, Sengupta T, Prince SE. Computational Analysis of Non-synonymous SNPs in ATM Kinase: Structural Insights, Functional Implications, and Inhibitor Discovery. Mol Biotechnol 2024:10.1007/s12033-024-01120-x. [PMID: 38489015 DOI: 10.1007/s12033-024-01120-x] [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: 10/11/2023] [Accepted: 02/13/2024] [Indexed: 03/17/2024]
Abstract
Ataxia telangiectasia-mutated (ATM) protein kinase, a key player in cellular integrity regulation, is known for its role in DNA damage response. This study investigates the broader impact of ATM on cellular processes and potential clinical manifestations arising from mutations, aiming to expand our understanding of ATM's diverse functions beyond conventional roles. The research employs a comprehensive set of computational techniques for a thorough analysis of ATM mutations. The mutation data are curated from dbSNP and HuVarBase databases. A meticulous assessment is conducted, considering factors such as deleterious effects, protein stability, oncogenic potential, and biophysical characteristics of the identified mutations. Conservation analysis, utilizing diverse computational tools, provides insights into the evolutionary significance of these mutations. Molecular docking and dynamic simulation analyses are carried out for selected mutations, investigating their interactions with Y2080D, AZD0156, and quercetin inhibitors to gauge potential therapeutic implications. Among the 419 mutations scrutinized, five (V1913C, Y2080D, L2656P, C2770G, and C2930G) are identified as both disease causing and protein destabilizing. The study reveals the oncogenic potential of these mutations, supported by findings from the COSMIC database. Notably, Y2080D is associated with haematopoietic and lymphoid cancers, while C2770G shows a correlation with squamous cell carcinomas. Molecular docking and dynamic simulation analyses highlight strong binding affinities of quercetin for Y2080D and AZD0156 for C2770G, suggesting potential therapeutic options. In summary, this computational analysis provides a comprehensive understanding of ATM mutations, revealing their potential implications in cellular integrity and cancer development. The study underscores the significance of Y2080D and C2770G mutations, offering valuable insights for future precision medicine targeting-specific ATM. Despite informative computational analyses, a significant research gap exists, necessitating essential in vitro and in vivo studies to validate the predicted effects of ATM mutations on protein structure and function.
Collapse
Affiliation(s)
- Nagesh Kishan Panchal
- Department of Biotechnology, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, 632 014, India
| | - Poorva Samdani
- Department of Biomedical Sciences, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, India
| | - Tiasa Sengupta
- Department of Biomedical Sciences, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, India
| | - Sabina Evan Prince
- Department of Biotechnology, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, 632 014, India.
| |
Collapse
|
24
|
Kang J, Wei S, Jia Z, Ma Y, Chen H, Sun C, Xu J, Tao J, Dong Y, Lv W, Tian H, Guo X, Bi S, Zhang C, Jiang Y, Lv H, Zhang M. Effects of genetic variation on the structure of RNA and protein. Proteomics 2024; 24:e2300235. [PMID: 38197532 DOI: 10.1002/pmic.202300235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Revised: 12/15/2023] [Accepted: 12/19/2023] [Indexed: 01/11/2024]
Abstract
Changes in the structure of RNA and protein, have an important impact on biological functions and are even important determinants of disease pathogenesis and treatment. Some genetic variations, including copy number variation, single nucleotide variation, and so on, can lead to changes in biological function and increased susceptibility to certain diseases by changing the structure of RNA or protein. With the development of structural biology and sequencing technology, a large amount of RNA and protein structure data and genetic variation data resources has emerged to be used to explain biological processes. Here, we reviewed the effects of genetic variation on the structure of RNAs and proteins, and investigated their impact on several diseases. An online resource (http://www.onethird-lab.com/gems/) to support convenient retrieval of common tools is also built. Finally, the challenges and future development of the effects of genetic variation on RNA and protein were discussed.
Collapse
Affiliation(s)
- Jingxuan Kang
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Siyu Wei
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Zhe Jia
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Yingnan Ma
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Haiyan Chen
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Chen Sun
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Jing Xu
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Junxian Tao
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Yu Dong
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Wenhua Lv
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
| | - Hongsheng Tian
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
| | - Xuying Guo
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
| | - Shuo Bi
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
| | - Chen Zhang
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
| | - Yongshuai Jiang
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Hongchao Lv
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| | - Mingming Zhang
- College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, China
- The Epigenome-Wide Association Study Project, Harbin, China
| |
Collapse
|
25
|
Zhang M, Gong C, Ge F, Yu DJ. FCMSTrans: Accurate Prediction of Disease-Associated nsSNPs by Utilizing Multiscale Convolution and Deep Feature Combination within a Transformer Framework. J Chem Inf Model 2024; 64:1394-1406. [PMID: 38349747 DOI: 10.1021/acs.jcim.3c02025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/27/2024]
Abstract
Nonsynonymous single-nucleotide polymorphisms (nsSNPs), implicated in over 6000 diseases, necessitate accurate prediction for expedited drug discovery and improved disease diagnosis. In this study, we propose FCMSTrans, a novel nsSNP predictor that innovatively combines the transformer framework and multiscale modules for comprehensive feature extraction. The distinctive attribute of FCMSTrans resides in a deep feature combination strategy. This strategy amalgamates evolutionary-scale modeling (ESM) and ProtTrans (PT) features, providing an understanding of protein biochemical properties, and position-specific scoring matrix, secondary structure, predicted relative solvent accessibility, and predicted disorder (PSPP) features, which are derived from four protein sequences and structure-oriented characteristics. This feature combination offers a comprehensive view of the molecular dynamics involving nsSNPs. Our model employs the transformer's self-attention mechanisms across multiple layers, extracting higher-level and abstract representations. Simultaneously, varied-level features are captured by multiscale convolutions, enriching feature abstraction at multiple echelons. Our comparative analyses with existing methodologies highlight significant improvements made possible by the integrated feature fusion approach adopted in FCMSTrans. This is further substantiated by performance assessments based on diverse data sets, such as PredictSNP, MMP, and PMD, with areas under the curve (AUCs) of 0.869, 0.819, and 0.693, respectively. Furthermore, FCMSTrans shows robustness and superiority by outperforming the current best predictor, PROVEAN, in a blind test conducted on a third-party data set, achieving an impressive AUC score of 0.7838. The Python code of FCMSTrans is available at https://github.com/gc212/FCMSTrans for academic usage.
Collapse
Affiliation(s)
- Ming Zhang
- School of Computer, Jiangsu University of Science and Technology, 666 Changhui Road, Zhenjiang 212100, China
| | - Chao Gong
- School of Computer, Jiangsu University of Science and Technology, 666 Changhui Road, Zhenjiang 212100, China
| | - Fang Ge
- State Key Laboratory of Organic Electronics and Information Displays & Institute of Advanced Materials (IAM), Nanjing University of Posts & Telecommunications, 9 Wenyuan Road, Nanjing 210023, China
| | - Dong-Jun Yu
- School of Computer Science and Engineering, Nanjing University of Science and Technology, 200 Xiaolingwei, Nanjing 210094, China
| |
Collapse
|
26
|
Gupta A, Jamal T, Rajbhar P, Gaur AS, Chauhan SS, Parthasarathi R. Cytokines inhibitory mechanism of Prunus domestica L. (Plum) peptides as potential immunomodulators against systemic lupus erythematosus: an in-silico screening. In Silico Pharmacol 2024; 12:12. [PMID: 38370860 PMCID: PMC10866836 DOI: 10.1007/s40203-023-00188-8] [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: 08/30/2023] [Accepted: 12/28/2023] [Indexed: 02/20/2024] Open
Abstract
Natural bioactive peptides exhibit various chemical and structural properties to enhance the immune response against multiple inflammatory and autoimmune related disorders. The immunomodulatory function and bioactivity of seed peptides show the capability for the development of biotherapeutics that could prevent autoimmune diseases. The aim of current study is to determine the immunomodulatory function of bioactive peptides derived from the seed of plum (Prunus domestica L.) by applying various immunoinformatic approaches. A thorough analysis of forty-one peptides was performed including drug likeliness, pharmacokinetic, and bioactivity profiling studies. Further, molecular docking and molecular dynamics (MD) simulations of screened peptides were carried out with the two interleukin targets (IL-17A and IL-23) of systemic lupus erythematosus (SLE). After the systematic screening, four peptides, namely HLLP, LPLL, LPAGV, and NLPL, were found as potential inhibitors against SLE. Additionally, site-directed mutagenesis analysis was conducted to explore the role of essential amino acid residues in the binding pattern/energy change. Computational alanine screening analysis found that CYS123, CYS121 of IL-17A and ASP270, and SER249 of IL-23 as hot spot residues that could play an important role in the inhibition property of screened peptides. Overall, the methodology described in the study can be utilized for developing unique peptide inhibitors that have a preventative role against SLE. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1007/s40203-023-00188-8.
Collapse
Affiliation(s)
- Anshika Gupta
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
| | - Tanya Jamal
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
| | - Priyanka Rajbhar
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
| | - Anamika Singh Gaur
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
| | - Shweta Singh Chauhan
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
- Academy of Scientific and Innovative Research (AcSIR), Uttar Pradesh, Ghaziabad, 201002 India
| | - Ramakrishnan Parthasarathi
- Computational Toxicology Facility, Toxicoinformatics Research Group, CSIR-Indian Institute of Toxicology Research, Vishvigyan Bhavan 31, Mahatma Gandhi Marg, Lucknow, Uttar Pradesh 226001 India
- Academy of Scientific and Innovative Research (AcSIR), Uttar Pradesh, Ghaziabad, 201002 India
| |
Collapse
|
27
|
Hristova SH, Zhivkov AM. Three-Dimensional Structural Stability and Local Electrostatic Potential at Point Mutations in Spike Protein of SARS-CoV-2 Coronavirus. Int J Mol Sci 2024; 25:2174. [PMID: 38396850 PMCID: PMC10889838 DOI: 10.3390/ijms25042174] [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: 01/17/2024] [Revised: 01/30/2024] [Accepted: 02/01/2024] [Indexed: 02/25/2024] Open
Abstract
The contagiousness of SARS-CoV-2 β-coronavirus is determined by the virus-receptor electrostatic association of its positively charged spike (S) protein with the negatively charged angiotensin converting enzyme-2 (ACE2 receptor) of the epithelial cells. If some mutations occur, the electrostatic potential on the surface of the receptor-binding domain (RBD) could be altered, and the S-ACE2 association could become stronger or weaker. The aim of the current research is to investigate whether point mutations can noticeably alter the electrostatic potential on the RBD and the 3D stability of the S1-subunit of the S-protein. For this purpose, 15 mutants with different hydrophilicity and electric charge (positive, negative, or uncharged) of the substituted and substituting amino acid residues, located on the RBD at the S1-ACE2 interface, are selected, and the 3D structure of the S1-subunit is reconstructed on the base of the crystallographic structure of the S-protein of the wild-type strain and the amino acid sequence of the unfolded polypeptide chain of the mutants. Then, the Gibbs free energy of folding, isoelectric point, and pH-dependent surface electrostatic potential of the S1-subunit are computed using programs for protein electrostatics. The results show alterations in the local electrostatic potential in the vicinity of the mutant amino acid residue, which can influence the S-ACE2 association. This approach allows prediction of the relative infectivity, transmissibility, and contagiousness (at equal social immune status) of new SARS-CoV-2 mutants by reconstruction of the 3D structure of the S1-subunit and calculation of the surface electrostatic potential.
Collapse
Affiliation(s)
- Svetlana H. Hristova
- Department of Medical Physics and Biophysics, Medical Faculty, Medical University—Sofia, Zdrave Street 2, 1431 Sofia, Bulgaria;
| | - Alexandar M. Zhivkov
- Scientific Research Center, “St. Kliment Ohridski” Sofia University, 8 Dragan Tsankov Blvd., 1164 Sofia, Bulgaria
| |
Collapse
|
28
|
Brandt E, Harjama L, Elomaa O, Saarela J, Donner K, Lappalainen K, Kivirikko S, Ranki A, Kere J, Kettunen K, Hannula-Jouppi K. A novel SERPINA12 variant and first European patients with diffuse palmoplantar keratoderma. J Eur Acad Dermatol Venereol 2024; 38:413-418. [PMID: 37684051 DOI: 10.1111/jdv.19498] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2023] [Accepted: 08/18/2023] [Indexed: 09/10/2023]
Abstract
BACKGROUND Hereditary palmoplantar keratodermas (hPPKs) comprise a heterogeneous group of skin disorders characterized by persistent palmoplantar hyperkeratosis. Loss-of-function variants in a serine peptidase inhibitor, SERPINA12, have recently been implicated in autosomal recessive diffuse hPPK. The disorder appears to share similarities with another hPPK associated with protease overactivity, namely Nagashima-type PPK (NPPK) caused by biallelic variants in SERPINB7. OBJECTIVES The aim of this study was to enhance the understanding of the clinical and genetic characteristics of serine protease-related hPPKs caused by variants in SERPINA12 and SERPINB7. METHODS Whole-exome sequencing (WES) was performed for hPPK patients. Haplotype analysis was completed for the patients with identified recessive SERPINA12 variants and their available family members. In addition, the current literature of SERPINA12- and SERPINB7-related hPPKs was summarized. RESULTS The phenotype of SERPINA12-related hPPK was confirmed by reporting three new SERPINA12 patients, the first of European origin. A novel SERPINA12 c.1100G>A p.(Gly367Glu) missense variant was identified confirming that the variant spectrum of SERPINA12 include both truncating and missense variants. The previously reported SERPINA12 c.631C>T p.(Arg211*) was indicated enriched in the Finnish population due to a plausible founder effect. In addition, SERPINA12 hPPK patients were shown to share a similar phenotype to patients with recessive variants in SERPINB7. The shared phenotype included diffuse transgradient PPK since birth or early childhood and frequent palmoplantar hyperhidrosis, aquagenic whitening and additional hyperkeratotic lesions in non-palmoplantar areas. SERPINA12 and SERPINB7 hPPK patients cannot be distinguished without genetic analysis. CONCLUSIONS Recessive variants in SERPINA12 and SERPINB7 leading to protease overactivity and hPPK produce a similar phenotype, indistinguishable without genetic analysis. SERPINA12 variants should be assessed also in non-Asian patients with diffuse transgradient PPK. Understanding the role of serine protease inhibitors will provide insights into the complex proteolytic network in epidermal homeostasis.
Collapse
Affiliation(s)
- E Brandt
- Department of Dermatology and Allergology, ERN-Skin Center, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
| | - L Harjama
- Department of Dermatology and Allergology, ERN-Skin Center, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
| | - O Elomaa
- Folkhälsan Research Center, Helsinki, Finland and Research Programs Unit, Stem Cells and Metabolism Research Program, University of Helsinki, Helsinki, Finland
| | - J Saarela
- Institute for Molecular Medicine Finland (FIMM), Helsinki Institute of Life Science (HiLIFE), University of Helsinki, Helsinki, Finland
- The Centre for Molecular Medicine Norway (NCMM), University of Oslo, Oslo, Norway
| | - K Donner
- Institute for Molecular Medicine Finland (FIMM), Helsinki Institute of Life Science (HiLIFE), University of Helsinki, Helsinki, Finland
| | - K Lappalainen
- Department of Dermatology and Allergology, ERN-Skin Center, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
| | - S Kivirikko
- Department of Clinical Genetics and Department of Medical and Clinical Genetics, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
| | - A Ranki
- Department of Dermatology and Allergology, ERN-Skin Center, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
| | - J Kere
- Folkhälsan Research Center, Helsinki, Finland and Research Programs Unit, Stem Cells and Metabolism Research Program, University of Helsinki, Helsinki, Finland
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| | - K Kettunen
- Institute for Molecular Medicine Finland (FIMM), Helsinki Institute of Life Science (HiLIFE), University of Helsinki, Helsinki, Finland
- HUS Diagnostic Center, Division of Genetics and Clinical Pharmacology, Laboratory of Genetics, University of Helsinki and Helsinki University Hospital, Helsinki, Finland
| | - K Hannula-Jouppi
- Department of Dermatology and Allergology, ERN-Skin Center, University of Helsinki and Helsinki University Central Hospital, Helsinki, Finland
- Folkhälsan Research Center, Helsinki, Finland and Research Programs Unit, Stem Cells and Metabolism Research Program, University of Helsinki, Helsinki, Finland
| |
Collapse
|
29
|
Kumar R, Jayaraman M, Ramadas K, Chandrasekaran A. Computational identification and analysis of deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene: a structural and functional impact. J Biomol Struct Dyn 2024; 42:1518-1532. [PMID: 37173831 DOI: 10.1080/07391102.2023.2211674] [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/02/2023] [Accepted: 04/02/2023] [Indexed: 05/15/2023]
Abstract
Cytochrome P450 oxidoreductase (POR) protein is essential for steroidogenesis, and POR gene mutations are frequently associated with P450 Oxidoreductase Deficiency (PORD), a disorder of hormone production. To our knowledge, no previous attempt has been made to identify and analyze the deleterious/pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene through an extensive computational approach. Computational algorithms and tools were employed to identify, characterize, and validate the pathogenic SNPs associated with certain diseases. To begin with, all the high-confidence SNPs were collected, and their structural and functional impacts on the protein structures were explored. The results of various in silico analyses affirm that the A287P and R457H variants of POR could destabilize the interactions between the amino acids and the hydrogen bond networks, resulting in functional deviations of POR. The literature study further confirms that the pathogenic mutations (A287P and R457H) are associated with the onset of PORD. Molecular dynamics simulations (MDS) and essential dynamics (ED) studies characterized the structural consequences of prioritized deleterious mutations, representing the structural destabilization that might disrupt POR biological function. The identified deleterious mutations at the cofactor's binding domains might interfere with the essential interactions between the protein and cofactors, thus inhibiting POR catalytic activity. The consolidated insights from the computational analyses can be used to predict potential deleterious mutants and understand the disease's pathological basis and the molecular mechanism of drug metabolism for the application of personalized medication. HIGHLIGHTSNADPH cytochrome P450 oxidoreductase (POR) mutations are associated with a broad spectrum of human diseasesIdentified and analyzed the most deleterious nsSNPs of POR through the sequence and structure-based prediction toolsInvestigated the structural and functional impacts of the most significant mutations (A287P and R457H) associated with PORDMolecular dynamics and PCA-based FEL analysis were utilized to probe the mutation-induced structural alterations in PORCommunicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Rajalakshmi Kumar
- Central Inter-Disciplinary Research Facility, Sri Balaji Vidyapeeth (Deemed to be University), Pillayarkuppam, Puducherry, India
| | - Manikandan Jayaraman
- Department of Bioinformatics, School of Life Sciences, Pondicherry University, Kalapet, Puducherry, India
| | - Krishna Ramadas
- Department of Bioinformatics, School of Life Sciences, Pondicherry University, Kalapet, Puducherry, India
| | - Adithan Chandrasekaran
- Central Inter-Disciplinary Research Facility, Sri Balaji Vidyapeeth (Deemed to be University), Pillayarkuppam, Puducherry, India
| |
Collapse
|
30
|
Martinez SE, Pandey AV, Perez Jimenez TE, Zhu Z, Court MH. Pharmacogenomics of poor drug metabolism in greyhounds: Canine P450 oxidoreductase genetic variation, breed heterogeneity, and functional characterization. PLoS One 2024; 19:e0297191. [PMID: 38300925 PMCID: PMC10833530 DOI: 10.1371/journal.pone.0297191] [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: 08/05/2023] [Accepted: 12/31/2023] [Indexed: 02/03/2024] Open
Abstract
Greyhounds metabolize cytochrome P450 (CYP) 2B11 substrates more slowly than other dog breeds. However, CYP2B11 gene variants associated with decreased CYP2B11 expression do not fully explain reduced CYP2B11 activity in this breed. P450 oxidoreductase (POR) is an essential redox partner for all CYPs. POR protein variants can enhance or repress CYP enzyme function in a CYP isoform and substrate dependent manner. The study objectives were to identify POR protein variants in greyhounds and determine their effect on coexpressed CYP2B11 and CYP2D15 enzyme function. Gene sequencing identified two missense variants (Glu315Gln and Asp570Glu) forming four alleles, POR-H1 (reference), POR-H2 (570Glu), POR-H3 (315Gln, 570Glu) and POR-H4 (315Gln). Out of 68 dog breeds surveyed, POR-H2 was widely distributed across multiple breeds, while POR-H3 was largely restricted to greyhounds and Scottish deerhounds (35% allele frequencies), and POR-H4 was rare. Three-dimensional protein structure modelling indicated significant effects of Glu315Gln (but not Asp570Glu) on protein flexibility through loss of a salt bridge between Glu315 and Arg519. Recombinant POR-H1 (reference) and each POR variant (H2-H4) were expressed alone or with CYP2B11 or CYP2D15 in insect cells. No substantial effects on POR protein expression or enzyme activity (cytochrome c reduction) were observed for any POR variant (versus POR-H1) when expressed alone or with CYP2B11 or CYP2D15. Furthermore, there were no effects on CYP2B11 or CYP2D15 protein expression, or on CYP2D15 enzyme kinetics by any POR variant (versus POR-H1). However, Vmax values for 7-benzyloxyresorufin, propofol and bupropion oxidation by CYP2B11 were significantly reduced by coexpression with POR-H3 (by 34-37%) and POR-H4 (by 65-72%) compared with POR-H1. Km values were unaffected. Our results indicate that the Glu315Gln mutation (common to POR-H3 and POR-H4) reduces CYP2B11 enzyme function without affecting at least one other major canine hepatic P450 (CYP2D15). Additional in vivo studies are warranted to confirm these findings.
Collapse
Affiliation(s)
- Stephanie E. Martinez
- Pharmacogenomics Laboratory, Program in Individualized Medicine (PrIMe), Department of Veterinary Clinical Sciences, College of Veterinary Medicine, Washington State University, Pullman, Washington, United States of America
| | - Amit V. Pandey
- Pediatric Endocrinology, Diabetology, and Metabolism, Department of Biomedical Research, University Children’s Hospital Bern, Switzerland and Translational Hormone Research Program, University of Bern, Bern, Switzerland
| | - Tania E. Perez Jimenez
- Pharmacogenomics Laboratory, Program in Individualized Medicine (PrIMe), Department of Veterinary Clinical Sciences, College of Veterinary Medicine, Washington State University, Pullman, Washington, United States of America
| | - Zhaohui Zhu
- Pharmacogenomics Laboratory, Program in Individualized Medicine (PrIMe), Department of Veterinary Clinical Sciences, College of Veterinary Medicine, Washington State University, Pullman, Washington, United States of America
| | - Michael H. Court
- Pharmacogenomics Laboratory, Program in Individualized Medicine (PrIMe), Department of Veterinary Clinical Sciences, College of Veterinary Medicine, Washington State University, Pullman, Washington, United States of America
| |
Collapse
|
31
|
Rodrigues CHM, Portelli S, Ascher DB. Exploring the effects of missense mutations on protein thermodynamics through structure-based approaches: findings from the CAGI6 challenges. Hum Genet 2024:10.1007/s00439-023-02623-4. [PMID: 38227011 DOI: 10.1007/s00439-023-02623-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Accepted: 11/18/2023] [Indexed: 01/17/2024]
Abstract
Missense mutations are known contributors to diverse genetic disorders, due to their subtle, single amino acid changes imparted on the resultant protein. Because of this, understanding the impact of these mutations on protein stability and function is crucial for unravelling disease mechanisms and developing targeted therapies. The Critical Assessment of Genome Interpretation (CAGI) provides a valuable platform for benchmarking state-of-the-art computational methods in predicting the impact of disease-related mutations on protein thermodynamics. Here we report the performance of our comprehensive platform of structure-based computational approaches to evaluate mutations impacting protein structure and function on 3 challenges from CAGI6: Calmodulin, MAPK1 and MAPK3. Our stability predictors have achieved correlations of up to 0.74 and AUCs of 1 when predicting changes in ΔΔG for MAPK1 and MAPK3, respectively, and AUC of up to 0.75 in the Calmodulin challenge. Overall, our study highlights the importance of structure-based approaches in understanding the effects of missense mutations on protein thermodynamics. The results obtained from the CAGI6 challenges contribute to the ongoing efforts to enhance our understanding of disease mechanisms and facilitate the development of personalised medicine approaches.
Collapse
Affiliation(s)
- Carlos H M Rodrigues
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, VIC, 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, QLD, 4072, Australia
| | - Stephanie Portelli
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, VIC, 3004, Australia
- School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, QLD, 4072, Australia
| | - David B Ascher
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, VIC, 3004, Australia.
- School of Chemistry and Molecular Biosciences, University of Queensland, St Lucia, QLD, 4072, Australia.
| |
Collapse
|
32
|
Liu B, Jiang Y, Yang Y, Chen JX. OmeDDG: Improved Protein Mutation Stability Prediction Based on Predicted 3D Structures. J Phys Chem B 2024; 128:67-76. [PMID: 38130113 DOI: 10.1021/acs.jpcb.3c05601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Determining changes in the protein's thermal stability following mutations is critical in protein engineering and understanding pathogenic missense mutations. Despite the development of various computational methods to predict the effects of single-point mutations, their accuracy remains limited. In this study, we propose a new computational method, OmeDDG, that more accurately predicts mutation-induced Gibbs free energy changes in protein folding (ΔΔG). OmeDDG takes the sequences of wild-type and mutant proteins as input, utilizes OmegaFold to obtain the 3D structure, employs a convolutional neural network to extract structural features, and combines them with protein mutation features and pretraining features to predict the stability of single-point mutations in proteins. We performed a comprehensive comparison between OmeDDG and other available prediction methods on four blind test datasets, confirming that OmeDDG can effectively enhance protein mutation prediction performance. Notably, on the antisymmetric dataset Ssym, OmeDDG achieves the best performance, demonstrating favorable antisymmetry with PCC = 0.79 and RMSE = 0.96 for forward mutations and PCC = 0.77 and RMSE = 0.97 for reverse mutant types.
Collapse
Affiliation(s)
- Baoying Liu
- School of Computing and Artificial Intelligence, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
| | - Yongquan Jiang
- School of Computing and Artificial Intelligence, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
- Artificial Intelligence Research Institute, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
| | - Yan Yang
- School of Computing and Artificial Intelligence, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
- Artificial Intelligence Research Institute, Southwest Jiaotong University, Chengdu 611756, Sichuan, China
| | - Jim X Chen
- Department of Computer Science, George Mason University, Fairfax, Virginia 22030-4444, United States
| |
Collapse
|
33
|
Annan A, Raiss N, Lemrabet S, Elomari N, Elmir EH, Filali-Maltouf A, Medraoui L, Oumzil H. Proposal of pharmacophore model for HIV reverse transcriptase inhibitors: Combined mutational effect analysis, molecular dynamics, molecular docking and pharmacophore modeling study. Int J Immunopathol Pharmacol 2024; 38:3946320241231465. [PMID: 38296818 PMCID: PMC10832406 DOI: 10.1177/03946320241231465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Accepted: 01/13/2024] [Indexed: 02/02/2024] Open
Abstract
OBJECTIVES Antiretroviral therapy (ART) efficacy is jeopardized by the emergence of drug resistance mutations in HIV, compromising treatment effectiveness. This study aims to propose novel analogs of Effavirenz (EFV) as potential direct inhibitors of HIV reverse transcriptase, employing computer-aided drug design methodologies. METHODS Three key approaches were applied: a mutational profile study, molecular dynamics simulations, and pharmacophore development. The impact of mutations on the stability, flexibility, function, and affinity of target proteins, especially those associated with NRTI, was assessed. Molecular dynamics analysis identified G190E as a mutation significantly altering protein properties, potentially leading to therapeutic failure. Comparative analysis revealed that among six first-line antiretroviral drugs, EFV exhibited notably low affinity with viral reverse transcriptase, further reduced by the G190E mutation. Subsequently, a search for EFV-similar inhibitors yielded 12 promising molecules based on their affinity, forming the basis for generating a pharmacophore model. RESULTS Mutational analysis pinpointed G190E as a crucial mutation impacting protein properties, potentially undermining therapeutic efficacy. EFV demonstrated diminished affinity with viral reverse transcriptase, exacerbated by the G190E mutation. The search for EFV analogs identified 12 high-affinity molecules, culminating in a pharmacophore model elucidating key structural features crucial for potent inhibition. CONCLUSION This study underscores the significance of EFV analogs as potential inhibitors of HIV reverse transcriptase. The findings highlight the impact of mutations on drug efficacy, particularly the detrimental effect of G190E. The generated pharmacophore model serves as a pivotal reference for future drug development efforts targeting HIV, providing essential structural insights for the design of potent inhibitors based on EFV analogs identified in vitro.
Collapse
Affiliation(s)
- Azzeddine Annan
- Research Center of Plant and Microbial Biotechnologies, Biodiversity and Environment, Faculty of Sciences, Mohammed V University, Rabat, Morocco
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
| | - Noureddine Raiss
- Research Center of Plant and Microbial Biotechnologies, Biodiversity and Environment, Faculty of Sciences, Mohammed V University, Rabat, Morocco
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
| | - Sanae Lemrabet
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
| | - Nezha Elomari
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
| | - El Harti Elmir
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
| | - Abdelkarim Filali-Maltouf
- Research Center of Plant and Microbial Biotechnologies, Biodiversity and Environment, Faculty of Sciences, Mohammed V University, Rabat, Morocco
| | - Leila Medraoui
- Research Center of Plant and Microbial Biotechnologies, Biodiversity and Environment, Faculty of Sciences, Mohammed V University, Rabat, Morocco
| | - Hicham Oumzil
- Virology Department, National Reference Laboratory for HIV, Institute National of Hygiene, Rabat, Morocco
- Pedagogy and Research Unit of Microbiology, and Genomic Center of Human Pathologies, School of Medicine and Pharmacy, Mohamed V University, Rabat, Morocco
| |
Collapse
|
34
|
Panchal NK, Mohanty S, Prince SE. Computational insights into NIMA-related kinase 6: unraveling mutational effects on structure and function. Mol Cell Biochem 2023:10.1007/s11010-023-04910-0. [PMID: 38117419 DOI: 10.1007/s11010-023-04910-0] [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: 09/09/2023] [Accepted: 11/25/2023] [Indexed: 12/21/2023]
Abstract
The NEK6 (NIMA-related kinase 6) serine/threonine kinase is a pivotal player in a multitude of cellular processes, including the regulation of the cell cycle and the response to DNA damage. Its significance extends to disease pathogenesis, as changes in NEK6 activity have been linked to the development of cancer. Non-synonymous single nucleotide polymorphisms (nsSNPs) in NEK6 have been linked to cancer as they alter the protein's native structure and function. The association between NEK6 activity and cancer development has prompted researchers to explore the effects of genetic variations within the NEK6 gene. Therefore, we utilized advanced computational tools to analyze 155 high-confidence nsSNPs in the NEK6 gene. From this analysis, 21 nsSNPs were identified as potentially harmful, raising concerns about their impact on NEK6 activity and cancer risk. These 21 mutations were then examined for structural alterations, and eight of nsSNPs (I51M, V76A, I134N, Y152D, R171Q, V186G, L237R, and C285S) were found to destabilize the protein. Among the destabilizing mutations screened, a specific mutation, R171Q, stood out due to its conserved nature. To understand its impact on the protein and conformation, all-atom molecular dynamics simulations (MDS) for 100 ns were performed for both Wildtype NEK6 (WT-NEK6) and R171Q. The simulations revealed that the R171Q variant was unstable and led to significant conformational changes in NEK6. This study provides valuable insights into NEK6 dysfunction caused by single amino acid alterations, offering a novel understanding of the molecular mechanisms underlying NEK6-related cancer progression.
Collapse
Affiliation(s)
- Nagesh Kishan Panchal
- Department of Biotechnology, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, 632 014, India
| | - Shruti Mohanty
- Department of Biomedical Sciences, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, 632014, India
| | - Sabina Evan Prince
- Department of Biotechnology, School of Biosciences and Technology, Vellore Institute of Technology, Vellore, 632 014, India.
| |
Collapse
|
35
|
Kolhe JA, Babu NL, Freeman BC. Protocol for establishing a protein interactome based on close physical proximity to a target protein within live budding yeast. STAR Protoc 2023; 4:102663. [PMID: 37883222 PMCID: PMC10630676 DOI: 10.1016/j.xpro.2023.102663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Revised: 06/28/2023] [Accepted: 10/03/2023] [Indexed: 10/28/2023] Open
Abstract
Here, we present a protocol for establishing a protein interactome based on close physical proximity to a target protein within live yeast cells. We describe steps for capturing both transient and stable binders by integrating a non-natural amino acid. We detail procedures for employing a site-directed method for labeling the surface that mediates protein associations and uncovers the binding sites on the interactors. Combined with mass spectrometry, our approach proves valuable in discovering binding partners and constructing a comprehensive protein-interaction network.
Collapse
Affiliation(s)
- Janhavi A Kolhe
- Department of Cell and Developmental Biology, School of Molecular and Cellular Biology, University of Illinois-Urbana-Champaign, Urbana, IL, USA
| | - Neethu L Babu
- Department of Cell and Developmental Biology, School of Molecular and Cellular Biology, University of Illinois-Urbana-Champaign, Urbana, IL, USA
| | - Brian C Freeman
- Department of Cell and Developmental Biology, School of Molecular and Cellular Biology, University of Illinois-Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
36
|
Forte G, Buonadonna AL, Pantaleo A, Fasano C, Capodiferro D, Grossi V, Sanese P, Cariola F, De Marco K, Lepore Signorile M, Manghisi A, Guglielmi AF, Simonetti S, Laforgia N, Disciglio V, Simone C. Classic Galactosemia: Clinical and Computational Characterization of a Novel GALT Missense Variant (p.A303D) and a Literature Review. Int J Mol Sci 2023; 24:17388. [PMID: 38139222 PMCID: PMC10744227 DOI: 10.3390/ijms242417388] [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: 10/20/2023] [Revised: 11/30/2023] [Accepted: 12/08/2023] [Indexed: 12/24/2023] Open
Abstract
Classic galactosemia is an autosomal recessive inherited liver disorder of carbohydrate metabolism caused by deficient activity of galactose-1-phosphate uridylyltransferase (GALT). While a galactose-restricted diet is lifesaving, most patients still develop long-term complications. In this study, we report on a two-week-old female patient who is a compound heterozygote for a known pathogenic variant (p.K285N) and a novel missense variant (p.A303D) in the GALT gene. Segregation analysis showed that the patient inherited the p.K285N pathogenic variant from her father and the p.A303D variant from her mother. A bioinformatics analysis to predict the impact of the p.A303D missense variant on the structure and stability of the GALT protein revealed that it may be pathogenic. Based on this finding, we performed a literature review of all GALT missense variants identified in homozygous and compound heterozygous galactosemia patients carrying the p.K285N pathogenic variant to explore their molecular effects on the clinical phenotype of the disease. Our analysis revealed that these missense variants are responsible for a wide range of molecular defects. This study expands the clinical and mutational spectrum in classic galactosemia and reinforces the importance of understanding the molecular consequences of genetic variants to incorporate genetic analysis into clinical care.
Collapse
Affiliation(s)
- Giovanna Forte
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Antonia Lucia Buonadonna
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Antonino Pantaleo
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Candida Fasano
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Donatella Capodiferro
- Section of Neonatology and Neonatal Intensive Care Unit, Department of Interdisciplinary Medicine, “Aldo Moro” University of Bari, 70121 Bari, Italy; (D.C.); (N.L.)
| | - Valentina Grossi
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Paola Sanese
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Filomena Cariola
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Katia De Marco
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Martina Lepore Signorile
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Andrea Manghisi
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Anna Filomena Guglielmi
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Simonetta Simonetti
- Clinical Pathology and Neonatal Screening, Azienda Ospedaliera Universitaria Policlinico-Giovanni XXIII, 70124 Bari, Italy;
| | - Nicola Laforgia
- Section of Neonatology and Neonatal Intensive Care Unit, Department of Interdisciplinary Medicine, “Aldo Moro” University of Bari, 70121 Bari, Italy; (D.C.); (N.L.)
| | - Vittoria Disciglio
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
| | - Cristiano Simone
- Medical Genetics, National Institute of Gastroenterology-IRCCS “Saverio de Bellis” Research Hospital, 70013 Castellana Grotte, Italy; (G.F.); (A.L.B.); (A.P.); (C.F.); (V.G.); (P.S.); (F.C.); (K.D.M.); (M.L.S.); (A.M.); (A.F.G.)
- Medical Genetics, Department of Precision and Regenerative Medicine and Jonic Area (DiMePRe-J), University of Bari Aldo Moro, 70124 Bari, Italy
| |
Collapse
|
37
|
Aljubouri TRS, Al-Shuhaib MBS. A missense SNP in the proopiomelanocortin (POMC) gene is associated with growth traits in Awassi and Karakul sheep. Anim Biotechnol 2023; 34:4837-4850. [PMID: 37071507 DOI: 10.1080/10495398.2023.2197469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2023]
Abstract
This study was conducted to assess the association between proopiomelanocortin (POMC) gene and growth traits in Awassi and Karakul sheep. PCR-single strand conformation polymorphism (SSCP) method was utilized to assess the polymorphism of POMC PCR amplicons with body weight and length, wither and rump height, chest and abdominal circumference measured at birth, 3, 6, 9, and 12 months intervals. Only one missense SNP (rs424417456:C > A) was detected in exon-2, in which glycine was converted to cysteine in the 65th position in POMC (p.65Gly > Cys). rs424417456 SNP showed significant associations with all growth traits in the third, sixth, ninth, and twelfth months. At the age of 3 months onward, lambs with CC genotype showed higher body weight, body length, wither and rump heights, and chest and abdominal circumferences than lambs with CA and AA genotypes, respectively. Prediction analyses indicated a deleterious effect of p.65Gly > Cys on POMC structure, function, and stability. Owing to the strong correlation between rs424417456:CC and better growth characteristics, this genotype is proposed as a promising marker to enhance growth traits in Awassi and Karakul sheep. The predicted damaging effects caused by rs424417456:CA and rs424417456:AA genotypes may entail a putative mechanism through which lambs with these genotypes exhibit lower growth traits.
Collapse
Affiliation(s)
- Thamer R S Aljubouri
- Department of Animal Production, College of Agriculture, Al-Qasim Green University, Babil, Iraq
| | | |
Collapse
|
38
|
Azmi MB, Jawed A, Ahmed SDH, Naeem U, Feroz N, Saleem A, Sardar K, Qureshi SA, Azim MK. Understanding the impact of structural modifications at the NNAT gene's post-translational acetylation site: in silico approach for predicting its drug-interaction role in anorexia nervosa. Eat Weight Disord 2023; 28:97. [PMID: 37987927 PMCID: PMC10663277 DOI: 10.1007/s40519-023-01618-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Accepted: 10/18/2023] [Indexed: 11/22/2023] Open
Abstract
PURPOSE Anorexia nervosa (AN) is a neuropsychological public health concern with a socially disabling routine and affects a person's healthy relationship with food. The role of the NNAT (Neuronatin) gene in AN is well established. The impact of mutation at the protein's post-translational modification (PTM) site has been exclusively associated with the worsening of the protein's biochemical dynamics. METHODS To understand the relationship between genotype and phenotype, it is essential to investigate the appropriate molecular stability of protein required for proper biological functioning. In this regard, we investigated the PTM-acetylation site of the NNAT gene in terms of 19 other specific amino acid probabilities in place of wild type (WT) through various in silico algorithms. Based on the highest pathogenic impact computed through the consensus classifier tool, we generated 3 residue-specific (K59D, P, W) structurally modified 3D models of NNAT. These models were further tested through the AutoDock Vina tool to compute the molecular drug binding affinities and inhibition constant (Ki) of structural variants and WT 3D models. RESULTS With trained in silico machine learning algorithms and consensus classifier; the three structural modifications (K59D, P, W), which were also the most deleterious substitution at the acetylation site of the NNAT gene, showed the highest structural destabilization and decreased molecular flexibility. The validation and quality assessment of the 3D model of these structural modifications and WT were performed. They were further docked with drugs used to manage AN, it was found that the ΔGbind (kcal/mol) values and the inhibition constants (Ki) were relatively lower in structurally modified models as compared to WT. CONCLUSION We concluded that any future structural variation(s) at the PTM-acetylation site of the NNAT gene due to possible mutational consequences, will serve as a basis to explore its relationship with the propensity of developing AN. LEVEL OF EVIDENCE No level of evidence-open access bioinformatics research.
Collapse
Affiliation(s)
- Muhammad Bilal Azmi
- Department of Biochemistry, Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan.
| | - Areesha Jawed
- Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan
| | - Syed Danish Haseen Ahmed
- Department of Biochemistry, Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan
| | - Unaiza Naeem
- Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan
| | - Nazia Feroz
- Department of Biochemistry, Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan
| | - Arisha Saleem
- Dow Medical College, Dow University of Health Sciences, Karachi, Pakistan
| | - Kainat Sardar
- Department of Biochemistry, University of Karachi, Karachi, Pakistan
- Department of Chemistry, Bahria College NORE-1, Karachi, Pakistan
| | | | - M Kamran Azim
- Department of Biosciences, Mohammad Ali Jinnah University, Karachi, Pakistan
| |
Collapse
|
39
|
Kurniawan J, Ishida T. Comparing Supervised Learning and Rigorous Approach for Predicting Protein Stability upon Point Mutations in Difficult Targets. J Chem Inf Model 2023; 63:6778-6788. [PMID: 37897811 DOI: 10.1021/acs.jcim.3c00750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/30/2023]
Abstract
Accurate prediction of protein stability upon a point mutation has important applications in drug discovery and personalized medicine. It remains a challenging issue in computational biology. Existing computational prediction methods, which range from mechanistic to supervised learning approaches, have experienced limited progress over the last few decades. This stagnation is largely due to their heavy reliance on both the quantity and quality of the training data. This is evident in recent state-of-the-art methods that continue to yield substantial errors on two challenging blind test sets: frataxin and p53, with average root-mean-square errors exceeding 3 and 1.5 kcal/mol, respectively, which is still above the theoretical 1 kcal/mol prediction barrier. Rigorous approaches, on the other hand, offer greater potential for accuracy without relying on training data but are computationally demanding and require both wild-type and mutant structure information. Although they showed high accuracy for conserving mutations, their performance is still limited for charge-changing mutation cases. This might be due to the lack of an available mutant structure, often represented by a simplified capped peptide. The recent advances in protein structure prediction methods now make it possible to obtain structures comparable to experimental ones, including complete mutant structure information. In this work, we compare the performance of supervised learning-based methods and rigorous approaches for predicting protein stability on point mutations in difficult targets: frataxin and p53. The rigorous alchemical method significantly surpasses state-of-the-art techniques in terms of both the root-mean-squared error and Pearson correlation coefficient in these two challenging blind test sets. Additionally, we propose an improved alchemical method that employs the pmx double-system/single-box approach to accurately predict the folding free energy change upon both conserving and charge-changing mutations. The enhanced protocol can accurately predict both types of mutations, thereby outperforming existing state-of-the-art methods in overall performance.
Collapse
Affiliation(s)
- Jason Kurniawan
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan
| | - Takashi Ishida
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan
| |
Collapse
|
40
|
Swain PS, Panda S, Pati S, Dehury B. Computational saturation mutagenesis to explore the effect of pathogenic mutations on extra-cellular domains of TREM2 associated with Alzheimer's and Nasu-Hakola disease. J Mol Model 2023; 29:360. [PMID: 37924367 DOI: 10.1007/s00894-023-05770-7] [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: 07/25/2023] [Accepted: 10/25/2023] [Indexed: 11/06/2023]
Abstract
CONTEXT The specialised family of triggering receptors expressed on myeloid cells (TREMs) plays a pivotal role in causing neurodegenerative disorders and activating microglial anti-inflammatory responses. Nasu-Hakola disease (NHD), a rare autosomal recessive disorder, has been associated with mutations in TREM2, which is also responsible for raising the risk of Alzheimer's disease (AD). Herein, we have made an endeavour to differentiate the confirmed pathogenic variants in TREM2 extra-cellular domain (ECD) linked with NHD and AD using mutation-induced fold stability change (∆∆G), with the computation of 12distinct structure-based methods through saturation mutagenesis. Correlation analysis between relative solvent accessibility (RSA) and ∆∆G expresses the discrete distributive behaviour of mutants associated with TREM2 in AD (R2 = 0.061) and NHD (R2 = 0.601). Our findings put an emphasis on W50 and V126 as major players in maintaining V-like domain in TREM2. Interestingly, we discern that both of them interact with a common residue Y108, which is dissolved upon mutation. This Y108 could have structural or functional role for TREM2 which can be an ideal candidate for further study. Furthermore, the residual interaction network highlights the importance of R47 and R62 in maintaining the CDR loops that are crucial for ligand binding. Future studies using biophysical characterisation of ligand interactions in TREM2-ECD would be helpful for the development of novel therapeutics for AD and NHD. METHODS ConSurf algorithm and ENDscript were used to determine the position and conservation of each residue in the wild-type ECD of TREM2. The mutation-induced fold stability change (∆∆G) of confirmed pathogenic mutants associated with NHD and AD was estimated using 12 state-of-the-art structure-based protein stability tools. Furthermore, we also computed the effect of random mutation on these sites using computational saturation mutagenesis. Linear regression analysis was performed using mutants ∆∆G and RSA through GraphPad software. In addition, a comprehensive non-bonded residual interaction network (RIN) of wild type and its mutants of TREM2-ECD was enumerated using RING3.0.
Collapse
Affiliation(s)
- Preety Sthutika Swain
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, 751023, Odisha, India
| | - Sunita Panda
- Mycology Laboratory, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, 751023, Odisha, India
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, 751023, Odisha, India.
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Nalco Square, Chandrasekharpur, Bhubaneswar, 751023, Odisha, India.
| |
Collapse
|
41
|
Gong J, Jiang L, Chen Y, Zhang Y, Li X, Ma Z, Fu Z, He F, Sun P, Ren Z, Tian M. THPLM: a sequence-based deep learning framework for protein stability changes prediction upon point variations using pretrained protein language model. Bioinformatics 2023; 39:btad646. [PMID: 37874953 PMCID: PMC10627365 DOI: 10.1093/bioinformatics/btad646] [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: 04/29/2023] [Revised: 09/25/2023] [Accepted: 10/22/2023] [Indexed: 10/26/2023] Open
Abstract
MOTIVATION Quantitative determination of protein thermodynamic stability is a critical step in protein and drug design. Reliable prediction of protein stability changes caused by point variations contributes to developing-related fields. Over the past decades, dozens of structure-based and sequence-based methods have been proposed, showing good prediction performance. Despite the impressive progress, it is necessary to explore wild-type and variant protein representations to address the problem of how to represent the protein stability change in view of global sequence. With the development of structure prediction using learning-based methods, protein language models (PLMs) have shown accurate and high-quality predictions of protein structure. Because PLM captures the atomic-level structural information, it can help to understand how single-point variations cause functional changes. RESULTS Here, we proposed THPLM, a sequence-based deep learning model for stability change prediction using Meta's ESM-2. With ESM-2 and a simple convolutional neural network, THPLM achieved comparable or even better performance than most methods, including sequence-based and structure-based methods. Furthermore, the experimental results indicate that the PLM's ability to generate representations of sequence can effectively improve the ability of protein function prediction. AVAILABILITY AND IMPLEMENTATION The source code of THPLM and the testing data can be accessible through the following links: https://github.com/FPPGroup/THPLM.
Collapse
Affiliation(s)
- Jianting Gong
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Lili Jiang
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Yongbing Chen
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Yixiang Zhang
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Xue Li
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Zhiqiang Ma
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Department of Computer Science, College of Humanities and Sciences of Northeast Normal University, Changchun 130117, China
| | - Zhiguo Fu
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Fei He
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Pingping Sun
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
| | - Zilin Ren
- School of Information Science and Technology, Institution of Computational Biology, Northeast Normal University, Changchun 130117, China
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| | - Mingyao Tian
- Changchun Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Changchun 130122, China
| |
Collapse
|
42
|
Asif M, Chiou CC, Hussain MF, Hussain M, Sajid Z, Gulsher M, Raheem A, Khan A, Nasreen N, Kloczkowski A, Hassan M, Iqbal F, Chen CC. Homozygous Mutations in GDAP1 and MFN2 Genes Resulted in Autosomal Recessive Forms of Charcot-Marie-Tooth Disease in Consanguineous Pakistani Families. DNA Cell Biol 2023; 42:697-708. [PMID: 37797217 PMCID: PMC11262584 DOI: 10.1089/dna.2023.0169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 07/09/2023] [Accepted: 08/23/2023] [Indexed: 10/07/2023] Open
Abstract
Charcot-Marie-Tooth disease (CMT) is a heritable neurodegenerative disease of peripheral nervous system diseases in which more than 100 genes and their mutations are associated. Two consanguineous families Dera Ghazi Khan (PAK-CMT1-DG KHAN) and Layyah (PAK-CMT2-LAYYAH) with multiple CMT-affected subjects were enrolled from Punjab province in Pakistan. Basic epidemiological data were collected for the subjects. Nerve conduction study (NCS) and electromyography (EMG) were performed for the patients. Whole-exome sequencing (WES) followed by Sanger sequencing was applied to report the genetic basic of CMT. The NCS findings revealed that sensory and motor nerve conduction velocities for both families were <38 m/s. EMG presented denervation, neuropathic motor unit potential, and reduced interference pattern of peripheral nerves. WES identified that a novel nonsense mutation (c. 226 G>T) in GADP1 gene and a previously known missense mutation in MFN2 gene (c. 334 G>A) cause CMT4A (Charcot-Marie-Tooth disease type 4A) in the PAK-CMT1-DG KHAN family and CMT2A (Charcot-Marie-Tooth disease type 2A) in the PAK-CMT2-LAYYAH family, respectively. Mutations followed Mendelian pattern with autosomal recessive mode of inheritance. Multiple sequence alignment by Clustal Omega indicated that mutation-containing domain in both genes is highly conserved, and in situ analysis revealed that both mutations are likely to be pathogenic. We reported that a novel nonsense mutation and a previously known missense mutation in GAPD1 gene and MFN2 gene, respectively, cause CMT in consanguineous Pakistani families.
Collapse
Affiliation(s)
- Muhammad Asif
- Institute of Molecular Biology and Biotechnology. Bahauddin Zakariya University, Multan, Pakistan
- Institute of Zoology, Bahauddin Zakariya University, Multan, Pakistan
| | - Chien-Chun Chiou
- Department of Dermatology, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi, Taiwan
| | | | - Manzoor Hussain
- Orthopedic Unit 1, Nishter Medical University Multan, Pakistan
| | - Zureesha Sajid
- Institute of Molecular Biology and Biotechnology. Bahauddin Zakariya University, Multan, Pakistan
- Department of Biotechnology, Institute of Biochemistry, Biotechnology and Bioinformatics, Baghdad-ul-Jadeed Campus, The Islamia University of Bahawalpur, Bahawalpur, Pakistan
| | - Muhammad Gulsher
- Children Hospital and Institute of Child Health, Multan, Pakistan
| | - Afifa Raheem
- Institute of Zoology, Bahauddin Zakariya University, Multan, Pakistan
| | - Adil Khan
- Department of Botany and Zoology, Bacha Khan University, Charsadda, Pakistan
| | - Nasreen Nasreen
- Department of Zoology, Abdul Wali Khan University, Mardan, Pakistan
| | - Andrzej Kloczkowski
- The Steve and Cindy Rasmussen Institute for Genomic Medicine, Nationwide Children's Hospital, Columbus, Ohio, USA
- Department of Pediatrics, The Ohio State University, Columbus, Ohio, USA
| | - Mubashir Hassan
- The Steve and Cindy Rasmussen Institute for Genomic Medicine, Nationwide Children's Hospital, Columbus, Ohio, USA
| | - Furhan Iqbal
- Institute of Zoology, Bahauddin Zakariya University, Multan, Pakistan
| | - Chien-Chin Chen
- Department of Pathology, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Chiayi, Taiwan
- Department of Cosmetic Science, Chia Nan University of Pharmacy and Science, Tainan, Taiwan
- Ph.D. Program in Translational Medicine, Rong Hsing Research Center for Translational Medicine, National Chung Hsing University, Taichung, Taiwan
- Department of Biotechnology and Bioindustry Sciences, College of Bioscience and Biotechnology, National Cheng Kung University, Tainan, Taiwan
| |
Collapse
|
43
|
Umerenkov D, Nikolaev F, Shashkova TI, Strashnov PV, Sindeeva M, Shevtsov A, Ivanisenko NV, Kardymon OL. PROSTATA: a framework for protein stability assessment using transformers. Bioinformatics 2023; 39:btad671. [PMID: 37935419 PMCID: PMC10651431 DOI: 10.1093/bioinformatics/btad671] [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: 06/26/2023] [Revised: 10/25/2023] [Accepted: 11/02/2023] [Indexed: 11/09/2023] Open
Abstract
MOTIVATION Accurate prediction of change in protein stability due to point mutations is an attractive goal that remains unachieved. Despite the high interest in this area, little consideration has been given to the transformer architecture, which is dominant in many fields of machine learning. RESULTS In this work, we introduce PROSTATA, a predictive model built in a knowledge-transfer fashion on a new curated dataset. PROSTATA demonstrates advantage over existing solutions based on neural networks. We show that the large improvement margin is due to both the architecture of the model and the quality of the new training dataset. This work opens up opportunities to develop new lightweight and accurate models for protein stability assessment. AVAILABILITY AND IMPLEMENTATION PROSTATA is available at https://github.com/AIRI-Institute/PROSTATA and https://prostata.airi.net.
Collapse
Affiliation(s)
| | | | | | - Pavel V Strashnov
- Bioinformatics Group, AIRI, Moscow 121170, Russia
- Department of Computer Design and Technology, Bauman Moscow State Technical University, Moscow 105005, Russia
| | | | - Andrey Shevtsov
- Bioinformatics Group, AIRI, Moscow 121170, Russia
- Regulatory Transcriptomics and Epigenomics Group, Institute of Bioengineering, Research Center of Biotechnology RAS, Moscow 117036, Russia
| | - Nikita V Ivanisenko
- Bioinformatics Group, AIRI, Moscow 121170, Russia
- Laboratory of Computational Proteomics, Institute of Cytology and Genetics SB RAS, Novosibirsk 630090, Russia
| | | |
Collapse
|
44
|
Joshi BP, Bhandare VV, Vankawala M, Patel P, Patel R, Vyas B, Krishnamurty R. Friedelin, a novel inhibitor of CYP17A1 in prostate cancer from Cassia tora. J Biomol Struct Dyn 2023; 41:9695-9720. [PMID: 36373336 DOI: 10.1080/07391102.2022.2145497] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Accepted: 11/04/2022] [Indexed: 11/16/2022]
Abstract
In prostate cancer (PC), drugs targeting CYP17A1 have shown great success in regulating PC progression. However, successful drug molecules show adverse side effects and therapeutic resistance in PC. Therefore, we proposed to discover the potent phytochemical-based inhibitor against CYP17A1 using virtual screening. In this study, a phytochemicals library of ∼13800 molecules was selected to screen the best possible inhibitors against CYP17A1. A molecular modelling approach investigated detailed intermolecular interactions, their structural stability, and binding affinity. Further, in vitro and in vivo studies were performed to confirm the anticancer activity of identified potential inhibitor against CYP17A1. Friedelin from Cassia tora (CT) is identified as the best possible inhibitor from the screened library. MD simulation study reveals stable binding of Friedelin to conserved binding pocket of CYP17A1 with higher binding affinity than studied control, that is, Orteronel. Friedelin was tested on hormone-sensitive (22Rv1) and insensitive (DU145) cell lines and the IC50 value was found to be 72.025 and 81.766 µg/ml, respectively. CT extract showed a 25.28% IC50 value against 22Rv1, ∼92.6% increase in late Apoptosis/Necrosis, and three folds decrease in early apoptosis in treated cells compared to untreated cells. Further, animal studies show a marked decrease in prostate weight by 39.6% and prostate index by 36.5%, along with a reduction in serum PSA level by 71.7% and testosterone level by 92.4% compared to the testosterone group, which was further validated with histopathological studies. Thus, we propose Friedelin and CT extract as potential leads, which could be taken further for drug development in PC.[Figure: see text]Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | | | - Mahima Vankawala
- Leicester Institute of Structural and Chemical Biology, Department of Molecular and Cell Biology, University of Leicester, Leicester, UK
| | - Prittesh Patel
- C. G. Bhakta Institute of Biotechnology, Uka Tarsadia University, Tarsadi, Surat, Gujarat, India
| | - Rajesh Patel
- Bioinformatics and Supercomputer Lab., Department of Biosciences (UGC-SAP-DRS-II & DST-FIST-I), Veer Narmad South Gujarat University, Surat, Gujarat, India
| | - Bhavin Vyas
- Department of Pharmacology, Maliba Pharmacy College, Uka Tarsadia University, Tarsadi, Surat, Gujarat, India
| | - Ramar Krishnamurty
- C. G. Bhakta Institute of Biotechnology, Uka Tarsadia University, Tarsadi, Surat, Gujarat, India
| |
Collapse
|
45
|
Chen X, Leyendecker S, van den Bedem H. SARS-CoV-2 main protease mutation analysis via a kinematic method. Proteins 2023; 91:1496-1509. [PMID: 37408369 DOI: 10.1002/prot.26543] [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: 11/30/2022] [Revised: 05/23/2023] [Accepted: 06/08/2023] [Indexed: 07/07/2023]
Abstract
The Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is the virus responsible for the COVID-19 pandemic. COVID-19 continues to cause millions of deaths globally in part due to immune-evading mutations. SARS-CoV-2 main protease (Mpro) is an important enzyme for viral replication and potentially an effective drug target. Mutations affect the dynamics of enzymes and thereby their activity and ability to bind ligands. Here, we use kinematic flexibility analysis (KFA) to identify how mutations and ligand binding changes the conformational flexibility of Mpro. KFA decomposes macromolecules into regions of different flexibility near-instantly from a static structure, allowing conformational dynamics analysis at scale. Altogether, we analyzed 47 mutation sites across 69 Mpro-ligand complexes resulting in more than 3300 different structures which includes 69 mutated structures with all 47 sites mutated simultaneously and 3243 single residue mutated structures. We found that mutations generally increased the conformational flexibility of the protein. Understanding the impact of mutations on the flexibility of Mpro is essential for identifying potential drug targets in the treatment of SARS-CoV-2. Further studies in this area can offer valuable insights into the mechanisms of molecular recognition.
Collapse
Affiliation(s)
- Xiyu Chen
- Department of Mechanical Engineering, Institute of Applied Dynamics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Sigrid Leyendecker
- Department of Mechanical Engineering, Institute of Applied Dynamics, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Henry van den Bedem
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| |
Collapse
|
46
|
Azmi MB, Ahmed A, Ahmed TF, Imtiaz F, Asif U, Zaman U, Khan KA, Sherwani AK. Transcript-Level In Silico Analysis of Alzheimer's Disease-Related Gene Biomarkers and Their Evaluation with Bioactive Flavonoids to Explore Therapeutic Interactions. ACS OMEGA 2023; 8:40695-40712. [PMID: 37929088 PMCID: PMC10621018 DOI: 10.1021/acsomega.3c05769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/06/2023] [Accepted: 09/25/2023] [Indexed: 11/07/2023]
Abstract
Alzheimer's disease (AD) is a progressive brain disorder that can significantly affect the quality of life. We used a variety of in silico tools to investigate the transcript-level mutational impact of exonic missense rare variations (single nucleotide polymorphisms, SNPs) on protein function and to identify potential druggable protein cavities that correspond to potential therapeutic targets for the management of AD. According to the NIA-AA (National Institute on Aging-Alzheimer's Association) framework, we selected three AD biomarker genes (APP, NEFL, and MAPT). We systematically screened transcript-level exonic rare SNPs from these genes with a minor allele frequency of 1% in 1KGD (1000 Genomes Project Database) and gnomAD (Genome Aggregation Database). With downstream functional effect predictions, a single variation (rs182024939: K > N) of the MAPT gene with nine transcript SNPs was identified as the most pathogenic variation from the large dataset of mutations. The machine learning consensus classifier predictor categorized these transcript-level SNPs as the most deleterious variations, resulting in a large decrease in protein structural stability (ΔΔG kcal/mol). The bioactive flavonoid library was screened for drug-likeness and toxicity risk. Virtual screening of eligible flavonoids was performed using the MAPT protein. Identification of druggable protein-binding cavities showed VAL305, GLU655, and LYS657 as consensus-interacting residues present in the MAPT-docked top-ranked flavonoid compounds. The MM/PB(GB)SA analysis indicated hesperetin (-5.64 kcal/mol), eriodictyol (-5.63 kcal/mol), and sakuranetin (-5.60 kcal/mol) as the best docked flavonoids with the near-native binding pose. The findings of this study provide important insights into the potential of hesperetin as a promising flavonoid that can be utilized for further rational drug design and lead optimization to open new gateways in the field of AD therapeutics.
Collapse
Affiliation(s)
- Muhammad Bilal Azmi
- Department
of Biochemistry, Dow Medical College, Dow
University of Health Sciences, Karachi 74400, Pakistan
| | - Affan Ahmed
- Dow
Medical College, Dow University of Health
Sciences, Karachi 74400, Pakistan
| | - Tehniat Faraz Ahmed
- Department
of Biochemistry, Dow International Dental College, Dow University of Health Sciences, Karachi 75460, Pakistan
| | - Fauzia Imtiaz
- Department
of Biochemistry, Dow Medical College, Dow
University of Health Sciences, Karachi 74400, Pakistan
| | - Uzma Asif
- Department
of Biochemistry, Medicine Program, Batterjee
Medical College, Jeddah 21442, Saudi Arabia
| | - Uzma Zaman
- Department
of Biochemistry, Dow International Medical College, Dow University of Health Sciences, Karachi 74200, Pakistan
| | - Khalid Ali Khan
- Unit of Bee
Research and Honey Production, Research Center for Advanced Materials
Science (RCAMS), King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
- Applied
College, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
| | - Asif Khan Sherwani
- Research
and Development Unit, Jamjoom Pharmaceuticals
Co. Ltd, Jeddah 21442, Saudi Arabia
| |
Collapse
|
47
|
Rothfuss MT, Becht DC, Zeng B, McClelland LJ, Yates-Hansen C, Bowler BE. High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST. J Am Chem Soc 2023; 145:22979-22992. [PMID: 37815921 PMCID: PMC10626973 DOI: 10.1021/jacs.3c04966] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/12/2023]
Abstract
The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.
Collapse
Affiliation(s)
- Michael T. Rothfuss
- Department of Chemistry and Biochemistry, University of Montana, Missoula, MT 59812, United States
| | - Dustin C. Becht
- Department of Chemistry and Biochemistry, University of Montana, Missoula, MT 59812, United States
| | - Baisen Zeng
- Center for Biomolecular Structure and Dynamics, University of Montana, Missoula, MT 59812, United States
| | - Levi J. McClelland
- Center for Biomolecular Structure and Dynamics, University of Montana, Missoula, MT 59812, United States
- Division of Biological Sciences, University of Montana, Missoula, MT 59812, United States
| | - Cindee Yates-Hansen
- Center for Biomolecular Structure and Dynamics, University of Montana, Missoula, MT 59812, United States
| | - Bruce E. Bowler
- Department of Chemistry and Biochemistry, University of Montana, Missoula, MT 59812, United States
- Center for Biomolecular Structure and Dynamics, University of Montana, Missoula, MT 59812, United States
| |
Collapse
|
48
|
Ajith A, Subbiah U. In silico screening of non-synonymous SNPs in human TUFT1 gene. J Genet Eng Biotechnol 2023; 21:95. [PMID: 37801178 PMCID: PMC10558407 DOI: 10.1186/s43141-023-00551-4] [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: 11/21/2022] [Accepted: 09/20/2023] [Indexed: 10/07/2023]
Abstract
BACKGROUND Tuftelin 1 (TUFT1) gene is important in the development and mineralization of dental enamel. The study aimed to identify potential functionally deleterious non-synonymous SNPs (nsSNPs) in the TUFT1 gene by using different in silico tools. The deleterious missense SNPs were identified from SIFT, PolyPhen-2, PROVEAN, SNPs & GO, PANTHER, and SNAP2. The stabilization, conservation, and three-dimensional modeling of mutant proteins were analyzed by I-Mutant 3.0, Consurf, and Project HOPE, respectively. The protein-protein interaction using STRING, GeneMANIA for gene-gene interaction, and DynaMut for evaluating the impact of the mutation on protein stability, conformation, and flexibility. RESULTS Eight deleterious nsSNPs (E242A, R303W, K182N, K123N, R117W, H289Q, R203W, and Q107R) out of 304 were found to have high-risk damaging effects using six in silico tools. Among them, K182N and K123N alone had increased stability, whereas E242A, R303W, R117W, H289Q, Q107R, and R203W exhibited a decrease in protein stability, based on DDG values. Meanwhile, all the eight deleterious nsSNPs altered the size, charge, hydrophobicity, and spatial organization of the amino acids and predominantly had alpha helix domains. These deleterious variants were located in highly conserved regions except R203W. Protein-protein interaction predicted that TUFT1 interacted with ten proteins that are involved in enamel mineralization and odontogenesis. Gene-gene interaction network showed that TUFT1 is involved in physical interactions, gene co-localization, and pathway interactions. DynaMut ΔΔG values predicted that five nsSNPs were destabilizing the protein, ΔΔG ENCoM values showed a destabilizing effect for all mutants, and seven nsSNPs increased the molecular flexibility of TUFT1. CONCLUSION Our study predicted eight functional SNPs that had detrimental effects on the structure and function of the TUFT1 gene. This will aid in the development of candidate deleterious markers as a potential target for disease diagnosis and therapeutic interventions.
Collapse
Affiliation(s)
- Athira Ajith
- Human Genetics Research Centre, Sree Balaji Dental College and Hospital, Bharath Institute of Higher Education and Research, Chennai, 600 100, Tamil Nadu, India
| | - Usha Subbiah
- Human Genetics Research Centre, Sree Balaji Dental College and Hospital, Bharath Institute of Higher Education and Research, Chennai, 600 100, Tamil Nadu, India.
| |
Collapse
|
49
|
Sarmiento-Mañús R, Fontcuberta-Cervera S, González-Bayón R, Hannah MA, Álvarez-Martínez FJ, Barrajón-Catalán E, Micol V, Quesada V, Ponce MR, Micol JL. Analysis of the Arabidopsis venosa4-0 mutant supports the role of VENOSA4 in dNTP metabolism. PLANT SCIENCE : AN INTERNATIONAL JOURNAL OF EXPERIMENTAL PLANT BIOLOGY 2023; 335:111819. [PMID: 37562732 DOI: 10.1016/j.plantsci.2023.111819] [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: 04/04/2023] [Revised: 07/07/2023] [Accepted: 08/04/2023] [Indexed: 08/12/2023]
Abstract
Human Sterile alpha motif and histidine-aspartate domain containing protein 1 (SAMHD1) functions as a dNTPase to maintain dNTP pool balance. In eukaryotes, the limiting step in de novo dNTP biosynthesis is catalyzed by RIBONUCLEOTIDE REDUCTASE (RNR). In Arabidopsis, the RNR1 subunit of RNR is encoded by CRINKLED LEAVES 8 (CLS8), and RNR2 by three paralogous genes, including TSO MEANING 'UGLY' IN CHINESE 2 (TSO2). In plants, DIFFERENTIAL DEVELOPMENT OF VASCULAR ASSOCIATED CELLS 1 (DOV1) catalyzes the first step of the de novo biosynthesis of purines. Here, to explore the role of VENOSA4 (VEN4), the most likely Arabidopsis ortholog of human SAMHD1, we studied the ven4-0 point mutation, whose leaf phenotype was stronger than those of its insertional alleles. Structural predictions suggested that the E249L substitution in the mutated VEN4-0 protein rigidifies its 3D structure. The morphological phenotypes of the ven4, cls8, and dov1 single mutants were similar, and those of the ven4 tso2 and ven4 dov1 double mutants were synergistic. The ven4-0 mutant had reduced levels of four amino acids related to dNTP biosynthesis, including glutamine and glycine, which are precursors in the de novo purine biosynthesis. Our results reveal high functional conservation between VEN4 and SAMHD1 in dNTP metabolism.
Collapse
Affiliation(s)
- Raquel Sarmiento-Mañús
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | | | - Rebeca González-Bayón
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | - Matthew A Hannah
- Max-Planck-Institut für Molekulare Pflanzenphysiologie, Am Mühlenberg 1, 14476 Potsdam, Germany
| | - Francisco Javier Álvarez-Martínez
- Instituto de Investigación, Desarrollo e Innovación en Biotecnología Sanitaria de Elche, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | - Enrique Barrajón-Catalán
- Instituto de Investigación, Desarrollo e Innovación en Biotecnología Sanitaria de Elche, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | - Vicente Micol
- Instituto de Investigación, Desarrollo e Innovación en Biotecnología Sanitaria de Elche, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | - Víctor Quesada
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain
| | - María Rosa Ponce
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain.
| | - José Luis Micol
- Instituto de Bioingeniería, Universidad Miguel Hernández, Campus de Elche, 03202 Elche, Spain.
| |
Collapse
|
50
|
Cankara F, Doğan T. ASCARIS: Positional feature annotation and protein structure-based representation of single amino acid variations. Comput Struct Biotechnol J 2023; 21:4743-4758. [PMID: 37822561 PMCID: PMC10562615 DOI: 10.1016/j.csbj.2023.09.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2023] [Revised: 09/15/2023] [Accepted: 09/15/2023] [Indexed: 10/13/2023] Open
Abstract
Background Genomic variations may cause deleterious effects on protein functionality and perturb biological processes. Elucidating the effects of variations is critical for developing novel treatment strategies for diseases of genetic origin. Computational approaches have been aiding the work in this field by modeling and analyzing the mutational landscape. However, new approaches are required, especially for accurate representation and data-centric analysis of sequence variations. Method In this study, we propose ASCARIS (Annotation and StruCture-bAsed RepresentatIon of Single amino acid variations), a method for the featurization (i.e., quantitative representation) of single amino acid variations (SAVs), which could be used for a variety of purposes, such as predicting their functional effects or building multi-omics-based integrative models. ASCARIS utilizes the direct and spatial correspondence between the location of the SAV on the sequence/structure and 30 different types of positional feature annotations (e.g., active/lipidation/glycosylation sites; calcium/metal/DNA binding, inter/transmembrane regions, etc.), along with structural features and physicochemical properties. The main novelty of this method lies in constructing reusable numerical representations of SAVs via functional annotations. Results We statistically analyzed the relationship between these features and the consequences of variations and found that each carries information in this regard. To investigate potential applications of ASCARIS, we trained variant effect prediction models that utilize our SAV representations as input. We carried out an ablation study and a comparison against the state-of-the-art methods and observed that ASCARIS has a competing and complementary performance against widely-used predictors. ASCARIS can be used alone or in combination with other approaches to represent SAVs from a functional perspective. ASCARIS is available as a programmatic tool at https://github.com/HUBioDataLab/ASCARIS and as a web-service at https://huggingface.co/spaces/HUBioDataLab/ASCARIS.
Collapse
Affiliation(s)
- Fatma Cankara
- Biological Data Science Laboratory, Dept. of Computer Engineering, Hacettepe University, Ankara, Turkey
- Department of Health Informatics, Graduate School of Informatics, METU, Ankara, Turkey
- Department of Computational Sciences and Engineering, Koc University, Istanbul, Turkey
| | - Tunca Doğan
- Biological Data Science Laboratory, Dept. of Computer Engineering, Hacettepe University, Ankara, Turkey
- Institute of Informatics, Hacettepe University, Ankara, Turkey
- Department of Bioinformatics, Graduate School of Health Sciences, Hacettepe University, Ankara, Turkey
| |
Collapse
|