Robustness and plasticity of metabolic pathway flux among uropathogenic isolates of Pseudomonas aeruginosa

Characterization of the Entner-Doudoroff pathway in Pseudomonas aeruginosa catheter-associated urinary tract infections

Pseudomonas aeruginosa is an opportunistic nosocomial pathogen responsible for a subset of catheter-associated urinary tract infections (CAUTI). In a murine model of P. aeruginosa CAUTI, we previously demonstrated that urea within urine suppresses quorum sensing and induces the Entner-Doudoroff (E-D) pathway. The E-D pathway consists of the genes zwf, pgl, edd, and eda. Zwf and Pgl convert glucose-6-phosphate into 6-phosphogluconate. Edd hydrolyzes 6-phosphogluconate to 2-keto-3-deoxy-6-phosphogluconate (KDPG). Finally, Eda cleaves KDPG to glyceraldehyde-3-phosphate and pyruvate, which enters the citric acid cycle. Here, we generated in-frame E-D mutants in the strain PA14 and assessed their growth phenotypes on chemically defined and complex media. These E-D mutants have a growth defect when grown on glucose or gluconate as the sole carbon source, which is similar to results previously reported for PAO1 mutants lacking E-D genes. RNA-sequencing following short exposure to urine revealed minimal gene regulation differences compared to the wild type. In a murine CAUTI model, virulence testing of E-D mutants revealed that two mutants lacking zwf and pgl showed minor fitness defects. Infection with the ∆pgl strain exhibited a 20% increase in host survival, and the ∆zwf strain displayed decreased colonization of the catheter and kidneys. Consequently, our findings suggest that the E-D pathway in P. aeruginosa is dispensable in this model of CAUTI. IMPORTANCE Prior studies have shown that the Entner-Doudoroff pathway is up-regulated when Pseudomonas aeruginosa is grown in urine. Pseudomonads use the Entner-Doudoroff (E-D) pathway to metabolize glucose instead of glycolysis, which led us to ask whether this pathway is required for urinary tract infection. Here, single-deletion mutants of each gene in the pathway were tested for growth on chemically defined media with single-carbon sources as well as complex media. The effect of each mutant on global gene expression in laboratory media and urine was characterized. The virulence of these mutants in a murine model of catheter-associated urinary tract infection revealed that these mutants had similar levels of colonization indicating that glucose is not the primary carbon source utilized in the urinary tract.

Diversity and function of viral AMGs associated with DNA biosynthesis in the Napahai plateau wetland

Viruses play an important role in microbial community structure and biodiversity by lysing host cells, and can also affect host metabolic pathways by expressing auxiliary metabolic genes (AMGs). As a unique low-latitude, high-altitude seasonal plateau wetland in China, Napahai has high research value. However, studies on the genetic diversity of AMGs and viruses associated with DNA biosynthesis have not been reported. Based on metagenomics, with the phylogenetic tree, PCoA, and α diversity analysis, we found that three DNA biosynthesis-related viral AMGs (cobS, mazG, and purM) in the Napahai plateau wetland were rich in genetic diversity, uniqueness, and differences compared with other habitats and host sources. Through the KEGG metabolic pathway and metabolic flow analysis of Pseudomonas mandelii (SW-3) and phage (VSW-3), the AMGs (cobS, mazG, and purM) genes of the three related viruses involved in DNA biosynthesis were upregulated and their expression increased significantly. In general, we systematically described the genetic diversity of AMGs associated with DNA biosynthesis in plateau wetland ecosystems and clarified the contribution of viral AMGs in the Napahai plateau wetland to DNA biosynthesis, as well as the changes of metabolites and genes. It further expands the understanding of phage-host interactions, which is of great significance for further revealing the role of viral AMGs in the biological evolution and biogeochemical cycle of wetland ecosystems.

Characterization of the Entner-Douderoff Pathway in Pseudomonas aeruginosa Catheter-associated Urinary Tract Infections

Pseudomonas aeruginosa is an opportunistic nosocomial pathogen responsible for catheter-associated urinary tract infections (CAUTI). In a murine model of P. aeruginosa CAUTI, we previously demonstrated that urea within urine suppresses quorum sensing and induces the Entner-Douderoff (E-D) pathway. The E-D pathway consists of the genes zwf, pgl, edd, and eda. Zwf and Pgl convert glucose-6-phosphate into 6-phosphogluconate. Edd hydrolyzes 6-phosphogluconate to 2-keto-3-deoxy-6-phosphogluconate (KDPG). Finally, Eda cleaves KDPG to glyceraldehyde-3-phosphate and pyruvate, which enters the citric acid cycle. Here, we generated in-frame E-D mutants in strain PA14 and assessed their growth phenotypes on chemically defined media. These E-D mutants have a growth defect when grown on glucose or gluconate as sole carbon source which are similar to results previously reported for PAO1 mutants lacking E-D genes. RNA-sequencing following short exposure to urine revealed minimal gene regulation differences compared to the wild type. In a murine CAUTI model, virulence testing of E-D mutants revealed that two mutants lacking zwf and pgl showed minor fitness defects. Infection with the ∆pgl strain exhibited a 20% increase in host survival, and the ∆zwf strain displayed decreased colonization of the catheter and kidneys. Consequently, our findings suggest that the E-D pathway in P. aeruginosa is dispensable in this model of CAUTI.

Systems-Wide Dissection of Organic Acid Assimilation in Pseudomonas aeruginosa Reveals a Novel Path To Underground Metabolism

The human pathogen Pseudomonas aeruginosa (Pa) is one of the most frequent and severe causes of nosocomial infection. This organism is also a major cause of airway infections in people with cystic fibrosis (CF). Pa is known to have a remarkable metabolic plasticity, allowing it to thrive under diverse environmental conditions and ecological niches; yet, little is known about the central metabolic pathways that sustain its growth during infection or precisely how these pathways operate. In this work, we used a combination of 'omics approaches (transcriptomics, proteomics, metabolomics, and ¹³C-fluxomics) and reverse genetics to provide systems-level insight into how the infection-relevant organic acids succinate and propionate are metabolized by Pa. Moreover, through structural and kinetic analysis of the 2-methylcitrate synthase (2-MCS; PrpC) and its paralogue citrate (CIT) synthase (GltA), we show how these two crucial enzymatic steps are interconnected in Pa organic acid assimilation. We found that Pa can rapidly adapt to the loss of GltA function by acquiring mutations in a transcriptional repressor, which then derepresses prpC expression. Our findings provide a clear example of how "underground metabolism," facilitated by enzyme substrate promiscuity, "rewires" Pa metabolism, allowing it to overcome the loss of a crucial enzyme. This pathogen-specific knowledge is critical for the advancement of a model-driven framework to target bacterial central metabolism. IMPORTANCE Pseudomonas aeruginosa is an opportunistic human pathogen that, due to its unrivalled resistance to antibiotics, ubiquity in the built environment, and aggressiveness in infection scenarios, has acquired the somewhat dubious accolade of being designated a "critical priority pathogen" by the WHO. In this work, we uncover the pathways and mechanisms used by P. aeruginosa to grow on a substrate that is abundant at many infection sites: propionate. We found that if the organism is prevented from metabolizing propionate, the substrate turns from being a convenient nutrient source into a potent poison, preventing bacterial growth. We further show that one of the enzymes involved in these reactions, 2-methylcitrate synthase (PrpC), is promiscuous and can moonlight for another essential enzyme in the cell (citrate synthase). Indeed, mutations that abolish citrate synthase activity (which would normally prevent the cell from growing) can be readily overcome if the cell acquires additional mutations that increase the expression of PrpC. This is a nice example of the evolutionary utility of so-called "underground metabolism."

Repressed Central Carbon Metabolism and Its Effect on Related Metabolic Pathways in Cefoperazone/Sulbactam-Resistant Pseudomonas aeruginosa

