Search

Blood metabolite levels are affected by numerous factors, including preanalytical factors such as collection methods and geographical sites. These perturbations have caused deleterious consequences for many metabolomics studies and represent a major challenge in the metabolomics field. It is important to understand these factors and develop models to reduce their perturbations. However, to date, the lack of suitable mathematical models for blood metabolite levels under homeostasis has hindered progress. In this study, we develop quantitative models of blood metabolite levels in healthy adults based on multisite sample cohorts that mimic the current challenge. Five cohorts of samples obtained across four geographically distinct sites were investigated, focusing on approximately 50 metabolites that were quantified using 1H NMR spectroscopy. More than one-third of the variation in these metabolite profiles is due to cross-cohort variation. A dramatic reduction in the variation of metabolite levels (90%), especially their site-to-site variation (95%), was achieved by modeling each metabolite using demographic and clinical factors and especially other metabolites, as observed in the top principal components. The results also reveal that several metabolites contribute disproportionately to such variation, which could be explained by their association with biological pathways including biosynthesis and degradation. The study demonstrates an intriguing network effect of metabolites that can be utilized to better define homeostatic metabolite levels, which may have implications for improved health monitoring. As an example of the potential utility of the approach, we show that modeling gender-related metabolic differences retains the interesting variance while reducing unwanted (site-related) variance.

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities. ONE SENTENCE SUMMARY: Computer-assisted biosynthetic capacities elucidation accelerates our ability to interpret microbial interactions, allowing us to understand better and establish a balance within ecosystems. [ABSTRACT FROM AUTHOR]

To explore the impact of microbial interactions on outcomes from three prevalent algorithms (Flux Balance Analysis (FBA), community FBA (cFBA), and SteadyCom) analyzing microbial community metabolic networks, five toy community models representing common microbial interactions were designed. These include commensalism, mutualism, competition, mutualism-competition, and commensalism-competition. Various scenarios, considering different biomass yields and substrate constraints, were examined for each type. In commensal communities, all algorithms consistently produced similar results. However, changes in biomass yields and substrate constraints led to variable abundances (0.33–0.8) and community growth rates (2–5 1/h) within a broad range. For competitive communities, all algorithms predicted growth of fastest-growing member. To comply with the natural coexistence of members, suboptimal solutions over optimal point are recommended. FBA faced challenges in modeling mutualism, consistently predicting growth of only one member. Although cFBA and SteadyCom resulted in a lower community growth rate, coexistence of both members were satisfied. In toy models with dual interactions, more realistic outcomes were achieved contrary to purely competitive model as the dependency fosters the coexistence which was missing in the competitive only scenarios. These findings emphasize the importance of algorithm choice based on specific microbial interaction types for reliable community behavior predictions.​ [ABSTRACT FROM AUTHOR]

Characterizing the phenotypic diversity and metabolic capabilities of industrially relevant manufacturing cell lines is critical to bioprocess optimization and cell line development. Metabolic capabilities of production hosts limit nutrient and resource channeling into desired cellular processes and can have a profound impact on productivity. These limitations cannot be directly inferred from measured data such as spent media concentrations or transcriptomics. Here, we present an integrated multi-omic analysis pipeline combining exo-metabolomics, transcriptomics, and genome-scale metabolic network analysis and apply it to three antibody-producing Chinese Hamster Ovary cell lines to identify reprogramming features associated with high-producing clones and metabolic bottlenecks limiting product formation in an industrial bioprocess. Analysis of individual datatypes revealed a decreased nitrogenous byproduct secretion in high-producing clones and the topological changes in peripheral metabolic pathway expression associated with phase shifts. An integrated omics analysis in the context of the genome-scale metabolic model elucidated the differences in central metabolism and identified amino acid utilization bottlenecks limiting cell growth and antibody production that were not evident from exo-metabolomics or transcriptomics alone. Thus, we demonstrate the utility of a multi-omics characterization in providing an in-depth understanding of cellular metabolism, which is critical to efforts in cell engineering and bioprocess optimization. • Multi-omics analysis permits synergistic evaluation of culture performance. • Overall phenotypic changes were traced back to phase-specific metabolic reprogramming. • Reprogramming nitrogen usage and metabolism was associated with increased productivity. • Amino acid usage and channeling towards desired products limited productivity. [ABSTRACT FROM AUTHOR]

Abstract Background Due to increasing ecological concerns, microbial production of biochemicals from sustainable carbon sources like acetate is rapidly gaining importance. However, to successfully establish large-scale production scenarios, a solid understanding of metabolic driving forces is required to inform bioprocess design. To generate such knowledge, we constructed isopropanol-producing Escherichia coli W strains. Results Based on strain screening and metabolic considerations, a 2-stage process was designed, incorporating a growth phase followed by a nitrogen-starvation phase. This process design yielded the highest isopropanol titers on acetate to date (13.3 g L−1). Additionally, we performed shotgun and acetylated proteomics, and identified several stress conditions in the bioreactor scenarios, such as acid stress and impaired sulfur uptake. Metabolic modeling allowed for an in-depth characterization of intracellular flux distributions, uncovering cellular demand for ATP and acetyl-CoA as limiting factors for routing carbon toward the isopropanol pathway. Moreover, we asserted the importance of a balance between fluxes of the NADPH-providing isocitrate dehydrogenase (ICDH) and the product pathway. Conclusions Using the newly gained system-level understanding for isopropanol production from acetate, we assessed possible engineering approaches and propose process designs to maximize production. Collectively, our work contributes to the establishment and optimization of acetate-based bioproduction systems. Graphical Abstract

Background: Due to increasing ecological concerns, microbial production of biochemicals from sustainable carbon sources like acetate is rapidly gaining importance. However, to successfully establish large-scale production scenarios, a solid understanding of metabolic driving forces is required to inform bioprocess design. To generate such knowledge, we constructed isopropanol-producing Escherichia coli W strains. Results: Based on strain screening and metabolic considerations, a 2-stage process was designed, incorporating a growth phase followed by a nitrogen-starvation phase. This process design yielded the highest isopropanol titers on acetate to date (13.3 g L−1). Additionally, we performed shotgun and acetylated proteomics, and identified several stress conditions in the bioreactor scenarios, such as acid stress and impaired sulfur uptake. Metabolic modeling allowed for an in-depth characterization of intracellular flux distributions, uncovering cellular demand for ATP and acetyl-CoA as limiting factors for routing carbon toward the isopropanol pathway. Moreover, we asserted the importance of a balance between fluxes of the NADPH-providing isocitrate dehydrogenase (ICDH) and the product pathway. Conclusions: Using the newly gained system-level understanding for isopropanol production from acetate, we assessed possible engineering approaches and propose process designs to maximize production. Collectively, our work contributes to the establishment and optimization of acetate-based bioproduction systems. [ABSTRACT FROM AUTHOR]

