Genome-scale reconstruction of metabolic network in Bacillus subtilis based on high-throughput phenotyping and gene essentiality data

Genome-scale metabolic modeling of P. thermoglucosidasius NCIMB 11955 reveals metabolic bottlenecks in anaerobic metabolism

Parageobacillus thermoglucosidasius represents a thermophilic, facultative anaerobic bacterial chassis, with several desirable traits for metabolic engineering and industrial production. To further optimize strain productivity, a systems level understanding of its metabolism is needed, which can be facilitated by a genome-scale metabolic model. Here, we present p-thermo, the most complete, curated and validated genome-scale model (to date) of Parageobacillus thermoglucosidasius NCIMB 11955. It spans a total of 890 metabolites, 1175 reactions and 917 metabolic genes, forming an extensive knowledge base for P. thermoglucosidasius NCIMB 11955 metabolism. The model accurately predicts aerobic utilization of 22 carbon sources, and the predictive quality of internal fluxes was validated with previously published ¹³C-fluxomics data. In an application case, p-thermo was used to facilitate more in-depth analysis of reported metabolic engineering efforts, giving additional insight into fermentative metabolism. Finally, p-thermo was used to resolve a previously uncharacterised bottleneck in anaerobic metabolism, by identifying the minimal required supplemented nutrients (thiamin, biotin and iron(III)) needed to sustain anaerobic growth. This highlights the usefulness of p-thermo for guiding the generation of experimental hypotheses and for facilitating data-driven metabolic engineering, expanding the use of P. thermoglucosidasius as a high yield production platform.

High-Quality Genome-Scale Models From Error-Prone, Long-Read Assemblies

Advances in nanopore-based sequencing techniques have enabled rapid characterization of genomes and transcriptomes. An emerging application of this sequencing technology is point-of-care characterization of pathogenic bacteria. However, genome assessments alone are unable to provide a complete understanding of the pathogenic phenotype. Genome-scale metabolic reconstruction and analysis is a bottom-up Systems Biology technique that has elucidated the phenotypic nuances of antimicrobial resistant (AMR) bacteria and other human pathogens. Combining these genome-scale models (GEMs) with point-of-care nanopore sequencing is a promising strategy for combating the emerging health challenge of AMR pathogens. However, the sequencing errors inherent to the nanopore technique may negatively affect the quality, and therefore the utility, of GEMs reconstructed from nanopore assemblies. Here we describe and validate a workflow for rapid construction of GEMs from nanopore (MinION) derived assemblies. Benchmarking the pipeline against a high-quality reference GEM of Escherichia coli K-12 resulted in nanopore-derived models that were >99% complete even at sequencing depths of less than 10× coverage. Applying the pipeline to clinical isolates of pathogenic bacteria resulted in strain-specific GEMs that identified canonical AMR genome content and enabled simulations of strain-specific microbial growth. Additionally, we show that treating the sequencing run as a mock metagenome did not degrade the quality of models derived from metagenome assemblies. Taken together, this study demonstrates that combining nanopore sequencing with GEM construction pipelines enables rapid, in situ characterization of microbial metabolism.

Essential gene prediction using limited gene essentiality information-An integrative semi-supervised machine learning strategy

Essential gene prediction helps to find minimal genes indispensable for the survival of any organism. Machine learning (ML) algorithms have been useful for the prediction of gene essentiality. However, currently available ML pipelines perform poorly for organisms with limited experimental data. The objective is the development of a new ML pipeline to help in the annotation of essential genes of less explored disease-causing organisms for which minimal experimental data is available. The proposed strategy combines unsupervised feature selection technique, dimension reduction using the Kamada-Kawai algorithm, and semi-supervised ML algorithm employing Laplacian Support Vector Machine (LapSVM) for prediction of essential and non-essential genes from genome-scale metabolic networks using very limited labeled dataset. A novel scoring technique, Semi-Supervised Model Selection Score, equivalent to area under the ROC curve (auROC), has been proposed for the selection of the best model when supervised performance metrics calculation is difficult due to lack of data. The unsupervised feature selection followed by dimension reduction helped to observe a distinct circular pattern in the clustering of essential and non-essential genes. LapSVM then created a curve that dissected this circle for the classification and prediction of essential genes with high accuracy (auROC > 0.85) even with 1% labeled data for model training. After successful validation of this ML pipeline on both Eukaryotes and Prokaryotes that show high accuracy even when the labeled dataset is very limited, this strategy is used for the prediction of essential genes of organisms with inadequate experimentally known data, such as Leishmania sp. Using a graph-based semi-supervised machine learning scheme, a novel integrative approach has been proposed for essential gene prediction that shows universality in application to both Prokaryotes and Eukaryotes with limited labeled data. The essential genes predicted using the pipeline provide an important lead for the prediction of gene essentiality and identification of novel therapeutic targets for antibiotic and vaccine development against disease-causing parasites.

Diagnosing and Predicting Mixed-Culture Fermentations with Unicellular and Guild-Based Metabolic Models

Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

Multispecies microbial communities determine the fate of materials in the environment and can be harnessed to produce beneficial products from renewable resources. In a recent example, fermentations by microbial communities have produced medium-chain fatty acids (MCFAs). Tools to predict, assess, and improve the performance of these communities, however, are limited. To provide such tools, we constructed two metabolic models of MCFA-producing microbial communities based on available genomic, transcriptomic, and metabolomic data. The first model is a unicellular model (iFermCell215), while the second model (iFermGuilds789) separates fermentation activities into functional guilds. Ethanol and lactate are fermentation products known to serve as substrates for MCFA production, while acetate is another common cometabolite during MCFA production. Simulations with iFermCell215 predict that low molar ratios of acetate to ethanol favor MCFA production, whereas the products of lactate and acetate coutilization are less dependent on the acetate-to-lactate ratio. In simulations of an MCFA-producing community fed a complex organic mixture derived from lignocellulose, iFermGuilds789 predicted that lactate was an extracellular cometabolite that served as a substrate for butyrate (C4) production. Extracellular hexanoic (C6) and octanoic (C8) acids were predicted by iFermGuilds789 to be from community members that directly metabolize sugars. Modeling results provide several hypotheses that can improve our understanding of microbial roles in an MCFA-producing microbiome and inform strategies to increase MCFA production. Further, these models represent novel tools for exploring the role of mixed microbial communities in carbon recycling in the environment, as well as in beneficial reuse of organic residues.

IMPORTANCE Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

A systematic strategy using a reconstructed genome-scale metabolic network for pathogen Streptococcuspneumoniae D39 to find novel potential drug targets

Streptococcus pneumoniae is a Gram-positive bacterium that is one of the major causes of various infections such as pneumonia, meningitis, otitis media and endocarditis. Since antibiotic resistance of S. pneumoniae is pointed out as a challenge in the treatment of these infections, more studies are required to focus on disease prevention. In this research, a first manually curated genome-scale metabolic network of the pathogen S. pneumoniae D39 was reconstructed based on its genome annotation data, and biochemical knowledge from literature and databases. The model was validated by amino acid auxotrophies, gene essentiality analysis, and different carbohydrate sources. Then, a two-stage strategy was developed to find target genes for growth reduction of the pathogen and their importance in the various infection sites. In the first stage, growth-associated genes were identified by integration of transcriptomic data with the model and in the second stage, the importance of each gene in the metabolism for growth was evaluated using principal component analysis. The reports presented in the literature confirm the effect of some found genes on the growth of S. pneumoniae.

