Search

Present in many industrial effluents and as common degradation product of organic matter, phenol is a widespread compound which may cause serious environmental problems, due to its toxicity to animals and humans. Degradation of phenol from the environment by mesophilic bacteria has been studied extensively over the past decades, but only little is known about phenol biodegradation at high temperatures or low pH. In this work we studied phenol degradation in the thermoacidophilic archaeon Saccharolobus solfataricus P2 (basonym: Sulfolobus solfataricus) under extreme conditions (80°C, pH 3.5). We combined metabolomics and transcriptomics together with metabolic modeling to elucidate the organism's response to growth with phenol as sole carbon source. Although S. solfataricus is able to utilize phenol for biomass production, the carbon source induces profound stress reactions, including genome rearrangement as well as a strong intracellular accumulation of polyamines. Furthermore, computational modeling revealed a 40% higher oxygen demand for substrate oxidation, compared to growth on glucose. However, only 16.5% of oxygen is used for oxidation of phenol to catechol, resulting in a less efficient integration of carbon into the biomass. Finally, our data underlines the importance of the phenol meta -degradation pathway in S. solfataricus and enables us to predict enzyme candidates involved in the degradation processes downstream of 2-hydroxymucconic acid. [ABSTRACT FROM AUTHOR]

Metabolism represents the biological level that is most related to the phenotypes of the cell and, alterations or reprogramming of this can (i) affect the production of primary or secondary metabolites in microorganisms of biotechnological interest, (ii) favor or not the inhibition of growth in pathogenic organisms and (iii) developing metabolic disorders such as obesity or diabetes, among others. The study of metabolism, redesign, and redirection of metabolic fluxes has become an important area of research (also known as Metabolic Engineering), as it has allowed the development and design of improved biological processes, the identification of therapeutic targets, the design of therapeutic strategies to cure metabolic disorders and the identification of biomarkers in cancer, among others. Currently, the development of computational methodologies is allowing to study the cell metabolism under different environmental conditions and directing experiments with the model's predictions. The purpose of this review is to highlight the importance of metabolic flux analysis as a general methodology to study metabolic reprogramming in different organisms of biotechnological, medical, and therapeutic interest. This work condenses the theoretical bases and key concepts to understand the analysis of metabolic fluxes, which will be a fundamental input for those who are entering to the world of systems biology or related areas. [ABSTRACT FROM AUTHOR]

The gut microbiota plays a crucial role in health and is significantly modulated by human diets. In addition to Western diets which are rich in proteins, high-protein diets are used for specific populations or indications, mainly weight loss. In this study, we investigated the effect of protein supplementation on Bacteroides caccae , a Gram-negative gut symbiont. The supplementation with whey proteins led to a significant increase in growth rate, final biomass, and short-chain fatty acids production. A comprehensive genomic analysis revealed that B. caccae possesses a set of 156 proteases with putative intracellular and extracellular localization and allowed to identify amino acid transporters and metabolic pathways. We developed a fully curated genome-scale metabolic model of B. caccae that incorporated its proteolytic activity and simulated its growth and production of fermentation-related metabolites in response to the different growth media. We validated the model by comparing the predicted phenotype to experimental data. The model accurately predicted B. caccae 's growth and metabolite production ( R 2 = 0.92 for the training set and R 2 = 0.89 for the validation set). We found that accounting for both ATP consumption related to proteolysis, and whey protein accessibility is necessary for accurate predictions of metabolites production. These results provide insights into B. caccae 's adaptation to a high-protein diet and its ability to utilize proteins as a source of nutrition. The proposed model provides a useful tool for understanding the feeding mechanism of B. caccae in the gut microbiome.IMPORTANCEMicrobial proteolysis is understudied despite the availability of dietary proteins for the gut microbiota. Here, the proteolytic potential of the gut symbiont Bacteroides caccae was analyzed for the first time using pan-genomics. This sketches a well-equipped bacteria for protein breakdown, capable of producing 156 different proteases with a broad spectrum of cleavage targets. This functional potential was confirmed by the enhancement of growth and metabolic activities at high protein levels. Proteolysis was included in a B. caccae metabolic model which was fitted with the experiments and validated on external data. This model pinpoints the links between protein availability and short-chain fatty acids production, and the importance for B. caccae to gain access to glutamate and asparagine to promote growth. This integrated approach can be generalized to other symbionts and upscaled to complex microbiota to get insights into the ecological impact of proteins on the gut microbiota., Competing Interests: The authors declare no conflict of interest.

Legumes are unique among plants as they can obtain nitrogen through symbiosis with nitrogen-fixing rhizobia that form root nodules in the host plants. Therefore they are valuable crops for sustainable agriculture. Increasing nitrogen fixation efficiency is not only important for achieving better plant growth and yield, but it is also crucial for reducing the use of nitrogen fertilizer. Arbuscular mycorrhizal fungi (AMF) are another group of important beneficial microorganisms that form symbiotic relationships with legumes. AMF can promote host plant growth by providing mineral nutrients and improving the soil ecosystem. The trilateral legume-rhizobia-AMF symbiotic relationships also enhance plant development and tolerance against biotic and abiotic stresses. It is known that domestication and agricultural activities have led to the reduced genetic diversity of cultivated germplasms and higher sensitivity to nutrient deficiencies in crop plants, but how domestication has impacted the capability of legumes to establish beneficial associations with rhizospheric microbes (including rhizobia and fungi) is not well-studied. In this review, we will discuss the impacts of domestication and agricultural practices on the interactions between legumes and soil microbes, focusing on the effects on AMF and rhizobial symbioses and hence nutrient acquisition by host legumes. In addition, we will summarize the genes involved in legume-microbe interactions and studies that have contributed to a better understanding of legume symbiotic associations using metabolic modeling.

