Metabolic flux responses to pyruvate kinase knockout in Escherichia coli

A parallel glycolysis provides a selective advantage through rapid growth acceleration

Glycolysis is a universal metabolic process that breaks down glucose to produce adenosine triphosphate (ATP) and biomass precursors. The Entner-Doudoroff (ED) pathway is a glycolytic pathway that parallels textbook glycolysis but yields half as much ATP. Accordingly, in organisms that possess both glycolytic pathways (for example, Escherichia coli), its raison d'être remains a mystery. In this study, we found that the ED pathway provides a selective advantage during growth acceleration. Upon carbon and nitrogen upshifts, E. coli accelerates growth faster with than without the ED pathway. Concurrent isotope tracing reveals that the ED pathway flux increases faster than that of textbook glycolysis. We attribute the fast response time of the ED pathway to its strong thermodynamic driving force and streamlining of glucose import. Intermittent nutrient supply manifests the evolutionary advantage of the parallel glycolysis; thus, the dynamic nature of an ostensibly redundant pathway's role in promoting rapid adaptation constitutes a metabolic design principle.

A reservoir of timescales emerges in recurrent circuits with heterogeneous neural assemblies

The temporal activity of many physical and biological systems, from complex networks to neural circuits, exhibits fluctuations simultaneously varying over a large range of timescales. Long-tailed distributions of intrinsic timescales have been observed across neurons simultaneously recorded within the same cortical circuit. The mechanisms leading to this striking temporal heterogeneity are yet unknown. Here, we show that neural circuits, endowed with heterogeneous neural assemblies of different sizes, naturally generate multiple timescales of activity spanning several orders of magnitude. We develop an analytical theory using rate networks, supported by simulations of spiking networks with cell-type specific connectivity, to explain how neural timescales depend on assembly size and show that our model can naturally explain the long-tailed timescale distribution observed in the awake primate cortex. When driving recurrent networks of heterogeneous neural assemblies by a time-dependent broadband input, we found that large and small assemblies preferentially entrain slow and fast spectral components of the input, respectively. Our results suggest that heterogeneous assemblies can provide a biologically plausible mechanism for neural circuits to demix complex temporal input signals by transforming temporal into spatial neural codes via frequency-selective neural assemblies.

Evaluating scaling of capillary photo-biofilm reactors for high cell density cultivation of mixed trophies artificial microbial consortia

Capillary biofilm reactors (CBRs) are attractive for growing photoautotrophic bacteria as they allow high cell-density cultivation. Here, we evaluated the CBR system's suitability to grow an artificial consortium composed of Synechocystis sp. PCC 6803 and Pseudomonas sp. VBL120. The impact of reactor material, flow rate, pH, O2, and medium composition on biomass development and long-term biofilm stability at different reactor scales was studied. Silicone was superior over other materials like glass or PVC due to its excellent O2 permeability. High flow rates of 520 μL min-1 prevented biofilm sloughing in 1 m capillary reactors, leading to a 54% higher biomass dry weight combined with the lowest O2 concentration inside the reactor compared to standard operating conditions. Further increase in reactor length to 5 m revealed a limitation in trace elements. Increasing trace elements by a factor of five allowed for complete surface coverage with a biomass dry weight of 36.8 g m-2 and, thus, a successful CBR scale-up by a factor of 25. Practical application: Cyanobacteria use light energy to upgrade CO2, thereby holding the potential for carbon-neutral production processes. One of the persisting challenges is low cell density due to light limitations and O2 accumulation often occurring in established flat panel or tubular photobioreactors. Compared to planktonic cultures, much higher cell densities (factor 10 to 100) can be obtained in cyanobacterial biofilms. The capillary biofilm reactor (CBR) offers good growth conditions for cyanobacterial biofilms, but its applicability has been shown only on the laboratory scale. Here, a first scale-up study based on sizing up was performed, testing the feasibility of this system for large-scale applications. We demonstrate that by optimizing nutrient supply and flow conditions, the system could be enlarged by factor 25 by enhancing the length of the reactor. This reactor concept, combined with cyanobacterial biofilms and numbering up, holds the potential to be applied as a flexible, carbon-neutral production platform for value-added compounds.

Mitigating biomass composition uncertainties in flux balance analysis using ensemble representations

The biomass equation is a critical component in genome-scale metabolic models (GEMs): it is used as the de facto objective function in flux balance analysis (FBA). This equation accounts for the quantities of all known biomass precursors that are required for cell growth based on the macromolecular and monomer compositions measured at certain conditions. However, it is often reported that the macromolecular composition of cells could change across different environmental conditions and thus the use of the same single biomass equation in FBA, under multiple conditions, is questionable. Herein, we first investigated the qualitative and quantitative variations of macromolecular compositions of three representative host organisms, Escherichia coli, Saccharomyces cerevisiae and Cricetulus griseus, across different environmental/genetic variations. While macromolecular building blocks such as RNA, protein, and lipid composition vary notably, changes in fundamental biomass monomer units such as nucleotides and amino acids are not appreciable. We also observed that flux predictions through FBA is quite sensitive to macromolecular compositions but not the monomer compositions. Based on these observations, we propose ensemble representations of biomass equation in FBA to account for the natural variation of cellular constituents. Such ensemble representations of biomass better predicted the flux through anabolic reactions as it allows for the flexibility in the biosynthetic demands of the cells. The current study clearly highlights that certain component of the biomass equation indeed vary across different conditions, and the ensemble representation of biomass equation in FBA by accounting for such natural variations could avoid inaccuracies that may arise from in silico simulations.

Metagenomic Sequencing Identified Specific Bacteriophage Signature Discriminating between Healthy and Diarrheal Neonatal Piglets

Neonatal diarrhea is one of the most severe diseases in human beings and pigs, leading to high mortality and growth faltering. Gut microbiome-related studies mostly focus on the relationship between bacteria and neonatal diarrhea onset, and no research study has investigated the role of the gut virome in neonatal diarrhea. Here, using metagenomic sequencing, we characterized the fecal viral community of diarrheal and healthy neonatal piglets. We found that the viral community of diarrheal piglets showed higher individual heterogeneity and elevated abundance of Myoviridae. By predicting the bacterial host of the identified viral genomes, phages infecting Proteobacteria, especially E. coli, were the dominant taxa in neonatal diarrheal piglets. Consistent with this, the antibiotic resistance gene of E. coli origin was also enriched in neonatal diarrheal piglets. Finally, we established a random forest model to accurately discriminate between neonatal diarrheal piglets and healthy controls and identified genus E. coli- and genus listeria-infecting bacteriophages, including psa and C5 viruses, as key biomarkers. In conclusion, we provide the first glance of viral community and function characteristics in diarrheal and healthy neonatal piglets. These findings expand our understanding of the relationship among phages, bacteria and diarrhea, and may facilitate the development of therapeutics for the prevention and treatment of neonatal diarrhea.

Neural mechanisms underlying the temporal organization of naturalistic animal behavior

Naturalistic animal behavior exhibits a strikingly complex organization in the temporal domain, with variability arising from at least three sources: hierarchical, contextual, and stochastic. What neural mechanisms and computational principles underlie such intricate temporal features? In this review, we provide a critical assessment of the existing behavioral and neurophysiological evidence for these sources of temporal variability in naturalistic behavior. Recent research converges on an emergent mechanistic theory of temporal variability based on attractor neural networks and metastable dynamics, arising via coordinated interactions between mesoscopic neural circuits. We highlight the crucial role played by structural heterogeneities as well as noise from mesoscopic feedback loops in regulating flexible behavior. We assess the shortcomings and missing links in the current theoretical and experimental literature and propose new directions of investigation to fill these gaps.