Charting the Metabolic Landscape of the Facultative Methylotroph Bacillus methanolicus

Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Bacillus methanolicus MGA3 is a thermotolerant and relatively fast-growing methylotroph able to secrete large quantities of glutamate and lysine. These natural characteristics make B. methanolicus a good candidate to become a new industrial chassis organism, especially in a methanol-based economy. Intriguingly, the only substrates known to support B. methanolicus growth as sole sources of carbon and energy are methanol, mannitol, and, to a lesser extent, glucose and arabitol. Because fluxomics provides the most direct readout of the cellular phenotype, we hypothesized that comparing methylotrophic and nonmethylotrophic metabolic states at the flux level would yield new insights into MGA3 metabolism. In this study, we designed and performed a ¹³C metabolic flux analysis (¹³C-MFA) of the facultative methylotroph B. methanolicus MGA3 growing on methanol, mannitol, and arabitol to compare the associated metabolic states. On methanol, results showed a greater flux in the ribulose monophosphate (RuMP) pathway than in the tricarboxylic acid (TCA) cycle, thus validating previous findings on the methylotrophy of B. methanolicus. New insights related to the utilization of cyclic RuMP versus linear dissimilation pathways and between the RuMP variants were generated. Importantly, we demonstrated that the linear detoxification pathways and the malic enzyme shared with the pentose phosphate pathway have an important role in cofactor regeneration. Finally, we identified, for the first time, the metabolic pathway used to assimilate arabitol. Overall, those data provide a better understanding of this strain under various environmental conditions.

IMPORTANCE Methanol is inexpensive, is easy to transport, and can be produced both from renewable and from fossil resources without mobilizing arable lands. As such, it is regarded as a potential carbon source to transition toward a greener industrial chemistry. Metabolic engineering of bacteria and yeast able to efficiently consume methanol is expected to provide cell factories that will transform methanol into higher-value chemicals in the so-called methanol economy. Toward that goal, the study of natural methylotrophs such as Bacillus methanolicus is critical to understand the origin of their efficient methylotrophy. This knowledge will then be leveraged to transform such natural strains into new cell factories or to design methylotrophic capability in other strains already used by the industry.

Assessment of transcriptomic constraint-based methods for central carbon flux inference

MOTIVATION

Determining intracellular metabolic flux through isotope labeling techniques such as 13C metabolic flux analysis (13C-MFA) incurs significant cost and effort. Previous studies have shown transcriptomic data coupled with constraint-based metabolic modeling can determine intracellular fluxes that correlate highly with 13C-MFA measured fluxes and can achieve higher accuracy than constraint-based metabolic modeling alone. These studies, however, used validation data limited to E. coli and S. cerevisiae grown on glucose, with significantly similar flux distribution for central metabolism. It is unclear whether those results apply to more diverse metabolisms, and therefore further, extensive validation is needed.

RESULTS

In this paper, we formed a dataset of transcriptomic data coupled with corresponding 13C-MFA flux data for 21 experimental conditions in different unicellular organisms grown on varying carbon substrates and conditions. Three computational flux-balance analysis (FBA) methods were comparatively assessed. The results show when uptake rates of carbon sources and key metabolites are known, transcriptomic data provides no significant advantage over constraint-based metabolic modeling (average correlation coefficients, transcriptomic E-Flux2 0.725 and SPOT 0.650 vs non-transcriptomic pFBA 0.768). When uptake rates are unknown, however, predictions obtained utilizing transcriptomic data are generally good and significantly better than those obtained using constraint-based metabolic modeling alone (E-Flux2 0.385 and SPOT 0.583 vs pFBA 0.237). Thus, transcriptomic data coupled with constraint-based metabolic modeling is a promising method to obtain intracellular flux estimates in microorganisms, particularly in cases where uptake rates of key metabolites cannot be easily determined, such as for growth in complex media or in vivo conditions.

Cell-free synthesis system-assisted pathway bottleneck diagnosis and engineering in Bacillus subtilis

Metabolic engineering is a key technology for cell factories construction by rewiring cellular resources to achieve efficient production of target chemicals. However, the existence of bottlenecks in synthetic pathway can seriously affect production efficiency, which is also one of the core issues for metabolic engineers to solve. Therefore, developing an approach for diagnosing potential metabolic bottlenecks in a faster and simpler manner is of great significance to accelerate cell factories construction. The cell-free reaction system based on cell lysates can transfer metabolic reactions from in vivo to in vitro, providing a flexible access to directly change protein and metabolite variables, thus provides a potential solution for rapid identification of bottlenecks. Here, bottleneck diagnosis of the N-acetylneuraminic acid (NeuAc) biosynthesis pathway in industrially important chassis microorganism Bacillus subtilis was performed using cell-free synthesis system. Specifically, a highly efficient B. subtilis cell-free system for NeuAc de novo synthesis was firstly constructed, which had a 305-fold NeuAc synthesis rate than that in vivo and enabled fast pathway dynamics analysis. Next, through the addition of all potential key intermediates in combination with substrate glucose respectively, it was found that insufficient phosphoenolpyruvate supply was one of the NeuAc pathway bottlenecks. Rational in vivo metabolic engineering of NeuAc-producing B. subtilis was further performed to eliminate the bottleneck. By down-regulating the expression level of pyruvate kinase throughout the growth phase or only in the stationary phase using inhibitory N-terminal coding sequences (NCSs) and growth-dependent regulatory NCSs respectively, the maximal NeuAc titer increased 2.0-fold. Our study provides a rapid method for bottleneck diagnosis, which may help to accelerate the cycle of design, build, test and learn cycle for metabolic engineering.

Genome scale metabolic models and analysis for evaluating probiotic potentials

Probiotics are live beneficial microorganisms that can be consumed in the form of dairy and food products as well as dietary supplements to promote a healthy balance of gut bacteria in humans. Practically, the main challenge is to identify and select promising strains and formulate multi-strain probiotic blends with consistent efficacy which is highly dependent on individual dietary regimes, gut environments, and health conditions. Limitations of current in vivo and in vitro methods for testing probiotic strains can be overcome by in silico model guided systems biology approaches where genome scale metabolic models (GEMs) can be used to describe their cellular behaviors and metabolic states of probiotic strains under various gut environments. Here, we summarize currently available GEMs of microbial strains with probiotic potentials and propose a knowledge-based framework to evaluate metabolic capabilities on the basis of six probiotic criteria. They include metabolic characteristics, stability, safety, colonization, postbiotics, and interaction with the gut microbiome which can be assessed by in silico approaches. As such, the most suitable strains can be identified to design personalized multi-strain probiotics in the future.

Genome-scale reconstruction of Paenarthrobacter aurescens TC1 metabolic model towards the study of atrazine bioremediation

Atrazine is an herbicide and a pollutant of great environmental concern that is naturally biodegraded by microbial communities. Paenarthrobacter aurescens TC1 is one of the most studied degraders of this herbicide. Here, we developed a genome scale metabolic model for P. aurescens TC1, iRZ1179, to study the atrazine degradation process at organism level. Constraint based flux balance analysis and time dependent simulations were used to explore the organism’s phenotypic landscape. Simulations aimed at designing media optimized for supporting growth and enhancing degradation, by passing the need in strain design via genetic modifications. Growth and degradation simulations were carried with more than 100 compounds consumed by P. aurescens TC1. In vitro validation confirmed the predicted classification of different compounds as efficient, moderate or poor stimulators of growth. Simulations successfully captured previous reports on the use of glucose and phosphate as bio-stimulators of atrazine degradation, supported by in vitro validation. Model predictions can go beyond supplementing the medium with a single compound and can predict the growth outcomes for higher complexity combinations. Hence, the analysis demonstrates that the exhaustive power of the genome scale metabolic reconstruction allows capturing complexities that are beyond common biochemical expertise and knowledge and further support the importance of computational platforms for the educated design of complex media. The model presented here can potentially serve as a predictive tool towards achieving optimal biodegradation efficiencies and for the development of ecologically friendly solutions for pollutant degradation.