Salt tolerance in the γ-proteobacterium Halomonas elongata is linked to its ability to produce the compatible solute ectoine. The metabolism of ectoine production is of great interest since it can shed light on the biochemical basis of halotolerance as well as pave the way for the improvement of the biotechnological production of such compatible solute. Ectoine belongs to the biosynthetic family of aspartate-derived amino-acids. Aspartate is formed from oxaloacetate, thereby connecting ectoine production to the anaplerotic reactions that refill carbon into the tricarboxylic acid cycle (TCA cycle). This places a high demand on these reactions and creates the need to regulate them not only in response to growth but also in response to extracellular salt concentration. In this work, we combine modeling and experiments to analyze how these different needs shape the anaplerotic reactions in H. elongata. First, the stoichiometric and thermodynamic factors that condition the flux distributions are analyzed, then the optimal patterns of operation for oxaloacetate production are calculated. Finally, the phenotype of two deletion mutants lacking potentially relevant anaplerotic enzymes: phosphoenolpyruvate carboxylase (Ppc) and oxaloacetate decarboxylase (Oad) are experimentally characterized. The results show that the anaplerotic reactions in H. elongata are indeed subject to evolutionary pressures that differ from those faced by other gram-negative bacteria. Ectoine producing halophiles must meet a higher metabolic demand for oxaloacetate and the reliance of many marine bacteria on the Entner-Doudoroff pathway compromises the anaplerotic efficiency of Ppc, which is usually one of the main enzymes fulfilling this role. The anaplerotic flux in H. elongata is contributed not only by Ppc but also by Oad, an enzyme that has not yet been shown to play this role in vivo. Ppc is necessary for H. elongata to grow normally at low salt concentrations but it is not required to achieve near maximal growth rates as long as there is a steep sodium gradient. On the other hand, the lack of Oad presents serious difficulties to grow at high salt concentrations. This points to a shared role of these two enzymes in guaranteeing the supply of oxaloacetate for biosynthetic reactions.

Microbial communities are pervasive in the natural environment, associated with many hosts, and of increasing importance in biotechnological applications. The complexity of these microbial systems makes the underlying mechanisms driving their dynamics difficult to identify. While experimental meta-OMICS techniques are routinely applied to record the inventory and activity of microbiomes over time, it remains difficult to obtain quantitative predictions based on such data. Mechanistic, quantitative mathematical modeling approaches hold the promise to both provide predictive power and shed light on cause-effect relationships driving these dynamic systems. We introduce μbialSim (pronounced “microbial sim”), a dynamic Flux-Balance-Analysis-based (dFBA) numerical simulator which is able to predict the time course in terms of composition and activity of microbiomes containing 100s of species in batch or chemostat mode. Activity of individual species is simulated by using separate FBA models which have access to a common pool of compounds, allowing for metabolite exchange. A novel augmented forward Euler method ensures numerical accuracy by temporarily reducing the time step size when compound concentrations decrease rapidly due to high compound affinities and/or the presence of many consuming species. We present three exemplary applications of μbialSim: a batch culture of a hydrogenotrophic archaeon, a syntrophic methanogenic biculture, and a 773-species human gut microbiome which exhibits a complex and dynamic pattern of metabolite exchange. Focusing on metabolite exchange as the main interaction type, μbialSim allows for the mechanistic simulation of microbiomes at their natural complexity. Simulated trajectories can be used to contextualize experimental meta-OMICS data and to derive hypotheses on cause-effect relationships driving community dynamics based on scenario simulations. μbialSim is implemented in Matlab and relies on the COBRA Toolbox or CellNetAnalyzer for FBA calculations. The source code is available under the GNU General Public License v3.0 at https://git.ufz.de/UMBSysBio/microbialsim.

ABSTRACT 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. IMPORTANCE 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.

Mesenchymal stem cells are a promising source for externally grown tissue replacements and patient-specific immunomodulatory treatments. This promise has not yet been fulfilled in part due to production scaling issues and the need to maintain the correct phenotype after re-implantation. One aspect of extracorporeal growth that may be manipulated to optimize cell growth and differentiation is metabolism. The metabolism of MSCs changes during and in response to differentiation and immunomodulatory changes. MSC metabolism may be linked to functional differences but how this occurs and influences MSC function remains unclear. Understanding how MSC metabolism relates to cell function is however important as metabolite availability and environmental circumstances in the body may affect the success of implantation. Genome-scale constraint based metabolic modeling can be used as a tool to fill gaps in knowledge of MSC metabolism, acting as a framework to integrate and understand various data types (e.g., genomic, transcriptomic and metabolomic). These approaches have long been used to optimize the growth and productivity of bacterial production systems and are being increasingly used to provide insights into human health research. Production of tissue for implantation using MSCs requires both optimized production of cell mass and the understanding of the patient and phenotype specific metabolic situation. This review considers the current knowledge of MSC metabolism and how it may be optimized along with the current and future uses of genome scale constraint based metabolic modeling to further this aim.

ABSTRACT Correlation-based analysis of paired microbiome-metabolome data sets is becoming a widespread research approach, aiming to comprehensively identify microbial drivers of metabolic variation. To date, however, the limitations of this approach and other microbiome-metabolome analysis methods have not been comprehensively evaluated. To address this challenge, we have introduced a mathematical framework to quantify the contribution of each taxon to metabolite variation based on uptake and secretion fluxes. We additionally used a multispecies metabolic model to simulate simplified gut communities, generating idealized microbiome-metabolome data sets. We then compared observed taxon-metabolite correlations in these data sets to calculated ground truth taxonomic contribution values. We found that in simulations of both a representative simple 10-species community and complex human gut microbiota, correlation-based analysis poorly identified key contributors, with an extremely low predictive value despite the idealized setting. We further demonstrate that the predictive value of correlation analysis is strongly influenced by both metabolite and taxon properties, as well as by exogenous environmental variation. We finally discuss the practical implications of our findings for interpreting microbiome-metabolome studies. IMPORTANCE Identifying the key microbial taxa responsible for metabolic differences between microbiomes is an important step toward understanding and manipulating microbiome metabolism. To achieve this goal, researchers commonly conduct microbiome-metabolome association studies, comprehensively measuring both the composition of species and the concentration of metabolites across a set of microbial community samples and then testing for correlations between microbes and metabolites. Here, we evaluated the utility of this general approach by first developing a rigorous mathematical definition of the contribution of each microbial taxon to metabolite variation and then examining these contributions in simulated data sets of microbial community metabolism. We found that standard correlation-based analysis of our simulated microbiome-metabolome data sets can identify true contributions with very low predictive value and that its performance depends strongly on specific properties of both metabolites and microbes, as well as on those of the surrounding environment. Combined, our findings can guide future interpretation and validation of microbiome-metabolome studies.

