Structural analysis of missense mutations in galactokinase 1 (GALK1) leading to galactosemia type-2

Computational identification and analysis of deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene: a structural and functional impact

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.

Whole-genome resequencing of Dorper and Hu sheep to reveal selection signatures associated with important traits

Dorper and Hu sheep exhibit different characteristics in terms of reproduction, growth, and meat quality. Comparison of the genomes of two breeds help to reveal important genomic information. In this study, whole genome resequencing of 30 individuals (Dorper, DB and Hu sheep, HY) identified 15,108,125 single nucleotide polymorphisms (SNPs). Population differentiation (Fst) and cross population extended haplotype homozygosity (XP-EHH) were performed for selective signal analysis. In total, 106 and 515 overlapped genes were present in both the Fst results and XP-EHH results in HY vs DB and in DB vs HY, respectively. In HY vs DB, 106 genes were enriched in 12 GO terms and 83 KEGG pathways, such as ATP binding (GO:0005524) and PI3K-Akt signaling pathway (oas04151). In DB vs HY, 515 genes were enriched in 109 GO terms and 215 KEGG pathways, such as skeletal muscle cell differentiation (GO:0035914) and MAPK signaling pathway (oas04010). According to the annotation results, we identified a series of candidate genes associated with reproduction (UNC5C, BMPR1B, and GLIS1), meat quality (MECOM, MEF2C, and MYF6), and immunity (GMDS, GALK1, and ITGB4). Our investigation has uncovered genomic information for important traits in sheep and provided a basis for subsequent studies of related traits.

Gene network interaction analysis to elucidate the antimicrobial resistance mechanisms in the Clostridiumdifficile

Antimicrobial resistance has caused chaos worldwide due to the depiction of multidrug-resistant (MDR) infective microorganisms. A thorough examination of antimicrobial resistance (AMR) genes and associated resistant mechanisms is vital to solving this problem. Clostridium difficile (C. difficile) is an opportunistic nosocomial bacterial strain that has acquired exogenous AMR genes that confer resistance to antimicrobials such as erythromycin, azithromycin, clarithromycin, rifampicin, moxifloxacin, fluoroquinolones, vancomycin, and others. A network of interactions, including 20 AMR genes, was created and analyzed. In functional enrichment analysis, Cellular components (CC), Molecular Functions (MF), and Biological Processes (BP) were discovered to have substantial involvement. Mutations in the rpl genes, which encode ribosomal proteins, confer resistance in Gram-positive bacteria. Full erythromycin and azithromycin cross-resistance can be conferred if more than one of the abovementioned genes is present. In the enriched BP, rps genes related to transcriptional regulation and biosynthesis were found. The genes belong to the rpoB gene family, which has previously been related to rifampicin resistance. The genes rpoB, gyrA, gyrB, rpoS, rpl genes, rps genes, and Van genes are thought to be the hub genes implicated in resistance in C. difficile. As a result, new medications could be developed using these genes. Overall, our observations provide a thorough understanding of C. difficile AMR mechanisms.

High risk genetic variants of human insulin receptor substrate 1(IRS1) infer structural instability and functional interference

Insulin receptor substrate 1(IRS1) is a signaling adapter protein encoded by the IRS1 gene. This protein delivers signals from insulin and insulin-like growth factor-1(IGF-1) receptors to the phosphatidylinositol 3-kinases (P13K)/protein kinase B (Akt) and Extracellular signal-regulated kinases (Erk) - Mitogen-activated protein (MAP) kinase pathways, which regulate particular cellular processes. Mutations in this gene have been linked to type 2 diabetes mellitus, a heightened risk of insulin resistance, and an increased likelihood of developing a number of different malignancies. The structure and function of IRS1 could be severely compromised as a result of single nucleotide polymorphism (SNP) type genetic variants. In this study, we focused on identification of the most harmful non-synonymous SNPs (nsSNPs) of the IRS1 gene as well as prediction of their structural and functional consequences. Six different algorithms made the initial prediction that 59 of the 1142 IRS1 nsSNPs would have a negative impact on the protein structure. In-depth evaluations detected 26 nsSNPs located inside the functional domains of IRS1. Following that, 16 nsSNPs were identified as more harmful based on conservation profile, hydrophobic interaction, surface accessibility, homology modelling, and inter-atomic interactions. Following an in-depth analysis of protein stability, M249T (rs373826433), I223T (rs1939785175) and V204G (rs1574667052) were identified as three most deleterious SNPs and were subjected to molecular dynamics simulation for further insights. These findings will help us understand the implications for disease susceptibility, cancer progression, and the efficacy of therapeutic development against IRS1 gene mutants.Communicated by Ramaswamy H. Sarma.