Synthetic microbial communities of heterotrophs and phototrophs facilitate sustainable growth

Microbial communities comprised of phototrophs and heterotrophs hold great promise for sustainable biotechnology. Successful application of these communities relies on the selection of appropriate partners. Here we construct four community metabolic models to guide strain selection, pairing phototrophic, sucrose-secreting Synechococcus elongatus with heterotrophic Escherichia coli K-12, Escherichia coli W, Yarrowia lipolytica, or Bacillus subtilis. Model simulations reveae metabolic exchanges that sustain the heterotrophs in minimal media devoid of any organic carbon source, pointing to S. elongatus-E. coli K-12 as the most active community. Experimental validation of flux predictions for this pair confirms metabolic interactions and potential production capabilities. Synthetic communities bypass member-specific metabolic bottlenecks (e.g. histidine- and transport-related reactions) and compensate for lethal genetic traits, achieving up to 27% recovery from lethal knockouts. The study provides a robust modelling framework for the rational design of synthetic communities with optimized growth sustainability using phototrophic partners.

A Genome-Scale Metabolic Model of 2,3-Butanediol Production by Thermophilic Bacteria Geobacillus icigianus

The thermophilic strain of the genus Geobacillus, Geobacillus icigianus is a promising bacterial chassis for a wide range of biotechnological applications. In this study, we explored the metabolic potential of Geobacillus icigianus for the production of 2,3-butanediol (2,3-BTD), one of the cost-effective commodity chemicals. Here we present a genome-scale metabolic model iMK1321 for Geobacillus icigianus constructed using an auto-generating pipeline with consequent thorough manual curation. The model contains 1321 genes and includes 1676 reactions and 1589 metabolites, representing the most-complete and publicly available model of the genus Geobacillus. The developed model provides new insights into thermophilic bacterial metabolism and highlights new strategies for biotechnological applications of the strain. Our analysis suggests that Geobacillus icigianus has a potential for 2,3-butanediol production from a variety of utilized carbon sources, including glycerine, a common byproduct of biofuel production. We identified a set of solutions for enhancing 2,3-BTD production, including cultivation under anaerobic or microaerophilic conditions and decreasing the TCA flux to succinate via reducing citrate synthase activity. Both in silico predicted metabolic alternatives have been previously experimentally verified for closely related strains including the genus Bacillus.

Metabolic cell communication within tumour microenvironment: models, methods and perspectives

Environmental cues are essential in defining tumour malignancy, by promoting tumour initiation, progression and metastatic spreading. Stromal cells may metabolically cooperate or compete with cancer cells, playing a mandatory role in defining cancer metabolic plasticity, potentially dictating the final tumour outcome. Assessing shared nutrients between different tumoural or stromal compartments is essential to understand the impact of environmental nutrients on the metabolic plasticity of tumours. Here, we review analytical and computational approaches for studying the tumour metabolic microenvironment, the destiny of nutrients shared among tumour and stromal populations, as well as the molecular modules of these metabolic relationships.

In vitro synergistic effects of three enzymes from Bacillus subtilis CH-1 on keratin decomposition

Extracellular protease Vpr (Vpr), gamma-glutamyltranspeptidase (GGT; EC 2.3.2.2) and glyoxal/methylglyoxal reductase (YvgN; EC 1.1.1.21) are extracellular enzymes involved in feather degradation, which were identified by secretome analyses from an efficient feather-degrading strain Bacillus subtilis CH-1. The encoding sequences corresponding to the three secretory enzymes were cloned into vector pET22b for recombinant expression in Escherichia coli strain BL21 (DE3). Afterward, the proteins containing the C-terminal His-tag were purified using a Ni-IDA column. The optimal temperatures and pH values for protease activity of recombinant Vpr, GGT, and YvgN were identified as 45 °C/pH 7.0, 40 °C/pH 8.0, and 50 °C/pH 6.0 respectively when casein is the substrate. Furthermore, the synergistic effects of the three enzymes were studied using feather powder as substrate. Vpr was the core enzyme to hydrolyze keratin, while GGT and YvgN were coenzymes providing reducing activities for keratin decomposition. The keratinolytic activity was enhanced to about 1.4-folds when YvgN and Vpr applied together in comparison to Vpr alone. And the keratinolytic activity almost reached to 1.5-folds when all the three enzymes were combined to use. The study provides a novel perspective of the mechanism of keratin degradation by microorganisms, and thereby may also be of relevance for the design of an industrial process for enzymatic keratin degradation; however, additional experiments must be done to substantiate this conclusion.

Genome-Scale Investigation of the Metabolic Determinants Generating Bacterial Fastidious Growth

Xylella fastidiosa is one of the most important threats to plant health worldwide, causing disease in the Americas on a range of agricultural crops and trees, and recently associated with a critical epidemic affecting olive trees in Europe. A main challenge for the detection of the pathogen and the development of physiological studies is its fastidious growth, as the generation time can vary from 10 to 100 h for some strains. This physiological peculiarity is shared with several human pathogens and is poorly understood. We performed an analysis of the metabolic capabilities of X. fastidiosa through a genome-scale metabolic model of the bacterium. This model was reconstructed and manually curated using experiments and bibliographical evidence. Our study revealed that fastidious growth most probably results from different metabolic specificities such as the absence of highly efficient enzymes or a global inefficiency in virulence factor production. These results support the idea that the fragility of the metabolic network may have been shaped during evolution to lead to the self-limiting behavior of X. fastidiosa.

High proliferation rate and robustness are vital characteristics of bacterial pathogens that successfully colonize their hosts. The observation of drastically slow growth in some pathogens is thus paradoxical and remains unexplained. In this study, we sought to understand the slow (fastidious) growth of the plant pathogen Xylella fastidiosa. Using genome-scale metabolic network reconstruction, modeling, and experimental validation, we explored its metabolic capabilities. Despite genome reduction and slow growth, the pathogen’s metabolic network is complete but strikingly minimalist and lacking in robustness. Most alternative reactions were missing, especially those favoring fast growth, and were replaced by less efficient paths. We also found that the production of some virulence factors imposes a heavy burden on growth. Interestingly, some specific determinants of fastidious growth were also found in other slow-growing pathogens, enriching the view that these metabolic peculiarities are a pathogenicity strategy to remain at a low population level.

IMPORTANCEXylella fastidiosa is one of the most important threats to plant health worldwide, causing disease in the Americas on a range of agricultural crops and trees, and recently associated with a critical epidemic affecting olive trees in Europe. A main challenge for the detection of the pathogen and the development of physiological studies is its fastidious growth, as the generation time can vary from 10 to 100 h for some strains. This physiological peculiarity is shared with several human pathogens and is poorly understood. We performed an analysis of the metabolic capabilities of X. fastidiosa through a genome-scale metabolic model of the bacterium. This model was reconstructed and manually curated using experiments and bibliographical evidence. Our study revealed that fastidious growth most probably results from different metabolic specificities such as the absence of highly efficient enzymes or a global inefficiency in virulence factor production. These results support the idea that the fragility of the metabolic network may have been shaped during evolution to lead to the self-limiting behavior of X. fastidiosa.

Suggestions for Standardized Identifiers for Fatty Acyl Compounds in Genome Scale Metabolic Models and Their Application to the WormJam Caenorhabditis elegans Model

