Genomic analyses of glycine decarboxylase neurogenic mutations yield a large-scale prediction model for prenatal disease.

Hundreds of mutations in a single gene result in rare diseases, but why mutations induce severe or attenuated states remains poorly understood. Defect in glycine decarboxylase (GLDC) causes Non-ketotic Hyperglycinemia (NKH), a neurological disease associated with elevation of plasma glycine. We unified a human multiparametric NKH mutation scale that separates severe from attenuated neurological disease with new in silico tools for murine and human genome level-analyses, gathered in vivo evidence from mice engineered with top-ranking attenuated and a highly pathogenic mutation, and integrated the data in a model of pre- and post-natal disease outcomes, relevant for over a hundred major and minor neurogenic mutations. Our findings suggest that highly severe neurogenic mutations predict fatal, prenatal disease that can be remedied by metabolic supplementation of dams, without amelioration of persistent plasma glycine. The work also provides a systems approach to identify functional consequences of mutations across hundreds of genetic diseases. Our studies provide a new framework for a large scale understanding of mutation functions and the prediction that severity of a neurogenic mutation is a direct measure of pre-natal disease in neurometabolic NKH mouse models. This framework can be extended to analyses of hundreds of monogenetic rare disorders where the underlying genes are known but understanding of the vast majority of mutations and why and how they cause disease, has yet to be realized.

Obesity increases hepatic glycine dehydrogenase and aminomethyltransferase expression while dietary glycine supplementation reduces white adipose tissue in Zucker diabetic fatty rats.

Obesity is associated with altered glycine metabolism in humans. This study investigated the mechanisms regulating glycine metabolism in obese rats. Eight-week-old Zucker diabetic fatty rats (ZDF; a type-II diabetic animal model) received either 1% glycine or 1.19% L-alanine (isonitrogenous control) in drinking water for 6 weeks. An additional group of lean Zucker rats also received 1.19% L-alanine as a lean control. Glycine concentrations in serum and liver were markedly lower in obese versus lean rats. Enteral glycine supplementation restored both serum and hepatic glycine levels, while reducing mesenteric and internal white fat mass compared with alanine-treated ZDF rats. Blood glucose and non-esterified fatty acid (NEFA) concentrations did not differ between the control and glycine-supplemented ZDF rats (P > 0.10). Both mRNA and protein expression of aminomethyltransferase (AMT) and glycine dehydrogenase, decarboxylating (GLDC) were increased in the livers of obese versus lean rats (P < 0.05). In contrast, glycine cleavage system H (GCSH) hepatic mRNA expression was downregulated in obese versus lean rats, although there was no change in protein expression. These findings indicate that reduced quantities of glycine observed in obese subjects likely results from an upregulation of the hepatic glycine cleavage system and that dietary glycine supplementation potentially reduces obesity in ZDF rats.

Impaired folate 1-carbon metabolism causes formate-preventable hydrocephalus in glycine decarboxylase-deficient mice.

Ventriculomegaly and hydrocephalus are associated with loss of function of glycine decarboxylase (Gldc) in mice and in humans suffering from non-ketotic hyperglycinemia (NKH), a neurometabolic disorder characterized by accumulation of excess glycine. Here, we showed that ventriculomegaly in Gldc-deficient mice is preceded by stenosis of the Sylvian aqueduct and malformation or absence of the subcommissural organ and pineal gland. Gldc functions in the glycine cleavage system, a mitochondrial component of folate metabolism, whose malfunction results in accumulation of glycine and diminished supply of glycine-derived 1-carbon units to the folate cycle. We showed that inadequate 1-carbon supply, as opposed to excess glycine, is the cause of hydrocephalus associated with loss of function of the glycine cleavage system. Maternal supplementation with formate prevented both ventriculomegaly, as assessed at prenatal stages, and postnatal development of hydrocephalus in Gldc-deficient mice. Furthermore, ventriculomegaly was rescued by genetic ablation of 5,10-methylene tetrahydrofolate reductase (Mthfr), which results in retention of 1-carbon groups in the folate cycle at the expense of transfer to the methylation cycle. In conclusion, a defect in folate metabolism can lead to prenatal aqueduct stenosis and resultant hydrocephalus. These defects are preventable by maternal supplementation with formate, which acts as a 1-carbon donor.