Computational refinement identifies functional destructive single nucleotide polymorphisms associated with human retinoid X receptor gene

Alterations in the nuclear retinoid X receptor (RXRs) signalling have been implicated in neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, stroke, multiple sclerosis and glaucoma. Single nucleotide polymorphisms (SNPs) are the main cause underlying single nucleic acid variations which in turn determine heterogeneity within various populations. These genetic polymorphisms have been suggested to associate with various degenerative disorders in population-wide analysis. This bioinformatics study was designed to investigate, search, retrieve and identify deleterious SNPs which may affect the structure and function of various RXR isoforms through a computational and molecular modelling approach. Amongst the 1,813 retrieved SNPs several were found to be deleterious with rs140464195_G139R, rs368400425_R358W and rs368586400_L383F RXRα mutant variants being the most detrimental ones causing changes in the interatomic interactions and decreasing the flexibility of the mutant proteins. Molecular genetics analysis identified seven missense mutations in RXRα/β/γ isoforms. Two novel mutations SNP IDs (rs1588299621 and rs1057519958) were identified in RXRα isoform. We used several in silico prediction tools such as SIFT, PolyPhen, I-Mutant, Protein Variation Effect Analyzer (PROVEAN), PANTHER, SNP&Go, PhD-SNP and SNPeffect to predict pathogenicity and protein stability associated with RXR mutations. The structural assessment by DynaMut tool revealed that hydrogen bonds were affected along with hydrophobic and carbonyl interactions resulting in reduced flexibility at the mutated residue positions but ultimately stabilizing the molecule as a whole. Summarizing, analysis of the missense mutations in RXR isoforms showed a mix of conclusive and inconclusive genotype-phenotype correlations suggesting the use of sophisticated computational analysis tools for studying RXR variants.Communicated by Ramaswamy H. Sarma.

Establishment of a risk score model for bladder urothelial carcinoma based on energy metabolism-related genes and their relationships with immune infiltration

Bladder urothelial carcinoma (BLCA) is a common malignant tumor of the human urinary system, and a large proportion of BLCA patients have a poor prognosis. Therefore, there is an urgent need to find more efficient and sensitive biomarkers for the prognosis of BLCA patients in clinical practice. RNA sequencing (RNA-seq) data and clinical information were obtained from The Cancer Genome Atlas, and 584 energy metabolism-related genes (EMRGs) were obtained from the Reactome pathway database. Cox regression analysis and least absolute shrinkage and selection operator analysis were applied to assess prognostic genes and build a risk score model. The estimate and cibersort algorithms were used to explore the immune microenvironment, immune infiltration, and checkpoints in BLCA patients. Furthermore, we used the Human Protein Atlas database and our single-cell RNA-seq datasets of BLCA patients to verify the expression of 13 EMRGs at the protein and single-cell levels. We constructed a risk score model; the area under the curve of the model at 5 years was 0.792. The risk score was significantly correlated with the immune markers M0 macrophages, M2 macrophages, CD8 T cells, follicular helper T cells, regulatory T cells, and dendritic activating cells. Furthermore, eight immune checkpoint genes were significantly upregulated in the high-risk group. The risk score model can accurately predict the prognosis of BLCA patients and has clinical application value. In addition, according to the differences in immune infiltration and checkpoints, BLCA patients with the most significant benefit can be selected for immune checkpoint inhibitor therapy.

Bioinformatics investigation on blood-based gene expressions of Alzheimer's disease revealed ORAI2 gene biomarker susceptibility: An explainable artificial intelligence-based approach