Metabolic shift and antibiotic resistance have been reported in Pseudomonas aeruginosa. However, the global metabolic characteristics remain largely unknown. The present study characterizes the central carbon metabolism and its effect on other metabolic pathways in cefoperazone-sulbactam (SCF)-resistant P. aeruginosa (PA-R_SCF). GC-MS-based metabolomics shows a repressed central carbon metabolism in PA-R_SCF, which is confirmed by measuring expression of genes and activity of enzymes in the metabolism. Furthermore, expression of the genes that encode the enzymes for the first step of fatty acid biosynthesis, glutamate metabolism, and electron transport chain is reduced, confirmed by their enzymatic activity assay, and the key enzyme for riboflavin metabolism is also reduced, indicating the decreased metabolic flux to the four related metabolic pathways. Moreover, the role of the reduced riboflavin metabolism, being related to ROS generation, in SCF resistance is explored. Exogenous H₂O₂ potentiates SCF-mediated killing in a dose-dependent manner, suggesting that the decreased ROS resulted from the reduced riboflavin metabolism that contributed to the resistance. These results indicate that the repressed central carbon metabolism and related riboflavin metabolism contribute to SCF resistance, but increasing ROS can restore SCF sensitivity. These findings characterize the repressed central carbon metabolism and its effect on other metabolic pathways as the global metabolic features in PA-R_SCF.

A study of the flexibility of the carbon catabolic pathways of extremophilic P. aeruginosa san ai exposed to benzoate versus glucose as sole carbon sources by multi omics analytical platform

Polyextremophilic, hydrocarbonoclastic Pseudomonas aeruginosa san ai can survive under extreme environmental challenges in the presence of a variety of pollutants such as organic solvents and hydrocarbons, particularly aromatics, heavy metals, and high pH. To date, the metabolic plasticity of the extremophilic P. aeruginosa, has not been sufficiently studied in regard to the effect of changing carbon sources. Therefore, the present study explores the carbon metabolic pathways of polyextremophilic P. aeruginosa san ai grown on sodium benzoate versus glucose and its potential for aromatic degradation. P. aeruginosa san ai removed/metabolised nearly 430 mg/L of benzoate for 48 h, demonstrating a high capacity for aromatic degradation. Comparative functional proteomics, targeted metabolomics and genomics analytical approaches were employed to study the carbon metabolism of the P. aeruginosa san ai. Functional proteomic study of selected enzymes participating in the β-ketoadipate and the Entner-Doudoroff pathways revealed a metabolic reconfiguration induced by benzoate compared to glucose. Metabolome analysis implied the existence of both catechol and protocatechuate branches of the β-ketoadipate pathway. Enzymatic study of benzoate grown cultures confirmed the activity of the ortho- catechol branch of the β-ketoadipate pathway. Even high concentrations of benzoate did not show increased stress protein synthesis, testifying to its extremophilic nature capable of surviving in harsh conditions. This ability of Pseudomonas aeruginosa san ai to efficiently degrade benzoate can provide a wide range of use of this strain in environmental and agricultural application.

Synthetic Control of Metabolic States in Pseudomonas putida by Tuning Polyhydroxyalkanoate Cycle

Polyhydroxyalkanoates (PHAs) are polyesters produced by numerous microorganisms for energy and carbon storage. Simultaneous synthesis and degradation of PHA drives a dynamic cycle linked to the central carbon metabolism, which modulates numerous and diverse bacterial processes, such as stress endurance, pathogenesis, and persistence. Here, we analyze the role of the PHA cycle in conferring robustness to the model bacterium P. putida KT2440. To assess the effect of this cycle in the cell, we began by constructing a PHA depolymerase (PhaZ) mutant strain that had its PHA cycle blocked. We then restored the flux through the cycle in the context of an engineered library of P. putida strains harboring differential levels of PhaZ. High-throughput phenotyping analyses of this collection of strains revealed significant changes in response to PHA cycle performance impacting cell number and size, PHA accumulation, and production of extracellular (R)-hydroxyalkanoic acids. To understand the metabolic changes at the system level due to PHA turnover, we contextualized these physiological data using the genome-scale metabolic model iJN1411. Model-based predictions suggest successive metabolic steady states during the growth curve and an important carbon flux rerouting driven by the activity of the PHA cycle. Overall, we demonstrate that modulating the activity of the PHA cycle gives us control over the carbon metabolism of P. putida, which in turn will give us the ability to tailor cellular mechanisms driving stress tolerance, e.g., defenses against oxidative stress, and any potential biotechnological applications. IMPORTANCE Despite large research efforts devoted to understanding the flexible metabolism of Pseudomonas beyond the role of key regulatory players, the metabolic basis powering the dynamic control of its biological fitness under disturbance conditions remains largely unknown. Among other metabolic hubs, the so-called PHA cycle, involving simultaneous synthesis and degradation of PHAs, is emerging as a pivotal metabolic trait powering metabolic robustness and resilience in this bacterial group. Here, we provide evidence suggesting that metabolic states in Pseudomonas can be anticipated, controlled, and engineered by tailoring the flux through the PHA cycle. Overall, our study suggests that the PHA cycle is a promising metabolic target toward achieving control over bacterial metabolic robustness. This is likely to open up a broad range of applications in areas as diverse as pathogenesis and biotechnology.

Glucose metabolism in Pseudomonas aeruginosa is cyclic when producing Polyhydroxyalkanoates and Rhamnolipids

Pseudomonas aeruginosa is an important chassis for production of polyhydroxyalkanoates (PHA) and rhamnolipids (RHL). Advances in the understanding of the biosynthesis metabolism of these biocompounds are crucial for increasing yield. ¹³C-Metabolic Flux Ratio Analysis (¹³C-MFA) is a technique to estimate in vivo metabolic fluxes ratios. PHA and RHL are essentially non-growth associated products of biotechnological interest and both contain hydroxyalkanoates (HAs), whose labeling patterns could be accessed by GC-MS. In this study, to reveal the relative contributions of the Entner-Doudoroff (ED) pathway and the non-oxidative Pentose Phosphate (PP) pathway to PHA and RHL production, ¹³C-MFA was performed in Pseudomonas aeruginosa LFM634 when supplied with labeled glucose. This bacterial strain lacks both functional EMP and the oxidative PP branch. Labeling patterns in HAs were measured. Experiments with [U-¹³C] glucose indicated a low flux though PP pathway. An optimal design of labeling experiment showed that [6-¹³C] glucose would be the best substrate to enable an estimation of the ED flux with high accuracy. Results of experiments performed with this isotope indicated that about two-thirds of glyceraldehyde 3-phosphate is recycled through a cyclic ED architecture, suggesting that P. aeruginosa utilizes that cycle to regulate the NADPH/Acetyl-CoA ratio for PHA and RHL biosynthesis.

The Pseudomonas aeruginosa whole genome sequence: A 20th anniversary celebration

Toward the end of August 2000, the 6.3 Mbp whole genome sequence of Pseudomonas aeruginosa strain PAO1 was published. With 5570 open reading frames (ORFs), PAO1 had the largest microbial genome sequenced up to that point in time-including a large proportion of metabolic, transport and antimicrobial resistance genes supporting its ability to colonize diverse environments. A remarkable 9% of its ORFs were predicted to encode proteins with regulatory functions, providing new insight into bacterial network complexity as a function of network size. In this celebratory article, we fast forward 20 years, and examine how access to this resource has transformed our understanding of P. aeruginosa. What follows is more than a simple review or commentary; we have specifically asked some of the leaders in the field to provide personal reflections on how the PAO1 genome sequence, along with the Pseudomonas Community Annotation Project (PseudoCAP) and Pseudomonas Genome Database (pseudomonas.com), have contributed to the many exciting discoveries in this field. In addition to bringing us all up to date with the latest developments, we also ask our contributors to speculate on how the next 20 years of Pseudomonas research might pan out.

Pseudomonas aeruginosa transcriptome adaptations from colonization to biofilm infection of skin wounds

In burn patients Pseudomonas aeruginosa infection is a major cause of morbidity. Analysis of the pathogen's gene expression as it transitions from colonization to acute and then biofilm wound infection may provide strategies for infection control. Toward this goal, we seeded log-phase P. aeruginosa (PAO1) into 3-day-old, full-thickness excision wounds (rabbit ear) and harvested the bacteria during colonization (Hrs 2 and 6), acute infection (Hr 24), and biofilm infection (Days 5 and 9) for transcriptome analysis (RNA-Seq). After 2-6 h in the wound, genes for metabolism and cell replication were down-regulated while wound-adaptation genes were up-regulated (vs. expression in log-phase culture). As the infection progressed from acute to biofilm infection, more genes became up-regulated than down-regulated, but the down-regulated genes enriched in more pathways, likely because the genes and pathways that bacteria already colonizing wounds up-regulate to establish biofilm infection are less known. Across the stages of infection, carbon-utilization pathways shifted. During acute infection, itaconate produced by myeloid cells appears to have been a carbon source because myeloid cell infiltration and the expression of the host gene, ACOD1, for itaconate production peaked coincidently with the expression of the PAO1 genes for itaconate transport and catabolism. Additionally, branched-chain amino acids are suggested to be a carbon source in acute infection and in biofilm infection. In biofilm infection, fatty acid degradation was also up-regulated. These carbon sources feed into the glyoxylate cycle that was coincidently up-regulated, suggesting it provided the precursors for P. aeruginosa to synthesize macromolecules in establishing wound infection.