Mutations in genes encoding the glycine cleavage system predispose to neural tube defects in mice and humans.

Neural tube defects (NTDs), including spina bifida and anencephaly, are common birth defects of the central nervous system. The complex multigenic causation of human NTDs, together with the large number of possible candidate genes, has hampered efforts to delineate their molecular basis. Function of folate one-carbon metabolism (FOCM) has been implicated as a key determinant of susceptibility to NTDs. The glycine cleavage system (GCS) is a multi-enzyme component of mitochondrial folate metabolism, and GCS-encoding genes therefore represent candidates for involvement in NTDs. To investigate this possibility, we sequenced the coding regions of the GCS genes: AMT, GCSH and GLDC in NTD patients and controls. Two unique non-synonymous changes were identified in the AMT gene that were absent from controls. We also identified a splice acceptor site mutation and five different non-synonymous variants in GLDC, which were found to significantly impair enzymatic activity and represent putative causative mutations. In order to functionally test the requirement for GCS activity in neural tube closure, we generated mice that lack GCS activity, through mutation of AMT. Homozygous Amt(-/-) mice developed NTDs at high frequency. Although these NTDs were not preventable by supplemental folic acid, there was a partial rescue by methionine. Overall, our findings suggest that loss-of-function mutations in GCS genes predispose to NTDs in mice and humans. These data highlight the importance of adequate function of mitochondrial folate metabolism in neural tube closure.

Computational analysis of storage synthesis in developing Brassica napus L. (oilseed rape) embryos: flux variability analysis in relation to ¹³C metabolic flux analysis.

Plant oils are an important renewable resource, and seed oil content is a key agronomical trait that is in part controlled by the metabolic processes within developing seeds. A large-scale model of cellular metabolism in developing embryos of Brassica napus (bna572) was used to predict biomass formation and to analyze metabolic steady states by flux variability analysis under different physiological conditions. Predicted flux patterns are highly correlated with results from prior ¹³C metabolic flux analysis of B. napus developing embryos. Minor differences from the experimental results arose because bna572 always selected only one sugar and one nitrogen source from the available alternatives, and failed to predict the use of the oxidative pentose phosphate pathway. Flux variability, indicative of alternative optimal solutions, revealed alternative pathways that can provide pyruvate and NADPH to plastidic fatty acid synthesis. The nutritional values of different medium substrates were compared based on the overall carbon conversion efficiency (CCE) for the biosynthesis of biomass. Although bna572 has a functional nitrogen assimilation pathway via glutamate synthase, the simulations predict an unexpected role of glycine decarboxylase operating in the direction of NH4⁺ assimilation. Analysis of the light-dependent improvement of carbon economy predicted two metabolic phases. At very low light levels small reductions in CO2 efflux can be attributed to enzymes of the tricarboxylic acid cycle (oxoglutarate dehydrogenase, isocitrate dehydrogenase) and glycine decarboxylase. At higher light levels relevant to the ¹³C flux studies, ribulose-1,5-bisphosphate carboxylase activity is predicted to account fully for the light-dependent changes in carbon balance.

[The isotope effect in the glycine dehydrogenase reaction is the cause of the intramolecular isotope inhomogeneity of glucose carbon of starch synthesized during photorespiration].

The isotope distribution of glucose-6-phosphate in the main pathways of its biosynthesis (in the processes of CO2 assimilation and photorespiration in the Calvin cycle and during resynthesis from the degradation products of lipids and proteins) was analyzed. For reconstructing the isotope distribution of glucoso-6-phosphate synthesized in the Calvin cycle during photorespiration, the functioning of the cycle with regard to its coupling with the glycolate chain, which together constitute the photorespiration chain, was considered. In the glycine dehydrogenase reaction of the glycolate cycle, there arises an isotope effect, which determines the distribution of isotopes in the glucose-6-phosphate and other photorespiration products. The isotope effect of the glycine dehydrogenase reaction increases at the expense of the exhaustion of glucose resources feeding the photorespiration chain. As a result, atoms C-3 and C-4 of glucose become enriched with the heavy isotope, and subsequent mixing of atoms and the specificity of interactions in the photorespiration chain lead to an isotope weighting of the other atoms and an uneven distribution of carbon isotopes in glucose-6-phosphate and other photorespiration products. A comparison of the glucose-6-phosphate isotope patterns in different pathways of the synthesis with the experimental data on the distribution of carbon isotopes in starch glucose of storing plant organs led to the conclusion that the starch resources are predominantly formed at the expense of glucose-6-phosphate of photorespiration. This is consistent with the earlier observed enhancement of photorespiration at the stage of plant maturation.