Heterogeneity of soil bacterial and bacteriophage communities in three rice agroecosystems and potential impacts of bacteriophage on nutrient cycling

BACKGROUND

As genetic entities infecting and replicating only in bacteria, bacteriophages can regulate the community structure and functions of their host bacteria. The ecological roles of bacteriophages in aquatic and forest environments have been widely explored, but those in agroecosystems remains limited. Here, we used metagenomic sequencing to analyze the diversity and interactions of bacteriophages and their host bacteria in soils from three typical rice agroecosystems in China: double cropping in Guangzhou, southern China, rice-wheat rotation cropping in Nanjing, eastern China and early maturing single cropping in Jiamusi, northeastern China. Enterobacter phage-NJ was isolated and its functions on soil nitrogen cycling and effect on soil bacterial community structure were verified in pot inoculation experiments and 16S rRNA gene sequencing.

RESULTS

Soil bacterial and viral diversity and predicted functions varied among the three agroecosystems. Genes detected in communities from the three agroecosystems were associated with typical functions: soil bacteria in Jiamusi were significantly enriched in genes related to carbohydrate metabolism, in Nanjing with xenobiotics biodegradation and metabolism, and in Guangzhou with virulence factors and scarce in secondary metabolite biosynthesis, which might lead to a significant occurrence of rice bacterial diseases. The virus community structure varies significantly among the three ecosystems, only 13.39% of the total viral species were shared by the three rice agroecosystems, 59.56% of the viral species were specific to one agroecosystem. Notably, over-represented auxiliary carbohydrate-active enzyme (CAZyme) genes were identified in the viruses, which might assist host bacteria in metabolizing carbon, and 67.43% of these genes were present in Jiamusi. In bacteriophage isolation and inoculation experiments, Enterobacter bacteriophage-NJ reduced the nitrogen fixation capacity of soil by lysing N-fixing host bacteria and changed the soil bacterial diversity and community structure.

CONCLUSION

Our results showed that diversity and function predicted of paddy soil bacteria and viruses varied in the three agroecosystems. Soil bacteriophages can affect nutrient cycling by boosting host metabolism through the carried auxiliary metabolic genes (AMGs) and lysing the host bacteria that are involved in biogeochemical cycles. These findings form a basis for better understanding bacterial and bacteriophage diversity in different rice agroecosystems, laying a solid foundation for further studies of soil microbial communities that support ecofriendly production of healthy rice.

Metabolic and process engineering for microbial production of protocatechuate with Corynebacterium glutamicum

3,4-Dihydroxybenzoate (protocatechuate, PCA) is a phenolic compound naturally found in edible vegetables and medicinal herbs. PCA is of high interest in the chemical industry and has wide potential for pharmaceutical applications. We designed and constructed a novel Corynebacterium glutamicum strain to enable the efficient utilization of d-xylose for microbial production of PCA. Shake flask cultivation of the engineered strain showed a maximum PCA titer of 62.1 ± 12.1 mM (9.6 ± 1.9 g L^-1 ) from d-xylose as the primary carbon and energy source. The corresponding yield was 0.33 C-mol PCA per C-mol d-xylose, which corresponds to 38% of the maximum theoretical yield. Under growth-decoupled bioreactor conditions, a comparable PCA titer and a total amount of 16.5 ± 1.1 g PCA could be achieved when d-glucose and d-xylose were combined as orthogonal carbon substrates for biocatalyst provision and product synthesis, respectively. Downstream processing of PCA was realized via electrochemically induced crystallization by taking advantage of the pH-dependent properties of PCA. This resulted in a maximum final purity of 95.4%. The established PCA production process represents a highly sustainable approach, which will serve as a blueprint for the bio-based production of other hydroxybenzoic acids from alternative sugar feedstocks.

Use of an optical sensor to directly monitor the metabolic activity of growing bacteria

The metabolic activity of growing bacteria was directly monitored by using an electro-optical (EO) sensor. The sensor enables examination of bacteria in batch and continuous cultures. As examples, we report studies with Еscherichia coli, a bacterium with an aerobic type of metabolism, and Lactobacillus plantarum, a bacterium with an anaerobic type of metabolism. Bacterial growth was accompanied by a simultaneous change in both the hydrodynamic mean size (HMS) of the bacteria and the concentration of ions in the cytoplasm (CIC). Both variables were associated with the regulation of cellular metabolic activity, which can be cyclic during intense bacterial growth. A simultaneous change in metabolic activity and osmotic regulation was also found. For СIC and HMS measurements, we used online results of the EO analysis of cells suspended in water. The measured results for the CIC and HMS can be used to directly monitor bacterial metabolism. The results of this study are of practical importance for the real-time EO monitoring of the metabolic activity of growing bacteria without preliminary sample preparation.

Metabolic flux configuration determination using information entropy

Constraint-based models use steady-state mass balances to define a solution space of flux configurations, which can be narrowed down by measuring as many fluxes as possible. Due to loops and redundant pathways, this process typically yields multiple alternative solutions. To address this ambiguity, flux sampling can estimate the probability distribution of each flux, or a flux configuration can be singled out by further minimizing the sum of fluxes according to the assumption that cellular metabolism favors states where enzyme-related costs are economized. However, flux sampling is susceptible to artifacts introduced by thermodynamically infeasible cycles and is it not clear if the economy of fluxes assumption (EFA) is universally valid. Here, we formulated a constraint-based approach, MaxEnt, based on the principle of maximum entropy, which in this context states that if more than one flux configuration is consistent with a set of experimentally measured fluxes, then the one with the minimum amount of unwarranted assumptions corresponds to the best estimation of the non-observed fluxes. We compared MaxEnt predictions to Escherichia coli and Saccharomyces cerevisiae publicly available flux data. We found that the mean square error (MSE) between experimental and predicted fluxes by MaxEnt and EFA-based methods are three orders of magnitude lower than the median of 1,350,000 MSE values obtained using flux sampling. However, only MaxEnt and flux sampling correctly predicted flux through E. coli’s glyoxylate cycle, whereas EFA-based methods, in general, predict no flux cycles. We also tested MaxEnt predictions at increasing levels of overflow metabolism. We found that MaxEnt accuracy is not affected by overflow metabolism levels, whereas the EFA-based methods show a decreasing performance. These results suggest that MaxEnt is less sensitive than flux sampling to artifacts introduced by thermodynamically infeasible cycles and that its predictions are less susceptible to overfitting than EFA-based methods.

Escherichia coli metabolism under short-term repetitive substrate dynamics: adaptation and trade-offs

Background

Microbial metabolism is highly dependent on the environmental conditions. Especially, the substrate concentration, as well as oxygen availability, determine the metabolic rates. In large-scale bioreactors, microorganisms encounter dynamic conditions in substrate and oxygen availability (mixing limitations), which influence their metabolism and subsequently their physiology. Earlier, single substrate pulse experiments were not able to explain the observed physiological changes generated under large-scale industrial fermentation conditions.

Results