Oxidative Stress Response in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative environmental and human opportunistic pathogen highly adapted to many different environmental conditions. It can cause a wide range of serious infections, including wounds, lungs, the urinary tract, and systemic infections. The high versatility and pathogenicity of this bacterium is attributed to its genomic complexity, the expression of several virulence factors, and its intrinsic resistance to various antimicrobials. However, to thrive and establish infection, P. aeruginosa must overcome several barriers. One of these barriers is the presence of oxidizing agents (e.g., hydrogen peroxide, superoxide, and hypochlorous acid) produced by the host immune system or that are commonly used as disinfectants in a variety of different environments including hospitals. These agents damage several cellular molecules and can cause cell death. Therefore, bacteria adapt to these harsh conditions by altering gene expression and eliciting several stress responses to survive under oxidative stress. Here, we used PubMed to evaluate the current knowledge on the oxidative stress responses adopted by P. aeruginosa. We will describe the genes that are often differently expressed under oxidative stress conditions, the pathways and proteins employed to sense and respond to oxidative stress, and how these changes in gene expression influence pathogenicity and the virulence of P. aeruginosa. Understanding these responses and changes in gene expression is critical to controlling bacterial pathogenicity and developing new therapeutic agents.

The three NADH dehydrogenases of Pseudomonas aeruginosa: Their roles in energy metabolism and links to virulence

Pseudomonas aeruginosa is a ubiquitous opportunistic pathogen which relies on a highly adaptable metabolism to achieve broad pathogenesis. In one example of this flexibility, to catalyze the NADH:quinone oxidoreductase step of the respiratory chain, P. aeruginosa has three different enzymes: NUO, NQR and NDH2, all of which carry out the same redox function but have different energy conservation and ion transport properties. In order to better understand the roles of these enzymes, we constructed two series of mutants: (i) three single deletion mutants, each of which lacks one NADH dehydrogenase and (ii) three double deletion mutants, each of which retains only one of the three enzymes. All of the mutants grew approximately as well as wild type, when tested in rich and minimal medium and in a range of pH and [Na+] conditions, except that the strain with only NUO (ΔnqrFΔndh) has an extended lag phase. During exponential phase, the NADH dehydrogenases contribute to total wild-type activity in the following order: NQR > NDH2 > NUO. Some mutants, including the strain without NQR (ΔnqrF) had increased biofilm formation, pyocyanin production, and killed more efficiently in both macrophage and mouse infection models. Consistent with this, ΔnqrF showed increased transcription of genes involved in pyocyanin production.

Pseudomonas aeruginosa reverse diauxie is a multidimensional, optimized, resource utilization strategy

Pseudomonas aeruginosa is a globally-distributed bacterium often found in medical infections. The opportunistic pathogen uses a different, carbon catabolite repression (CCR) strategy than many, model microorganisms. It does not utilize a classic diauxie phenotype, nor does it follow common systems biology assumptions including preferential consumption of glucose with an 'overflow' metabolism. Despite these contradictions, P. aeruginosa is competitive in many, disparate environments underscoring knowledge gaps in microbial ecology and systems biology. Physiological, omics, and in silico analyses were used to quantify the P. aeruginosa CCR strategy known as 'reverse diauxie'. An ecological basis of reverse diauxie was identified using a genome-scale, metabolic model interrogated with in vitro omics data. Reverse diauxie preference for lower energy, nonfermentable carbon sources, such as acetate or succinate over glucose, was predicted using a multidimensional strategy which minimized resource investment into central metabolism while completely oxidizing substrates. Application of a common, in silico optimization criterion, which maximizes growth rate, did not predict the reverse diauxie phenotypes. This study quantifies P. aeruginosa metabolic strategies foundational to its wide distribution and virulence including its potentially, mutualistic interactions with microorganisms found commonly in the environment and in medical infections.

Bacterial metabolism-inspired molecules to modulate antibiotic efficacy

The decreasing antibiotic susceptibility of bacterial pathogens calls for novel antimicrobial therapies. Traditional screening pathways based on drug-target interaction have gradually reached the stage of diminishing returns. Thus, novel strategies are urgently needed in the fight against antibiotic-refractory bacteria, particularly for tolerant bacteria. Recently, evidence has accumulated demonstrating that microbial changes caused by bacterial metabolic processes significantly modulate antibiotic killing. A better understanding of these bacterial metabolic processes is indicating a need to screen novel metabolic modulators as potential antibiotic adjuvants. In this review, we describe the state of our current knowledge about how these bacterial metabolism-inspired molecules affect antibiotic efficacy, including potentiation and inhibition activity. In addition, the challenges faced and prospects for bringing them into clinic are also discussed. These examples may provide candidates or targets for the development of novel antibiotic adjuvants.

Wax ester production in nitrogen-rich conditions by metabolically engineered Acinetobacter baylyi ADP1

Metabolic engineering can be used as a powerful tool to redirect cell resources towards product synthesis, also in conditions that are not optimal for the production. An example of synthesis strongly dependent on external conditions is the production of storage lipids, which typically requires a high carbon/nitrogen ratio. This requirement also limits the use of abundant nitrogen-rich materials, such as industrial protein by-products, as substrates for lipid production. Acinetobacter baylyi ADP1 is known for its ability to produce industrially interesting storage lipids, namely wax esters (WEs). Here, we engineered A. baylyi ADP1 by deleting the gene aceA encoding for isocitrate lyase and overexpressing fatty acyl-CoA reductase Acr1 in the wax ester production pathway to allow redirection of carbon towards WEs. This strategy led to 3-fold improvement in yield (0.075 ​g/g glucose) and 3.15-fold improvement in titer (1.82 ​g/L) and productivity (0.038 ​g/L/h) by a simple one-stage batch cultivation with glucose as carbon source. The engineered strain accumulated up to 27% WEs of cell dry weight. The titer and cellular WE content are the highest reported to date among microbes. We further showed that the engineering strategy alleviated the inherent requirement for high carbon/nitrogen ratio and demonstrated the production of wax esters using nitrogen-rich substrates including casamino acids, yeast extract, and baker's yeast hydrolysate, which support biomass production but not WE production in wild-type cells. The study demonstrates the power of metabolic engineering in overcoming natural limitations in the production of storage lipids.

The aerobic respiratory chain of Pseudomonas aeruginosa cultured in artificial urine media: Role of NQR and terminal oxidases

Pseudomonas aeruginosa is a Gram-negative γ-proteobacterium that forms part of the normal human microbiota and it is also an opportunistic pathogen, responsible for 30% of all nosocomial urinary tract infections. P. aeruginosa carries a highly branched respiratory chain that allows the colonization of many environments, such as the urinary tract, catheters and other medical devices. P. aeruginosa respiratory chain contains three different NADH dehydrogenases (complex I, NQR and NDH-2), whose physiologic roles have not been elucidated, and up to five terminal oxidases: three cytochrome c oxidases (COx), a cytochrome bo₃ oxidase (CYO) and a cyanide-insensitive cytochrome bd-like oxidase (CIO). In this work, we studied the composition of the respiratory chain of P. aeruginosa cells cultured in Luria Broth (LB) and modified artificial urine media (mAUM), to understand the metabolic adaptations of this microorganism to the growth in urine. Our results show that the COx oxidases play major roles in mAUM, while P. aeruginosa relies on CYO when growing in LB medium. Moreover, our data demonstrate that the proton-pumping NQR complex is the main NADH dehydrogenase in both LB and mAUM. This enzyme is resistant to HQNO, an inhibitory molecule produced by P. aeruginosa, and may provide an advantage against the natural antibacterial agents produced by this organism. This work offers a clear picture of the composition of this pathogen’s aerobic respiratory chain and the main roles that NQR and terminal oxidases play in urine, which is essential to understand its physiology and could be used to develop new antibiotics against this notorious multidrug-resistant microorganism.

Metabolic manipulation by Pseudomonas fluorescens: a powerful stratagem against oxidative and metal stress