Interaction between photorespiration and respiration in transgenic potato plants with antisense reduction in glycine decarboxylase.

Potato (Solanum tuberosum L. cv. Désirée) plants with an antisense reduction in the P-protein of the glycine decarboxylase complex (GDC) were used to study the interaction between respiration and photorespiration. Mitochondria isolated from transgenic plants had a decreased capacity for glycine oxidation and glycine accumulated in the leaves. Malate consumption increased in leaves of GDC deficient plants and the capacity for malate and NADH oxidation increased in isolated mitochondria. A lower level of alternative oxidase protein and decreased partitioning of electrons to the alternative pathway was found in these plants. The adenylate status was altered in protoplasts from transgenic plants, most notably the chloroplastic ATP/ADP ratio increased. The lower capacity for photorespiration in leaves of GDC deficient plants was compensated for by increased respiratory decarboxylations in the light. This is interpreted as a decreased light suppression of the tricarboxylic acid cycle in GDC deficient plants in comparison to wild-type plants. The results support the view that respiratory decarboxylations in the light are restricted at the level of the pyruvate dehydrogenase complex and/or isocitrate dehydrogenase and that this effect is likely to be mediated by mitochondrial photorespiratory products.

Identification of a novel one-carbon metabolism regulon in Saccharomyces cerevisiae.

Glycine specifically induces genes encoding subunits of the glycine decarboxylase complex (GCV1, GCV2, and GCV3), and this is mediated by a fall in cytoplasmic levels of 5,10-methylenetetrahydrofolate caused by inhibition of cytoplasmic serine hydroxymethyltransferase. Here it is shown that this control system extends to genes for other enzymes of one-carbon metabolism and de novo purine biosynthesis. Northern analysis of the response to glycine demonstrated that the induction of the GCV genes and the induction of other amino acid metabolism genes are temporally distinct. The genome-wide response to glycine revealed that several other genes are rapidly co-induced with the GCV genes, including SHM2, which encodes cytoplasmic serine hydroxymethyltransferase. These results were refined by examining transcript levels in an shm2Delta strain (in which cytoplasmic 5,10-methylenetetrahydrofolate levels are reduced) and a met13Delta strain, which lacks the main methylenetetrahydrofolate reductase activity of yeast and is effectively blocked at consumption of 5,10-methylene tetrahydrofolate for methionine synthesis. Glycine addition also caused a substantial transient disturbance to metabolism, including a sequence of changes in induction of amino acid biosynthesis and respiratory chain genes. Analysis of the glycine response in the shm2Delta strain demonstrated that apart from the one-carbon regulon, most of these transient responses were not contingent on a disturbance to one-carbon metabolism. The one-carbon response is distinct from the Bas1p purine biosynthesis regulon and thus represents the first example of transcriptional regulation in response to activated one-carbon status.

Structure and expression of the glycine cleavage system in rat central nervous system.

The glycine cleavage system (GCS) is a mitochondrial multienzyme system consisting of four individual proteins, three specific components (P-, T-, and H-proteins) and one house-keeping enzyme, dihydrolipoamide dehydrogenase. Inherited deficiency of the GCS causes nonketotic hyperglycinemia (NKH), an inborn error of glycine metabolism. NKH is characterized by massive accumulation of glycine in serum and cerebrospinal fluids and severe neuronal dysfunction in neonates. To elucidate the neuropathogenesis of NKH, we cloned cDNAs encoding three specific components of the GCS and studied the gene expression in rat central nervous system. P-, T-, and H-protein cDNAs encoded 1024, 403, and 170 amino acids, respectively. In situ hybridization analysis revealed that P-protein mRNA was expressed mainly in glial-like cells, including Bergmann glias in the cerebellum, while T- and H-protein mRNAs were detected in both glial-like cells and neurons. T- and H-protein mRNAs, but not P-protein mRNA, were expressed in the spinal cord. Primary astrocyte cultures established from cerebral cortex had higher GCS activities than hepatocytes whereas those from spinal cord expressed only H-protein mRNA and had no enzymatic activity. An important role of glycine as inhibitory neurotransmitter has been established in the brainstem and spinal cord and another role of glycine as an excitation modulator of N-methyl-D-aspartate receptor is suggested in the hippocampus, cerebral cortex, olfactory bulbus, and cerebellum. Our results suggest that the GCS plays a major role in the forebrain and cerebellum rather than in the spinal cord, and that N-methyl-D-aspartate receptor may participate in neuropathogenesis of NKH.

