Simulation of the purine nucleotide cycle as an anaplerotic process in skeletal muscle

A computer model of purine nucleotide and citric acid cycles joined through fumarate is given. Steady-state equations corresponding to metabolic enzymes are written based on the information from the literature about their kinetic behavior. Numerical integration of this set of equations is performed and in order to maintain an overall stabilization between the two cycles, enzymatic activities, in the form of V, have been calculated. Sensitivity coefficients for enzymes indicate that the control is exerted, depending upon the intermediate concentrations, and furthermore, it is demonstrated that AMP concentration in muscle should be very low. From stabilization, simulation of exercise conditions has been performed by diminishing [ATP] and increasing accordingly [ADP] and [AMP]. In such conditions the operation of purine nucleotide cycle leads to a considerable increase in the level of citric acid cycle intermediates. Disruption of purine nucleotide cycle by altering some of the three enzymatic steps leads to a lesser increase of these intermediates. The set of results presented seems to confirm the hypothesis that purine nucleotide cycle acts as an anaplerotic process in muscle, as the experimental results of Aragon and Lowenstein (Aragon, J.J., and Lowenstein, J.M. (1980) Eur. J. Biochem. 110, 371-377) suggest.     

Alteration of purine metabolism by AICA-riboside in human B lymphoblasts.

The effect of 5-amino-4-imidazole-carboximide (AI-CA)-riboside on different pathways of purine metabolism (biosynthesis de novo, salvage pathways, adenosine metabolism, ATP catabolism) was studied in human B lymphoblasts (WI-L2). AICA-Riboside markedly decreased intracellular levels of 5-phosphoribosyl-1-pyrophosphate and in consequence affected purine biosynthesis de novo and purine salvage pathways. AICA-riboside inhibited incorporation of glycine into purine nucleotides, but when formate was used as the precursor of purine biosynthesis de novo, a biphasic effect was observed. The incorporation of formate into purine nucleotides was increased by AICA-riboside at concentrations up to 2 mM but decreased at higher concentrations. Salvage of the purine bases adenine, hypoxanthine, and guanine was markedly inhibited and utilization of extracellular adenosine in B lymphoblasts was reduced by AICA-riboside. AICA-riboside increased ribose 1-phosphate concentrations and increased degradation of prelabeled ATP. No effect on the intracellular levels of orthophosphate was found. Proliferation of WI-L2 lymphoblasts was only slightly affected at concentrations of AICA-riboside below 500 microM but markedly inhibited by higher concentrations.     

The purine metabolism of human erythrocytes

This review summarizes currently available information about a crucial part of erythrocyte metabolism, that is, purine nucleotide conversions and their relationships with other conversion pathways. We describe the cellular resynthesis, interconversion, and degradation of purine compounds, and also the regulatory mechanisms in the conversion pathways. We also mention purine metabolism disorders and their clinical consequences. The literature is fragmentary because studies have concentrated only on selected aspects of purine metabolism; hence the need for a synthetic approach. Published in Russian in Biokhimiya, 2006, Vol. 71, No. 5, pp. 581–591.     

Mitochondrial one-carbon metabolism is adapted to the specific needs of yeast, plants and mammals

In eukaryotes, folate metabolism is compartmentalized between the cytoplasm and organelles. The folate pathways of mitochondria are adapted to serve the metabolism of the organism. In yeast, mitochondria support cytoplasmic purine synthesis through the generation of formate. This pathway is important but not essential for survival, consistent with the flexibility of yeast metabolism. In plants, the mitochondrial pathways support photorespiration by generating serine from glycine. This pathway is essential under photosynthetic conditions and the enzyme expression varies with photosynthetic activity. In mammals, the expression of the mitochondrial enzymes varies in tissues and during development. In embryos, mitochondria supply formate and glycine for purine synthesis, a process essential for survival; in adult tissues, flux through mitochondria can favor serine production. The differences in the folate pathways of mitochondria depending on species, tissues and developmental stages, profoundly alter the nature of their metabolic contribution.     

