Homology modeling studies on human galactose-1-phosphate uridylyltransferase and on its galactosemia-related mutant Q188R provide an explanation of molecular effects of the mutation on homo- and heterodimers

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.

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.

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.

Structure and Ligands Interactions of Citrus Tryptophan Decarboxylase by Molecular Modeling and Docking Simulations

In a previous work, we in silico annotated protein sequences of Citrus genus plants as putative tryptophan decarboxylase (pTDC). Here, we investigated the structural properties of Citrus pTDCs by using the TDC sequence of Catharanthus roseus as an experimentally annotated reference to carry out comparative modeling and substrate docking analyses. The functional annotation as TDC was verified by combining 3D molecular modeling and docking simulations, evidencing the peculiarities and the structural similarities with C. roseus TDC. Docking with l-tryptophan as a ligand showed specificity of pTDC for this substrate. These combined results confirm our previous in silico annotation of the examined protein sequences of Citrus as TDC and provide support for TDC activity in this plant genus.

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.

Uncommon functional properties of the first piscine 26S proteasome from the Antarctic notothenioid Trematomus bernacchii

The antioxidant defense mechanisms have a great impact on the life of Antarctic organisms. The present study could represent the first evidence of a direct involvement of the 26S proteasome in the antioxidant defense systems of fish adapted to cold.

Protein homoeostasis is a fundamental process allowing the preservation of functional proteins and it has a great impact on the life of the Antarctic organisms. However, the effect of low temperatures on protein turnover is poorly understood and the cold-adaptation of the degradation machinery remains an unresolved issue. As the 26S proteasome represents the main proteolytic system devoted to the controlled degradation of intracellular proteins, the purpose of the present study was to investigate the functions of this complex in the notothenioid Trematomus bernacchii, in order to better understand its role in the physiology of Antarctic fish. To this aim, we purified and characterized the 26S proteasome from T. bernacchii and isolated the cDNAs codifying seven of the 14 subunits belonging to the proteasome 20S core particle. Results provided evidences of the high resistance of the piscine 26S proteasome to oxidative agents and of its ‘uncommon’ ability to efficiently hydrolyse oxidized bovine serum albumin (BSA), suggesting that this enzymatic complex could play a key role in the antioxidant defense systems in fish inhabiting permanently cold marine environments. These unique properties were also reflected by the 3D model analysis, which revealed a higher structural stability of the piscine complex respect to the murine template. Finally, a comparative analysis, performed in a variety of tissues collected from T. bernacchii and the temperate fish Dicentrarchus labrax, showed a lower protein retention in the cold-adapted fish, possibly due to a better efficiency of its degradation machinery.

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.

A New APEH Cluster with Antioxidant Functions in the Antarctic Hemoglobinless Icefish Chionodraco hamatus

Acylpeptide hydrolase (APEH) is a ubiquitous cytosolic protease that plays an important role in the detoxification of oxidised proteins. In this work, to further explore the physiological role of this enzyme, two apeh cDNAs were isolated from the Chionodraco hamatus icefish, which lives in the highly oxygenated Antarctic marine environment. The encoded proteins (APEH-1(Ch) and APEH-2(Ch)) were characterised in comparison with the uniquely expressed isoform from the temperate fish Dicentrarchus labrax (APEH-1Dl) and the two APEHs from the red-blooded Antarctic fish Trematomus bernacchii (APEH-1(Tb) and APEH-2(Tb)). Homology modelling and kinetic characterisation of the APEH isoforms provided new insights into their structure/function properties. APEH-2 isoforms were the only ones capable of hydrolysing oxidised proteins, with APEH-2(Ch) being more efficient than APEH-2(Tb) at this specific function. Therefore, this ability of APEH-2 isoforms is the result of an evolutionary adaptation due to the pressure of a richly oxygenated environment. The lack of expression of APEH-2 in the tissues of the temperate fish used as the controls further supported this hypothesis. In addition, analysis of gene expression showed a significant discrepancy between the levels of transcripts and those of proteins, especially for apeh-2 genes, which suggests the presence of post-transcriptional regulation mechanisms, triggered by cold-induced oxidative stress, that produce high basal levels of APEH-2 mRNA as a reserve that is ready to use in case of the accumulation of oxidised proteins.

