Analysis of galactosemia-linked mutations of GALT enzyme using a computational biology approach

Sam-Sam Association Between EphA2 and SASH1: In Silico Studies of Cancer-Linked Mutations

Recently, SASH1 has emerged as a novel protein interactor of a few Eph tyrosine kinase receptors like EphA2. These interactions involve the first N-terminal Sam (sterile alpha motif) domain of SASH1 (SASH1-Sam1) and the Sam domain of Eph receptors. Currently, the functional meaning of the SASH1-Sam1/EphA2-Sam complex is unknown, but EphA2 is a well-established and crucial player in cancer onset and progression. Thus, herein, to investigate a possible correlation between the formation of the SASH1-Sam1/EphA2-Sam complex and EphA2 activity in cancer, cancer-linked mutations in SASH1-Sam1 were deeply analyzed. Our research plan relied first on searching the COSMIC database for cancer-related SASH1 variants carrying missense mutations in the Sam1 domain and then, through a variety of bioinformatic tools and molecular dynamic simulations, studying how these mutations could affect the stability of SASH1-Sam1 alone, leading eventually to a defective fold. Next, through docking studies, with the support of AlphaFold2 structure predictions, we investigated if/how mutations in SASH1-Sam1 could affect binding to EphA2-Sam. Our study, apart from presenting a solid multistep research protocol to analyze structural consequences related to cancer-associated protein variants with the support of cutting-edge artificial intelligence tools, suggests a few mutations that could more likely modulate the interaction between SASH1-Sam1 and EphA2-Sam.

Clinical and biochemical phenotypes, genotypes, and long-term outcomes of individuals with galactosemia type I from a single metabolic genetics center in Alberta

Background

Galactosemia type I is an autosomal recessive disorder of galactose metabolism due to galactose-1-phosphate uridyltransferase deficiency, encoded by GALT. To investigate the phenotypes, genotypes and long-term outcomes of galactosemia, we performed a retrospective cohort study in our center.

Methods

All individuals with galactosemia type I were included. We divided individuals into two groups to compare the outcomes of those treated symptomatically (SymX) and asymptomatically (AsymX). We reviewed electronic patient charts for clinical features, biochemical investigations, molecular genetic investigations, treatments, and outcomes.

Results

There were 25 individuals including classic (n = 17), clinical variant (n = 4), and biochemical variant (Duarte) galactosemia (n = 4). Twelve individuals were diagnosed symptomatically (SymX), and 9 individuals were diagnosed asymptomatically (AsymX). We did not include individuals with biochemical variant (Duarte) galactosemia into any of these groups. At the time of the diagnosis, conjugated hyperbilirubinemia was present in 83.3% of SymX group, whereas only 22% of AsymX group. SymX group had hepatomegaly (25%), failure to thrive (33.3%), cataract (16.7%) and sepsis (25%), whereas none of the individuals in the AsymX group had these clinical features. Fourteen variants in GALT were identified including pathogenic/likely pathogenic (n = 12), and likely benign/benign (n = 2) variants. The vast majority of individuals with classic and clinical variant galactosemia were treated with a galactose-lactose-free diet for life (n = 20/21). Intellectual disability was present in 54.5% of the SymX group, and in 37.5% of the AsymX group as a long-term outcome. Tremors were present 50% of the SymX group, and in 22% of the AsymX group as a long-term outcome. Although, intellectual disability and tremors seem to be less common in the AsymX group, there was no statistically significant difference between both groups. Primary ovarian insufficiency was present 50% of the SymX group, whereas in 20% of the AsymX group in post-pubertal females. We report a novel hypomorphic GALT variant (p.Ala303Ser) in one individual with clinical variant galactosemia. We also report an individual with clinical variant galactosemia with normal urine galactitol levels on a normal diet.

Conclusion

It seems that newborn screening and early administration of a galactose-lactose-free diet decreases the long-term galactosemia-associated complications but does not prevent them completely. It may be that not all individuals with clinical variant galactosemia may need a galactose-lactose-free diet. It is timely to find new therapeutic strategies that can reduce the frequency of late-onset complications in galactosemia.

Cancer-Related Mutations in the Sam Domains of EphA2 Receptor and Ship2 Lipid Phosphatase: A Computational Study

