Search

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding betweenmicrobiologists andmodelers/bioinformaticians, stemming froma lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored formicrobiologists, researchers newtomicrobiomemodeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling. [ABSTRACT FROM AUTHOR]

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed. [ABSTRACT FROM AUTHOR]

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

ABSTRACTGenome-scale metabolic modeling is a powerful framework for predicting metabolic phenotypes of any organism with an annotated genome. For two decades, this framework has been used for the rational design of microbial cell factories. In the last decade, the range of applications has exploded, and new frontiers have emerged, including the study of the gut microbiome and its health implications and the role of microbial communities in global ecosystems. However, all the critical steps in this framework, from model construction to simulation, require the use of powerful linear optimization solvers, with the choice often relying on commercial solvers for their well-known computational efficiency. In this work, I benchmark a total of six solvers (two commercial and four open source) and measure their performance to solve linear and mixed-integer linear problems of increasing complexity. Although commercial solvers are still the fastest, at least two open-source solvers show comparable performance. These results show that genome-scale metabolic modeling does not need to be hindered by commercial licensing schemes and can become a truly open science framework for solving urgent societal challenges.IMPORTANCEModeling the metabolism of organisms and communities allows for computational exploration of their metabolic capabilities and testing their response to genetic and environmental perturbations. This holds the potential to address multiple societal issues related to human health and the environment. One of the current limitations is the use of commercial optimization solvers with restrictive licenses for academic and non-academic use. This work compares the performance of several commercial and open-source solvers to solve some of the most complex problems in the field. Benchmarking results show that, although commercial solvers are indeed faster, some of the open-source options can also efficiently tackle the hardest problems, showing great promise for the development of open science applications.

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed.

A genome-scale metabolic model, encompassing a total of 623 genes, 727 reactions, and 865 metabolites, was developed for Pyrococcus furiosus, an archaeon that grows optimally at 100°C by carbohydrate and peptide fermentation. The model uses subsystem-based genome annotation, along with extensive manual curation of 237 genereaction associations including those involved in central carbon metabolism, amino acid metabolism, and energy metabolism. The redox and energy balance of P. furiosus was investigated through random sampling of flux distributions in the model during growth on disaccharides. The core energy balance of the model was shown to depend on high acetate production and the coupling of a sodium-dependent ATP synthase and membrane-bound hydrogenase, which generates a sodium gradient in a ferredoxin-dependent manner, aligning with existing understanding of P. furiosus metabolism. The model was utilized to inform genetic engineering designs that favor the production of ethanol over acetate by implementing an NADPH and CO-dependent energy economy. The P. furiosus model is a powerful tool for understanding the relationship between generation of end products and redox/energy balance at a systems-level that will aid in the design of optimal engineering strategies for production of bio-based chemicals and fuels. IMPORTANCE The bio-based production of organic chemicals provides a sustainable alternative to fossil-based production in the face of today's climate challenges. In this work, we present a genome-scale metabolic reconstruction of Pyrococcus furiosus, a wellestablished platform organism that has been engineered to produce a variety of chemicals and fuels. The metabolic model was used to design optimal engineering strategies to produce ethanol. The redox and energy balance of P. furiosus was examined in detail, which provided useful insights that will guide future engineering designs. [ABSTRACT FROM AUTHOR]

Evidence suggests that novel enzyme functions evolved from low-level promiscuous activities in ancestral enzymes. Yet, the evolutionary dynamics and physiological mechanisms of how such side activities contribute to systems-level adaptations are not well characterized. Furthermore, it remains untested whether knowledge of an organism's promiscuous reaction set, or underground metabolism, can aid in forecasting the genetic basis of metabolic adaptations. Here, we employ a computational model of underground metabolism and laboratory evolution experiments to examine the role of enzyme promiscuity in the acquisition and optimization of growth on predicted non-native substrates in Escherichia coli K-12 MG1655. After as few as approximately 20 generations, evolved populations repeatedly acquired the capacity to grow on five predicted non-native substrates-D-lyxose, D-2-deoxyribose, D-arabinose, m-tartrate, and monomethyl succinate. Altered promiscuous activities were shown to be directly involved in establishing high-efficiency pathways. Structural mutations shifted enzyme substrate turnover rates toward the new substrate while retaining a preference for the primary substrate. Finally, genes underlying the phenotypic innovations were accurately predicted by genome-scale model simulations of metabolism with enzyme promiscuity.

There is great interest in engineering photoautotrophic metabolism to generate bioproducts of societal importance. Despite the success in employing genome-scale modeling coupled with flux balance analysis to engineer heterotrophic metabolism, the lack of proper constraints necessary to generate biologically realistic predictions has hindered broad application of this methodology to phototrophic metabolism. Here we describe a methodology for constraining genome-scale models of photoautotrophy in the cyanobacteria Synechococcus elongatus PCC 7942. Experimental photophysiology parameters coupled to genome-scale flux balance analysis resulted in accurate predictions of growth rates and metabolic reaction fluxes at low and high light conditions. Additionally, by constraining photon uptake fluxes, we characterized the metabolic cost of excess excitation energy. The predicted energy fluxes were consistent with known light-adapted phenotypes in cyanobacteria. Finally, we leveraged the modeling framework to characterize existing photoautotrophic and photomixtotrophic engineering strategies for 2,3-butanediol production in S. elongatus. This methodology, applicable to genome-scale modeling of all phototrophic microorganisms, can facilitate the use of flux balance analysis in the engineering of light-driven metabolism.