Simple Summary: This study describes a large-scale analysis of cancer metabolism across The Cancer Genome Atlas and Multiple Myeloma Research Foundation MMRF-COMMPASS study encompassing 10,915 patients from 34 different cancer types. To the best of our knowledge, this is the first such large-scale metabolic analysis performed using these data. Our approach identified biologically interpretable metabolic biomarkers of overall survival in 31 out of 34 cancer types evaluated. Additionally, in a subset of 7 cancer types, we predict chemosensitivity to antifolate-based therapies and a biomarker of response based on our patient-specific metabolic models. Finally, this work introduces the concept of "interpretation-driven biomarker identification". In this framework, biological interpretability of biomarker signatures is ensured first by using network-based simulations to characterize tumor biology, then using survival analysis to identify any biological features associated with overall survival—essentially decoupling biomarker discovery and biological interpretation of biomarker scores. Despite the high variability in cancer biology, cancers nevertheless exhibit cohesive hallmarks across multiple cancer types, notably dysregulated metabolism. Metabolism plays a central role in cancer biology, and shifts in metabolic pathways have been linked to tumor aggressiveness and likelihood of response to therapy. We therefore sought to interrogate metabolism across cancer types and understand how intrinsic modes of metabolism vary within and across indications and how they relate to patient prognosis. We used context specific genome-scale metabolic modeling to simulate metabolism across 10,915 patients from 34 cancer types from The Cancer Genome Atlas and the MMRF-COMMPASS study. We found that cancer metabolism clustered into modes characterized by differential glycolysis, oxidative phosphorylation, and growth rate. We also found that the simulated activities of metabolic pathways are intrinsically prognostic across cancer types, especially tumor growth rate, fatty acid biosynthesis, folate metabolism, oxidative phosphorylation, steroid metabolism, and glutathione metabolism. This work shows the prognostic power of individual patient metabolic modeling across multiple cancer types. Additionally, it shows that analyzing large-scale models of cancer metabolism with survival information provides unique insights into underlying relationships across cancer types and suggests how therapies designed for one cancer type may be repurposed for use in others. [ABSTRACT FROM AUTHOR]

Metabolic engineering for high productivity and increased robustness is needed to enable sustainable biomanufacturing of lactic acid from lignocellulosic biomass. Lactic acid is an important commodity chemical used for instance as a monomer for production of polylactic acid, a biodegradable polymer. Here, rational and model-based optimization was used to engineer a diploid, xylose fermenting Saccharomyces cerevisiae strain to produce L-lactic acid. The metabolic flux was steered towards lactic acid through the introduction of multiple lactate dehydrogenase encoding genes while deleting ERF2 , GPD1 , and CYB2. A production of 93 g/L of lactic acid with a yield of 0.84 g/g was achieved using xylose as the carbon source. To increase xylose utilization and reduce acetic acid synthesis, PHO13 and ALD6 were also deleted from the strain. Finally, CDC19 encoding a pyruvate kinase was overexpressed, resulting in a yield of 0.75 g lactic acid/g sugars consumed, when the substrate used was a synthetic lignocellulosic hydrolysate medium, containing hexoses, pentoses and inhibitors such as acetate and furfural. Notably, modeling also provided leads for understanding the influence of oxygen in lactic acid production. High lactic acid production from xylose, at oxygen-limitation could be explained by a reduced flux through the oxidative phosphorylation pathway. On the contrast, higher oxygen levels were beneficial for lactic acid production with the synthetic hydrolysate medium, likely as higher ATP concentrations are needed for tolerating the inhibitors therein. The work highlights the potential of S. cerevisiae for industrial production of lactic acid from lignocellulosic biomass. [Display omitted] • A xylose fermenting S. cerevisiae strain was engineered to produce L-lactic acid. • Rational and model-guided design was applied to increase lactic acid production. • 93 g/L of lactic acid at a yield of 0.84 g/g was produced from xylose. • A yield of 0.75 g/g was achieved in a synthetic lignocellulosic hydrolysate. • Modeling was used to understand the influence of oxygen in lactic acid production. [ABSTRACT FROM AUTHOR]

The current study investigated the impact of Saccharomyces cerevisiae fermentation on apricot kernel oil extraction and the resulting quality characteristics of the oil. The seed kernels underwent physicochemical and antioxidant studies after being separated, gelatinized, and fermented. The results revealed 43.6 % of oil yield by fermentation. The extracted oil had a lower acid value (3.74 mg KOH/g), and peroxide value (4.25 meq/kg). The investigation of additional oil properties was included, like iodine value (100.26 mg KOH/g), saponification value (190.5 mg KOH/g), specific gravity (0.919), viscosity (4.62), moisture content (0.592 %), and refractive index (1.4836), which signifies optimum physicochemical parameters of oil. As per the FTIR data, some distinct functional groups like C-H stretch and C O, C C, CH3, OH, and R-O-R esters were found in the extracted oil. The oil's fatty acid composition, with 50.47 % saturated, 22.23 % mono-unsaturated, and 35.13 % polyunsaturated fatty acids, enhances stability, oxidative resistance, and nutritional value, justifying its favorable physicochemical properties. The DPPH radical scavenging activity of the oil showed to have good antioxidant activity. Metabolic modeling of S. cerevisiae provided information about compounds received after fermentation. Hence, it can be concluded that fermentation is a suitable approach to extract apricot seed kernel oil with beneficial properties. [Display omitted] • Apricot seed kernel is a good source of oil, we can use it for multiple purposes. • Saccharomyces cerevisiae -based fermentation, we can get 43.6 % of oil. • Different beneficial fatty acids were obtained from this oil: Caprylic, Lauric, Margaric, Paullinic, Phthalic acid, etc. [ABSTRACT FROM AUTHOR]

This study presents a novel approach for developing generic metabolic Raman calibration models for in‐line cell culture analysis using glucose and lactate stock solution titration in an aqueous phase and data augmentation techniques. First, a successful set‐up of the titration method was achieved by adding glucose or lactate solution at several different constant rates into the aqueous phase of a bench‐top bioreactor. Subsequently, the in‐line glucose and lactate concentration were calculated and interpolated based on the rate of glucose and lactate addition, enabling data augmentation and enhancing the robustness of the metabolic calibration model. Nine different combinations of spectra pretreatment, wavenumber range selection, and number of latent variables were evaluated and optimized using aqueous titration data as training set and a historical cell culture data set as validation and prediction set. Finally, Raman spectroscopy data collected from 11 historical cell culture batches (spanning four culture modes and scales ranging from 3 to 200 L) were utilized to predict the corresponding glucose and lactate values. The results demonstrated a high prediction accuracy, with an average root mean square errors of prediction of 0.65 g/L for glucose, and 0.48 g/L for lactate. This innovative method establishes a generic metabolic calibration model, and its applicability can be extended to other metabolites, reducing the cost of deploying real‐time cell culture monitoring using Raman spectroscopy in bioprocesses. [ABSTRACT FROM AUTHOR]