Metabolic response of potato plants to an antisense reduction of the P-protein of glycine decarboxylase.

Potato (Solanum tuberosum L. cv. Desiré) plants with reduced amounts of P-protein, one of the subunits of glycine decarboxylase (GDC), have been generated by introduction of an antisense transgene. Two transgenic lines, containing about 60-70% less P-protein in the leaves compared to wild-type potato, were analysed in more detail. The reduction in P-protein amount led to a decrease in the ability of leaf mitochondria to decarboxylate glycine. Photosynthetic and growth rates were reduced but the plants were viable under ambient air and produced tubers. Glycine concentrations within the leaves were elevated up to about 100-fold during illumination. Effects on other amino acids and on sucrose and hexoses were minor. Nearly all of the glycine accumulated during the day was metabolised during the following night. The data suggest that the GDC operates far below substrate saturation under normal conditions thus allowing a flexible and fast response to changes in the environment.

The cytotoxic lipid peroxidation product, 4-hydroxy-2-nonenal, specifically inhibits decarboxylating dehydrogenases in the matrix of plant mitochondria.

4-Hydroxy-2-nonenal (HNE), a cytotoxic product of lipid peroxidation, inhibits O(2) consumption by potato tuber mitochondria. 2-Oxoglutarate dehydrogenase (OGDC), pyruvate dehydrogenase complex (PDC) (both 80% inhibited) and NAD-malic enzyme (50% inhibited) are its major targets. Mitochondrial proteins identified by reaction with antibodies raised to lipoic acid lost this antigenicity following HNE treatment. These proteins were identified as acetyltransferases of PDC (78 kDa and 55 kDa), succinyltransferases of OGDC (50 kDa and 48 kDa) and glycine decarboxylase H protein (17 kDa). The significance of the effect of these inhibitions on the impact of lipid peroxidation and plant respiratory functions is discussed.

Regulation of the balance of one-carbon metabolism in Saccharomyces cerevisiae.

One-carbon metabolism in yeast is an essential process that relies on at least one of three one-carbon donor molecules: serine, glycine, or formate. By a combination of genetics and biochemistry we have shown how cells regulate the balance of one-carbon flow between the donors by regulating cytoplasmic serine hydroxymethyltransferase activity in a side reaction occurring in the presence of excess glycine. This control governs the level of 5,10-methylene tetrahydrofolate (5,10-CH(2)-H(4)folate) in the cytoplasm, which has a direct role in signaling transcriptional control of the expression of key genes, particularly those encoding the unique components of the glycine decarboxylase complex (GCV1, GCV2, and GCV3). Based on these and other observations, we propose a model for how cells balance the need to supplement their one-carbon pools when charged folates are limiting or when glycine is in excess. We also propose that under normal conditions, cytoplasmic 5,10-CH(2)-H(4)folate is mainly directed to generating methyl groups via methionine, whereas one-carbon units generated from glycine in mitochondria are more directed to purine biosynthesis. When glycine is in excess, 5, 10-CH(2)-H(4)folate is decreased, and the regulation loop shifts the balance of generation of one-carbon units into the mitochondrion.

Synthesis and characterization of selenolipoylated H-protein of the glycine cleavage system.

