Why metabolic enzymes are essential or nonessential for growth of Escherichia coli K12 on glucose

Machine learning identifies key metabolic reactions in bacterial growth on different carbon sources

Carbon source-dependent control of bacterial growth is fundamental to bacterial physiology and survival. However, pinpointing the metabolic steps important for cell growth is challenging due to the complexity of cellular networks. Here, the elastic net model and multilayer perception model that integrated genome-wide gene-deletion data and simulated flux distributions were constructed to identify metabolic reactions beneficial or detrimental to Escherichia coli grown on 30 different carbon sources. Both models outperformed traditional in silico methods by identifying not just essential reactions but also nonessential ones that promote growth. They successfully predicted metabolic reactions beneficial to cell growth, with high convergence between the models. The models revealed that biosynthetic pathways generally promote growth across various carbon sources, whereas the impact of energy-generating pathways varies with the carbon source. Intriguing predictions were experimentally validated for findings beyond experimental training data and the impact of various carbon sources on the glyoxylate shunt, pyruvate dehydrogenase reaction, and redundant purine biosynthesis reactions. These highlight the practical significance and predictive power of the models for understanding and engineering microbial metabolism.

Acinetobacter baumannii coordinates central metabolism, plasmid dissemination, and virulence by sensing nutrient availability

IMPORTANCE

Plasmid conjugation is known to be an energy-expensive process, but our understanding of the molecular linkage between conjugation and metabolism is limited. Our finding reveals that Acinetobacter baumannii utilizes a two-component system to co-regulate metabolism, plasmid transfer, and virulence by sensing reaction intermediates of key metabolic pathways, which suggests that nutrient availability dictates not only bacterial proliferation but also horizontal gene transfer. The identification of Dot/Icm-like proteins as components of a conjugation system involved in the dissemination of antibiotic-resistance genes by A. baumannii has provided important targets for the development of agents capable of inhibiting virulence and the spread of anti-microbial-resistance genes in bacterial communities.

Exchange of Vitamin B1 and Its Biosynthesis Intermediates Shapes the Composition of Synthetic Microbial Cocultures and Reveals Complexities of Nutrient Sharing

Microbial communities occupy diverse niches in nature, and community members routinely exchange a variety of nutrients among themselves. While large-scale metagenomic and metabolomic studies shed some light on these exchanges, the contribution of individual species and the molecular details of specific interactions are difficult to track. In this study, we follow the exchange of vitamin B₁ (thiamin) and its intermediates between microbes within synthetic cocultures of Escherichia coli and Vibrio anguillarum. Thiamin contains two moieties, 4-amino-5-hydroxymethyl-2-methylpyrimidine (HMP) and 4-methyl-5-(2-hydroxyethyl)thiazole (THZ), which are synthesized by distinct pathways using enzymes ThiC and ThiG, respectively, and then coupled by ThiE to form thiamin. Even though E. coli ΔthiC, ΔthiE, and ΔthiG mutants are thiamin auxotrophs, we observed that cocultures of ΔthiC-ΔthiE and ΔthiC-ΔthiG mutants are able to grow in a thiamin-deficient medium, whereas the ΔthiE-ΔthiG coculture does not. Further, the exchange of thiamin and its intermediates in V. anguillarum cocultures and in mixed cocultures of V. anguillarum and E. coli revealed that there exist specific patterns for thiamin metabolism and exchange among these microbes. Our findings show that HMP is shared more frequently than THZ, concurrent with previous observations that free HMP and HMP auxotrophy is commonly found in various environments. Furthermore, we observe that the availability of exogenous thiamin in the media affects whether these strains interact with each other or grow independently. These findings collectively underscore the importance of the exchange of essential metabolites as a defining factor in building and modulating synthetic or natural microbial communities. IMPORTANCE Vitamin B₁ (thiamin) is an essential nutrient for cellular metabolism. Microorganisms that are unable to synthesize thiamin either fully or in part exogenously obtain it from their environment or via exchanges with other microbial members in their community. In this study, we created synthetic microbial cocultures that rely on sharing thiamin and its biosynthesis intermediates and observed that some of them are preferentially exchanged. We also observed that the coculture composition is dictated by the production and/or availability of thiamin and its intermediates. Our studies with synthetic cocultures provide the molecular basis for understanding thiamin sharing among microorganisms and lay out broad guidelines for setting up synthetic microbial cocultures by using the exchange of an essential metabolite as their foundation.

Biochemical and functional characterization of Mycobacterium tuberculosis ketol-acid reductoisomerase

Branched-chain amino acids (BCAAs) are essential amino acids, but their biosynthetic pathway is absent in mammals. Ketol-acid reductoisomerase (IlvC) is a BCAA biosynthetic enzyme that is coded by Rv3001c in Mycobacterium tuberculosis H37Rv (Mtb-Rv) and MRA_3031 in M. tuberculosis H37Ra (Mtb-Ra). IlvCs are essential in Mtb-Rv as well as in Escherichia coli. Compared to wild-type and IlvC-complemented Mtb-Ra strains, IlvC knockdown strain showed reduced survival at low pH and under low pH+starvation stress conditions. Further, increased expression of IlvC was observed under low pH and starvation stress conditions. Confirmation of a role for IlvC in pH and starvation stress was achieved by developing E. coli BL21(DE3) IlvC knockout, which was defective for growth in M9 minimal medium, but growth could be rescued by isoleucine and valine supplementation. Growth was also restored by complementing with over-expressing constructs of Mtb-Ra and E. coli IlvCs. The E. coli knockout also had a survival deficit at pH=5.5 and 4.5 and was more susceptible to killing at pH=3.0. The biochemical characterization of Mtb-Ra and E. coli IlvCs confirmed that both have NADPH-dependent activity. In conclusion, this study demonstrates the functional complementation of E. coli IlvC by Mtb-Ra IlvC and also suggests that IlvC has a role in tolerance to low pH and starvation stress.

Preferential use of carbon central metabolism and anaerobic respiratory chains in porcine extraintestinal pathogenic Escherichia coli during bloodstream infection

Porcine extraintestinal pathogenic Escherichia coli (ExPEC) is occurring with increasing frequency in China, and leads to significant economic and welfare costs in the swine industry. The underlying mechanisms of porcine ExPEC in blood colonization during systematic infection is poorly understood. Here we measured the gene expression of porcine ExPEC in infected animal bloodstream in vivo and fresh swine blood in vitro. Using comparisons with P values of ≤ 0.01, we identified 354 and 313 genes as being significantly up- or down-regulated at least 2-fold change during bloodstream infection, respectively. Excepting for an array of iron acquisition systems, numerous genes involved in carbon central metabolism and anaerobic respiratory chains were upregulated here. These genes were categorized into several clusters including the TCA-cycle (frdABCD, citCEFXG), d-ribose transporter (rbsDACB), nickel transporter (nikABCDER), NiFe hydrogenase (hybOABCDEF, hycBCDEFG), Hyp-complex (hypABCDE), DMSO reductase (dmsABC and ynfEFGHI), format dehydrogenase (fdnGHI) and NADH dehydrogenase I (nuoA-N). The mutant with simultaneous inactivation of ribose and citrate imports showed significant reduced fitness in host blood, suggesting these two carbohydrates are utilized by central metabolism network as important carbon-source during bloodstream infection. Similar deficiency was also observed in the mutant double deleted NiFe hydrogenase 2 and 3 anaerobic respiratory chains. Further study found that FNR (a global regulator facilitating bacterial adaptation to anaerobic conditions) is an important regulator in response to bloodstream to activate center metabolism and anaerobic respiratory chains, thus contribute to the full-virulence of porcine ExPEC. These findings provide compelling evidence to support the notion that carbon central metabolism network and anaerobic respiratory chains play key roles for porcine ExPEC fitness within host bloodstream.