Simultaneous determination of purine and pyrimidine metabolites in HPRT-deficient cell lines

Genetic mutations in the purine salvage enzyme, hypoxanthine-guanine phosphoribosyltransferase (HPRT), are known to cause Lesch-Nyhan syndrome and Kelley-Seegmiller syndrome. In patients, purine metabolism is different from that of normal persons. We have previously developed a method for simultaneously determining the concentration of purine and pyrimidine nucleosides and nucleotides. This system was applied to determine the concentrations of nucleosides and nucleotides in HPRT-deficient cell lines. The amount of inosine 5'-monophosphate (IMP) was different in Lesch-Nyhan syndrome, Kelley-Seegmiller syndrome, and control cell lines. The difference in the amount of IMP confirmed the mutation of the enzyme.     

Effect of enzyme organization on the stability of Yates-Pardee pathways

Yates-Pardee-type metabolic pathways in a heterogenous cell milieu are modeled as a system of coupled non-linear partial differential equations. A numerical solution to this systmm is described and some properties of such a physiological system are studied. Confinement with and without a membrane is considered and it is shown how confinement results in an increase in the stability of the metabolite concentrations. These results suggest that the enzyme organization may contribute to the stability of the cellular metabolism.     

Untargeted plasma metabolite profiling reveals the broad systemic consequences of xanthine oxidoreductase inactivation in mice

A major challenge in systems biology is integration of molecular findings for individual enzyme activities into a cohesive high-level understanding of cellular metabolism and physiology/pathophysiology. However, meaningful prediction for how a perturbed enzyme activity will globally impact metabolism in a cell, tissue or intact organisms is precluded by multiple unknowns, including in vivo enzymatic rates, subcellular distribution and pathway interactions. To address this challenge, metabolomics offers the potential to simultaneously survey changes in thousands of structurally diverse metabolites within complex biological matrices. The present study assessed the capability of untargeted plasma metabolite profiling to discover systemic changes arising from inactivation of xanthine oxidoreductase (XOR), an enzyme that catalyzes the final steps in purine degradation. Using LC-MS coupled with a multivariate statistical data analysis platform, we confidently surveyed >3,700 plasma metabolites (50-1,000 Da) for differential expression in XOR wildtype vs. mice with inactivated XOR, arising from gene deletion or pharmacological inhibition. Results confirmed the predicted derangements in purine metabolism, but also revealed unanticipated perturbations in metabolism of pyrimidines, nicotinamides, tryptophan, phospholipids, Krebs and urea cycles, and revealed kidney dysfunction biomarkers. Histochemical studies confirmed and characterized kidney failure in xor-nullizygous mice. These findings provide new insight into XOR functions and demonstrate the power of untargeted metabolite profiling for systemic discovery of direct and indirect consequences of gene mutations and drug treatments.     

Supercomputer analysis of purine and pyrimidine metabolism leading to DNA synthesis

A model-system is established to analyze purine and pyrimidine metabolism leading to DNA synthesis. The principal aim is to explore the flow and regulation of terminal deoxynucleoside triophosphates (dNTPs) in various input and parametric conditions. A series of flow equations are established, which are subsequently converted to differential equations. These are programmed (Fortran) and analyzed on a Cray X-MP/48 supercomputer. The pool concentrations are presented as a function of time in conditions in which various pertinent parameters of the system are modified. The system is formulated by 100 differential equations.     

Simultaneous Determination of Purine and Pyrimidine Metabolites in Hprt-Deficient Cell Lines

Genetic mutations in the purine salvage enzyme, hypoxanthine-guanine phosphoribosyltransferase (HPRT), are known to cause Lesch–Nyhan syndrome and Kelley–Seegmiller syndrome. In patients, purine metabolism is different from that of normal persons. We have previously developed a method for simultaneously determining the concentration of purine and pyrimidine nucleosides and nucleotides. This system was applied to determine the concentrations of nucleosides and nucleotides in HPRT-deficient cell lines. The amount of inosine 5′-monophosphate (IMP) was different in Lesch–Nyhan syndrome, Kelley–Seegmiller syndrome, and control cell lines. The difference in the amount of IMP confirmed the mutation of the enzyme.     