Situations of extremely low substrate availability, resulting in slow growth, are common in natural environments. To mimic these conditions, Lactobacillus plantarum was grown in a carbon-limited retentostat with complete biomass retention. The physiology of extremely slow-growing L. plantarum--as studied by genome-scale modeling and transcriptomics--was fundamentally different from that of stationary-phase cells. Stress resistance mechanisms were not massively induced during transition to extremely slow growth. The energy-generating metabolism was remarkably stable and remained largely based on the conversion of glucose to lactate. The combination of metabolic and transcriptomic analyses revealed behaviors involved in interactions with the environment, more particularly with plants: production of plant hormones or precursors thereof, and preparedness for the utilization of plant-derived substrates. Accordingly, the production of compounds interfering with plant root development was demonstrated in slow-growing L. plantarum. Thus, conditions of slow growth and limited substrate availability seem to trigger a plant environment-like response, even in the absence of plant-derived material, suggesting that this might constitute an intrinsic behavior in L. plantarum.

Abstract Background Microbial bioengineering has the potential to become a key contributor to the future development of human society by providing sustainable, novel, and cost-effective production pipelines. However, the sustained productivity of genetically engineered strains is often a challenge, as spontaneous non-producing mutants tend to grow faster and take over the population. Novel strategies to prevent this issue of strain instability are urgently needed. Results In this study, we propose a novel strategy applicable to all microbial production systems for which a genome-scale metabolic model is available that aligns the production of native metabolites to the formation of biomass. Based on well-established constraint-based analysis techniques such as OptKnock and FVA, we developed an in silico pipeline—FRUITS—that specifically ‘Finds Reactions Usable in Tapping Side-products’. It analyses a metabolic network to identify compounds produced in anabolism that are suitable to be coupled to growth by deletion of their re-utilization pathway(s), and computes their respective biomass and product formation rates. When applied to Synechocystis sp. PCC6803, a model cyanobacterium explored for sustainable bioproduction, a total of nine target metabolites were identified. We tested our approach for one of these compounds, acetate, which is used in a wide range of industrial applications. The model-guided engineered strain shows an obligatory coupling between acetate production and photoautotrophic growth as predicted. Furthermore, the stability of acetate productivity in this strain was confirmed by performing prolonged turbidostat cultivations. Conclusions This work demonstrates a novel approach to stabilize the production of target compounds in cyanobacteria that culminated in the first report of a photoautotrophic growth-coupled cell factory. The method developed is generic and can easily be extended to any other modeled microbial production system.

Methyl ketones present a group of highly reduced platform chemicals industrially produced from petroleum-derived hydrocarbons. They find applications in the fragrance, flavor, pharmacological, and agrochemical industries, and are further discussed as biodiesel blends. In recent years, intense research has been carried out to achieve sustainable production of these molecules by re-arranging the fatty acid metabolism of various microbes. One challenge in the development of a highly productive microbe is the high demand for reducing power. Here, we engineered Pseudomonas taiwanensis VLB120 for methyl ketone production as this microbe has been shown to sustain exceptionally high NAD(P)H regeneration rates. The implementation of published strategies resulted in 2.1 g L aq −1 methyl ketones in fed-batch fermentation. We further increased the production by eliminating competing reactions suggested by metabolic analyses. These efforts resulted in the production of 9.8 g L aq −1 methyl ketones (corresponding to 69.3 g L org −1 in the in situ extraction phase) at 53% of the maximum theoretical yield. This represents a 4-fold improvement in product titer compared to the initial production strain and the highest titer of recombinantly produced methyl ketones reported to date. Accordingly, this study underlines the high potential of P. taiwanensis VLB120 to produce methyl ketones and emphasizes model-driven metabolic engineering to rationalize and accelerate strain optimization efforts. Image 1 • Metabolic engineering of Pseudomonas taiwanensis VLB120 enables high-level methyl ketone production. • Model-supported strain design improved the product titer by 4-fold. • Carbon-excess conditions favor carbon flux towards methyl ketone synthesis. • Highest reported methyl ketone titer and yield in a recombinant microbial host. [ABSTRACT FROM AUTHOR]

Legumes are unique among plants as they can obtain nitrogen through symbiosis with nitrogen-fixing rhizobia that form root nodules in the host plants. Therefore they are valuable crops for sustainable agriculture. Increasing nitrogen fixation efficiency is not only important for achieving better plant growth and yield, but it is also crucial for reducing the use of nitrogen fertilizer. Arbuscular mycorrhizal fungi (AMF) are another group of important beneficial microorganisms that form symbiotic relationships with legumes. AMF can promote host plant growth by providing mineral nutrients and improving the soil ecosystem. The trilateral legume-rhizobia-AMF symbiotic relationships also enhance plant development and tolerance against biotic and abiotic stresses. It is known that domestication and agricultural activities have led to the reduced genetic diversity of cultivated germplasms and higher sensitivity to nutrient deficiencies in crop plants, but how domestication has impacted the capability of legumes to establish beneficial associations with rhizospheric microbes (including rhizobia and fungi) is not well-studied. In this review, we will discuss the impacts of domestication and agricultural practices on the interactions between legumes and soil microbes, focusing on the effects on AMF and rhizobial symbioses and hence nutrient acquisition by host legumes. In addition, we will summarize the genes involved in legume-microbe interactions and studies that have contributed to a better understanding of legume symbiotic associations using metabolic modeling. [ABSTRACT FROM AUTHOR]

Salt tolerance in the γ-proteobacterium Halomonas elongata is linked to its ability to produce the compatible solute ectoine. The metabolism of ectoine production is of great interest since it can shed light on the biochemical basis of halotolerance as well as pave the way for the improvement of the biotechnological production of such compatible solute. Ectoine belongs to the biosynthetic family of aspartate-derived amino-acids. Aspartate is formed from oxaloacetate, thereby connecting ectoine production to the anaplerotic reactions that refill carbon into the tricarboxylic acid cycle (TCA cycle). This places a high demand on these reactions and creates the need to regulate them not only in response to growth but also in response to extracellular salt concentration. In this work, we combine modeling and experiments to analyze how these different needs shape the anaplerotic reactions in H. elongata. First, the stoichiometric and thermodynamic factors that condition the flux distributions are analyzed, then the optimal patterns of operation for oxaloacetate production are calculated. Finally, the phenotype of two deletion mutants lacking potentially relevant anaplerotic enzymes: phosphoenolpyruvate carboxylase (Ppc) and oxaloacetate decarboxylase (Oad) are experimentally characterized. The results show that the anaplerotic reactions in H. elongata are indeed subject to evolutionary pressures that differ from those faced by other gram-negative bacteria. Ectoine producing halophiles must meet a higher metabolic demand for oxaloacetate and the reliance of many marine bacteria on the Entner-Doudoroff pathway compromises the anaplerotic efficiency of Ppc, which is usually one of the main enzymes fulfilling this role. The anaplerotic flux in H. elongata is contributed not only by Ppc but also by Oad, an enzyme that has not yet been shown to play this role in vivo. Ppc is necessary for H. elongata to grow normally at low salt concentrations but it is not required to achieve near maximal growth rates as long as there is a steep sodium gradient. On the other hand, the lack of Oad presents serious difficulties to grow at high salt concentrations. This points to a shared role of these two enzymes in guaranteeing the supply of oxaloacetate for biosynthetic reactions. [ABSTRACT FROM AUTHOR]