In silico co-factor balance estimation using constraint-based modelling informs metabolic engineering in Escherichia coli

In the growing field of metabolic engineering, where cells are treated as ‘factories’ that synthesize industrial compounds, it is essential to consider the ability of the cells’ native metabolism to accommodate the demands of synthetic pathways, as these pathways will alter the homeostasis of cellular energy and electron metabolism. From the breakdown of substrate, microorganisms activate and reduce key co-factors such as ATP and NAD(P)H, which subsequently need to be hydrolysed and oxidized, respectively, in order to restore cellular balance. A balanced supply and consumption of such co-factors, here termed co-factor balance, will influence biotechnological performance. To aid the strain selection and design process, we used stoichiometric modelling (FBA, pFBA, FVA and MOMA) and the Escherichia coli (E.coli) core stoichiometric model to investigate the network-wide effect of butanol and butanol precursor production pathways differing in energy and electron demand on product yield. An FBA-based co-factor balance assessment (CBA) algorithm was developed to track and categorise how ATP and NAD(P)H pools are affected in the presence of a new pathway. CBA was compared to the balance calculations proposed by Dugar et al. (Nature Biotechnol. 29 (12), 1074–1078). Predicted solutions were compromised by excessively underdetermined systems, displaying greater flexibility in the range of reaction fluxes than experimentally measured by ¹³C-metabolic flux analysis (MFA) and the appearance of unrealistic futile co-factor cycles. With the assumption that futile cycles are tightly regulated in reality, the FBA models were manually constrained in a step-wise manner. Solutions with minimal futile cycling diverted surplus energy and electrons towards biomass formation. As an alternative, the use of loopless FBA or constraining the models with measured flux ranges were tried but did not prevent futile co-factor cycles. The results highlight the need to account for co-factor imbalance and confirm that better-balanced pathways with minimal diversion of surplus towards biomass formation present the highest theoretical yield. The analysis also suggests that ATP and NAD(P)H balancing cannot be assessed in isolation from each other, or even from the balance of additional co-factors such as AMP and ADP. We conclude that, through revealing the source of co-factor imbalance CBA can facilitate pathway and host selection when designing new biocatalysts for implementation by metabolic engineering.

The chemicals industry is a major contributor to greenhouse gas emissions and desperately requires more sustainable alternatives. Genetically engineered microorganisms can be used as ‘bio-factories’ to manufacture chemicals, replacing those currently sourced from fossil fuels or unsustainable tropical plant agriculture. However, due to the complexity of biology, the features that render one bio-factory design more efficient than others are difficult to identify. Computational modelling of such designs can enable the selection of optimally performing designs, but it remains challenging as biology is complex and not fully understood. Microorganisms require energy for their own growth and maintenance, but also to convert molecules into desired target products. The supply and consumption of such energy is through co-factors, and the balance of such co-factors influences the performance of the engineered bio-factories. This study developed a computer-aided approach for quantification of the co-factor balance of bio-factories. Using the chemical n-butanol as a case study, our study explores the impact of variant bio-factory designs with differing co-factor balance on the potential efficiency of biomanufacturing. We provide insights into the relative balance of different designs and provide a computational framework to select the best-performing designs.

Activity evaluation of glycolytic promoters from Escherichia coli and application for mevalonate biosynthesis

Promoters are the most important tools to control and regulate gene expression in synthetic biology and metabolic engineering. The expression of target genes in Escherichia coli is usually controlled by inducible promoters which may cause excessive metabolic load on the host and may be uneconomical due to inducer cost. Therefore, it is important to identify more constitutive promoters that are capable of avoiding these limitations. In this study, using the monomeric red fluorescent protein (mRFP) as the reporter gene, ten promoters from the glycolytic pathway of E. coli were cloned and characterized. We found that glycolytic promoters exhibited the advantages of constitutive promoters and higher strength compared with the commonly used inducible promoter Plac. We further introduced glycolytic promoters into the mevalonate biosynthesis system. The maximum mevalonate titer produced by engineered E. coli under the control of glycolytic promoters was obviously higher than that under the control of Plac, indicating the superiority of glycolytic promoter for the metabolic engineering of E. coli. This set of glycolytic promoters significantly expands the range of engineering tools available for E. coli and can be applied in future metabolic engineering studies.

sRNA-Based Screening Chromosomal Gene Targets and Modular Designing Escherichia coli for High-Titer Production of Aglycosylated Immunoglobulin G

The production of the aglycosylated immunoglobulin G (IgG) in Escherichia coli has received wide interest for its analytical and therapeutic applications. To enhance the production titer of IgG, we first used synthetic sRNAs to perform a systematical analysis of the gene expression in the translational level in the glycolytic pathway (module 1) and the tricarboxylic acid (TCA) cycle (module 2) to reveal the critical genes for the efficient IgG production. Second, to provide sufficient amino acid precursors for the protein biosynthesis, amino acid biosynthesis pathways (module 3) were enhanced to facilitate the IgG production. Upon integrated engineering of these genes in the three modules (module 1, aceF; module 2, gltA and acnA; module 3, serB) and optimization of fermentation conditions, the recombinant E. coli enabled a titer of the full-assembled IgG of 4.5 ± 0.6 mg/L in flask cultures and 184 ± 9.2 mg/L in the 5 L high cell density fed-batch fermenter, which is, as far as we know, the highest reported titer of IgG production in recombinant E. coli.

Cladogenesis and Genomic Streamlining in Extracellular Endosymbionts of Tropical Stink Bugs

Phytophagous stink bugs are globally distributed and many harbor vertically inherited bacterial symbionts that are extracellular, yet little is known about how the symbiont’s genomes have evolved under this transmission strategy. Genome reduction is common in insect intracellular symbionts but limited genome sampling of the extracellular symbionts of distantly related stink bugs has precluded inferring patterns of extracellular symbiont genome evolution. To address this knowledge gap, we completely sequenced the genomes of the uncultivable bacterial symbionts of four neotropical stink bugs of the Edessa genus. Phylogenetic and comparative analyses indicated that the symbionts form a clade within the Pantoea genus and their genomes are highly reduced (∼0.8 Mb). Furthermore, genome synteny analysis and a jackknife approach for phylogenetic reconstruction, which corrected for long branch attraction artifacts, indicated that the Edessa symbionts were the result of a single symbiotic event that was distinct from the symbiosis event giving rise to Candidatus “Pantoea carbekii,” the extracellular symbiont of the invasive pentatomid stink bug, Halyomorpha halys. Metabolic functions inferred from the Edessa symbiont genomes suggests a shift in genomic composition characteristic of its lifestyle in that they retained many host-supportive functions while undergoing dramatic gene loss and establishing a stable relationship with their host insects. Given the undersampled nature of extracellular insect symbionts, this study is the first comparative analysis of these symbiont genomes from four distinct Edessa stink bug species. Finally, we propose the candidate name “Candidatus Pantoea edessiphila” for the species of these symbionts with strain designations according to their host species.

Overexpression of thermostable meso-diaminopimelate dehydrogenase to redirect diaminopimelate pathway for increasing L-lysine production in Escherichia coli