Metabolism is the foundation of all living organisms and is at the core of numerous if not all biological processes. The ability of an organism to modulate its metabolism is a central characteristic needed to proliferate, to be dormant and to survive any assault. Pseudomonas fluorescens is bestowed with a uniquely versatile metabolic framework that enables the microbe to adapt to a wide range of conditions including disparate nutrients and toxins. In this mini-review we elaborate on the various metabolic reconfigurations evoked by this microbial system to combat reactive oxygen/nitrogen species and metal stress. The fine-tuning of the NADH/NADPH homeostasis coupled with the production of α-keto-acids and ATP allows for the maintenance of a reductive intracellular milieu. The metabolic networks propelling the synthesis of metabolites like oxalate and aspartate are critical to keep toxic metals at bay. The biochemical processes resulting from these defensive mechanisms provide molecular clues to thwart infectious microbes and reveal elegant pathways to generate value-added products.

Contextual Flexibility in Pseudomonas aeruginosa Central Carbon Metabolism during Growth in Single Carbon Sources

Pseudomonas aeruginosa is an opportunistic human pathogen that is well known for causing infections in the airways of people with cystic fibrosis. Although it is clear that P. aeruginosa is metabolically well adapted to life in the CF lung, little is currently known about how the organism metabolizes the nutrients available in the airways. In this work, we used a combination of gene expression and isotope tracer (“fluxomic”) analyses to find out exactly where the input carbon goes during growth on two CF-relevant carbon sources, acetate and glycerol (derived from the breakdown of lung surfactant). We found that carbon is routed (“fluxed”) through very different pathways during growth on these substrates and that this is accompanied by an unexpected remodeling of the cell’s electron transfer pathways. Having access to this “blueprint” is important because the metabolism of P. aeruginosa is increasingly being recognized as a target for the development of much-needed antimicrobial agents.

Pseudomonas aeruginosa is an opportunistic human pathogen, particularly noted for causing infections in the lungs of people with cystic fibrosis (CF). Previous studies have shown that the gene expression profile of P. aeruginosa appears to converge toward a common metabolic program as the organism adapts to the CF airway environment. However, we still have only a limited understanding of how these transcriptional changes impact metabolic flux at the systems level. To address this, we analyzed the transcriptome, proteome, and fluxome of P. aeruginosa grown on glycerol or acetate. These carbon sources were chosen because they are the primary breakdown products of an airway surfactant, phosphatidylcholine, which is known to be a major carbon source for P. aeruginosa in CF airways. We show that the fluxes of carbon throughout central metabolism are radically different among carbon sources. For example, the newly recognized “EDEMP cycle” (which incorporates elements of the Entner-Doudoroff [ED] pathway, the Embden-Meyerhof-Parnas [EMP] pathway, and the pentose phosphate [PP] pathway) plays an important role in supplying NADPH during growth on glycerol. In contrast, the EDEMP cycle is attenuated during growth on acetate, and instead, NADPH is primarily supplied by the reaction catalyzed by isocitrate dehydrogenase(s). Perhaps more importantly, our proteomic and transcriptomic analyses revealed a global remodeling of gene expression during growth on the different carbon sources, with unanticipated impacts on aerobic denitrification, electron transport chain architecture, and the redox economy of the cell. Collectively, these data highlight the remarkable metabolic plasticity of P. aeruginosa; that plasticity allows the organism to seamlessly segue between different carbon sources, maximizing the energetic yield from each.

Pseudomonad reverse carbon catabolite repression, interspecies metabolite exchange, and consortial division of labor

Microorganisms acquire energy and nutrients from dynamic environments, where substrates vary in both type and abundance. The regulatory system responsible for prioritizing preferred substrates is known as carbon catabolite repression (CCR). Two broad classes of CCR have been documented in the literature. The best described CCR strategy, referred to here as classic CCR (cCCR), has been experimentally and theoretically studied using model organisms such as Escherichia coli. cCCR phenotypes are often used to generalize universal strategies for fitness, sometimes incorrectly. For instance, extremely competitive microorganisms, such as Pseudomonads, which arguably have broader global distributions than E. coli, have achieved their success using metabolic strategies that are nearly opposite of cCCR. These organisms utilize a CCR strategy termed 'reverse CCR' (rCCR), because the order of preferred substrates is nearly reverse that of cCCR. rCCR phenotypes prefer organic acids over glucose, may or may not select preferred substrates to optimize growth rates, and do not allocate intracellular resources in a manner that produces an overflow metabolism. cCCR and rCCR have traditionally been interpreted from the perspective of monocultures, even though most microorganisms live in consortia. Here, we review the basic tenets of the two CCR strategies and consider these phenotypes from the perspective of resource acquisition in consortia, a scenario that surely influenced the evolution of cCCR and rCCR. For instance, cCCR and rCCR metabolism are near mirror images of each other; when considered from a consortium basis, the complementary properties of the two strategies can mitigate direct competition for energy and nutrients and instead establish cooperative division of labor.

High-quality genome-scale metabolic modelling of Pseudomonas putida highlights its broad metabolic capabilities

Genome-scale reconstructions of metabolism are computational species-specific knowledge bases able to compute systemic metabolic properties. We present a comprehensive and validated reconstruction of the biotechnologically relevant bacterium Pseudomonas putida KT2440 that greatly expands computable predictions of its metabolic states. The reconstruction represents a significant reactome expansion over available reconstructed bacterial metabolic networks. Specifically, iJN1462 (i) incorporates several hundred additional genes and associated reactions resulting in new predictive capabilities, including new nutrients supporting growth; (ii) was validated by in vivo growth screens that included previously untested carbon (48) and nitrogen (41) sources; (iii) yielded gene essentiality predictions showing large accuracy when compared with a knock-out library and Bar-seq data; and (iv) allowed mapping of its network to 82 P. putida sequenced strains revealing functional core that reflect the large metabolic versatility of this species, including aromatic compounds derived from lignin. Thus, this study provides a thoroughly updated metabolic reconstruction and new computable phenotypes for P. putida, which can be leveraged as a first step toward understanding the pan metabolic capabilities of Pseudomonas.

Early Detection of Cystic Fibrosis Acute Pulmonary Exacerbations by Exhaled Breath Condensate Metabolomics

The most common cause of death in cystic fibrosis (CF) patients is progressive lung function decline, which is punctuated by acute pulmonary exacerbations (APEs). A major challenge is to discover biomarkers for detecting an oncoming APE and allow for pre-emptive clinical interventions. Metabolic profiling of exhaled breath condensate (EBC) samples collected from CF patients before, during, and after APEs and under stable conditions (n = 210) was performed using ultraperformance liquid chromatography (UPLC) coupled to Orbitrap mass spectrometry (MS). Negative ion mode MS data showed that classification between metabolic profiles from "pre-APE" (pending APE before the CF patient had any signs of illness) and stable CF samples was possible with good sensitivities (85.7 and 89.5%), specificities (88.4 and 84.1%), and accuracies (87.7 and 85.7%) for pediatric and adult patients, respectively. Improved classification performance was achieved by combining positive with negative ion mode data. Discriminant metabolites included two potential biomarkers identified in a previous pilot study: lactic acid and 4-hydroxycyclohexylcarboxylic acid. Some of the discriminant metabolites had microbial origins, indicating a possible role of bacterial metabolism in APE progression. The results show promise for detecting an oncoming APE using EBC metabolites, thus permitting early intervention to abort such an event.

Reverse diauxie phenotype in Pseudomonas aeruginosa biofilm revealed by exometabolomics and label-free proteomics

Microorganisms enhance fitness by prioritizing catabolism of available carbon sources using a process known as carbon catabolite repression (CCR). Planktonically grown Pseudomonas aeruginosa is known to prioritize the consumption of organic acids including lactic acid over catabolism of glucose using a CCR strategy termed "reverse diauxie." P. aeruginosa is an opportunistic pathogen with well-documented biofilm phenotypes that are distinct from its planktonic phenotypes. Reverse diauxie has been described in planktonic cultures, but it has not been documented explicitly in P. aeruginosa biofilms. Here a combination of exometabolomics and label-free proteomics was used to analyze planktonic and biofilm phenotypes for reverse diauxie. P. aeruginosa biofilm cultures preferentially consumed lactic acid over glucose, and in addition, the cultures catabolized the substrates completely and did not exhibit the acetate secreting "overflow" metabolism that is typical of many model microorganisms. The biofilm phenotype was enabled by changes in protein abundances, including lactate dehydrogenase, fumarate hydratase, GTP cyclohydrolase, L-ornithine N(5)-monooxygenase, and superoxide dismutase. These results are noteworthy because reverse diauxie-mediated catabolism of organic acids necessitates a terminal electron acceptor like O2, which is typically in low supply in biofilms due to diffusion limitation. Label-free proteomics identified dozens of proteins associated with biofilm formation including 16 that have not been previously reported, highlighting both the advantages of the methodology utilized here and the complexity of the proteomic adaptation for P. aeruginosa biofilms. Documenting the reverse diauxic phenotype in P. aeruginosa biofilms is foundational for understanding cellular nutrient and energy fluxes, which ultimately control growth and virulence.

