Purification, characterization and mode of action of PdhR, the transcriptional repressor of the pdhR-aceEF-lpd operon of Escherichia coli

Predictive Biophysical Neural Network Modeling of a Compendium of in vivo Transcription Factor DNA Binding Profiles for Escherichia coli

The DNA binding of most Escherichia coli Transcription Factors (TFs) has not been comprehensively mapped, and few have models that can quantitatively predict binding affinity. We report the global mapping of in vivo DNA binding for 139 E. coli TFs using ChIP-Seq. We used these data to train BoltzNet, a novel neural network that predicts TF binding energy from DNA sequence. BoltzNet mirrors a quantitative biophysical model and provides directly interpretable predictions genome-wide at nucleotide resolution. We used BoltzNet to quantitatively design novel binding sites, which we validated with biophysical experiments on purified protein. We have generated models for 125 TFs that provide insight into global features of TF binding, including clustering of sites, the role of accessory bases, the relevance of weak sites, and the background affinity of the genome. Our paper provides new paradigms for studying TF-DNA binding and for the development of biophysically motivated neural networks.

Lighting up Pyruvate Metabolism in Saccharomyces cerevisiae by a Genetically Encoded Fluorescent Biosensor

Monitoring intracellular pyruvate is useful for the exploration of fundamental metabolism and for guiding the construction of yeast cell factories for chemical production. Here, we employed a genetically encoded fluorescent Pyronic biosensor to light up the pyruvate metabolic state in the cytoplasm, nucleus, and mitochondria of Saccharomyces cerevisiae BY4741. A strong correlation was observed between the pyruvate fluctuation in mitochondria and cytoplasm when exposed to different metabolites. Further metabolic analysis of pyruvate uptake and glycolytic dynamics showed that glucose and fructose dose-dependently activated cytoplasmic pyruvate levels more effectively than direct exposure to pyruvate. Meanwhile, the Pyronic biosensor could visually distinguish phenotypes of the wild-type S. cerevisiae BY4741 and the pyruvate-hyperproducing S. cerevisiae TAM at a single-cell resolution, having the potential for high-throughput screening. Overall, Pyronic biosensors targeting different suborganelles contribute to mapping and studying the central carbon metabolism in-depth and guide the design and construction of yeast cell factories.

Transcriptome Analysis Reveals the Effect of PdhR in Plesiomonas shigelloides

The pyruvate dehydrogenase complex regulator (PdhR) was originally identified as a repressor of the pdhR-aceEF-lpd operon, which encodes the pyruvate dehydrogenase complex (PDHc) and PdhR itself. According to previous reports, PdhR plays a regulatory role in the physiological and metabolic pathways of bacteria. At present, the function of PdhR in Plesiomonas shigelloides is still poorly understood. In this study, RNA sequencing (RNA-Seq) of the wild-type strain and the ΔpdhR mutant strains was performed for comparison to identify the PdhR-controlled pathways, revealing that PdhR regulates ~7.38% of the P. shigelloides transcriptome. We found that the deletion of pdhR resulted in the downregulation of practically all polar and lateral flagella genes in P. shigelloides; meanwhile, motility assay and transmission electron microscopy (TEM) confirmed that the ΔpdhR mutant was non-motile and lacked flagella. Moreover, the results of RNA-seq and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) showed that PdhR positively regulated the expression of the T3SS cluster, and the ΔpdhR mutant significantly reduced the ability of P. shigelloides to infect Caco-2 cells compared with the WT. Consistent with previous research, pyruvate-sensing PdhR directly binds to its promoter and inhibits pdhR-aceEF-lpd operon expression. In addition, we identified two additional downstream genes, metR and nuoA, that are directly negatively regulated by PdhR. Furthermore, we also demonstrated that ArcA was identified as being located upstream of pdhR and lpdA and directly negatively regulating their expression. Overall, we revealed the function and regulatory pathway of PdhR, which will allow for a more in-depth investigation into P. shigelloides pathogenicity as well as the complex regulatory network.

Escherichia coli aceE variants coding pyruvate dehydrogenase improve the generation of pyruvate-derived acetoin

Several chromosomally expressed AceE variants were constructed in Escherichia coli ΔldhA ΔpoxB ΔppsA and compared using glucose as the sole carbon source. These variants were examined in shake flask cultures for growth rate, pyruvate accumulation, and acetoin production via heterologous expression of the budA and budB genes from Enterobacter cloacae ssp. dissolvens. The best acetoin-producing strains were subsequently studied in controlled batch culture at the one-liter scale. PDH variant strains attained up to four-fold greater acetoin than the strain expressing the wild-type PDH. In a repeated batch process, the H106V PDH variant strain attained over 43 g/L of pyruvate-derived products, acetoin (38.5 g/L) and 2R,3R-butanediol (5.0 g/L), corresponding to an effective concentration of 59 g/L considering the dilution. The acetoin yield from glucose was 0.29 g/g with a volumetric productivity of 0.9 g/L·h (0.34 g/g and 1.0 g/L·h total products). The results demonstrate a new tool in pathway engineering, the modification of a key metabolic enzyme to improve the formation of a product via a kinetically slow, introduced pathway. Direct modification of the pathway enzyme offers an alternative to promoter engineering in cases where the promoter is involved in a complex regulatory network.

Optimizing E. coli as a formatotrophic platform for bioproduction via the reductive glycine pathway

Microbial C1 fixation has a vast potential to support a sustainable circular economy. Hence, several biotechnologically important microorganisms have been recently engineered for fixing C1 substrates. However, reports about C1-based bioproduction with these organisms are scarce. Here, we describe the optimization of a previously engineered formatotrophic Escherichia coli strain. Short-term adaptive laboratory evolution enhanced biomass yield and accelerated growth of formatotrophic E. coli to 3.3 g-CDW/mol-formate and 6 h doubling time, respectively. Genome sequence analysis revealed that manipulation of acetate metabolism is the reason for better growth performance, verified by subsequent reverse engineering of the parental E. coli strain. Moreover, the improved strain is capable of growing to an OD₆₀₀ of 22 in bioreactor fed-batch experiments, highlighting its potential use for industrial bioprocesses. Finally, demonstrating the strain's potential to support a sustainable, formate-based bioeconomy, lactate production from formate was engineered. The optimized strain generated 1.2 mM lactate -10% of the theoretical maximum- providing the first proof-of-concept application of the reductive glycine pathway for bioproduction.

Is energy excess the initial trigger of carbon overflow metabolism? Transcriptional network response of carbon-limited Escherichia coli to transient carbon excess

Background

Escherichia coli adapted to carbon-limiting conditions is generally geared for energy-efficient carbon utilization. This includes also the efficient utilization of glucose, which serves as a source for cellular building blocks as well as energy. Thus, catabolic and anabolic functions are balanced under these conditions to minimize wasteful carbon utilization. Exposure to glucose excess interferes with the fine-tuned coupling of anabolism and catabolism leading to the so-called carbon overflow metabolism noticeable through acetate formation and eventually growth inhibition.

Results

Cellular adaptations towards sudden but timely limited carbon excess conditions were analyzed by exposing slow-growing cells in steady state glucose-limited continuous culture to a single glucose pulse. Concentrations of metabolites as well as time-dependent transcriptome alterations were analyzed and a transcriptional network analysis performed to determine the most relevant transcription and sigma factor combinations which govern these adaptations. Down-regulation of genes related to carbon catabolism is observed mainly at the level of substrate uptake and downstream of pyruvate and not in between in the glycolytic pathway. It is mainly accomplished through the reduced activity of CRP-cAMP and through an increased influence of phosphorylated ArcA. The initiated transcriptomic change is directed towards down-regulation of genes, which contribute to active movement, carbon uptake and catabolic carbon processing, in particular to down-regulation of genes which contribute to efficient energy generation. Long-term changes persisting after glucose depletion and consumption of acetete encompassed reduced expression of genes related to active cell movement and enhanced expression of genes related to acid resistance, in particular acid resistance system 2 (GABA shunt) which can be also considered as an inefficient bypass of the TCA cycle.

Conclusions

Our analysis revealed that the major part of the trancriptomic response towards the glucose pulse is not directed towards enhanced cell proliferation but towards protection against excessive intracellular accumulation of potentially harmful concentration of metabolites including among others energy rich compounds such as ATP. Thus, resources are mainly utilized to cope with “overfeeding” and not for growth including long-lasting changes which may compromise the cells future ability to perform optimally under carbon-limiting conditions (reduced motility and ineffective substrate utilization).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01787-4.

Acetate overflow metabolism regulates a major metabolic shift after glucose depletion in Escherichia coli

Acetate overflow refers to the metabolism by which a large part of carbon incorporated as glucose into Escherichia coli cells is catabolized and excreted as acetate into the medium. We previously found that mutants for the acetate overflow pathway enzymes phosphoacetyltransferase (Pta) and acetate kinase (AckA) showed significant diauxic growth after glucose depletion in E. coli. Here, we analyzed the underlying mechanism in the pta mutant. Proteomic and other analyses revealed an increase in pyruvate dehydrogenase complex subunits and a decrease in glyoxylate shunt enzymes, which resulted from pyruvate accumulation. Since restoration of these enzyme levels by overexpressing PdhR (pyruvate-sensing transcription factor) or deleting iclR (gene encoding a pyruvate- and glyoxylate-sensing transcription factor) alleviated the growth lag of the pta mutant after glucose depletion, these changes were considered as the reason for the phenotype. Given the evidence for decreased coenzyme A (HS-CoA) levels in the pta mutant, the growth inhibition after glucose depletion was partly explained by limited availability of HS-CoA in the cell. The findings provide insights into the role of acetate overflow in metabolic regulation, which may be useful for biotechnological applications.