Dehydrogenase pathway, one of diaminopimelate pathway, is important to the biosynthesis of L-lysine and peptidoglycan via one single reaction catalyzed by meso-diaminopimelate dehydrogenase (DapDH). In this study, the thermostable DapDH was introduced into diaminopimelate pathway that increased the final titer (from 71.8 to 119.5 g/L), carbon yield (from 35.3% to 49.1%) and productivity (from 1.80 to 2.99 g/(L∙h)) of L-lysine by LATR12-2∆rpiB::ddh_St in fed-batch fermentation. To do this, the kinetic properties and the effects of different DapDHs on L-lysine production were investigated, and the results indicated that overexpression of StDapDH in LATR12-2 was beneficial to construct an L-lysine producer with good productive performance because it exhibited the best of kinetic characteristics and optimal temperature as well as thermostability in reductive amination. Furthermore, ammonium availability was optimized, and found that 20 g/L of (NH₄)₂SO₄ was the optimal ammonium concentration for improving the efficiency of L-lysine production by LATR12-2∆rpiB::ddh_St. Metabolomics analysis showed that introducing the StDapDH significantly enhanced carbon flux into pentose phosphate pathway and L-lysine biosynthetic pathway, thus increasing the levels of NADPH and precursors for L-lysine biosynthesis. This is the first report of a rational modification of diaminopimelate pathway that improves the efficiency of L-lysine production through overexpression of thermostable DapDH in E. coli.

The antimicrobial action of polyaniline involves production of oxidative stress while functionalisation of polyaniline introduces additional mechanisms

Polyaniline (PANI) and functionalised polyanilines (fPANI) are novel antimicrobial agents whose mechanism of action was investigated. Escherichia coli single gene deletion mutants revealed that the antimicrobial mechanism of PANI likely involves production of hydrogen peroxide while homopolymer poly(3-aminobenzoic acid), P3ABA, used as an example of a fPANI, disrupts metabolic and respiratory machinery, by targeting ATP synthase and causes acid stress. PANI was more active against E. coli in aerobic, compared to anaerobic, conditions, while this was apparent for P3ABA only in rich media. Greater activity in aerobic conditions suggests involvement of reactive oxygen species. P3ABA treatment causes an increase in intracellular free iron, which is linked to perturbation of metabolic enzymes and could promote reactive oxygen species production. Addition of exogenous catalase protected E. coli from PANI antimicrobial action; however, this was not apparent for P3ABA treated cells. The results presented suggest that PANI induces production of hydrogen peroxide, which can promote formation of hydroxyl radicals causing biomolecule damage and potentially cell death. P3ABA is thought to act as an uncoupler by targeting ATP synthase resulting in a futile cycle, which precipitates dysregulation of iron homeostasis, oxidative stress, acid stress, and potentially the fatal loss of proton motive force.

Growth medium-dependent antimicrobial activity of early stage MEP pathway inhibitors

The in vivo microenvironment of bacterial pathogens is often characterized by nutrient limitation. Consequently, conventional rich in vitro culture conditions used widely to evaluate antibacterial agents are often poorly predictive of in vivo activity, especially for agents targeting metabolic pathways. In one such pathway, the methylerythritol phosphate (MEP) pathway, which is essential for production of isoprenoids in bacterial pathogens, relatively little is known about the influence of growth environment on antibacterial properties of inhibitors targeting enzymes in this pathway. The early steps of the MEP pathway are catalyzed by 1-deoxy-d-xylulose 5-phosphate (DXP) synthase and reductoisomerase (IspC). The in vitro antibacterial efficacy of the DXP synthase inhibitor butylacetylphosphonate (BAP) was recently reported to be strongly dependent upon growth medium, with high potency observed under nutrient limitation and exceedingly weak activity in nutrient-rich conditions. In contrast, the well-known IspC inhibitor fosmidomycin has potent antibacterial activity in nutrient-rich conditions, but to date, its efficacy had not been explored under more relevant nutrient-limited conditions. The goal of this work was to thoroughly characterize the effects of BAP and fosmidomycin on bacterial cells under varied growth conditions. In this work, we show that activities of both inhibitors, alone and in combination, are strongly dependent upon growth medium, with differences in cellular uptake contributing to variance in potency of both agents. Fosmidomycin is dissimilar to BAP in that it displays relatively weaker activity in nutrient-limited compared to nutrient-rich conditions. Interestingly, while it has been generally accepted that fosmidomycin activity depends upon expression of the GlpT transporter, our results indicate for the first time that fosmidomycin can enter cells by an alternative mechanism under nutrient limitation. Finally, we show that the potency and relationship of the BAP-fosmidomycin combination also depends upon the growth medium, revealing a striking loss of BAP-fosmidomycin synergy under nutrient limitation. This change in BAP-fosmidomycin relationship suggests a shift in the metabolic and/or regulatory networks surrounding DXP accompanying the change in growth medium, the understanding of which could significantly impact targeting strategies against this pathway. More generally, our findings emphasize the importance of considering physiologically relevant growth conditions for predicting the antibacterial potential MEP pathway inhibitors and for studies of their intracellular targets.

Re-Factoring Glycolytic Genes for Targeted Engineering of Catabolism in Gram-Negative Bacteria

The Embden-Meyerhof-Parnas (EMP) pathway is widely accepted to be the biochemical standard of glucose catabolism. The well-characterized glycolytic route of Escherichia coli, based on the EMP catabolism, is an example of an intricate pathway in terms of genomic organization of the genes involved and patterns of gene expression and regulation. This intrinsic genetic and metabolic complexity renders it difficult to engineer glycolytic activities and transfer them onto other microbial cell factories, thus limiting the biotechnological potential of bacterial hosts that lack the route. Taking into account the potential applications of such a portable tool for targeted pathway engineering, in the present protocol we describe how 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 (i.e., GlucoBrick) that facilitates their grouping in the form of functional modules that can be combined at the user's will. This novel genetic tool allows for the à la carte implementation or boosting of EMP pathway activities into different Gram-negative bacteria. The potential of the GlucoBrick platform is further illustrated by engineering novel glycolytic activities in the most representative members of the Pseudomonas genus (Pseudomonas putida and Pseudomonas aeruginosa).

Analysis of metabolic networks of Streptomyces leeuwenhoekii C34 by means of a genome scale model: Prediction of modifications that enhance the production of specialized metabolites

The first genome scale model (GSM) for Streptomyces leeuwenhoekii C34 was developed to study the biosynthesis pathways of specialized metabolites and to find metabolic engineering targets for enhancing their production. The model, iVR1007, consists of 1,722 reactions, 1,463 metabolites, and 1,007 genes, it includes the biosynthesis pathways of chaxamycins, chaxalactins, desferrioxamines, ectoine, and other specialized metabolites. iVR1007 was validated using experimental information of growth on 166 different sources of carbon, nitrogen and phosphorous, showing an 83.7% accuracy. The model was used to predict metabolic engineering targets for enhancing the biosynthesis of chaxamycins and chaxalactins. Gene knockouts, such as sle03600 (L-homoserine O-acetyltransferase), and sle39090 (trehalose-phosphate synthase), that enhance the production of the specialized metabolites by increasing the pool of precursors were identified. Using the algorithm of flux scanning based on enforced objective flux (FSEOF) implemented in python, 35 and 25 over-expression targets for increasing the production of chaxamycin A and chaxalactin A, respectively, that were not directly associated with their biosynthesis routes were identified. Nineteen over-expression targets that were common to the two specialized metabolites studied, like the over-expression of the acetyl carboxylase complex (sle47660 (accA) and any of the following genes: sle44630 (accA_1) or sle39830 (accA_2) or sle27560 (bccA) or sle59710) were identified. The predicted knockouts and over-expression targets will be used to perform metabolic engineering of S. leeuwenhoekii C34 and obtain overproducer strains.

Gracob: a novel graph-based constant-column biclustering method for mining growth phenotype data

Motivation

Growth phenotype profiling of genome-wide gene-deletion strains over stress conditions can offer a clear picture that the essentiality of genes depends on environmental conditions. Systematically identifying groups of genes from such high-throughput data that share similar patterns of conditional essentiality and dispensability under various environmental conditions can elucidate how genetic interactions of the growth phenotype are regulated in response to the environment.

Results

