Roadmap to determine the point mutations involved in cardiomyopathy disorder: A Bayesian approach

In silico comprehensive analysis of coding and non-coding SNPs in human mTOR protein

The mammalian/mechanistic target of rapamycin (mTOR) protein is an important growth regulator and has been linked with multiple diseases including cancer and diabetes. Non-synonymous mutations of this gene have already been found in patients with renal clear cell carcinoma, melanoma, and acute lymphoid leukemia among many others. Such mutations can potentially affect a protein’s structure and hence its functions. In this study, therefore, the most deleterious SNPs of mTOR protein have been determined to identify potential biomarkers for various disease treatments. The aim is to generate a structured dataset of the mTOR gene’s SNPs that may prove to be an asset for the identification and treatment of multiple diseases associated with the target gene. Both sequence and structure-based approaches were adopted and a wide variety of bioinformatics tools were applied to analyze the SNPs of mTOR protein. In total 11 nsSNPs have been filtered out of 2178 nsSNPs along with two non-coding variations. All of the nsSNPs were found to destabilize the protein structure and disrupt its function. While R619C, A1513D, and T1977R mutations were shown to alter C alpha distances and bond angles of the mTOR protein, L509Q, R619C and N2043S were predicted to disrupt the mTOR protein’s interaction with NBS1 protein and FKBP1A/rapamycin complex. In addition, one of the non-coding SNPs was shown to alter miRNA binding sites. Characterizing nsSNPs and non-coding SNPs and their harmful effects on a protein’s structure and functions will enable researchers to understand the critical impact of mutations on the molecular mechanisms of various diseases. This will ultimately lead to the identification of potential targets for disease diagnosis and therapeutic interventions.

High-Risk Polymorphisms Associated with the Molecular Function of Human HMGCR Gene Infer the Inhibition of Cholesterol Biosynthesis

HMG-CoA reductase or HMGCR (3-hydroxy-3-methylglutaryl-CoA reductase) is a rate-limiting enzyme involved in cholesterol biosynthesis. HMGCR plays an important role in the possible occurrence of hypercholesterolemia leading to atherosclerosis and coronary heart disease. This enzyme is a major target for cholesterol-lowering drugs such as "statin" which blocks the synthesis of mevalonate, a precursor for cholesterol biosynthesis. This study is aimed at characterizing deleterious mutations and classifying functional single nucleotide polymorphisms (SNPs) of the HMGCR gene through analysis of functional and structural evaluation, domain association, solvent accessibility, and energy minimization studies. The functional and characterization tools such as SIFT, PolyPhen, SNPs and GO, Panther, I-Mutant, and Pfam along with programming were employed to explore all the available SNPs in the HMGCR gene in the database. Among 6815 SNP entries from different databases, approximately 388 SNPs were found to be missense. Analysis showed that seven missense SNPs are more likely to have deleterious effects. A tertiary model of the mutant protein was constructed to determine the functional and structural effects of the HMGCR mutation. In addition, the location of the mutations suggests that they may have deleterious effects because most of the mutations are residing in the functional domain of the protein. The findings from the analysis predicted that rs147043821 and rs193026499 missense SNPs could cause significant structural and functional instability in the mutated proteins of the HMGCR gene. The findings of the current study will likely be useful in future efforts to uncover the mechanism and cause of hypercholesterolemia. In addition, the identified SNPs of HMGCR gene could set up a strong foundation for further therapeutic discovery.

Isolation and structure elucidation of a steroidal moiety from Withania somnifera and in silico evaluation of antimalarial efficacy against artemisinin resistance Plasmodium falciparum kelch 13 protein

According to the 2021 Malaria report, 241 million clinical episodes with 627000 deaths penalty was estimated across the worldwide. However, mutation in the propeller domain of Plasmodium falciparum kelch 13 protein resulted in longer parasite clearance time following an artemisinin-based treatment and had a greater survival rate of ring-stage parasites even after a brief exposure to a high dose of artesunate. Clinical manifestations become more complex and worse with the emerging trend of drug resistance against artemisinin derivatives and the poor effectiveness of malaria vaccination drive. Steroidal lactone (withanolide) moiety (C-28) isolated from methanolic leaf extract Withania somnifera show a greater affinity towards Pfkelch 13 protein in comparison to the artemisinin derivatives (artesunate, artemether). The isolated compound was characterized to be withaferin A with a percentage yield of 29.01% w/w in chloroform fraction, 1.75% w/w in methanolic extract, and 0.29% w/w in raw leaf powder. Structure-based analysis shows that withaferin A (docking score -8.253, -9.802) has a higher affinity for two distinct binding pockets I and II of the Plasmodium falciparum kelch 13 protein than artesunate (docking score -4.470, -3.656). Further, Gibbs binding free energy signifies thermodynamic stability of the docked complex of withaferin A (-43.25, -43.76 Kcal/mol) in comparison to artesunate docked complex (-8.49, -5.75 Kcal/mol). The pharmacokinetic profile of withaferin A shows more drug-likeness characteristics without violating Jorgensen's rule of three, and Lipinski's rule of five. Hence above experimental findings suggest withaferin A could be a suitable therapeutic adjunct for preclinical evaluation of antimalarial potentiality in artemisinin-resistant malaria.