In this study we applied a repetitive feast–famine regime in an aerobic Escherichia coli culture in a time-scale of seconds. The regime was applied for several generations, allowing cells to adapt to the (repetitive) dynamic environment. The observed response was highly reproducible over the cycles, indicating that cells were indeed fully adapted to the regime. We observed an increase of the specific substrate and oxygen consumption (average) rates during the feast–famine regime, compared to a steady-state (chemostat) reference environment. The increased rates at same (average) growth rate led to a reduced biomass yield (30% lower). Interestingly, this drop was not followed by increased by-product formation, pointing to the existence of energy-spilling reactions. During the feast–famine cycle, the cells rapidly increased their uptake rate. Within 10 s after the beginning of the feeding, the substrate uptake rate was higher (4.68 μmol/g_CDW/s) than reported during batch growth (3.3 μmol/g_CDW/s). The high uptake led to an accumulation of several intracellular metabolites, during the feast phase, accounting for up to 34% of the carbon supplied. Although the metabolite concentrations changed rapidly, the cellular energy charge remained unaffected, suggesting well-controlled balance between ATP producing and ATP consuming reactions.

Conclusions

The adaptation of the physiology and metabolism of E. coli under substrate dynamics, representative for large-scale fermenters, revealed the existence of several cellular mechanisms coping with stress. Changes in the substrate uptake system, storage potential and energy-spilling processes resulted to be of great importance. These metabolic strategies consist a meaningful step to further tackle reduced microbial performance, observed under large-scale cultivations.

Applications of Chromatography-Ultra High-Resolution MS for Stable Isotope-Resolved Metabolomics (SIRM) Reconstruction of Metabolic Networks

Metabolism is a complex network of compartmentalized and coupled chemical reactions, which often involve transfers of substructures of biomolecules, thus requiring metabolite substructures to be tracked. Stable isotope resolved metabolomics (SIRM) enables pathways reconstruction, even among chemically identical metabolites, by tracking the provenance of stable isotope-labeled substructures using NMR and ultrahigh resolution (UHR) MS. The latter can resolve and count isotopic labels in metabolites and can identify isotopic enrichment in substructures when operated in tandem MS mode. However, MS2 is difficult to implement with chromatography-based UHR-MS due to lengthy MS1 acquisition time that is required to obtain the molecular isotopologue count, which is further exacerbated by the numerous isotopologue source ions to fragment. We review here recent developments in tandem MS applications of SIRM to obtain more detailed information about isotopologue distributions in metabolites and their substructures.

Microbial energy management-A product of three broad tradeoffs

Wherever thermodynamics allows, microbial life has evolved to transform and harness energy. Microbial life thus abounds in the most unexpected places, enabled by profound metabolic diversity. Within this diversity, energy is transformed primarily through variations on a few core mechanisms. Energy is further managed by the physiological processes of cell growth and maintenance that use energy. Some aspects of microbial physiology are streamlined for energetic efficiency while other aspects seem suboptimal or even wasteful. We propose that the energy that a microbe harnesses and devotes to growth and maintenance is a product of three broad tradeoffs: (i) economic, trading enzyme synthesis or operational cost for functional benefit, (ii) environmental, trading optimization for a single environment for adaptability to multiple environments, and (iii) thermodynamic, trading energetic yield for forward metabolic flux. Consideration of these tradeoffs allows one to reconcile features of microbial physiology that seem to opposingly promote either energetic efficiency or waste.

Central metabolic pathway modification to improve L-tryptophan production in Escherichia coli

Tryptophan, an aromatic amino acid, has been widely used in food industry because it participates in the regulation of protein synthesis and metabolic network in vivo. In this study, we obtained a strain named TRP03 by enhancing the tryptophan synthesis pathway, which could accumulate tryptophan at approximately 35 g/L in a 5 L bioreactor. We then modified the central metabolic pathway of TRP03, to increase the supply of the precursor phosphoenolpyruvate (PEP), the genes related to PEP were modified. Furthermore, citric acid transport system and TCA were upregulated to effectively increase cell growth. We observed that strain TRP07 that could accumulate tryptophan at approximately 49 g/L with a yield of 0.186 g tryptophan/g glucose in a 5 L bioreactor. By-products such as glutamate and acetic acid were reduced to 0.8 g/L and 2.2 g/L, respectively.

The Absence of Pyruvate Kinase Affects Glucose-Dependent Carbon Catabolite Repression in Bacillus subtilis

Pyruvate is a key intermediate of diverse metabolic pathways of central carbon metabolism. In addition to being the end product of glycolysis, pyruvate is an essential carbon distribution point to oxidative metabolism, amino acid and fatty acid syntheses, and overflow metabolite production. Hence, a tight regulation of pyruvate kinase (Pyk) activity is of great importance. This study aimed to analyze targeted metabolites from several pathways and possible changes in Bacillus subtilis lacking Pyk. Wild type and Δpyk cells were cultivated in chemically defined medium with glucose and pyruvate as carbon sources, and the extracted metabolites were analyzed by ¹H-NMR, GC-MS, HPLC-MS, and LC-MS/MS. The results showed that the perturbation created in the pyruvate node drove an adaptation to new conditions by altering the nutritional compounds’ consumption. In Δpyk, pyruvate, which is subject to glucose-dependent carbon catabolite repression, did not comply with the hierarchy in carbon source utilization. Other metabolic alterations were observed such as the higher secretion of the overflow metabolites acetoin and 2,3-butanediol by Δpyk. Our results help to elucidate the regulatory transport of glucose and pyruvate in B. subtilis and possible metabolic reroute to alternative pathways in the absence of Pyk.

Redox rebalance against genetic perturbations and modulation of central carbon metabolism by the oxidative stress regulation

The micro-aerophilic organisms and aerobes as well as yeast and higher organisms have evolved to gain energy through respiration (via oxidative phosphorylation), thereby enabling them to grow much faster than anaerobes. However, during respiration, reactive oxygen species (ROSs) are inherently (inevitably) generated, and threaten the cell's survival. Therefore, living organisms (or cells) must furnish the potent defense systems to keep such ROSs at harmless level, where the cofactor balance plays crucial roles. Namely, NADH is the source of energy generation (catabolism) in the respiratory chain reactions, through which ROSs are generated, while NADPH plays important roles not only for the cell synthesis (anabolism) but also for detoxifying ROSs. Therefore, the cell must rebalance the redox ratio by modulating the fluxes of the central carbon metabolism (CCM) by regulating the multi-level regulation machinery upon genetic perturbations and the change in the growth conditions. Here, we discuss about how aerobes accomplish such cofactor homeostasis against redox perturbations. In particular, we consider how single-gene mutants (including pgi, pfk, zwf, gnd and pyk mutants) modulate their metabolisms in relation to cofactor rebalance (and also by adaptive laboratory evolution). We also discuss about how the overproduction of NADPH (by the pathway gene mutation) can be utilized for the efficient production of useful value-added chemicals such as medicinal compounds, polyhydroxyalkanoates, and amino acids, all of which require NADPH in their synthetic pathways. We then discuss about the metabolic responses against oxidative stress, where αketoacids play important roles not only for the coordination between catabolism and anabolism, but also for detoxifying ROSs by non-enzymatic reactions, as well as for reducing the production of ROSs by repressing the activities of the TCA cycle and respiration (via carbon catabolite repression). Thus, we discuss about the mechanisms (basic strategies) that modulate the metabolism from respiration to respiro-fermentative metabolism causing overflow, based on the role of Pyk activity, affecting the NADPH production at the oxidative pentose phosphate (PP) pathway, and the roles of αketoacids for the change in the source of energy generation from the oxidative phosphorylation to the substrate level phosphorylation.