The production of monoclonal antibodies (mAbs) using Chinese Hamster Ovary (CHO) cells has revolutionized the treatment of numerous diseases, solidifying their position as a cornerstone of the biopharmaceutical industry. However, achieving maximum mAb production while upholding strict product quality standards remains a significant hurdle. Optimizing cell culture media emerges as a critical factor in this endeavor, requiring a nuanced understanding of the complex interplay of nutrients, growth factors, and other components that profoundly influence cellular growth, productivity, and product quality. Significant strides have been made in media optimization, including techniques such as media blending, one factor at a time, and statistical design of experiments approaches. The present review provides a comprehensive analysis of the recent advancements in culture media design strategies, focusing on the comparative application of systems biology (SB) and machine learning (ML) approaches. The applications of SB and ML in optimizing CHO cell culture medium and successful examples of their use are summarized. Finally, we highlight the immense potential of integrating SB and ML, emphasizing the development of hybrid models that leverage the strengths of both approaches for robust, efficient, and scalable optimization of mAb production in CHO cells. This review provides a roadmap for researchers and industry professionals to navigate the complex landscape of mAb production optimization, paving the way for developing next-generation CHO cell culture media that drive significant improvements in yield and productivity. [ABSTRACT FROM AUTHOR]

ABSTRACT Microbial communities experience continuous environmental changes, with temperature fluctuations being the most impacting. This is particularly important considering the ongoing global warming but also in the “simpler” context of seasonal variability of sea-surface temperature. Understanding how microorganisms react at the cellular level can improve our understanding of their possible adaptations to a changing environment. In this work, we investigated the mechanisms through which metabolic homeostasis is maintained in a cold-adapted marine bacterium during growth at temperatures that differ widely (15 and 0°C). We have quantified its intracellular and extracellular central metabolomes together with changes occurring at the transcriptomic level in the same growth conditions. This information was then used to contextualize a genome-scale metabolic reconstruction, and to provide a systemic understanding of cellular adaptation to growth at 2 different temperatures. Our findings indicate a strong metabolic robustness at the level of the main central metabolites, counteracted by a relatively deep transcriptomic reprogramming that includes changes in gene expression of hundreds of metabolic genes. We interpret this as a transcriptomic buffering of cellular metabolism, able to produce overlapping metabolic phenotypes, despite the wide temperature gap. Moreover, we show that metabolic adaptation seems to be mostly played at the level of few key intermediates (e.g., phosphoenolpyruvate) and in the cross talk between the main central metabolic pathways. Overall, our findings reveal a complex interplay at gene expression level that contributes to the robustness/resilience of core metabolism, also promoting the leveraging of state-of-the-art multi-disciplinary approaches to fully comprehend molecular adaptations to environmental fluctuations. IMPORTANCE This manuscript addresses a central and broad interest topic in environmental microbiology, i.e. the effect of growth temperature on microbial cell physiology. We investigated if and how metabolic homeostasis is maintained in a cold-adapted bacterium during growth at temperatures that differ widely and that match measured changes on the field. Our integrative approach revealed an extraordinary robustness of the central metabolome to growth temperature. However, this was counteracted by deep changes at the transcriptional level, and especially in the metabolic part of the transcriptome. This conflictual scenario was interpreted as a transcriptomic buffering of cellular metabolism, and was investigated using genome-scale metabolic modeling. Overall, our findings reveal a complex interplay at gene expression level that contributes to the robustness/resilience of core metabolism, also promoting the use of state-of-the-art multi-disciplinary approaches to fully comprehend molecular adaptations to environmental fluctuations.

SUMMARY: Characterizing photosynthetic productivity is necessary to understand the ecological contributions and biotechnology potential of plants, algae, and cyanobacteria. Light capture efficiency and photophysiology have long been characterized by measurements of chlorophyll fluorescence dynamics. However, these investigations typically do not consider the metabolic network downstream of light harvesting. By contrast, genome‐scale metabolic models capture species‐specific metabolic capabilities but have yet to incorporate the rapid regulation of the light harvesting apparatus. Here, we combine chlorophyll fluorescence parameters defining photosynthetic and non‐photosynthetic yield of absorbed light energy with a metabolic model of the pennate diatom Phaeodactylum tricornutum. This integration increases the model predictive accuracy regarding growth rate, intracellular oxygen production and consumption, and metabolic pathway usage. Through the quantification of excess electron transport, we uncover the sequential activation of non‐radiative energy dissipation processes, cross‐compartment electron shuttling, and non‐photochemical quenching as the rapid photoacclimation strategy in P. tricornutum. Interestingly, the photon absorption thresholds that trigger the transition between these mechanisms were consistent at low and high incident photon fluxes. We use this understanding to explore engineering strategies for rerouting cellular resources and excess light energy towards bioproducts in silico. Overall, we present a methodology for incorporating a common, informative data type into computational models of light‐driven metabolism and show its utilization within the design–build–test–learn cycle for engineering of photosynthetic organisms. Significance Statement: Whole genome models of metabolism are improved by the addition of parameters associated with measured photo‐physiology. The approach is demonstrated in the photosynthetic alga, Phaeodactylum tricornutum, but this approach could be applied to any plant system. [ABSTRACT FROM AUTHOR]