The progressive, chronic nature of Alzheimer's disease (AD), a form of dementia, defaces the adulthood of elderly individuals. The pathogenesis of the condition is primarily unascertained, turning the treatment efficacy more arduous. Therefore, understanding the genetic etiology of AD is essential to identifying targeted therapeutics. This study aimed to use machine-learning techniques of expressed genes in patients with AD to identify potential biomarkers that can be used for future therapy. The dataset is accessed from the Gene Expression Omnibus (GEO) database (Accession Number: GSE36980). The subgroups (AD blood samples from frontal, hippocampal, and temporal regions) are individually investigated against non-AD models. Prioritized gene cluster analyses are conducted with the STRING database. The candidate gene biomarkers were trained with various supervised machine-learning (ML) classification algorithms. The interpretation of the model prediction is perpetrated with explainable artificial intelligence (AI) techniques. This experiment revealed 34, 60, and 28 genes as target biomarkers of AD mapped from the frontal, hippocampal, and temporal regions. It is identified ORAI2 as a shared biomarker in all three areas strongly associated with AD's progression. The pathway analysis showed that STIM1 and TRPC3 are strongly associated with ORAI2. We found three hub genes, TPI1, STIM1, and TRPC3, in the network of the ORAI2 gene that might be involved in the molecular pathogenesis of AD. Naive Bayes classified the samples of different groups by fivefold cross-validation with 100% accuracy. AI and ML are promising tools in identifying disease-associated genes that will advance the field of targeted therapeutics against genetic diseases.

Current Understanding on the Genetic Basis of Key Metabolic Disorders: A Review

Simple Summary

Metabolic disorders (MD) are a challenge to healthcare systems; the emergence of the modern socio-economic system has led to a profound change in lifestyles in terms of dietary habits, exercise regimens, and behavior, all of which complement the genetic factors associated with MD. Diabetes Mellitus and Familial hypercholesterolemia are two of the 14 most widely researched MD, as they pose the greatest challenge to the public healthcare system and have an impact on productivity and the economy. Research findings have led to the development of new therapeutic molecules for the mitigation of MD as well as the invention of experimental strategies, which target the genes themselves via gene editing and RNA interference. Although these approaches may herald the emergence of a new toolbox to treat MD, the current therapeutic approaches still heavily depend on substrate reduction, dietary restrictions based on genetic factors, exercise, and the maintenance of good mental health. The development of orphan drugs for the less common MD such as Krabbe, Farber, Fabry, and Gaucher diseases, remains in its infancy, owing to the lack of investment in research and development, and this has driven the development of personalized therapeutics based on gene silencing and related technologies.

Abstract

Advances in data acquisition via high resolution genomic, transcriptomic, proteomic and metabolomic platforms have driven the discovery of the underlying factors associated with metabolic disorders (MD) and led to interventions that target the underlying genetic causes as well as lifestyle changes and dietary regulation. The review focuses on fourteen of the most widely studied inherited MD, which are familial hypercholesterolemia, Gaucher disease, Hunter syndrome, Krabbe disease, Maple syrup urine disease, Metachromatic leukodystrophy, Mitochondrial encephalopathy lactic acidosis stroke-like episodes (MELAS), Niemann-Pick disease, Phenylketonuria (PKU), Porphyria, Tay-Sachs disease, Wilson’s disease, Familial hypertriglyceridemia (F-HTG) and Galactosemia based on genome wide association studies, epigenetic factors, transcript regulation, post-translational genetic modifications and biomarker discovery through metabolomic studies. We will delve into the current approaches being undertaken to analyze metadata using bioinformatic approaches and the emerging interventions using genome editing platforms as applied to animal models.

Galactosemia: Biochemistry, Molecular Genetics, Newborn Screening, and Treatment

Galactosemia is an inborn disorder of carbohydrate metabolism characterized by the inability to metabolize galactose, a sugar contained in milk (the main source of nourishment for infants), and convert it into glucose, the sugar used by the body as the primary source of energy. Galactosemia is an autosomal recessive genetic disease that can be diagnosed at birth, even in the absence of symptoms, with newborn screening by assessing the level of galactose and the GALT enzyme activity, as GALT defect constitutes the most frequent cause of galactosemia. Currently, galactosemia cannot be cured, but only treated by means of a diet with a reduced content of galactose and lactose. Although the diet is able to reverse the neonatal clinical picture, it does not prevent the development of long-term complications. This review provides an overview of galactose metabolism, molecular genetics, newborn screening and therapy of galactosemia. Novel treatments for galactosemia currently being investigated in (pre)clinical studies and potentially able to prevent long-term complications are also presented.

Reduced off-target effect of NG-BE4max by using NG-HiFi system