Summary: Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST‐MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast‐growing algae, and produced a synergistic multi‐organ flux map that is a first in MFA for any biological system. These insights could not be elucidated through other approaches and show the potential of INST‐MFA to correct an oversimplified understanding of plant metabolism. [ABSTRACT FROM AUTHOR]

We review the process of reconstructing genome-scale metabolic models (GEMs) for single species and microbial communities. We present a detailed compendium of model-building platforms for experts and non-experts. We review and compare the approaches for testing GEMs towards engineering microbial systems. We discuss current challenges and future directions in model development, including the need for standards in microbial community models. Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels. [ABSTRACT FROM AUTHOR]

Itaconic acid is a platform chemical with a range of applications in polymer synthesis and is also discussed for biofuel production. While produced in industry from glucose or sucrose, co‐feeding of glucose and acetate was recently discussed to increase itaconic acid production by the smut fungus Ustilago maydis. In this study, we investigate the optimal co‐feeding conditions by interlocking experimental and computational methods. Flux balance analysis indicates that acetate improves the itaconic acid yield up to a share of 40% acetate on a carbon molar basis. A design of experiment results in the maximum yield of 0.14 itaconic acid per carbon source from 100 g L−1 $\,\text{g L}{}^{-1}$ glucose and 12 g L−1 $\,\text{g L}{}^{-1}$ acetate. The yield is improved by around 22% when compared to feeding of glucose as sole carbon source. To further improve the yield, gene deletion targets are discussed that were identified using the metabolic optimization tool OptKnock. The study contributes ideas to reduce land use for biotechnology by incorporating acetate as co‐substrate, a C2‐carbon source that is potentially derived from carbon dioxide. [ABSTRACT FROM AUTHOR]

ABSTRACT Streptococcus pyogenes is responsible for a range of diseases in humans contributing significantly to morbidity and mortality. Among more than 200 serotypes of S. pyogenes, serotype M1 strains hold the greatest clinical relevance due to their high prevalence in severe human infections. To enhance our understanding of pathogenesis and discovery of potential therapeutic approaches, we have developed the first genome-scale metabolic model (GEM) for a serotype M1 S. pyogenes strain, which we name iYH543. The curation of iYH543 involved cross-referencing a draft GEM of S. pyogenes serotype M1 from the AGORA2 database with gene essentiality and autotrophy data obtained from transposon mutagenesis-based and growth screens. We achieved a 92.6% (503/543 genes) accuracy in predicting gene essentiality and a 95% (19/20 amino acids) accuracy in predicting amino acid auxotrophy. Additionally, Biolog Phenotype microarrays were employed to examine the growth phenotypes of S. pyogenes, which further contributed to the refinement of iYH543. Notably, iYH543 demonstrated 88% accuracy (168/190 carbon sources) in predicting growth on various sole carbon sources. Discrepancies observed between iYH543 and the actual behavior of living S. pyogenes highlighted areas of uncertainty in the current understanding of S. pyogenes metabolism. iYH543 offers novel insights and hypotheses that can guide future research efforts and ultimately inform novel therapeutic strategies.IMPORTANCEGenome-scale models (GEMs) play a crucial role in investigating bacterial metabolism, predicting the effects of inhibiting specific metabolic genes and pathways, and aiding in the identification of potential drug targets. Here, we have developed the first GEM for the S. pyogenes highly virulent serotype, M1, which we name iYH543. The iYH543 achieved high accuracy in predicting gene essentiality. We also show that the knowledge obtained by substituting actual measurement values for iYH543 helps us gain insights that connect metabolism and virulence. iYH543 will serve as a useful tool for rational drug design targeting S. pyogenes metabolism and computational screening to investigate the interplay between inhibiting virulence factor synthesis and growth.

While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model's biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

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 Δsdh4Δgpd1, Δsdh5Δgpd1, Δsdh6Δgpd1, Δsdh4Δgpd2, Δsdh5Δgpd2, and Δsdh6Δgpd2, 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.

Genome-scale metabolic models comprehensively describe an organisms metabolism and can be tailored using omics data to model condition-specific physiology. The quality of context-specific models is impacted by (i) choice of algorithm and parameters and (ii) alternate context-specific models that equally explain the -omics data. Here we quantify the influence of alternate optima on microbial and mammalian model extraction using GIMME, iMAT, MBA, and mCADRE. We find that metabolic tasks defining an organisms phenotype must be explicitly and quantitatively protected. The scope of alternate models is strongly influenced by algorithm choice and the topological properties of the parent genome-scale model with fatty acid metabolism and intracellular metabolite transport contributing much to alternate solutions in all models. mCADRE extracted the most reproducible context-specific models and models generated using MBA had the most alternate solutions. There were fewer qualitatively different solutions generated by GIMME in E. coli, but these increased substantially in the mammalian models. Screening ensembles using a receiver operating characteristic plot identified the best-performing models. A comprehensive evaluation of models extracted using combinations of extraction methods and expression thresholds revealed that GIMME generated the best-performing models in E. coli, whereas mCADRE is better suited for complex mammalian models. These findings suggest guidelines for benchmarking -omics integration algorithms and motivate the development of a systematic workflow to enumerate alternate models and extract biologically relevant context-specific models.

Gliomas are the most common type of malignant brain tumors, with glioblastoma multiforme (GBM) having a median survival of 15 months due to drug resistance and relapse. The treatment of gliomas relies on surgery, radiotherapy and chemotherapy. Only 12 anti-brain tumor chemotherapies (AntiBCs), mostly alkylating agents, have been approved so far. Glioma subtype–specific metabolic models were reconstructed to simulate metabolite exchanges, in silico knockouts and the prediction of drug and drug combinations for all three subtypes. The simulations were confronted with literature, high-throughput screenings (HTSs), xenograft and clinical trial data to validate the workflow and further prioritize the drug candidates. The three subtype models accurately displayed different degrees of dependencies toward glutamine and glutamate. Furthermore, 33 single drugs, mainly antimetabolites and TXNRD1-inhibitors, as well as 17 drug combinations were predicted as potential candidates for gliomas. Half of these drug candidates have been previously tested in HTSs. Half of the tested drug candidates reduce proliferation in cell lines and two-thirds in xenografts. Most combinations were predicted to be efficient for all three glioma types. However, eflornithine/rifamycin and cannabidiol/adapalene were predicted specifically for GBM and low-grade glioma, respectively. Most drug candidates had comparable efficiency in preclinical tests, cerebrospinal fluid bioavailability and mode-of-action to AntiBCs. However, fotemustine and valganciclovir alone and eflornithine and celecoxib in combination with AntiBCs improved the survival compared to AntiBCs in two-arms, phase I/II and higher glioma clinical trials. Our work highlights the potential of metabolic modeling in advancing glioma drug discovery, which accurately predicted metabolic vulnerabilities, repurposable drugs and combinations for the glioma subtypes. [ABSTRACT FROM AUTHOR]