In-silico phenotype prediction by normal mode variant analysis in TUBB4A-related disease

TUBB4A-associated disorder is a rare condition affecting the central nervous system. It displays a wide phenotypic spectrum, ranging from isolated late-onset torsion dystonia to a severe early-onset disease with developmental delay, neurological deficits, and atrophy of the basal ganglia and cerebellum, therefore complicating variant interpretation and phenotype prediction in patients carrying TUBB4A variants. We applied entropy-based normal mode analysis (NMA) to investigate genotype–phenotype correlations in TUBB4A-releated disease and to develop an in-silico approach to assist in variant interpretation and phenotype prediction in this disorder. Variants included in our analysis were those reported prior to the conclusion of data collection for this study in October 2019. All TUBB4A pathogenic missense variants reported in ClinVar and Pubmed, for which associated clinical information was available, and all benign/likely benign TUBB4A missense variants reported in ClinVar, were included in the analysis. Pathogenic variants were divided into five phenotypic subgroups. In-silico point mutagenesis in the wild-type modeled protein structure was performed for each variant. Wild-type and mutated structures were analyzed by coarse-grained NMA to quantify protein stability as entropy difference value (ΔG) for each variant. Pairwise ΔG differences between all variant pairs in each structural cluster were calculated and clustered into dendrograms. Our search yielded 41 TUBB4A pathogenic variants in 126 patients, divided into 11 partially overlapping structural clusters across the TUBB4A protein. ΔG-based cluster analysis of the NMA results revealed a continuum of genotype–phenotype correlation across each structural cluster, as well as in transition areas of partially overlapping structural clusters. Benign/likely benign variants were integrated into the genotype–phenotype continuum as expected and were clearly separated from pathogenic variants. We conclude that our results support the incorporation of the NMA-based approach used in this study in the interpretation of variant pathogenicity and phenotype prediction in TUBB4A-related disease. Moreover, our results suggest that NMA may be of value in variant interpretation in additional monogenic conditions.

Computational Analysis of Gly482Ser Single-Nucleotide Polymorphism in PPARGC1A Gene Associated with CAD, NAFLD, T2DM, Obesity, Hypertension, and Metabolic Diseases

Peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PPARGC1A) regulates the expression of energy metabolism's genes and mitochondrial biogenesis. The essential roles of PPARGC1A encouraged the researchers to assess the relation between metabolism-related diseases and its variants. To study Gly482Ser (+1564G/A) single-nucleotide polymorphism (SNP) after PPARGC1A modeling, we substitute Gly482 for Ser482. Stability prediction tools showed that this substitution decreases the stability of PPARGC1A or has a destabilizing effect on this protein. We then utilized molecular dynamics simulation of both the Gly482Ser variant and wild type of the PPARGC1A protein to analyze the structural changes and to reveal the conformational flexibility of the PPARGC1A protein. We observed loss flexibility in the RMSD plot of the Gly482Ser variant, which was further supported by a decrease in the SASA value in the Gly482Ser variant structure of PPARGC1A and an increase of H-bond with the increase of β-sheet and coil and decrease of turn in the DSSP plot of the Gly482Ser variant. Such alterations may significantly impact the structural conformation of the PPARGC1A protein, and it might also affect its function. It showed that the Gly482Ser variant affects the PPARGC1A structure and makes the backbone less flexible to move. In general, molecular dynamics simulation (MDS) showed more flexibility in the native PPARGC1A structure. Essential dynamics (ED) also revealed that the range of eigenvectors in the conformational space has lower extension of motion in the Gly482Ser variant compared with WT. The Gly482Ser variant also disrupts PPARGC1A interaction. Due to this single-nucleotide polymorphism in PPARGC1A, it became more rigid and might disarray the structural conformation and catalytic function of the protein and might also induce type 2 diabetes mellitus (T2DM), coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). The results obtained from this study will assist wet lab research in expanding potent treatment on T2DM.