We first demonstrate that detecting such ‘co-fit’ gene groups can be cast as a less well-studied problem in biclustering, i.e. constant-column biclustering. Despite significant advances in biclustering techniques, very few were designed for mining in growth phenotype data. Here, we propose Gracob, a novel, efficient graph-based method that casts and solves the constant-column biclustering problem as a maximal clique finding problem in a multipartite graph. We compared Gracob with a large collection of widely used biclustering methods that cover different types of algorithms designed to detect different types of biclusters. Gracob showed superior performance on finding co-fit genes over all the existing methods on both a variety of synthetic data sets with a wide range of settings, and three real growth phenotype datasets for E. coli, proteobacteria and yeast.

Availability and Implementation

Our program is freely available for download at http://sfb.kaust.edu.sa/Pages/Software.aspx.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Challenges and Hallmarks of Establishing Alkylacetylphosphonates as Probes of Bacterial 1-Deoxy-d-xylulose 5-Phosphate Synthase

1-Deoxy-d-xylulose 5-phosphate (DXP) synthase catalyzes the thiamin diphosphate (ThDP)-dependent formation of DXP from pyruvate and d-glyceraldehyde 3-phosphate. DXP is at a metabolic branch point in bacteria, feeding into the methylerythritol phosphate pathway to indispensable isoprenoids and acting as a precursor for biosynthesis of essential cofactors in central metabolism, pyridoxal phosphate and ThDP, the latter of which is also required for DXP synthase catalysis. DXP synthase follows a unique random sequential mechanism and possesses an unusually large active site. These features have guided the design of sterically demanding alkylacetylphosphonates (alkylAPs) toward the development of selective DXP synthase inhibitors. alkylAPs studied here display selective, low μM inhibitory activity against DXP synthase. They are weak inhibitors of bacterial growth in standard nutrient rich conditions. However, bacteria are significantly sensitized to most alkylAPs in defined minimal growth medium, with minimal inhibitory concentrations (MICs) ranging from low μM to low mM and influenced by alkyl-chain length. The longest analog (C₈) displays the weakest antimicrobial activity and is a substrate for efflux via AcrAB-TolC. The dependence of inhibitor potency on growth environment emphasizes the need for antimicrobial screening conditions that are relevant to the in vivo microbial microenvironment during infection. DXP synthase expression and thiamin supplementation studies offer support for DXP synthase as an intracellular target for some alkylAPs and reveal both the challenges and intriguing aspects of these approaches to study target engagement.

The In Vitro Redundant Enzymes PurN and PurT Are Both Essential for Systemic Infection of Mice in Salmonella enterica Serovar Typhimurium

Metabolic enzymes show a high degree of redundancy, and for that reason they are generally ignored in searches for novel targets for anti-infective substances. The enzymes PurN and PurT are redundant in vitro in Salmonella enterica serovar Typhimurium, in which they perform the third step of purine synthesis. Surprisingly, the results of the current study demonstrated that single-gene deletions of each of the genes encoding these enzymes caused attenuation (competitive infection indexes [CI] of <0.03) in mouse infections. While the ΔpurT mutant multiplied as fast as the wild-type strain in cultured J774A.1 macrophages, net multiplication of the ΔpurN mutant was reduced approximately 50% in 20 h. The attenuation of the ΔpurT mutant was abolished by simultaneous removal of the enzyme PurU, responsible for the formation of formate, indicating that the attenuation was related to formate accumulation or wasteful consumption of formyl tetrahydrofolate by PurU. In the process of further characterization, we disclosed that the glycine cleavage system (GCV) was the most important for formation of C1 units in vivo (CI = 0.03 ± 0.03). In contrast, GlyA was the only important enzyme for the formation of C1 units in vitro The results with the ΔgcvT mutant further revealed that formation of serine by SerA and further conversion of serine into C1 units and glycine by GlyA were not sufficient to ensure C1 formation in S Typhimurium in vivo The results of the present study call for reinvestigations of the concept of metabolic redundancy in S Typhimurium in vivo.

Energy Utilization for Survival and Fertilization-Parsimonious Quiescent Sperm Turn Extravagant on Motility Activation in Rat

Quiescent sperm survive in cauda epididymis for long periods of time under extreme crowding conditions and with a very limited energy substrate, while after ejaculation, motile sperm live for a much shorter period with an unlimited energy resource and without crowding. Thus, the energy metabolism in relation to the energy requirement of the two may be quite different. A simple physiological technique was evolved to collect viable quiescent sperm from rat cauda epididymis to compare its energy metabolism with motile sperm. Quiescent sperm exhibited 40%-60% higher activities of mitochondrial electron transport chain complexes I-IV and ATP synthase in comparison to motile sperm and accumulated Ca(2+) in the midpiece mitochondria to enhance oxidative phosphorylation (OxPhos). In contrast, motile sperm displayed up to 75% higher activities of key glycolytic enzymes and secreted more than two times the lactate than quiescent sperm. Quiescent sperm phosphorylated AMPK and MAPK-p38, while motile sperm phosphorylated AKT and MAPK/ERK. Glycolytic inhibitor iodoacetamide prevented motility activation of quiescent rat sperm and inhibited conception in rabbits more effectively than OxPhos uncoupler 2,4-dinitrophenol. Apparently, quiescent sperm employ the most energy efficient OxPhos to survive for extended periods of time under extreme conditions of nutrition and crowding. However, on motility initiation, sperm switch predominantly to glycolysis to cater to their high- and quick-energy requirement of much shorter periods. This study also presents a proof of concept for targeting sperm energy metabolism for contraception.

Escherichia coli Pyruvate Dehydrogenase Complex Is an Important Component of CXCL10-Mediated Antimicrobial Activity

Chemokines are best recognized for their role within the innate immune system as chemotactic cytokines, signaling and recruiting host immune cells to sites of infection. Certain chemokines, such as CXCL10, have been found to play an additional role in innate immunity, mediating CXCR3-independent killing of a diverse array of pathogenic microorganisms. While this is still not clearly understood, elucidating the mechanisms underlying chemokine-mediated antimicrobial activity may facilitate the development of novel therapeutic strategies effective against antibiotic-resistant Gram-negative pathogens. Here, we show that CXCL10 exerts antibacterial effects on clinical and laboratory strains of Escherichia coli and report that disruption of pyruvate dehydrogenase complex (PDHc), which converts pyruvate to acetyl coenzyme A, enables E. coli to resist these antimicrobial effects. Through generation and screening of a transposon mutant library, we identified two mutants with increased resistance to CXCL10, both with unique disruptions of the gene encoding the E1 subunit of PDHc, aceE. Resistance to CXCL10 also occurred following deletion of either aceF or lpdA, genes that encode the remaining two subunits of PDHc. Although PDHc resides within the bacterial cytosol, electron microscopy revealed localization of immunogold-labeled CXCL10 to the bacterial cell surface in both the E. coli parent and aceE deletion mutant strains. Taken together, our findings suggest that while CXCL10 interacts with an as-yet-unidentified component on the cell surface, PDHc is an important mediator of killing by CXCL10. To our knowledge, this is the first description of PDHc as a key bacterial component involved in the antibacterial effect of a chemokine.

Enhancement of 3-hydroxypropionic acid production from glycerol by using a metabolic toggle switch

Background

3-hydroxypropionic acid (3-HP) is an important platform for the production of C3 chemicals, including acrylic acid, methyl acrylate, and acrylamide. Microbial production of 3-HP is mainly due to glycerol metabolism. In this study, in order to improve microbial 3-HP production, we applied a metabolic toggle switch for controlling the glycerol metabolism to redirect the excess metabolic flux of central metabolic pathway toward an exogenous 3-HP producing pathway in Escherichia coli.

Results