Cheese taste and flavour properties result from complex metabolic processes occurring in microbial communities. A deeper understanding of such mechanisms makes it possible to improve both industrial production processes and end-product quality through the design of microbial consortia. In this work, we caracterise the metabolism of a three-species community consisting of Lactococcus lactis , Lactobacillus plantarum and Propionibacterium freudenreichii during a seven-week cheese production process. Using genome-scale metabolic models and omics data integration, we modeled and calibrated individual dynamics using monoculture experiments, and coupled these models to capture the metabolism of the community. This model accurately predicts the dynamics of the community, enlightening the contribution of each microbial species to organoleptic compound production. Further metabolic exploration revealed additional possible interactions between the bacterial species. This work provides a methodological framework for the prediction of community-wide metabolism and highlights the added value of dynamic metabolic modeling for the comprehension of fermented food processes. [ABSTRACT FROM AUTHOR]

Population heterogeneity is due to intrinsic and extrinsic noises; cell to cell variability is related to cell partitioning at division and external forcing by the environment. In this article, a Monte‐Carlo calculation of integral source terms is coupled to a compartment approach. The biological phase is represented by an ensemble of discrete particles which are equipped with a coarse‐grained model that describes the intracellular response to extracellular conditions. The redistribution of cell properties at division generates unbalanced cells and introduces additional time scales into the population dynamics. Using this fully two‐way coupled approach, the consequences of repeated exposure to feast and famine events on both the reactor and the population heterogenities are investigated in a CSTR+PFR scale‐down bioreactor. It is shown that the integration of intrinsic noise due to cell division continuously produces cells that are not at equilibrium with their environment and this contributes to population heterogeneity, even in homogeneous environments. [ABSTRACT FROM AUTHOR]

Mechanistic, constraint-based models of microbial isolates or communities are a staple in the metabolic analysis toolbox, but predictions about microbe-microbe and microbe-environment interactions are only as good as the accuracy of transporter annotations. A number of hurdles stand in the way of comprehensive functional assignments for membrane transporters. These include general or non-specific substrate assignments, ambiguity in the localization, directionality and reversibility of a transporter, and the many-to-many mapping of substrates, transporters and genes. In this perspective, we summarize progress in both experimental and computational approaches used to determine the function of transporters and consider paths forward that integrate both. Investment in accurate, high-throughput functional characterization is needed to train the next-generation of predictive tools toward genome-scale metabolic network reconstructions that better predict phenotypes and interactions. More reliable predictions in this domain will benefit fields ranging from personalized medicine to metabolic engineering to microbial ecology.

Fabry disease (FD) is an X-linked lysosomal disease caused by an enzyme deficiency of alpha-galactosidase A (α-gal A). This deficiency leads to the accumulation of glycosphingolipids in lysosomes, resulting in a range of clinical symptoms. The complex pathogenesis of FD involves lysosomal dysfunction, altered autophagy, and mitochondrial abnormalities. Omics sciences, particularly transcriptomic analysis, comprehensively understand molecular mechanisms underlying diseases. This study focuses on genome-wide expression analysis in an FD human podocyte model to gain insights into the underlying mechanisms of podocyte dysfunction. Human control and GLA-edited podocytes were used. Gene expression data was generated using RNA-seq analysis, and differentially expressed genes were identified using DESeq2. Principal component analysis and Spearman correlation have explored gene expression trends. Functional enrichment and Reporter metabolite analyses were conducted to identify significantly affected metabolites and metabolic pathways. Differential expression analysis revealed 247 genes with altered expression levels in GLA-edited podocytes compared to control podocytes. Among these genes, 136 were underexpressed, and 111 were overexpressed in GLA-edited cells. Functional analysis of differentially expressed genes showed their involvement in various pathways related to oxidative stress, inflammation, fatty acid metabolism, collagen and extracellular matrix homeostasis, kidney injury, apoptosis, autophagy, and cellular stress response. The study provides insights into molecular mechanisms underlying Fabry podocyte dysfunction. Integrating transcriptomics data with genome-scale metabolic modeling further unveiled metabolic alterations in GLA-edited podocytes. This comprehensive approach contributes to a better understanding of Fabry disease and may lead to identifying new biomarkers and therapeutic targets for this rare lysosomal disorder.

ABSTRACT Pseudomonas aeruginosa is a ubiquitous, opportunistic human pathogen. Since it often expresses multidrug resistance, new treatment options are urgently required. Such new treatments are usually assessed with one of the canonical laboratory strains, PAO1 or PA14. However, these two strains are unlikely representative of the strains infecting patients, because they have adapted to laboratory conditions and do not capture the enormous genomic diversity of the species. Here, we characterized the major P. aeruginosa clone type (mPact) panel. This panel consists of 20 strains, which reflect the species’ genomic diversity, cover all major clone types, and have both patient and environmental origins. We found significant strain variation in distinct responses toward antibiotics and general growth characteristics. Only few of the measured traits are related, suggesting independent trait optimization across strains. High resistance levels were only identified for clinical mPact isolates and could be linked to known antimicrobial resistance (AMR) genes. One strain, H01, produced highly unstable AMR combined with reduced growth under drug-free conditions, indicating an evolutionary cost to resistance. The expression of microcolonies was common among strains, especially for strain H15, which also showed reduced growth, possibly indicating another type of evolutionary trade-off. By linking isolation source, growth, and virulence to life history traits, we further identified specific adaptive strategies for individual mPact strains toward either host processes or degradation pathways. Overall, the mPact panel provides a reasonably sized set of distinct strains, enabling in-depth analysis of new treatment designs or evolutionary dynamics in consideration of the species’ genomic diversity.IMPORTANCENew treatment strategies are urgently needed for high-risk pathogens such as the opportunistic and often multidrug-resistant pathogen Pseudomonas aeruginosa. Here, we characterize the major P. aeruginosa clone type (mPact) panel. It consists of 20 strains with different origins that cover the major clone types of the species as well as its genomic diversity. This mPact panel shows significant variation in (i) resistance against distinct antibiotics, including several last resort antibiotics; (ii) related traits associated with the response to antibiotics; and (iii) general growth characteristics. We further developed a novel approach that integrates information on resistance, growth, virulence, and life-history characteristics, allowing us to demonstrate the presence of distinct adaptive strategies of the strains that focus either on host interaction or resource processing. In conclusion, the mPact panel provides a manageable number of representative strains for this important pathogen for further in-depth analyses of treatment options and evolutionary dynamics.