Developing a pyruvate-driven metabolic scenario for growth-coupled microbial production

Microbial-based chemical synthesis serves as a promising approach for sustainable production of industrially important products. However, limited production performance caused by metabolic burden or genetic variations poses one of the major challenges in achieving an economically viable biomanufacturing process. To address this issue, one superior strategy is to couple the product synthesis with cellular growth, which renders production obligatory for cell survival. Here we create a pyruvate-driven metabolic scenario in engineered Escherichia coli for growth-coupled bioproduction, with which we demonstrate its application in boosting production of anthranilate and its derivatives. Deletion of a minimal set of endogenous pyruvate-releasing pathways engenders anthranilate synthesis as the salvage route for pyruvate generation to support cell growth, concomitant with simultaneous anthranilate production. Further introduction of native and non-native downstream pathways affords production enhancement of two anthranilate-derived high-value products including L-tryptophan and cis, cis-muconic acid from different carbon sources. The work reported here presents a new growth-coupled strategy with demonstrated feasibility for promoting microbial production.

Mixed-species biofilms for high-cell-density application of Synechocystis sp. PCC 6803 in capillary reactors for continuous cyclohexane oxidation to cyclohexanol

Photosynthetic microorganisms have enormous potential to produce fuels and value-added compounds sustainably. Efficient cultivation concepts that enable optimal light and CO₂ supply are necessary for the realization of high cell densities (HCDs), and subsequently for process implementation. We introduce capillary biofilm reactors with a high surface to volume ratio, and thus enhanced light availability, enabling HCDs of photo-autotrophic microorganisms. However, oxygenic photosynthesis leads to O₂ accumulation in such systems, impairing biofilm growth. We combined O₂ producing Synechocystis with O₂ respiring Pseudomonas using proto-cooperation to achieve HCDs of up to 51.8 g_BDW L^-1. This concept was coupled to the challenging C-H oxyfunctionalization of cyclohexane to cyclohexanol with a remarkable conversion of >98% and selectivity of 100% (KA oil). High photoautotrophic biocatalyst concentrations were established and resulted in a productivity of 3.76 g_cyclohexanol m^-2 day^-1, which was maintained for at least one month.

The Manganese-Dependent Pyruvate Kinase PykM Is Required for Wild-Type Glucose Utilization by Brucella abortus 2308 and Its Virulence in C57BL/6 Mice

Pyruvate kinase plays a central role in glucose catabolism in bacteria, and efficient utilization of this hexose has been linked to the virulence of Brucella strains in mice. The brucellae produce a single pyruvate kinase which is an ortholog of the Bradyrhizobium manganese (Mn)-dependent pyruvate kinase PykM. A biochemical analysis of the Brucella pyruvate kinase and phenotypic analysis of a Brucella abortus mutant defective in high-affinity Mn import indicate that this enzyme is an authentic PykM ortholog which functions as a Mn-dependent enzyme in vivo The loss of PykM has a negative impact on the capacity of the parental 2308 strain to utilize glucose, fructose, and galactose but not on its ability to utilize ribose, xylose, arabinose, or erythritol, and a pykM mutant displays significant attenuation in C57BL/6 mice. Although the enzyme pyruvate phosphate dikinase (PpdK) can substitute for the loss of pyruvate kinase in some bacteria and is also an important virulence determinant in Brucella, a phenotypic analysis of B. abortus 2308 and isogenic pykM, ppdK, and pykM ppdK mutants indicates that PykM and PpdK make distinctly different contributions to carbon metabolism and virulence in these bacteria.IMPORTANCE Mn plays a critical role in the physiology and virulence of Brucella strains, and the results presented here suggest that one of the important roles that the high-affinity Mn importer MntH plays in the pathogenesis of these strains is supporting the function of the Mn-dependent kinase PykM. A better understanding of how the brucellae adapt their physiology and metabolism to sustain their intracellular persistence in host macrophages will provide knowledge that can be used to design improved strategies for preventing and treating brucellosis, a disease that has a significant impact on both the veterinary and public health communities worldwide.

Advances in metabolic flux analysis toward genome-scale profiling of higher organisms

Methodological and technological advances have recently paved the way for metabolic flux profiling in higher organisms, like plants. However, in comparison with omics technologies, flux profiling has yet to provide comprehensive differential flux maps at a genome-scale and in different cell types, tissues, and organs. Here we highlight the recent advances in technologies to gather metabolic labeling patterns and flux profiling approaches. We provide an opinion of how recent local flux profiling approaches can be used in conjunction with the constraint-based modeling framework to arrive at genome-scale flux maps. In addition, we point at approaches which use metabolomics data without introduction of label to predict either non-steady state fluxes in a time-series experiment or flux changes in different experimental scenarios. The combination of these developments allows an experimentally feasible approach for flux-based large-scale systems biology studies.

Evolutionary Convergence of Pathway-Specific Enzyme Expression Stoichiometry

Coexpression of proteins in response to pathway-inducing signals is the founding paradigm of gene regulation. However, it remains unexplored whether the relative abundance of co-regulated proteins requires precise tuning. Here, we present large-scale analyses of protein stoichiometry and corresponding regulatory strategies for 21 pathways and 67-224 operons in divergent bacteria separated by 0.6-2 billion years. Using end-enriched RNA-sequencing (Rend-seq) with single-nucleotide resolution, we found that many bacterial gene clusters encoding conserved pathways have undergone massive divergence in transcript abundance and architectures via remodeling of internal promoters and terminators. Remarkably, these evolutionary changes are compensated post-transcriptionally to maintain preferred stoichiometry of protein synthesis rates. Even more strikingly, in eukaryotic budding yeast, functionally analogous proteins that arose independently from bacterial counterparts also evolved to convergent in-pathway expression. The broad requirement for exact protein stoichiometries despite regulatory divergence provides an unexpected principle for building biological pathways both in nature and for synthetic activities.

A Network Biology Approach to Decipher Stress Response in Bacteria Using Escherichia coli As a Model

The development of drug-resistant pathogenic bacteria poses challenges to global health for their treatment and control. In this context, stress response enables bacterial populations to survive extreme perturbations in the environment but remains poorly understood. Specific modules are activated for unique stressors with few recognized global regulators. The phenomenon of cross-stress protection strongly suggests the presence of central proteins that control the diverse stress responses. In this work, Escherichia coli was used to model the bacterial stress response. A Protein-Protein Interaction Network was generated by integrating differentially expressed genes in eight stress conditions of pH, temperature, and antibiotics with relevant gene ontology terms. Topological analysis identified 24 central proteins. The well-documented role of 16 central proteins in stress indicates central control of the response, while the remaining eight proteins may have a novel role in stress response. Cluster analysis of the generated network implicated RNA binding, flagellar assembly, ABC transporters, and DNA repair as important processes during response to stress. Pathway analysis showed crosstalk of Two Component Systems with metabolic processes, oxidative phosphorylation, and ABC transporters. The results were further validated by analysis of an independent cross-stress protection dataset. This study also reports on the ways in which bacterial stress response can progress to biofilm formation. In conclusion, we suggest that drug targets or pathways disrupting bacterial stress responses can potentially be exploited to combat antibiotic tolerance and multidrug resistance in the future.