A study of deoxyribonucleotide metabolism and its relation to DNA synthesis. Supercomputer simulation and model-system analysis.

A model system (1) was established to analyze purine and pyrimidine metabolism. This system has been expanded to include macrosimulation of DNA synthesis and the study of its regulation by terminal deoxynucleoside triphosphates (dNTPs) via a complex set of interactions. Computer experiments reveal that our model exhibits adequate and reasonable sensitivity in terms of dNTP pool levels and rates of DNA synthesis when inputs to the system are varied. These simulation experiments reveal that in order to achieve maximum DNA synthesis (in terms of purine metabolism), a proper balance is required in guanine and adenine input into this metabolic system. Excessive inputs will become inhibitory to DNA synthesis. In addition, studies are carried out on rates of DNA synthesis when various parameters are changed quantitatively. The current system is formulated by 110 differential equations.     

Supercomputer analysis of purine and pyrimidine metabolism leading to DNA synthesis

A model-system is established to analyze purine and pyrimidine metabolism leading to DNA synthesis. The principal aim is to explore the flow and regulation of terminal deoxynucleoside triophosphates (dNTPs) in various input and parametric conditions. A series of flow equations are established, which are subsequently converted to differential equations. These are programmed (Fortran) and analyzed on a Cray chi-MP/48 supercomputer. The pool concentrations are presented as a function of time in conditions in which various pertinent parameters of the system are modified. The system is formulated by 100 differential equations.     

The origin and evolution of modern metabolism

One fundamental goal of current research is to understand how complex biomolecular networks took the form that we observe today. Cellular metabolism is probably one of the most ancient biological networks and constitutes a good model system for the study of network evolution. While many evolutionary models have been proposed, a substantial body of work suggests metabolic pathways evolve fundamentally by recruitment, in which enzymes are drawn from close or distant regions of the network to perform novel chemistries or use different substrates. Here we review how structural and functional genomics has impacted our knowledge of evolution of modern metabolism and describe some approaches that merge evolutionary and structural genomics with advances in bioinformatics. These include mining the data on structure and function of enzymes for salient patterns of enzyme recruitment. Initial studies suggest modern metabolism originated in enzymes of nucleotide metabolism harboring the P-loop hydrolase fold, probably in pathways linked to the purine metabolic subnetwork. This gateway of recruitment gave rise to pathways related to the synthesis of nucleotides and cofactors for an ancient RNA world. Once the TIM beta/alpha-barrel fold architecture was discovered, it appears metabolic activities were recruited explosively giving rise to subnetworks related to carbohydrate and then amino acid metabolism. Remarkably, recruitment occurred in a layered system reminiscent of Morowitz's prebiotic shells, supporting the notion that modern metabolism represents a palimpsest of ancient metabolic chemistries.     

Sequential Activation of Metabolic Pathways: a Dynamic Optimization Approach

The regulation of cellular metabolism facilitates robust cellular operation in the face of changing external conditions. The cellular response to this varying environment may include the activation or inactivation of appropriate metabolic pathways. Experimental and numerical observations of sequential timing in pathway activation have been reported in the literature. It has been argued that such patterns can be rationalized by means of an underlying optimal metabolic design. In this paper we pose a dynamic optimization problem that accounts for time-resource minimization in pathway activation under constrained total enzyme abundance. The optimized variables are time-dependent enzyme concentrations that drive the pathway to a steady state characterized by a prescribed metabolic flux. The problem formulation addresses unbranched pathways with irreversible kinetics. Neither specific reaction kinetics nor fixed pathway length are assumed. In the optimal solution, each enzyme follows a switching profile between zero and maximum concentration, following a temporal sequence that matches the pathway topology. This result provides an analytic justification of the sequential activation previously described in the literature. In contrast with the existent numerical approaches, the activation sequence is proven to be optimal for a generic class of monomolecular kinetics. This class includes, but is not limited to, Mass Action, Michaelis–Menten, Hill, and some Power-law models. This suggests that sequential enzyme expression may be a common feature of metabolic regulation, as it is a robust property of optimal pathway activation.     