The metabolic toggle switch enables conditional repression of the expression of a target gene during the fermentation. We individually performed conditional repression of glpK, tpiA, and gapA, which are involved in glycerol metabolism. The conditional repression of glpK and tpiA was not effective for 3-HP production under our experimental conditions. However, gapA conditional repression contributed to improve 3-HP production (titer, 54.2 ± 1.5 mM; yield, 32.1 ± 1.3 %) compared with that for the wild type strain. Additional deletion of endogenous yqhD, which is responsible for the production of a major byproduct, 1,3-propandiol, further increased 3-HP production (titer, 67.3 ± 2.1 mM; yield, 51.5 ± 3.2 %). The titer and yield were 80 and 94 % higher than those of the wild type strain, respectively. The obtained 3-HP yield from glycerol is comparable with the highest yield ever reported for microbial 3-HP production using glycerol as a sole carbon source. The measurement of intracellular metabolites showed the metabolic toggle switch successfully controlled the metabolic flux.

Conclusion

The conditional repression of gapA by using the metabolic toggle switch combined with deletion of endogeneous yqhD increased 3-HP production approximately twofold from glycerol. This result indicates the metabolic toggle switch can be applied in various bio-production using diverse substrates.

Assessing the Metabolic Diversity of Streptococcus from a Protein Domain Point of View

Understanding the diversity and robustness of the metabolism of bacteria is fundamental for understanding how bacteria evolve and adapt to different environments. In this study, we characterised 121 Streptococcus strains and studied metabolic diversity from a protein domain perspective. Metabolic pathways were described in terms of the promiscuity of domains participating in metabolic pathways that were inferred to be functional. Promiscuity was defined by adapting existing measures based on domain abundance and versatility. The approach proved to be successful in capturing bacterial metabolic flexibility and species diversity, indicating that it can be described in terms of reuse and sharing functional domains in different proteins involved in metabolic activity. Additionally, we showed striking differences among metabolic organisation of the pathogenic serotype 2 Streptococcus suis and other strains.

Habitat visualization and genomic analysis of "Candidatus Pantoea carbekii," the primary symbiont of the brown marmorated stink bug

Phytophagous pentatomid insects can negatively impact agricultural productivity and the brown marmorated stink bug (Halyomorpha halys) is an emerging invasive pest responsible for damage to many fruit crops and ornamental plants in North America. Many phytophagous stink bugs, including H. halys, harbor gammaproteobacterial symbionts that likely contribute to host development, and characterization of symbiont transmission/acquisition and their contribution to host fitness may offer alternative strategies for managing pest species. “Candidatus Pantoea carbekii” is the primary occupant of gastric ceca lumina flanking the distal midgut of H. halys insects and it is acquired each generation when nymphs feed on maternal extrachorion secretions following hatching. Insects prevented from symbiont uptake exhibit developmental delays and aberrant behaviors. To infer contributions of Ca. P. carbekii to H. halys, the complete genome was sequenced and annotated from a North American H. halys population. Overall, the Ca. P. carbekii genome is nearly one-fourth (1.2 Mb) that of free-living congenerics, and retains genes encoding many functions that are potentially host-supportive. Gene content reflects patterns of gene loss/retention typical of intracellular mutualists of plant-feeding insects. Electron and fluorescence in situ microscopic imaging of H. halys egg surfaces revealed that maternal extrachorion secretions were populated with Ca. P. carbekii cells. The reported findings detail a transgenerational mode of symbiont transmission distinct from that observed for intracellular insect mutualists and illustrate the potential additive functions contributed by the bacterial symbiont to this important agricultural pest.

Tailoring specialized metabolite production in streptomyces

Streptomycetes are prolific producers of a plethora of medically useful metabolites. These compounds are made by complex secondary (specialized) metabolic pathways, which utilize primary metabolic intermediates as building blocks. In this review we discuss the evolution of specialized metabolites and how expansion of gene families in primary metabolism has lead to the evolution of diversity in these specialized metabolic pathways and how developing a better understanding of expanded primary metabolic pathways can help enhance synthetic biology approaches to industrial pathway engineering.

Pyranose dehydrogenase ligand promiscuity: a generalized approach to simulate monosaccharide solvation, binding, and product formation

The flavoenzyme pyranose dehydrogenase (PDH) from the litter decomposing fungus Agaricus meleagris oxidizes many different carbohydrates occurring during lignin degradation. This promiscuous substrate specificity makes PDH a promising catalyst for bioelectrochemical applications. A generalized approach to simulate all 32 possible aldohexopyranoses in the course of one or a few molecular dynamics (MD) simulations is reported. Free energy calculations according to the one-step perturbation (OSP) method revealed the solvation free energies (ΔG_solv) of all 32 aldohexopyranoses in water, which have not yet been reported in the literature. The free energy difference between β- and α-anomers (ΔG_β-α) of all d-stereoisomers in water were compared to experimental values with a good agreement. Moreover, the free-energy differences (ΔG) of the 32 stereoisomers bound to PDH in two different poses were calculated from MD simulations. The relative binding free energies (ΔΔG_bind) were calculated and, where available, compared to experimental values, approximated from K_m values. The agreement was very good for one of the poses, in which the sugars are positioned in the active site for oxidation at C1 or C2. Distance analysis between hydrogens of the monosaccharide and the reactive N5-atom of the flavin adenine dinucleotide (FAD) revealed that oxidation is possible at HC1 or HC2 for pose A, and at HC3 or HC4 for pose B. Experimentally detected oxidation products could be rationalized for the majority of monosaccharides by combining ΔΔG_bind and a reweighted distance analysis. Furthermore, several oxidation products were predicted for sugars that have not yet been tested experimentally, directing further analyses. This study rationalizes the relationship between binding free energies and substrate promiscuity in PDH, providing novel insights for its applicability in bioelectrochemistry. The results suggest that a similar approach could be applied to study promiscuity of other enzymes.

Generally, enzymes are perceived as being specific for both their substrates and the reaction they catalyze. This standard paradigm started to shift and currently enzyme promiscuity towards various substrates is perceived rather as the rule than the exception. Enzyme promiscuity seems to be vital for proteins to acquire new functions, and therefore for evolution itself. The driving forces for promiscuity are manifold and consequently challenging to study. Binding free energies, which can be calculated from computer simulations, represent a convenient measure for them. Here, we investigate the binding free energies between the enzyme pyranose dehydrogenase (PDH) and a sugar-substrate by computational means. PDH has an extraordinarily promiscuous substrate-specificity, making it interesting for e.g. bioelectrochemical applications. By introducing modifications to the sugar-structure used for the molecular dynamics simulations, we could simultaneously study all 32 possible aldohexopyranoses from a single simulation. This saves costly computational resources and time for setting up and analyzing the simulations. We could nicely reproduce experimental results and predict so far undetected sugar-oxidation products, directing further experiments. This study gives novel insights into PDH's substrate promiscuity and the enzyme's applicability. A similar approach could be applied to study the promiscuity of other enzymes.

Molecular annotation of ketol-acid reductoisomerases from Streptomyces reveals a novel amino acid biosynthesis interlock mediated by enzyme promiscuity