Flux balance analysis (FBA) is currently the standard method to compute metabolic fluxes in genome-scale networks. Several FBA extensions employing diverse objective functions and/or constraints have been published. Here we propose a hybrid semi-parametric FBA extension that combines mechanistic-level constraints (parametric) with empirical constraints (non-parametric) in the same linear program. A CHO dataset with 27 measured exchange fluxes obtained from 21 reactor experiments served to evaluate the method. The mechanistic constraints were deduced from a reduced CHO-K1 genome-scale network with 686 metabolites, 788 reactions and 210 degrees of freedom. The non-parametric constraints were obtained by principal component analysis of the flux dataset. The two types of constraints were integrated in the same linear program showing comparable computational cost to standard FBA. The hybrid FBA is shown to significantly improve the specific growth rate prediction under different constraints scenarios. A metabolically efficient cell growth feed targeting minimal byproducts accumulation was designed by hybrid FBA. It is concluded that integrating parametric and nonparametric constraints in the same linear program may be an efficient approach to reduce the solution space and to improve the predictive power of FBA methods when critical mechanistic information is missing. [ABSTRACT FROM AUTHOR]

Heme is an oxygen carrier and a cofactor of both industrial enzymes and food additives. The intracellular level of free heme is low, which limits the synthesis of heme proteins. Therefore, increasing heme synthesis allows an increased production of heme proteins. Using the genome-scale metabolic model (GEM) Yeast8 for the yeast Saccharomyces cerevisiae, we identified fluxes potentially important to heme synthesis. With this model, in silico simulations highlighted 84 gene targets for balancing biomass and increasing heme production. Of those identified, 76 genes were individually deleted or overexpressed in experiments. Empirically, 40 genes individually increased heme production (up to threefold). Heme was increased by modifying target genes, which not only included the genes involved in heme biosynthesis, but also those involved in glycolysis, pyruvate, Fe-S clusters, glycine, and succinyl-coenzyme A (CoA) metabolism. Next, we developed an algorithmic method for predicting an optimal combination of these genes by using the enzyme-constrained extension of the Yeast8 model, ecYeast8. The computationally identified combination for enhanced heme production was evaluated using the heme ligand-binding biosensor (Heme-LBB). The positive targets were combined using CRISPR-Cas9 in the yeast strain (IMX581-HEM15-HEM14-HEM3-Δshm1-HEM2-Δhmx1-FET4-Δgcv2-HEM1-Δgcv1-HEM13), which produces 70-fold-higher levels of intracellular heme. [ABSTRACT FROM AUTHOR]

Escherichia coli strains are widely used in academic research and biotechnology. New technologies for quantifying strain-specific differences and their underlying contributing factors promise greater understanding of how these differences significantly impact physiology, synthetic biology, metabolic engineering, and process design. Here, we quantified strain-specific differences in seven widely used strains of E. coli (BL21, C, Crooks, DH5a, K-12 MG1655, K-12 W3110, and W) using genomics, phenomics, transcriptomics, and genome-scale modeling. Metabolic physiology and gene expression varied widely with downstream implications for productivity, product yield, and titer. These differences could be linked to differential regulatory structure. Analyzing high-flux reactions and expression of encoding genes resulted in a correlated and quantitative link between these sets, with strain-specific caveats. Integrated modeling revealed that certain strains are better suited to produce given compounds or express desired constructs considering native expression states of pathways that enable high-production phenotypes. This study yields a framework for quantitatively comparing strains in a species with implications for strain selection.

The availability of multi-omics data sets and genome-scale metabolic models for various organisms provide a platform for modeling and analyzing genotype-to-phenotype relationships. Flux balance analysis is the main tool for predicting flux distributions in genome-scale metabolic models and various data-integrative approaches enable modeling context-specific network behavior. Due to its linear nature, this optimization framework is readily scalable to multi-tissue or -organ and even multi-organism models. However, both data and model size can hamper a straightforward biological interpretation of the estimated fluxes. Moreover, flux balance analysis simulates metabolism at steady-state and thus, in its most basic form, does not consider kinetics or regulatory events. The integration of flux balance analysis with complementary data analysis and modeling techniques offers the potential to overcome these challenges. In particular machine learning approaches have emerged as the tool of choice for data reduction and selection of most important variables in big data sets. Kinetic models and formal languages can be used to simulate dynamic behavior. This review article provides an overview of integrative studies that combine flux balance analysis with machine learning approaches, kinetic models, such as physiology-based pharmacokinetic models, and formal graphical modeling languages, such as Petri nets. We discuss the mathematical aspects and biological applications of these integrated approaches and outline challenges and future perspectives.

The human upper respiratory tract is the reservoir of a diverse community of commensals and potential pathogens (pathobionts), including Streptococcus pneumoniae (pneumococcus), Haemophilus influenzae, Moraxella catarrhalis, and Staphylococcus aureus, which occasionally turn into pathogens causing infectious diseases, while the contribution of many nasal microorganisms to human health remains undiscovered. To better understand the composition of the nasal microbiome community, we create a workflow of the community model, which mimics the human nasal environment. To address this challenge, constraint-based reconstruction of biochemically accurate genome-scale metabolic models (GEMs) networks of microorganisms is mandatory. Our workflow applies constraint-based modeling (CBM), simulates the metabolism between species in a given microbiome, and facilitates generating novel hypotheses on microbial interactions. Utilizing this workflow, we hope to gain a better understanding of interactions from the metabolic modeling perspective. This article presents nasal community modeling workflow (NCMW)—a python package based on GEMs of species as a starting point for understanding the composition of the nasal microbiome community. The package is constructed as a step-by-step mathematical framework for metabolic modeling and analysis of the nasal microbial community. Using constraint-based models reduces the need for culturing species in vitro, a process that is not convenient in the environment of human noses.Availability: NCMW is freely available on the Python Package Index (PIP) via pip install NCMW. The source code, documentation, and usage examples (Jupyter Notebook and example files) are available at https://github.com/manuelgloeckler/ncmw.

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. The advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