Bio-detoxification Bacteria Isolated from Dye-Polluted Soils Promote Lactic Acid Production from Ammonia Pretreated Corn Stover

Agro-stovers are the most abundant substrates for producing lactic acid, which has great potential application in the production of biodegradable and biocompatible polylactic acid polymers. However, chemical pretreatments on agro-stovers generate inhibitors that repress the subsequent lactic acid fermentation. In this study, three bacterial strains (Enterococcus faecalis B101, Acinetobacter calcoaceticus C1, and Pseudomonas aeruginosa CS) isolated from dye-polluted soils could utilize phenolic inhibitor mimics (vanillin, 4- hydroxybenzaldehyde, or syringaldehyde) from alkaline pretreated corn stovers as a sole carbon source. Lactic acid titer increased from 27.42 g/L (Bacillus coagulans LA204 alone) to 44.76 g/L (CS and LA204) using 50 g/L glucose with 1 g/L 4-hydroxybenzaldehyde added. Lactic acid production from 50 g/L ammonia pretreated corn stover was increased nearly twofold by inoculating phenolic degradation bacteria and lactic acid bacteria (C1& Lactobacillus pentosus FL0421). In the control (FL0421 alone), only 16.98 g/L of lactic acid was produced. The isolated and identified strains degraded the phenolic compounds and increased the lactic acid production from glucose and ammonia pretreated corn stover. These characteristics of the strains support industrial application with efficient in situ detoxification of phenolic compounds during lactic acid production from agro-stovers using simultaneous saccharification and fermentation (SSF).

A Cyclic Metabolic Network in Pseudomonas protegens Pf-5 Prioritizes the Entner-Doudoroff Pathway and Exhibits Substrate Hierarchy during Carbohydrate Co-Utilization

The genetic characterization of Pseudomonas protegens Pf-5 was recently completed. However, the inferred metabolic network structure has not yet been evaluated experimentally. Here, we employed 13C-tracers and quantitative flux analysis to investigate the intracellular network for carbohydrate metabolism. In lieu of the direct phosphorylation of glucose by glucose kinase, glucose catabolism was characterized primarily by the oxidation of glucose to gluconate and 2-ketogluconate before the phosphorylation of these metabolites to feed the Entner-Doudoroff (ED) pathway. In the absence of phosphofructokinase activity, a cyclic flux from the ED pathway to the upper Embden-Meyerhof-Parnas (EMP) pathway was responsible for routing glucose-derived carbons to the non-oxidative pentose phosphate (PP) pathway. Consistent with the lack of annotated genes in P. protegens Pf-5 for the transport or initial catabolism of pentoses and galactose, only glucose was assimilated into intracellular metabolites in the presence of xylose, arabinose, or galactose. However, when glucose was fed simultaneously with fructose or mannose, co-uptake of these hexoses was evident, but glucose was preferred over fructose (3 to 1) and over mannose (4 to 1). Despite gene annotation of mannose catabolism to fructose-6-phosphate, metabolite labeling patterns revealed that mannose was assimilated into fructose-1,6-bisphosphate, similarly to fructose catabolism. Remarkably, carbons from mannose and fructose were also found to cycle backward through the upper EMP pathway toward the ED pathway. Therefore, the operational metabolic network for processing carbohydrates in P. protegens Pf-5 prioritizes flux through the ED pathway to channel carbons to EMP, PP, and downstream pathways.IMPORTANCE Species of the Pseudomonas genus thrive in various nutritional environments and have strong biocatalytic potential due to their diverse metabolic capabilities. Carbohydrate substrates are ubiquitous both in environmental matrices and in feedstocks for engineered bioconversion. Here, we investigated the metabolic network for carbohydrate metabolism in Pseudomonas protegens Pf-5. Metabolic flux quantitation revealed the relative involvement of different catabolic routes in channeling carbohydrate carbons through a cyclic metabolic network. We also uncovered that mannose catabolism was similar to fructose catabolism, despite the annotation of a different pathway in the genome. Elucidation of the constitutive metabolic network in P. protegens is important for understanding its innate carbohydrate processing, thus laying the foundation for targeting metabolic engineering of this untapped Pseudomonas species.

Regulation of glycolytic flux and overflow metabolism depending on the source of energy generation for energy demand

Overflow metabolism is a common phenomenon observed at higher glycolytic flux in many bacteria, yeast (known as Crabtree effect), and mammalian cells including cancer cells (known as Warburg effect). This phenomenon has recently been characterized as the trade-offs between protein costs and enzyme efficiencies based on coarse-graining approaches. Moreover, it has been recognized that the glycolytic flux increases as the source of energy generation changes from energetically efficient respiration to inefficient respiro-fermentative or fermentative metabolism causing overflow metabolism. It is highly desired to clarify the metabolic regulation mechanisms behind such phenomena. Metabolic fluxes are located on top of the hierarchical regulation systems, and represent the outcome of the integrated response of all levels of cellular regulation systems. In the present article, we discuss about the different levels of regulation systems for the modulation of fluxes depending on the growth rate, growth condition such as oxygen limitation that alters the metabolism towards fermentation, and genetic perturbation affecting the source of energy generation from respiration to respiro-fermentative metabolism in relation to overflow metabolism. The intracellular metabolite of the upper glycolysis such as fructose 1,6-bisphosphate (FBP) plays an important role not only for flux sensing, but also for the regulation of the respiratory activity either directly or indirectly (via transcription factors) at higher growth rate. The glycolytic flux regulation is backed up (enhanced) by unphosphorylated EIIA and HPr of the phosphotransferase system (PTS) components, together with the sugar-phosphate stress regulation, where the transcriptional regulation is further modulated by post-transcriptional regulation via the degradation of mRNA (stability of mRNA) in Escherichia coli. Moreover, the channeling may also play some role in modulating the glycolytic cascade reactions.

Bacterial Metabolism During Biofilm Growth Investigated by 13C Tracing

This study investigated the metabolism of Pseudomonas aeruginosa PAO1 during its biofilm development via microscopy imaging, gene expression analysis, and ¹³C-labeling. First, dynamic labeling was employed to investigate glucose utilization rate in fresh biofilms (thickness 40∼60 micrometer). The labeling turnover time of glucose-6-P indicated biofilm metabolism was substantially slower than planktonic cells. Second, PAO1 was cultured in continuous tubular biofilm reactors or shake flasks. Then ¹³C-metabolic flux analysis of PAO1 was performed based on the isotopomer patterns of proteinogenic amino acids. The results showed that PAO1 biofilm cells during growth conserved the flux features as their planktonic mode. (1) Glucose could be degraded by two cyclic routes (the TCA cycle and the Entner-Doudoroff-Embden-Meyerhof-Parnas loop) that facilitated NAD(P)H supplies. (2) Anaplerotic pathways (including pyruvate shunt) increased flux plasticity. (3) Biofilm growth phenotype did not require significant intracellular flux rewiring (variations between biofilm and planktonic flux network, normalized by glucose uptake rate as 100%, were less than 20%). (4) Transcription analysis indicated that key catabolic genes in fresh biofilm cells had expression levels comparable to planktonic cells. Finally, PAO1, Shewanella oneidensis (as the comparing group), and their c-di-GMP transconjugants (with different biofilm formation capabilities) were ¹³C-labeled under biofilm reactors or planktonic conditions. Analysis of amino acid labeling variances from different cultures indicated Shewanella flux network was more flexibly changed than PAO1 during its biofilm formation.

From General to Specific: Can Pseudomonas Primary Metabolism Be Exploited for Narrow-Spectrum Antibiotics?