ABSTRACT In a healthy colon, the stratified mucus layer serves as a crucial innate immune barrier to protect the epithelium from microbes. Mucins are complex glycoproteins that serve as a nutrient source for resident microflora and can be exploited by pathogens. We aimed to understand how the intestinal pathogen, Clostridioides difficile, independently uses or manipulates mucus to its benefit, without contributions from members of the microbiota. Using a 2-D primary human intestinal epithelial cell model to generate physiologic mucus, we assessed C. difficile-mucus interactions through growth assays, RNA-Seq, biophysical characterization of mucus, and contextualized metabolic modeling. We found that host-derived mucus promotes C. difficile growth both in vitro and in an infection model. RNA-Seq revealed significant upregulation of genes related to central metabolism in response to mucus, including genes involved in sugar uptake, the Wood-Ljungdahl pathway, and the glycine cleavage system. In addition, we identified differential expression of genes related to sensing and transcriptional control. Analysis of mutants with deletions in highly upregulated genes reflected the complexity of C. difficile-mucus interactions, with potential interplay between sensing and growth. Mucus also stimulated biofilm formation in vitro, which may in turn alter the viscoelastic properties of mucus. Context-specific metabolic modeling confirmed differential metabolism and the predicted importance of enzymes related to serine and glycine catabolism with mucus. Subsequent growth experiments supported these findings, indicating mucus is an important source of serine. Our results better define responses of C. difficile to human gastrointestinal mucus and highlight flexibility in metabolism that may influence pathogenesis.IMPORTANCEClostridioides difficile results in upward of 250,000 infections and 12,000 deaths annually in the United States. Community-acquired infections continue to rise, and recurrent disease is common, emphasizing a vital need to understand C. difficile pathogenesis. C. difficile undoubtedly interacts with colonic mucus, but the extent to which the pathogen can independently respond to and take advantage of this niche has not been explored extensively. Moreover, the metabolic complexity of C. difficile remains poorly understood but likely impacts its capacity to grow and persist in the host. Here, we demonstrate that C. difficile uses native colonic mucus for growth, indicating C. difficile possesses mechanisms to exploit the mucosal niche. Furthermore, mucus induces metabolic shifts and biofilm formation in C. difficile, which has potential ramifications for intestinal colonization. Overall, our work is crucial to better understand the dynamics of C. difficile-mucus interactions in the context of the human gut.

Blood metabolite levels are affected by numerous factors, including preanalytical factors such as collection methods and geographical sites. These perturbations have caused deleterious consequences for many metabolomics studies and represent a major challenge in the metabolomics field. It is important to understand these factors and develop models to reduce their perturbations. However, to date, the lack of suitable mathematical models for blood metabolite levels under homeostasis has hindered progress. In this study, we develop quantitative models of blood metabolite levels in healthy adults based on multisite sample cohorts that mimic the current challenge. Five cohorts of samples obtained across four geographically distinct sites were investigated, focusing on approximately 50 metabolites that were quantified using 1H NMR spectroscopy. More than one-third of the variation in these metabolite profiles is due to cross-cohort variation. A dramatic reduction in the variation of metabolite levels (90%), especially their site-to-site variation (95%), was achieved by modeling each metabolite using demographic and clinical factors and especially other metabolites, as observed in the top principal components. The results also reveal that several metabolites contribute disproportionately to such variation, which could be explained by their association with biological pathways including biosynthesis and degradation. The study demonstrates an intriguing network effect of metabolites that can be utilized to better define homeostatic metabolite levels, which may have implications for improved health monitoring. As an example of the potential utility of the approach, we show that modeling gender-related metabolic differences retains the interesting variance while reducing unwanted (site-related) variance. [ABSTRACT FROM AUTHOR]

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/BjaR1 QS system in Bradyrhizobium diazoe fficiens 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 ΔbjaR1 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 BjaR1 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, BjaR1 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/BjaR1 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. [ABSTRACT FROM AUTHOR]

Non-conventional yeasts are promising cell factories to produce lipids and oleochemicals, metabolites of industrial interest (e.g., organics acids, esters, and alcohols), and enzymes. They can also use different agro-industrial by-products as substrates within the context of a circular economy. Some of these yeasts can also comprise economic and health burdens as pathogens. Genome-scale metabolic models (GEMs), networks reconstructed based on the genomic and metabolic information of one or more organisms, are great tools to understand metabolic functions and landscapes, as well as propose engineering targets to improve metabolite production or propose novel drug targets. Previous reviews on yeast GEMs have mainly focused on the history and the evaluation of Saccharomyces cerevisiae modeling paradigms or the accessibility and usability of yeast GEMs. However, they did not describe the reconstruction strategies, limitations, validations, challenges, and research gaps of non-conventional yeast GEMs. Herein, we focused on the reconstruction of available non-Saccharomyces GEMs, their validation, underscoring the physiological insights, as well as the identification of both metabolic engineering and drug targets. We also discuss the challenges and knowledge gaps and propose strategies to boost their use and novel reconstructions. [ABSTRACT FROM AUTHOR]

Parkinson's disease (PD) is a chronic and progressive neurodegenerative disease with the major symptoms comprising loss of movement coordination (motor dysfunction) and non-motor dysfunction, including gastrointestinal symptoms. Alterations in the gut microbiota composition have been reported in PD patients vs. controls. However, it is still unclear how these compositional changes contribute to disease etiology and progression. Furthermore, most of the available studies have focused on European, Asian, and North American cohorts, but the microbiomes of PD patients in Latin America have not been characterized. To address this problem, we obtained fecal samples from Colombian participants (n = 25 controls, n = 25 PD idiopathic cases) to characterize the taxonomical community changes during disease via 16S rRNA gene sequencing. An analysis of differential composition, diversity, and personalized computational modeling was carried out, given the fecal bacterial composition and diet of each participant. We found three metabolites that differed in dietary habits between PD patients and controls: carbohydrates, trans fatty acids, and potassium. We identified six genera that changed significantly in their relative abundance between PD patients and controls, belonging to the families Lachnospiraceae, Lactobacillaceae, Verrucomicrobioaceae, Peptostreptococcaceae, and Streptococcaceae. Furthermore, personalized metabolic modeling of the gut microbiome revealed changes in the predicted production of seven metabolites (Indole, tryptophan, fructose, phenylacetic acid, myristic acid, 3-Methyl-2-oxovaleric acid, and N-Acetylneuraminic acid). These metabolites are associated with the metabolism of aromatic amino acids and their consumption in the diet. Therefore, this research suggests that each individual's diet and intestinal composition could affect host metabolism. Furthermore, these findings open the door to the study of microbiome–host interactions and allow us to contribute to personalized medicine. [ABSTRACT FROM AUTHOR]