H-protein of the glycine cleavage system has a lipoic acid prosthetic group. Selenolipoic acid is a lipoic acid analog in which both sulfur atoms are replaced by selenium atoms. Two isoforms of bovine lipoyltransferase that are responsible for the attachment of lipoic acid to H-protein had an affinity for selenolipoyl-AMP and transferred the selenolipoyl moiety to bovine apoH-protein comparable to lipoyl-AMP. Selenolipoylated H-protein was overexpressed in Escherichia coli and purified. Selenolipoylated H-protein was 26% as effective as lipoylated H-protein in the glycine decarboxylation reaction, in which reduction of the diselenide bond of selenolipoylated H-protein is catalyzed by P-protein. The diselenide form of selenolipoylated H-protein was a poor substrate for L-protein, and the rate of reduction was 0.5% of that of lipoylated H-protein. The rate of the overall glycine cleavage reaction with selenolipoylated H-protein was <1% of that with lipoylated H-protein. These results are consistent with the difference in the redox potential between the diselenide and disulfide bonds. In contrast, selenolipoylated H-protein showed three times as high glycine-14CO2 exchange activity as lipoylated H-protein, presumably because the rate of reoxidation of reduced selenolipoylated H-protein is much higher than that of lipoylated H-protein.

Molecular characterization of GCV3, the Saccharomyces cerevisiae gene coding for the glycine cleavage system hydrogen carrier protein.

YAL044, a gene on the left arm of Saccharomyces cerevisiae chromosome one, is shown to code for the H-protein subunit of the multienzyme glycine cleavage system. The gene designation has therefore been changed to GCV3, reflecting its role in the glycine cleavage system. GCV3 encodes a 177-residue protein with a putative mitochondrial targeting signal at its amino terminus. Targeted gene replacement shows that GCV3 is not required for growth on minimal medium; however, it is essential when glycine serves as the sole nitrogen source. Studies of GCV3 expression revealed that it is highly regulated. Supplementation of minimal medium with glycine, the glycine cleavage system's substrate, induced expression at least 30-fold. In contrast, and consistent with the cleavage of glycine providing activated single-carbon units, the addition of the metabolic end products that require activated single-carbon units repressed expression about 10-fold. Finally, like many amino acid biosynthetic genes, GCV3 is subject to regulation by the general amino acid control system.

Interaction between glycine decarboxylase, serine hydroxymethyltransferase and tetrahydrofolate polyglutamates in pea leaf mitochondria.

The aim of the present work was to further determine how the T-protein of the glycine-cleavage system and serine hydroxy-methyltransferase (SHMT), two folate-dependent enzymes from pea leaf mitochondria, interact through a common pool of tetrahydrofolate polyglutamates (H4PteGlun). It was observed that the binding affinity of tetrahydrofolate polyglutamates for these proteins continuously increased with increasing number of glutamates up to six residues. It was also established that, once bound to the proteins, tetrahydrofolate, a very O2-sensitive molecule, was protected from oxidative degradation. The dissociation constants (Kd) of H4PteGlu5, the most predominant form of polyglutamate in the mitochondria, were approximately 0.5 microM for both T-protein and SHMT, whereas the Kd values of CH2-H4PteGlu5 were higher, 2.7 and 7 microM respectively. In a matrix extract from pea leaf mitochondria, the maximal activity of the glycine-cleavage system was about 2.5 times higher than the maximal activity of SHMT. This resulted in a permanent disequilibrium of the SHMT-catalysed reaction which was therefore driven toward the production of serine and H4PteGlun, the thermodynamically unfavourable direction. Indeed, measurements of the steady-state ratio of CH2-H4PteGlun/H4PteGlun (n = 1 or n = 5) during the course of glycine oxidation demonstrated that the methylene form accounted for 65-80% of the folate pool. This indicates that, in our in vitro experiments, CH2-H4PteGlun with long polyglutamate chains accumulated in the bulk medium. This observation suggests that, in these in vitro experiments at least, there was no channelling of CH2-H4PteGlu5 between the T-protein and SHMT.

X-ray structure determination at 2.6-A resolution of a lipoate-containing protein: the H-protein of the glycine decarboxylase complex from pea leaves.

H-protein, a lipoic acid-containing protein of the glycine decarboxylase (EC 1.4.4.2) complex from pea (Pisum sativum) was crystallized from ammonium sulfate solution at pH 5.2 in space group P3(1)21. The x-ray crystal structure was determined to 2.6-A resolution by multiple isomorphous replacement techniques. The structure was refined to an R value of 23% for reflections between 15- and 2.6-A resolution (F > 2 sigma), including the lipoate moiety and 50 water molecules, for the two protein molecules of the asymmetric unit. The 131-amino acid residues form seven beta-strands arranged into two antiparallel beta-sheets forming a "sandwich" structure. One alpha-helix is observed at the C-terminal end. The lipoate cofactor attached to Lys-63 is located in the loop of a hairpin configuration. The lipoate moiety points toward the residues His-34 and Asp-128 and is situated at the surface of the H-protein. This allows the flexibility of the lipoate arm. This is the first x-ray determination of a lipoic acid-containing protein, and the present results are in agreement with previous theoretical predictions and NMR studies of the catalytic domains of lipoic acid- and biotin-containing proteins.