A genetic screen to identify factors affected by undecaprenyl phosphate recycling uncovers novel connections to morphogenesis in Escherichia coli

Undecaprenyl phosphate (Und-P) is an essential lipid carrier that ferries cell wall intermediates across the cytoplasmic membrane in bacteria. Und-P is generated by dephosphorylating undecaprenyl pyrophosphate (Und-PP). In Escherichia coli, BacA, PgpB, YbjG, and LpxT dephosphorylate Und-PP and are conditionally essential. To identify vulnerabilities that arise when Und-P metabolism is defective, we developed a genetic screen for synthetic interactions which, in combination with ΔybjG ΔlpxT ΔbacA, are lethal or reduce fitness. The screen uncovered novel connections to cell division, DNA replication/repair, signal transduction, and glutathione metabolism. Further analysis revealed several new morphogenes; loss of one of these, qseC, caused cells to enlarge and lyse. QseC is the sensor kinase component of the QseBC two-component system. Loss of QseC causes overactivation of the QseB response regulator by PmrB cross-phosphorylation. Here, we show that deleting qseB completely reverses the shape defect of ΔqseC cells, as does overexpressing rprA (a small RNA). Surprisingly, deleting pmrB only partially suppressed qseC-related shape defects. Thus, QseB is activated by multiple factors in QseC's absence and prior functions ascribed to QseBC may originate from cell wall defects. Altogether, our findings provide a framework for identifying new determinants of cell integrity that could be targeted in future therapies.

Expanded roles of pyruvate-sensing PdhR in transcription regulation of the Escherichia coli K-12 genome: fatty acid catabolism and cell motility

The transcription factor PdhR has been recognized as the master regulator of the pyruvate catabolism pathway in Escherichia coli, including both NAD-linked oxidative decarboxylation of pyruvate to acetyl-CoA by PDHc (pyruvate dehydrogenase complex) and respiratory electron transport of NADH to oxygen by Ndh-CyoABCD enzymes. To identify the whole set of regulatory targets under the control of pyruvate-sensing PdhR, we performed genomic SELEX (gSELEX) screening in vitro. A total of 35 PdhR-binding sites were identified along the E. coli K-12 genome, including previously identified targets. Possible involvement of PdhR in regulation of the newly identified target genes was analysed in detail by gel shift assay, RT-qPCR and Northern blot analysis. The results indicated the participation of PdhR in positive regulation of fatty acid degradation genes and negative regulation of cell mobility genes. In fact, GC analysis indicated an increase in free fatty acids in the mutant lacking PdhR. We propose that PdhR is a bifunctional global regulator for control of a total of 16–23 targets, including not only the genes involved in central carbon metabolism but also some genes for the surrounding pyruvate-sensing cellular pathways such as fatty acid degradation and flagella formation. The activity of PdhR is controlled by pyruvate, the key node between a wide variety of metabolic pathways, including generation of metabolic energy and cell building blocks.

A highly responsive pyruvate sensor reveals pathway-regulatory role of the mitochondrial pyruvate carrier MPC

Mitochondria generate ATP and building blocks for cell growth and regeneration, using pyruvate as the main substrate. Here we introduce PyronicSF, a user-friendly GFP-based sensor of improved dynamic range that enables real-time subcellular quantitation of mitochondrial pyruvate transport, concentration and flux. We report that cultured mouse astrocytes maintain mitochondrial pyruvate in the low micromolar range, below cytosolic pyruvate, which means that the mitochondrial pyruvate carrier MPC is poised to exert ultrasensitive control on the balance between respiration and anaplerosis/gluconeogenesis. The functionality of the sensor in living tissue is demonstrated in the brain of Drosophila melanogaster larvae. Mitochondrial subpopulations are known to coexist within a given cell, which differ in their morphology, mobility, membrane potential, and vicinity to other organelles. The present tool can be used to investigate how mitochondrial diversity relates to metabolism, to study the role of MPC in disease, and to screen for small-molecule MPC modulators.

Systems assessment of transcriptional regulation on central carbon metabolism by Cra and CRP

Two major transcriptional regulators of carbon metabolism in bacteria are Cra and CRP. CRP is considered to be the main mediator of catabolite repression. Unlike for CRP, in vivo DNA binding information of Cra is scarce. Here we generate and integrate ChIP-exo and RNA-seq data to identify 39 binding sites for Cra and 97 regulon genes that are regulated by Cra in Escherichia coli. An integrated metabolic-regulatory network was formed by including experimentally-derived regulatory information and a genome-scale metabolic network reconstruction. Applying analysis methods of systems biology to this integrated network showed that Cra enables optimal bacterial growth on poor carbon sources by redirecting and repressing glycolysis flux, by activating the glyoxylate shunt pathway, and by activating the respiratory pathway. In these regulatory mechanisms, the overriding regulatory activity of Cra over CRP is fundamental. Thus, elucidation of interacting transcriptional regulation of core carbon metabolism in bacteria by two key transcription factors was possible by combining genome-wide experimental measurement and simulation with a genome-scale metabolic model.

Elucidating the Role and Regulation of a Lactate Permease as Lactate Transporter in Bacillus coagulans DSM1

A key feature of Bacillus coagulans is its ability to produce l-lactate via homofermentative metabolism. A putative lactate permease-encoding gene (lutP) and the gene encoding its regulator (lutR) were identified in one operon in B. coagulans strains. LutP orthologs are highly conserved and located adjacent to the gene cluster related to lactate utilization in most lactate-utilizing microorganisms. However, no lactate utilization genes were found adjacent to lutP in all sequenced B. coagulans strains. The stand-alone presence of lutP in l-lactate producers indicates that it may have functions in lactate production. In this study, B. coagulans DSM1 was used as a representative strain, and the critical roles of LutP and its regulation were described. Transport property assays showed that LutP was essential for lactate uptake. Its regulator LutR directly interacted with the lutP-lutR intergenic region, and lutP transcription was activated by l-lactate via regulation by LutR. A biolayer interferometry assay further confirmed that LutR bound to an 11-bp inverted repeat in the intergenic region, and lutP transcription began when the binding of LutR to the lutP upstream sequence was inhibited. We conclusively showed that lutP encodes a functional lactate permease in B. coagulans IMPORTANCE Lactate-utilizing strains require lactate permease (LutP) to transport lactate into cells. Bacillus coagulans LutP is a previously uncharacterized lactate permease with no lactate utilization genes situated either adjacent to or remotely from it. In this study, an active lactate permease in an l-lactate producer, B. coagulans DSM1, was identified. Lactate supplementation regulated the expression of lactate permease. This study presents physiological evidence of the presence of a lactate transporter in B. coagulans Our findings indicate a potential target for the engineering of strains in order to improve their fermentation characteristics.

Construction and evolution of an Escherichia coli strain relying on nonoxidative glycolysis for sugar catabolism

We constructed an Escherichia coli strain that does not use glycolysis for sugar catabolism. Instead, it uses the synthetic nonoxidative glycolysis cycle to directly synthesize stoichiometric amounts of the two-carbon building block (acetyl-CoA), which is then converted to three-carbon metabolites to support growth. The resulting strain grows aerobically in glucose minimal medium and can achieve near-complete carbon conservation in the production of acetyl-CoA–derived products during anaerobic fermentation. This strain improves the theoretical carbon yield from 66.7% to 100% in acetyl-CoA–derived product formation.

The Embden–Meyerhoff–Parnas (EMP) pathway, commonly known as glycolysis, represents the fundamental biochemical infrastructure for sugar catabolism in almost all organisms, as it provides key components for biosynthesis, energy metabolism, and global regulation. EMP-based metabolism synthesizes three-carbon (C3) metabolites before two-carbon (C2) metabolites and must emit one CO₂ in the synthesis of the C2 building block, acetyl-CoA, a precursor for many industrially important products. Using rational design, genome editing, and evolution, here we replaced the native glycolytic pathways in Escherichia coli with the previously designed nonoxidative glycolysis (NOG), which bypasses initial C3 formation and directly generates stoichiometric amounts of C2 metabolites. The resulting strain, which contains 11 gene overexpressions, 10 gene deletions by design, and more than 50 genomic mutations (including 3 global regulators) through evolution, grows aerobically in glucose minimal medium but can ferment anaerobically to products with nearly complete carbon conservation. We confirmed that the strain metabolizes glucose through NOG by ¹³C tracer experiments. This redesigned E. coli strain represents a different approach for carbon catabolism and may serve as a useful platform for bioproduction.

Carbon Monoxide Gas Is Not Inert, but Global, in Its Consequences for Bacterial Gene Expression, Iron Acquisition, and Antibiotic Resistance

AIMS

Carbon monoxide is a respiratory poison and gaseous signaling molecule. Although CO-releasing molecules (CORMs) deliver CO with temporal and spatial specificity in mammals, and are proven antimicrobial agents, we do not understand the modes of CO toxicity. Our aim was to explore the impact of CO gas per se, without intervention of CORMs, on bacterial physiology and gene expression.

RESULTS