The 6-phosphogluconate dehydrogenase superfamily oxidize and reduce a wide range of substrates, making their functional annotation challenging. Ketol-acid reductoisomerase (KARI), encoded by the ilvC gene in branched-chain amino acids biosynthesis, is a promiscuous reductase enzyme within this superfamily. Here, we obtain steady-state enzyme kinetic parameters for 10 IlvC homologues from the genera Streptomyces and Corynebacterium, upon eight selected chemically diverse substrates, including some not normally recognized by enzymes of this superfamily. This biochemical data suggested a Streptomyces biosynthetic interlock between proline and the branched-chain amino acids, mediated by enzyme substrate promiscuity, which was confirmed via mutagenesis and complementation analyses of the proC, ilvC1 and ilvC2 genes in Streptomyces coelicolor. Moreover, both ilvC orthologues and paralogues were analysed, such that the relationship between gene duplication and functional diversification could be explored. The KARI paralogues present in S. coelicolor and Streptomyces lividans, despite their conserved high sequence identity (97%), were shown to be more promiscuous, suggesting a recent functional diversification. In contrast, the KARI paralogue from Streptomyces viridifaciens showed selectivity towards the synthesis of valine precursors, explaining its recruitment within the biosynthetic gene cluster of valanimycin. These results allowed us to assess substrate promiscuity indices as a tool to annotate new molecular functions with metabolic implications.

A genome-scale metabolic flux model of Escherichia coli K-12 derived from the EcoCyc database

BACKGROUND

Constraint-based models of Escherichia coli metabolic flux have played a key role in computational studies of cellular metabolism at the genome scale. We sought to develop a next-generation constraint-based E. coli model that achieved improved phenotypic prediction accuracy while being frequently updated and easy to use. We also sought to compare model predictions with experimental data to highlight open questions in E. coli biology.

RESULTS

We present EcoCyc-18.0-GEM, a genome-scale model of the E. coli K-12 MG1655 metabolic network. The model is automatically generated from the current state of EcoCyc using the MetaFlux software, enabling the release of multiple model updates per year. EcoCyc-18.0-GEM encompasses 1445 genes, 2286 unique metabolic reactions, and 1453 unique metabolites. We demonstrate a three-part validation of the model that breaks new ground in breadth and accuracy: (i) Comparison of simulated growth in aerobic and anaerobic glucose culture with experimental results from chemostat culture and simulation results from the E. coli modeling literature. (ii) Essentiality prediction for the 1445 genes represented in the model, in which EcoCyc-18.0-GEM achieves an improved accuracy of 95.2% in predicting the growth phenotype of experimental gene knockouts. (iii) Nutrient utilization predictions under 431 different media conditions, for which the model achieves an overall accuracy of 80.7%. The model's derivation from EcoCyc enables query and visualization via the EcoCyc website, facilitating model reuse and validation by inspection. We present an extensive investigation of disagreements between EcoCyc-18.0-GEM predictions and experimental data to highlight areas of interest to E. coli modelers and experimentalists, including 70 incorrect predictions of gene essentiality on glucose, 80 incorrect predictions of gene essentiality on glycerol, and 83 incorrect predictions of nutrient utilization.

CONCLUSION

Significant advantages can be derived from the combination of model organism databases and flux balance modeling represented by MetaFlux. Interpretation of the EcoCyc database as a flux balance model results in a highly accurate metabolic model and provides a rigorous consistency check for information stored in the database.

Nitrogen assimilation in Escherichia coli: putting molecular data into a systems perspective

We present a comprehensive overview of the hierarchical network of intracellular processes revolving around central nitrogen metabolism in Escherichia coli. The hierarchy intertwines transport, metabolism, signaling leading to posttranslational modification, and transcription. The protein components of the network include an ammonium transporter (AmtB), a glutamine transporter (GlnHPQ), two ammonium assimilation pathways (glutamine synthetase [GS]-glutamate synthase [glutamine 2-oxoglutarate amidotransferase {GOGAT}] and glutamate dehydrogenase [GDH]), the two bifunctional enzymes adenylyl transferase/adenylyl-removing enzyme (ATase) and uridylyl transferase/uridylyl-removing enzyme (UTase), the two trimeric signal transduction proteins (GlnB and GlnK), the two-component regulatory system composed of the histidine protein kinase nitrogen regulator II (NRII) and the response nitrogen regulator I (NRI), three global transcriptional regulators called nitrogen assimilation control (Nac) protein, leucine-responsive regulatory protein (Lrp), and cyclic AMP (cAMP) receptor protein (Crp), the glutaminases, and the nitrogen-phosphotransferase system. First, the structural and molecular knowledge on these proteins is reviewed. Thereafter, the activities of the components as they engage together in transport, metabolism, signal transduction, and transcription and their regulation are discussed. Next, old and new molecular data and physiological data are put into a common perspective on integral cellular functioning, especially with the aim of resolving counterintuitive or paradoxical processes featured in nitrogen assimilation. Finally, we articulate what still remains to be discovered and what general lessons can be learned from the vast amounts of data that are available now.

A vector library for silencing central carbon metabolism genes with antisense RNAs in Escherichia coli

We describe here the construction of a series of 71 vectors to silence central carbon metabolism genes in Escherichia coli. The vectors inducibly express antisense RNAs called paired-terminus antisense RNAs, which have a higher silencing efficacy than ordinary antisense RNAs. By measuring mRNA amounts, measuring activities of target proteins, or observing specific phenotypes, it was confirmed that all the vectors were able to silence the expression of target genes efficiently. Using this vector set, each of the central carbon metabolism genes was silenced individually, and the accumulation of metabolites was investigated. We were able to obtain accurate information on ways to increase the production of pyruvate, an industrially valuable compound, from the silencing results. Furthermore, the experimental results of pyruvate accumulation were compared to in silico predictions, and both sets of results were consistent. Compared to the gene disruption approach, the silencing approach has an advantage in that any E. coli strain can be used and multiple gene silencing is easily possible in any combination.

A comprehensive survey of small-molecule binding pockets in proteins

Many biological activities originate from interactions between small-molecule ligands and their protein targets. A detailed structural and physico-chemical characterization of these interactions could significantly deepen our understanding of protein function and facilitate drug design. Here, we present a large-scale study on a non-redundant set of about 20,000 known ligand-binding sites, or pockets, of proteins. We find that the structural space of protein pockets is crowded, likely complete, and may be represented by about 1,000 pocket shapes. Correspondingly, the growth rate of novel pockets deposited in the Protein Data Bank has been decreasing steadily over the recent years. Moreover, many protein pockets are promiscuous and interact with ligands of diverse scaffolds. Conversely, many ligands are promiscuous and interact with structurally different pockets. Through a physico-chemical and structural analysis, we provide insights into understanding both pocket promiscuity and ligand promiscuity. Finally, we discuss the implications of our study for the prediction of protein-ligand interactions based on pocket comparison.

The life of a living cell relies on many distinct proteins to carry out their functions. Most of these functions are rooted in interactions between the proteins and metabolites, small-molecules essential for life. By targeting specific proteins relevant to a disease, drug molecules may provide a cure. A deep understanding of the nature of interactions between proteins and small-molecules (or ligands) through analyzing their structures may help predict protein function or improve drug design. In this contribution, we present a large-scale analysis of a non-redundant set of over 20,000 experimental protein-ligand complex structures available in the current Protein Data Bank. We seek answers to several fundamental questions: How many representative pockets are there that serve as ligand-binding sites in proteins? To what extent can we infer a similar protein-ligand interaction by matching the structures of protein pockets? How different are the ligands found in the same pocket? For a promiscuous protein pocket, how does a pocket maintain favorable interactions with very different ligands? Conversely, how different are those pockets that interact with the same ligand? We find the structural space of protein pocket is small and that both protein promiscuity and ligand promiscuity are very common in Nature.

Determination of the structure of the catabolic N-succinylornithine transaminase (AstC) from Escherichia coli