Characterization of the primary structure of H-protein from Pisum sativum and location of a lipoic acid residue by combined liquid chromatography/mass spectrometry and liquid chromatography/tandem mass spectrometry.

A purified extract of H-protein, a subunit of the glycine cleavage complex of the pea leaf mitochondria, was investigated by liquid chromatography/mass spectrometry (LC/MS) and liquid chromatography/tandem mass spectrometry (LC/MS/MS), using both continuous flow fast atom bombardment (CF-FAB) and electrospray ionization (ESI) mass spectrometry. Determination of the molecular weight of the entire protein, a 14 kDa subunit of the glycine decarboxylase complex, was achieved by ESI mass spectrometry and revealed covalent binding of the protein to the stabilizing agent beta-mercapto-ethanol. On-line LC/MS analysis of peptides arising from the endoproteinase Glu-C digestion of the H-protein was achieved using capillary columns (0.25 mm i.d.), and permitted confirmation of the previously reported sequence deduced from cDNA cloning experiments. The detailed interpretation of data extracted from these LC/MS experiments facilitated identification of peptides containing modified amino acid residues. In particular the identification of a lipoic acid cofactor, a rather unusual modified lysine residue which interacts with different active sites in the enzyme complex, was achieved using both LC/CF-FAB-MS and LC/ESI-MS. The exact location of this modified lysine residue was determined by obtaining fragment spectra of multiply protonated precursor ions of selected peptides, using on-line LC/MS/MS techniques.

Effects of tetrahydrofolate polyglutamates on the kinetic parameters of serine hydroxymethyltransferase and glycine decarboxylase from pea leaf mitochondria.

Plant tissues contain highly conjugated forms of folate. Despite this, the ability of plant folate-dependent enzymes to utilize tetrahydrofolate polyglutamates has not been examined in detail. In leaf mitochondria, the glycine-cleavage system and serine hydroxymethyltransferase, present in large amounts in the matrix space and involved in the photorespiratory cycle, necessitate the presence of tetrahydrofolate as a cofactor. The aim of the present work was to determine whether glutamate chain length (one to six glutamate residues) influenced the affinity constant for tetrahydrofolate and the maximal velocities displayed by these two enzymes. The results show that the affinity constant decreased by at least one order of magnitude when the tetrahydrofolate substrate contained three or more glutamate residues. In contrast, maximal velocities were not altered in the presence of these substrates. These results are consistent with analyses of mitochondrial folates which revealed a pool of polyglutamates dominated by tetra and pentaglutamates. The equilibrium constant of the serine hydroxymethyltransferase suggests that, during photorespiration, the reaction must be permanently pushed toward the formation of serine (the unfavourable direction) to allow the recycling of tetrahydrofolate necessary for the operation of the glycine decarboxylase T-protein.

Prediction of the three-dimensional structures of the biotinylated domain from yeast pyruvate carboxylase and of the lipoylated H-protein from the pea leaf glycine cleavage system: a new automated method for the prediction of protein tertiary structure.

A new, automated, knowledge-based method for the construction of three-dimensional models of proteins is described. Geometric restraints on target structures are calculated from a consideration of homologous template structures and the wider knowledge base of unrelated protein structures. Three-dimensional structures are calculated from initial partly folded states by high-temperature molecular dynamics simulations followed slow cooling of the system (simulated annealing) using nonphysical potentials. Three-dimensional models for the biotinylated domain from the pyruvate carboxylase of yeast and the lipoylated H-protein from the glycine cleavage system of pea leaf were constructed, based on the known structures of two lipoylated domains of 2-oxo acid dehydrogenase multienzyme complexes. Despite their weak sequence similarity, the three proteins are predicted to have similar three-dimensional structures, representative of a new protein module. Implications for the mechanisms of posttranslational modification of these proteins and their catalytic function are discussed.