The promising novel biomarkers and candidate small molecule drugs in lower-grade glioma: Evidence from bioinformatics analysis of high-throughput data

Overall survival of patients with low-grade glioma (LGG) has shown no significant improvement over the past 30 years, with survival averaging approximately 7 years. This study aimed to identify novel promising biomarkers of LGG and reveal its potential molecular mechanisms by integrated bioinformatics analysis. The microarray datasets of GSE68848 and GSE4290 were selected from GEO database for integrated analysis. In total, 293 overlapping differentially expressed genes (DEGs) were detected using the limma package. One hundred and eighty-eight nodes with 603 interactions were obtained from the establishment of protein-protein interaction (PPI) network. Functional and signaling pathway enriched were significantly correlated with the synapse and calcium signaling pathway, respectively. Module analysis revealed eight hub genes with high connectivity, which included CHRM1, DLG2, GABRD, GRIN1, HTR2A, KCNJ3, KCNJ9, and NUSAP1, and they were markedly correlated with patients' prognosis. The mining of the Gene Expression Profiling Interactive Analysis database and qPCR further confirmed the abnormal expression of these key genes with their prognostic value in LGG. We eventually predicted the 20 most vital small molecule drugs, which potentially reverse the carcinogenic state of LGG, as per the CMap (connectivity map) database and these DEGs, and MS-275 (enrichment score = -0.939) was considered as the most promising small molecule to treat LGG. In conclusion, our study provided eight reliable novel molecular biomarkers for diagnosis, prognosis prediction, and treatment targets for LGG. These conclusions will contribute to a better comprehension of molecular mechanisms fundamental to LGG occurrence and progression, and providing new insights for future development of genomic individualized treatment in LGG.

Prediction of a highly deleterious mutation E17K in AKT-1 gene: An in silico approach

The AKT1 (v-akt murine thymoma viral oncogene homologue 1) kinase is a member of most frequently activated proliferation and survival signaling pathway in cancer. Recently, hyperactivation of AKT1, due to functional point mutation in the pleckstrin homology (PH) domain of AKT1 gene, has been found to be associated with human colorectal, breast and ovarian cancer. Thus, considering its crucial role in cellular signaling pathway, a functional analysis of missense mutations of AKT1 gene was undertaken in this study. Twenty nine nsSNPs (non-synonymous single nucleotide polymorphism) within coding region of AKT1 gene were selected for our investigation and six SNPs were found to be deleterious by combinatorial predictions of various computational tools. RMSD values were calculated for the mutant models which predicted four substitutions (E17K, E319G, D32E and A255T) to be highly deleterious. The insight of the structural attribute was gained through analysis of, secondary structures, solvent accessibility and intermolecular hydrogen bond analysis which confirmed one missense mutation (E17K) to be highly deleterious nsSNPs. In conclusion, the investigated gene AKT1 has twenty nine SNPs in the coding region and through progressive analysis using different bioinformatics tools one highly deleterious SNP with rs121434592 was profiled. Thus, results of this study can pave a new platform to sort nsSNPs for several important regulatory genes that can be undertaken for the confirmation of their phenotype and their correlation with diseased status in case control studies.

•

We have added a small portion of text in introduction part as per reviewers comment.

•

We have added a schematic representation of methodology used (Fig. 1).

•

We have added text in the discussion portion as per the comment of reviewer.

•

We have also corrected the conclusion as per reviewer's comments.

In silico approach to explore the disruption in the molecular mechanism of human hyaluronidase 1 by mutant E268K that directs Natowicz syndrome

Natowicz syndrome (mucopolysaccharidoses type 9) is a lysosomal storage disorder caused by deficient or defective human hyaluronidase 1. The disorder is not well studied at the molecular level. Therefore, a new in silico approach was proposed to study the molecular basis on which one clinically observed mutation, Glu268Lys, results in a defective enzyme. The native and mutant structures were subjected to comparative analyses using a conformational sampling approach for geometrical variables viz, RMSF, RMSD, and Ramachandran plot. In addition, the strength of a Cys207-Cys221 disulfide bond and electrostatic interaction between Arg265 and Asp206 were studied, as they are known to be involved in the catalytic activity of the enzyme. Native and mutant E268K showed statistically significant variations with p < 0.05 in RMSD, Ramachandran plot, strengths of disulfide bond, and electrostatic interactions. Further, single model analysis showed variations between native and mutant structures in terms of intra-protein interactions, hydrogen bond dilution, secondary structure, and dihedral angles. Docking analysis predicted the mutant to have a less favorable substrate binding energy compared to the native protein. Additionally, steered MD analysis indicated that the substrate should have more affinity to the native than mutant enzymes. The observed changes theoretically explain the less favorable binding energy of substrate towards mutant E268K, thereby providing a structural basis for its reduced catalytic activity. Hence, our study provides a basis for understanding the disruption in the molecular mechanism of human hyaluronidase 1 by mutation E268K, which may prove useful for the development of synthetic chaperones as a treatment option for Natowicz syndrome.

