Search

Abstract The bacterium Bacillus subtilis undergoes asymmetric cell division during sporulation, producing a mother cell and a smaller forespore connected by the SpoIIQ-SpoIIIA (or Q-A) channel. The two cells differentiate metabolically, and the forespore becomes dependent on the mother cell for essential building blocks. Here, we investigate the metabolic interactions between mother cell and forespore using genome-scale metabolic and expression models as well as experiments. Our results indicate that nucleotides are synthesized in the mother cell and transported in the form of nucleoside di- or tri-phosphates to the forespore via the Q-A channel. However, if the Q-A channel is inactivated later in sporulation, then glycolytic enzymes can form an ATP and NADH shuttle, providing the forespore with energy and reducing power. Our integrated in silico and in vivo approach sheds light into the intricate metabolic interactions underlying cell differentiation in B. subtilis, and provides a foundation for future studies of metabolic differentiation.

Abstract Bacillus subtilis is a well-characterized microorganism and a model for the study of Gram-positive bacteria. The bacterium can produce proteins at high densities and yields, which has made it valuable for industrial bioproduction. Like other cell factories, metabolic modeling of B. subtilis has discovered ways to optimize its metabolism toward various applications. The first genome-scale metabolic model (M-model) of B. subtilis was published more than a decade ago and has been applied extensively to understand metabolism, to predict growth phenotypes, and served as a template to reconstruct models for other Gram-positive bacteria. However, M-models are ill-suited to simulate the production and secretion of proteins as well as their proteomic response to stress. Thus, a new generation of metabolic models, known as metabolism and gene expression models (ME-models), has been initiated. Here, we describe the reconstruction and validation of a ME model of B. subtilis, iJT964-ME. This model achieved higher performance scores on the prediction of gene essentiality as compared to the M-model. We successfully validated the model by integrating physiological and omics data associated with gene expression responses to ethanol and salt stress. The model further identified the mechanism by which tryptophan synthesis is upregulated under ethanol stress. Further, we employed iJT964-ME to predict amylase production rates under two different growth conditions. We analyzed these flux distributions and identified key metabolic pathways that permitted the increase in amylase production. Models like iJT964-ME enable the study of proteomic response to stress and the illustrate the potential for optimizing protein production in bacteria.

Nitrogen fixation is an important metabolic process carried out by microorganisms, which converts molecular nitrogen into inorganic nitrogenous compounds such as ammonia (NH3). These nitrogenous compounds are crucial for biogeochemical cycles and for the synthesis of essential biomolecules, i.e. nucleic acids, amino acids and proteins. Azotobacter vinelandii is a bacterial non-photosynthetic model organism to study aerobic nitrogen fixation (diazotrophy) and hydrogen production. Moreover, the diazotroph can produce biopolymers like alginate and polyhydroxybutyrate (PHB) that have important industrial applications. However, many metabolic processes such as partitioning of carbon and nitrogen metabolism in A. vinelandii remain unknown to date.Genome-scale metabolic models (M-models) represent reliable tools to unravel and optimize metabolic functions at genome-scale. M-models are mathematical representations that contain information about genes, reactions, metabolites and their associations. M-models can simulate optimal reaction fluxes under a wide variety of conditions using experimentally determined constraints. Here we report on the development of a M-model of the wild type bacterium A. vinelandii DJ (iDT1278) which consists of 2,003 metabolites, 2,469 reactions, and 1,278 genes. We validated the model using high-throughput phenotypic and physiological data, testing 180 carbon sources and 95 nitrogen sources. iDT1278 was able to achieve an accuracy of 89% and 91% for growth with carbon sources and nitrogen source, respectively. This comprehensive M-model will help to comprehend metabolic processes associated with nitrogen fixation, ammonium assimilation, and production of organic nitrogen in an environmentally important microorganism.

Abstract Production of biofuels and bioenergy precursors by phototrophic microorganisms, such as microalgae and cyanobacteria, is a promising alternative to conventional fuels obtained from non-renewable resources. Several species of microalgae have been investigated as potential candidates for the production of biofuels, for the most part due to their exceptional metabolic capability to accumulate large quantities of lipids. Constraint-based modeling, a systems biology approach that accurately predicts the metabolic phenotype of phototrophs, has been deployed to identify suitable culture conditions as well as to explore genetic enhancement strategies for bioproduction. Core metabolic models were employed to gain insight into the central carbon metabolism in photosynthetic microorganisms. More recently, comprehensive genome-scale models, including organelle-specific information at high resolution, have been developed to gain new insight into the metabolism of phototrophic cell factories. Here, we review the current state of the art of constraint-based modeling and computational method development and discuss how advanced models led to increased prediction accuracy and thus improved lipid production in microalgae.

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

The denitrifying bacterium Thauera sp. MZ1T, a common member of microbial communities in wastewater treatment facilities, can produce different compounds from a range of carbon (C) and nitrogen (N) sources under aerobic and anaerobic conditions. In these different conditions, Thauera modifies its metabolism to produce different compounds that influence the microbial community. In particular, Thauera sp. MZ1T produces different exopolysaccharides with floc-forming properties, impacting the physical disposition of wastewater consortia and the efficiency of nutrient assimilation by the microbial community. Under N-limiting conditions, Thauera sp. MZ1T decreases its growth rate and accelerates the accumulation of polyhydroxyalkanoate-related (PHA) compounds including polyhydroxybutyrate (PHB), which plays a fundamental role as C and energy storage in this β-proteobacterium. However, the metabolic mechanisms employed by Thauera sp. MZ1T to assimilate and catabolize many of the different C and N sources under aerobic and anaerobic conditions remain unknown. Systems biology approaches such as genome-scale metabolic modeling have been successfully used to unveil complex metabolic mechanisms for various microorganisms. Here, we developed a comprehensive metabolic model (M-model) for Thauera sp. MZ1T (iThauera861), consisting of 1,744 metabolites, 2,384 reactions, and 861 genes. We validated the model experimentally using over 70 different C and N sources under both aerobic and anaerobic conditions. iThauera861 achieved a prediction accuracy of 95% for growth on various C and N sources and close to 85% for assimilation of aromatic compounds under denitrifying conditions. The M-model was subsequently deployed to determine the effects of substrates, oxygen presence, and the C:N ratio on the production of PHB and exopolysaccharides (EPS), showing the highest polymer yields are achieved with nucleotides and amino acids under aerobic conditions. This comprehensive M-model will help reveal the metabolic processes by which this ubiquitous species influences communities in wastewater treatment systems and natural environments.