Genome scale metabolic models (GSMs) are a representation of the current knowledge on the metabolism of a given organism or superorganism. They group metabolites, genes, enzymes and reactions together to form a mathematical model and representation that can be used to analyze metabolic networks in silico or used for analysis of omics data. Beside correct mass and charge balance, correct structural annotation of metabolites represents an important factor for analysis of these metabolic networks. However, several metabolites in different GSMs have no or only partial structural information associated with them. Here, a new systematic nomenclature for acyl-based metabolites such as fatty acids, acyl-carnitines, acyl-coenzymes A or acyl-carrier proteins is presented. This nomenclature enables one to encode structural details in the metabolite identifiers and improves human readability of reactions. As proof of principle, it was applied to the fatty acid biosynthesis and degradation in the Caenorhabditis elegans consensus model WormJam.

Genome-scale metabolic models of Microbacterium species isolated from a high altitude desert environment

The Atacama Desert is the most arid desert on Earth, focus of important research activities related to microbial biodiversity studies. In this context, metabolic characterization of arid soil bacteria is crucial to understand their survival strategies under extreme environmental stress. We investigated whether strain-specific features of two Microbacterium species were involved in the metabolic ability to tolerate/adapt to local variations within an extreme desert environment. Using an integrative systems biology approach we have carried out construction and comparison of genome-scale metabolic models (GEMs) of two Microbacterium sp., CGR1 and CGR2, previously isolated from physicochemically contrasting soil sites in the Atacama Desert. Despite CGR1 and CGR2 belong to different phylogenetic clades, metabolic pathways and attributes are highly conserved in both strains. However, comparison of the GEMs showed significant differences in the connectivity of specific metabolites related to pH tolerance and CO₂ production. The latter is most likely required to handle acidic stress through decarboxylation reactions. We observed greater GEM connectivity within Microbacterium sp. CGR1 compared to CGR2, which is correlated with the capacity of CGR1 to tolerate a wider pH tolerance range. Both metabolic models predict the synthesis of pigment metabolites (β-carotene), observation validated by HPLC experiments. Our study provides a valuable resource to further investigate global metabolic adaptations of bacterial species to grow in soils with different abiotic factors within an extreme environment.

Comparison of Isomerase and Weimberg Pathway for γ-PGA Production From Xylose by Engineered Bacillus subtilis

The production of poly-γ-glutamic acid (γ-PGA), a biopolymer consisting of D- and L-glutamic acid monomers, currently relies on L-glutamate, or citrate as carbon substrates. Here we aimed at using plant biomass-derived substrates such as xylose. γ-PGA producing microorganisms including Bacillus subtilis natively metabolize xylose via the isomerase pathway. The Weimberg pathway, a xylose utilization pathway first described for Caulobacter crescentus, offers a carbon-efficient alternative converting xylose to 2-oxoglutarate without carbon loss. We engineered a recombinant B. subtilis strain that was able to grow on xylose with a growth rate of 0.43 h^-1 using a recombinant Weimberg pathway. Although ion-pair reversed-phase LC/MS/MS metabolome analysis revealed lower concentrations of γ-PGA precursors such as 2-oxoglutarate, the γ-PGA titer was increased 6-fold compared to the native xylose isomerase strain. Further metabolome analysis indicates a metabolic bottleneck in the phosphoenolpyruvate-pyruvate-oxaloacetate node causing bi-phasic (diauxic) growth of the recombinant Weimberg strain. Flux balance analysis (FBA) of the γ-PGA producing B. subtilis indicated that a maximal theoretical γ-PGA yield is achieved on D-xylose/ D-glucose mixtures. The results of the B. subtilis strain harboring the Weimberg pathway on such D-xylose/ D-glucose mixtures demonstrate indeed resource efficient, high yield γ-PGA production from biomass-derived substrates.

Reinforcement Learning for Bioretrosynthesis

Metabolic engineering aims to produce chemicals of interest from living organisms, to advance toward greener chemistry. Despite efforts, the research and development process is still long and costly, and efficient computational design tools are required to explore the chemical biosynthetic space. Here, we propose to explore the bioretrosynthesis space using an artificial intelligence based approach relying on the Monte Carlo Tree Search reinforcement learning method, guided by chemical similarity. We implement this method in RetroPath RL, an open-source and modular command line tool. We validate it on a golden data set of 20 manually curated experimental pathways as well as on a larger data set of 152 successful metabolic engineering projects. Moreover, we provide a novel feature that suggests potential media supplements to complement the enzymatic synthesis plan.

The functional repertoire contained within the native microbiota of the model nematode Caenorhabditis elegans

The microbiota is generally assumed to have a substantial influence on the biology of multicellular organisms. The exact functional contributions of the microbes are often unclear and cannot be inferred easily from 16S rRNA genotyping, which is commonly used for taxonomic characterization of bacterial associates. In order to bridge this knowledge gap, we here analyzed the metabolic competences of the native microbiota of the model nematode Caenorhabditis elegans. We integrated whole-genome sequences of 77 bacterial microbiota members with metabolic modeling and experimental characterization of bacterial physiology. We found that, as a community, the microbiota can synthesize all essential nutrients for C. elegans. Both metabolic models and experimental analyses revealed that nutrient context can influence how bacteria interact within the microbiota. We identified key bacterial traits that are likely to influence the microbe's ability to colonize C. elegans (i.e., the ability of bacteria for pyruvate fermentation to acetoin) and affect nematode fitness (i.e., bacterial competence for hydroxyproline degradation). Considering that the microbiota is usually neglected in C. elegans research, the resource presented here will help our understanding of this nematode's biology in a more natural context. Our integrative approach moreover provides a novel, general framework to characterize microbiota-mediated functions.

Manually curated genome-scale reconstruction of the metabolic network of Bacillus megaterium DSM319

Bacillus megaterium is a microorganism widely used in industrial biotechnology for production of enzymes and recombinant proteins, as well as in bioleaching processes. Precise understanding of its metabolism is essential for designing engineering strategies to further optimize B. megaterium for biotechnology applications. Here, we present a genome-scale metabolic model for B. megaterium DSM319, iJA1121, which is a result of a metabolic network reconciliation process. The model includes 1709 reactions, 1349 metabolites, and 1121 genes. Based on multiple-genome alignments and available genome-scale metabolic models for other Bacillus species, we constructed a draft network using an automated approach followed by manual curation. The refinements were performed using a gap-filling process. Constraint-based modeling was used to scrutinize network features. Phenotyping assays were performed in order to validate the growth behavior of the model using different substrates. To verify the model accuracy, experimental data reported in the literature (growth behavior patterns, metabolite production capabilities, metabolic flux analysis using ¹³C glucose and formaldehyde inhibitory effect) were confronted with model predictions. This indicated a very good agreement between in silico results and experimental data. For example, our in silico study of fatty acid biosynthesis and lipid accumulation in B. megaterium highlighted the importance of adopting appropriate carbon sources for fermentation purposes. We conclude that the genome-scale metabolic model iJA1121 represents a useful tool for systems analysis and furthers our understanding of the metabolism of B. megaterium.

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.

Influence of the Crc global regulator on substrate uptake rates and the distribution of metabolic fluxes in Pseudomonas putida KT2440 growing in a complete medium

When the soil bacterium Pseudomonas putida grows in a complete medium, it prioritizes the assimilation of preferred carbon sources, optimizing its metabolism and growth. This regulatory process is orchestrated by the Crc and Hfq proteins. The present work examines the changes that occur in metabolic fluxes when the crc gene is inactivated and cells grow exponentially in LB complete medium. Analyses were performed at three different moments during exponential growth, examining the assimilation rates for the compounds present in LB, changes in the proteome, and the changes in metabolic fluxes predicted by the iJN1411 metabolic model for P. putida KT2440. During the early exponential phase, consumption rates for sugars, many organic acids and most amino acids were higher in a Crc-null strain than in the wild type, leading to an overflow of the metabolic pathways and the leakage of pyruvate and acetate. These accelerated consumption rates decreased during the mid-exponential phase, when cells mostly used sugars and alanine. At later times, pyruvate was recovered from the medium and utilized. The higher consumption rates of the Crc-null strain reduced the growth rate. The lack of the Crc/Hfq regulatory system thus led to unbalanced metabolism with poorly optimized metabolic fluxes.