5-Amino-4-imidazolecarboxamide riboside (Z-riboside) metabolism in eukaryotic cells

Metabolites of 5-amino-4-imidazolecarboxamide riboside (Z-riboside) have potential roles in the regulation of cellular metabolism and as pharmacological agents in several pathological situations. Before studying Z-riboside metabolism it was necessary to develop methods for identifying and quantitating 5(4)-amino-4(5)-imidazolecarboxamide metabolites. These studies utilized Chinese hamster ovary fibroblast auxotrophic mutants to identify and isolate compounds relevant to Z-riboside metabolism by a combination of high performance liquid chromatographic procedures. In order to study Z-riboside metabolism wild-type and mutant cells were cultured in Z-riboside. This ribosyl precursor to a purine de novo intermediate does not undergo any detectable phosphorolysis but rather is phosphorylated by adenosine kinase in an unregulated manner. This results in the intracellular accumulation of 5-amino-4-imidazolecarboxamide ribotide (ZMP), the levels of which control flow from Z-riboside to the following metabolites: 1) IMP and other purine nucleotides, 2) 5-amino-4-imidazole-N-succinocarboxamide ribotide (sZMP), and 3) 5-amino-4-imidazolecarboxamide riboside 5'-triphosphate (ZTP). At low ZMP concentrations, the predominant metabolic fate is IMP. Initially, IMP enters the adenylate and guanylate pools, but subsequently is hydrolyzed to inosine and this phosphorolyzed to hypoxanthine. At intermediate ZMP concentrations there is net retrograde flux through the bifunctional enzyme adenylosuccinate AMP lyase resulting in sZMP synthesis and antegrade flux leads to the accumulation of adenylosuccinate. At high ZMP concentrations, ZTP accumulates. In addition to these effects on purine metabolism, pyrimidine nucleotide pools are depleted when ZMP accumulates. These results are discussed in relation to the regulation of purine nucleotide synthesis and the use of Z-riboside as a pharmacological intervention in pathophysiological situations.     

Erythrocytic adenosine monophosphate as an alternative purine source in Plasmodium falciparum

Plasmodium falciparum is a purine auxotroph, salvaging purines from erythrocytes for synthesis of RNA and DNA. Hypoxanthine is the key precursor for purine metabolism in Plasmodium. Inhibition of hypoxanthine-forming reactions in both erythrocytes and parasites is lethal to cultured P. falciparum. We observed that high concentrations of adenosine can rescue cultured parasites from purine nucleoside phosphorylase and adenosine deaminase blockade but not when erythrocyte adenosine kinase is also inhibited. P. falciparum lacks adenosine kinase but can salvage AMP synthesized in the erythrocyte cytoplasm to provide purines when both human and Plasmodium purine nucleoside phosphorylases and adenosine deaminases are inhibited. Transport studies in Xenopus laevis oocytes expressing the P. falciparum nucleoside transporter PfNT1 established that this transporter does not transport AMP. These metabolic patterns establish the existence of a novel nucleoside monophosphate transport pathway in P. falciparum.     

Phylogenetic Analysis and Molecular Evolution of Guanine Deaminases: From Guanine to Dendrites