We used tightly controlled chemostat conditions and integrated transcriptomic datasets with statistical modeling to reveal the global effects of CO. CO is known to inhibit bacterial respiration, and we found expression of genes encoding energy-transducing pathways to be significantly affected via the global regulators, Fnr, Arc, and PdhR. Aerobically, ArcA-the response regulator-is transiently phosphorylated and pyruvate accumulates, mimicking anaerobiosis. Genes implicated in iron acquisition, and the metabolism of sulfur amino acids and arginine, are all perturbed. The global iron-related changes, confirmed by modulation of activity of the transcription factor Fur, may underlie enhanced siderophore excretion, diminished intracellular iron pools, and the sensitivity of CO-challenged bacteria to metal chelators. Although CO gas (unlike H2S and NO) offers little protection from antibiotics, a ruthenium CORM is a potent adjuvant of antibiotic activity.

INNOVATION

This is the first detailed exploration of global bacterial responses to CO, revealing unexpected targets with implications for employing CORMs therapeutically.

CONCLUSION

This work reveals the complexity of bacterial responses to CO and provides a basis for understanding the impacts of CO from CORMs, heme oxygenase activity, or environmental sources. Antioxid. Redox Signal. 24, 1013-1028.

Biotin and Lipoic Acid: Synthesis, Attachment, and Regulation

Two vitamins, biotin and lipoic acid, are essential in all three domains of life. Both coenzymes function only when covalently attached to key metabolic enzymes. There they act as "swinging arms" that shuttle intermediates between two active sites (= covalent substrate channeling) of key metabolic enzymes. Although biotin was discovered over 100 years ago and lipoic acid 60 years ago, it was not known how either coenzyme is made until recently. In Escherichia coli the synthetic pathways for both coenzymes have now been worked out for the first time. The late steps of biotin synthesis, those involved in assembling the fused rings, were well described biochemically years ago, although recent progress has been made on the BioB reaction, the last step of the pathway in which the biotin sulfur moiety is inserted. In contrast, the early steps of biotin synthesis, assembly of the fatty acid-like "arm" of biotin were unknown. It has now been demonstrated that the arm is made by using disguised substrates to gain entry into the fatty acid synthesis pathway followed by removal of the disguise when the proper chain length is attained. The BioC methyltransferase is responsible for introducing the disguise, and the BioH esterase is responsible for its removal. In contrast to biotin, which is attached to its cognate proteins as a finished molecule, lipoic acid is assembled on its cognate proteins. An octanoyl moiety is transferred from the octanoyl acyl carrier protein of fatty acid synthesis to a specific lysine residue of a cognate protein by the LipB octanoyltransferase followed by sulfur insertion at carbons C-6 and C-8 by the LipA lipoyl synthetase. Assembly on the cognate proteins regulates the amount of lipoic acid synthesized, and, thus, there is no transcriptional control of the synthetic genes. In contrast, transcriptional control of the biotin synthetic genes is wielded by a remarkably sophisticated, yet simple, system, exerted through BirA, a dual-function protein that both represses biotin operon transcription and ligates biotin to its cognate proteins.

Tricarboxylic Acid Cycle and Glyoxylate Bypass

The tricarboxylic acid (TCA) cycle plays two essential roles in metabolism. First, under aerobic conditions the cycle is responsible for the total oxidation of acetyl-CoA that is derived mainly from the pyruvate produced by glycolysis. Second, TCA cycle intermediates are required in the biosynthesis of several amino acids. Although the TCA cycle has long been considered a "housekeeping" pathway in Escherichia coli and Salmonella enterica, the pathway is highly regulated at the transcriptional level. Much of this control is exerted in response to respiratory conditions. The TCA cycle gene-protein relationship and mutant phenotypes have been well studied, although a few loose ends remain. The realization that a "shadow" TCA cycle exists that proceeds through methylcitrate has cleared up prior ambiguities. The glyoxylate bypass has long been known to be essential for growth on carbon sources such as acetate or fatty acids because this pathway allowsnet conversion of acetyl-CoA to metabolic intermediates. Strains lacking this pathway fail to grow on these carbon sources, since acetate carbon entering the TCA cycle is quantitatively lost as CO2 resulting in the lack of a means to replenish the dicarboxylic acids consumed in amino acid biosynthesis. The TCA cycle gene-protein relationship and mutant phenotypes have been well studied, although the identity of the small molecule ligand that modulates transcriptional control of the glyoxylate cycle genes by binding to the IclR repressor remains unknown. The activity of the cycle is also exerted at the enzyme level by the reversible phosphorylation of the TCA cycle enzyme isocitrate dehydrogenase catalyzed by a specific kinase/phosphatase to allow isocitratelyase to compete for isocitrate and cleave this intermediate to glyoxylate and succinate.

Biotin and Lipoic Acid: Synthesis, Attachment, and Regulation

Two vitamins, biotin and lipoic acid, are essential in all three domains of life. Both coenzymes function only when covalently attached to key metabolic enzymes. There they act as "swinging arms" that shuttle intermediates between two active sites (= covalent substrate channeling) of key metabolic enzymes. Although biotin was discovered over 100 years ago and lipoic acid was discovered 60 years ago, it was not known how either coenzyme is made until recently. In Escherichia coli the synthetic pathways for both coenzymes have now been worked out for the first time. The late steps of biotin synthesis, those involved in assembling the fused rings, were well described biochemically years ago, although recent progress has been made on the BioB reaction, the last step of the pathway, in which the biotin sulfur moiety is inserted. In contrast, the early steps of biotin synthesis, assembly of the fatty acid-like "arm" of biotin, were unknown. It has now been demonstrated that the arm is made by using disguised substrates to gain entry into the fatty acid synthesis pathway followed by removal of the disguise when the proper chain length is attained. The BioC methyltransferase is responsible for introducing the disguise and the BioH esterase for its removal. In contrast to biotin, which is attached to its cognate proteins as a finished molecule, lipoic acid is assembled on its cognate proteins. An octanoyl moiety is transferred from the octanoyl-ACP of fatty acid synthesis to a specific lysine residue of a cognate protein by the LipB octanoyl transferase, followed by sulfur insertion at carbons C6 and C8 by the LipA lipoyl synthetase. Assembly on the cognate proteins regulates the amount of lipoic acid synthesized, and thus there is no transcriptional control of the synthetic genes. In contrast, transcriptional control of the biotin synthetic genes is wielded by a remarkably sophisticated, yet simple, system exerted through BirA, a dual-function protein that both represses biotin operon transcription and ligates biotin to its cognate protein.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

Rapid, randomized development of genetically encoded FRET sensors for small molecules

A generally applicable protocol for random and yet efficient construction of genetically encoded FRET sensors for small molecules was established.

Benzoate mediates the simultaneous repression of anaerobic 4-methylbenzoate and succinate utilization in Magnetospirillum sp. strain pMbN1

Background

At high concentrations of organic substrates, microbial utilization of preferred substrates (i.e., supporting fast growth) often results in diauxic growth with hierarchical substrate depletion. Unlike the carbon catabolite repression-mediated discriminative utilization of carbohydrates, the substrate preferences of non-carbohydrate-utilizing bacteria for environmentally relevant compound classes (e.g., aliphatic or aromatic acids) are rarely investigated. The denitrifying alphaproteobacterium Magnetospirillum sp. strain pMbN1 anaerobically degrades a wide variety of aliphatic and aromatic compounds and is unique for anaerobic degradation of 4-methylbenzoate. The latter proceeds via a distinct reaction sequence analogous to the central anaerobic benzoyl-CoA pathway to intermediates of central metabolism. Considering the presence of these two different anaerobic “aromatic ring degrading” pathways, substrate preferences of Magnetospirillum sp. strain pMbN1 were investigated. Anaerobic growth and substrate consumption were monitored in binary and ternary mixtures of 4-methylbenzoate, benzoate and succinate, in conjuction with time-resolved abundance profiling of selected transcripts and/or proteins related to substrate uptake and catabolism.

Results

Diauxic growth with benzoate preference was observed for binary and ternary substrate mixtures containing 4-methylbenzoate and succinate (despite adaptation of Magnetospirillum sp. strain pMbN1 to one of the latter two substrates). On the contrary, 4-methylbenzoate and succinate were utilized simultaneously from a binary mixture, as well as after benzoate depletion from the ternary mixture. Apparently, simultaneous repression of 4-methylbenzoate and succinate utilization from the ternary substrate mixture resulted from (i) inhibition of 4-methylbenzoate uptake, and (ii) combined inhibition of succinate uptake (via the two transporters DctPQM and DctA) and succinate conversion to acetyl-CoA (via pyruvate dehydrogenase). The benzoate-mediated repression of C₄-dicarboxylate utilization in Magnetospirillum sp. strain pMbN1 differs from that recently described for “Aromatoleum aromaticum” EbN1 (involving only DctPQM).

Conclusions

The preferential or simultaneous utilization of benzoate and other aromatic acids from mixtures with aliphatic acids may represent a more common nutritional behavior among (anaerobic) degradation specialist than previously thought. Preference of Magnetospirillum sp. strain pMbN1 for benzoate from mixtures with 4-methylbenzoate, and thus temporal separation of the benzoyl-CoA (first) and 4-methylbenzoyl-CoA (second) pathway, may reflect a catabolic tuning towards metabolic efficiency and the markedly broader range of aromatic substrates feeding into the central anaerobic benzoyl-CoA pathway.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-014-0269-4) contains supplementary material, which is available to authorized users.