Escherichia coli possesses two acyl ornithine aminotransferases, one catabolic (AstC) and the other anabolic (ArgD), that participate in L-arginine metabolism. Although only 58% identical, the enzymes have been shown to be functionally interchangeable. Here we have purified AstC and have obtained X-ray crystal structures of apo and holo-AstC and of the enzyme complexed with its physiological substrate, succinylornithine. We compare the structures obtained in this study with those of ArgD from Salmonella typhimurium obtained elsewhere, finding several notable differences. Docking studies were used to explore the docking modes of several substrates (ornithine, succinylornithine and acetylornithine) and the co-substrate glutamate/α-ketogluterate. The docking studies support our observations that AstC has a strong preference for acylated ornithine species over ornithine itself, and suggest that the increase in specificity associated with acylation is caused by steric and desolvation effects rather than specific interactions between the substrate and enzyme.

Three independent determinants of protein evolutionary rate

One of the most widely accepted ideas related to the evolutionary rates of proteins is that functionally important residues or regions evolve slower than other regions, a reasonable outcome of which should be a slower evolutionary rate of the proteins with a higher density of functionally important sites. Oddly, the role of functional importance, mainly measured by essentiality, in determining evolutionary rate has been challenged in recent studies. Several variables other than protein essentiality, such as expression level, gene compactness, protein-protein interactions, etc., have been suggested to affect protein evolutionary rate. In the present review, we try to refine the concept of functional importance of a gene, and consider three factors-functional importance, expression level, and gene compactness, as independent determinants of evolutionary rate of a protein, based not only on their known correlation with evolutionary rate but also on a reasonable mechanistic model. We suggest a framework based on these mechanistic models to correctly interpret the correlations between evolutionary rates and the various variables as well as the interrelationships among the variables.

The Entner-Doudoroff pathway empowers Pseudomonas putida KT2440 with a high tolerance to oxidative stress

Glucose catabolism of Pseudomonas putida is carried out exclusively through the Entner-Doudoroff (ED) pathway due to the absence of 6-phosphofructokinase. In order to activate the Embden-Meyerhof-Parnas (EMP) route we transferred the pfkA gene from Escherichia coli to a P. putida wild-type strain as well as to an eda mutant, i.e. lacking 2-keto-3-deoxy-6-phosphogluconate aldolase. PfkA(E. coli) failed to redirect the carbon flow from the ED route towards the EMP pathway, suggesting that ED was essential for sugar catabolism. The presence of PfkA(E. coli) was detrimental for growth, which could be traced to the reduction of ATP and NAD(P)H pools along with alteration of the NAD(P)H/NADP(+) ratio. Pseudomonas putida cells carrying PfkA(E. coli) became highly sensitive to diamide and hydrogen peroxide, the response to which is very demanding of NADPH. The inhibitory effect of PfkA(E. coli) could in part be relieved by methionine, the synthesis of which relies much on NADPH. These results expose the role of the ED pathway for generating the redox currency (NADPH) that is required for counteracting oxidative stress. It is thus likely that environmental bacteria that favour the ED pathway over the EMP pathway do so in order to gear their aerobic metabolism to endure oxidative-related insults.

Nutrient dependence of RNase E essentiality in Escherichia coli

Escherichia coli cells normally require RNase E activity to form colonies (colony-forming ability [CFA]). The CFA-defective phenotype of cells lacking RNase E is partly reversed by overexpression of the related endoribonuclease RNase G or by mutation of the gene encoding the RNA helicase DeaD. We found that the carbon source utilization by rne deaD doubly mutant bacteria differs from that of rne(+) cells and from that of cells mutated in deaD alone and that the loss of rne function in these bacteria limits conversion of the glycolytic pathway product phosphoenolpyruvate to the tricarboxylic acid (TCA) cycle intermediate oxaloacetic acid. We show that the mechanism underlying this effect is reduced production of the enzyme phosphoenolpyruvate carboxylase (PPC) and that adventitious overexpression of PPC, which facilitates phosphoenolpyruvate utilization and connects the glycolytic pathway with the TCA cycle, restored CFA to rne deaD mutant bacteria cultured on carbon sources that otherwise were unable to sustain growth. We further show that bacteria producing full-length RNase E, which allows formation of degradosomes, have nutritional requirements different from those of cells supplied with only the N-terminal catalytic region of RNase E and that mitigation of RNase E deficiency by overexpression of a related RNase, RNase G, is also affected by carbon source. Our results reveal previously unsuspected effects of RNase E deficiency and degradosome formation on nutrient utilization by E. coli cells.

Probing the mutational interplay between primary and promiscuous protein functions: a computational-experimental approach

Protein promiscuity is of considerable interest due its role in adaptive metabolic plasticity, its fundamental connection with molecular evolution and also because of its biotechnological applications. Current views on the relation between primary and promiscuous protein activities stem largely from laboratory evolution experiments aimed at increasing promiscuous activity levels. Here, on the other hand, we attempt to assess the main features of the simultaneous modulation of the primary and promiscuous functions during the course of natural evolution. The computational/experimental approach we propose for this task involves the following steps: a function-targeted, statistical coupling analysis of evolutionary data is used to determine a set of positions likely linked to the recruitment of a promiscuous activity for a new function; a combinatorial library of mutations on this set of positions is prepared and screened for both, the primary and the promiscuous activities; a partial-least-squares reconstruction of the full combinatorial space is carried out; finally, an approximation to the Pareto set of variants with optimal primary/promiscuous activities is derived. Application of the approach to the emergence of folding catalysis in thioredoxin scaffolds reveals an unanticipated scenario: diverse patterns of primary/promiscuous activity modulation are possible, including a moderate (but likely significant in a biological context) simultaneous enhancement of both activities. We show that this scenario can be most simply explained on the basis of the conformational diversity hypothesis, although alternative interpretations cannot be ruled out. Overall, the results reported may help clarify the mechanisms of the evolution of new functions. From a different viewpoint, the partial-least-squares-reconstruction/Pareto-set-prediction approach we have introduced provides the computational basis for an efficient directed-evolution protocol aimed at the simultaneous enhancement of several protein features and should therefore open new possibilities in the engineering of multi-functional enzymes.

Interpretations of evolutionary processes at the molecular level have been determined to a significant extent by the concept of “trade-off”, the idea that improving a given feature of a protein molecule by mutation will likely bring about deterioration in other features. For instance, if a protein is able to carry out two different molecular tasks based on the same functional site (competing tasks), optimization for one task could be naively expected to impair its performance for the other task. In this work, we report a computational/experimental approach to assess the potential patterns of modulation of two competing molecular tasks in the course of natural evolution. Contrary to the naïve expectation, we find that diverse modulation patterns are possible, including the simultaneous optimization of the two tasks. We show, however, that this simultaneous optimization is not in conflict with the trade-offs expected for two competing tasks: using the language of the theory of economic efficiency, trade-offs are realized in the Pareto set of optimal variants for the two tasks, while most protein variants do not belong to such Pareto set. That is, most protein variants are not Pareto-efficient and can potentially be improved in terms of several features.

High-throughput workflow for monitoring and mining bioprocess data and its application to inferring the physiological response of Escherichia coli to perturbations

Miniaturization and high-throughput screening are currently the focus of emerging research areas such as systems biology and systems biotechnology. A fluorescence-based screening assay for the online monitoring of oxygen and pH and a numerical method to mine the resulting online process data are described. The assay employs commercial phosphorescent oxygen- and pH-sensitive probes in standard 48- or 96-well plates on a plate reader equipped with a shaker. In addition to dual parametric analysis of both pH and oxygen in a single well, the assay allows monitoring of growth, as measured by absorbance. Validation of the assay is presented and compared with commercially available plates equipped with optical sensors for oxygen and pH. By using model-free fitting to the readily available online measurements, the length and rate of each phase such as the duration of lag and transition phase or acidification, growth, and oxygen consumption rates are automatically detected. In total, nine physiological descriptors, which can be used for further statistical and comparison analysis, are extracted from the pH, oxygen partial pressure (pO(2)), and optical density (OD) profiles. The combination of a simple mix-and-measure procedure with an automatic data mining method allows high sample throughput and good reproducibility while providing a physiological state identification and characterization of test cells. As a proof of concept, the utility of the workflow in assessing the physiological response of Escherichia coli to environmental and genetic perturbations is demonstrated.

