Simulating Metabolic Processes Using an Architecture Based on Networks of Bio-inspired Processors.

In this work, we propose the Networks of Evolutionary Processors (NEP) [2] as a computational model to solve problems related with biological phenomena. In our first approximation, we simulate biological processes related with cellular signaling and their implications in the metabolism, by using an architecture based on NEP (NEP architecture) and their specializations: Networks of Polarized Evolutionary Processors (NPEP) [1] and NEP Transducers (NEPT) [3]. In particular, we use this architecture to simulate the interplay between cellular processes related with the metabolism as the Krebs cycle and the malate-aspartate shuttle pathway (MAS) both being altered by signaling by calcium.

Identification of a human nuclear receptor defines a new signaling pathway for CYP3A induction

Nuclear receptors regulate metabolic pathways in response to changes in the environment by appropriate alterations in gene expression of key metabolic enzymes. Here, a computational search approach based on iteratively built hidden Markov models of nuclear receptors was used to identify a human nuclear receptor, termed hPAR, that is expressed in liver and intestines. hPAR was found to be efficiently activated by pregnanes and by clinically used drugs including rifampicin, an antibiotic known to selectively induce human but not murine CYP3A expression. The CYP3A drug-metabolizing enzymes are expressed in gut and liver in response to environmental chemicals and clinically used drugs. Interestingly, hPAR is not activated by pregnenolone 16α-carbonitrile, which is a potent inducer of murine CYP3A genes and an activator of the mouse receptor PXR.1. Furthermore, hPAR was found to bind to and trans-activate through a conserved regulatory sequence present in human but not murine CYP3A genes. These results provide evidence that hPAR and PXR.1 may represent orthologous genes from different species that have evolved to regulate overlapping target genes in response to pharmacologically distinct CYP3A activators, and have potential implications for the in vitro identification of drug interactions important to humans.

Deficient DNA-ligase activity in the metabolic disease tyrosinemia type I

Hereditary tyrosinemia type I (HT1) is an autosomal recessive inborn error of metabolism caused by the deficiency of fumarylacetoacetate hydrolase, the last enzyme in the tyrosine catabolism pathway. This defect results in accumulation of succinylacetone (SA) that reacts with amino acids and proteins to form stable adducts via Schiff base formation, lysine being the most reactive amino acid. HT1 patients surviving beyond infancy are at considerable risk for the development of hepatocellular carcinoma, and a high level of chromosomal breakage is observed in HT1 cells, suggesting a defect in the processing of DNA. In this paper we show that the overall DNA-ligase activity is low in HT1 cells (about 20% of the normal value) and that Okazaki fragments are rejoined at a reduced rate compared with normal fibroblasts. No mutation was found by sequencing the ligase I cDNA from HT1 cells, and the level of expression of the ligase I mRNA was similar in normal and HT1 fibroblasts, suggesting the presence of a ligase inhibitor. SA was shown to inhibit in vitro the overall DNA-ligase activity present in normal cell extracts. The activity of purified T4 DNA-ligase, whose active site is also a lysine residue, was inhibited by SA in a dose-dependent manner. These results suggest that accumulation of SA reduces the overall ligase activity in HT1 cells and indicate that metabolism errors may play a role in regulating enzymatic activities involved in DNA replication and repair.

Interactions among enzymes of the Arabidopsis flavonoid biosynthetic pathway

Flavonoids are secondary metabolites derived from phenylalanine and acetate metabolism that perform a variety of essential functions in higher plants. Studies over the past 30 years have supported a model in which flavonoid metabolism is catalyzed by an enzyme complex localized to the endoplasmic reticulum [Hrazdina, G. & Wagner, G. J. (1985) Arch. Biochem. Biophys. 237, 88–100]. To test this model further we assayed for direct interactions between several key flavonoid biosynthetic enzymes in developing Arabidopsis seedlings. Two-hybrid assays indicated that chalcone synthase, chalcone isomerase (CHI), and dihydroflavonol 4-reductase interact in an orientation-dependent manner. Affinity chromatography and immunoprecipitation assays further demonstrated interactions between chalcone synthase, CHI, and flavonol 3-hydroxylase in lysates from Arabidopsis seedlings. These results support the hypothesis that the flavonoid enzymes assemble as a macromolecular complex with contacts between multiple proteins. Evidence was also found for posttranslational modification of CHI. The importance of understanding the subcellular organization of elaborate enzyme systems is discussed in the context of metabolic engineering.