Cell line-specific, genome-scale metabolic models enable rigorous and systematic in silico investigation of cellular metabolism. Such models have recently become available for Chinese hamster ovary (CHO) cells. However, a key ingredient, namely an experimentally validated biomass function that summarizes the cellular composition, was so far missing. Here, we close this gap by providing extensive experimental data on the biomass composition of 13 parental and producer CHO cell lines under various conditions. We report total protein, lipid, DNA, RNA and carbohydrate content, cell dry mass, and detailed protein and lipid composition. Furthermore, we present meticulous data on exchange rates between cells and environment and provide detailed experimental protocols on how to determine all of the above. The biomass composition is converted into cell line- and condition-specific biomass functions for use in cell line-specific, genome-scale metabolic models of CHO. Finally, flux balance analysis (FBA) is used to demonstrate consistency between in silico predictions and experimental analysis. Our study reveals a strong variability of the total protein content and cell dry mass across cell lines. However, the relative amino acid composition is independent of the cell line and condition and thus needs not be explicitly measured for each new cell line. In contrast, the lipid composition is strongly influenced by the growth media and thus will have to be determined in each case. These cell line-specific variations in biomass composition have a small impact on growth rate predictions with FBA, as inaccuracies in the predictions are rather dominated by inaccuracies in the exchange rate spectra. Cell-specific biomass variations only become important if the experimental errors in the exchange rate spectra drop below twenty percent. Image 1 [ABSTRACT FROM AUTHOR]

Metabolic modeling and machine learning are key components in the emerging next generation of systems and synthetic biology tools, targeting the genotype-phenotype-environment relationship. Rather than being used in isolation, it is becoming clear that their value is maximized when they are combined. However, the potential of integrating these two frameworks for omic data augmentation and integration is largely unexplored. We propose, rigorously assess, and compare machine-learning-based data integration techniques, combining gene expression profiles with computationally generated metabolic flux data to predict yeast cell growth. To this end, we create strain-specific metabolic models for 1,143 Saccharomyces cerevisiae mutants and we test 27 machine-learning methods, incorporating state-of-the-art feature selection and multiview learning approaches. We propose a multiview neural network using fluxomic and transcriptomic data, showing that the former increases the predictive accuracy of the latter and reveals functional patterns that are not directly deducible from gene expression alone. We test the proposed neural network on a further 86 strains generated in a different experiment, therefore verifying its robustness to an additional independent dataset. Finally, we show that introducing mechanistic flux features improves the predictions also for knockout strains whose genes were not modeled in the metabolic reconstruction. Our results thus demonstrate that fusing experimental cues with in silico models, based on known biochemistry, can contribute with disjoint information toward biologically informed and interpretable machine learning. Overall, this study provides tools for understanding and manipulating complex phenotypes, increasing both the prediction accuracy and the extent of discernible mechanistic biological insights. [ABSTRACT FROM AUTHOR]

Comprehensive molecular‐level models of human metabolism have been generated on a cellular level. However, models of whole‐body metabolism have not been established as they require new methodological approaches to integrate molecular and physiological data. We developed a new metabolic network reconstruction approach that used organ‐specific information from literature and omics data to generate two sex‐specific whole‐body metabolic (WBM) reconstructions. These reconstructions capture the metabolism of 26 organs and six blood cell types. Each WBM reconstruction represents whole‐body organ‐resolved metabolism with over 80,000 biochemical reactions in an anatomically and physiologically consistent manner. We parameterized the WBM reconstructions with physiological, dietary, and metabolomic data. The resulting WBM models could recapitulate known inter‐organ metabolic cycles and energy use. We also illustrate that the WBM models can predict known biomarkers of inherited metabolic diseases in different biofluids. Predictions of basal metabolic rates, by WBM models personalized with physiological data, outperformed current phenomenological models. Finally, integrating microbiome data allowed the exploration of host–microbiome co‐metabolism. Overall, the WBM reconstructions, and their derived computational models, represent an important step toward virtual physiological humans. Synopsis: Sex‐specific, whole‐body human metabolic models were developed and constrained with physiological, dietary, and metabolomic data. They recapitulate known whole‐body metabolic functions and enable mechanistic exploration of host‐microbiome co‐metabolism. Sex‐specific whole‐body metabolic reconstructions represent the integrated function of 26 organs and six blood cell types.Stoichiometric reconstructions of metabolism can be constrained with whole‐body physiological and metabolomic data to generate personalized models.Whole‐body metabolic models recapitulate known inter‐organ metabolic cycles and energy use, successfully predict known biomarkers of inherited metabolic diseases, and explore potential host‐microbiome co‐metabolism. [ABSTRACT FROM AUTHOR]

Background: The systems-scale analysis of cellular metabolites, "metabolomics," provides data ideal for applications in metabolic engineering. However, many of the computational tools for strain design are built around Flux Balance Analysis (FBA), which makes assumptions that preclude direct integration of metabolomics data into the underlying models. Finding a way to retain the advantages of FBA's linear structure while relaxing some of its assumptions could allow us to account for metabolite levels and metabolite-dependent regulation in strain design tools built from FBA, improving the accuracy of predictions made by these tools. We designed, implemented, and characterized a modeling strategy based on Dynamic FBA (DFBA), called Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), to satisfy these specifications. Our strategy adds constraints describing the dynamics and regulation of metabolism that are strictly linear. We evaluated LK-DFBA against alternative modeling frameworks using simulated noisy data from a small in silico model and a larger model of central carbon metabolism in E. coli, and compared each framework's ability to recapitulate the original system. Results: In the smaller model, we found that we could use regression from a dynamic flux estimation (DFE) with an optional non-linear parameter optimization to reproduce metabolite concentration dynamic trends more effectively than an ordinary differential equation model with generalized mass action rate laws when tested under realistic data sampling frequency and noise levels. We observed detrimental effects across all tested modeling approaches when metabolite time course data were missing, but found these effects to be smaller for LK-DFBA in most cases. With the E. coli model, we produced qualitatively reasonable results with similar properties to the smaller model and explored two different parameterization structures that yield trade-offs in computation time and accuracy. Conclusions: LK-DFBA allows for calculation of metabolite concentrations and considers metabolite-dependent regulation while still retaining many computational advantages of FBA. This provides the proof-of-principle for a new metabolic modeling framework with the potential to create genome-scale dynamic models and the potential to be applied in strain engineering tools that currently use FBA. [ABSTRACT FROM AUTHOR]