Recently, a rationally engineered SpCas9 variant (SpCas9-NG) that can recognize a minimal NG protospacer adjacent motif (PAM) was reported to expand the targeting scope in genome editing. However, increased genome-wide off-target mutations with this variant compared with SpCas9 were reported in previous studies. In addition, lower base editing frequencies and higher unintended off-target mutations were also found in Hoxc13-ablated rabbits generated by NG-BE4max in our study. Here, a high-fidelity base editor, NG-HiFi, in comparison to NG-BE4max, showed retention of on-target activity while exhibiting significantly decreased off-target activity in Hoxc13-ablated rabbits. Collectively, the improved specificity and reduced off-target effect of SpCas9-NG assisted in cytidine base editing with the NG-HiFi system, providing a promising tool to precisely model human diseases in rabbits.

[13C]-galactose breath test in a patient with galactokinase deficiency and spastic diparesis

Galactokinase deficiency is an inborn error of carbohydrate metabolism due to a block in the formation of galactose-1-phosphate from galactose. Although the association of galactokinase deficiency with formation of cataracts is well established, the extent of the clinical phenotype is still under investigation. We describe a 6-year-old female who was diagnosed with galactokinase deficiency due to cataract formation when she was 10 months of age and initially started on galactose-restricted diet at that time for 5 months. She developed gait abnormality at 4 years of age. Breath tests via measurement of 13C isotope in exhaled carbon dioxide following 13C-labeled galactose administration at carbon-1 and carbon-2 positions revealed oxidation rates within the normal range. The results in this patient strikingly contrast with the results of another patient with GALK1 deficiency that underwent breath testing with [1-14C]-galactose and [2-14C]-galactose. Extension of in vivo breath tests to other galactokinase patients is needed to better understand the pathophysiology of this disease.

Molecular dynamics, residue network analysis, and cross-correlation matrix to characterize the deleterious missense mutations in GALE causing galactosemia III

Epimerase-deficiency galactosemia (EDG) is caused by mutations in the UDP-galactose 4'-epimerase enzyme, encoded by gene GALE. Catalyzing the last reaction in the Leloir pathway, UDP-galactose-4-epimerase catalyzes the interconversion of UDP-galactose and UDP-glucose. This study aimed to use in-depth computational strategies to prioritize the pathogenic missense mutations in GALE protein and investigate the systemic behavior, conformational spaces, atomic motions, and cross-correlation matrix of the GALE protein. We searched four databases (dbSNP, ClinVar, UniProt, and HGMD) and major biological literature databases (PubMed, Science Direct, and Google Scholar), for missense mutations that are associated with EDG patients, our search yielded 190 missense mutations. We applied a systematic computational prediction pipeline, including pathogenicity, stability, biochemical, conservational, protein residue contacts, and structural analysis, to predict the pathogenicity of these mutations. We found three mutations (p.K161N, p.R239W, and p.G302D) with a severe phenotype in patients with EDG that correlated with our computational prediction analysis; thus, they were selected for further structural and simulation analyses to compute the flexibility and stability of the mutant GALE proteins. The three mutants were subjected to molecular dynamics simulation (MDS) with native protein for 200 ns using GROMACS. The MDS demonstrated that these mutations affected the beta-sheets and helical region that are responsible for the catalytic activity; subsequently, affects the stability and flexibility of the mutant proteins along with a decrease and more deviations in compactness when compared to that of a native. Also, three mutations created major variations in the combined atomic motions of the catalytic and C-terminal regions. The network analysis of the residues in the native and three mutant protein structures showed disturbed residue contacts occurred owing to the missense mutations. Our findings help to understand the structural behavior of a protein owing to mutation and are intended to serve as a platform for prioritizing mutations, which could be potential targets for drug discovery and development of targeted therapeutics.

Galactosemia: Towards Pharmacological Chaperones