Expanding Primary Metabolism Helps Generate the Metabolic Robustness To Facilitate Antibiotic Biosynthesis in Streptomyces

The expansion of the genetic repertoire of an organism by gene duplication or horizontal gene transfer (HGT) can aid adaptation. Streptomyces bacteria are prolific producers of bioactive specialized metabolites that have adaptive functions in nature and have found extensive utility in human medicine. While the biosynthesis of these specialized metabolites is directed by dedicated biosynthetic gene clusters, little attention has been focused on how these organisms have evolved robustness in their genomes to facilitate the metabolic plasticity required to provide chemical precursors for biosynthesis during the complex metabolic transitions from vegetative growth to specialized metabolite production and sporulation. Here, we examine genetic redundancy in actinobacteria and show that specialized metabolite-producing bacterial families exhibit gene family expansion in primary metabolism. Focusing on a gene duplication event, we show that the two pyruvate kinases in the genome of Streptomyces coelicolor arose by an ancient duplication event and that each has evolved altered enzymatic kinetics, with Pyk1 having a 20-fold-higher kcat than Pyk2 (4,703 s-1 compared to 215 s-1, respectively), and yet both are constitutively expressed. The pyruvate kinase mutants were also found to be compromised in terms of fitness compared to wild-type Streptomyces These data suggest that expanding gene families can help maintain cell functionality during metabolic perturbation such as nutrient limitation and/or specialized metabolite production.IMPORTANCE The rise of antimicrobial-resistant infections has prompted a resurgence in interest in understanding the production of specialized metabolites, such as antibiotics, by Streptomyces The presence of multiple genes encoding the same enzymatic function is an aspect of Streptomyces biology that has received little attention; however, understanding how the metabolic expansion influences these organisms can help enhance production of clinically useful molecules. Here, we show that expanding the number of pyruvate kinases enables metabolic adaptation, increases strain fitness, and represents an excellent target for metabolic engineering of industrial specialized metabolite-producing bacteria and the activation of cryptic specialized metabolites.

Dietary energy drives the dynamic response of bovine rumen viral communities

BACKGROUND

Rumen microbes play a greater role in host energy acquisition than that of gut-associated microbes in monogastric animals. Although genome-enabled advancements are providing access to the vast diversity of uncultivated microbes, our understanding of variables shaping rumen microbial communities is in its infancy. Viruses have been shown to impact microbial populations through a myriad of processes, including cell lysis and reprogramming of host metabolism. However, little is known about the processes shaping the distribution of rumen viruses or how viruses may modulate microbial-driven processes in the rumen. To this end, we investigated how rumen bacterial and viral community structure and function responded in five steers fed four randomized dietary treatments in a crossover design.

RESULTS

Total digestible nutrients (TDN), a measure of dietary energy, best explained the variation in bacterial and viral communities. Additional ecological drivers of viral communities included dietary zinc content and microbial functional diversity. Using partial least squares regression, we demonstrate significant associations between the abundances of 267 viral populations and variables driving the variation in rumen viral communities. While rumen viruses were dynamic, 14 near ubiquitous viral populations were identified, suggesting the presence of a core rumen virome largely comprised of novel viruses. Moreover, analysis of virally encoded auxiliary metabolic genes (AMGs) indicates rumen viruses have glycosidic hydrolases to potentially augment the breakdown of complex carbohydrates to increase energy production. Other AMGs identified have a role in redirecting carbon to the pentose phosphate pathway and one carbon pools by folate to boost viral replication.

CONCLUSIONS

We demonstrate that rumen bacteria and viruses have differing responses and ecological drivers to dietary perturbation. Our results show that rumen viruses have implications for understanding the structuring of the previously identified core rumen microbiota and impacting microbial metabolism through a vast array of AMGs. AMGs in the rumen appear to have consequences for microbial metabolism that are largely in congruence with the current paradigm established in marine systems. This study provides a foundation for future hypotheses regarding the dynamics of viral-mediated processes in the rumen.

Biocatalytic conversion of cycloalkanes to lactones using an in-vivo cascade in Pseudomonas taiwanensis VLB120

Chemical synthesis of lactones from cycloalkanes is a multi-step process challenged by limitations in reaction efficiency (conversion and yield), atom economy (by-products) and environmental performance. A heterologous pathway comprising novel enzymes with compatible kinetics was designed in Pseudomonas taiwanensis VLB120 enabling in-vivo cascade for synthesizing lactones from cycloalkanes. The respective pathway included cytochrome P450 monooxygenase (CHX), cyclohexanol dehydrogenase (CDH), and cyclohexanone monooxygenase (CHXON) from Acidovorax sp. CHX100. Resting (non-growing) cells of the recombinant host P. taiwanensis VLB120 converted cyclohexane, cyclohexanol, and cyclohexanone to ϵ-caprolactone at 22, 80-100, and 170 U g_CDW^-1 , respectively. Cyclohexane (5 mM) was completely converted with a selectivity of 65% for ϵ-caprolactone formation in 2 hr without accumulation of intermediate products. Promiscuity of the whole-cell biocatalyst gave access to analogous lactones from cyclooctane and cyclodecane. A total product concentration of 2.3 g L^-1 and a total turnover number of 36,720 was achieved over 5 hr with a biocatalyst concentration of 6.8 g_CDW  L^-1 .

Proteomics and 1H NMR-based metabolomics analysis of pathogenic Vibrio vulnificus aquacultures isolated from sewage drains

Vibrio bacteria live in both marine and freshwater habitats and are associated with aquatic animals. Vibrio vulnificus is a pathogenic bacterium that infects people and livestock. It is usually found in offshore waters or within fish and shellfish. This study presents a comparative proteomic analysis of the outer membrane protein (OMP) changes in V. vulnificus proteins after stimulation with sewage from sewage drains. Using two-dimensional electrophoresis followed by MALDI-TOF MS/MS, 32 protein spots with significant differences in abundance were identified and characterized. These identified proteins were found to be involved in various functional categories, including catalysis, transport, membrane proteins progresses, receptor activity, energy metabolism, cytokine activity, and protein metabolism. The mRNA expression levels of 12 differential proteins were further assessed by qRT-PCR. Seven genes including carboxypeptidase, hemoglobin receptor, succinate dehydrogenase iron-sulfur subunit, ATP synthase subunit alpha, thioredoxin, succinyl-CoA synthetase subunit, and alanine dehydrogenase were downregulated upon stimulation, whereas the protein expression levels HupA receptor, type I secretion outer membrane protein, glutamine synthetase, superoxide dismutase, OmpU, and VuuA were upregulated. ¹H NMR spectra showed 18 dysregulated metabolites from V. vulnificus after the sewage stimulation and the pathogenicity was enhanced after that.

13 C metabolic flux profiling of Pichia pastoris grown in aerobic batch cultures on glucose revealed high relative anabolic use of TCA cycle and limited incorporation of provided precursors of branched-chain amino acids