Enzyme promiscuity toward substrates has been discussed in evolutionary terms as providing the flexibility to adapt to novel environments. In the present work, we describe an approach toward exploring such enzyme promiscuity in the space of a metabolic network. This approach leverages genome-scale models, which have been widely used for predicting growth phenotypes in various environments or following a genetic perturbation; however, these predictions occasionally fail. Failed predictions of gene essentiality offer an opportunity for targeting biological discovery, suggesting the presence of unknown underground pathways stemming from enzymatic cross-reactivity. We demonstrate a workflow that couples constraint-based modeling and bioinformatic tools with KO strain analysis and adaptive laboratory evolution for the purpose of predicting promiscuity at the genome scale. Three cases of genes that are incorrectly predicted as essential in Escherichia coli--aspC, argD, and gltA--are examined, and isozyme functions are uncovered for each to a different extent. Seven isozyme functions based on genetic and transcriptional evidence are suggested between the genes aspC and tyrB, argD and astC, gabT and puuE, and gltA and prpC. This study demonstrates how a targeted model-driven approach to discovery can systematically fill knowledge gaps, characterize underground metabolism, and elucidate regulatory mechanisms of adaptation in response to gene KO perturbations.

Cordyceps militaris is an entomopathogenic fungus which is often used in Asia as a traditional medicine developed from age-old wisdom. Presently, cordycepin from C. militaris is a great interest in medicinal applications. However, cellular growth of C. militaris and the association with cordycepin production remain poorly understood. To explore the metabolism of C. militaris as potential cell factories in medical and biotechnology applications, this study developed a high-quality genome-scale metabolic model of C. militaris, iNR1329, based on its genomic content and physiological data. The model included a total of 1329 genes, 1821 biochemical reactions, and 1171 metabolites among 4 different cellular compartments. Its in silico growth simulation results agreed well with experimental data on different carbon sources. iNR1329 was further used for optimizing the growth and cordycepin overproduction using a novel approach, POPCORN, for rational design of synthetic media. In addition to the high-quality GEM iNR1329, the presented POPCORN approach was successfully used to rationally design an optimal synthetic medium with C:N ratio of 8:1 for enhancing 3.5-fold increase in cordycepin production. This study thus provides a novel insight into C. militaris physiology and highlights a potential GEM-driven method for synthetic media design and metabolic engineering application. The iNR1329 and the POPCORN approach are available at the GitHub repository: https://github.com/sysbiomics/Cordyceps_militaris-GEM.

Abstract Background Yarrowia lipolytica, an oleaginous yeast, is a promising platform strain for production of biofuels and oleochemicals as it can accumulate a high level of lipids in response to nitrogen limitation. Accordingly, many metabolic engineering efforts have been made to develop engineered strains of Y. lipolytica with higher lipid yields. Genome-scale model of metabolism (GEM) is a powerful tool for identifying novel genetic designs for metabolic engineering. Several GEMs for Y. lipolytica have recently been developed; however, not many applications of the GEMs have been reported for actual metabolic engineering of Y. lipolytica. The major obstacle impeding the application of Y. lipolytica GEMs is the lack of proper methods for predicting phenotypes of the cells in the nitrogen-limited condition, or more specifically in the stationary phase of a batch culture. Results In this study, we showed that environmental version of minimization of metabolic adjustment (eMOMA) can be used for predicting metabolic flux distribution of Y. lipolytica under the nitrogen-limited condition and identifying metabolic engineering strategies to improve lipid production in Y. lipolytica. Several well-characterized overexpression targets, such as diglyceride acyltransferase, acetyl-CoA carboxylase, and stearoyl-CoA desaturase, were successfully rediscovered by our eMOMA-based design method, showing the relevance of prediction results. Interestingly, the eMOMA-based design method also suggested non-intuitive knockout targets, and we experimentally validated the prediction with a mutant lacking YALI0F30745g, one of the predicted targets involved in one-carbon/methionine metabolism. The mutant accumulated 45% more lipids compared to the wild-type. Conclusion This study demonstrated that eMOMA is a powerful computational method for understanding and engineering the metabolism of Y. lipolytica and potentially other oleaginous microorganisms.

Abstract Background The maturation of the female germ cell, the oocyte, requires the synthesis and storing of all the necessary metabolites to support multiple divisions after fertilization. Oocyte maturation is only possible in the presence of surrounding, diverse, and changing layers of somatic cells. Our understanding of metabolic interactions between the oocyte and somatic cells has been limited due to dynamic nature of ovarian follicle development, thus warranting a systems approach. Results Here, we developed a genome-scale metabolic model of the mouse ovarian follicle. This model was constructed using an updated mouse general metabolic model (Mouse Recon 2) and contains several key ovarian follicle development metabolic pathways. We used this model to characterize the changes in the metabolism of each follicular cell type (i.e., oocyte, granulosa cells, including cumulus and mural cells), during ovarian follicle development in vivo. Using this model, we predicted major metabolic pathways that are differentially active across multiple follicle stages. We identified a set of possible secreted and consumed metabolites that could potentially serve as biomarkers for monitoring follicle development, as well as metabolites for addition to in vitro culture media that support the growth and maturation of primordial follicles. Conclusions Our systems approach to model follicle metabolism can guide future experimental studies to validate the model results and improve oocyte maturation approaches and support growth of primordial follicles in vitro.