The lipid phosphatase Ship2 interacts with the EphA2 receptor by forming a heterotypic Sam (sterile alpha motif)-Sam complex. Ship2 works as a negative regulator of receptor endocytosis and consequent degradation, and anti-oncogenic effects in cancer cells should be induced by hindering its association with EphA2. Herein, a computational approach is presented to investigate the relationship between Ship2-Sam/EphA2-Sam interaction and cancer onset and further progression. A search was first conducted through the COSMIC (Catalogue of Somatic Mutations in Cancer) database to identify cancer-related missense mutations positioned inside or close to the EphA2-Sam and Ship2-Sam reciprocal binding interfaces. Next, potential differences in the chemical-physical properties of mutant and wild-type Sam domains were evaluated by bioinformatics tools based on analyses of primary sequences. Three-dimensional (3D) structural models of mutated EphA2-Sam and Ship2-Sam domains were built as well and deeply analysed with diverse computational instruments, including molecular dynamics, to classify potentially stabilizing and destabilizing mutations. In the end, the influence of mutations on the EphA2-Sam/Ship2-Sam interaction was studied through docking techniques. This in silico approach contributes to understanding, at the molecular level, the mutation/cancer relationship by predicting if amino acid substitutions could modulate EphA2 receptor endocytosis.

Classic Galactosemia: Clinical and Computational Characterization of a Novel GALT Missense Variant (p.A303D) and a Literature Review

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.

Comparison of In Vitro and In Silico Assessments of Human Galactose-1-Phosphate Uridylyltransferase Coding Variants

Introduction Human pathogenic coding variations of the galactose-1-phosphate uridylyltransferase (GALT) gene cause classic galactosemia, a recessive disease of galactose metabolism. Unfortunately, there are many variants of uncertain significance (VUS) that need to be characterized in order to be able to predict the likelihood of classic galactosemia for all possible genotypes. There are many bioinformatic resources available that attempt to predict the pathogenicity of a human variant, but it is unclear if these methods realistically predict the consequence of these variants. To determine the clinical application of these resources, we compared the results of in vitro enzymatic assays with in silico predictive models. Methods In all assays, we compared the activity of the three human GALT VUS (Alanine81Threonine, Histidine47Aspartate, Glutamate58Lysine) to native GALT (nGALT) and to a variant of known pathogenic clinical significance (Glutamine188Arginine). The enzymatic activities of VUS recombinant proteins were compared to the results of in silico analytical methods. The in silico methods included the comparison of molecular dynamic simulation root-mean-square deviation (RMSD) results and the results from predictive programs PredictSNP, evolutionary model of variant effect (EVE), ConSurf, and sorting intolerant from tolerant (SIFT). Results The enzymatic assays showed that the variants tested had diminished Vmax relative to the native protein. The VUS RMSD data for both the whole protein and individual residues in the molecular dynamics simulations were not significantly different when compared to nGALT. The other predictive programs had mixed results for each VUS and were not consistent with the enzyme activity or simulation results. Conclusions Our experiments indicated a statistically significant decrease in enzymatic activity of the VUS when compared to nGALT. These experiments also demonstrated significant differences between in silico predictions and in vitro results. These results suggest that the in silico tools used may not be beneficial in determining the pathogenicity of GALT VUS.

Pan-phylum In Silico Analyses of Nematode Endocannabinoid Signalling Systems Highlight Novel Opportunities for Parasite Drug Target Discovery

The endocannabinoid signalling (ECS) system is a complex lipid signalling pathway that modulates diverse physiological processes in both vertebrate and invertebrate systems. In nematodes, knowledge of endocannabinoid (EC) biology is derived primarily from the free-living model species Caenorhabditis elegans, where ECS has been linked to key aspects of nematode biology. The conservation and complexity of nematode ECS beyond C. elegans is largely uncharacterised, undermining the understanding of ECS biology in nematodes including species with key importance to human, veterinary and plant health. In this study we exploited publicly available omics datasets, in silico bioinformatics and phylogenetic analyses to examine the presence, conservation and life stage expression profiles of EC-effectors across phylum Nematoda. Our data demonstrate that: (i) ECS is broadly conserved across phylum Nematoda, including in therapeutically and agriculturally relevant species; (ii) EC-effectors appear to display clade and lifestyle-specific conservation patterns; (iii) filarial species possess a reduced EC-effector complement; (iv) there are key differences between nematode and vertebrate EC-effectors; (v) life stage-, tissue- and sex-specific EC-effector expression profiles suggest a role for ECS in therapeutically relevant parasitic nematodes. To our knowledge, this study represents the most comprehensive characterisation of ECS pathways in phylum Nematoda and inform our understanding of nematode ECS complexity. Fundamental knowledge of nematode ECS systems will seed follow-on functional studies in key nematode parasites to underpin novel drug target discovery efforts.