Metabolic Modeling of Streptococcus mutans Reveals Complex Nutrient Requirements of an Oral Pathogen

Streptococcus mutans is a Gram-positive bacterium that thrives under acidic conditions and is a primary cause of tooth decay (dental caries). To better understand the metabolism of S. mutans on a systematic level, we manually constructed a genome-scale metabolic model of the S. mutans type strain UA159. The model, called iSMU, contains 675 reactions involving 429 metabolites and the products of 493 genes. We validated iSMU by comparing simulations with growth experiments in defined medium. The model simulations matched experimental results for 17 of 18 carbon source utilization assays and 47 of 49 nutrient depletion assays. We also simulated the effects of single gene deletions. The model's predictions agreed with 78.1% and 84.4% of the gene essentiality predictions from two experimental data sets. Our manually curated model is more accurate than S. mutans models generated from automated reconstruction pipelines and more complete than other manually curated models. We used iSMU to generate hypotheses about the S. mutans metabolic network. Subsequent genetic experiments confirmed that (i) S. mutans catabolizes sorbitol via a sorbitol-6-phosphate 2-dehydrogenase (SMU_308) and (ii) the Leloir pathway is required for growth on complex carbohydrates such as raffinose. We believe the iSMU model is an important resource for understanding the metabolism of S. mutans and guiding future experiments.IMPORTANCE Tooth decay is the most prevalent chronic disease in the United States. Decay is caused by the bacterium Streptococcus mutans, an oral pathogen that ferments sugars into tooth-destroying lactic acid. We constructed a complete metabolic model of S. mutans to systematically investigate how the bacterium grows. The model provides a valuable resource for understanding and targeting S. mutans' ability to outcompete other species in the oral microbiome.

Metabolic engineering of Bacillus amyloliquefaciens LL3 for enhanced poly-γ-glutamic acid synthesis

Poly-γ-glutamic acid (γ-PGA) is a biocompatible and biodegradable polypeptide with wide-ranging applications in foods, cosmetics, medicine, agriculture and wastewater treatment. Bacillus amyloliquefaciens LL3 can produce γ-PGA from sucrose that can be obtained easily from sugarcane and sugar beet. In our previous work, it was found that low intracellular glutamate concentration was the limiting factor for γ-PGA production by LL3. In this study, the γ-PGA synthesis by strain LL3 was enhanced by chromosomally engineering its glutamate metabolism-relevant networks. First, the downstream metabolic pathways were partly blocked by deleting fadR, lysC, aspB, pckA, proAB, rocG and gudB. The resulting strain NK-A6 synthesized 4.84 g l-1 γ-PGA, with a 31.5% increase compared with strain LL3. Second, a strong promoter PC 2up was inserted into the upstream of icd gene, to generate strain NK-A7, which further led to a 33.5% improvement in the γ-PGA titre, achieving 6.46 g l-1 . The NADPH level was improved by regulating the expression of pgi and gndA. Third, metabolic evolution was carried out to generate strain NK-A9E, which showed a comparable γ-PGA titre with strain NK-A7. Finally, the srf and itu operons were deleted respectively, from the original strains NK-A7 and NK-A9E. The resulting strain NK-A11 exhibited the highest γ-PGA titre (7.53 g l-1 ), with a 2.05-fold improvement compared with LL3. The results demonstrated that the approaches described here efficiently enhanced γ-PGA production in B. amyloliquefaciens fermentation.

Genome-Scale Metabolic Network Models of Bacillus Species Suggest that Model Improvement is Necessary for Biotechnological Applications

Background

A genome-scale metabolic network model (GEM) is a mathematical representation of an organism’s metabolism. Today, GEMs are popular tools for computationally simulating the biotechnological processes and for predicting biochemical properties of (engineered) strains.

Objectives

In the present study, we have evaluated the predictive power of two GEMs, namely iBsu1103 (for Bacillus subtilis 168) and iMZ1055 (for Bacillus megaterium WSH002).

Materials and Methods

For comparing the predictive power of Bacillus subtilis and Bacillus megaterium GEMs, experimental data were obtained from previous wet-lab studies included in PubMed. By using these data, we set the environmental, stoichiometric and thermodynamic constraints on the models, and FBA is performed to predict the biomass production rate, and the values of other fluxes. For simulating experimental conditions in this study, COBRA toolbox was used.

Results

By using the wealth of data in the literature, we evaluated the accuracy of in silico simulations of these GEMs. Our results suggest that there are some errors in these two models which make them unreliable for predicting the biochemical capabilities of these species. The inconsistencies between experimental and computational data are even greater where B. subtilis and B. megaterium do not have similar phenotypes.

Conclusions

Our analysis suggests that literature-based improvement of genome-scale metabolic network models of the two Bacillus species is essential if these models are to be successfully applied in biotechnology and metabolic engineering.

An unconventional uptake rate objective function approach enhances applicability of genome-scale models for mammalian cells

Constraint-based modeling has been applied to analyze metabolism of numerous organisms via flux balance analysis and genome-scale metabolic models, including mammalian cells such as the Chinese hamster ovary (CHO) cells-the principal cell factory platform for therapeutic protein production. Unfortunately, the application of genome-scale model methodologies using the conventional biomass objective function is challenged by the presence of overly-restrictive constraints, including essential amino acid exchange fluxes that can lead to improper predictions of growth rates and intracellular flux distributions. In this study, these constraints are found to be reliably predicted by an "essential nutrient minimization" approach. After modifying these constraints with the predicted minimal uptake values, a series of unconventional objective functions are applied to minimize each individual non-essential nutrient uptake rate, revealing useful insights about metabolic exchange rates and flows across different cell lines and culture conditions. This unconventional uptake-rate objective functions (UOFs) approach is able to distinguish metabolic differences between three distinct CHO cell lines (CHO-K1, -DG44, and -S) not directly observed using the conventional biomass growth maximization solutions. Further, a comparison of model predictions with experimental data from literature correctly correlates with the specific CHO-DG44-derived cell line used experimentally, and the corresponding dual prices provide fruitful information concerning coupling relationships between nutrients. The UOFs approach is likely to be particularly suited for mammalian cells and other complex organisms which contain multiple distinct essential nutrient inputs, and may offer enhanced applicability for characterizing cell metabolism and physiology as well as media optimization and biomanufacturing control.

Genome-driven cell engineering review: in vivo and in silico metabolic and genome engineering

Producing 'designer cells' with specific functions is potentially feasible in the near future. Recent developments, including whole-cell models, genome design algorithms and gene editing tools, have advanced the possibility of combining biological research and mathematical modelling to further understand and better design cellular processes. In this review, we will explore computational and experimental approaches used for metabolic and genome design. We will highlight the relevance of modelling in this process, and challenges associated with the generation of quantitative predictions about cell behaviour as a whole: although many cellular processes are well understood at the subsystem level, it has proved a hugely complex task to integrate separate components together to model and study an entire cell. We explore these developments, highlighting where computational design algorithms compensate for missing cellular information and underlining where computational models can complement and reduce lab experimentation. We will examine issues and illuminate the next steps for genome engineering.