Temperature plays a fundamental role in biology, influencing cellular function, chemical reaction rates, molecular structures, and interactions. While the temperature dependence of many biochemical reactions is well defined in vitro, the effect of temperature on metabolic function at the network level is poorly understood, and it remains an important challenge in optimizing the storage of cells and tissues at lower temperatures. Here, we used time-course metabolomic data and systems biology approaches to characterize the effects of storage temperature on human platelets (PLTs) in a platelet additive solution. We observed that changes to the metabolome with storage time do not simply scale with temperature but instead display complex temperature dependence, with only a small subset of metabolites following an Arrhenius-type relationship. Investigation of PLT energy metabolism through integration with computational modeling revealed that oxidative metabolism is more sensitive to temperature changes than glycolysis. The increased contribution of glycolysis to ATP turnover at lower temperatures indicates a stronger glycolytic phenotype with decreasing storage temperature. More broadly, these results demonstrate that the temperature dependence of the PLT metabolic network is not uniform, suggesting that efforts to improve the health of stored PLTs could be targeted at specific pathways. [ABSTRACT FROM AUTHOR]

SUMMARY: Potato (Solanum tuberosum) is a significant non‐grain food crop in terms of global production. However, its yield potential might be raised by identifying means to release bottlenecks within photosynthetic metabolism, from the capture of solar energy to the synthesis of carbohydrates. Recently, engineered increases in photosynthetic rates in other crops have been directly related to increased yield – how might such increases be achieved in potato? To answer this question, we derived the photosynthetic parameters Vcmax and Jmax to calibrate a kinetic model of leaf metabolism (e‐Photosynthesis) for potato. This model was then used to simulate the impact of manipulating the expression of genes and their protein products on carbon assimilation rates in silico through optimizing resource investment among 23 photosynthetic enzymes, predicting increases in photosynthetic CO2 uptake of up to 67%. However, this number of manipulations would not be practical with current technologies. Given a limited practical number of manipulations, the optimization indicated that an increase in amounts of three enzymes – Rubisco, FBP aldolase, and SBPase – would increase net assimilation. Increasing these alone to the levels predicted necessary for optimization increased photosynthetic rate by 28% in potato. [ABSTRACT FROM AUTHOR]

13C‐Metabolic Flux Analysis (13C‐MFA) and Flux Balance Analysis (FBA) are widely used to investigate the operation of biochemical networks in both biological and biotechnological research. Both methods use metabolic reaction network models of metabolism operating at steady state so that reaction rates (fluxes) and the levels of metabolic intermediates are constrained to be invariant. They provide estimated (MFA) or predicted (FBA) values of the fluxes through the network in vivo, which cannot be measured directly. These fluxes can shed light on basic biology and have been successfully used to inform metabolic engineering strategies. Several approaches have been taken to test the reliability of estimates and predictions from constraint‐based methods and to compare alternative model architectures. Despite advances in other areas of the statistical evaluation of metabolic models, such as the quantification of flux estimate uncertainty, validation and model selection methods have been underappreciated and underexplored. We review the history and state‐of‐the‐art in constraint‐based metabolic model validation and model selection. Applications and limitations of the χ2‐test of goodness‐of‐fit, the most widely used quantitative validation and selection approach in 13C‐MFA, are discussed, and complementary and alternative forms of validation and selection are proposed. A combined model validation and selection framework for 13C‐MFA incorporating metabolite pool size information that leverages new developments in the field is presented and advocated for. Finally, we discuss how adopting robust validation and selection procedures can enhance confidence in constraint‐based modeling as a whole and ultimately facilitate more widespread use of FBA in biotechnology. [ABSTRACT FROM AUTHOR]