Investigating the Effects of Amino Acid Variations in Human Menin

Human menin is a nuclear protein that participates in many cellular processes, as transcriptional regulation, DNA damage repair, cell signaling, cell division, proliferation, and migration, by interacting with many other proteins. Mutations of the gene encoding menin cause multiple endocrine neoplasia type 1 (MEN1), a rare autosomal dominant disorder associated with tumors of the endocrine glands. In order to characterize the structural and functional effects at protein level of the hundreds of missense variations, we investigated by computational methods the wild-type menin and more than 200 variants, predicting the amino acid variations that change secondary structure, solvent accessibility, salt-bridge and H-bond interactions, protein thermostability, and altering the capability to bind known protein interactors. The structural analyses are freely accessible online by means of a web interface that integrates also a 3D visualization of the structure of the wild-type and variant proteins. The results of the study offer insight into the effects of the amino acid variations in view of a more complete understanding of their pathological role.

Simulation of the Interactions of Arginine with Wild-Type GALT Enzyme and the Classic Galactosemia-Related Mutant p.Q188R by a Computational Approach

Classic galactosemia is an inborn error of metabolism associated with mutations that impair the activity and the stability of galactose-1-phosphate uridylyltransferase (GALT), catalyzing the third step in galactose metabolism. To date, no treatments (including dietary galactose deprivation) are able to prevent or alleviate the long-term complications affecting galactosemic patients. Evidence that arginine is able to improve the activity of the human enzyme expressed in a prokaryotic model of classic galactosemia has induced researchers to suppose that this amino acid could act as a pharmacochaperone, but no effects were detected in four galactosemic patients treated with this amino acid. Given that no molecular characterizations of the possible effects of arginine on GALT have been performed, and given that the samples of patients treated with arginine are extremely limited for drawing definitive conclusions at the clinical level, we performed computational simulations in order to predict the interactions (if any) between this amino acid and the enzyme. Our results do not support the possibility that arginine could function as a pharmacochaperone for GALT, but information obtained by this study could be useful for identifying, in the future, possible pharmacochaperones for this enzyme.

Analysis of the Structure-Function-Dynamics Relationships of GALT Enzyme and of Its Pathogenic Mutant p.Q188R: A Molecular Dynamics Simulation Study in Different Experimental Conditions

The third step of the catabolism of galactose in mammals is catalyzed by the enzyme galactose-1-phosphate uridylyltransferase (GALT), a homodimeric enzyme with two active sites located in the proximity of the intersubunit interface. Mutations of this enzyme are associated to the rare inborn error of metabolism known as classic galactosemia; in particular, the most common mutation, associated with the most severe phenotype, is the one that replaces Gln188 in the active site of the enzyme with Arg (p.Gln188Arg). In the past, and more recently, the structural effects of this mutation were deduced on the static structure of the wild-type human enzyme; however, we feel that a dynamic view of the proteins is necessary to deeply understand their behavior and obtain tips for possible therapeutic interventions. Thus, we performed molecular dynamics simulations of both wild-type and p.Gln188Arg GALT proteins in the absence or in the presence of the substrates in different conditions of temperature. Our results suggest the importance of the intersubunit interactions for a correct activity of this enzyme and can be used as a starting point for the search of drugs able to rescue the activity of this enzyme in galactosemic patients.

Structural Dissection of Viral Spike-Protein Binding of SARS-CoV-2 and SARS-CoV-1 to the Human Angiotensin-Converting Enzyme 2 (ACE2) as Cellular Receptor

An outbreak by a new severe acute respiratory syndrome betacoronavirus (SARS-CoV-2) has spread CoronaVirus Disease 2019 (COVID-19) all over the world. Immediately, following studies have confirmed the human Angiotensin-Converting Enzyme 2 (ACE2) as a cellular receptor of viral Spike-Protein (Sp) that mediates the CoV-2 invasion into the pulmonary host cells. Here, we compared the molecular interactions of the viral Sp from previous SARS-CoV-1 of 2002 and SARS-CoV-2 with the host ACE2 protein by in silico analysis of the available experimental structures of Sp-ACE2 complexes. The K417 amino acid residue, located in the region of Sp Receptor-Binding Domain (RBD) of the new coronavirus SARS-CoV-2, showed to have a key role for the binding to the ACE2 N-terminal region. The R426 residue of SARS-CoV-1 Sp-RBD also plays a key role, although by interacting with the central region of the ACE2 sequence. Therefore, our study evidenced peculiarities in the interactions of the two Sp-ACE2 complexes. Our outcomes were consistent with previously reported mutagenesis studies on SARS-CoV-1 and support the idea that a new and different RBD was acquired by SARS-CoV-2. These results have interesting implications and suggest further investigations.