Abstract Evidence suggests that novel enzyme functions evolved from low‐level promiscuous activities in ancestral enzymes. Yet, the evolutionary dynamics and physiological mechanisms of how such side activities contribute to systems‐level adaptations are not well characterized. Furthermore, it remains untested whether knowledge of an organism's promiscuous reaction set, or underground metabolism, can aid in forecasting the genetic basis of metabolic adaptations. Here, we employ a computational model of underground metabolism and laboratory evolution experiments to examine the role of enzyme promiscuity in the acquisition and optimization of growth on predicted non‐native substrates in Escherichia coli K‐12 MG1655. After as few as approximately 20 generations, evolved populations repeatedly acquired the capacity to grow on five predicted non‐native substrates—D‐lyxose, D‐2‐deoxyribose, D‐arabinose, m‐tartrate, and monomethyl succinate. Altered promiscuous activities were shown to be directly involved in establishing high‐efficiency pathways. Structural mutations shifted enzyme substrate turnover rates toward the new substrate while retaining a preference for the primary substrate. Finally, genes underlying the phenotypic innovations were accurately predicted by genome‐scale model simulations of metabolism with enzyme promiscuity.

The gut microbiome has emerged as a contributing factor in non-communicable disease, rendering it a target of health-promoting interventions. Yet current understanding of the host-microbiome dynamic is insufficient to predict the variation in intervention outcomes across individuals. We explore the mechanisms that underpin the gut bacterial ecosystem and highlight how a more complete understanding of this ecology will enable improved intervention outcomes. This ecology varies within the gut over space and time. Interventions disrupt these processes, with cascading consequences throughout the ecosystem. In vivo studies cannot isolate and probe these processes at the required spatiotemporal resolutions, and in vitro studies lack the representative complexity required. However, we highlight that, together, both approaches can inform in silico models that integrate cellular-level dynamics, can extrapolate to explain bacterial community outcomes, permit experimentation and observation over ecological processes at high spatiotemporal resolution, and can serve as predictive platforms on which to prototype interventions. Thus, it is a concerted integration of these techniques that will enable rational targeted manipulations of the gut ecosystem.

Many cancer cells consume glutamine at high rates; counterintuitively, they simultaneously excrete glutamate, the first intermediate in glutamine metabolism. Glutamine consumption has been linked to replenishment of tricarboxylic acid cycle (TCA) intermediates and synthesis of adenosine triphosphate (ATP), but the reason for glutamate excretion is unclear. Here, we dynamically profile the uptake and excretion fluxes of a liver cancer cell line (HepG2) and use genome-scale metabolic modeling for in-depth analysis. We find that up to 30% of the glutamine is metabolized in the cytosol, primarily for nucleotide synthesis, producing cytosolic glutamate. We hypothesize that excreting glutamate helps the cell to increase the nucleotide synthesis rate to sustain growth. Indeed, we show experimentally that partial inhibition of glutamate excretion reduces cell growth. Our integrative approach thus links glutamine addiction to glutamate excretion in cancer and points toward potential drug targets. [ABSTRACT FROM AUTHOR]

Living cells can express different metabolic pathways that support growth. The criteria that determine which pathways are selected in which environment remain unclear. One recurrent selection is overflow metabolism: the simultaneous usage of an ATP-efficient and -inefficient pathway, shown for example in Escherichia coli, Saccharomyces cerevisiae and cancer cells. Many models, based on different assumptions, can reproduce this observation. Therefore, they provide no conclusive evidence which mechanism is causing overflow metabolism. We compare the mathematical structure of these models. Although ranging from flux balance analyses to self-fabricating metabolism and expression models, we can rewrite all models into one standard form. We conclude that all models predict overflow metabolism when two, model-specific, growth-limiting constraints are hit. This is consistent with recent theory. Thus, identifying these two constraints is essential for understanding overflow metabolism. We list all imposed constraints by these models, so that they can hopefully be tested in future experiments. [ABSTRACT FROM AUTHOR]

ABSTRACT Quantitative understanding and prediction of microbial community dynamics are an outstanding challenge. We test the hypothesis that metabolic mechanisms provide a foundation for accurate prediction of dynamics in microbial systems. In our research, metabolic models have been able to accurately predict species interactions, evolutionary trajectories, and response to perturbation in simple synthetic consortia. However, metabolic models have many constraints and often serve best as null models to identify additional processes at play. We anticipate that major advances in metabolic systems biology will involve scaling bottom-up approaches to complex communities and expanding the processes that are incorporated in a metabolic perspective. Ultimately, cellular metabolism will inform predictive ecology that enables precision management of microbial systems.

The emergence and spread of metallo-beta-lactamase–producing multidrug-resistant (MDR) Klebsiella pneumoniae is a serious public health threat, which is further complicated by the increased prevalence of colistin resistance. The link between antimicrobial resistance acquired by strains of Klebsiella and their unique metabolic capabilities has not been determined. Here, we reconstruct genome-scale metabolic models for 22 K. pneumoniae strains with various resistance profiles to different antibiotics, including two strains exhibiting colistin resistance isolated from Cairo, Egypt. We use the models to predict growth capabilities on 265 different sole carbon, nitrogen, sulfur, and phosphorus sources for all 22 strains. Alternate nitrogen source utilization of glutamate, arginine, histidine, and ethanolamine among others provided discriminatory power for identifying resistance to amikacin, tetracycline, and gentamicin. Thus, genome-scale model based predictions of growth capabilities on alternative substrates may lead to construction of classification trees that are indicative of antibiotic resistance in Klebsiella isolates.