Galactosemia is a rare inherited metabolic disease resulting from mutations in the four genes which encode enzymes involved in the metabolism of galactose. The current therapy, the removal of galactose from the diet, is inadequate. Consequently, many patients suffer lifelong physical and cognitive disability. The phenotype varies from almost asymptomatic to life-threatening disability. The fundamental biochemical cause of the disease is a decrease in enzymatic activity due to failure of the affected protein to fold and/or function correctly. Many novel therapies have been proposed for the treatment of galactosemia. Often, these are designed to treat the symptoms and not the fundamental cause. Pharmacological chaperones (PC) (small molecules which correct the folding of misfolded proteins) represent an exciting potential therapy for galactosemia. In theory, they would restore enzyme function, thus preventing downstream pathological consequences. In practice, no PCs have been identified for potential application in galactosemia. Here, we review the biochemical basis of the disease, identify opportunities for the application of PCs and describe how these might be discovered. We will conclude by considering some of the clinical issues which will affect the future use of PCs in the treatment of galactosemia.

Identification and prognostic value of a glycolysis-related gene signature in patients with bladder cancer

ABSTRACT

Bladder cancer (BC) is one of the most common malignancies worldwide. Several biomarkers related to the prognosis of patients with BC have previously been identified. However, these prognostic models use only one gene and are thus not reliable or accurate enough. The purpose of our study was to develop an innovative gene signature that has greater prognostic value in BC. So, in this study, we performed mRNA expression profiling of glycolysis-related genes in BC (n = 407) cohorts by mining data from The Cancer Genome Atlas (TCGA) database. The glycolysis-related gene sets were confirmed using the Gene Set Enrichment Analysis (GSEA). Using Cox regression analysis, a risk score staging model was built based on the genes that were determined to be significantly associated with BC outcome. Eventually, the system of risk score was structured to predict a patient's survival, and we identified four genes (CHPF, AK3, GALK1, and NUP188) that were associated with the outcomes of BC patients. According to the above-mentioned gene signature, patients were divided into two risk subgroups. The analysis showed that our constructed risk model was independent of clinical features and that the risk score was a highly powerful tool for predicting the overall survival (OS) of BC patients. Taking together, we identified a gene signature associated with glycolysis that could effectively predict the prognosis of BC patients. Our findings offer a new perspective for the clinical research and treatment of BC.

TK, Ahmad MN, Okashah SS, Kamaraj B, Al-Subaie AM, C. GPD, Zayed H. Deciphering the Role of Filamin B Calponin-Homology Domain in Causing the Larsen Syndrome, Boomerang Dysplasia, and Atelosteogenesis Type I Spectrum Disorders via a Computational Approach

Filamins (FLN) are a family of actin-binding proteins involved in regulating the cytoskeleton and signaling phenomenon by developing a network with F-actin and FLN-binding partners. The FLN family comprises three conserved isoforms in mammals: FLNA, FLNB, and FLNC. FLNB is a multidomain monomer protein with domains containing an actin-binding N-terminal domain (ABD 1-242), encompassing two calponin-homology domains (assigned CH1 and CH2). Primary variants in FLNB mostly occur in the domain (CH2) and surrounding the hinge-1 region. The four autosomal dominant disorders that are associated with FLNB variants are Larsen syndrome, atelosteogenesis type I (AOI), atelosteogenesis type III (AOIII), and boomerang dysplasia (BD). Despite the intense clustering of FLNB variants contributing to the LS-AO-BD disorders, the genotype-phenotype correlation is still enigmatic. In silico prediction tools and molecular dynamics simulation (MDS) approaches have offered the potential for variant classification and pathogenicity predictions. We retrieved 285 FLNB missense variants from the UniProt, ClinVar, and HGMD databases in the current study. Of these, five and 39 variants were located in the CH1 and CH2 domains, respectively. These variants were subjected to various pathogenicity and stability prediction tools, evolutionary and conservation analyses, and biophysical and physicochemical properties analyses. Molecular dynamics simulation (MDS) was performed on the three candidate variants in the CH2 domain (W148R, F161C, and L171R) that were predicted to be the most pathogenic. The MDS analysis results showed that these three variants are highly compact compared to the native protein, suggesting that they could affect the protein on the structural and functional levels. The computational approach demonstrates the differences between the FLNB mutants and the wild type in a structural and functional context. Our findings expand our knowledge on the genotype-phenotype correlation in FLNB-related LS-AO-BD disorders on the molecular level, which may pave the way for optimizing drug therapy by integrating precision medicine.

Comprehensive in silico screening and molecular dynamics studies of missense mutations in Sjogren-Larsson syndrome associated with the ALDH3A2 gene

Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.

Bioinformatics classification of mutations in patients with Mucopolysaccharidosis IIIA