The spread of antimicrobial resistance is a major threat to human health, and patients requiring prolonged antibiotic exposure are in desperate need of new therapeutic strategies. It has been hypothesized that tailoring our antibiotics to inhibit molecular targets specific to pathogens might stem the spread of resistance. A prime candidate for such a strategy is Pseudomonas aeruginosa, which can be found in the lungs of nearly all adult cystic fibrosis patients and, due to chronic exposure to antibiotics, has a high rate of multidrug-resistant strains. Although much research has been done on P. aeruginosa virulence factors as narrow-spectrum targets, less attention has been paid to primary carbon metabolism being leveraged for pathogen-specific mechanisms. However, early studies show that primary metabolic pathways, although shared amongst all organisms, contain intricacies specific to Pseudomonas species that have potential for antibiotic exploitation. Here we lay out some of this work in the hopes that it inspires researchers to continue developing a knowledge base for future antibiotic discovery to build upon and include a case study of a Pseudomonas primary metabolic pathway that has been targeted by small molecules in a species-specific manner.

Gluconeogenic precursor availability regulates flux through the glyoxylate shunt in Pseudomonas aeruginosa

The glyoxylate shunt bypasses the oxidative decarboxylation steps of the tricarboxylic acid (TCA) cycle, thereby conserving carbon skeletons for gluconeogenesis and biomass production. In Escherichia coli, carbon flux is redirected through the first enzyme of the glyoxylate shunt, isocitrate lyase (ICL), following phosphorylation and inactivation of the TCA cycle enzyme, isocitrate dehydrogenase (ICD), by the kinase/phosphatase, AceK. In contrast, mycobacterial species lack AceK and employ a phosphorylation-insensitive isocitrate dehydrogenase (IDH), which is allosterically activated by the product of ICL activity, glyoxylate. However, Pseudomonas aeruginosa expresses IDH, ICD, ICL, and AceK, raising the question of how these enzymes are regulated to ensure proper flux distribution between the competing pathways. Here, we present the structure, kinetics, and regulation of ICL, IDH, and ICD from P. aeruginosa We found that flux partitioning is coordinated through reciprocal regulation of these enzymes, linking distribution of carbon flux to the availability of the key gluconeogenic precursors, oxaloacetate and pyruvate. Specifically, a greater abundance of these metabolites activated IDH and inhibited ICL, leading to increased TCA cycle flux. Regulation was also exerted through AceK-dependent phosphorylation of ICD; high levels of acetyl-CoA (which would be expected to accumulate when oxaloacetate is limiting) stimulated the kinase activity of AceK, whereas high levels of oxaloacetate stimulated its phosphatase activity. In summary, the TCA cycle-glyoxylate shunt branch point in P. aeruginosa has a complex enzymology that is profoundly different from those in other species characterized to date. Presumably, this reflects its predilection for consuming fatty acids, especially during infection scenarios.

Transcriptome analysis of a Pseudomonas aeruginosasn-glycerol-3-phosphate dehydrogenase mutant reveals a disruption in bioenergetics

Pseudomonas aeruginosa causes acute and chronic human infections and is the major cause of morbidity and mortality in cystic fibrosis (CF) patients. We previously determined that the sn-glycerol-3-phosphate dehydrogenase encoded by glpD plays a larger role in P. aeruginosa physiology beyond its role in glycerol metabolism. To better understand the effect of a glpD mutation on P. aeruginosa physiology we compared the transcriptomes of P. aeruginosa strain PAO1 and the PAO1ΔglpD mutant using RNA-seq analysis. We determined that a null mutation of glpD significantly altered amino acid metabolism in P. aeruginosa and affected the production of intermediates that are channelled into the tricarboxylic acid cycle. Moreover, the loss of glpD induced a general stress response mediated by RpoS in P. aeruginosa. Several other phenotypes observed for the P. aeruginosa glpD mutant include increased persister cell formation, reduced extracellular ATP accumulation and increased heat output. Taken together, these findings implicate sn-glycerol-3-phosphate dehydrogenase as a key player in energy metabolism in P. aeruginosa.

Computational Systems Biology of Metabolism in Infection

A systems approach to elucidate the effect of infection on cell metabolism provides several opportunities from a better understanding of molecular mechanisms to the identification of potential biomarkers and drug targets. This is obvious from the fact that we have witnessed the accelerated use of computational systems biology in the last five years to study metabolic changes in pathogen and/or host cells in response to infection. In this chapter, we aim to present a comprehensive review of the recent research by focusing on genome-scale metabolic network models of pathogen-host systems and genome-wide metabolomics and fluxomics analysis of infected cells.

Cloning, overexpression, and purification of glucose-6-phosphate dehydrogenase of Pseudomonas aeruginosa

Glucose-6-phosphate dehydrogenase (G6PDH) (EC 1.1.1.363) plays an important role in the human pathogen Pseudomonas aeruginosa because it generates NADPH, an essential cofactor for several biosynthetic pathways and antioxidant enzymes. P. aeruginosa G6PDH is also a key enzyme in the metabolism of various carbon sources, such as glucose, glycerol, fructose, and mannitol. Understanding the kinetic characteristics and mechanisms that control the activity of this enzyme is crucial for future studies in this context. However, one of the impediments to achieving this goal is the limited amount of protein obtained when current purification protocols are implemented, a factor curtailing its biochemical characterization. In this study, we report a fast, efficient and reproducible procedure for the purification of P. aeruginosa G6PDH that can be implemented in a short period (2 days). In order to establish this protocol, the zwf gene, which encodes for this enzyme, was cloned and overexpressed in Escherichia coli cells. In contrast to other procedures, our method is based on protein precipitation with CaCl₂ and further purification by ion exchange chromatography. Using this protocol, we were able to obtain 31 mg/L of pure protein that manifested specific activity of 145.7 U/mg. The recombinant enzyme obtained in this study manifested similar physicochemical and kinetic properties to those reported in previous works for this molecule. The large quantities of active enzyme obtained using this procedure will facilitate its structural characterization and identify differences between P. aeruginosa- and human G6PDH, thus contributing to the search for selective inhibitors against the bacterial enzyme.

Refactoring the Embden-Meyerhof-Parnas Pathway as a Whole of Portable GlucoBricks for Implantation of Glycolytic Modules in Gram-Negative Bacteria

The Embden–Meyerhof–Parnas (EMP) pathway is generally considered to be the biochemical standard for glucose catabolism. Alas, its native genomic organization and the control of gene expression in Escherichia coli are both very intricate, which limits the portability of the EMP pathway to other biotechnologically important bacterial hosts that lack the route. In this work, the genes encoding all the enzymes of the linear EMP route have been individually recruited from the genome of E. coli K-12, edited in silico to remove their endogenous regulatory signals, and synthesized de novo following a standard (GlucoBrick) that enables their grouping in the form of functional modules at the user’s will. After verifying their activity in several glycolytic mutants of E. coli, the versatility of these GlucoBricks was demonstrated in quantitative physiology tests and biochemical assays carried out in Pseudomonas putida KT2440 and P. aeruginosa PAO1 as the heterologous hosts. Specific configurations of GlucoBricks were also adopted to streamline the downward circulation of carbon from hexoses to pyruvate in E. coli recombinants, thereby resulting in a 3-fold increase of poly(3-hydroxybutyrate) synthesis from glucose. Refactoring whole metabolic blocks in the fashion described in this work thus eases the engineering of biochemical processes where the optimization of carbon traffic is facilitated by the operation of the EMP pathway—which yields more ATP than other glycolytic routes such as the Entner–Doudoroff pathway.

Carbon Sources Tune Antibiotic Susceptibility in Pseudomonas aeruginosa via Tricarboxylic Acid Cycle Control

Metabolically dormant bacteria present a critical challenge to effective antimicrobial therapy because these bacteria are genetically susceptible to antibiotic treatment but phenotypically tolerant. Such tolerance has been attributed to impaired drug uptake, which can be reversed by metabolic stimulation. Here, we evaluate the effects of central carbon metabolite stimulations on aminoglycoside sensitivity in the pathogen Pseudomonas aeruginosa. We identify fumarate as a tobramycin potentiator that activates cellular respiration and generates a proton motive force by stimulating the tricarboxylic acid (TCA) cycle. In contrast, we find that glyoxylate induces phenotypic tolerance by inhibiting cellular respiration with acetyl-coenzyme A diversion through the glyoxylate shunt, despite drug import. Collectively, this work demonstrates that TCA cycle activity is important for both aminoglycoside uptake and downstream lethality and identifies a potential strategy for potentiating aminoglycoside treatment of P. aeruginosa infections.

Recent advances in high-throughput 13C-fluxomics