Carbon metabolism of Crabtree-negative yeast Pichia pastoris was profiled using ¹³ C nuclear magnetic resonance (NMR) to delineate regulation during exponential growth and to study the import of two precursors for branched-chain amino acid biosynthesis, α-ketoisovalerate and α-ketobutyrate. Cells were grown in aerobic batch cultures containing (a) only glucose, (b) glucose along with the precursors, or (c) glucose and Val. The study provided the following new insights. First, ¹³ C flux ratio analyses of central metabolism reveal an unexpectedly high anaplerotic supply of the tricarboxylic acid cycle for a Crabtree-negative yeast, and show that a substantial fraction of glucose catabolism proceeds through the pentose phosphate pathway. A comparison with previous flux ratio analyses for batch cultures of Crabtree-negative Pichia stipitis and Crabtree-positive Saccharomyces cerevisiae indicate that the overall regulation of central carbon metabolism in P. pastoris is intermediate in between P. stipitis and S. cerevisiae. Second, excess α-ketoisovalerate in the medium is not transported into the cytoplasm indicating that P. pastoris lacks a suitable transporter. In contrast, excess Val is efficiently taken up and largely fulfills demands for both Val and Leu for protein synthesis. Third, excess α-ketobutyrate is transported into the mitochondria for Ile biosynthesis. However, the import does not efficiently inhibit the synthesis of α-ketobutyrate from pyruvate indicating that P. pastoris has not been optimized evolutionarily to take full advantage of this carbon source. These findings have direct implications for preparing uniformly ² H,¹³ C,¹⁵ N-labeled proteins containing protonated Ile, Val, and Leu methyl groups in P. pastoris for NMR-based structural biology.

ENZYMES

Acetohydroxy acid isomeroreductase (EC 1.1.1.86), branched-chain amino acid aminotransferase (BCAT, EC 2.6.1.42), fumarase (EC 4.2.1.2), malic enzyme (EC 1.1.1.39/1.1.1.40), phosphoenolpyruvate carboxykinase (EC 4.1.1.49), pyruvate carboxylase (EC 6.4.1.1), pyruvate kinase (EC 2.7.1.40), l-serine hydroxymethyltransferase (EC 2.1.2.1), threonine aldolase (EC 4.1.2.5), threonine dehydratase (EC 4.3.1.19); transketolase (EC 2.2.1.1), transaldolase (EC 2.2.1.2).

Beyond topology: coevolution of structure and flux in metabolic networks

Interactions between the structure of a metabolic network and its functional properties underlie its evolutionary diversification, but the mechanism by which such interactions arise remains elusive. Particularly unclear is whether metabolic fluxes that determine the concentrations of compounds produced by a metabolic network, are causally linked to a network's structure or emerge independently of it. A direct empirical study of populations where both structural and functional properties vary among individuals' metabolic networks is required to establish whether changes in structure affect the distribution of metabolic flux. In a population of house finches (Haemorhous mexicanus), we reconstructed full carotenoid metabolic networks for 442 individuals and uncovered 11 structural variants of this network with different compounds and reactions. We examined the consequences of this structural diversity for the concentrations of plumage-bound carotenoids produced by flux in these networks. We found that concentrations of metabolically derived, but not dietary carotenoids, depended on network structure. Flux was partitioned similarly among compounds in individuals of the same network structure: within each network, compound concentrations were closely correlated. The highest among-individual variation in flux occurred in networks with the strongest among-compound correlations, suggesting that changes in the magnitude, but not the distribution of flux, underlie individual differences in compound concentrations on a static network structure. These findings indicate that the distribution of flux in carotenoid metabolism closely follows network structure. Thus, evolutionary diversification and local adaptations in carotenoid metabolism may depend more on the gain or loss of enzymatic reactions than on changes in flux within a network structure.

Carbonyl Stress in Bacteria: Causes and Consequences

Pathways of synthesis of the α-reactive carbonyl compound methylglyoxal (MG) in prokaryotes are described in this review. Accumulation of MG leads to development of carbonyl stress. Some pathways of MG formation are similar for both pro- and eukaryotes, but there are reactions specific for prokaryotes, e.g. the methylglyoxal synthase reaction. This reaction and the glyoxalase system constitute an alternative pathway of glucose catabolism - the MG shunt not associated with the synthesis of ATP. In violation of the regulation of metabolism, the cell uses MG shunt as well as other glycolysis shunting pathways and futile cycles enabling stabilization of its energetic status. MG was first examined as a biologically active metabolic factor participating in the formation of phenotypic polymorphism and hyperpersistent potential of bacterial populations. The study of carbonyl stress is interesting for evolutionary biology and can be useful for constructing highly effective producer strains.

Growth of Pseudomonas taiwanensis VLB120∆C biofilms in the presence of n-butanol

Biocatalytic processes often encounter problems due to toxic reactants and products, which reduce biocatalyst viability. Thus, robust organisms capable of tolerating or adapting towards such compounds are of high importance. This study systematically investigated the physiological response of Pseudomonas taiwanensisVLB120∆C biofilms when exposed to n‐butanol, one of the potential next generation biofuels as well as a toxic substance using microscopic and biochemical methods. Initially P. taiwanensisVLB120∆C biofilms did not show any observable growth in the presence of 3% butanol. Prolonged cultivation of 10 days led to biofilm adaptation, glucose and oxygen uptake doubled and consequently it was possible to quantify biomass. Complementing the medium with yeast extract and presumably reducing the metabolic burden caused by butanol exposure further increased the biomass yield. In course of cultivation cells reduced their size in the presence of n‐butanol which results in an enlarged surface‐to‐volume ratio and thus increased nutrient uptake. Finally, biofilm enhanced its extracellular polymeric substances (EPS) production when exposed to n‐butanol. The predominant response of these biofilms under n‐butanol stress are higher energy demand, increased biomass yield upon medium complements, larger surface‐to‐volume ratio and enhanced EPS production. Although we observed a distinct increase in biomass in the presence of 3% butanol it was not possible to cultivate P. taiwanensisVLB120∆C biofilms at higher n‐butanol concentrations. Thereby this study shows that biofilms are not per se tolerant against solvents, and need to adapt to toxic n‐butanol concentrations.

SUMOFLUX: A Generalized Method for Targeted 13C Metabolic Flux Ratio Analysis

Metabolic fluxes are a cornerstone of cellular physiology that emerge from a complex interplay of enzymes, carriers, and nutrients. The experimental assessment of in vivo intracellular fluxes using stable isotopic tracers is essential if we are to understand metabolic function and regulation. Flux estimation based on 13C or 2H labeling relies on complex simulation and iterative fitting; processes that necessitate a level of expertise that ordinarily preclude the non-expert user. To overcome this, we have developed SUMOFLUX, a methodology that is broadly applicable to the targeted analysis of 13C-metabolic fluxes. By combining surrogate modeling and machine learning, we trained a predictor to specialize in estimating flux ratios from measurable 13C-data. SUMOFLUX targets specific flux features individually, which makes it fast, user-friendly, applicable to experimental design and robust in terms of experimental noise and exchange flux magnitude. Collectively, we predict that SUMOFLUX's properties realistically pave the way to high-throughput flux analyses.

PFA toolbox: a MATLAB tool for Metabolic Flux Analysis

BACKGROUND

Metabolic Flux Analysis (MFA) is a methodology that has been successfully applied to estimate metabolic fluxes in living cells. However, traditional frameworks based on this approach have some limitations, particularly when measurements are scarce and imprecise. This is very common in industrial environments. The PFA Toolbox can be used to face those scenarios.