Mucopolysaccharidosis (MPS) IIIA, also known as Sanfilippo syndrome type A, is a severe, progressive disease that affects the central nervous system (CNS). MPS IIIA is inherited in an autosomal recessive manner and is caused by a deficiency in the lysosomal enzyme sulfamidase, which is required for the degradation of heparan sulfate. The sulfamidase is produced by the N-sulphoglucosamine sulphohydrolase (SGSH) gene. In MPS IIIA patients, the excess of lysosomal storage of heparan sulfate often leads to mental retardation, hyperactive behavior, and connective tissue impairments, which occur due to various known missense mutations in the SGSH, leading to protein dysfunction. In this study, we focused on three mutations (R74C, S66W, and R245H) based on in silico pathogenic, conservation, and stability prediction tool studies. The three mutations were further subjected to molecular dynamic simulation (MDS) analysis using GROMACS simulation software to observe the structural changes they induced, and all the mutants exhibited maximum deviation patterns compared with the native protein. Conformational changes were observed in the mutants based on various geometrical parameters, such as conformational stability, fluctuation, and compactness, followed by hydrogen bonding, physicochemical properties, principal component analysis (PCA), and salt bridge analyses, which further validated the underlying cause of the protein instability. Additionally, secondary structure and surrounding amino acid analyses further confirmed the above results indicating the loss of protein function in the mutants compared with the native protein. The present results reveal the effects of three mutations on the enzymatic activity of sulfamidase, providing a molecular explanation for the cause of the disease. Thus, this study allows for a better understanding of the effect of SGSH mutations through the use of various computational approaches in terms of both structure and functions and provides a platform for the development of therapeutic drugs and potential disease treatments.

Dynamic origins of substrate promiscuity in bacterial galactokinases

Galactokinase catalyses the ATP-dependent phosphorylation of galactose and structurally related sugars. The enzyme has attracted interest as a potential biocatalyst for the production of sugar 1-phosphates and several attempts have been made to broaden its specificity. In general, bacterial galactokinases have wider substrate ranges than mammalian ones. The enzymes from Escherichia coli and Lactococcus lactis have received particular attention and a number of variants with increased promiscuity have been identified. Here, we present a molecular dynamics study designed to investigate the molecular causes of the wider substrate ranges of these enzymes and their variants with particular reference to protein mobility. Some regions close to the active site of the enzyme have different structures in the bacterial enzymes compared to the human one. Alterations known to increase the substrate range (e.g. Y371H in the E. coli enzyme), tend to alter the conformation of a key α-helical region (residues 216-232 in the E. coli enzyme). The equivalent helix in the human enzyme has previously been predicted to be altered in variants which affect catalytic activity or protein stability. This helix appears to be a key region in galactokinases from a range of species and may represent an interesting target for future attempts to broaden the specificity of galactokinases.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.

Computational modelling approaches as a potential platform to understand the molecular genetics association between Parkinson's and Gaucher diseases

Gaucher's disease (GD) is a genetic disorder in which glucocerebroside accumulates in cells and specific organs. It is broadly classified into type I, type II and type III. Patients with GD are at high risk of Parkinson's disease (PD), and the clinical and pathological presentation of GD patients with PD is almost identical to idiopathic PD. Several experimental models like cell culture, animal models, and transgenic mice models were used to understand the molecular mechanism behind GD and PD association; however, such mechanism remains unclear. In this context, based on literature reports, we identified the most common mutations K198T, E326K, T369M, N370S, V394L, D409H, L444P, and R496H, in the Glucosylceramidase (GBA) protein that are known to cause GD1, and represent a risk of developing PD. However, to date, no computational analyses have designed to elucidate the potential functional role of GD mutations with increased risk of PD. The present computational pipeline allows us to understand the structural and functional significance of these GBA mutations with PD. Based on the published data, the most common and severe mutations were E326K, N370S, and L444P, which further selected for our computational analysis. PredictSNP and iStable servers predicted L444P mutant to be the most deleterious and responsible for the protein destabilization, followed by the N370S mutation. Further, we used the structural analysis and molecular dynamics approach to compare the most frequent deleterious mutations (N370S and L444P) with the mild mutation E326K. The structural analysis demonstrated that the location of E326K and N370S in the alpha helix region of the protein whereas the mutant L444P was in the starting region of the beta sheet, which might explain the predicted pathogenicity level and destabilization effect of the L444P mutant. Finally, Molecular Dynamics (MD) at 50 ns showed the highest deviation and fluctuation pattern in the L444P mutant compared to the two mutants E326K and N370S and the native protein. This was consistent with more loss of intramolecular hydrogen bonds and less compaction of the radius of gyration in the L444P mutant. The proposed study is anticipated to serve as a potential platform to understand the mechanism of the association between GD and PD, and might facilitate the process of drug discovery against both GD and PD.