The metabolic potential of Escherichia coli BL21 in defined and rich medium

BACKGROUND

The proteome reflects the available cellular machinery to deal with nutrients and environmental challenges. The most common E. coli strain BL21 growing in different, commonly employed media was evaluated using a detailed quantitative proteome analysis.

RESULTS

The presence of preformed biomass precursor molecules in rich media such as Luria Bertani supported rapid growth concomitant to acetate formation and apparently unbalanced abundances of central metabolic pathway enzymes, e.g. high levels of lower glycolytic pathway enzymes as well as pyruvate dehydrogenase, and low levels of TCA cycle and high levels of the acetate forming enzymes Pta and AckA. The proteome of cells growing exponentially in glucose-supplemented mineral salt medium was dominated by enzymes of amino acid synthesis pathways, contained more balanced abundances of central metabolic pathway enzymes, and a lower portion of ribosomal and other translational proteins. Entry into stationary phase led to a reconstruction of the bacterial proteome by increasing e.g. the portion of proteins required for scavenging rare nutrients and general cell protection. This proteomic reconstruction during entry into stationary phase was more noticeable in cells growing in rich medium as they have a greater reservoir of recyclable proteins from the translational machinery.

CONCLUSIONS

The proteomic comparison of cells growing exponentially in different media reflected the antagonistic and competitive regulation of central metabolic pathways through the global transcriptional regulators Cra, Crp, and ArcA. For example, the proteome of cells growing exponentially in rich medium was consistent with a dominating role of phosphorylated ArcA most likely a result from limitations in reoxidizing reduced quinones in the respiratory chain under these growth conditions. The proteomic alterations of exponentially growing cells into stationary phase cells were consistent with stringent-like and stationary phase responses and a dominating control through DksA-ppGpp and RpoS.

Imaging mitochondrial flux in single cells with a FRET sensor for pyruvate

Mitochondrial flux is currently accessible at low resolution. Here we introduce a genetically-encoded FRET sensor for pyruvate, and methods for quantitative measurement of pyruvate transport, pyruvate production and mitochondrial pyruvate consumption in intact individual cells at high temporal resolution. In HEK293 cells, neurons and astrocytes, mitochondrial pyruvate uptake was saturated at physiological levels, showing that the metabolic rate is determined by intrinsic properties of the organelle and not by substrate availability. The potential of the sensor was further demonstrated in neurons, where mitochondrial flux was found to rise by 300% within seconds of a calcium transient triggered by a short theta burst, while glucose levels remained unaltered. In contrast, astrocytic mitochondria were insensitive to a similar calcium transient elicited by extracellular ATP. We expect the improved resolution provided by the pyruvate sensor will be of practical interest for basic and applied researchers interested in mitochondrial function.

Genome-scale analysis of escherichia coli FNR reveals complex features of transcription factor binding

FNR is a well-studied global regulator of anaerobiosis, which is widely conserved across bacteria. Despite the importance of FNR and anaerobiosis in microbial lifestyles, the factors that influence its function on a genome-wide scale are poorly understood. Here, we report a functional genomic analysis of FNR action. We find that FNR occupancy at many target sites is strongly influenced by nucleoid-associated proteins (NAPs) that restrict access to many FNR binding sites. At a genome-wide level, only a subset of predicted FNR binding sites were bound under anaerobic fermentative conditions and many appeared to be masked by the NAPs H-NS, IHF and Fis. Similar assays in cells lacking H-NS and its paralog StpA showed increased FNR occupancy at sites bound by H-NS in WT strains, indicating that large regions of the genome are not readily accessible for FNR binding. Genome accessibility may also explain our finding that genome-wide FNR occupancy did not correlate with the match to consensus at binding sites, suggesting that significant variation in ChIP signal was attributable to cross-linking or immunoprecipitation efficiency rather than differences in binding affinities for FNR sites. Correlation of FNR ChIP-seq peaks with transcriptomic data showed that less than half of the FNR-regulated operons could be attributed to direct FNR binding. Conversely, FNR bound some promoters without regulating expression presumably requiring changes in activity of condition-specific transcription factors. Such combinatorial regulation may allow Escherichia coli to respond rapidly to environmental changes and confer an ecological advantage in the anaerobic but nutrient-fluctuating environment of the mammalian gut.

Regulation of gene expression by transcription factors (TFs) is key to adaptation to environmental changes. Our comprehensive, genome-scale analysis of a prototypical global TF, the anaerobic regulator FNR from Escherichia coli, leads to several novel and unanticipated insights into the influences on FNR binding genome-wide and the complex structure of bacterial regulons. We found that binding of NAPs restricts FNR binding at a subset of sites, suggesting that the bacterial genome is not freely accessible for FNR binding. Our finding that less than half of the predicted FNR binding sites were occupied in vivo further challenges the utility of using bioinformatic searches alone to predict regulon structure, reinforcing the need for experimental determination of TF binding. By correlating the occupancy data with transcriptomic data, we confirm that FNR serves as a global signal of anaerobiosis but expression of some operons in the FNR regulon require other regulators sensitive to alternative environmental stimuli. Thus, FNR binding and regulation appear to depend on both the nucleoprotein structure of the chromosome and on combinatorial binding of FNR with other regulators. Both of these phenomena are typical of TF binding in eukaryotes; our results establish that they are also features of bacterial TF binding.

More than just a metabolic regulator--elucidation and validation of new targets of PdhR in Escherichia coli

BACKGROUND

The pyruvate dehydrogenase regulator protein (PdhR) of Escherichia coli acts as a transcriptional regulator in a pyruvate dependent manner to control central metabolic fluxes. However, the complete PdhR regulon has not yet been uncovered. To achieve an extended understanding of its gene regulatory network, we combined large-scale network inference and experimental verification of results obtained by a systems biology approach.

RESULTS

22 new genes contained in two operons controlled by PdhR (previously only 20 regulatory targets in eight operons were known) were identified by analysing a large-scale dataset of E. coli from the Many Microbes Microarray Database and novel expression data from a pdhR knockout strain, as well as a PdhR overproducing strain. We identified a regulation of the glycolate utilization operon glcDEFGBA using chromatin immunoprecipitation and gel shift assays. We show that this regulation could be part of a cross-induction between genes necessary for acetate and pyruvate utilisation controlled through PdhR. Moreover, a link of PdhR regulation to the replication machinery of the cell via control of the transcription of the dcw-cluster was verified in experiments. This augments our knowledge of the functions of the PdhR-regulon and demonstrates its central importance for further cellular processes in E. coli.

CONCLUSIONS

We extended the PdhR regulon by 22 new genes contained in two operons and validated the regulation of the glcDEFGBA operon for glycolate utilisation and the dcw-cluster for cell division proteins experimentally. Our results provide, for the first time, a plausible regulatory link between the nutritional status of the cell and cell replication mediated by PdhR.

Stationary-Phase Gene Regulation in Escherichia coli §

In their stressful natural environments, bacteria often are in stationary phase and use their limited resources for maintenance and stress survival. Underlying this activity is the general stress response, which in Escherichia coli depends on the σS (RpoS) subunit of RNA polymerase. σS is closely related to the vegetative sigma factor σ70 (RpoD), and these two sigmas recognize similar but not identical promoter sequences. During the postexponential phase and entry into stationary phase, σS is induced by a fine-tuned combination of transcriptional, translational, and proteolytic control. In addition, regulatory "short-cuts" to high cellular σS levels, which mainly rely on the rapid inhibition of σS proteolysis, are triggered by sudden starvation for various nutrients and other stressful shift conditons. σS directly or indirectly activates more than 500 genes. Additional signal input is integrated by σS cooperating with various transcription factors in complex cascades and feedforward loops. Target gene products have stress-protective functions, redirect metabolism, affect cell envelope and cell shape, are involved in biofilm formation or pathogenesis, or can increased stationary phase and stress-induced mutagenesis. This review summarizes these diverse functions and the amazingly complex regulation of σS. At the molecular level, these processes are integrated with the partitioning of global transcription space by sigma factor competition for RNA polymerase core enzyme and signaling by nucleotide second messengers that include cAMP, (p)ppGpp, and c-di-GMP. Physiologically, σS is the key player in choosing between a lifestyle associated with postexponential growth based on nutrient scavenging and motility and a lifestyle focused on maintenance, strong stress resistance, and increased adhesiveness. Finally, research with other proteobacteria is beginning to reveal how evolution has further adapted function and regulation of σS to specific environmental niches.

Reprogramming of Escherichia coli K-12 metabolism during the initial phase of transition from an anaerobic to a micro-aerobic environment

Background

Many bacteria undergo transitions between environments with differing O₂ availabilities as part of their natural lifestyles and during biotechnological processes. However, the dynamics of adaptation when bacteria experience changes in O₂ availability are understudied. The model bacterium and facultative anaerobe Escherichia coli K-12 provides an ideal system for exploring this process.

Methods and Findings

Time-resolved transcript profiles of E. coli K-12 during the initial phase of transition from anaerobic to micro-aerobic conditions revealed a reprogramming of gene expression consistent with a switch from fermentative to respiratory metabolism. The changes in transcript abundance were matched by changes in the abundances of selected central metabolic proteins. A probabilistic state space model was used to infer the activities of two key regulators, FNR (O₂ sensing) and PdhR (pyruvate sensing). The model implied that both regulators were rapidly inactivated during the transition from an anaerobic to a micro-aerobic environment. Analysis of the external metabolome and protein levels suggested that the cultures transit through different physiological states during the process of adaptation, characterized by the rapid inactivation of pyruvate formate-lyase (PFL), a slower induction of pyruvate dehydrogenase complex (PDHC) activity and transient excretion of pyruvate, consistent with the predicted inactivation of PdhR and FNR.