Biochemical and computational analyses of two phenotypically related GALT mutations (S222N and S135L) that lead to atypical galactosemia

Galactosemia is a metabolic disorder caused by mutations in the GALT gene [1,2]. We encountered a patient heterozygous for a known pathogenic H132Q mutation and a novel S222N variant of unknown significance [3]. Reminiscent of patients with the S135L mutation, our patient had loss of GALT enzyme activity in erythrocytes but a very mild clinical phenotype [3–8]. We performed splicing experiments and computational structural analyses to investigate the role of the novel S222N variant. Alamut software data predicted loss of splicing enhancers for the S222N and S135L mutations [9,10]. A cDNA library was generated from our patient׳s RNA to investigate for splicing errors, but no change in transcript length was seen [3]. In silico structural analysis was performed to investigate enzyme stability and attempt to understand the mechanism of the atypical galactosemia phenotype. Stability results are publicly available in the GALT Protein Database 2.0 [11–14]. Animations were created to give the reader a dynamic view of the enzyme structure and mutation locations. Protein database files and python scripts are included for further investigation.

Compound heterozygosity with a novel S222N GALT mutation leads to atypical galactosemia with loss of GALT activity in erythrocytes but little evidence of clinical disease

Galactosemia is an inborn error of galactose metabolism caused by mutations in the GALT gene. Though early detection and galactose restriction prevent severe liver disease, affected individuals have persistently elevated biomarkers and often neuro-developmental symptoms. We present a teenage compound heterozygote for a known pathogenic mutation (H132Q) and a novel variant of unknown significance (S222N), with nearly absent erythrocyte GALT enzyme activity but normal biomarkers and only mild anxiety despite diet non-adherence. This case is similar to a previously reported S135L mutation. In this report we investigate the novel S222N variant and critically evaluate a clinically puzzling case.

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.

A frequent splicing mutation and novel missense mutations color the updated mutational spectrum of classic galactosemia in Portugal

Classic galactosemia is an autosomal recessive disorder caused by deficient galactose-1-phosphate uridylyltransferase (GALT) activity. Patients develop symptoms in the neonatal period, which can be ameliorated by dietary restriction of galactose. Many patients develop long-term complications, with a broad range of clinical symptoms whose pathophysiology is poorly understood. The high allelic heterogeneity of GALT gene that characterizes this disorder is thought to play a determinant role in biochemical and clinical phenotypes. We aimed to characterize the mutational spectrum of GALT deficiency in Portugal and to assess potential genotype-phenotype correlations. Direct sequencing of the GALT gene and in silico analyses were employed to evaluate the impact of uncharacterized mutations upon GALT functionality. Molecular characterization of 42 galactosemic Portuguese patients revealed a mutational spectrum comprising 14 nucleotide substitutions: ten missense, two nonsense and two putative splicing mutations. Sixteen different genotypic combinations were detected, half of the patients being p.Q188R homozygotes. Notably, the second most frequent variation is a splicing mutation. In silico predictions complemented by a close-up on the mutations in the protein structure suggest that uncharacterized missense mutations have cumulative point effects on protein stability, oligomeric state, or substrate binding. One splicing mutation is predicted to cause an alternative splicing event. This study reinforces the difficulty in establishing a genotype-phenotype correlation in classic galactosemia, a monogenic disease whose complex pathogenesis and clinical features emphasize the need to expand the knowledge on this "cloudy" disorder.

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.

Homology modeling a fast tool for drug discovery: current perspectives