Genome-scale metabolic modeling is a powerful framework for predicting metabolic phenotypes of any organism with an annotated genome. For two decades, this framework has been used for the rational design of microbial cell factories. In the last decade, the range of applications has exploded, and new frontiers have emerged, including the study of the gut microbiome and its health implications and the role of microbial communities in global ecosystems. However, all the critical steps in this framework, from model construction to simulation, require the use of powerful linear optimization solvers, with the choice often relying on commercial solvers for their well-known computational efficiency. In this work, I benchmark a total of six solvers (two commercial and four open source) and measure their performance to solve linear and mixed-integer linear problems of increasing complexity. Although commercial solvers are still the fastest, at least two open-source solvers show comparable performance. These results show that genome-scale metabolic modeling does not need to be hindered by commercial licensing schemes and can become a truly open science framework for solving urgent societal challenges.IMPORTANCEModeling the metabolism of organisms and communities allows for computational exploration of their metabolic capabilities and testing their response to genetic and environmental perturbations. This holds the potential to address multiple societal issues related to human health and the environment. One of the current limitations is the use of commercial optimization solvers with restrictive licenses for academic and non-academic use. This work compares the performance of several commercial and open-source solvers to solve some of the most complex problems in the field. Benchmarking results show that, although commercial solvers are indeed faster, some of the open-source options can also efficiently tackle the hardest problems, showing great promise for the development of open science applications., Competing Interests: The author declares no conflict of interest.

Metabolic reprogramming is a hallmark of cancer. The main aim of this paper is to integrate a genome-scale metabolic description of tumor cells into a tumor growth model that accounts for the spatiotemporally heterogeneous tumor microenvironment, in order to study the effects of microscopic characteristics on tumor evolution. A lactate maximization metabolic strategy that allows near-optimal growth solution, while maximizing lactate secretion, is assumed. The proposed sub-cellular metabolic model is then incorporated into a hybrid discrete-continuous model of tumor growth. We produced several phenotypes by applying different constraints and optimization criteria in the metabolic model and explored the tumor evolution of the various phenotypes in different vasculature conditions and extracellular matrix densities. At first, we showed that the metabolic capabilities of phenotypes depending on resource availability can vary in a counter-intuitive manner. We then showed that: first, tumor population, morphology, and spread are affected differently in different conditions, allowing thus phenotypes to be superior than others in different conditions; and second, polyclonal tumors consisting of different phenotypes can exploit their different metabolic capabilities to enhance further tumor evolution. The proposed framework comprises a proof-of-concept demonstration showing the importance of considering the metabolic capabilities of phenotypes on predicting tumor evolution. The proposed framework allows the incorporation of context-specific and patient-specific data for the study of personalized tumor evolution and therapy efficacy, linking genome to metabolic capabilities and tumor dynamics. [ABSTRACT FROM AUTHOR]

Background: The maturation of the female germ cell, the oocyte, requires the synthesis and storing of all the necessary metabolites to support multiple divisions after fertilization. Oocyte maturation is only possible in the presence of surrounding, diverse, and changing layers of somatic cells. Our understanding of metabolic interactions between the oocyte and somatic cells has been limited due to dynamic nature of ovarian follicle development, thus warranting a systems approach. Results: Here, we developed a genome-scale metabolic model of the mouse ovarian follicle. This model was constructed using an updated mouse general metabolic model (Mouse Recon 2) and contains several key ovarian follicle development metabolic pathways. We used this model to characterize the changes in the metabolism of each follicular cell type (i.e., oocyte, granulosa cells, including cumulus and mural cells), during ovarian follicle development in vivo. Using this model, we predicted major metabolic pathways that are differentially active across multiple follicle stages. We identified a set of possible secreted and consumed metabolites that could potentially serve as biomarkers for monitoring follicle development, as well as metabolites for addition to in vitro culture media that support the growth and maturation of primordial follicles. Conclusions: Our systems approach to model follicle metabolism can guide future experimental studies to validate the model results and improve oocyte maturation approaches and support growth of primordial follicles in vitro. [ABSTRACT FROM AUTHOR]

Candida glabrata is a major cause of candidiasis and the second most frequent opportunistic yeast pathogen. Its infectious and antifungal mechanisms are globally regulated by the transcription systems of pathogenic fungi. In this study, we reconstructed the genome-scale transcriptional regulatory network (TRN) of C. glabrata, consisting of 6634 interactive relationships between 145 transcription factors and 3230 target genes, based on genomic and transcriptomic data. The C. glabrata TRN was found to have a typical topological structure and significant network cohesiveness. Moreover, this network could be functionally divided into several sub-networks, including networks involving carbon, nitrogen, growth-associated metabolic profiles, stress response to acidity, hyperosmosis, peroxidation, hypoxia and virulence. Furthermore, by integrating the genome-scale metabolic model of C. glabrata, six essential metabolites and eight related enzymes were systematically selected as drug targets. Overall, elucidation of the genome-scale TRN of C. glabrata has expanded our knowledge of the contents and structures of microbial regulatory networks and improved our understanding of the regulatory behaviors of growth, metabolism and gene expression programs in response to environmental stimuli. Reconstruction of genome-scale transcriptional regulatory network for non-model microbe Candida glabrata and integration with metabolic model for the in silico drug targets. [ABSTRACT FROM AUTHOR]