Background: Constraint-based metabolic modeling has been applied to understand metabolism related disease mechanisms, to predict potential new drug targets and anti-metabolites, and to identify biomarkers of complex diseases. Although the state-of-art modeling toolbox, COBRA 3.0, is powerful, it requires substantial computing time conducting flux balance analysis, knockout analysis, and Markov Chain Monte Carlo (MCMC) sampling, which may limit its application in large scale genome-wide analysis. Results: Here, we rewrote the underlying code of COBRA 3.0 using C/C++, and developed a toolbox, termed FastMM, to effectively conduct constraint-based metabolic modeling. The results showed that FastMM is 2~400 times faster than COBRA 3.0 in performing flux balance analysis and knockout analysis and returns consistent outputs. When applied to MCMC sampling, FastMM is 8 times faster than COBRA 3.0. FastMM is also faster than some efficient metabolic modeling applications, such as Cobrapy and Fast-SL. In addition, we developed a Matlab/Octave interface for fast metabolic modeling. This interface was fully compatible with COBRA 3.0, enabling users to easily perform complex applications for metabolic modeling. For example, users who do not have deep constraint-based metabolic model knowledge can just type one command in Matlab/Octave to perform personalized metabolic modeling. Users can also use the advance and multiple threading parameters for complex metabolic modeling. Thus, we provided an efficient and user-friendly solution to perform large scale genome-wide metabolic modeling. For example, FastMM can be applied to the modeling of individual cancer metabolic profiles of hundreds to thousands of samples in the Cancer Genome Atlas (TCGA). Conclusion: FastMM is an efficient and user-friendly toolbox for large-scale personalized constraint-based metabolic modeling. It can serve as a complementary and invaluable improvement to the existing functionalities in COBRA 3.0. FastMM is under GPL license and can be freely available at GitHub site: https://github.com/GonghuaLi/FastMM. [ABSTRACT FROM AUTHOR]

Background: In order to improve the accuracy of constraint-based metabolic models, several approaches have been developed which intend to integrate additional biological information. Two of these methods, MOMENT and GECKO, incorporate enzymatic (kcat) parameters and enzyme mass constraints to further constrain the space of feasible metabolic flux distributions. While both methods have been proven to deliver useful extensions of metabolic models, they may considerably increase size and complexity of the models and there is currently no tool available to fully automate generation and calibration of such enzyme-constrained models from given stoichiometric models. Results: In this work we present three major developments. We first conceived short MOMENT (sMOMENT), a simplified version of the MOMENT approach, which yields the same predictions as MOMENT but requires significantly fewer variables and enables direct inclusion of the relevant enzyme constraints in the standard representation of a constraint-based model. When measurements of enzyme concentrations are available, these can be included as well leading in the extreme case, where all enzyme concentrations are known, to a model representation that is analogous to the GECKO approach. Second, we developed the AutoPACMEN toolbox which allows an almost fully automated creation of sMOMENT-enhanced stoichiometric metabolic models. In particular, this includes the automatic read-out and processing of relevant enzymatic data from different databases and the reconfiguration of the stoichiometric model with embedded enzymatic constraints. Additionally, tools have been developed to adjust (kcat and enzyme pool) parameters of sMOMENT models based on given flux data. We finally applied the new sMOMENT approach and the AutoPACMEN toolbox to generate an enzyme-constrained version of the E. coli genome-scale model iJO1366 and analyze its key properties and differences with the standard model. In particular, we show that the enzyme constraints improve flux predictions (e.g., explaining overflow metabolism and other metabolic switches) and demonstrate, for the first time, that these constraints can markedly change the spectrum of metabolic engineering strategies for different target products. Conclusions: The methodological and tool developments presented herein pave the way for a simplified and routine construction and analysis of enzyme-constrained metabolic models. [ABSTRACT FROM AUTHOR]

Proteinaceous organic wastes are suitable substrates to produce high added‐value products in anaerobic mixed‐culture fermentations. In these processes, the stoichiometry of the biotransformation depends highly on operational conditions such as pH or feeding characteristics and there are still no tools that allow the process to be directed toward those products of interest. Indeed, the lack of product selectivity strongly limits the potential industrial development of these bioprocesses. In this work, we developed a mathematical metabolic model for the production of volatile fatty acids from protein‐rich wastes. In particular, the effect of pH on the product yields is analyzed and, for the first time, the observed changes are mechanistically explained. The model reproduces experimental results at both neutral and acidic pH and it is also capable of predicting the tendencies in product yields observed with a pH drop. It also offers mechanistic insights into the interaction among the different amino acids (AAs) of a particular protein and how an AA might yield different products depending on the relative abundance of other AAs. Particular emphasis is placed on the utility of this mathematical model as a process design tool and different examples are given on how to use the model for this purpose. [ABSTRACT FROM AUTHOR]

Cupriavidus necator is a bacterium with a high phenotypic diversity and versatile metabolic capabilities. It has been extensively studied as a model hydrogen oxidizer, as well as a producer of polyhydroxyalkanoates (PHA), plastic-like biopolymers with a high potential to substitute petroleum-based materials. Thanks to its adaptability to diverse metabolic lifestyles and to the ability to accumulate large amounts of PHA, C. necator is employed in many biotechnological processes, with particular focus on PHA production from waste carbon sources. The large availability of genomic information has enabled a characterization of C. necator 's metabolism, leading to the establishment of metabolic models which are used to devise and optimize culture conditions and genetic engineering approaches. In this work, the characteristics of available C. necator strains and genomes are reviewed, underlining how a thorough comprehension of the genetic variability of C. necator is lacking and it could be instrumental for wider application of this microorganism. The metabolic paradigms of C. necator and how they are connected to PHA production and accumulation are described, also recapitulating the variety of carbon substrates used for PHA accumulation, highlighting the most promising strategies to increase the yield. Finally, the review describes and critically analyzes currently available genome-scale metabolic models and reduced metabolic network applications commonly employed in the optimization of PHA production. Overall, it appears that the capacity of C. necator of performing CO 2 bioconversion to PHA is still underexplored, both in biotechnological applications and in metabolic modeling. However, the accurate characterization of this organism and the efforts in using it for gas fermentation can help tackle this challenging perspective in the future. • Taxonomy of Cupriavidus necator strains is examined by comparative genomics. • The metabolism of C. necator is suited to diverse aerobic and anaerobic lifestyles. • C. necator can convert many waste products to polyhydroxyalkanoates. • Metabolic modeling of C. necator is a powerful tool for bioprocess optimization. [ABSTRACT FROM AUTHOR]