Conclusion

Perturbation of anaerobic steady-state cultures by introduction of a limited supply of O₂ combined with time-resolved transcript, protein and metabolite profiling, and probabilistic modeling has revealed that pyruvate (sensed by PdhR) is a key metabolic signal in coordinating the reprogramming of E. coli K-12 gene expression by working alongside the O₂ sensor FNR during transition from anaerobic to micro-aerobic conditions.

Large-scale 13C-flux analysis reveals distinct transcriptional control of respiratory and fermentative metabolism in Escherichia coli

The authors analyze the role transcription plays in regulating bacterial metabolic flux. Of 91 transcriptional regulators studied, 2/3 affect absolute fluxes, but only a small number of regulators control the partitioning of flux between different metabolic pathways.

In contrast to the canonical respiro-fermentative glucose metabolism, fully respiratory galactose metabolism depends exclusively on the PEP-glyoxylate cycle.

Of 91 transcription factors, 2/3 affect absolute fluxes, but only one controls the distribution of fluxes on galactose and nine on glucose.

Transcriptional control of hexose flux distributions is confined to the acetyl-CoA branch point.

The PEP-glyoxylate cycle is controlled by cAMP-Crp in a hexose uptake rate-dependent manner.

Focusing on central carbon metabolism of Escherichia coli, we aim here to systematically identify transcriptional regulators that control the distribution of metabolic fluxes during aerobic growth on hexoses. To assess the condition dependence of transcriptional control of flux, we selected glucose and galactose as two substrates that are highly similar, yet lead to distinct growth rates (Soupene et al, 2003), overall metabolic rates (De Anda et al, 2006; Samir El et al, 2009) and levels of catabolite repression (Hogema et al, 1998; Bettenbrock et al, 2007).

Experimentally determined fluxes (Fischer and Sauer, 2003a) during growth on glucose and galactose reveal two distinct metabolic states. On glucose, high metabolic rates lead to high overflow metabolism and respiratory fluxes through the TCA cycle. On galactose, in contrast, metabolism was much slower without overflow metabolism and respiratory fluxes exclusively through the PEP-glyoxylate cycle (Fischer and Sauer, 2003b).

To determine which transcriptional events controlled these two distinct metabolic states, we determined intracellular fluxes in 91 transcription factor mutants. These genetic perturbations primarily affected absolute fluxes but not the distribution of fluxes. The distribution of flux between glycolysis and pentose–phosphate pathway in upper metabolism, e.g., remained constant in all mutants under all conditions. Transcriptional control of the flux distribution was exclusively seen at the acetyl-CoA branch point. On glucose, nine transcription factors controlled the distribution of fluxes at this branch point, five of which (ArcA, IHFA, IHFB, PdhR, Fur) did so presumably directly through their known targets in TCA cycle and/or respiration. Without known targets in the relevant pathways, the remaining four transcription factors (GlpR, QseB, HdfR, GlcC) may act either indirectly or directly through unknown targets. On galactose, transcriptional control focused exclusively on the PEP-glyoxylate cycle. While deletion of six transcription factors (Cra, Crp, IHFA, IHFB, Mlc, NagC) abolished or reduced the PEP-glyoxylate cycle flux, we demonstrate by substrate-limited chemostat experiments, derepression of galactose uptake and show by metabolomics that five of these transcription factors act indirectly through increased cAMP concentrations that allosterically activate Crp, the only direct transcription factors that controls the PEP-glyoxylate cycle (Nanchen et al, 2008).

Overall, our absolute flux data demonstrate that control of flux splitting during growth on hexoses was confined to the acetyl-CoA branch point in E. coli. Of the 36 transcription factors known to target genes in pathways that diverge from the acetyl-CoA branch point, only one transcription factor on galactose and five plus potentially four others on glucose showed altered flux splitting. The primary focus of steady state transcriptional control on the acetyl-CoA branch point, and thus the metabolic decision between the energetically efficient respiration and the less efficient but more rapid fermentation, was recently also demonstrated with only relative flux data for Saccharomyces cerevisiae (Fendt et al, 2010).

In contrast to glucose-grown Bacillus subtilis (Fischer et al, 2005), for batch glucose-grown E. coli, none of the investigated transcription factor mutants exhibited improved biomass productivity. However, the mutants Cra, IHF A, IHF B and NagC with increased uptake rates grew much faster at almost unaltered biomass yields. As the removal of the glucose PTS-based repression with a Crr mutant also resulted in increased galactose uptake, we provide evidence that E. coli actively represses its galactose uptake at the expense of otherwise possible rapid growth.

Despite our increasing topological knowledge on regulation networks in model bacteria, it is largely unknown which of the many co-occurring regulatory events actually control metabolic function and the distribution of intracellular fluxes. Here, we unravel condition-dependent transcriptional control of Escherichia coli metabolism by large-scale ¹³C-flux analysis in 91 transcriptional regulator mutants on glucose and galactose. In contrast to the canonical respiro-fermentative glucose metabolism, fully respiratory galactose metabolism depends exclusively on the phosphoenol-pyruvate (PEP)-glyoxylate cycle. While 2/3 of the regulators directly or indirectly affected absolute flux rates, the partitioning between different pathways remained largely stable with transcriptional control focusing primarily on the acetyl-CoA branch point. Flux distribution control was achieved by nine transcription factors on glucose, including ArcA, Fur, PdhR, IHF A and IHF B, but was exclusively mediated by the cAMP-dependent Crp regulation of the PEP-glyoxylate cycle flux on galactose. Five further transcription factors affected this flux only indirectly through cAMP and Crp by increasing the galactose uptake rate. Thus, E. coli actively limits its galactose catabolism at the expense of otherwise possible faster growth.

The operonic location of auto-transcriptional repressors is highly conserved in bacteria

Bacterial genes are commonly encoded in clusters, known as operons, which share transcriptional regulatory control and often encode functionally related proteins that take part in certain biological pathways. Operons that are coregulated are known to colocalize in the genome, suggesting that their spatial organization is under selection for efficient expression regulation. However, the internal order of genes within operons is believed to be poorly conserved, and hence expression requirements are claimed to be too weak to oppose gene rearrangements. In light of these opposing views, we set out to investigate whether the internal location of the regulatory genes within operons is under selection. Our analysis shows that transcription factors (TFs) are preferentially encoded as either first or last in their operons, in the two diverged model bacteria Escherichia coli and Bacillus subtilis. In a higher resolution, we find that TFs that repress transcription of the operon in which they are encoded (autorepressors), contribute most of this signal by specific preference of the first operon position. We show that this trend is strikingly conserved throughout highly diverged bacterial phyla. Moreover, these autorepressors regulate operons that carry out highly diverse biological functions. We propose a model according to which autorepressors are selected to be located first in their operons in order to optimize transcription regulation. Specifically, the first operon position helps autorepressors to minimize leaky transcription of the operon structural genes, thus minimizing energy waste. Our analysis provides statistically robust evidence for a paradigm of bacterial autorepressor preferential operonic location. Corroborated with our suggested model, an additional layer of operon expression control that is common throughout the bacterial domain is revealed.

Novel members of the Cra regulon involved in carbon metabolism in Escherichia coli

Cra (catabolite repressor activator) is a global regulator of the genes for carbon metabolism in Escherichia coli. To gain insights into the regulatory roles of Cra, attempts were made to identify the whole set of regulation targets using an improved genomic SELEX (systematic evolution of ligands by exponential enrichment) system. Surprisingly, a total of 164 binding sites were identified for Cra, 144 (88%) of which were newly identified. The majority of known targets were included in the SELEX chip pattern. The promoters examined by the lacZ reporter assay in vivo were all regulated by Cra. These two lines of evidence indicate that a total of as many as 178 promoters are under the control of Cra. The majority of Cra targets are the genes coding for the enzymes involved in central carbon metabolism, covering all the genes for the enzymes involved in glycolysis and metabolism downstream of glycolysis, including the tricarboxylic acid (TCA) cycle and aerobic respiration. Taken together, we propose that Cra plays a key role in balancing the levels of the enzymes for carbon metabolism.

Doubling the catabolic reducing power (NADH) output of Escherichia coli fermentation for production of reduced products

Homofermentative production of reduced products requires additional reducing power output (NADH) from glucose catabolism. Anaerobic expression of the pyruvate dehydrogenase complex (PDH, encoded by aceEF-lpd, a normal aerobic operon) is able to provide the additional NADH required for production of reduced products in Escherichia coli fermentation. The multiple promoters (pflBp(1-7)) of pyruvate formate lyase (pflB) were evaluated for anaerobic expression of the aceEF-lpd operon. Four chromosomal constructs, pflBp(1-7)-aceEF-lpd, pflBp(1-6)-aceEF-lpd, pflBp(6,7)-aceEF-lpd, and pflBp6-aceEF-lpd efficiently expressed the PDH complex in anaerobically grown cells. Doubling the reducing power output was achieved when glucose was oxidized to acetyl-CoA through glycolysis and pyruvate oxidation by the anaerobically expressed PDH complex (glucose -->2 acetyl-CoA + 4 NADH). This additional reducing power output can be used for production of reduced products in anaerobic E. coli fermentation.