Impact of missense mutations in survival motor neuron protein (SMN1) leading to Spinal Muscular Atrophy (SMA): A computational approach

Spinal muscular atrophy (SMA) is a neuromuscular disorder caused by the mutations in survival motor neuron 1 gene (SMN1). The molecular pathology of missense mutations in SMN1 is not thoroughly investigated so far. Therefore, we collected all missense mutations in the SMN1 protein, using all possible search terms, from three databases (PubMed, PMC and Google Scholar). All missense mutations were subjected to in silico pathogenicity, conservation, and stability analysis tools. We used statistical analysis as a QC measure for validating the specificity and accuracy of these tools. PolyPhen-2 demonstrated the highest specificity and accuracy. While PolyPhen-1 showed the highest sensitivity; overall, PolyPhen2 showed better measures in comparison to other in silico tools. Three mutations (D44V, Y272C, and Y277C) were identified as the most pathogenic and destabilizing. Further, we compared the physiochemical properties of the native and the mutant amino acids and observed loss of H-bonds and aromatic stacking upon the cysteine to tyrosine substitution, which led to the loss of aromatic rings and may reduce protein stability. The three mutations were further subjected to Molecular Dynamics Simulation (MDS) analysis using GROMACS to understand the structural changes. The Y272C and Y277C mutants exhibited maximum deviation pattern from the native protein as compared to D44V mutant. Further MDS analysis predicted changes in the stability that may have been contributed due to the loss of hydrogen bonds as observed in intramolecular hydrogen bond analysis and physiochemical analysis. A loss of function/structural impact was found to be severe in the case of Y272C and Y277C mutants in comparison to D44V mutation. Correlating the results from in silico predictions, physiochemical analysis, and MDS, we were able to observe a loss of stability in all the three mutants. This combinatorial approach could serve as a platform for variant interpretation and drug design for spinal muscular dystrophy resulting from missense mutations.

Computational approach to unravel the impact of missense mutations of proteins (D2HGDH and IDH2) causing D-2-hydroxyglutaric aciduria 2

The 2-hydroxyglutaric aciduria (2-HGA) is a rare neurometabolic disorder that leads to the development of brain damage. It is classified into three categories: D-2-HGA, L-2-HGA, and combined D,L-2-HGA. The D-2-HGA includes two subtypes: type I and type II caused by the mutations in D2HGDH and IDH2 proteins, respectively. In this study, we studied six mutations, four in the D2HGDH (I147S, D375Y, N439D, and V444A) and two in the IDH2 proteins (R140G, R140Q). We performed in silico analysis to investigate the pathogenicity and stability changes of the mutant proteins using pathogenicity (PANTHER, PhD-SNP, SIFT, SNAP, and META-SNP) and stability (i-Mutant, MUpro, and iStable) predictors. All the mutations of both D2HGDH and IDH2 proteins were predicted as disease causing except V444A, which was predicted as neutral by SIFT. All the mutants were also predicted to be destabilizing the protein except the mutants D375Y and N439D. DSSP plugin of the PyMOL and Molecular Dynamics Simulations (MDS) were used to study the structural changes in the mutant proteins. In the case of D2HGDH protein, the mutations I147S and V444A that are positioned in the beta sheet region exhibited higher Root Mean Square Deviation (RMSD), decrease in compactness and number of intramolecular hydrogen bonds compared to the mutations N439D and D375Y that are positioned in the turn and loop region, respectively. While the mutants R140Q and R140QG that are positioned in the alpha helix region of the protein. MDS results revealed the mutation R140Q to be more destabilizing (higher RMSD values, decrease in compactness and number of intramolecular hydrogen bonds) compared to the mutation R140G of the IDH2 protein. This study is expected to serve as a platform for drug development against 2-HGA and pave the way for more accurate variant assessment and classification for patients with genetic diseases.