ABSTRACTThe 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 (R2 = 0.92 for the training set and R2 = 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.

Genetic background and environmental conditions affect the production of sensory impact compounds by Saccharomyces cerevisiae. The relative importance of the strain-specific metabolic capabilities for the production of volatile organic compounds (VOCs) remains unclear. We investigated which amino acids contribute to VOC production and whether amino acid-VOC relations are conserved among yeast strains. Amino acid consumption and production of VOCs during grape juice fermentation was investigated using four commercial wine yeast strains: Elixir, Opale, R2, and Uvaferm. Principal component analysis of the VOC data demonstrated that Uvaferm correlated with ethyl acetate and ethyl hexanoate production, R2 negatively correlated with the acetate esters, and Opale positively correlated with fusel alcohols. Biomass formation was similar for all strains, pointing to metabolic differences in the utilization of nutrients to form VOCs. Partial least-squares linear regression showed that total aroma production is a function of nitrogen utilization (R2 = 0.87). We found that glycine, tyrosine, leucine, and lysine utilization were positively correlated with fusel alcohols and acetate esters. Mechanistic modeling of the yeast metabolic network via parsimonious flux balance analysis and flux enrichment analysis revealed enzymes with crucial roles, such as transaminases and decarboxylases. Our work provides insights in VOC production in wine yeasts. IMPORTANCE Saccharomyces cerevisiae is widely used in grape juice fermentation to produce wines. Along with the genetic background, the nitrogen in the environment in which S. cerevisiae grows impacts its regulation of metabolism. Also, commercial S. cerevisiae strains exhibit immense diversity in their formation of aromas, and a desirable aroma bouquet is an essential characteristic for wines. Since nitrogen affects aroma formation in wines, it is essential to know the extent of this connection and how it leads to strain-dependent aroma profiles in wines. We evaluated the differences in the production of key aroma compounds among four commercial wine strains. Moreover, we analyzed the role of nitrogen utilization on the formation of various aroma compounds. This work illustrates the unique aroma-producing differences among industrial yeast strains and suggests more intricate, nitrogen-associated routes influencing those aroma-producing differences.

ABSTRACT The Lyme disease bacterium Borrelia burgdorferi is extremely host dependent, with a small genome and correspondingly limited metabolism. As such, it is an excellent candidate for the development of targeted, narrow-spectrum antimicrobials. To accelerate drug discovery in this fastidious bacterium, in silico genome-scale metabolic modeling was used to construct a map of B. burgdorferi’s metabolism. This map was used to predict essential genes and enzymes; experimental data validated these predicted hits as viable drug targets. Repurposing existing small-molecule inhibitors, it is shown that inhibition of two predicted essential enzymes (pyridoxal kinase and serine hydroxymethyltransferase) selectively kills B. burgdorferi in culture. Thus, the essential processes identified here represent targets for the development of narrow-spectrum antimicrobials. This pipeline, pairing in silico discovery with validation in culture, may be useful for other genetically intractable pathogens.IMPORTANCELyme 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.

ABSTRACT Microbes typically reside in communities containing multiple species, whose interactions have considerable impacts on the robustness and functionality of such communities. To manage microbial communities, it is essential to understand the factors driving their assemblage and maintenance. Even though the community composition could be easily assessed, interspecies interactions during community establishment remain poorly understood. Here, we combined co-occurrence network analysis with quantitative PCR to examine the importance of each species within synthetic communities (SynComs) of pellicle biofilms. Genome-scale metabolic models and in vitro experiments indicated that the biomass of SynComs was primarily affected by keystone species that are acting either as metabolic facilitators or as competitors. Our study sets an example of how to construct a SynCom and investigate interspecies interactions.IMPORTANCECo-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.

The advent of the first synthetic cell has triggered significant interest in synthetic biology research and its potential impact on various fields, including industry, agriculture, health, and the environment. With its unique blend of molecular biology, genetics, engineering, and computational modeling, synthetic biology is becoming a crucial tool in addressing challenges in food science. For example, the development of synthetic genetic circuits is underway for the improvement of food safety and quality, such as creating biosensors for detecting toxins and heavy metals, enhancing the nutritional value of foods, and preserving the environment using sustainable technologies. Regardless of its potential, the use of synthetic biology in food production raises concerns regarding biosafety, asking for stringent control and comprehensive legislation to prevent misuse of the technology. This review provides an overview of synthetic biology in food science, discussing the current state-of-art of its toolbox achievements for cost-effective and nutritional improvement, likewise, the legal, ethical, and biosecurity considerations surrounding the use of synthetic circuits for the advancement of food science.

The modeling of rates of biochemical reactions—fluxes—in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint‐based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to‐date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C‐metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands‐on computer laboratory component of a four‐day metabolic modeling workshop and participant survey results showed improvements in learners' self‐assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level. [ABSTRACT FROM AUTHOR]

Many studies have demonstrated that the gut microbiota is associated with human health and disease. Manipulation of the gut microbiota, e.g. supplementation of probiotics, has been suggested to be feasible, but subject to limited therapeutic efficacy. To develop efficient microbiota-targeted diagnostic and therapeutic strategies, metabolic engineering has been applied to construct genetically modified probiotics and synthetic microbial consortia. This review mainly discusses commonly adopted strategies for metabolic engineering in the human gut microbiome, including the use of in silico , in vitro , or in vivo approaches for iterative design and construction of engineered probiotics or microbial consortia. Especially, we highlight how genome-scale metabolic models can be applied to advance our understanding of the gut microbiota. Also, we review the recent applications of metabolic engineering in gut microbiome studies as well as discuss important challenges and opportunities. • An adapted iterative cycle of Learn-Design-Build-Test (LDBT) for engineering the human gut microbiome was proposed. • How genome-scale metabolic models and community models can be applied to understand the gut microbiota was elaborated. • Recent progress of metabolic engineering in gut microbiome field was summarized. • Critical challenges and opportunities for metabolic engineering of gut microbiota were discussed. [ABSTRACT FROM AUTHOR]

ABSTRACT Haemophilus influenzae is a major cause of community-acquired pneumonia. While studied extensively in various laboratory models, less is known about the cell function while inside the human lung. We present the first analysis of the global gene expression of H. influenzae while the bacteria are in the lung during pneumonia (in vivo conditions) and contrast it with bacterial isolates that have been cultured under standard laboratory conditions (in vitro conditions). Patients with pneumonia were recruited from emergency departments and intensive care units during 2018–2020 (n = 102). Lower respiratory samples were collected for bacterial culture and RNA extraction. Patient samples with H. influenzae (n = 8) and colonies from bacterial cultures (n = 6) underwent RNA sequencing. The reads were then pseudo-aligned to core and pan genomes created from 15 reference strains. While bacteria cultured in vitro clustered tightly by principal component analysis of core genome (n = 1067) gene expression, bacteria in the patient samples had more diverse transcriptomic signatures and did not group with their lab-cultured counterparts. In total, 328 core genes were significantly differentially expressed between in vitro and in vivo conditions. The most highly upregulated genes in vivo included tbpA and fbpA, which are involved in the acquisition of iron from transferrin, and the stress response gene msrAB. The biosynthesis of nucleotides/purines and molybdopterin-scavenging processes were also significantly enriched in vivo. In contrast, major metabolic pathways and iron-sequestering genes were downregulated under this condition. In conclusion, extensive transcriptomic differences were found between bacteria while in the human lung and bacteria that were cultured in vitro. IMPORTANCE The human-specific pathogen Haemophilus influenzae is generally not well suited for studying in animal models, and most laboratory models are unlikely to approximate the diverse environments encountered by bacteria in the human airways accurately. Thus, we have examined the global gene expression of H. influenzae during pneumonia. Extensive differences in the global gene expression profiles were found in H. influenzae while in the human lung compared to bacteria that were grown in the laboratory. In contrast, the gene expression profiles of isolates collected from different patients were found to cluster together when grown under the same laboratory conditions. Interesting observations were made of how H. influenzae acquires and uses iron and molybdate, endures oxidative stress, and regulates central metabolism while in the lung. Our results indicate important processes during infection and can guide future research on genes and pathways that are relevant in the pathogenesis of H. influenzae pneumonia.

Fermentation employing Saccharomyces cerevisiae has produced alcoholic beverages and bread for millennia. More recently, S. cerevisiae has been used to manufacture specific metabolites for the food, pharmaceutical, and cosmetic industries. Among the most important of these metabolites are compounds associated with desirable aromas and flavors, including higher alcohols and esters. Although the physiology of yeast has been well‐studied, its metabolic modulation leading to aroma production in relevant industrial scenarios such as winemaking is still unclear. Here we ask what are the underlying metabolic mechanisms that explain the conserved and varying behavior of different yeasts regarding aroma formation under enological conditions? We employed dynamic flux balance analysis (dFBA) to answer this key question using the latest genome‐scale metabolic model (GEM) of S. cerevisiae. The model revealed several conserved mechanisms among wine yeasts, for example, acetate ester formation is dependent on intracellular metabolic acetyl‐CoA/CoA levels, and the formation of ethyl esters facilitates the removal of toxic fatty acids from cells using CoA. Species‐specific mechanisms were also found, such as a preference for the shikimate pathway leading to more 2‐phenylethanol production in the Opale strain as well as strain behavior varying notably during the carbohydrate accumulation phase and carbohydrate accumulation inducing redox restrictions during a later cell growth phase for strain Uvaferm. In conclusion, our new metabolic model of yeast under enological conditions revealed key metabolic mechanisms in wine yeasts, which will aid future research strategies to optimize their behavior in industrial settings. [ABSTRACT FROM AUTHOR]

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease. [ABSTRACT FROM AUTHOR]

Cyanobacteria are a promising platform for the production of the triterpene squalene (C30), a precursor for all plant and animal sterols, and a highly attractive intermediate towards triterpenoids, a large group of secondary plant metabolites. Synechocystis sp. PCC 6803 natively produces squalene from CO2 through the MEP pathway. Based on the predictions of a constraint-based metabolic model, we took a systematic overexpression approach to quantify native Synechocystis gene's impact on squalene production in a squalene-hopene cyclase gene knock-out strain (Δshc). Our in silico analysis revealed an increased flux through the Calvin-Benson-Bassham cycle in the Δshc mutant compared to the wildtype, including the pentose phosphate pathway, as well as lower glycolysis, while the tricarboxylic acid cycle predicted to be downregulated. Further, all enzymes of the MEP pathway and terpenoid synthesis, as well as enzymes from the central carbon metabolism, Gap2, Tpi and PyrK, were predicted to positively contribute to squalene production upon their overexpression. Each identified target gene was integrated into the genome of Synechocystis Δshc under the control of the rhamnose-inducible promoter Prha. Squalene production was increased in an inducer concentration dependent manner through the overexpression of most predicted genes, which are genes of the MEP pathway, ispH, ispE, and idi, leading to the greatest improvements. Moreover, we were able to overexpress the native squalene synthase gene (sqs) in Synechocystis Δshc, which reached the highest production titer of 13.72 mg l-1 reported for squalene in Synechocystis sp. PCC 6803 so far, thereby providing a promising and sustainable platform for triterpene production. [ABSTRACT FROM AUTHOR]

Genetic engineering of hyperthermophilic organisms for the production of fuels and other useful chemicals is an emerging biotechnological opportunity. In particular, for volatile organic compounds such as ethanol, fermentation at high temperatures could allow for straightforward separation by direct distillation. Currently, the upper growth temperature limit for native ethanol producers is 72°C in the bacterium Thermoanaerobacter ethanolicus JW200, and the highest temperature for heterologouslyengineered bioethanol production was recently demonstrated at 85°C in the archaeon Pyrococcus furiosus. Here, we describe an engineered strain of P. furiosus that synthesizes ethanol at 95°C, utilizing a homologously-expressed native alcohol dehydrogenase, termed AdhF. Ethanol biosynthesis was compared at 75°C and 95°C with various engineered strains. At lower temperatures, the acetaldehyde substrate for AdhF is most likely produced from acetate by aldehyde ferredoxin oxidoreductase (AOR). At higher temperatures, the effect of AOR on ethanol production is negligible, suggesting that acetaldehyde is produced by pyruvate ferredoxin oxidoreductase (POR) via oxidative decarboxylation of pyruvate, a reaction known to occur only at higher temperatures. Heterologous expression of a carbon monoxide dehydrogenase complex in the AdhF overexpression strain enabled it to use CO as a source of energy, leading to increased ethanol production. A genome reconstruction model for P. furiosus was developed to guide metabolic engineering strategies and understand outcomes. This work opens the door to the potential for 'bioreactive distillation' since fermentation can be performed well above the normal boiling point of ethanol. IMPORTANCE Previously, the highest temperature for biological ethanol production was 85°C. Here, we have engineered ethanol production at 95°C by the hyperthermophilic archaeon Pyrococcus furiosus. Using mutant strains, we showed that ethanol production occurs by different pathways at 75°C and 95°C. In addition, by heterologous expression of a carbon monoxide dehydrogenase complex, ethanol production by this organism was driven by the oxidation of carbon monoxide. A genome reconstruction model for P. furiosus was developed to guide metabolic engineering strategies and understand outcomes. [ABSTRACT FROM AUTHOR]

All plant and animal kingdom organisms use highly connected biochemical networks to facilitate sustaining, proliferation, and growth functions. While the biochemical network details are well known, the understanding of the intense regulation principles is still limited. We chose to investigate the Hermetia illucens fly at the larval stage because this stage is a crucial period for the successful accumulation and allocation of resources for the subsequent organism's developmental stages. We combined iterative wet lab experiments and innovative metabolic modeling design approaches to simulate and explain the H. illucens larval stage resource allocation processes and biotechnology potential. We performed time-based growth and high-value chemical compound accumulation wet lab chemical analysis experiments on larvae and the Gainesville diet composition. We built and validated the first H. illucens medium-size, stoichiometric metabolic model to predict the effects of diet-based alterations on fatty acid allocation potential. Using optimization methods such as flux balance and flux variability analysis on the novel insect metabolic model, we predicted that doubled essential amino acid consumption increased the growth rate by 32%, but pure glucose consumption had no positive impact on growth. In the case of doubled pure valine consumption, the model predicted a 2% higher growth rate. In this study, we describe a new framework for researching the impact of dietary alterations on the metabolism of multi-cellular organisms at different developmental stages for improved, sustainable, and directed high-value chemicals. [ABSTRACT FROM AUTHOR]

Clarifying the general rules behind microbial community assembly will foster the development of microbiome-based technological solutions. Here, we study microbial community assembly through a computational analysis of phylogenetic core groups (PCGs): discrete portions of the bacterial phylogeny with high prevalence in the ecosystem under study. We first show that the existence of PCGs was a predominant feature of the varied set of microbial ecosystems studied. Then, we re-analyzed an in vitro experimental dataset using a PCG-based approach, drawing only from its community composition data and from publicly available genomic databases. Using mainly genome scale metabolic models and population dynamics modeling, we obtained ecological insights on metabolic niche structure and population dynamics comparable to those gained after canonical experimentation. Thus, leveraging phylogenetic signal to help unravel microbiome function and assembly rules offers a potential avenue to gain further insight on Earth’s microbial ecosystems.

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.

The urgent necessity to fight antimicrobial resistance is universally recognized. In the search of new targets and strategies to face this global challenge, a promising approach resides in the study of the cellular response to antimicrobial exposure and on the impact of global cellular reprogramming on antimicrobial drugs' efficacy. The metabolic state of microbial cells has been shown to undergo several antimicrobial-induced modifications and, at the same time, to be a good predictor of the outcome of an antimicrobial treatment. Metabolism is a promising reservoir of potential drug targets/adjuvants that has not been fully exploited to date. One of the main problems in unraveling the metabolic response of cells to the environment resides in the complexity of such metabolic networks. To solve this problem, modeling approaches have been developed, and they are progressively gaining in popularity due to the huge availability of genomic information and the ease at which a genome sequence can be converted into models to run basic phenotype predictions. Here, we review the use of computational modeling to study the relationship between microbial metabolism and antimicrobials and the recent advances in the application of genome-scale metabolic modeling to the study of microbial responses to antimicrobial exposure. [ABSTRACT FROM AUTHOR]