Bacterial adaptation through distributed sensing of metabolic fluxes

We present a large-scale differential equation model of E. coli's central metabolism and its enzymatic, transcriptional, and posttranslational regulation. This model reproduces E. coli's known physiological behavior.

We found that the interplay of known interactions in E. coli's central metabolism can indirectly recognize the presence of extracellular carbon sources through measuring intracellular metabolic flux patterns.

We found that E. coli's system-level adaptations between glycolytic and gluconeogenic carbon sources are realized on the molecular level by global feedback architectures that overarch the enzymatic and transcriptional regulatory layers.

We found that the capability for closed-loop self-regulation can emerge within metabolism itself and therefore, metabolic operation may adapt itself autonomously to changing carbon sources (not requiring upstream sensing and signaling).

Adaptations to fluctuating carbon source availability are of particular importance for bacteria. To understand these adaptations, it needs to be understood how a system's behavior emerges from the interactions between the characterized molecules (Kitano, 2002b). To attain such a system understanding of bacterial metabolic adaptations to carbon source availability, the coupling between the recognition and adjustment aspects and between the enzymatic and genetic regulatory layers must be understood. For many carbon sources, neither transmembrane sensors nor regulatory proteins with sensing function have been identified. Also, it remains unclear how multiple local regulations work together to accomplish a coherent adjustment on the systems level. In this paper, we show that (1) the interplay of the known interactions in E. coli's central metabolism is capable of recognizing carbon sources indirectly, and that (2) these molecular interactions can adjust E. coli's metabolic operation between growth on glycolytic and gluconeogenic carbon sources, and that (3) this adaptation is governed by general principles.

We hypothesized that the system-level adaptations between growth on glycolytic and gluconeogenic carbon sources are accomplished by a system-wide regulation architecture that emerges when the known enzymatic and transcriptional regulations become coupled through five transcription factor (TF)–metabolite interactions. To (1) assess whether such coupled molecular interactions can indeed work together to adapt metabolic operation, and if yes, (2) to understand this system-level adaptation in molecular-level detail, we constructed a large-scale differential equation model. The model topology comprises the Embden–Meyerhoff pathway, the tricarboxylic acid (TCA) cycle, the glyoxylate (GLX) shunt, the anaplerotic reactions, the diversion of carbon flux to the GLX shunt, the uptake of glucose, the uptake and excretion of acetate, enzymatic regulation, transcriptional regulation by four TFs, and the regulation of these TFs' activities through TF–metabolite interactions. We translated the topology into differential equations by assigning the most appropriate rate law to each interaction. The kinetic model comprises 47 ordinary differential equations and 193 parameters. Parameter values were estimated through application of the ‘divide-and-conquer approach' (Kotte and Heinemann, 2009) on published experimental steady state-omics data sets.

Model simulations reproduce E. coli's known physiological behavior in an environment with fluctuating carbon source availability. But how does the in silico cell recognize acetate without a transmembrane sensor for extracellular acetate or a TF binding to intracellular acetate? Similarly, it is unclear whether the glucose sensing function of the phosphotransferase system is the exclusive mechanism to recognize glucose, or whether this sensing function is integrated into a larger sensing architecture. The model suggests that the recognition is performed indirectly through a mechanism we termed distributed sensing of intracellular metabolic fluxes. This mechanism uses two distinct motifs, which we termed pathway usage and flux direction, to establish defined correlations between metabolic fluxes and the levels of certain, here termed flux-signaling metabolites. The binding of these metabolites to TFs propagates the flux information to the transcriptional regulatory layer. A molecular sensor for intracellular metabolic flux is thus defined as a system of regulations and enzyme kinetics, comprising (1) either of the two motifs pathway usage or flux direction and (2) the binding of the thus established flux-signaling metabolites to TF(s).

As the in silico cell establishes and uses sensors for several intracellular metabolic fluxes, the overall sensing architecture infers the present carbon sources from a pattern of metabolic fluxes and is as such of a distributed nature. The core of this sensing architecture is formed not by transmembrane sensors but by four flux sensors, which establish flux-signaling metabolites according to the two proposed general motifs. These flux sensors use intracellular metabolic flux as a means to correlate the presence of extracellular carbon sources with the levels of intracellular metabolites. The recognition of glucose through the PTS transmembrane complex is embedded as one flux sensor in this distributed sensing architecture; the other three flux sensors function without the help of transmembrane complexes.

The in silico cell achieves the coupling between recognition and adjustment through its TFs, whose activities respond to the available carbon sources and at the same time regulate the expression of target genes. This combined recognition and adjustment, centered on the four TFs, closes four global feedback loops that overarch the metabolic and genetic layers as illustrated in Figure 6. The adaptation of the in silico cell arises from the global feedback loop-embedded, flux sensor-adjusted transcriptional regulation of the four TFs, with each TF performing one part of the overall adaptation. This adaptation incorporates both the influence of the metabolic on the genetic layer, achieved through TF–metabolite interactions, and of the genetic on the metabolic layer, achieved through the impact of adjusted enzyme levels on metabolic fluxes.

The existence of the global feedback architectures challenges the conventional view that top-level regulatory proteins recognize environmental conditions and adjust downstream metabolic operation. It suggests that the capability for closed-loop self-regulation can emerge within metabolism itself and therefore, metabolic operation may adapt itself autonomously (not requiring upstream sensing and regulation) to changing carbon sources.

To conclude, the presented differential equation model of E. coli's central metabolism offers a consistent explanation of how a multitude of known molecular interactions fit into a coherent systems picture; the interactions work together like gear wheels that mesh with one another to adapt central metabolism between growth on the glycolytic substrate glucose and the gluconeogenic substrate acetate. The deduced general functional principles provide the missing link to understand system-level adaptations to carbon sources in molecular-level detail. The proposed principles fall under the umbrella of distributed flux sensing. The flux sensing mechanism entails the binding of TFs to flux-signaling metabolites, which are established through the motifs signaling of pathway usage and signaling of flux direction, and are embedded in global feedback loop architectures. These principles allow an autonomous adaptation of metabolic operation to growth in fluctuating environments.

The recognition of carbon sources and the regulatory adjustments to recognized changes are of particular importance for bacterial survival in fluctuating environments. Despite a thorough knowledge base of Escherichia coli's central metabolism and its regulation, fundamental aspects of the employed sensing and regulatory adjustment mechanisms remain unclear. In this paper, using a differential equation model that couples enzymatic and transcriptional regulation of E. coli's central metabolism, we show that the interplay of known interactions explains in molecular-level detail the system-wide adjustments of metabolic operation between glycolytic and gluconeogenic carbon sources. We show that these adaptations are enabled by an indirect recognition of carbon sources through a mechanism we termed distributed sensing of intracellular metabolic fluxes. This mechanism uses two general motifs to establish flux-signaling metabolites, whose bindings to transcription factors form flux sensors. As these sensors are embedded in global feedback loop architectures, closed-loop self-regulation can emerge within metabolism itself and therefore, metabolic operation may adapt itself autonomously (not requiring upstream sensing and signaling) to fluctuating carbon sources.

Metabolic flux control at the pyruvate node in an anaerobic Escherichia coli strain with an active pyruvate dehydrogenase

During anaerobic growth of Escherichia coli, pyruvate formate-lyase (PFL) and lactate dehydrogenase (LDH) channel pyruvate toward a mixture of fermentation products. We have introduced a third branch at the pyruvate node in a mutant of E. coli with a mutation in pyruvate dehydrogenase (PDH*) that renders the enzyme less sensitive to inhibition by NADH. The key starting enzymes of the three branches at the pyruvate node in such a mutant, PDH*, PFL, and LDH, have different metabolic potentials and kinetic properties. In such a mutant (strain QZ2), pyruvate flux through LDH was about 30%, with the remainder of the flux occurring through PFL, indicating that LDH is a preferred route of pyruvate conversion over PDH*. In a pfl mutant (strain YK167) with both PDH* and LDH activities, flux through PDH* was about 33% of the total, confirming the ability of LDH to outcompete the PDH pathway for pyruvate in vivo. Only in the absence of LDH (strain QZ3) was pyruvate carbon equally distributed between the PDH* and PFL pathways. A pfl mutant with LDH and PDH* activities, as well as a pfl ldh double mutant with PDH* activity, had a surprisingly low cell yield per mole of ATP (Y(ATP)) (about 7.0 g of cells per mol of ATP) compared to 10.9 g of cells per mol of ATP for the wild type. The lower Y(ATP) suggests the operation of a futile energy cycle in the absence of PFL in this strain. An understanding of the controls at the pyruvate node during anaerobic growth is expected to provide unique insights into rational metabolic engineering of E. coli and related bacteria for the production of various biobased products at high rates and yields.

TonB-dependent transporters: regulation, structure, and function