The rise of high throughput (HT) strain engineering tools accompanying the area of synthetic biology is supporting the generation of a large number of microbial cell factories. A current bottleneck in process development is our limited capacity to rapidly analyze the metabolic state of the engineered strains, and in particular their intracellular fluxes. HT ¹³C-fluxomics workflows have not yet become commonplace, despite the existence of several HT tools at each of the required stages. This includes cultivation and sampling systems, analytics for isotopic analysis, and software for data processing and flux calculation. Here, we review recent advances in the field and highlight bottlenecks that must be overcome to allow the emergence of true HT ¹³C-fluxomics workflows.

Metabolic flux analyses of Pseudomonas aeruginosa cystic fibrosis isolates

Pseudomonas aeruginosa is a metabolically versatile wide-ranging opportunistic pathogen. In humans P. aeruginosa causes infections of the skin, urinary tract, blood, and the lungs of Cystic Fibrosis patients. In addition, P. aeruginosa's broad environmental distribution, relatedness to biotechnologically useful species, and ability to form biofilms have made it the focus of considerable interest. We used ¹³C metabolic flux analysis (MFA) and flux balance analysis to understand energy and redox production and consumption and to explore the metabolic phenotypes of one reference strain and five strains isolated from the lungs of cystic fibrosis patients. Our results highlight the importance of the oxidative pentose phosphate and Entner-Doudoroff pathways in P. aeruginosa growth. Among clinical strains we report two divergent metabolic strategies and identify changes between genetically related strains that have emerged during a chronic infection of the same patient. MFA revealed that the magnitude of fluxes through the glyoxylate cycle correlates with growth rates.

Metabolic modeling of a chronic wound biofilm consortium predicts spatial partitioning of bacterial species

Background

Chronic wounds are often colonized by consortia comprised of different bacterial species growing as biofilms on a complex mixture of wound exudate. Bacteria growing in biofilms exhibit phenotypes distinct from planktonic growth, often rendering the application of antibacterial compounds ineffective. Computational modeling represents a complementary tool to experimentation for generating fundamental knowledge and developing more effective treatment strategies for chronic wound biofilm consortia.

Results

We developed spatiotemporal models to investigate the multispecies metabolism of a biofilm consortium comprised of two common chronic wound isolates: the aerobe Pseudomonas aeruginosa and the facultative anaerobe Staphylococcus aureus. By combining genome-scale metabolic reconstructions with partial differential equations for metabolite diffusion, the models were able to provide both temporal and spatial predictions with genome-scale resolution. The models were used to analyze the metabolic differences between single species and two species biofilms and to demonstrate the tendency of the two bacteria to spatially partition in the multispecies biofilm as observed experimentally. Nutrient gradients imposed by supplying glucose at the bottom and oxygen at the top of the biofilm induced spatial partitioning of the two species, with S. aureus most concentrated in the anaerobic region and P. aeruginosa present only in the aerobic region. The two species system was predicted to support a maximum biofilm thickness much greater than P. aeruginosa alone but slightly less than S. aureus alone, suggesting an antagonistic metabolic effect of P. aeruginosa on S. aureus. When each species was allowed to enhance its growth through consumption of secreted metabolic byproducts assuming identical uptake kinetics, the competitiveness of P. aeruginosa was further reduced due primarily to the more efficient lactate metabolism of S. aureus. Lysis of S. aureus by a small molecule inhibitor secreted from P. aeruginosa and/or P. aeruginosa aerotaxis were predicted to substantially increase P. aeruginosa competitiveness in the aerobic region, consistent with in vitro experimental studies.

Conclusions

Our biofilm modeling approach allows the prediction of individual species metabolism and interspecies interactions in both time and space with genome-scale resolution. This study yielded new insights into the multispecies metabolism of a chronic wound biofilm, in particular metabolic factors that may lead to spatial partitioning of the two bacterial species. We believe that P. aeruginosa lysis of S. aureus combined with nutrient competition is a particularly relevant scenario for which model predictions could be tested experimentally.

Electronic supplementary material

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

Towards stable kinetics of large metabolic networks: Nonequilibrium potential function approach

While the biochemistry of metabolism in many organisms is well studied, details of the metabolic dynamics are not fully explored yet. Acquiring adequate in vivo kinetic parameters experimentally has always been an obstacle. Unless the parameters of a vast number of enzyme-catalyzed reactions happened to fall into very special ranges, a kinetic model for a large metabolic network would fail to reach a steady state. In this work we show that a stable metabolic network can be systematically established via a biologically motivated regulatory process. The regulation is constructed in terms of a potential landscape description of stochastic and nongradient systems. The constructed process draws enzymatic parameters towards stable metabolism by reducing the change in the Lyapunov function tied to the stochastic fluctuations. Biologically it can be viewed as interplay between the flux balance and the spread of workloads on the network. Our approach allows further constraints such as thermodynamics and optimal efficiency. We choose the central metabolism of Methylobacterium extorquens AM1 as a case study to demonstrate the effectiveness of the approach. Growth efficiency on carbon conversion rate versus cell viability and futile cycles is investigated in depth.

Rapid Prediction of Bacterial Heterotrophic Fluxomics Using Machine Learning and Constraint Programming