Ste6p Mutants Defective in Exit from the Endoplasmic Reticulum (ER) Reveal Aspects of an ER Quality Control Pathway in Saccharomyces cerevisiae

We are studying the intracellular trafficking of the multispanning membrane protein Ste6p, the a-factor transporter in Saccharomyces cerevisiae and a member of the ATP-binding cassette superfamily of proteins. In the present study, we have used Ste6p as model for studying the process of endoplasmic reticulum (ER) quality control, about which relatively little is known in yeast. We have identified three mutant forms of Ste6p that are aberrantly ER retained, as determined by immunofluorescence and subcellular fractionation. By pulse-chase metabolic labeling, we demonstrate that these mutants define two distinct classes. The single member of Class I, Ste6–166p, is highly unstable. We show that its degradation involves the ubiquitin–proteasome system, as indicated by its in vivo stabilization in certain ubiquitin–proteasome mutants or when cells are treated with the proteasome inhibitor drug MG132. The two Class II mutant proteins, Ste6–13p and Ste6–90p, are hyperstable relative to wild-type Ste6p and accumulate in the ER membrane. This represents the first report of a single protein in yeast for which distinct mutant forms can be channeled to different outcomes by the ER quality control system. We propose that these two classes of ER-retained Ste6p mutants may define distinct checkpoint steps in a linear pathway of ER quality control in yeast. In addition, a screen for high-copy suppressors of the mating defect of one of the ER-retained ste6 mutants has identified a proteasome subunit, Hrd2p/p97, previously implicated in the regulated degradation of wild-type hydroxymethylglutaryl-CoA reductase in the ER membrane.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.

Integrated control of appetite and fat metabolism by the leptin-proopiomelanocortin pathway

Leptin deficiency results in a complex obesity phenotype comprising both hyperphagia and lowered metabolism. The hyperphagia results, at least in part, from the absence of induction by leptin of melanocyte stimulating hormone (MSH) secretion in the hypothalamus; the MSH normally then binds to melanocortin-4 receptor expressing neurons and inhibits food intake. The basis for the reduced metabolic rate has been unknown. Here we show that leptin administered to leptin-deficient (ob/ob) mice results in a large increase in peripheral MSH levels; further, peripheral administration of an MSH analogue results in a reversal of their abnormally low metabolic rate, in an acceleration of weight loss during a fast, in partial restoration of thermoregulation in a cold challenge, and in inducing serum free fatty acid levels. These results support an important peripheral role for MSH in the integration of metabolism with appetite in response to perceived fat stores indicated by leptin levels.

Induction of male sterility in plants by metabolic engineering of the carbohydrate supply

Extracellular invertase mediates phloem unloading via an apoplastic pathway. The gene encoding isoenzyme Nin88 from tobacco was cloned and shown to be characterized by a specific spatial and temporal expression pattern. Tissue-specific antisense repression of Nin88 under control of the corresponding promoter in tobacco results in a block during early stages of pollen development, thus, causing male sterility. This result demonstrates a critical role of extracellular invertase in pollen development and strongly supports the essential function of extracellular sucrose cleavage for supplying carbohydrates to sink tissues via the apoplast. The specific interference with phloem unloading, the sugar status, and metabolic signaling during pollen formation will be a potentially valuable approach to induce male sterility in various crop species for hybrid seed production.

Quantitative trait loci and metabolic pathways

The interpretation of quantitative trait locus (QTL) studies is limited by the lack of information on metabolic pathways leading to most economic traits. Inferences about the roles of the underlying genes with a pathway or the nature of their interaction with other loci are generally not possible. An exception is resistance to the corn earworm Helicoverpa zea (Boddie) in maize (Zea mays L.) because of maysin, a C-glycosyl flavone synthesized in silks via a branch of the well characterized flavonoid pathway. Our results using flavone synthesis as a model QTL system indicate: (i) the importance of regulatory loci as QTLs, (ii) the importance of interconnecting biochemical pathways on product levels, (iii) evidence for “channeling” of intermediates, allowing independent synthesis of related compounds, (iv) the utility of QTL analysis in clarifying the role of specific genes in a biochemical pathway, and (v) identification of a previously unknown locus on chromosome 9S affecting flavone level. A greater understanding of the genetic basis of maysin synthesis and associated corn earworm resistance should lead to improved breeding strategies. More broadly, the insights gained in relating a defined genetic and biochemical pathway affecting a quantitative trait should enhance interpretation of the biological basis of variation for other quantitative traits.