Pseudomonas putida KT2440 metabolism undergoes sequential modifications during exponential growth in a complete medium as compounds are gradually consumed

Pseudomonas putida is a soil bacterium with a versatile and robust metabolism. When confronted with mixtures of carbon sources, it prioritizes the utilization of the preferred compounds, optimizing metabolism and growth. This response is particularly strong when growing in a complex medium such as LB. This work examines the changes occurring in P. putida KT2440 metabolic fluxes, while it grows exponentially in LB medium and sequentially consumes the compounds available. Integrating the uptake rates for each compound at three different moments during the exponential growth with the changes observed in the proteome, and with the metabolic fluxes predicted by the iJN1411 metabolic model for this strain, allowed the metabolic rearrangements that occurred to be determined. The results indicate that the bacterium changes significantly the configuration of its metabolism during the early, mid and late exponential phases of growth. Sugars served as an energy source during the early phase and later as energy and carbon source. The configuration of the tricarboxylic acids cycle varied during growth, providing no energy in the early phase, and turning to a reductive mode in the mid phase and to an oxidative mode later on. This work highlights the dynamism and flexibility of P. putida metabolism.

iDS372, a Phenotypically Reconciled Model for the Metabolism of Streptococcus pneumoniae Strain R6

A high-quality GSM model for Streptococcus pneumoniae R6 model strain (iDS372), comprising 372 genes and 529 reactions, was developed. The construction of this model involved performing a genome-wide reannotation to identify the metabolic capacity of the bacterium. A reaction representing the abstraction of the biomass composition was reconciled from several studies reported in the literature and previous models, and included in the model. The final model comprises two compartments and manifold automatically generated gene rules. The validation was performed with experimental data from recent studies, regarding the usability of carbon sources, the effect of the presence of oxygen, and the requirement of amino acids for growth. This model can be used to better understand the metabolism of this major pathogen, provide clues regarding new drug targets, and eventually design strategies for fighting infections by these bacteria.

Identification of critical connectors in the directed reaction-centric graphs of microbial metabolic networks

BACKGROUND

Detection of central nodes in asymmetrically directed biological networks depends on centrality metrics quantifying individual nodes' importance in a network. In topological analyses on metabolic networks, various centrality metrics have been mostly applied to metabolite-centric graphs. However, centrality metrics including those not depending on high connections are largely unexplored for directed reaction-centric graphs.

RESULTS

We applied directed versions of centrality metrics to directed reaction-centric graphs of microbial metabolic networks. To investigate the local role of a node, we developed a novel metric, cascade number, considering how many nodes are closed off from information flow when a particular node is removed. High modularity and scale-freeness were found in the directed reaction-centric graphs and betweenness centrality tended to belong to densely connected modules. Cascade number and bridging centrality identified cascade subnetworks controlling local information flow and irreplaceable bridging nodes between functional modules, respectively. Reactions highly ranked with bridging centrality and cascade number tended to be essential, compared to reactions that other central metrics detected.

CONCLUSIONS

We demonstrate that cascade number and bridging centrality are useful to identify key reactions controlling local information flow in directed reaction-centric graphs of microbial metabolic networks. Knowledge about the local flow connectivity and connections between local modules will contribute to understand how metabolic pathways are assembled.

Current status and applications of genome-scale metabolic models

Genome-scale metabolic models (GEMs) computationally describe gene-protein-reaction associations for entire metabolic genes in an organism, and can be simulated to predict metabolic fluxes for various systems-level metabolic studies. Since the first GEM for Haemophilus influenzae was reported in 1999, advances have been made to develop and simulate GEMs for an increasing number of organisms across bacteria, archaea, and eukarya. Here, we review current reconstructed GEMs and discuss their applications, including strain development for chemicals and materials production, drug targeting in pathogens, prediction of enzyme functions, pan-reactome analysis, modeling interactions among multiple cells or organisms, and understanding human diseases.

Metabolic characterization of the chitinolytic bacterium Serratia marcescens using a genome-scale metabolic model

Background

Serratia marcescens is a chitinolytic bacterium that can potentially be used for consolidated bioprocessing to convert chitin to value-added chemicals. Currently, S. marcescens is poorly characterized and studies on intracellular metabolic and regulatory mechanisms would expedite development of bioprocessing applications.

Results

In this study, our goal was to characterize the metabolic profile of S. marcescens to provide insight for metabolic engineering applications and fundamental biological studies. Hereby, we constructed a constraint-based genome-scale metabolic model (iSR929) including 929 genes, 1185 reactions and 1164 metabolites based on genomic annotation of S. marcescens Db11. The model was tested by comparing model predictions with experimental data and analyzed to identify essential aspects of the metabolic network (e.g. 138 essential genes predicted). The model iSR929 was refined by integrating RNAseq data of S. marcescens growth on three different carbon sources (glucose, N-acetylglucosamine, and glycerol). Significant differences in TCA cycle utilization were found for growth on the different carbon substrates, For example, for growth on N-acetylglucosamine, S. marcescens exhibits high pentose phosphate pathway activity and nucleotide synthesis but low activity of the TCA cycle.

Conclusions

Our results show that S. marcescens model iSR929 can provide reasonable predictions and can be constrained to fit with experimental values. Thus, our model may be used to guide strain designs for metabolic engineering to produce chemicals such as 2,3-butanediol, N-acetylneuraminic acid, and n-butanol using S. marcescens.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2826-1) contains supplementary material, which is available to authorized users.

Genome-Scale Metabolic Model of a Microbial Cell Factory (Brevibacillus thermoruber 423) with Multi-Industry Potentials for Exopolysaccharide Production

Brevibacillus thermoruber 423 is a thermophilic bacterium capable of producing high levels of exopolysaccharide (EPS) that has broad applications in nutrition, feed, cosmetics, pharmaceutical, and chemical industries, not to mention in health and bionanotechnology sectors. EPS is a natural, nontoxic, and biodegradable polymer of sugar residues and plays pivotal roles in cell-to-cell interactions, adhesion, biofilm formation, and protection of cell against environmental extremes. This bacterium is a thermophilic EPS producer while exceeding other thermophilic producers by virtue of high level of polymer synthesis. Recently, B. thermoruber 423 was noted for relevance to multiple industry sectors because of its capacity to use xylose, and produce EPS, isoprenoids, ethanol/butanol, lipases, proteases, cellulase, and glucoamylase enzymes as well as its resistance to arsenic. A key step in understanding EPS production with a systems-based approach is the knowledge of microbial genome sequence. To speed biotechnology and industrial applications, this study reports on a genome-scale metabolic model (GSMM) of B. thermoruber 423, constructed using the recently available high-quality genome sequence that we have subsequently validated using physiological data on batch growth and EPS production on seven different carbon sources. The model developed contains 1454 reactions (of which 1127 are assigned an enzyme commission number) and 1410 metabolites from 925 genes. This GSMM offers the promise to enable and accelerate further systems biology and industrial scale studies, not to mention the ability to calculate metabolic flux distribution in large networks and multiomic data integration.

Bacillus velezensis: phylogeny, useful applications, and avenues for exploitation

Some members of the Bacillus velezensis (Bv) group (e.g., Bv FZB42T and AS3.43) were previously assigned grouping with B. subtilis and B. amyloliquefaciens, based on the fact that they shared a 99% DNA-DNA percentage phylogenetic similarity. However, hinging on current assessments of the pan-genomic reassignments, the differing phylogenomic characteristics of Bv from B. subtilis and B. amyloliquefaciens are now better understood. Within this re-grouping/reassignment, the various strains within the Bv share a close phylogenomic resemblance, and a number of these strains have received a lot of attention in recent years, due to their genomic robustness, and the growing evidence for their possible utilization in the agricultural industry for managing plant diseases. Only a few applications for their use medicinally/pharmaceutically, environmentally, and in the food industry have been reported, and this may be due to the fact that the majority of those strains investigated are those typically occurring in soil. Although the intracellular unique biomolecules of Bv strains have been revealed via in silico genome modeling and investigated using transcriptomics and proteomics, a further inquisition into the Bv metabolome using newer technologies such as metabolomics could elucidate additional applications of this economically relevant Bacillus species, beyond that of primarily the agricultural sector.