Umbilical cord blood-based gene signatures related to prenatal major depressive disorder

BACKGROUND

Prenatal exposure to depression has been considered as a risk factor for adverse childhood, while it is accompanied by unknown molecular mechanisms. The aim of this study was to identify differentially expressed genes (DEGs) and associated biological processes between cord blood samples from neonates born to mothers who exposed to major depressive disorder (MDD) and healthy mothers.

METHODS

The microarray data GSE114852 were downloaded to analyze the mRNA expression profiles of umbilical cord blood with 31 samples exposed to prenatal MDD and 62 samples with healthy mothers. Kyoto Encyclopedia of Genes and Genomes pathway and Gene ontology enrichment analyses were conducted to identify associated biochemical pathways and functional categories of the DEGs. The protein-protein interaction network was constructed and the top 10 hub genes in the network were predicted.

RESULTS

The results showed several immunity related processes, such as "phagosome", "Epstein-Barr virus infection", "proteasome", "positive regulation of I-kappaB kinase/NF-kappaB signaling", "interferon-gamma-mediated signaling pathway", and "tumor necrosis factor" presented significant differences between two groups. Most of the hub genes (for example PSMD2, PSMD6, PSMB8, PSMB9) were also associated with immune pathways.

CONCLUSION

This bioinformatic analysis demonstrated immune-mediated mechanisms might play a fatal role in abnormalities in fetal gene expression profiles caused by prenatal MDD.

Effect of genotype on galactose-1-phosphate in classic galactosemia patients

Impaired activity of galactose-1-phosphate uridyltransferase (GALT) causes classic galactosemia (OMIM 230400), characterized by the accumulation of galactose-1-phosphate (GAL1P) in patients' red blood cells (RBCs). Our recent study demonstrated a correlation between RBC GAL1P and long-term outcomes in galactosemia patients. Here, we analyze biochemical and molecular results in 77 classic galactosemia patients to evaluate the association between GALT genotypes and GAL1P concentration in RBCs. Experimental data from model organisms were also included to assess the correlation between GAL1P and predicted residual activity of each genotype. Although all individuals in this study showed markedly reduced RBC GALT activity, we observed significant differences in RBC GAL1P concentrations among galactosemia genotypes. While levels of GAL1P on treatment did not correlate with RBC GALT activities (p = 0.166), there was a negative nonlinear correlation between mean GAL1P concentrations and predicted residual enzyme activity of genotype (p = 0.004). These studies suggest that GAL1P levels in RBCs on treatment likely reflect the overall functional impairment of GALT in patients with galactosemia.

Galactose-1-phosphate uridyltransferase deficiency: A literature review of the putative mechanisms of short and long-term complications and allelic variants

Galactosemia type 1 is an autosomal recessive disorder of galactose metabolism, determined by a deficiency in the enzyme galactose-1-phosphate uridyltransferase (GALT). GALT deficiency is classified as severe or variant depending on biochemical phenotype, genotype and potential to develop acute and long-term complications. Neonatal symptoms usually resolve after galactose-restricted diet; however, some patients, despite the diet, can develop long-term complications, in particular when the GALT enzyme activity results absent or severely decreased. The mechanisms of acute and long-term complications are still discussed and several hypotheses are presented in the literature like enzymatic inhibition, osmotic stress, endoplasmic reticulum stress, oxidative stress, defects of glycosylation or epigenetic modification. This review summarizes the current knowledge of galactosemia, in particular the putative mechanisms of neonatal and long-term complications and the molecular genetics of GALT deficiency.

Sweet and sour: an update on classic galactosemia