The emergence and spread of metallo-beta-lactamase–producing multidrug-resistant (MDR) Klebsiella pneumoniae is a serious public health threat, which is further complicated by the increased prevalence of colistin resistance. The link between antimicrobial resistance acquired by strains of Klebsiella and their unique metabolic capabilities has not been determined. Here, we reconstruct genome-scale metabolic models for 22 K. pneumoniae strains with various resistance profiles to different antibiotics, including two strains exhibiting colistin resistance isolated from Cairo, Egypt. We use the models to predict growth capabilities on 265 different sole carbon, nitrogen, sulfur, and phosphorus sources for all 22 strains. Alternate nitrogen source utilization of glutamate, arginine, histidine, and ethanolamine among others provided discriminatory power for identifying resistance to amikacin, tetracycline, and gentamicin. Thus, genome-scale model based predictions of growth capabilities on alternative substrates may lead to construction of classification trees that are indicative of antibiotic resistance in Klebsiella isolates. [ABSTRACT FROM AUTHOR]

Summary: Photoacclimation consists of short‐ and long‐term strategies used by photosynthetic organisms to adapt to dynamic light environments. Observable photophysiology changes resulting from these strategies have been used in coarse‐grained models to predict light‐dependent growth and photosynthetic rates. However, the contribution of the broader metabolic network, relevant to species‐specific strategies and fitness, is not accounted for in these simple models.We incorporated photophysiology experimental data with genome‐scale modeling to characterize organism‐level, light‐dependent metabolic changes in the model diatom Phaeodactylum tricornutum. Oxygen evolution and photon absorption rates were combined with condition‐specific biomass compositions to predict metabolic pathway usage for cells acclimated to four different light intensities.Photorespiration, an ornithine‐glutamine shunt, and branched‐chain amino acid metabolism were hypothesized as the primary intercompartment reductant shuttles for mediating excess light energy dissipation. Additionally, simulations suggested that carbon shunted through photorespiration is recycled back to the chloroplast as pyruvate, a mechanism distinct from known strategies in photosynthetic organisms.Our results suggest a flexible metabolic network in P. tricornutum that tunes intercompartment metabolism to optimize energy transport between the organelles, consuming excess energy as needed. Characterization of these intercompartment reductant shuttles broadens our understanding of energy partitioning strategies in this clade of ecologically important primary producers. [ABSTRACT FROM AUTHOR]

Abstract There is great interest in engineering photoautotrophic metabolism to generate bioproducts of societal importance. Despite the success in employing genome-scale modeling coupled with flux balance analysis to engineer heterotrophic metabolism, the lack of proper constraints necessary to generate biologically realistic predictions has hindered broad application of this methodology to phototrophic metabolism. Here we describe a methodology for constraining genome-scale models of photoautotrophy in the cyanobacteria Synechococcus elongatus PCC 7942. Experimental photophysiology parameters coupled to genome-scale flux balance analysis resulted in accurate predictions of growth rates and metabolic reaction fluxes at low and high light conditions. Additionally, by constraining photon uptake fluxes, we characterized the metabolic cost of excess excitation energy. The predicted energy fluxes were consistent with known light-adapted phenotypes in cyanobacteria. Finally, we leveraged the modeling framework to characterize existing photoautotrophic and photomixtotrophic engineering strategies for 2,3-butanediol production in S. elongatus. This methodology, applicable to genome-scale modeling of all phototrophic microorganisms, can facilitate the use of flux balance analysis in the engineering of light-driven metabolism. Highlights • Light uptake and oxygen evolution accurately constrain phototrophic metabolism. • Metabolic fluxes predicted by these constraints recapitulate experimental results. • Alternative electron transport and photodamage repair rates are quantified. • Photo- and photomixtotrophic 2,3-butanediol engineering designs are characterized. • The constraints can be extended to genome-scale models of any phototroph. [ABSTRACT FROM AUTHOR]

Genome-scale stoichiometric modeling methods, in particular Flux Balance Analysis (FBA) and variations thereof, are widely used to investigate cell metabolism and to optimize biotechnological processes. Given (1) a metabolic network, which can be reconstructed from an organism’s genome sequence, and (2) constraints on reaction rates, which may be based on measured nutrient uptake rates, FBA predicts which reactions maximize an objective flux, usually the production of cell components. Although FBA solutions may accurately predict the metabolic behavior of a cell, the actual flux predictions are often hard to interpret. This is especially the case for conditions with many constraints, such as for organisms growing in rich nutrient environments: it remains unclear why a certain solution was optimal. Here, we rationalize FBA solutions by explaining for which properties the optimal combination of metabolic strategies is selected. We provide a graphical formalism in which the selection of solutions can be visualized; we illustrate how this perspective provides a glimpse of the logic that underlies genome-scale modeling by applying our formalism to models of various sizes.