Bacterial genes that change in expression upon environmental disturbance have commonly been seen as those that must also phenotypically matter. However, several studies suggest that differentially expressed genes are rarely phenotypically important. We demonstrate, for Gram-positive and Gram-negative bacteria, that these seemingly uncoordinated gene sets are involved in responses that can be linked through topological network analysis. However, the level of coordination is stress dependent. While a well-coordinated response is triggered in response to nutrient stress, antibiotics trigger an uncoordinated response in which transcriptionally and phenotypically important genes are neither linked spatially nor in their magnitude. Moreover, a gene expression meta-analysis reveals that genes with large fitness changes during stress have low transcriptional variation across hundreds of other conditions, and vice versa. Our work suggests that cellular responses can be understood through network models that incorporate regulatory and genetic relationships, which could aid drug target predictions and genetic network engineering.

Antibiotic‐associated infections with Clostridioides difficile are a severe and often lethal risk for hospitalized patients, and can also affect populations without these classical risk factors. For a rational design of therapeutical concepts, a better knowledge of the metabolism of the pathogen is crucial. Metabolic modeling can provide a simulation of quantitative growth and usage of metabolic pathways, leading to a deeper understanding of the organism. Here, we present an elaborate genome‐scale metabolic model of C. difficile 630Δerm. The model iHD992 includes experimentally determined product and substrate uptake rates and is able to simulate the energy metabolism and quantitative growth of C. difficile. Dynamic flux balance analysis was used for time‐resolved simulations of the quantitative growth in two different media. The model predicts oxidative Stickland reactions and glucose degradation as main sources of energy, while the resulting reduction potential is mostly used for acetogenesis via the Wood–Ljungdahl pathway. Initial modeling experiments did not reproduce the observed growth behavior before the production of large quantities of a previously unknown polysaccharide was detected. Combined genome analysis and laboratory experiments indicated that the polysaccharide is an acetylated glucose polymer. Time‐resolved simulations showed that polysaccharide secretion was coupled to growth even during unstable glucose uptake in minimal medium. This is accomplished by metabolic shifts between active glycolysis and gluconeogenesis which were also observed in laboratory experiments.

The genomic era has revolutionized research on secondary metabolites and bioinformatics methods have in recent years revived the antibiotic discovery process after decades with only few new active molecules being identified. New computational tools are driven by genomics and metabolomics analysis, and enables rapid identification of novel secondary metabolites. To translate this increased discovery rate into industrial exploitation, it is necessary to integrate secondary metabolite pathways in the metabolic engineering process. In this review, we will describe the novel advances in discovery of secondary metabolites produced by filamentous fungi, highlight the utilization of genome-scale metabolic models (GEMs) in the design of fungal cell factories for the production of secondary metabolites and review strategies for optimizing secondary metabolite production through the construction of high yielding platform cell factories.

ABSTRACT 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 adaptation to the host cell is important for obligate intracellular pathogens such as Chlamydia trachomatis (Ct). Here we infer the flux differences for Ct from proteome and qRT-PCR data by comprehensive pathway modeling. We compare the comparatively inert infectious elementary body (EB) and the active replicative reticulate body (RB) systematically using a genome-scale metabolic model with 321 metabolites and 277 reactions. This did yield 84 extreme pathways based on a published proteomics dataset at three different time points of infection. Validation of predictions was done by quantitative RT-PCR of enzyme mRNA expression at three time points. Ct’s major active pathways are glycolysis, gluconeogenesis, glycerol-phospholipid (GPL) biosynthesis (support from host acetyl-CoA) and pentose phosphate pathway (PPP), while its incomplete TCA and fatty acid biosynthesis are less active. The modeled metabolic pathways are much more active in RB than in EB. Our in silico model suggests that EB and RB utilize folate to generate NAD(P)H using independent pathways. The only low metabolic flux inferred for EB involves mainly carbohydrate metabolism. RB utilizes energy -rich compounds to generate ATP in nucleic acid metabolism. Validation data for the modeling include proteomics experiments (model basis) as well as qRT-PCR confirmation of selected metabolic enzyme mRNA expression differences. The metabolic modeling is made fully available here. Its detailed insights and models on Ct metabolic adaptations during infection are a useful modeling basis for future studies.

ABSTRACT The oomycete pathogen Phytophthora infestans causes potato and tomato late blight, a disease that is a serious threat to agriculture. P. infestans is a hemibiotrophic pathogen, and during infection, it scavenges nutrients from living host cells for its own proliferation. To date, the nutrient flux from host to pathogen during infection has hardly been studied, and the interlinked metabolisms of the pathogen and host remain poorly understood. Here, we reconstructed an integrated metabolic model of P. infestans and tomato (Solanum lycopersicum) by integrating two previously published models for both species. We used this integrated model to simulate metabolic fluxes from host to pathogen and explored the topology of the model to study the dependencies of the metabolism of P. infestans on that of tomato. This showed, for example, that P. infestans, a thiamine auxotroph, depends on certain metabolic reactions of the tomato thiamine biosynthesis. We also exploited dual-transcriptome data of a time course of a full late blight infection cycle on tomato leaves and integrated the expression of metabolic enzymes in the model. This revealed profound changes in pathogen-host metabolism during infection. As infection progresses, P. infestans performs less de novo synthesis of metabolites and scavenges more metabolites from tomato. This integrated metabolic model for the P. infestans-tomato interaction provides a framework to integrate data and generate hypotheses about in planta nutrition of P. infestans throughout its infection cycle. IMPORTANCE Late blight disease caused by the oomycete pathogen Phytophthora infestans leads to extensive yield losses in tomato and potato cultivation worldwide. To effectively control this pathogen, a thorough understanding of the mechanisms shaping the interaction with its hosts is paramount. While considerable work has focused on exploring host defense mechanisms and identifying P. infestans proteins contributing to virulence and pathogenicity, the nutritional strategies of the pathogen are mostly unresolved. Genome-scale metabolic models (GEMs) can be used to simulate metabolic fluxes and help in unravelling the complex nature of metabolism. We integrated a GEM of tomato with a GEM of P. infestans to simulate the metabolic fluxes that occur during infection. This yields insights into the nutrients that P. infestans obtains during different phases of the infection cycle and helps in generating hypotheses about nutrition in planta.