RESULTS

Here we present the PFA (Possibilistic Flux Analysis) Toolbox for MATLAB, which simplifies the use of Interval and Possibilistic Metabolic Flux Analysis. The main features of the PFA Toolbox are the following: (a) It provides reliable MFA estimations in scenarios where only a few fluxes can be measured or those available are imprecise. (b) It provides tools to easily plot the results as interval estimates or flux distributions. (c) It is composed of simple functions that MATLAB users can apply in flexible ways. (d) It includes a Graphical User Interface (GUI), which provides a visual representation of the measurements and their uncertainty. (e) It can use stoichiometric models in COBRA format. In addition, the PFA Toolbox includes a User's Guide with a thorough description of its functions and several examples.

CONCLUSIONS

The PFA Toolbox for MATLAB is a freely available Toolbox that is able to perform Interval and Possibilistic MFA estimations.

Dynamics of benzoate metabolism in Pseudomonas putida KT2440

Soil microorganisms mineralize lignin-derived aromatic carbon sources using oxidative catabolic pathways, such as the β-ketoadipate pathway. Although this aromatic pathway is one of the best-studied pathways in biochemistry, the complete pathway, including its regulation by aromatic carbon sources, has not been integrated into the metabolic network. In particular, information about the in vivo operation (e.g., kinetics and flux capacity) of the pathway is lacking. In this contribution, we use kinetic modeling and thermodynamic analysis to evaluate the in vivo operation of this key aromatic multi-step pathway. The resulting ab initio deterministic model of benzoate degradation via the β-ketoadipate (ortho-cleavage) pathway in Pseudomonas putida KT2440 is presented. The kinetic model includes mechanistic rate expressions for the enzymes and transport processes. The design and experimental validation of the model are driven by data generated from short-term perturbation experiments in a benzoate-limited continuous culture. The results of rigorous modeling of the in vivo dynamics provide strong support for flux regulation by the benzoate transporter and the enzymes forming and cleaving catechol. Revisiting the β-ketoadipate pathway might be valuable for applications in different fields, such as biochemistry and metabolic engineering, that use lignin monomers as a carbon source.

•

We describe a kinetic model for the β-ketoadipate pathway.

•

Short term metabolic responses were tracked on metabolite level by rapid sampling.

•

The model captures steady state and dynamic conditions of the β-ketoadipate pathway.

•

Thermodynamic analysis revealed regulation points of the pathway.

•

The results are discussed in the context of metabolic network operation.

Web application for genetic modification flux with database to estimate metabolic fluxes of genetic mutants

Computational analysis of metabolic fluxes is essential in understanding the structure and function of a metabolic network and in rationally designing genetically modified mutants for an engineering purpose. We had presented the genetic modification flux (GMF) that predicts the flux distribution of a broad range of genetically modified mutants. To enhance the feasibility and usability of GMF, we have developed a web application with a metabolic network database to predict a flux distribution of genetically modified mutants. One hundred and twelve data sets of Escherichia coli, Corynebacterium glutamicum, Saccharomyces cerevisiae, and Chinese hamster ovary were registered as standard models.

Indispensability of Horizontally Transferred Genes and Its Impact on Bacterial Genome Streamlining

Why are certain bacterial genomes so small and compact? The adaptive genome streamlining hypothesis posits that selection acts to reduce genome size because of the metabolic burden of replicating DNA. To reveal the impact of genome streamlining on cellular traits, we reduced the Escherichia coli genome by up to 20% by deleting regions which have been repeatedly subjects of horizontal transfer in nature. Unexpectedly, horizontally transferred genes not only confer utilization of specific nutrients and elevate tolerance to stresses, but also allow efficient usage of resources to build new cells, and hence influence fitness in routine and stressful environments alike. Genome reduction affected fitness not only by gene loss, but also by induction of a general stress response. Finally, we failed to find evidence that the advantage of smaller genomes would be due to a reduced metabolic burden of replicating DNA or a link with smaller cell size. We conclude that as the potential energetic benefit gained by deletion of short genomic segments is vanishingly small compared with the deleterious side effects of these deletions, selection for reduced DNA synthesis costs is unlikely to shape the evolution of small genomes.

Continuous cyclohexane oxidation to cyclohexanol using a novel cytochrome P450 monooxygenase from Acidovorax sp. CHX100 in recombinant P. taiwanensis VLB120 biofilms

The applications of biocatalysts in chemical industries are characterized by activity, selectivity, and stability. One key strategy to achieve high biocatalytic activity is the identification of novel enzymes with kinetics optimized for organic synthesis by Nature. The isolation of novel cytochrome P450 monooxygenase genes from Acidovorax sp. CHX100 and their functional expression in recombinant Pseudomonas taiwanensis VLB120 enabled efficient oxidation of cyclohexane to cyclohexanol. Although initial resting cell activities of 20 U gCDW (-1) were achieved, the rapid decrease in catalytic activity due to the toxicity of cyclohexane prevented synthetic applications. Cyclohexane toxicity was reduced and cellular activities stabilized over the reaction time by delivering the toxic substrate through the vapor phase and by balancing the aqueous phase mass transfer with the cellular conversion rate. The potential of this novel CYP enzyme was exploited by transferring the shake flask reaction to an aqueous-air segmented flow biofilm membrane reactor for maximizing productivity. Cyclohexane was continuously delivered via the silicone membrane. This ensured lower reactant toxicity and continuous product formation at an average volumetric productivity of 0.4 g L tube (-1) h(-1) for several days. This highlights the potential of combining a powerful catalyst with a beneficial reactor design to overcome critical issues of cyclohexane oxidation to cyclohexanol. It opens new opportunities for biocatalytic transformations of compounds which are toxic, volatile, and have low solubility in water.

Fluxome study of Pseudomonas fluorescens reveals major reorganisation of carbon flux through central metabolic pathways in response to inactivation of the anti-sigma factor MucA

Background

The bacterium Pseudomonas fluorescens switches to an alginate-producing phenotype when the pleiotropic anti-sigma factor MucA is inactivated. The inactivation is accompanied by an increased biomass yield on carbon sources when grown under nitrogen-limited chemostat conditions. A previous metabolome study showed significant changes in the intracellular metabolite concentrations, especially of the nucleotides, in mucA deletion mutants compared to the wild-type. In this study, the P. fluorescens SBW25 wild-type and an alginate non-producing mucA- ΔalgC double-knockout mutant are investigated through model-based ¹³C-metabolic flux analysis (¹³C-MFA) to explore the physiological consequences of MucA inactivation at the metabolic flux level. Intracellular metabolite extracts from three carbon labelling experiments using fructose as the sole carbon source are analysed for ¹³C-label incorporation in primary metabolites by gas and liquid chromatography tandem mass spectrometry.

Results

From mass isotopomer distribution datasets, absolute intracellular metabolic reaction rates for the wild type and the mutant are determined, revealing extensive reorganisation of carbon flux through central metabolic pathways in response to MucA inactivation. The carbon flux through the Entner-Doudoroff pathway was reduced in the mucA- ΔalgC mutant, while flux through the pentose phosphate pathway was increased. Our findings also indicated flexibility of the anaplerotic reactions through down-regulation of the pyruvate shunt in the mucA- ΔalgC mutant and up-regulation of the glyoxylate shunt.