Costless metabolic secretions as drivers of interspecies interactions in microbial ecosystems

Metabolic exchange mediates interactions among microbes, helping explain diversity in microbial communities. As these interactions often involve a fitness cost, it is unclear how stable cooperation can emerge. Here we use genome-scale metabolic models to investigate whether the release of “costless” metabolites (i.e. those that cause no fitness cost to the producer), can be a prominent driver of intermicrobial interactions. By performing over 2 million pairwise growth simulations of 24 species in a combinatorial assortment of environments, we identify a large space of metabolites that can be secreted without cost, thus generating ample cross-feeding opportunities. In addition to providing an atlas of putative interactions, we show that anoxic conditions can promote mutualisms by providing more opportunities for exchange of costless metabolites, resulting in an overrepresentation of stable ecological network motifs. These results may help identify interaction patterns in natural communities and inform the design of synthetic microbial consortia.

In considering cross-feeding among microbes within communities, it is typically assumed that metabolic secretions are costly to produce. However, Pacheco et al. use metabolic models to show that ‘costless’ secretions could be common in some environments and important for structuring interactions among microbes.

Integration of enzymatic data in Bacillus subtilis genome-scale metabolic model improves phenotype predictions and enables in silico design of poly-γ-glutamic acid production strains

Background

Genome-scale metabolic models (GEMs) allow predicting metabolic phenotypes from limited data on uptake and secretion fluxes by defining the space of all the feasible solutions and excluding physio-chemically and biologically unfeasible behaviors. The integration of additional biological information in genome-scale models, e.g., transcriptomic or proteomic profiles, has the potential to improve phenotype prediction accuracy. This is particularly important for metabolic engineering applications where more accurate model predictions can translate to more reliable model-based strain design.

Results

Here we present a GEM with Enzymatic Constraints using Kinetic and Omics data (GECKO) model of Bacillus subtilis, which uses publicly available proteomic data and enzyme kinetic parameters for central carbon (CC) metabolic reactions to constrain the flux solution space. This model allows more accurate prediction of the flux distribution and growth rate of wild-type and single-gene/operon deletion strains compared to a standard genome-scale metabolic model. The flux prediction error decreased by 43% and 36% for wild-type and mutants respectively. The model additionally increased the number of correctly predicted essential genes in CC pathways by 2.5-fold and significantly decreased flux variability in more than 80% of the reactions with variable flux. Finally, the model was used to find new gene deletion targets to optimize the flux toward the biosynthesis of poly-γ-glutamic acid (γ-PGA) polymer in engineered B. subtilis. We implemented the single-reaction deletion targets identified by the model experimentally and showed that the new strains have a twofold higher γ-PGA concentration and production rate compared to the ancestral strain.

Conclusions

This work confirms that integration of enzyme constraints is a powerful tool to improve existing genome-scale models, and demonstrates the successful use of enzyme-constrained models in B. subtilis metabolic engineering. We expect that the new model can be used to guide future metabolic engineering efforts in the important industrial production host B. subtilis.

Chasing bacterial chassis for metabolic engineering: a perspective review from classical to non-traditional microorganisms

The last few years have witnessed an unprecedented increase in the number of novel bacterial species that hold potential to be used for metabolic engineering. Historically, however, only a handful of bacteria have attained the acceptance and widespread use that are needed to fulfil the needs of industrial bioproduction - and only for the synthesis of very few, structurally simple compounds. One of the reasons for this unfortunate circumstance has been the dearth of tools for targeted genome engineering of bacterial chassis, and, nowadays, synthetic biology is significantly helping to bridge such knowledge gap. Against this background, in this review, we discuss the state of the art in the rational design and construction of robust bacterial chassis for metabolic engineering, presenting key examples of bacterial species that have secured a place in industrial bioproduction. The emergence of novel bacterial chassis is also considered at the light of the unique properties of their physiology and metabolism, and the practical applications in which they are expected to outperform other microbial platforms. Emerging opportunities, essential strategies to enable successful development of industrial phenotypes, and major challenges in the field of bacterial chassis development are also discussed, outlining the solutions that contemporary synthetic biology-guided metabolic engineering offers to tackle these issues.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

IN SILICO GENOME-SCALE RECONSTRUCTION AND ANALYSIS OF THE SHEWANELLA LOIHICA PV-4 METABOLIC NETWORK

Microbial Fuel Cells (MFCs) are devices that generate electricity directly from organic compounds with microbes (electricigens) serving as anodic catalysts. As a novel environment-friendly energy source, MFCs have extensive practical value. Since the biological features and metabolic mechanism of electricigens have a great effect on the electricity production of MFCs, it is a big deal to screen strains with high electricity productivity for improving the power output of MFC. Reconstructions and simulations of metabolic networks are of significant help in studying the metabolism of microorganisms so as to guide gene engineering and metabolic engineering to improve their power-generating efficiency. Herein, we reconstructed a genome-scale constraint-based metabolic network model of Shewanella loihica PV-4, an important electricigen, based on its genomic functional annotations, reaction databases and published metabolic network models of seven microorganisms. The resulting network model iGX790 consists of 902 reactions (including 71 exchange reactions), 798 metabolites and 790 genes, covering the main pathways such as carbon metabolism, energy metabolism, amino acid metabolism, nucleic acid metabolism and lipid metabolism. Using the model, we simulated the growth rate, the maximal synthetic rate of ATP, the flux variability analysis of metabolic network, gene deletion and so on to examine the metabolism of S. loihica PV-4.

Fast automated reconstruction of genome-scale metabolic models for microbial species and communities

Genome-scale metabolic models are instrumental in uncovering operating principles of cellular metabolism, for model-guided re-engineering, and unraveling cross-feeding in microbial communities. Yet, the application of genome-scale models, especially to microbial communities, is lagging behind the availability of sequenced genomes. This is largely due to the time-consuming steps of manual curation required to obtain good quality models. Here, we present an automated tool, CarveMe, for reconstruction of species and community level metabolic models. We introduce the concept of a universal model, which is manually curated and simulation ready. Starting with this universal model and annotated genome sequences, CarveMe uses a top-down approach to build single-species and community models in a fast and scalable manner. We show that CarveMe models perform closely to manually curated models in reproducing experimental phenotypes (substrate utilization and gene essentiality). Additionally, we build a collection of 74 models for human gut bacteria and test their ability to reproduce growth on a set of experimentally defined media. Finally, we create a database of 5587 bacterial models and demonstrate its potential for fast generation of microbial community models. Overall, CarveMe provides an open-source and user-friendly tool towards broadening the use of metabolic modeling in studying microbial species and communities.

Sequence Specific Modeling of E. coli Cell-Free Protein Synthesis