ABSTRACT A better understanding of essential cellular functions in pathogenic bacteria is important for the development of more effective antimicrobial agents. We performed a comprehensive identification of essential genes in Mycobacterium tuberculosis, the major causative agent of tuberculosis, using a combination of transposon insertion sequencing (Tn-seq) and comparative genomic analysis. To identify conditionally essential genes by Tn-seq, we used media with different nutrient compositions. Although many conditional gene essentialities were affected by the presence of relevant nutrient sources, we also found that the essentiality of genes in a subset of metabolic pathways was unaffected by metabolite availability. Comparative genomic analysis revealed that not all essential genes identified by Tn-seq were fully conserved within the M. tuberculosis complex, including some existing antitubercular drug target genes. In addition, we utilized an available M. tuberculosis genome-scale metabolic model, iSM810, to predict M. tuberculosis gene essentiality in silico. Comparing the sets of essential genes experimentally identified by Tn-seq to those predicted in silico reveals the capabilities and limitations of gene essentiality predictions, highlighting the complexity of M. tuberculosis essential metabolic functions. This study provides a promising platform to study essential cellular functions in M. tuberculosis. IMPORTANCE Mycobacterium tuberculosis causes 10 million cases of tuberculosis (TB), resulting in over 1 million deaths each year. TB therapy is challenging because it requires a minimum of 6 months of treatment with multiple drugs. Protracted treatment times and the emergent spread of drug-resistant M. tuberculosis necessitate the identification of novel targets for drug discovery to curb this global health threat. Essential functions, defined as those indispensable for growth and/or survival, are potential targets for new antimicrobial drugs. In this study, we aimed to define gene essentialities of M. tuberculosis on a genomewide scale to comprehensively identify potential targets for drug discovery. We utilized a combination of experimental (functional genomics) and in silico approaches (comparative genomics and flux balance analysis). Our functional genomics approach identified sets of genes whose essentiality was affected by nutrient availability. Comparative genomics revealed that not all essential genes were fully conserved within the M. tuberculosis complex. Comparing sets of essential genes identified by functional genomics to those predicted by flux balance analysis highlighted gaps in current knowledge regarding M. tuberculosis metabolic capabilities. Thus, our study identifies numerous potential antitubercular drug targets and provides a comprehensive picture of the complexity of M. tuberculosis essential cellular functions.

The biosynthetic capabilities of microbes underlie their growth and interactions, playing a prominent role in microbial community structure. For large, diverse microbial communities, prediction of these capabilities is limited by uncertainty about metabolic functions and environmental conditions. To address this challenge, we propose a probabilistic method, inspired by percolation theory, to computationally quantify how robustly a genome-derived metabolic network produces a given set of metabolites under an ensemble of variable environments. We used this method to compile an atlas of predicted biosynthetic capabilities for 97 metabolites across 456 human oral microbes. This atlas captures taxonomically-related trends in biomass composition, and makes it possible to estimate inter-microbial metabolic distances that correlate with microbial co-occurrences. We also found a distinct cluster of fastidious/uncultivated taxa, including several Saccharibacteria (TM7) species, characterized by their abundant metabolic deficiencies. By embracing uncertainty, our approach can be broadly applied to understanding metabolic interactions in complex microbial ecosystems.

ABSTRACT Cystic fibrosis (CF) is a fatal genetic disease characterized by chronic lung infections due to aberrant mucus production and the inability to clear invading pathogens. The traditional view that CF infections are caused by a single pathogen has been replaced by the realization that the CF lung usually is colonized by a complex community of bacteria, fungi, and viruses. To help unravel the complex interplay between the CF lung environment and the infecting microbial community, we developed a community metabolic model comprised of the 17 most abundant bacterial taxa, which account for >95% of reads across samples, from three published studies in which 75 sputum samples from 46 adult CF patients were analyzed by 16S rRNA gene sequencing. The community model was able to correctly predict high abundances of the “rare” pathogens Enterobacteriaceae, Burkholderia, and Achromobacter in three patients whose polymicrobial infections were dominated by these pathogens. With these three pathogens removed, the model correctly predicted that the remaining 43 patients would be dominated by Pseudomonas and/or Streptococcus. This dominance was predicted to be driven by relatively high monoculture growth rates of Pseudomonas and Streptococcus as well as their ability to efficiently consume amino acids, organic acids, and alcohols secreted by other community members. Sample-by-sample heterogeneity of community composition could be qualitatively captured through random variation of the simulated metabolic environment, suggesting that experimental studies directly linking CF lung metabolomics and 16S sequencing could provide important insights into disease progression and treatment efficacy. IMPORTANCE Cystic fibrosis (CF) is a genetic disease in which chronic airway infections and lung inflammation result in respiratory failure. CF airway infections are usually caused by bacterial communities that are difficult to eradicate with available antibiotics. Using species abundance data for clinically stable adult CF patients assimilated from three published studies, we developed a metabolic model of CF airway communities to better understand the interactions between bacterial species and between the bacterial community and the lung environment. Our model predicted that clinically observed CF pathogens could establish dominance over other community members across a range of lung nutrient conditions. Heterogeneity of species abundances across 75 patient samples could be predicted by assuming that sample-to-sample heterogeneity was attributable to random variations in the CF nutrient environment. Our model predictions provide new insights into the metabolic determinants of pathogen dominance in the CF lung and could facilitate the development of improved treatment strategies.

Metabolic networks have been used to successfully predict phenotypes based on optimization principles. However, a general framework that would extend to situations not governed by simple optimization, such as multispecies communities, is still lacking. Concepts from evolutionary game theory have been proposed to amend the situation. Alternative metabolic states can be seen as strategies in a “metabolic game,” and phenotypes can be predicted based on the equilibria of this game. In this survey, we review the literature on applying game theory to the study of metabolism, present the general idea of a metabolic game, and discuss open questions and future challenges.