Patterns of activities of root phosphomonoesterase and phosphodiesterase in wetland plants as a function of macrophyte species and ambient phosphorus regime

Phosphorus (P)-limited plants produce higher amounts of root phosphatases, but research has mostly focused on phosphomonoesterases (PMEs). Because phosphate diesters can form a significant proportion of organic P in wetlands, we aimed to determine whether wetland plants produce both root PMEs and root phosphodiesterases (PDEs), and, if so, what factors influence activities of these enzymes. We measured the activities of root PMEs and PDEs colorimetrically in a wide range of macrophytes from natural and P-enriched wetlands. Hydrolyzable P in sediments was analyzed using commercially available PMEs and PDEs. In all species, both root PMEs and PDEs were always present, and their activities were closely correlated. Sedges and broadleaved emergents had the highest activity of both enzymes, while those of floating-leaved plants were the lowest. Redundancy analysis revealed close association between root enzymes and the proportion of monoesterase- and diesterase-hydrolyzable dissolved unreactive P. Both enzymes were positively correlated with root tissue N : P ratio. Both plant and sediment traits were important when explaining differences in enzyme activities. Although the activities are related to ambient P regime, the relationship was not close enough to use root enzymes as reliable predictors of dissolved unreactive P that is hydrolyzed by sediment phosphomono- and diesterases.

Properties and identification of antibiotic drug targets

BACKGROUND

We analysed 48 non-redundant antibiotic target proteins from all bacteria, 22 antibiotic target proteins from E. coli only and 4243 non-drug targets from E. coli to identify differences in their properties and to predict new potential drug targets.

RESULTS

When compared to non-targets, bacterial antibiotic targets tend to be long, have high beta-sheet and low alpha-helix contents, are polar, are found in the cytoplasm rather than in membranes, and are usually enzymes, with ligases particularly favoured. Sequence features were used to build a support vector machine model for E. coli proteins, allowing the assignment of any sequence to the drug target or non-target classes, with an accuracy in the training set of 94%. We identified 319 proteins (7%) in the non-target set that have target-like properties, many of which have unknown function. 63 of these proteins have significant and undesirable similarity to a human protein, leaving 256 target like proteins that are not present in humans.

CONCLUSIONS

We suggest that antibiotic discovery programs would be more likely to succeed if new targets are chosen from this set of target like proteins or their homologues. In particular, 64 are essential genes where the cell is not able to recover from a random insertion disruption.

Goldbeter-Koshland model for open signaling cascades: a mathematical study

A modified Goldbeter-Koshland (GK) type model is proposed to study the open signaling cascade with forward activation. Three different regimes related to the ultrasensitivity are clearly defined. When the cascade is sufficiently long, the steady states in the downstream cycles exhibit a limiting behavior. Meanwhile, a temporal switch-like behavior is revealed in the pre-ultrasensitivity regime without any feedback circuitry. A detailed analysis of this phenomenon is provided.

Mapping human metabolic pathways in the small molecule chemical space

The work presented here is a study of human metabolic pathways, as projected in the chemical space of the small molecules they comprise, and it is composed of three parts: a) a study of the extent of clustering and overlap of these pathways in chemical space, b) the development and assessment of a statistical model for estimating the proximity to a given pathway of any small molecule, and c) the use of the above model in estimating the proximity of marketed drugs to human metabolic pathways. The distribution, overlap, and relationships of human metabolic pathways in this space are revealed using both visual and quantitative approaches. A set of selected physicochemical and topological descriptors is used to build a classifier, whose aim is to predict metabolic class and pathway membership of any small molecule. The classifier performs well for tightly clustered, isolated pathways but is, naturally, much less accurate for strongly overlapping pathways. Finally, the extent of overlap of a set of known drugs with the human metabolome is examined, and the classifier is used to predict likely cross-interactions between drugs and the major metabolic pathways in humans.

Systems-level approaches for identifying and analyzing genetic interaction networks in Escherichia coli and extensions to other prokaryotes

Molecular interactions define the functional organization of the cell. Epistatic (genetic, or gene-gene) interactions, one of the most informative and commonly encountered forms of functional relationships, are increasingly being used to map process architecture in model eukaryotic organisms. In particular, 'systems-level' screens in yeast and worm aimed at elucidating genetic interaction networks have led to the generation of models describing the global modular organization of gene products and protein complexes within a cell. However, comparable data for prokaryotic organisms have not been available. Given its ease of growth and genetic manipulation, the Gram-negative bacterium Escherichia coli appears to be an ideal model system for performing comprehensive genome-scale examinations of genetic redundancy in bacteria. In this review, we highlight emerging experimental and computational techniques that have been developed recently to examine functional relationships and redundancy in E. coli at a systems-level, and their potential application to prokaryotes in general. Additionally, we have scanned PubMed abstracts and full-text published articles to manually curate a list of approximately 200 previously reported synthetic sick or lethal genetic interactions in E. coli derived from small-scale experimental studies.

Fitness of Escherichia coli during urinary tract infection requires gluconeogenesis and the TCA cycle

Microbial pathogenesis studies traditionally encompass dissection of virulence properties such as the bacterium's ability to elaborate toxins, adhere to and invade host cells, cause tissue damage, or otherwise disrupt normal host immune and cellular functions. In contrast, bacterial metabolism during infection has only been recently appreciated to contribute to persistence as much as their virulence properties. In this study, we used comparative proteomics to investigate the expression of uropathogenic Escherichia coli (UPEC) cytoplasmic proteins during growth in the urinary tract environment and systematic disruption of central metabolic pathways to better understand bacterial metabolism during infection. Using two-dimensional fluorescence difference in gel electrophoresis (2D-DIGE) and tandem mass spectrometry, it was found that UPEC differentially expresses 84 cytoplasmic proteins between growth in LB medium and growth in human urine (P<0.005). Proteins induced during growth in urine included those involved in the import of short peptides and enzymes required for the transport and catabolism of sialic acid, gluconate, and the pentose sugars xylose and arabinose. Proteins required for the biosynthesis of arginine and serine along with the enzyme agmatinase that is used to produce the polyamine putrescine were also up-regulated in urine. To complement these data, we constructed mutants in these genes and created mutants defective in each central metabolic pathway and tested the relative fitness of these UPEC mutants in vivo in an infection model. Import of peptides, gluconeogenesis, and the tricarboxylic acid cycle are required for E. coli fitness during urinary tract infection while glycolysis, both the non-oxidative and oxidative branches of the pentose phosphate pathway, and the Entner-Doudoroff pathway were dispensable in vivo. These findings suggest that peptides and amino acids are the primary carbon source for E. coli during infection of the urinary tract. Because anaplerosis, or using central pathways to replenish metabolic intermediates, is required for UPEC fitness in vivo, we propose that central metabolic pathways of bacteria could be considered critical components of virulence for pathogenic microbes.

Bacteria that cause infections often have genes known as virulence factors that are required for bacteria to cause disease. Studying virulence factors such as toxins, adhesins, and secretion and iron-acquisition systems is a fundamental part of understanding infectious disease mechanisms. In contrast, little is known about the contribution of bacterial metabolism to infectious disease. This study shows that E. coli, which cause most urinary tract infections, utilize peptides as a preferred carbon source in vivo and requires some, but not all, of the central metabolic pathways to infect the urinary tract. Specifically, pathways that can be used to replenish metabolites, known as anaplerotic reactions, are important for uropathogenic E. coli infections. These findings help explain how metabolism can contribute to the ability of bacteria to cause a common infection.