Classic galactosemia is a rare inherited disorder of galactose metabolism caused by deficient activity of galactose-1-phosphate uridylyltransferase (GALT), the second enzyme of the Leloir pathway. It presents in the newborn period as a life-threatening disease, whose clinical picture can be resolved by a galactose-restricted diet. The dietary treatment proves, however, insufficient in preventing severe long-term complications, such as cognitive, social and reproductive impairments. Classic galactosemia represents a heavy burden on patients' and their families' lives. After its first description in 1908 and despite intense research in the past century, the exact pathogenic mechanisms underlying galactosemia are still not fully understood. Recently, new important insights on molecular and cellular aspects of galactosemia have been gained, and should open new avenues for the development of novel therapeutic strategies. Moreover, an international galactosemia network has been established, which shall act as a platform for expertise and research in galactosemia. Herein are reviewed some of the latest developments in clinical practice and research findings on classic galactosemia, an enigmatic disorder with many unanswered questions warranting dedicated research.

In silico prediction of the effects of mutations in the human triose phosphate isomerase gene: Towards a predictive framework for TPI deficiency

Triose phosphate isomerase (TPI) deficiency is a rare, but highly debilitating, inherited metabolic disease. Almost all patients suffer severe neurological effects and the most severely affected are unlikely to live beyond early childhood. Here, we describe an in silico study into well-characterised variants which are associated with the disease alongside an investigation into 79 currently uncharacterised TPI variants which are known to occur in the human population. The majority of the disease-associated mutations affected amino acid residues close to the dimer interface or the active site. However, the location of the altered amino acid residue did not predict the severity of the resulting disease. Prediction of the effect on protein stability using a range of different programs suggested a relationship between the degree of instability caused by the sequence variation and the severity of the resulting disease. Disease-associated variations tended to affect well-conserved residues in the protein's sequence. However, the degree of conservation of the residue was not predictive of disease severity. The majority of the 79 uncharacterised variants are potentially associated with disease since they were predicted to destabilise the protein and often occur in well-conserved residues. We predict that individuals homozygous for the corresponding mutations would be likely to suffer from TPI deficiency.

Molecular basis of classic galactosemia from the structure of human galactose 1-phosphate uridylyltransferase

Classic galactosemia is a potentially lethal disease caused by the dysfunction of galactose 1-phosphate uridylyltransferase (GALT). Over 300 disease-associated GALT mutations have been reported, with the majority being missense changes, although a better understanding of their underlying molecular effects has been hindered by the lack of structural information for the human enzyme. Here, we present the 1.9 Å resolution crystal structure of human GALT (hGALT) ternary complex, revealing a homodimer arrangement that contains a covalent uridylylated intermediate and glucose-1-phosphate in the active site, as well as a structural zinc-binding site, per monomer. hGALT reveals significant structural differences from bacterial GALT homologues in metal ligation and dimer interactions, and therefore is a zbetter model for understanding the molecular consequences of disease mutations. Both uridylylation and zinc binding influence the stability and aggregation tendency of hGALT. This has implications for disease-associated variants where p.Gln188Arg, the most commonly detected, increases the rate of aggregation in the absence of zinc likely due to its reduced ability to form the uridylylated intermediate. As such our structure serves as a template in the future design of pharmacological chaperone therapies and opens new concepts about the roles of metal binding and activity in protein misfolding by disease-associated mutants.

Value of predictive bioinformatics in inherited metabolic diseases

Typically, inherited metabolic diseases arise from point mutations in genes encoding metabolic enzymes. Although some of these mutations directly affect amino acid residues in the active sites of these enzymes, the majority do not. It is now well accepted that the majority of these disease-associated mutations exert their effects through alteration of protein stability, which causes a reduction in enzymatic activity. This finding suggests a way to predict the severity of newly discovered mutations. In silico prediction of the effects of amino acid sequence alterations on protein stability often correlates with disease severity. However, no stability prediction tool is perfect and, in general, better results are obtained if the predictions from a variety of tools are combined and then interpreted. In addition to predicted alterations to stability, the degree of conservation of a particular residue can also be a factor which needs to be taken into account: alterations to highly conserved residues are more likely to be associated with severe forms of the disease. The approach has been successfully applied in a variety of inherited metabolic diseases, but further improvements are necessary to enable robust translation into clinically useful tools.

Clinical and molecular spectra in galactosemic patients from neonatal screening in northeastern Italy: structural and functional characterization of new variations in the galactose-1-phosphate uridyltransferase (GALT) gene