Conclusions

Absolute metabolic fluxes and metabolite levels give detailed, integrated insight into the physiology of this industrially, medically and agriculturally important bacterial species and suggest that the most efficient way of using a mucA- mutant as a cell factory for alginate production would be to use non-growing conditions and nitrogen deprivation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0148-0) contains supplementary material, which is available to authorized users.

Tailoring strain construction strategies for muconic acid production in S. cerevisiae and E. coli

There is currently a strong interest to derive the biological precursor cis,cis-muconic acid from shikimate pathway-branches to develop a biological replacement for adipic acid. Pioneered by the Frost laboratory this concept has regained interest: Recent approaches (Boles, Alper, Yan) however suffer from low product titres. Here an in silico comparison of all strain construction strategies was conducted to highlight stoichiometric optimizations. Using elementary mode analysis new knock-out strategies were determined in Saccharomyces cerevisiae and Escherichia coli. The strain construction strategies are unique to each pathway-branch and organism, allowing significantly different maximum and minimum yields. The maximum theoretical product carbon yields on glucose ranged from 86% (dehydroshikimate-branch) to 69% (anthranilate-branch). In most cases a coupling of product formation to growth was possible. Especially in S. cerevisiae chorismate-routes a minimum yield constraint of 46.9% could be reached. The knock-out targets are non-obvious, and not-transferable, highlighting the importance of tailored strain construction strategies.

•

Maximum theoretical product carbon yields for biological production of cis,cis-muconic acid from glucose were calculated.

•

Existing strain construction strategies were compared in silico and weaknesses were identified.

•

New strain construction strategies were determined to allow coupling of product formation to growth for E. coli and S. cerevisiae.

•

In particular abolishing one target reaction introduces a minimum yield constraint of 28% in S. cerevisiae.

•

Determined knock-out targets are unique to each pathway and organism, highlighting the importance of tailored strain construction strategies.

Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C-metabolic flux analysis

We engineered a fatty acid overproducing Escherichia coli strain through overexpressing tesA (“pull”) and fadR (“push”) and knocking out fadE (“block”). This “pull-push-block” strategy yielded 0.17 g of fatty acids (C12–C18) per gram of glucose (equivalent to 48% of the maximum theoretical yield) in batch cultures during the exponential growth phase under aerobic conditions. Metabolic fluxes were determined for the engineered E. coli and its control strain using tracer ([1,2-13C]glucose) experiments and 13C-metabolic flux analysis. Cofactor (NADPH) and energy (ATP) balances were also investigated for both strains based on estimated fluxes. Compared to the control strain, fatty acid overproduction led to significant metabolic responses in the central metabolism: (1) Acetic acid secretion flux decreased 10-fold; (2) Pentose phosphate pathway and Entner–Doudoroff pathway fluxes increased 1.5- and 2.0-fold, respectively; (3) Biomass synthesis flux was reduced 1.9-fold; (4) Anaplerotic phosphoenolpyruvate carboxylation flux decreased 1.7-fold; (5) Transhydrogenation flux converting NADH to NADPH increased by 1.7-fold. Real-time quantitative RT-PCR analysis revealed the engineered strain increased the transcription levels of pntA (encoding the membrane-bound transhydrogenase) by 2.1-fold and udhA (encoding the soluble transhydrogenase) by 1.4-fold, which is in agreement with the increased transhydrogenation flux. Cofactor and energy balances analyses showed that the fatty acid overproducing E. coli consumed significantly higher cellular maintenance energy than the control strain. We discussed the strategies to future strain development and process improvements for fatty acid production in E. coli.

The pyruvate-tricarboxylic acid cycle node: a focal point of virulence control in the enteric pathogen Yersinia pseudotuberculosis

Despite our increasing knowledge of the specific pathogenicity factors in bacteria, the contribution of metabolic processes to virulence is largely unknown. Here, we elucidate a tight connection between pathogenicity and core metabolism in the enteric pathogen Yersinia pseudotuberculosis by integrated transcriptome and [(13)C]fluxome analysis of the wild type and virulence-regulator mutants. During aerobic growth on glucose, Y. pseudotuberculosis reveals an unusual flux distribution with a high level of secreted pyruvate. The absence of the transcriptional and post-transcriptional regulators RovA, CsrA, and Crp strongly perturbs the fluxes of carbon core metabolism at the level of pyruvate metabolism and the tricarboxylic acid (TCA) cycle, and these perturbations are accompanied by transcriptional changes in the corresponding enzymes. Knock-outs of regulators of this metabolic branch point and of its central enzyme, pyruvate kinase (ΔpykF), result in mutants with significantly reduced virulence in an oral mouse infection model. In summary, our work identifies the pyruvate-TCA cycle node as a focal point for controlling the host colonization and virulence of Yersinia.

Comparison and analysis of objective functions in flux balance analysis

Flux balance analysis (FBA) is currently one of the most important and used techniques for estimation of metabolic reaction rates (fluxes). This mathematical approach utilizes an optimization criterion in order to select a distribution of fluxes from the feasible space delimited by the metabolic reactions and some restrictions imposed over them, assuming that cellular metabolism is in steady state. Therefore, the obtained flux distribution depends on the specific objective function used. Multiple studies have been aimed to compare distinct objective functions at given conditions, in order to determine which of those functions produces values of fluxes closer to real data when used as objective in the FBA; in other words, what is the best objective function for modeling cell metabolism at a determined environmental condition. However, these comparative studies have been designed in very dissimilar ways, and in general, several factors that can change the ideal objective function in a cellular condition have not been adequately considered. Additionally, most of them have used only one dataset for representing one condition of cell growth, and different measuring techniques have been used. For these reasons, a rigorous study on the effect of factors such as the quantity of used data, the number and type of fluxes utilized as input data, and the selected classification of growth conditions, are required in order to obtain useful conclusions for these comparative studies, allowing limiting clearly the application range on any of those results.

The functional structure of central carbon metabolism in Pseudomonas putida KT2440

What defines central carbon metabolism? The classic textbook scheme of central metabolism includes the Embden-Meyerhof-Parnas (EMP) pathway of glycolysis, the pentose phosphate pathway, and the citric acid cycle. The prevalence of this definition of central metabolism is, however, equivocal without experimental validation. We address this issue using a general experimental approach that combines the monitoring of transcriptional and metabolic flux changes between steady states on alternative carbon sources. This approach is investigated by using the model bacterium Pseudomonas putida with glucose, fructose, and benzoate as carbon sources. The catabolic reactions involved in the initial uptake and metabolism of these substrates are expected to show a correlated change in gene expressions and metabolic fluxes. However, there was no correlation for the reactions linking the 12 biomass precursor molecules, indicating a regulation mechanism other than mRNA synthesis for central metabolism. This result substantiates evidence for a (re)definition of central carbon metabolism including all reactions that are bound to tight regulation and transcriptional invariance. Contrary to expectations, the canonical Entner-Doudoroff and EMP pathways sensu stricto are not a part of central carbon metabolism in P. putida, as they are not regulated differently from the aromatic degradation pathway. The regulatory analyses presented here provide leads on a qualitative basis to address the use of alternative carbon sources by deregulation and overexpression at the transcriptional level, while rate improvements in central carbon metabolism require careful adjustment of metabolite concentrations, as regulation resides to a large extent in posttranslational and/or metabolic regulation.