The metabolic implications of intracellular circulation

Two views currently dominate research into cell function and regulation. Model I assumes that cell behavior is quite similar to that expected for a watery bag of enzymes and ligands. Model II assumes that three-dimensional order and structure constrain and determine metabolite behavior. A major problem in cell metabolism is determining why essentially all metabolite concentrations are remarkably stable (are homeostatic) over large changes in pathway fluxes—for convenience, this is termed the [s] stability paradox. For muscle cells, ATP and O2 are the most perfectly homeostatic, even though O2 delivery and metabolic rate correlate in a 1:1 fashion. In total, more than 60 metabolites are known to be remarkably homeostatic in differing metabolic states. Several explanations of [s] stability are usually given by traditional model I studies—none of which apply to all enzymes in a pathway, and all of which require diffusion as the means for changing enzyme–substrate encounter rates. In contrast, recent developments in our understanding of intracellular myosin, kinesin, and dyenin motors running on actin and tubulin tracks or cables supply a mechanistic basis for regulated intracellular circulation systems with cytoplasmic streaming rates varying over an approximately 80-fold range (from 1 to >80 μm × sec−1). These new studies raise a model II hypothesis of intracellular perfusion or convection as a primary means for bringing enzymes and substrates together under variable metabolic conditions. In this view, change in intracellular perfusion rates cause change in enzyme–substrate encounter rates and thus change in pathway fluxes with no requirement for large simultaneous changes in substrate concentrations. The ease with which this hypothesis explains the [s] stability paradox is one of its most compelling features.

3-Methylcrotonyl-Coenzyme A Carboxylase Is a Component of the Mitochondrial Leucine Catabolic Pathway in Plants1

3-Methylcrotonyl-coenzyme A carboxylase (MCCase) is a mitochondrial biotin-containing enzyme whose metabolic function is not well understood in plants. In soybean (Glycine max) seedlings the organ-specific and developmentally induced changes in MCCase expression are regulated by mechanisms that control the accumulation of MCCase mRNA and the activity of the enzyme. During soybean cotyledon development, when seed-storage proteins are degraded, leucine (Leu) accumulation peaks transiently at 8 d after planting. The coincidence between peak MCCase expression and the decline in Leu content provides correlative evidence that MCCase is involved in the mitochondrial catabolism of Leu. Direct evidence for this conclusion was obtained from radiotracer metabolic studies using extracts from isolated mitochondria. These experiments traced the metabolic fate of [U-14C]Leu and NaH14CO3, the latter of which was incorporated into methylglutaconyl-coenzyme A (CoA) via MCCase. These studies directly demonstrate that plant mitochondria can catabolize Leu via the following scheme: Leu → α-ketoisocaproate → isovaleryl-CoA → 3-methylcrotonyl-CoA → 3-methylglutaconyl-CoA → 3-hydroxy-3-methylglutaryl-CoA → acetoacetate + acetyl-CoA. These findings demonstrate for the first time, to our knowledge, that the enzymes responsible for Leu catabolism are present in plant mitochondria. We conclude that a primary metabolic role of MCCase in plants is the catabolism of Leu.

Ethylene-Mediated Phospholipid Catabolic Pathway in Glucose-Starved Carrot Suspension Cells1

Glucose (Glc) starvation of suspension-cultured carrot (Daucus carota L.) cells resulted in sequential activation of phospholipid catabolic enzymes. Among the assayed enzymes involved in the degradation, phospholipase D (PLD) and lipolytic acyl hydrolase were activated at the early part of starvation, and these activities were followed by β-oxidation and the glyoxylate cycle enzymes in order. The activity of PLD and lipolytic acyl hydrolase was further confirmed by in vivo-labeling experiments. It was demonstrated that Glc added to a medium containing starving cells inhibited the phospholipid catabolic activities, indicating that phospholipid catabolism is negatively regulated by Glc. There was a burst of ethylene production 6 h after starvation. Ethylene added exogeneously to a Glc-sufficient medium activated PLD, indicating that ethylene acts as an element in the signal transduction pathway leading from Glc depletion to PLD activation. Activation of lipid peroxidation, suggestive of cell death, occurred immediately after the decrease of the phospholipid degradation, suggesting that the observed phospholipid catabolic pathway is part of the metabolic strategies by which cells effectively survive under Glc starvation.