Major goal of structural biology involve formation of protein-ligand complexes; in which the protein molecules act energetically in the course of binding. Therefore, perceptive of protein-ligand interaction will be very important for structure based drug design. Lack of knowledge of 3D structures has hindered efforts to understand the binding specificities of ligands with protein. With increasing in modeling software and the growing number of known protein structures, homology modeling is rapidly becoming the method of choice for obtaining 3D coordinates of proteins. Homology modeling is a representation of the similarity of environmental residues at topologically corresponding positions in the reference proteins. In the absence of experimental data, model building on the basis of a known 3D structure of a homologous protein is at present the only reliable method to obtain the structural information. Knowledge of the 3D structures of proteins provides invaluable insights into the molecular basis of their functions. The recent advances in homology modeling, particularly in detecting and aligning sequences with template structures, distant homologues, modeling of loops and side chains as well as detecting errors in a model contributed to consistent prediction of protein structure, which was not possible even several years ago. This review focused on the features and a role of homology modeling in predicting protein structure and described current developments in this field with victorious applications at the different stages of the drug design and discovery.

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.

Structural and molecular biology of type I galactosemia: disease-associated mutations

Type I galactosemia results from reduced galactose 1-phosphate uridylyltransferase (GALT) activity. Signs of disease include damage to the eyes, brain, liver, and ovaries. However, the exact nature and severity of the pathology depends on the mutation(s) in the patient's genes and his/her environment. Considerable enzymological and structural knowledge has been accumulated and this provides a basis to explain, at a biochemical level, impairment in the enzyme in the more than 230 disease-associated variants, which have been described. The most common variant, Q188R, occurs close to the active site and the dimer interface. The substitution probably disrupts both UDP-sugar binding and homodimer stability. Other alterations, for example K285N, occur close to the surface of the enzyme and most likely affect the folding and stability of the enzyme. There are a number of unanswered questions in the field, which require resolution. These include the possibility that the main enzymes of galactose metabolism form a supramolecular complex and the need for a high resolution crystal structure of human GALT.

Galactose toxicity in animals

In most organisms, productive utilization of galactose requires the highly conserved Leloir pathway of galactose metabolism. Yet, if this metabolic pathway is perturbed due to congenital deficiencies of the three associated enzymes, or an overwhelming presence of galactose, this monosaccharide which is abundantly present in milk and many non-dairy foodstuffs, will become highly toxic to humans and animals. Despite more than four decades of intense research, little is known about the molecular mechanisms of galactose toxicity in human patients and animal models. In this contemporary review, we take a unique approach to present an overview of galactose toxicity resulting from the three known congenital disorders of galactose metabolism and from experimental hypergalactosemia. Additionally, we update the reader about research progress on animal models, as well as advances in clinical management and therapies of these disorders.

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

We describe the prediction of the structural and functional effects of mutations on the enzyme galactose-1-phosphate uridyltransferase related to the genetic disease galactosemia, using a fully computational approach. One hundred and seven single-point mutants were simulated starting from the structural model of the enzyme obtained by homology modeling methods. Several bioinformatics programs were then applied to each resulting mutant protein to analyze the effect of the mutations. The mutations have a direct effect on the active site, or on the dimer assembly and stability, or on the monomer stability. We describe how mutations may exert their effect at a molecular level by altering H-bonds, salt bridges, secondary structure or surface features. The alteration of protein stability, at level of monomer and/or dimer, is the main effect observed. We found an agreement between our results and the functional experimental data available in literature for some mutants. The data and analyses for all the mutants are fully available in the web-accessible database hosted at http://bioinformatica.isa.cnr.it/GALT.

GALT protein database, a bioinformatics resource for the management and analysis of structural features of a galactosemia-related protein and its mutants

We describe the GALT-Prot database and its related web-based application that have been developed to collect information about the structural and functional effects of mutations on the human enzyme galactose-1-phosphate uridyltransferase (GALT) involved in the genetic disease named galactosemia type I. Besides a list of missense mutations at gene and protein sequence levels, GALT-Prot reports the analysis results of mutant GALT structures. In addition to the structural information about the wild-type enzyme, the database also includes structures of over 100 single point mutants simulated by means of a computational procedure, and the analysis to each mutant was made with several bioinformatics programs in order to investigate the effect of the mutations. The web-based interface allows querying of the database, and several links are also provided in order to guarantee a high integration with other resources already present on the web. Moreover, the architecture of the database and the web application is flexible and can be easily adapted to store data related to other proteins with point mutations. GALT-Prot is freely available at http://bioinformatica.isa.cnr.it/GALT/.