Flux balance analysis (FBA) is a computational methodology to model and analyze the metabolic behavior of cells. In this chapter, we break down the key steps for formulating an FBA model and other FBA-derived methodologies in the context of mammalian cell biology, including strain design, developing cell line-specific models, and conducting flux sampling. We provide annotated COBRApy code for each step to show how it would work in practice., (© 2024. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Photosynthetic cyanobacteria are important primary producers and model organisms for studying photosynthesis and elements cycling on earth. Due to the ability to absorb sunlight and utilize carbon dioxide, cyanobacteria have also been proposed as renewable chassis for carbon-neutral “microbial cell factories”. Recent progresses on cyanobacterial synthetic biology have led to the successful production of more than two dozen of fuels and fine chemicals directly from CO 2 , demonstrating their potential for scale-up application in the future. However, compared with popular heterotrophic chassis like Escherichia coli and Saccharomyces cerevisiae , where abundant genetic tools are available for manipulations at levels from single gene, pathway to whole genome, limited genetic tools are accessible to cyanobacteria. Consequently, this significant technical hurdle restricts both the basic biological researches and further development and application of these renewable systems. Though still lagging the heterotrophic chassis, the vital roles of genetic tools in tuning of gene expression, carbon flux re-direction as well as genome-wide manipulations have been increasingly recognized in cyanobacteria. In recent years, significant progresses on developing and introducing new and efficient genetic tools have been made for cyanobacteria, including promoters, riboswitches, ribosome binding site engineering, clustered regularly interspaced short palindromic repeats/CRISPR-associated nuclease (CRISPR/Cas) systems, small RNA regulatory tools and genome-scale modeling strategies. In this review, we critically summarize recent advances on development and applications as well as technical limitations and future directions of the genetic tools in cyanobacteria. In addition, toolboxes feasible for using in large-scale cultivation are also briefly discussed. [ABSTRACT FROM AUTHOR]

Neurology research and clinical practice are transforming toward postgenomics integrative biology. One such example is the study of human brain metabolism that is highly sophisticated due to reactions occurring in and between the astrocytes and neurons. Because of the inherent difficulty of performing experimental studies in human brain, metabolic network modeling has grown in importance to decipher the contribution of brain metabolite kinetics to human health and disease. Multiomics system science-driven metabolic models, using genome-scale and transcriptomics Big Data, offer the promise of new insights on metabolic networks in human brain. Added to this, the availability of omics technologies in both developed and developing world, neurology research, and clinical practice ought to be repositioned with a view to systems medicine. In this expert analysis, we present a critical and in-depth overview of the basic tenets of human brain metabolism, together with the most recent metabolic modeling strategies and computational studies of brain in health and neurological diseases. Human genome-scale metabolic models developed in a both global and brain-specific manner and multiomics synthesis of knowledge are highlighted in particular. We conclude by underscoring the value of multiomics modeling for metabolic diseases and computational investigations of the brain networks, with a view to unlocking the pathophysiology of Alzheimer's disease, Parkinson's disease, migraine, stroke, epilepsy, and multiple sclerosis, among other neurological disorders of importance for global health. [ABSTRACT FROM AUTHOR]

The primary aims and challenges associated with microbial fermentation include achieving faster cell growth, higher productivity, and more robust production processes. Genome-scale biological models, predicting the formation of an interaction among genetic materials, enzymes, and metabolites, constitute a systematic and comprehensive platform to analyze and optimize the microbial growth and production of biological products. Genome-scale biological models can help optimize microbial growth-associated traits by simulating biomass formation, predicting growth rates, and identifying the requirements for cell growth. With regard to microbial product biosynthesis, genome-scale biological models can be used to design product biosynthetic pathways, accelerate production efficiency, and reduce metabolic side effects, leading to improved production performance. The present review discusses the development of microbial genome-scale biological models since their emergence and emphasizes their pertinent application in improving industrial microbial fermentation of biological products. [ABSTRACT FROM AUTHOR]

Metabolic engineering modifies cellular function to address various biochemical applications. Underlying metabolic engineering efforts are a host of tools and knowledge that are integrated to enable successful outcomes. Concurrent development of computational and experimental tools has enabled different approaches to metabolic engineering. One approach is to leverage knowledge and computational tools to prospectively predict designs to achieve the desired outcome. An alternative approach is to utilize combinatorial experimental tools to empirically explore the range of cellular function and to screen for desired traits. This mini-review focuses on computational systems biology and synthetic biology tools that can be used in combination for prospective in silico strain design.

In recent years, metabolic engineering of microorganisms has attained much research interest to produce biofuels and industrially pertinent chemicals. Owing to the relatively fast growth rate, genetic malleability, and carbon neutral production process, cyanobacteria has been recognized as a specialized microorganism with a significant biotechnological perspective. Metabolically engineering cyanobacterial strains have shown great potential for the photosynthetic production of an array of valuable native or non-native chemicals and metabolites with profound agricultural and pharmaceutical significance using CO2 as a building block. In recent years, substantial improvements in developing and introducing novel and efficient genetic tools such as genome-scale modeling, high throughput omics analyses, synthetic/system biology tools, metabolic flux analysis and clustered regularly interspaced short palindromic repeats (CRISPR)-associated nuclease (CRISPR/cas) systems have been made for engineering cyanobacterial strains. Use of these tools and technologies has led to a greater understanding of the host metabolism, as well as endogenous and heterologous carbon regulation mechanisms which consequently results in the expansion of maximum productive ability and biochemical diversity. This review summarizes recent advances in engineering cyanobacteria to produce biofuel and industrially relevant fine chemicals of high interest. Moreover, the development and applications of cutting-edge toolboxes such as the CRISPR-cas9 system, synthetic biology, high-throughput “omics”, and metabolic flux analysis to engineer cyanobacteria for large-scale cultivation are also discussed.