TonB-dependent transporters (TBDTs) are bacterial outer membrane proteins that bind and transport ferric chelates, called siderophores, as well as vitamin B(12), nickel complexes, and carbohydrates. The transport process requires energy in the form of proton motive force and a complex of three inner membrane proteins, TonB-ExbB-ExbD, to transduce this energy to the outer membrane. The siderophore substrates range in complexity from simple small molecules such as citrate to large proteins such as serum transferrin and hemoglobin. Because iron uptake is vital for almost all bacteria, expression of TBDTs is regulated in a number of ways that include metal-dependent regulators, σ/anti-σ factor systems, small RNAs, and even a riboswitch. In recent years, many new structures of TBDTs have been solved in various states, resulting in a more complete understanding of siderophore selectivity and binding, signal transduction across the outer membrane, and interaction with the TonB-ExbB-ExbD complex. However, the transport mechanism is still unclear. In this review, we summarize recent progress in understanding regulation, structure, and function in TBDTs and questions remaining to be answered.

An insight into the role of phosphotransacetylase (pta) and the acetate/acetyl-CoA node in Escherichia coli

Background

Acetate metabolism in Escherichia coli plays an important role in the control of the central metabolism and in bioprocess performance. The main problems related to the use of E. coli as cellular factory are i) the deficient utilization of carbon source due to the excretion of acetate during aerobic growth, ii) the inhibition of cellular growth and protein production by acetate and iii) the need for cofactor recycling (namely redox coenzymes and free CoASH) to sustain balanced growth and cellular homeostasis.

Results

This work analyzes the effect of mutations in the acetate excretion/assimilation pathways, acetyl-CoA synthethase (acs) and phosphotransacetylase (pta), in E. coli BW25113 grown on glucose or acetate minimal media. Biomass and metabolite production, redox (NADH/NAD⁺) and energy (ATP) state, enzyme activities and gene expression profiles related to the central metabolism were analyzed. The knock-out of pta led to a more altered phenotype than that of acs. Deletion of pta reduced the ability to grow on acetate as carbon source and strongly affected the expression of several genes related to central metabolic pathways.

Conclusion

Results showed that pta limits biomass yield in aerobic glucose cultures, due to acetate production (overflow metabolism) and its inefficient use during glucose starvation. Deletion of pta severely impaired growth on acetate minimal medium and under anaerobiosis due to decreased acetyl-coenzyme A synthethase, glyoxylate shunt and gluconeogenic activities, leading to lower growth rate. When acetate is used as carbon source, the joint expression of pta and acs is crucial for growth and substrate assimilation, while pta deletion severely impaired anaerobic growth. Finally, at an adaptive level, pta deficiency makes the strain more sensitive to environmental changes and de-regulates the central metabolism.

Dual role of LldR in regulation of the lldPRD operon, involved in L-lactate metabolism in Escherichia coli

The lldPRD operon of Escherichia coli, involved in L-lactate metabolism, is induced by growth in this compound. We experimentally identified that this system is transcribed from a single promoter with an initiation site located 110 nucleotides upstream of the ATG start codon. On the basis of computational data, it had been proposed that LldR and its homologue PdhR act as regulators of the lldPRD operon. Nevertheless, no experimental data on the function of these regulators have been reported so far. Here we show that induction of an lldP-lacZ fusion by L-lactate is lost in an Delta lldR mutant, indicating the role of LldR in this induction. Expression analysis of this construct in a pdhR mutant ruled out the participation of PdhR in the control of lldPRD. Gel shift experiments showed that LldR binds to two operator sites, O1 (positions -105 to -89) and O2 (positions +22 to +38), with O1 being filled at a lower concentration of LldR. L-Lactate induced a conformational change in LldR that did not modify its DNA binding activity. Mutations in O1 and O2 enhanced the basal transcriptional level. However, only mutations in O1 abolished induction by L-lactate. Mutants with a change in helical phasing between O1 and O2 behaved like O2 mutants. These results were consistent with the hypothesis that LldR has a dual role, acting as a repressor or an activator of lldPRD. We propose that in the absence of L-lactate, LldR binds to both O1 and O2, probably leading to DNA looping and the repression of transcription. Binding of L-lactate to LldR promotes a conformational change that may disrupt the DNA loop, allowing the formation of the transcription open complex.

Regulation of L-lactate utilization by the FadR-type regulator LldR of Corynebacterium glutamicum

Corynebacterium glutamicum can grow on L-lactate as a sole carbon and energy source. The NCgl2816-lldD operon encoding a putative transporter (NCgl2816) and a quinone-dependent L-lactate dehydrogenase (LldD) is required for L-lactate utilization. DNA affinity chromatography revealed that the FadR-type regulator LldR (encoded by NCgl2814) binds to the upstream region of NCgl2816-lldD. Overexpression of lldR resulted in strongly reduced NCgl2816-lldD mRNA levels and strongly reduced LldD activity, and as a consequence, a severe growth defect was observed in cells grown on L-lactate as the sole carbon and energy source, but not in cells grown on glucose, ribose, or acetate. Deletion of lldR had no effect on growth on these carbon sources but resulted in high NCgl2816-lldD mRNA levels and high LldD activity in the presence and absence of L-lactate. Purified His-tagged LldR bound to a 54-bp fragment of the NCgl2816-lldD promoter, which overlaps with the transcriptional start site determined by random amplification of cDNA ends-PCR and contains a putative operator motif typical of FadR-type regulators, which is -1TNGTNNNACNA10. Mutational analysis revealed that this motif with hyphenated dyad symmetry is essential for binding of LldD to the NCgl2816-lldD promoter. L-Lactate, but not D-lactate, interfered with binding of LldRHis to the NCgl2816-lldD promoter. Thus, during growth on media lacking L-lactate, LldR represses expression of NCgl2816-lldD. In the presence of L-lactate in the growth medium or under conditions leading to intracellular L-lactate accumulation, the L-lactate utilization operon is induced.

Rethinking proteasome evolution: two novel bacterial proteasomes

The proteasome is a multisubunit structure that degrades proteins. Protein degradation is an essential component of regulation because proteins can become misfolded, damaged, or unnecessary. Proteasomes and their homologues vary greatly in complexity: from HslV (heat shock locus v), which is encoded by 1 gene in bacteria, to the eukaryotic 20S proteasome, which is encoded by more than 14 genes. Despite this variation in complexity, all the proteasomes are composed of homologous subunits. We searched 238 complete bacterial genomes for structures related to the proteasome and found evidence of two novel groups of bacterial proteasomes. The first, which we name Anbu, is sparsely distributed among cyanobacteria and proteobacteria. We hypothesize that Anbu must be very ancient because of its distribution within the cyanobacteria, and that it has been lost in many more recent species. We also present evidence for a fourth type of bacterial proteasome found in a few beta-proteobacteria, which we call beta-proteobacteria proteasome homologue (BPH). Sequence and structural analyses show that Anbu and BPH are both distinct from known bacterial proteasomes but have homologous structures. Anbu is encoded by one gene, so we postulate a duplication of Anbu created the 20S proteasome. Anbu's function appears to be related to transglutaminase activity, not the general stress response associated with HslV. We have found different combinations of Anbu, BPH, and HslV within these bacterial genomes, which raises questions about specialized protein degradation systems.

GntR family of regulators in Mycobacterium smegmatis: a sequence and structure based characterization

Background

Mycobacterium smegmatis is fast growing non-pathogenic mycobacteria. This organism has been widely used as a model organism to study the biology of other virulent and extremely slow growing species like Mycobacterium tuberculosis. Based on the homology of the N-terminal DNA binding domain, the recently sequenced genome of M. smegmatis has been shown to possess several putative GntR regulators. A striking characteristic feature of this family of regulators is that they possess a conserved N-terminal DNA binding domain and a diverse C-terminal domain involved in the effector binding and/or oligomerization. Since the physiological role of these regulators is critically dependent upon effector binding and operator sites, we have analysed and classified these regulators into their specific subfamilies and identified their potential binding sites.

Results

The sequence analysis of M. smegmatis putative GntRs has revealed that FadR, HutC, MocR and the YtrA-like regulators are encoded by 45, 8, 8 and 1 genes respectively. Further out of 45 FadR-like regulators, 19 were classified into the FadR group and 26 into the VanR group. All these proteins showed similar secondary structural elements specific to their respective subfamilies except MSMEG_3959, which showed additional secondary structural elements. Using the reciprocal BLAST searches, we further identified the orthologs of these regulators in Bacillus subtilis and other mycobacteria. Since the expression of many regulators is auto-regulatory, we have identified potential operator sites for a number of these GntR regulators by analyzing the upstream sequences.

Conclusion

This study helps in extending the annotation of M. smegmatis GntR proteins. It identifies the GntR regulators of M. smegmatis that could serve as a model for studying orthologous regulators from virulent as well as other saprophytic mycobacteria. This study also sheds some light on the nucleotide preferences in the target-motifs of GntRs thus providing important leads for initiating the experimental characterization of these proteins, construction of the gene regulatory network for these regulators and an understanding of the influence of these proteins on the physiology of the mycobacteria.

PdhR (pyruvate dehydrogenase complex regulator) controls the respiratory electron transport system in Escherichia coli