Dominant versus recessive: molecular mechanisms in metabolic disease

Inborn errors of metabolism used to be regarded as simple monogenic traits, but a closer look at how different alleles of a gene determine different phenotypes shows that the molecular mechanisms in the individual case are often complicated. Most metabolic disorders represent a spectrum of phenotypes from normal via attenuated to severe (and sometimes prenatally fatal), and disease manifestation is often influenced by other specific genetic or exogenous factors. The terms 'dominant' or 'recessive' relate to the functional consequences of differing alleles in the (compound) heterozygous individual; the terms are irrelevant for homozygous individuals and inappropriate for X-linked disorders. Mutations affecting the same amino acid residue may be associated with different inheritance patterns. True dominant inheritance in metabolism is rare; it may be found e.g. in tightly regulated biosynthetic pathways or when minor changes in metabolite concentrations have a functional effect. Some disorders such as erythropoietic protoporphyria show pseudodominant inheritance due to prevalent loss-of-function polymorphisms in the general population and are better acknowledged as recessive traits. The term 'variable expressivity' is not helpful with regard to autosomal recessive disorders when variable phenotypes are explained by different mutations in the respective gene. Clonal unmasking of a heterozygous mutation through somatic loss of the second allele, the main pathomechanism in inherited tumour predisposition syndromes, is rare in metabolic disorders, but focal congenital hyperinsulinism is a notable exception. Somatic mosaicism for an OTC gene mutation is given as an example of an apparently heterozygous mutation pattern in a boy with an X-linked disease.

Simulation of conformational changes occurring when a protein interacts with its receptor

In order to simulate the conformational changes occurring when a protein interacts with its receptor, we firstly evaluated the structural differences between the experimental unbound and bound conformations for selected proteins and created theoretical complexes by replacing, in each experimental complex, the protein-bound with the protein-unbound chain. The theoretical models were then subjected to additional modeling refinements to improve the side chain geometry. Comparing the theoretical and experimental complexes in term of structural and energetic factors is resulted that the refined theoretical complexes became more similar to the experimental ones. We applied the same procedure within an homology modeling experiment, using as templates the experimental structures of human interleukin-1beta (IL-1beta) unbound and bound with its receptor, to build models of the homologous proteins from mouse and trout in unbound and bound conformations and to simulate the interaction with the related receptors. Our results suggest that homology modeling techniques are sensitive to differences between bound and unbound conformations, and that modeling with accuracy the side chains in the complex improves the interaction and molecular recognition. Moreover, our refinement procedure could be used in protein-protein interaction studies and, also, applied in conjunction with rigid-body docking when is not available the protein-bound conformation.

Modeling the 3D structure of wheat subtilisin/chymotrypsin inhibitor (WSCI). Probing the reactive site with two susceptible proteinases by time-course analysis and molecular dynamics simulations

Comparative modeling and time-course hydrolysis experiments have been applied to investigate two enzyme-inhibitor complexes formed between the wheat subtilisin-chymotrypsin inhibitor (WSCI) and two susceptible proteinases. WSCI represents the first case of a wheat protein inhibitor active against animal chymotrypsins and bacterial subtilisins. The model was created using as template structure that of the CI-2A inhibitor from barley (PDB code: 2CI2), which shares 87% sequence identity with WSCI. Under these conditions of high similarity, the comparative modeling approach can be successfully applied. We predicted the WSCI 3D model and used it to investigate enzyme-inhibitor complex systems. Experimental observations indicated that chymotrypsin, but not subtilisin, in addition to cleavage at the primary reactive site Met48-Glu49, is able to hydrolyze a second peptide bond between Phe58 and Val59. Here, we report on cleavage of the peptide bond at the inhibitor's reactive site (Met48-Glu49) determined using time-course hydrolysis experiments; the same event was investigated for both subtilisin/WSCI and chymotrypsin/WSCI complexes using molecular dynamics simulations. The molecular details of the initial inhibitor-enzyme interactions, as well as of the changes observed during the simulations, allow us to speculate on the different fates of the two WSCI-proteinase complexes.