Classical galactosemia is an autosomal recessive inborn error of metabolism due to mutations of the GALT gene leading to toxic accumulation of galactose and derived metabolites. With the benefit of early diagnosis by neonatal screening and early therapy, the acute presentation of classical galactosemia can be prevented. However, despite early diagnosis and treatment, the long term outcome for these patients is still unpredictable because they may go on to develop cognitive disability, speech problems, neurological and/or movement disorders and, in females, ovarian dysfunction. The objectives of the current study were to report our experience with a group of galactosemic patients identified through the neonatal screening programs in northeastern Italy during the last 30years. No neonatal deaths due to galactosemia complications occurred after the introduction of the neonatal screening program. However, despite the early diagnosis and dietary treatment, the patients with classical galactosemia showed one or more long-term complications. A total of 18 different variations in the GALT gene were found in the patient cohort: 12 missense, 2 frameshift, 1 nonsense, 1 deletion, 1 silent variation, and 1 intronic. Six (p.R33P, p.G83V, p.P244S, p.L267R, p.L267V, p.E271D) were new variations. The most common variation was p.Q188R (12 alleles, 31.5%), followed by p.K285N (6 alleles, 15.7%) and p.N314D (6 alleles, 15.7%). The other variations comprised 1 or 2 alleles. In the patients carrying a new mutation, the biochemical analysis of GALT activity in erythrocytes showed an activity of <1%. In silico analysis (SIFT, PolyPhen-2 and the computational analysis on the static protein structure) showed potentially damaging effects of the six new variations on the GALT protein, thus expanding the genetic spectrum of GALT variations in Italy. The study emphasizes the difficulty in establishing a genotype-phenotype correlation in classical galactosemia and underlines the importance of molecular diagnostic testing prior to making any treatment.

GALT protein database: querying structural and functional features of GALT enzyme