Cell-free protein synthesis (CFPS) is a widely used research tool in systems and synthetic biology. However, if CFPS is to become a mainstream technology for applications such as point of care manufacturing, we must understand the performance limits and costs of these systems. Toward this question, we used sequence specific constraint based modeling to evaluate the performance of E. coli cell-free protein synthesis. A core E. coli metabolic network, describing glycolysis, the pentose phosphate pathway, energy metabolism, amino acid biosynthesis, and degradation was augmented with sequence specific descriptions of transcription and translation and effective models of promoter function. Model parameters were largely taken from literature; thus the constraint based approach coupled the transcription and translation of the protein product, and the regulation of gene expression, with the availability of metabolic resources using only a limited number of adjustable model parameters. We tested this approach by simulating the expression of two model proteins: chloramphenicol acetyltransferase and dual emission green fluorescent protein, for which we have data sets; we then expanded the simulations to a range of additional proteins. Protein expression simulations were consistent with measurements for a variety of cases. The constraint based simulations confirmed that oxidative phosphorylation was active in the CAT cell-free extract, as without it there was no feasible solution within the experimental constraints of the system. We then compared the metabolism of theoretically optimal and experimentally constrained CFPS reactions, and developed parameter free correlations which could be used to estimate productivity as a function of carbon number and promoter type. Lastly, global sensitivity analysis identified the key metabolic processes that controlled CFPS productivity and energy efficiency. In summary, sequence specific constraint based modeling of CFPS offered a novel means to a priori estimate the performance of a cell-free system, using only a limited number of adjustable parameters. While we modeled the production of a single protein in this study, the approach could easily be extended to multiprotein synthetic circuits, RNA circuits, or the cell-free production of small molecule products.

An enormous potential for niche construction through bacterial cross-feeding in a homogeneous environment

Microorganisms modify their environment by excreting by-products of metabolism, which can create new ecological niches that can help microbial populations diversify. A striking example comes from experimental evolution of genetically identical Escherichia coli populations that are grown in a homogeneous environment with the single carbon source glucose. In such experiments, stable communities of genetically diverse cross-feeding E. coli cells readily emerge. Some cells that consume the primary carbon source glucose excrete a secondary carbon source, such as acetate, that sustains other community members. Few such cross-feeding polymorphisms are known experimentally, because they are difficult to screen for. We studied the potential of bacterial metabolism to create new ecological niches based on cross-feeding. To do so, we used genome scale models of the metabolism of E. coli and metabolisms of similar complexity, to identify unique pairs of primary and secondary carbon sources in these metabolisms. We then combined dynamic flux balance analysis with analytical calculations to identify which pair of carbon sources can sustain a polymorphic cross-feeding community. We identified almost 10,000 such pairs of carbon sources, each of them corresponding to a unique ecological niche. Bacterial metabolism shows an immense potential for the construction of new ecological niches through cross feeding.

Biodiversity can emerge in a completely homogeneous environment from populations with initially genetically identical individuals. This striking observation comes from experimental evolution of bacteria, which create new ecological niches when they excrete nutrient-rich waste products that can sustain the life of other bacteria. It is difficult to estimate the potential of any one organism for such metabolic niche construction experimentally, because it is challenging to screen for novel metabolic abilities on a large scale. We therefore used experimentally validated models of bacterial metabolism to predict how many novel niches organisms like Escherichia coli can construct, if a novel niche must be able to sustain a stable community of microbes that differ in the nutrients they consume. We identify thousands of such niches. They differ in their primary carbon source and a secondary carbon source that is excreted by some microbes and used by others. Because we restricted ourselves to chemically simple environments, we may even have underestimated the enormous potential of microbes for niche construction.

Genome-scale metabolic reconstruction and analysis for Clostridium kluyveri

Clostridium kluyveri is an anaerobic microorganism that is well-known for producing butyrate and hexanoate using ethanol and acetate. It is also an important bacterium in the production of Chinese strong flavour baijiu (SFB). To obtain a comprehensive understanding of its metabolism, a curated genome-scale metabolic model (GSMM) of C. kluyveri, including 708 genes, 994 reactions, and 804 metabolites, was constructed and named iCKL708. This model was used to simulate the growth of C. kluyveri on different carbon substrates and the results agreed well with the experimental data. The butyrate, pentanoate, and hexanoate biosynthesis pathways were also elucidated. Flux balance analysis indicated that the ratio of ethanol to acetate, as well as the uptake rate of carbon dioxide, affected hexanoate production. The GSMM iCKL708 described here provides a platform to further our understanding and exploration of the metabolic potential of C. kluyveri.

Increased flux through the TCA cycle enhances bacitracin production by Bacillus licheniformis DW2

The dodecapeptide antibiotic bacitracin, produced by several strains of Bacillus licheniformis and Bacillus subtilis, is widely used as an antibacterial animal feed additive. Several genetic strategies were explored to enhance its production. The availability of building block amino acids for bacitracin production was found to play an important role in its synthesis. In this study, the TCA cycle in the industrial strain B. licheniformis DW2 was strengthened by overexpression of the key enzymes citrate synthase and isocitrate dehydrogenase (ICDH). As the central metabolic pathway, the TCA cycle is a major source for energy supply and intermediates for anabolism. By enhancing flux through the TCA cycle, more energy and precursors were generated for amino acid biosynthesis and uptake, resulting in enlarged intracellular pool of bacitracin-containing amino acids for bacitracin production. This study unveiled the metabolic responses of the increased TCA cycle flux in B. licheniformis and provided a novel strategy for enhancing bacitracin production.

Using genome-scale metabolic models to compare serovars of the foodborne pathogen Listeria monocytogenes

Listeria monocytogenes is a microorganism of great concern for the food industry and the cause of human foodborne disease. Therefore, novel methods of control are needed, and systems biology is one such approach to identify them. Using a combination of computational techniques and laboratory methods, genome-scale metabolic models (GEMs) can be created, validated, and used to simulate growth environments and discern metabolic capabilities of microbes of interest, including L. monocytogenes. The objective of the work presented here was to generate GEMs for six different strains of L. monocytogenes, and to both qualitatively and quantitatively validate these GEMs with experimental data to examine the diversity of metabolic capabilities of numerous strains from the three different serovar groups most associated with foodborne outbreaks and human disease. Following qualitative validation, 57 of the 95 carbon sources tested experimentally were present in the GEMs, and; therefore, these were the compounds from which comparisons could be drawn. Of these 57 compounds, agreement between in silico predictions and in vitro results for carbon source utilization ranged from 80.7% to 91.2% between strains. Nutrient utilization agreement between in silico predictions and in vitro results were also conducted for numerous nitrogen, phosphorous, and sulfur sources. Additionally, quantitative validation showed that the L. monocytogenes GEMs were able to generate in silico predictions for growth rate and growth yield that were strongly and significantly (p < 0.0013 and p < 0.0015, respectively) correlated with experimental results. These findings are significant because they show that these GEMs for L. monocytogenes are comparable to published GEMs of other organisms for agreement between in silico predictions and in vitro results. Therefore, as with the other GEMs, namely those for Escherichia coli, Staphylococcus aureus, Vibrio vulnificus, and Salmonella spp., they can be used to determine new methods of growth control and disease treatment.

On an algorithmic definition for the components of the minimal cell

Living cells are highly complex systems comprising a multitude of elements that are engaged in the many convoluted processes observed during the cell cycle. However, not all elements and processes are essential for cell survival and reproduction under steady-state environmental conditions. To distinguish between essential from expendable cell components and thus define the ‘minimal cell’ and the corresponding ‘minimal genome’, we postulate that the synthesis of all cell elements can be represented as a finite set of binary operators, and within this framework we show that cell elements that depend on their previous existence to be synthesized are those that are essential for cell survival. An algorithm to distinguish essential cell elements is presented and demonstrated within an interactome. Data and functions implementing the algorithm are given as supporting information. We expect that this algorithmic approach will lead to the determination of the complete interactome of the minimal cell, which could then be experimentally validated. The assumptions behind this hypothesis as well as its consequences for experimental and theoretical biology are discussed.