The acyl-homoserine lactones (AHLs)-mediated LuxI/LuxR quorum sensing (QS) system orchestrates diverse bacterial behaviors in response to changes in population density. The role of the BjaI/BjaR 1 QS system in Bradyrhizobium diazoefficiens USDA 110, which shares homology with LuxI/LuxR, remains elusive during symbiotic interaction with soybean. Here this genetic system in wild-type (WT) bacteria residing inside nodules exhibited significantly reduced activity compared to free-living cells, potentially attributed to soybean-mediated suppression. The deletion mutant strain ΔbjaR 1 showed significantly enhanced nodulation induction and nitrogen fixation ability. Nevertheless, its ultimate symbiotic outcome (plant dry weight) in soybeans was compromised. Furthermore, comparative analysis of the transcriptome, proteome, and promoter activity revealed that the inactivation of BjaR 1 systematically activated and inhibited genomic modules associated with nodulation and nitrogen metabolism. The former appeared to be linked to a significant decrease in the expression of NodD2, a key cell-density-dependent repressor of nodulation genes, while the latter conferred bacterial growth and nitrogen fixation insensitivity to environmental nitrogen. In addition, BjaR 1 exerted a positive influence on the transcription of multiple genes involved in a so-called central intermediate metabolism within the nodule. In conclusion, our findings highlight the crucial role of the BjaI/BjaR 1 QS circuit in positively regulating bacterial nitrogen metabolism and emphasize the significance of the soybean-mediated suppression of this genetic system for promoting efficient symbiotic nitrogen fixation by B. diazoefficiens .IMPORTANCEThe present study demonstrates, for the first time, that the BjaI/BjaR1 QS system of Bradyrhizobium diazoefficiens has a significant impact on its nodulation and nitrogen fixation capability in soybean by positively regulating NodD2 expression and bacterial nitrogen metabolism. Moreover, it provides novel insights into the importance of suppressing the activity of this QS circuit by the soybean host plant in establishing an efficient mutual relationship between the two symbiotic partners. This research expands our understanding of legumes' role in modulating symbiotic nitrogen fixation through rhizobial QS-mediated metabolic functioning, thereby deepening our comprehension of symbiotic coevolution theory. In addition, these findings may hold great promise for developing quorum quenching technology in agriculture., Competing Interests: The authors declare no conflict of interest.

Understanding the ancestral transition from anaerobic to aerobic lifestyles is essential for comprehending life's early evolution. However, the biological adaptations occurring during this crucial transition remain largely unexplored. Thiamine is an important cofactor involved in central carbon metabolism and aerobic respiration. Here, we explored the phylogenetic and global distribution of thiamine-auxotrophic and thiamine-prototrophic bacteria based on the thiamine biosynthetic pathway in 154 838 bacterial genomes. We observed strong coincidences of the origin of thiamine-synthetic bacteria with the "Great Oxygenation Event," indicating that thiamine biosynthesis in bacteria emerged as an adaptation to aerobic respiration. Furthermore, we demonstrated that thiamine-mediated metabolic interactions are fundamental factors influencing the assembly and diversity of bacterial communities by a global survey across 4245 soil samples. Through our newly established stable isotope probing-metabolic modeling method, we uncovered the active utilization of thiamine-mediated metabolic interactions by bacterial communities in response to changing environments, thus revealing an environmental adaptation strategy employed by bacteria at the community level. Our study demonstrates the widespread thiamine-mediated metabolic interactions in bacterial communities and their crucial roles in setting the stage for an evolutionary transition from anaerobic to aerobic lifestyles and subsequent environmental adaptation. These findings provide new insights into early bacterial evolution and their subsequent growth and adaptations to environments., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Society for Microbial Ecology.)

Metabolic modeling and machine learning (ML) are crucial components of the evolving next-generation tools in systems and synthetic biology, aiming to unravel the intricate relationship between genotype, phenotype, and the environment. Nonetheless, the comprehensive exploration of integrating these two frameworks, and fully harnessing the potential of fluxomic data, remains an unexplored territory. In this study, we present, rigorously evaluate, and compare ML-based techniques for data integration. The hybrid model revealed that the overexpression of six target genes and the knockout of seven target genes contribute to enhanced ethanol production. Specifically, we investigated the influence of succinate dehydrogenase (SDH) on ethanol biosynthesis in Saccharomyces cerevisiae through shake flask experiments. The findings indicate a noticeable increase in ethanol yield, ranging from 6 % to 10 %, in SDH subunit gene knockout strains compared to the wild-type strain. Moreover, in pursuit of a high-yielding strain for ethanol production, dual-gene deletion experiments were conducted targeting glycerol-3-phosphate dehydrogenase (GPD) and SDH. The results unequivocally demonstrate significant enhancements in ethanol production for the engineered strains Δsdh 4 Δgpd 1, Δsdh 5 Δgpd 1, Δsdh 6 Δgpd 1, Δsdh 4 Δgpd 2, Δsdh 5 Δgpd 2, and Δsdh 6 Δgpd 2, with improvements of 21.6 %, 27.9 %, and 22.7 %, respectively. Overall, the results highlighted that integrating mechanistic flux features substantially improves the prediction of gene knockout strains not accounted for in metabolic reconstructions. In addition, the finding in this study delivers valuable tools for comprehending and manipulating intricate phenotypes, thereby enhancing prediction accuracy and facilitating deeper insights into mechanistic aspects within the field of synthetic biology., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

Importance: Co-occurrence network analysis is an effective tool for predicting complex networks of microbial interactions in the natural environment. Using isolates from a rhizosphere, we constructed multi-species biofilm communities and investigated co-occurrence patterns between microbial species in genome-scale metabolic models and in vitro experiments. According to our results, metabolic exchanges and resource competition may partially explain the co-occurrence network analysis results found in synthetic bacterial biofilm communities., Competing Interests: The authors declare no conflict of interest.

Importance: Lyme disease is often treated using long courses of antibiotics, which can cause side effects for patients and risks the evolution of antimicrobial resistance. Narrow-spectrum antimicrobials would reduce these risks, but their development has been slow because the Lyme disease bacterium, Borrelia burgdorferi , is difficult to work with in the laboratory. To accelerate the drug discovery pipeline, we developed a computational model of B. burgdorferi 's metabolism and used it to predict essential enzymatic reactions whose inhibition prevented growth in silico . These predictions were validated using small-molecule enzyme inhibitors, several of which were shown to have specific activity against B. burgdorferi . Although the specific compounds used are not suitable for clinical use, we aim to use them as lead compounds to develop optimized drugs targeting the pathways discovered here., Competing Interests: The authors declare no conflict of interest.