Knowledge of the impact of variations on protein structure can enhance the comprehension of the mechanisms of genetic diseases related to that protein. Here, we present a new version of GALT Protein Database, a Web-accessible data repository for the storage and interrogation of structural effects of variations of the enzyme galactose-1-phosphate uridylyltransferase (GALT), the impairment of which leads to classic Galactosemia, a rare genetic disease. This new version of this database now contains the models of 201 missense variants of GALT enzyme, including heterozygous variants, and it allows users not only to retrieve information about the missense variations affecting this protein, but also to investigate their impact on substrate binding, intersubunit interactions, stability, and other structural features. In addition, it allows the interactive visualization of the models of variants collected into the database. We have developed additional tools to improve the use of the database by nonspecialized users. This Web-accessible database (http://bioinformatica.isa.cnr.it/GALT/GALT2.0) represents a model of tools potentially suitable for application to other proteins that are involved in human pathologies and that are subjected to genetic variations.

Functional and structural impact of the most prevalent missense mutations in classic galactosemia

Galactose-1-phosphate uridylyltransferase (GALT) is a key enzyme in galactose metabolism, particularly important in the neonatal period due to ingestion of galactose-containing milk. GALT deficiency results in the genetic disorder classic galactosemia, whose pathophysiology is still not fully elucidated. Whereas classic galactosemia has been hypothesized to result from GALT misfolding, a thorough functional-structural characterization of GALT most prevalent variants was still lacking, hampering the development of alternative therapeutic approaches. The aim of this study was to investigate the structural-functional effects of nine GALT mutations, four of which account for the vast majority of the mutations identified in galactosemic patients. Several methodologies were employed to evaluate the mutations' impact on GALT function, on the protein secondary and tertiary structures, and on the aggregation propensity. The major structural effect concerns disturbed propensity for aggregation, particularly striking for the p.Q188R variant, resulting from the most frequent (∼60%) allele at a worldwide scale. The absence of major effects at the secondary and tertiary structure levels suggests that the disturbed aggregation results from subtle perturbations causing a higher and/or longer exposure of hydrophobic residues in the variants as compared to WT GALT. The results herein described indicate a possible benefit from introducing proteostasis regulators and/or chemical/pharmacological chaperones to prevent the accumulation of protein aggregates, in new avenues of therapeutic research for classic galactosemia.

In silico prediction of the effects of mutations in the human UDP-galactose 4'-epimerase gene: towards a predictive framework for type III galactosemia

The enzyme UDP-galactose 4'-epimerase (GALE) catalyses the reversible epimerisation of both UDP-galactose and UDP-N-acetyl-galactosamine. Deficiency of the human enzyme (hGALE) is associated with type III galactosemia. The majority of known mutations in hGALE are missense and private thus making clinical guidance difficult. In this study a bioinformatics approach was employed to analyse the structural effects due to each mutation using both the UDP-glucose and UDP-N-acetylglucosamine bound structures of the wild-type protein. Changes to the enzyme's overall stability, substrate/cofactor binding and propensity to aggregate were also predicted. These predictions were found to be in good agreement with previous in vitro and in vivo studies when data was available and allowed for the differentiation of those mutants that severely impair the enzyme's activity against UDP-galactose. Next this combination of techniques were applied to another twenty-six reported variants from the NCBI dbSNP database that have yet to be studied to predict their effects. This identified p.I14T, p.R184H and p.G302R as likely severely impairing mutations. Although severely impaired mutants were predicted to decrease the protein's stability, overall predicted stability changes only weakly correlated with residual activity against UDP-galactose. This suggests other protein functions such as changes in cofactor and substrate binding may also contribute to the mechanism of impairment. Finally this investigation shows that this combination of different in silico approaches is useful in predicting the effects of mutations and that it could be the basis of an initial prediction of likely clinical severity when new hGALE mutants are discovered.

Misfolding of galactose 1-phosphate uridylyltransferase can result in type I galactosemia

Type I galactosemia is a genetic disorder that is caused by the impairment of galactose-1-phosphate uridylyltransferase (GALT; EC 2.7.7.12). Although a large number of mutations have been detected through genetic screening of the human GALT (hGALT) locus, for many it is not known how they cause their effects. The majority of these mutations are missense, with predicted substitutions scattered throughout the enzyme structure and thus causing impairment by other means rather than direct alterations to the active site. To clarify the fundamental, molecular basis of hGALT impairment we studied five disease-associated variants p.D28Y, p.L74P, p.F171S, p.F194L and p.R333G using both a yeast model and purified, recombinant proteins. In a yeast expression system there was a correlation between lysate activity and the ability to rescue growth in the presence of galactose, except for p.R333G. Kinetic analysis of the purified proteins quantified each variant's level of enzymatic impairment and demonstrated that this was largely due to altered substrate binding. Increased surface hydrophobicity, altered thermal stability and changes in proteolytic sensitivity were also detected. Our results demonstrate that hGALT requires a level of flexibility to function optimally and that altered folding is the underlying reason of impairment in all the variants tested here. This indicates that misfolding is a common, molecular basis of hGALT deficiency and suggests the potential of pharmacological chaperones and proteostasis regulators as novel therapeutic approaches for type I galactosemia.

Myosin-binding protein C DNA variants in domestic cats (A31P, A74T, R820W) and their association with hypertrophic cardiomyopathy

BACKGROUND

Two mutations in the MYBPC3 gene have been identified in Maine Coon (MCO) and Ragdoll (RD) cats with hypertrophic cardiomyopathy (HCM).

OBJECTIVE

This study examined the frequency of these mutations and of the A74T polymorphism to describe their worldwide distribution and correlation with echocardiography.

ANIMALS

1855 cats representing 28 breeds and random-bred cats worldwide, of which 446 underwent echocardiographic examination.

METHODS

This is a prospective cross-sectional study. Polymorphisms were genotyped by Illumina VeraCode GoldenGate or by direct sequencing. The disease status was defined by echocardiography according to established guidelines. Odds ratios for the joint probability of having HCM and the alleles were calculated by meta-analysis. Functional analysis was simulated.

RESULTS

The MYBPC3 A31P and R820W were restricted to MCO and RD, respectively. Both purebred and random-bred cats had HCM and the incidence increased with age. The A74T polymorphism was not associated with any phenotype. HCM was most prevalent in MCO homozygote for the A31P mutation and the penetrance increased with age. The penetrance of the heterozygote genotype was lower (0.08) compared with the P/P genotype (0.58) in MCO.

CONCLUSIONS AND CLINICAL IMPORTANCE

A31P mutation occurs frequently in MCO cats. The high incidence of HCM in homozygotes for the mutation supports the causal nature of the A31P mutation. Penetrance is incomplete for heterozygotes at A31P locus, at least at a young age. The A74T variant does not appear to be correlated with HCM.

Mutation spectrum in the French cohort of galactosemic patients and structural simulation of 27 novel missense variations

BACKGROUND

Classic galactosemia refers to galactose-1-phosphate uridyltransferase (GALT) deficiency and is characterized by long-term complications of unknown mechanism and high allelic heterogeneity of GALT gene.

AIM

To report molecular characterization of GALT variations in 210 French families, to analyze the structural effects of novel missense variations and to assess informativity of structural data in predicting outcome.

METHODS

Sequencing of exons and intron-exon junctions of GALT gene was completed in unsolved cases by analysis of a long range PCR product. Structural consequences of novel missense variations were predicted using a homology model of GALT protein and a semi-automated analysis which integrates simulation of variations, structural analyses and two web servers dedicated to identify mutation-induced change of protein stability.

RESULTS

Forty four novel variations were identified, among them 27 nucleotide substitutions. In silico modeling of these missense variations showed that 12 variations are predicted to impair subunit interactions and/or active site conformation and that 23 variations modify H-bond or salt-bridge networks. Twenty variations decrease the global stability of the protein. Five variations had apparently no structural effect.

CONCLUSION

Our results expand the mutation spectrum in GALT gene and the list of GALT variations analyzed at the structural level, providing new data to assess the pathophysiology of galactosemia.

Correlation assessment among clinical phenotypes, expression analysis and molecular modeling of 14 novel variations in the human galactose-1-phosphate uridylyltransferase gene

Galactose-1-phosphate uridylyltransferase (GALT) catalyzes the conversion of galactose-1-phosphate to UDP-galactose, a key step in the galactose metabolism. Deficiency of GALT activity in humans caused by deleterious variations in the GALT gene can cause a potentially lethal disease called classic galactosemia. In this study, we selected 14 novel nucleotide sequence changes in the GALT genes found in galactosemic patients for expression analysis and molecular modeling. Several variants showed decreased levels of expression and decreased abundance in the soluble fraction of the Escherichia coli cell extracts, suggesting altered stability and solubility. Only six variant GALT enzymes had detectable enzymatic activities. Kinetic studies showed that their V(max) decreased significantly. To further characterize the variants at molecular level, we performed static and dynamic molecular modeling studies. Effects of variations on local and/or global structural features of the enzyme were anticipated for the majority of variants. In-depth studies with molecular dynamic simulations on selected variants predicted the alteration of the protein structure even though static models apparently did not highlight any perturbation. Overall, these studies offered new insights on the molecular properties of GALT enzyme, with the aim of correlating them with the clinical outcome. Hum Mutat 33:1107-1115, 2012. © 2012 Wiley Periodicals, Inc.

PoPMuSiC 2.1: a web server for the estimation of protein stability changes upon mutation and sequence optimality

Background

The rational design of modified proteins with controlled stability is of extreme importance in a whole range of applications, notably in the biotechnological and environmental areas, where proteins are used for their catalytic or other functional activities. Future breakthroughs in medical research may also be expected from an improved understanding of the effect of naturally occurring disease-causing mutations on the molecular level.

Results

PoPMuSiC-2.1 is a web server that predicts the thermodynamic stability changes caused by single site mutations in proteins, using a linear combination of statistical potentials whose coefficients depend on the solvent accessibility of the mutated residue. PoPMuSiC presents good prediction performances (correlation coefficient of 0.8 between predicted and measured stability changes, in cross validation, after exclusion of 10% outliers). It is moreover very fast, allowing the prediction of the stability changes resulting from all possible mutations in a medium size protein in less than a minute. This unique functionality is user-friendly implemented in PoPMuSiC and is particularly easy to exploit. Another new functionality of our server concerns the estimation of the optimality of each amino acid in the sequence, with respect to the stability of the structure. It may be used to detect structural weaknesses, i.e. clusters of non-optimal residues, which represent particularly interesting sites for introducing targeted mutations. This sequence optimality data is also expected to have significant implications in the prediction and the analysis of particular structural or functional protein regions. To illustrate the interest of this new functionality, we apply it to a dataset of known catalytic sites, and show that a much larger than average concentration of structural weaknesses is detected, quantifying how these sites have been optimized for function rather than stability.

Conclusion

The freely available PoPMuSiC-2.1 web server is highly useful for identifying very rapidly a list of possibly relevant mutations with the desired stability properties, on which subsequent experimental studies can be focused. It can also be used to detect sequence regions corresponding to structural weaknesses, which could be functionally important or structurally delicate regions, with obvious applications in rational protein design.

Structural assessment of the effects of amino acid substitutions on protein stability and protein protein interaction

A structure-based approach is described for predicting the effects of amino acid substitutions on protein function. Structures were predicted using a homology modelling method. Folding and binding energy differences between wild-type and mutant structures were computed to quantitatively assess the effects of amino acid substitutions on protein stability and protein protein interaction, respectively. We demonstrated that pathogenic mutations at the interaction interface could affect binding energy and destabilise protein complex, whereas mutations at the non-interface might reduce folding energy and destabilise monomer structure. The results suggest that the structure-based analysis can provide useful information for understanding the molecular mechanisms of diseases.