The structural effects of mutations can aid in differential phenotype prediction of beta-myosin heavy chain (Myosin-7) missense variants

MOTIVATION

High-throughput sequencing platforms are increasingly used to screen patients with genetic disease for pathogenic mutations, but prediction of the effects of mutations remains challenging. Previously we developed SAAPdap (Single Amino Acid Polymorphism Data Analysis Pipeline) and SAAPpred (Single Amino Acid Polymorphism Predictor) that use a combination of rule-based structural measures to predict whether a missense genetic variant is pathogenic. Here we investigate whether the same methodology can be used to develop a differential phenotype predictor, which, once a mutation has been predicted as pathogenic, is able to distinguish between phenotypes-in this case the two major clinical phenotypes (hypertrophic cardiomyopathy, HCM and dilated cardiomyopathy, DCM) associated with mutations in the beta-myosin heavy chain (MYH7) gene product (Myosin-7).

RESULTS

A random forest predictor trained on rule-based structural analyses together with structural clustering data gave a Matthews' correlation coefficient (MCC) of 0.53 (accuracy, 75%). A post hoc removal of machine learning models that performed particularly badly, increased the performance (MCC = 0.61, Acc = 79%). This proof of concept suggests that methods used for pathogenicity prediction can be extended for use in differential phenotype prediction.

AVAILABILITY AND IMPLEMENTATION

Analyses were implemented in Perl and C and used the Java-based Weka machine learning environment. Please contact the authors for availability.

CONTACTS

andrew@bioinf.org.uk or andrew.martin@ucl.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Novel Mutation of SMAD3 Identified in a Chinese Family with Aneurysms-Osteoarthritis Syndrome

Aneurysms-osteoarthritis syndrome (AOS) is a recently delineated autosomal dominant disorder characterized by aneurysms, dissections, and tortuosity throughout the arterial tree in association with early onset osteoarthritis, mild craniofacial features, and skeletal and cutaneous anomalies. Previous studies have demonstrated that mutations in SMAD3, a key regulator of TGF-β signal transduction, contribute to AOS. Here, we investigated a family of three generations affected by AOS. A novel SMAD3 mutation, c.266G>A (p.C89Y), was identified and cosegregated with the affected individuals in this family. Our finding expands the mutation spectrum of SMAD3 gene and further strengthens the connection between the presence of aneurysms-osteoarthritis phenotype and SMAD3 mutations, which facilitates the understanding of the genotype-phenotype correlation of AOS.

In-silico screening of cancer associated mutation on PLK1 protein and its structural consequences

The Polo-like kinases (Plks) are a conserved subfamily of serine-threonine protein kinases that have significant roles in cell proliferation. The serine/threonine protein kinases or polo-like kinase 1 (PLK1) exist in centrosome during interphase and is an important regulatory enzyme in cell cycle progression during M phase. Mutations in mammalian PLK1 were found to be over expressed in various human cancers and it is disrupting the binding ability of polo box domain with target peptide. In this analysis we implemented a computational approach to filter the most deleterious and cancer associated mutation on PLK1 protein. We found W414F as the most deleterious and cancer associated by Polyphen 2.0, SIFT, I-mutant 3.0, PANTHER, PhD-SNP, SNP&GO, Mutpred and Dr Cancer tools. Molecular docking and molecular dynamics simulation (MDS) approach was used to investigate the structural and functional behavior of PLK1 protein upon mutation. MDS and docking results showed stability loss in mutant PLK1 protein. Due to mutation, PLK1 protein became more flexible and alters the dynamic property of protein which might affect the interaction with target peptide and leads to cell proliferation. Our study provided a well designed computational methodology to examine the cancer associated nsSNPs and their molecular mechanism. It further helps scientists to develop a drug therapy against PLK1 cancer-associated diseases.