13C metabolic flux analysis (13C-MFA) has been widely used to measure in vivo enzyme reaction rates (i.e., metabolic flux) in microorganisms. Mining the relationship between environmental and genetic factors and metabolic fluxes hidden in existing fluxomic data will lead to predictive models that can significantly accelerate flux quantification. In this paper, we present a web-based platform MFlux (http://mflux.org) that predicts the bacterial central metabolism via machine learning, leveraging data from approximately 100 13C-MFA papers on heterotrophic bacterial metabolisms. Three machine learning methods, namely Support Vector Machine (SVM), k-Nearest Neighbors (k-NN), and Decision Tree, were employed to study the sophisticated relationship between influential factors and metabolic fluxes. We performed a grid search of the best parameter set for each algorithm and verified their performance through 10-fold cross validations. SVM yields the highest accuracy among all three algorithms. Further, we employed quadratic programming to adjust flux profiles to satisfy stoichiometric constraints. Multiple case studies have shown that MFlux can reasonably predict fluxomes as a function of bacterial species, substrate types, growth rate, oxygen conditions, and cultivation methods. Due to the interest of studying model organism under particular carbon sources, bias of fluxome in the dataset may limit the applicability of machine learning models. This problem can be resolved after more papers on 13C-MFA are published for non-model species.

Proteome and carbon flux analysis of Pseudomonas aeruginosa clinical isolates from different infection sites

Pseudomonas aeruginosa is known as opportunistic pathogen frequently isolated from different infection sites. To investigate the expression rates of P. aeruginosa proteins commonly expressed by different clinical isolates, absolute protein quantities were determined employing a gel-free and data-independent LC-IMS(E) approach. Moreover, the metabolic diversity of these isolates was investigated by (13) C-metabolic flux analyses. 812 proteins were reproducibly identified and absolutely quantified for the reference strain P. aeruginosa PAO1, 363 of which were also identified and relatively quantified in all isolates. Whilst the majority of these proteins were expressed in constant amounts, expression rates of 42 proteins were highly variable between the isolates. Notably, the outer membrane protein OprH and the response regulator PhoP were strongly expressed in burned wounds isolates compared to lung/urinary tract isolates. Moreover, proteins involved in iron/amino acids uptake were found to be highly abundant in urinary tract isolates. The fluxome data revealed a conserved glycolysis, and a niche-specific divergence in fluxes through the glyoxylate shunt and the TCA cycle among the isolates. The integrated proteome/fluxome analysis did not indicate straightforward correlation between the protein amount and flux, but rather points to additional layers of regulation that mediate metabolic adaption of P. aeruginosa to different host environments. All MS data have been deposited in the ProteomeXchange with identifier PXD002373 (http://proteomecentral.proteomexchange.org/dataset/PXD002373).

The Basics of Bacteriuria: Strategies of Microbes for Persistence in Urine

Bacteriuria, the presence of bacteria in urine, is associated with asymptomatic, as well as symptomatic, urinary tract infection (UTI). Bacteriuria underpins some of the dynamics of microbial colonization of the urinary tract, and probably impacts the progression and persistence of infection in some individuals. Recent molecular discoveries in vitro have elucidated how some key bacterial traits can enable organisms to survive and grow in human urine as a means of microbial fitness adaptation for UTI. Several microbial characteristics that confer bacteruric potential have been identified including de novo synthesis of guanine, relative resistance to D-serine, and catabolism of malic acid. Microbial characteristics such as these are increasingly being defined through the use of synthetic human urine (SHU) in vitro as a model to mimic the in vivo environment that bacteria encounter in the bladder. There is considerable variation in the SHU model systems that have been used to study bacteriuria to date, and this influences the utility of these models. In this review, we discuss recent advances in our understanding of bacteruric potential with a focus on the specific mechanisms underlying traits that promote the growth of bacteria in urine. We also review the application of SHU in research studies modeling UTI and discuss the chemical makeup, and benefits and limitations that are encountered in utilizing SHU to study bacterial growth in urine in vitro.

Which microbial factors really are important in Pseudomonas aeruginosa infections?

Over the last two decades, tens of millions of dollars have been invested in understanding virulence in the human pathogen, Pseudomonas aeruginosa. However, the top ‘hits’ obtained in a recent TnSeq analysis aimed at identifying those genes that are conditionally essential for infection did not include most of the known virulence factors identified in these earlier studies. Instead, it seems that P. aeruginosa faces metabolic challenges in vivo, and unless it can overcome these, it fails to thrive and is cleared from the host. In this review, we look at the kinds of metabolic pathways that the pathogen seems to find essential, and comment on how this knowledge might be therapeutically exploited.

Pseudomonas putida KT2440 Strain Metabolizes Glucose through a Cycle Formed by Enzymes of the Entner-Doudoroff, Embden-Meyerhof-Parnas, and Pentose Phosphate Pathways

The soil bacterium Pseudomonas putida KT2440 lacks a functional Embden-Meyerhof-Parnas (EMP) pathway, and glycolysis is known to proceed almost exclusively through the Entner-Doudoroff (ED) route. To investigate the raison d'être of this metabolic arrangement, the distribution of periplasmic and cytoplasmic carbon fluxes was studied in glucose cultures of this bacterium by using (13)C-labeled substrates, combined with quantitative physiology experiments, metabolite quantification, and in vitro enzymatic assays under both saturating and non-saturating, quasi in vivo conditions. Metabolic flux analysis demonstrated that 90% of the consumed sugar was converted into gluconate, entering central carbon metabolism as 6-phosphogluconate and further channeled into the ED pathway. Remarkably, about 10% of the triose phosphates were found to be recycled back to form hexose phosphates. This set of reactions merges activities belonging to the ED, the EMP (operating in a gluconeogenic fashion), and the pentose phosphate pathways to form an unforeseen metabolic architecture (EDEMP cycle). Determination of the NADPH balance revealed that the default metabolic state of P. putida KT2440 is characterized by a slight catabolic overproduction of reducing power. Cells growing on glucose thus run a biochemical cycle that favors NADPH formation. Because NADPH is required not only for anabolic functions but also for counteracting different types of environmental stress, such a cyclic operation may contribute to the physiological heftiness of this bacterium in its natural habitats.

Elementary Flux Mode Analysis Revealed Cyclization Pathway as a Powerful Way for NADPH Regeneration of Central Carbon Metabolism

NADPH regeneration capacity is attracting growing research attention due to its important role in resisting oxidative stress. Besides, NADPH availability has been regarded as a limiting factor in production of industrially valuable compounds. The central carbon metabolism carries the carbon skeleton flux supporting the operation of NADPH-regenerating enzyme and offers flexibility in coping with NADPH demand for varied intracellular environment. To acquire an insightful understanding of its NADPH regeneration capacity, the elementary mode method was employed to compute all elementary flux modes (EFMs) of a network representative of central carbon metabolism. Based on the metabolic flux distributions of these modes, a cluster analysis of EFMs with high NADPH regeneration rate was conducted using the self-organizing map clustering algorithm. The clustering results were used to study the relationship between the flux of total NADPH regeneration and the flux in each NADPH producing enzyme. The results identified several reaction combinations supporting high NADPH regeneration, which are proven to be feasible in cells via thermodynamic analysis and coincident with a great deal of previous experimental report. Meanwhile, the reaction combinations showed some common characteristics: there were one or two decarboxylation oxidation reactions in the combinations that produced NADPH and the combination constitution included certain gluconeogenesis pathways. These findings suggested cyclization pathways as a powerful way for NADPH regeneration capacity of bacterial central carbon metabolism.

Comparative genomics of Roseobacter clade bacteria isolated from the accessory nidamental gland of Euprymna scolopes

The accessory nidamental gland (ANG) of the female Hawaiian bobtail squid, Euprymna scolopes, houses a consortium of bacteria including members of the Flavobacteriales, Rhizobiales, and Verrucomicrobia but is dominated by members of the Roseobacter clade (Rhodobacterales) within the Alphaproteobacteria. These bacteria are deposited into the jelly coat of the squid’s eggs, however, the function of the ANG and its bacterial symbionts has yet to be elucidated. In order to gain insight into this consortium and its potential role in host reproduction, we cultured 12 Rhodobacterales isolates from ANGs of sexually mature female squid and sequenced their genomes with Illumina sequencing technology. For taxonomic analyses, the ribosomal proteins of 79 genomes representing both roseobacters and non-roseobacters along with a separate MLSA analysis of 33 housekeeping genes from Roseobacter organisms placed all 12 isolates from the ANG within two groups of a single Roseobacter clade. Average nucelotide identity analysis suggests the ANG isolates represent three genera (Leisingera, Ruegeria, and Tateyamaria) comprised of seven putative species groups. All but one of the isolates contains a predicted Type VI secretion system, which has been shown to be important in secreting signaling and/or effector molecules in host–microbe associations and in bacteria–bacteria interactions. All sequenced genomes also show potential for secondary metabolite production, and are predicted to be involved with the production of acyl homoserine lactones (AHLs) and/or siderophores. An AHL bioassay confirmed AHL production in three tested isolates and from whole ANG homogenates. The dominant symbiont, Leisingera sp. ANG1, showed greater viability in iron-limiting conditions compared to other roseobacters, possibly due to higher levels of siderophore production. Future comparisons will try to elucidate novel metabolic pathways of the ANG symbionts to understand their putative role in host development.

Large-Scale 13C flux profiling reveals conservation of the Entner-Doudoroff pathway as a glycolytic strategy among marine bacteria that use glucose

Marine bacteria form one of the largest living surfaces on Earth, and their metabolic activity is of fundamental importance for global nutrient cycling. Here, we explored the largely unknown intracellular pathways in 25 microbes representing different classes of marine bacteria that use glucose: Alphaproteobacteria, Gammaproteobacteria, and Flavobacteriia of the Bacteriodetes phylum. We used (13)C isotope experiments to infer metabolic fluxes through their carbon core pathways. Notably, 90% of all strains studied use the Entner-Doudoroff (ED) pathway for glucose catabolism, whereas only 10% rely on the Embden-Meyerhof-Parnas (EMP) pathway. This result differed dramatically from the terrestrial model strains studied, which preferentially used the EMP pathway yielding high levels of ATP. Strains using the ED pathway exhibited a more robust resistance against the oxidative stress typically found in this environment. An important feature contributing to the preferential use of the ED pathway in the oceans could therefore be enhanced supply of NADPH through this pathway. The marine bacteria studied did not specifically rely on a distinct anaplerotic route, but the carboxylation of phosphoenolpyruvate (PEP) or pyruvate for fueling of the tricarboxylic acid (TCA) cycle was evenly distributed. The marine isolates studied belong to clades that dominate the uptake of glucose, a major carbon source for bacteria in seawater. Therefore, the ED pathway may play a significant role in the cycling of mono- and polysaccharides by bacterial communities in marine ecosystems.

Multi -omics and metabolic modelling pipelines: challenges and tools for systems microbiology

Integrated -omics approaches are quickly spreading across microbiology research labs, leading to (i) the possibility of detecting previously hidden features of microbial cells like multi-scale spatial organization and (ii) tracing molecular components across multiple cellular functional states. This promises to reduce the knowledge gap between genotype and phenotype and poses new challenges for computational microbiologists. We underline how the capability to unravel the complexity of microbial life will strongly depend on the integration of the huge and diverse amount of information that can be derived today from -omics experiments. In this work, we present opportunities and challenges of multi -omics data integration in current systems biology pipelines. We here discuss which layers of biological information are important for biotechnological and clinical purposes, with a special focus on bacterial metabolism and modelling procedures. A general review of the most recent computational tools for performing large-scale datasets integration is also presented, together with a possible framework to guide the design of systems biology experiments by microbiologists.