Guanine deaminase (GDA; guanase) is a ubiquitous enzyme that catalyzes the first step of purine metabolism by hydrolytic deamination of guanine, resulting in the production of xanthine. This hydrolase subfamily member plays an essential role in maintaining homeostasis of cellular triphosphate nucleotides for energy, signal transduction pathways, and nitrogen sources. In mammals, GDA protein levels can play a role in neuronal development by regulating dendritic arborization. We previously demonstrated that the most abundant alternative splice form of GDA in mammals, termed cypin (cytosolic PSD-95 interactor), interacts with postsynaptic density proteins, regulates microtubule polymerization, and increases dendrite number. Since purine metabolism and dendrite development were previously thought to be independent cellular processes, this multifunctional protein serves as a new target for the treatment of cognitive disorders characterized by aberrant neuronal morphology and purine metabolism. Although the enzymatic activity of GDA has been conserved during evolution from prokaryotes to higher eukaryotes, a detailed evolutionary assessment of the principal domains in GDA proteins has not yet been put forward. In this study, we perform a complete evolutionary analysis of the full-length sequences and the principal domains in guanine deaminases. Furthermore, we reconstruct the molecular phylogeny of guanine deaminases with neighbor-joining, maximum-likelihood, and UPGMA methods of phylogenetic inference. This study can act as a model whereby a universal housekeeping enzyme may be adapted to act also as a key regulator of a developmental process.     

Characterization of purine nucleotide metabolism in primary rat cardiomyocyte cultures

Primary rat cardiomyocyte cultures were utilized as a model for the study of purine nucleotide metabolism in the heart muscle, especially in connection with the mechanisms operating for the conservation of adenine nucleotides. The cultures exhibited capacity to produce purine nucleotides from nonpurine molecules (de novo synthesis), as well as from preformed purines (salvage synthesis). The conversion of adenosine to AMP, catalyzed by adenosine kinase, appears to be the most important physiological salvage pathway of adenine nucleotide synthesis in the cardiomyocytes. The study of the metabolic fate of IMP formed from [14C]formate or [14C]hypoxanthine and that of AMP formed from [14C]adenine or [14C]adenosine revealed that in the cardiomyocyte the main flow in the nucleotide interconversion pathways is from IMP to AMP, whereas the flux from AMP to IMP appeared to be markedly slower. Following synthesis from labeled precursors by either de novo or salvage pathways, most of the radioactivity in purine nucleotides accumulated in adenine nucleotides, and only a small proportion of it resided in IMP. The results suggest that the main pathway of AMP degradation in the cardiomyocyte proceeds through adenosine rather than through IMP. About 90% of the total radioactivity in purines effluxed from the cells during de novo synthesis from [14C]formate or following prelabeling of adenine nucleotides with [14C]adenine were found to reside in hypoxanthine. The activities in cell extracts of AMP 5'-nucleotidase and IMP 5'-nucleotidase, which catalyze nucleotide degradation, and of AMP deaminase, a key enzyme in the purine nucleotide cycle, were low. The nucleotidase activity resembles, and that of the AMP deaminase contrasts the respective enzyme activities in extracts of cultured skeletal-muscle myotubes. The results indicate that in the cardiomyocyte, in contrast to the myotube, the main mechanism operating for conservation of nucleotides is prompt phosphorylation of AMP, rather than operation of the purine nucleotide cycle. The primary cardiomyocyte cultures are a plausible model for the study of purine nucleotide metabolism in the heart muscle.     

Genome-based kinetic modeling of cytosolic glucose metabolism in industrially relevant cell lines: Saccharomyces cerevisiae and Chinese hamster ovary cells

Model-based analysis of cellular metabolism can facilitate our understanding of intracellular kinetics and aid the improvement of cell growth and biological product manufacturing. In this paper, a model-based kinetic study of cytosolic glucose metabolism for two industrially relevant cell lines, Saccharomyces cerevisiae and Chinese hamster ovary (CHO) cells, based on enzyme genetic presence and expression information is described. We have reconstructed the cytosolic glucose metabolism map for S. cerevisiae and CHO cells, containing 18 metabolites and 18 enzymes using information from the Kyoto Encyclopedia of Genes and Genomes (KEGG). Based on this map, we have developed akinetic mathematical model for the pathways involved,considering regulation and/or inhibition by products orco-substrates. The values of the maximum rates of reactions(V(max)) were estimated based on kinetic parameter information and metabolic flux analysis results available in literature, and the resulting simulation results for steady-state metabolite concentrations are in good agreement with published experimental data. Finally, the model was used to analyse how the production of DHAP, an important intermediate in fine chemicals synthesis, could be increased using gene knockout.