TAT peptide on the surface of liposomes affords their efficient intracellular delivery even at low temperature and in the presence of metabolic inhibitors

To achieve an efficient intracellular drug and DNA delivery, attempts were made to target microparticulate drug carriers into cytoplasm bypassing the endocytotic pathway. TAT peptides derived from the HIV-1 TAT protein facilitate intracellular delivery of proteins and small colloidal particles. We demonstrated that relatively large drug carriers, such as 200-nm liposomes, can also be delivered into cells by TAT peptide attached to the liposome surface. Liposomes were fluorescently labeled with membranotropic rhodamine-phosphatidylethanolamine or by entrapping FITC-dextran. Incubation of fluorescent TAT liposomes with mouse Lewis lung carcinoma cells, human breast tumor BT20 cells, and rat cardiac myocyte H9C2 results in intracellular localization of certain liposomes. Steric hindrances for TAT peptide⋅cell interaction (attachment of TAT directly to the liposome surface without spacer or the presence of a high MW polyethylene glycol on the liposome surface) abolish liposome internalization, evidencing the importance of direct contact of TAT peptide with the cell surface. Low temperature or metabolic inhibitors, sodium azide or iodoacetamide, have little influence on the translocation of TAT liposomes into cells, confirming the energy-independent character of this process. The approach may have important implications for drug delivery directly into cell cytoplasm.

Metabolic engineering of essential oil yield and composition in mint by altering expression of deoxyxylulose phosphate reductoisomerase and menthofuran synthase

Peppermint (Mentha × piperita L.) was independently transformed with a homologous sense version of the 1-deoxy-d-xylulose-5-phosphate reductoisomerase cDNA and with a homologous antisense version of the menthofuran synthase cDNA, both driven by the CaMV 35S promoter. Two groups of transgenic plants were regenerated in the reductoisomerase experiments, one of which remained normal in appearance and development; another was deficient in chlorophyll production and grew slowly. Transgenic plants of normal appearance and growth habit expressed the reductoisomerase transgene strongly and constitutively, as determined by RNA blot analysis and direct enzyme assay, and these plants accumulated substantially more essential oil (about 50% yield increase) without change in monoterpene composition compared with wild-type. Chlorophyll-deficient plants did not afford detectable reductoisomerase mRNA or enzyme activity and yielded less essential oil than did wild-type plants, indicating cosuppression of the reductoisomerase gene. Plants transformed with the antisense version of the menthofuran synthase cDNA were normal in appearance but produced less than half of this undesirable monoterpene oil component than did wild-type mint grown under unstressed or stressed conditions. These experiments demonstrate that essential oil quantity and quality can be regulated by metabolic engineering. Thus, alteration of the committed step of the mevalonate-independent pathway for supply of terpenoid precursors improves flux through the pathway that leads to increased monoterpene production, and antisense manipulation of a selected downstream monoterpene biosynthetic step leads to improved oil composition.

The dwarf-1 (dt) Mutant of Zea mays blocks three steps in the gibberellin-biosynthetic pathway.

In plants, gibberellin (GA)-responding mutants have been used as tools to identify the genes that control specific steps in the GA-biosynthetic pathway. They have also been used to determine which native GAs are active per se, i.e., further metabolism is not necessary for bioactivity. We present metabolic evidence that the D1 gene of maize (Zea mays L.) controls the three biosynthetic steps: GA20 to GA1, Ga20 to GA5, and GA5 to GA3. We also present evidence that three gibberellins, GA1, GA5, and GA3, have per se activity in stimulating shoot elongation in maize. The metabolic evidence comes from the injection of [17-13C,3H]GA20 and [17-13C,3H]GA5 into seedlings of d1 and controls (normal and d5), followed by isolation and identification of the 13C-labeled metabolites by full-scan GC-MS and Kovats retention index. For the controls, GA20 was metabolized to GA1,GA3, and GA5; GA5 was metabolized to GA3. For the d1 mutant, GA20 was not metabolized to GA1, GA3, or to GA5, and GA5 was not metabolized to GA3. The bioassay evidence is based on dosage response curves using d1 seedlings for assay. GA1, GA3, and GA5 had similar bioactivities, and they were 10-times more active than GA20.