The pyruvate dehydrogenase (PDH) multienzyme complex plays a key role in the metabolic interconnection between glycolysis and the citric acid cycle. Transcription of the Escherichia coli genes for all three components of the PDH complex in the pdhR-aceEF-lpdA operon is repressed by the pyruvate-sensing PdhR, a GntR family transcription regulator, and derepressed by pyruvate. After a systematic search for the regulation targets of PdhR using genomic systematic evolution of ligands by exponential enrichment (SELEX), we have identified two novel targets, ndh, encoding NADH dehydrogenase II, and cyoABCDE, encoding the cytochrome bo-type oxidase, both together forming the pathway of respiratory electron transport downstream from the PDH cycle. PDH generates NADH, while Ndh and CyoABCDE together transport electrons from NADH to oxygen. Using gel shift and DNase I footprinting assays, the PdhR-binding site (PdhR box) was defined, which includes a palindromic consensus sequence, ATTGGTNNNACCAAT. The binding in vitro of PdhR to the PdhR box decreased in the presence of pyruvate. Promoter assays in vivo using a two-fluorescent-protein vector also indicated that the newly identified operons are repressed by PdhR and derepressed by the addition of pyruvate. Taken together, we propose that PdhR is a master regulator for controlling the formation of not only the PDH complex but also the respiratory electron transport system.

In silico analysis and characterization of GntR family of regulators from Mycobacterium tuberculosis

The genome of Mycobacterium tuberculosis contains a large number of hypothetical and poorly characterized proteins including the proteins belonging to the GntR family. The regulators of this family show a conserved N-terminal DNA-binding domain but have a highly diverse C-terminal domain involved in the effector-binding and/or oligomerization. This heterogeneity has led to a further classification of this family into various subfamilies. The sequence analysis of the M. tuberculosis genome revealed that five genes encode for FadR-like regulators, one gene for HutC-like regulator and one for YtrA-like regulator. This classification was also consistent with specific secondary structural features known to be associated with FadR, HutC and YtrA subfamilies. Out of the five FadR-like regulators three of the regulators were further subclassified into FadR group and two of them into the VanR group. Interestingly Rv3060c, a FadR-like regulator, was shown to have an unusual size which led us to demonstrate it as a product of a gene duplication and fusion event. Thus this study extends the genome annotation of M. tuberculosis and provides important leads for initiating experimental characterization of these proteins, which in turn will enrich our knowledge of their role in cellular physiology.

Function, attachment and synthesis of lipoic acid in Escherichia coli

A series of genetic, biochemical, and physiological studies in Escherichia coli have elucidated the unusual pathway whereby lipoic acid is synthesized. Here we describe the results of these investigations as well as the functions of enzyme proteins that are modified by covalent attachment of lipoic acid and the enzymes that catalyze the modification reactions. Some aspects of the synthesis and attachment mechanisms have strong parallels in the pathways used in synthesis and attachment of biotin and these are compared and contrasted. Homologues of the lipoic acid metabolism proteins are found in all branches of life, save the Archea, and thus these findings seem to have wide biological relevance.

Functional analysis of 1440 Escherichia coli genes using the combination of knock-out library and phenotype microarrays

Escherichia coli is one of the best elucidated organisms. However, about 40% of E. coli genes have not been assigned to their function yet. We analyzed 1440 single gene knock-out mutants using the GN2-MicroPlate, which permits assay of 95 carbon-source utilizations simultaneously. In the knock-out library there are 1044 of so called y-genes with no apparent function. The raw dataset was analyzed and genes were interrelated by the clustering method of the GeneSpring software. In the resulted dendrogram of genes, a group of genes with known and related function tended to be assembled into a cluster. Our clustering method would be useful for functional assignment of so called y-genes with no apparent function, since the resulted dendrogram could connect y-genes to phenotype and function of well-studied genes.

E1 enzyme of the pyruvate dehydrogenase complex in Corynebacterium glutamicum: molecular analysis of the gene and phylogenetic aspects

The E1p enzyme is an essential part of the pyruvate dehydrogenase complex (PDHC) and catalyzes the oxidative decarboxylation of pyruvate with concomitant acetylation of the E2p enzyme within the complex. We analyzed the Corynebacterium glutamicum aceE gene, encoding the E1p enzyme, and constructed and characterized an E1p-deficient mutant. Sequence analysis of the C. glutamicum aceE gene and adjacent regions revealed that aceE is not flanked by genes encoding other enzymes of the PDHC. Transcriptional analysis revealed that aceE from C. glutamicum is monocistronic and that its transcription is initiated 121 nucleotides upstream of the translational start site. Inactivation of the chromosomal aceE gene led to the inability to grow on glucose and to the absence of PDHC and E1p activities, indicating that only a single E1p enzyme is present in C. glutamicum and that the PDHC is essential for the growth of this organism on carbohydrate substrates. Surprisingly, the E1p enzyme of C. glutamicum showed up to 51% identity to homodimeric E1p proteins from gram-negative bacteria but no similarity to E1 alpha- or beta-subunits of heterotetrameric E1p enzymes which are generally assumed to be typical for gram-positives. To investigate the distribution of E1p enzymes in bacteria, we compiled and analyzed the phylogeny of 46 homodimeric E1p proteins and of 58 alpha-subunits of heterotetrameric E1p proteins deposited in public databases. The results revealed that the distribution of homodimeric and heterotetrameric E1p subunits in bacteria is not in accordance with the rRNA-based phylogeny of bacteria and is more heterogeneous than previously assumed.

Characterization of a Corynebacterium glutamicum lactate utilization operon induced during temperature-triggered glutamate production

Gene expression changes of glutamate-producing Corynebacterium glutamicum were identified in transcriptome comparisons by DNA microarray analysis. During glutamate production induced by a temperature shift, C. glutamicum strain 2262 showed significantly higher mRNA levels of the NCgl2816 and NCgl2817 genes than its non-glutamate-producing derivative 2262NP. Reverse transcription-PCR analysis showed that the two genes together constitute an operon. NCgl2816 putatively codes for a lactate permease, while NCgl2817 was demonstrated to encode quinone-dependent l-lactate dehydrogenase, which was named LldD. C. glutamicum LldD displayed Michaelis-Menten kinetics for the substrate l-lactate with a K(m) of about 0.51 mM. The specific activity of LldD was about 10-fold higher during growth on l-lactate or on an l-lactate-glucose mixture than during growth on glucose, d-lactate, or pyruvate, while the specific activity of quinone-dependent d-lactate dehydrogenase differed little with the carbon source. RNA levels of NCgl2816 and lldD were about 18-fold higher during growth on l-lactate than on pyruvate. Disruption of the NCgl2816-lldD operon resulted in loss of the ability to utilize l-lactate as the sole carbon source. Expression of lldD restored l-lactate utilization, indicating that the function of the permease gene NCgl2816 is dispensable, while LldD is essential, for growth of C. glutamicum on l-lactate.

Integrated regulatory responses of fimB to N-acetylneuraminic (sialic) acid and GlcNAc in Escherichia coli K-12

Bacterial-host attachment by means of bacterial adhesins is a key step in host colonization. Phase variation (reversible on-off switching) of the type 1 fimbrial adhesin of Escherichia coli involves a DNA inversion catalyzed by FimB (switching in either direction) or FimE (mainly on-to-off switching). fimB is separated from the divergent yjhATS operon by a large (1.4 kbp) intergenic region. Short ( approximately 28 bp) cis-active elements (regions 1 and 2) close to yjhA stimulate fimB expression and are required for sialic acid (Neu(5)Ac) sensitivity of its expression [El-Labany, S., Sohanpal, B. K., Lahooti, M., Akerman, R. & Blomfield, I. C. (2003) Mol. Microbiol. 49, 1109-1118]. Here, we show that whereas NanR, a sialic acid-response regulator, binds to region 1, NagC, a GlcNAc-6P-responsive protein, binds to region 2 instead. The NanR- and NagC-binding sites lie adjacent to deoxyadenosine methylase (Dam) methylation sites (5'-GATC) that are protected from modification, and the two regulators are shown to be required for methylation protection at regions 1 and 2, respectively. Mutations in nanR and nagC diminish fimB expression, and both fimB expression and FimB recombination are inhibited by GlcNAc (3- and >35-fold, respectively). Sialic acid catabolism generates GlcNAc-6-P, and whereas GlcNAc disrupts methylation protection by NagC alone, Neu(5)Ac inhibits the protection mediated by both NanR and NagC as expected. Type 1 fimbriae are proinflammatory, and host defenses enhance the release of both Neu(5)Ac and GlcNAc by a variety of mechanisms. Inhibition of type 1 fimbriation by these amino sugars may thus help balance the interaction between E. coli and its hosts.

Bacillus subtilis GabR, a protein with DNA-binding and aminotransferase domains, is a PLP-dependent transcriptional regulator

Bacillus subtilis GabR is a member of a poorly characterized but widespread family of chimeric bacterial proteins that have apparent DNA binding and aminotransferase domains. GabR positively regulates expression of the gabTD operon responsible for utilization of gamma-aminobutyric acid (GABA) and represses the divergently transcribed gabR gene. Purified GabR bound specifically to the DNA region overlapping the -35 region of the gabT promoter and the -10 and +1 regions of the gabR promoter. Two 6 bp direct repeats located at the ends of this region appeared to be essential for GabR binding. In transcription reactions in vitro, GabR alone repressed expression from the gabR promoter but activated expression from the gabT promoter only in the presence of GABA and pyridoxal 5'-phosphate, an essential cofactor of aminotransferases. A similar requirement for pyridoxal 5'-phosphate and GABA for GabR-mediated transcription activation was shown in vivo. In vitro this requirement could be partially satisfied with pyridoxamine 5'-phosphate and succinic semialdehyde, the products of a GABA-dependent aminotransferase half-reaction. We hypothesize that the GabR-catalyzed aminotransferase-like reaction between GABA and pyridoxal 5'-phosphate is essential for GabR action as a transcriptional activator.