Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0

Exploring metabolic effects of dipeptide feed media on CHO cell cultures by in silico model-guided flux analysis

There is a growing interest in perfusion or continuous processes to achieve higher productivity of biopharmaceuticals in mammalian cell culture, specifically Chinese hamster ovary (CHO) cells, towards advanced biomanufacturing. These intensified bioprocesses highly require concentrated feed media in order to counteract their dilution effects. However, designing such condensed media formulation poses several challenges, particularly regarding the stability and solubility of specific amino acids. To address the difficulty and complexity in relevant media development, the biopharmaceutical industry has recently suggested forming dipeptides by combining one from problematic amino acids with selected pairs to compensate for limitations. In this study, we combined one of the lead amino acids, L-tyrosine, which is known for its poor solubility in water due to its aromatic ring and hydroxyl group, with glycine as the partner, thus forming glycyl-L-tyrosine (GY) dipeptide. Subsequently, we investigated the utilization of GY dipeptide during fed-batch cultures of IgG-producing CHO cells, by changing its concentrations (0.125 × , 0.25 × , 0.5 × , 1.0 × , and 2.0 ×). Multivariate statistical analysis of culture profiles was then conducted to identify and correlate the most significant nutrients with the production, followed by in silico model-guided analysis to systematically evaluate their effects on the culture performance, and elucidate metabolic states and cellular behaviors. As such, it allowed us to explain how the cells can more efficiently utilize GY dipeptide with respect to the balance of cofactor regeneration and energy distribution for the required biomass and protein synthesis. For example, our analysis results uncovered specific amino acids (Asn and Gln) and the 0.5 × GY dipeptide in the feed medium synergistically alleviated the metabolic bottleneck, resulting in enhanced IgG titer and productivity. In the validation experiments, we tested and observed that lower levels of Asn and Gln led to decreased secretion of toxic metabolites, enhanced longevity, and elevated specific cell growth and titer. KEY POINTS: • Explored the optimal Tyr dipeptide for the enhanced CHO cell culture performance • Systematically analyzed effects of dipeptide media by model-guided approach • Uncovered synergistic metabolic utilization of amino acids with dipeptide.

Harnessing meta-omics to unveil and mitigate methane emissions in ruminants: Integrative approaches and future directions

Methane emissions from enteric fermentation present a dual challenge globally: they not only contribute significantly to atmospheric greenhouse gases but also represent a considerable energy loss for ruminant animals. Utilizing high-throughput omics technologies to analyze rumen microbiome samples (meta-omics, i.e., metagenomics, metatranscriptomics, metaproteomics, metabolomics) holds vast potential for uncovering the intricate interplay between diet, microbiota, and methane emissions in these animals. The primary obstacle is the effective integration of diverse meta-omic approaches and their broader application across different ruminant species. Genetic variability significantly impacts methane production in ruminants, suggesting that genomic selection could be a viable strategy to reduce emissions. While substantial research has been conducted on the microbiological aspects of methane production, there remains a critical need to delineate the specific genetic interactions between the host and its microbiome. Advancements in meta-omics technologies are poised to shed light on these interactions, enhancing our understanding of the genetic factors that govern methane output. This review explores the potential of meta-omics to accelerate genetic advancements that could lead to reduced methane emissions in ruminants. By employing a systems biology approach, the integration of various omics technologies allows for the identification of key genomic regions and genetic markers linked to methane production. These markers can then be leveraged in selective breeding programs to cultivate traits associated with lower emissions. Moreover, the review addresses current challenges in applying genomic selection for this purpose and discusses how omics technologies can overcome these obstacles. The systematic integration and analysis of diverse biological data provide deeper insights into the genetic underpinnings and overall biology of methane production traits in ruminants. Ultimately, this comprehensive approach not only aids in reducing the environmental impact of agriculture but also contributes to the sustainability and efficiency of livestock management.

Applications of genome-scale metabolic models to the study of human diseases: A systematic review

BACKGROUND AND OBJECTIVES

Genome-scale metabolic networks (GEMs) represent a valuable modeling and computational tool in the broad field of systems biology. Their ability to integrate constraints and high-throughput biological data enables the study of intricate metabolic aspects and processes of different cell types and conditions. The past decade has witnessed an increasing number and variety of applications of GEMs for the study of human diseases, along with a huge effort aimed at the reconstruction, integration and analysis of a high number of organisms. This paper presents a systematic review of the scientific literature, to pursue several important questions about the application of constraint-based modeling in the investigation of human diseases. Hopefully, this paper will provide a useful reference for researchers interested in the application of modeling and computational tools for the investigation of metabolic-related human diseases.

METHODS

This systematic review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Elsevier Scopus®, National Library of Medicine PubMed® and Clarivate Web of Science™ databases were enquired, resulting in 566 scientific articles. After applying exclusion and eligibility criteria, a total of 169 papers were selected and individually examined.

RESULTS

The reviewed papers offer a thorough and up-to-date picture of the latest modeling and computational approaches, based on genome-scale metabolic models, that can be leveraged for the investigation of a large variety of human diseases. The numerous studies have been categorized according to the clinical research area involved in the examined disease. Furthermore, the paper discusses the most typical approaches employed to derive clinically-relevant information using the computational models.

CONCLUSIONS

The number of scientific papers, utilizing GEM-based approaches for the investigation of human diseases, suggests an increasing interest in these types of approaches; hopefully, the present review will represent a useful reference for scientists interested in applying computational modeling approaches to investigate the aetiopathology of human diseases; we also hope that this work will foster the development of novel applications and methods for the discovery of clinically-relevant insights on metabolic-related diseases.

Integration of proteomic data with genome-scale metabolic models: A methodological overview

The integration of proteomics data with constraint-based reconstruction and analysis (COBRA) models plays a pivotal role in understanding the relationship between genotype and phenotype and bridges the gap between genome-level phenomena and functional adaptations. Integrating a generic genome-scale model with information on proteins enables generation of a context-specific metabolic model which improves the accuracy of model prediction. This review explores methodologies for incorporating proteomics data into genome-scale models. Available methods are grouped into four distinct categories based on their approach to integrate proteomics data and their depth of modeling. Within each category section various methods are introduced in chronological order of publication demonstrating the progress of this field. Furthermore, challenges and potential solutions to further progress are outlined, including the limited availability of appropriate in vitro data, experimental enzyme turnover rates, and the trade-off between model accuracy, computational tractability, and data scarcity. In conclusion, methods employing simpler approaches demand fewer kinetic and omics data, consequently leading to a less complex mathematical problem and reduced computational expenses. On the other hand, approaches that delve deeper into cellular mechanisms and aim to create detailed mathematical models necessitate more extensive kinetic and omics data, resulting in a more complex and computationally demanding problem. However, in some cases, this increased cost can be justified by the potential for more precise predictions.

Computational Synthetic Biology Enabled through JAX: A Showcase

Mathematical modeling is indispensable in synthetic biology but remains underutilized. Tackling problems, from optimizing gene networks to simulating intracellular dynamics, can be facilitated by the ever-growing body of modeling approaches, be they mechanistic, stochastic, data-driven, or AI-enabled. Thanks to progress in the AI community, robust frameworks have emerged to enable researchers to access complex computational hardware and compilation. Previously, these frameworks focused solely on deep learning, but they have been developed to the point where running different forms of computation is relatively simple, as made possible, notably, by the JAX library. Running simulations at scale on GPUs speeds up research, which compounds enable larger-scale experiments and greater usability of code. As JAX remains underexplored in computational biology, we demonstrate its utility in three example projects ranging from synthetic biology to directed evolution, each with an accompanying demonstrative Jupyter notebook. We hope that these tutorials serve to democratize the flexible scaling, faster run-times, easy GPU portability, and mathematical enhancements (such as automatic differentiation) that JAX brings, all with only minor restructuring of code.

Protocol for genome-scale differential flux analysis to interrogate metabolic differences from gene expression data

Deciphering the functional differences between diseased and healthy cells requires understanding the alterations in biochemical flux patterns. We present a genome-scale differential flux analysis (GS-DFA) protocol to elucidate these metabolic disparities by integrating condition-specific gene expression data into the human genome-scale metabolic model (humanGEM). In this protocol, we describe the steps to normalize and integrate data into the humanGEM and analyze differential flux across the biochemical network between diseased and healthy cells. For complete details on the use and execution of this protocol, please refer to Nanda et al.1.

Enhanced growth of wheat in contaminated fields via synthetic microbiome as revealed by genome-scale metabolic modeling

The relationship between plants and soil microbial communities is complex and subtle, with microbes playing a crucial role in plant growth. Autochthonous bioaugmentation and nutrient biostimulation are promising bioremediation methods for herbicides in contaminated agricultural soils, but how microbes interact to promote biodegradation and plant growth on barren fields, especially in response to the treatment of the herbicide bromoxynil after wheat seedlings, remains poorly understood. In this study, we explored the microbial community reassembly process from the three-leaf stage to the tillering stage of wheat and put forward the idea of using the overlapping results of three methods (network Zi-Pi analysis, LEfSe analysis, and Random Forest analysis) as keystones for the simplification and optimization of key microbial species in the soil. Then we used genome-scale metabolic models (GSMMs) to design a targeted synthetic microbiome for promoting wheat seedling growing. The results showed that carbon source was more helpful in enriching soil microbial diversity and promoting the role of functional microbial communities, which facilitated the degradation of bromoxynil. Designed a multifunctional synthetic consortium consisting of seven non-degraders which unexpectedly assisted in the degradation of indigenous bacteria, which increased the degradation rate of bromoxynil by 2.05 times, and when adding nutritional supplementation, it increased the degradation rate by 3.65 times. In summary, this study provides important insights for rational fertilization and precise microbial consortium management to improve plant seedling growth in contaminated fields.

Multi-omics analyses highlight molecular differences between clinical and neuropathological diagnoses in Alzheimer's disease

Both clinical diagnosis and neuropathological diagnosis are commonly used in literature to categorize individuals as Alzheimer's disease (AD) or non-AD in omics analyses. Whether these diagnostic strategies result in distinct profiles of molecular abnormalities is poorly understood. Here, we analysed one of the most commonly used AD omics datasets in the literature from the Religious Orders Study and Memory and Aging Project (ROSMAP) cohort and compared the two diagnosis strategies using brain transcriptome and metabolome by grouping individuals as non-AD and AD according to clinical or neuropathological diagnosis separately. Differentially expressed genes, associated pathways related with AD hallmarks and AD-related genes showed that the categorization based on neuropathological diagnosis more accurately reflects the disease state at the molecular level than the categorization based on clinical diagnosis. We further identified consensus biomarker candidates between the two diagnosis strategies such as 5-hydroxylysine, sphingomyelin and 1-myristoyl-2-palmitoyl-GPC as metabolite biomarkers and sphingolipid metabolism as a pathway biomarker, which could be robust AD biomarkers since they are independent of diagnosis strategies. We also used consensus AD and consensus non-AD individuals between the two diagnostic strategies to train a machine-learning based model, which we used to classify the individuals who were cognitively normal but diagnosed as AD based on neuropathological diagnosis (asymptomatic AD individuals). The majority of these individuals were classified as consensus AD patients for both omics data types. Our study provides a detailed characterization of both diagnostic strategies in terms of the association of the corresponding multi-omics profiles with AD.

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

Reconstruction and Analysis of a Genome-Scale Metabolic Model of Acinetobacter lwoffii

Acinetobacter lwoffii is widely considered to be a harmful bacterium that is resistant to medicines and disinfectants. A. lwoffii NL1 degrades phenols efficiently and shows promise as an aromatic compound degrader in antibiotic-contaminated environments. To gain a comprehensive understanding of A. lwoffii, the first genome-scale metabolic model of A. lwoffii was constructed using semi-automated and manual methods. The iNX811 model, which includes 811 genes, 1071 metabolites, and 1155 reactions, was validated using 39 unique carbon and nitrogen sources. Genes and metabolites critical for cell growth were analyzed, and 12 essential metabolites (mainly in the biosynthesis and metabolism of glycan, lysine, and cofactors) were identified as antibacterial drug targets. Moreover, to explore the metabolic response to phenols, metabolic flux was simulated by integrating transcriptomics, and the significantly changed metabolism mainly included central carbon metabolism, along with some transport reactions. In addition, the addition of substances that effectively improved phenol degradation was predicted and validated using the model. Overall, the reconstruction and analysis of model iNX811 helped to study the antimicrobial systems and biodegradation behavior of A. lwoffii.

Origin of biogeographically distinct ecotypes during laboratory evolution

Resource partitioning is central to the incredible productivity of microbial communities, including gigatons in annual methane emissions through syntrophic interactions. Previous work revealed how a sulfate reducer (Desulfovibrio vulgaris, Dv) and a methanogen (Methanococcus maripaludis, Mm) underwent evolutionary diversification in a planktonic context, improving stability, cooperativity, and productivity within 300-1000 generations. Here, we show that mutations in just 15 Dv and 7 Mm genes within a minimal assemblage of this evolved community gave rise to co-existing ecotypes that were spatially enriched within a few days of culturing in a fluidized bed reactor. The spatially segregated communities partitioned resources in the simulated subsurface environment, with greater lactate utilization by attached Dv but partial utilization of resulting H2 by low affinity hydrogenases of Mm in the same phase. The unutilized H2 was scavenged by high affinity hydrogenases of planktonic Mm, producing copious amounts of methane. Our findings show how a few mutations can drive resource partitioning amongst niche-differentiated ecotypes, whose interplay synergistically improves productivity of the entire mutualistic community.

Insights into yeast response to chemotherapeutic agent through time series genome-scale metabolic models

Organism-specific genome-scale metabolic models (GSMMs) can unveil molecular mechanisms within cells and are commonly used in diverse applications, from synthetic biology, biotechnology, and systems biology to metabolic engineering. There are limited studies incorporating time-series transcriptomics in GSMM simulations. Yeast is an easy-to-manipulate model organism for tumor research. Here, a novel approach (TS-GSMM) was proposed to integrate time-series transcriptomics with GSMMs to narrow down the feasible solution space of all possible flux distributions and attain time-series flux samples. The flux samples were clustered using machine learning techniques, and the clusters' functional analysis was performed using reaction set enrichment analysis. A time series transcriptomics response of Yeast cells to a chemotherapeutic reagent-doxorubicin-was mapped onto a Yeast GSMM. Eleven flux clusters were obtained with our approach, and pathway dynamics were displayed. Induction of fluxes related to bicarbonate formation and transport, ergosterol and spermidine transport, and ATP production were captured. Integrating time-series transcriptomics data with GSMMs is a promising approach to reveal pathway dynamics without any kinetic modeling and detects pathways that cannot be identified through transcriptomics-only analysis. The codes are available at https://github.com/karabekmez/TS-GSMM.

Interspecies synergistic interactions mediated by cofactor exchange enhance stress tolerance by inducing biofilm formation

Metabolic exchange plays a crucial role in shaping microbial community interactions and functions, including the exchange of small molecules such as cofactors. Cofactors are fundamental to enzyme catalytic activities; however, the role of cofactors in microbial stress tolerance is unclear. Here, we constructed a synergistic consortium containing two strains that could efficiently mineralize di-(2-ethylhexyl) phthalate under hyperosmotic stress. Integration of transcriptomic analysis, metabolic profiling, and a genome-scale metabolic model (GEM) facilitated the discovery of the potential mechanism of microbial interactions. Multi-omics analysis revealed that the vitamin B12-dependent methionine-folate cycle could be a key pathway for enhancing the hyperosmotic stress tolerance of synergistic consortium. Further GEM simulations revealed interspecies exchange of S-adenosyl-L-methionine and riboflavin, cofactors needed for vitamin B12 biosynthesis, which was confirmed by in vitro experiments. Overall, we proposed a new mechanism of bacterial hyperosmotic stress tolerance: bacteria might promote the production of vitamin B12 to enhance biofilm formation, and the species collaborate with each other by exchanging cofactors to improve consortium hyperosmotic stress tolerance. These findings offer new insights into the role of cofactors in microbial interactions and stress tolerance and are potentially exploitable for environmental remediation.

IMPORTANCE

Metabolic interactions (also known as cross-feeding) are thought to be ubiquitous in microbial communities. Cross-feeding is the basis for many positive interactions (e.g., mutualism) and is a primary driver of microbial community assembly. In this study, a combination of multi-omics analysis and metabolic modeling simulation was used to reveal the metabolic interactions of a synthetic consortium under hyperosmotic stress. Interspecies cofactor exchange was found to promote biofilm formation under hyperosmotic stress. This provides a new perspective for understanding the role of metabolic interactions in microbial communities to enhance environmental adaptation, which is significant for improving the efficiency of production activities and environmental bioremediation.

Cell Geometry and Membrane Protein Crowding Constrain Growth Rate, Overflow Metabolism, Respiration, and Maintenance Energy

A metabolic theory is presented for predicting maximum growth rate, overflow metabolism, respiration efficiency, and maintenance energy flux based on the intersection of cell geometry, membrane protein crowding, and metabolism. The importance of cytosolic macromolecular crowding on phenotype has been established in the literature but the importance of surface area has been largely overlooked due to incomplete knowledge of membrane properties. We demonstrate that the capacity of the membrane to host proteins increases with growth rate offsetting decreases in surface area-to-volume ratios (SA:V). This increase in membrane protein is hypothesized to be essential to competitive Escherichia coli phenotypes. The presented membrane-centric theory uses biophysical properties and metabolic systems analysis to successfully predict the phenotypes of E. coli K-12 strains, MG1655 and NCM3722, which are genetically similar but have SA:V ratios that differ up to 30%, maximum growth rates on glucose media that differ by 40%, and overflow phenotypes that start at growth rates that differ by 80%. These analyses did not consider cytosolic macromolecular crowding, highlighting the distinct properties of the presented theory. Cell geometry and membrane protein crowding are significant biophysical constraints on phenotype and provide a theoretical framework for improved understanding and control of cell biology.

Systems Metabolic Engineering to Elucidate and Enhance Intestinal Metabolic Activities of Escherichia coli Nissle 1917

Escherichia coli Nissle 1917 (EcN) is one of the most widely used probiotics to treat gastrointestinal diseases. Recently, many studies have engineered EcN to release therapeutic proteins to treat specific diseases. However, because EcN exhibits intestinal metabolic activities, it is difficult to predict outcomes after administration. In silico and fermentation profiles revealed mucin metabolism of EcN. Multiomics revealed that fucose metabolism contributes to the intestinal colonization of EcN by enhancing the synthesis of flagella and nutrient uptake. The multiomics results also revealed that excessive intracellular trehalose synthesis in EcN, which is responsible for galactose metabolism, acts as a metabolic bottleneck, adversely affecting growth. To improve the ability of EcN to metabolize galactose, otsAB genes for trehalose synthesis were deleted, resulting in the ΔotsAB strain; the ΔotsAB strain exhibited a 1.47-fold increase in the growth rate and a 1.37-fold increase in the substrate consumption rate relative to wild-type EcN.

Personalized metabolic whole-body models for newborns and infants predict growth and biomarkers of inherited metabolic diseases

Comprehensive whole-body models (WBMs) accounting for organ-specific dynamics have been developed to simulate adult metabolism, but such models do not exist for infants. Here, we present a resource of 360 organ-resolved, sex-specific models of newborn and infant metabolism (infant-WBMs) spanning the first 180 days of life. These infant-WBMs were parameterized to represent the distinct metabolic characteristics of newborns and infants, including nutrition, energy requirements, and thermoregulation. We demonstrate that the predicted infant growth was consistent with the recommendation by the World Health Organization. We assessed the infant-WBMs' reliability and capabilities for personalization by simulating 10,000 newborns based on their blood metabolome and birth weight. Furthermore, the infant-WBMs accurately predicted changes in known biomarkers over time and metabolic responses to treatment strategies for inherited metabolic diseases. The infant-WBM resource holds promise for personalized medicine, as the infant-WBMs could be a first step to digital metabolic twins for newborn and infant metabolism.

Epigenetic control of metabolic identity across cell types

Constraint-based network modelling is a powerful tool for analysing cellular metabolism at genomic scale. Here, we conducted an integrative analysis of metabolic networks reconstructed from RNA-seq data with paired epigenomic data from the EpiATLAS resource of the International Human Epigenome Consortium (IHEC). Applying a state-of-the-art contextualisation algorithm, we reconstructed metabolic networks across 1,555 samples corresponding to 58 tissues and cell types. Analysis of these networks revealed the distribution of metabolic functionalities across human cell types and provides a compendium of human metabolic activity. This integrative approach allowed us to define, across tissues and cell types, i) reactions that fulfil the basic metabolic processes (core metabolism), and ii) cell type-specific functions (unique metabolism), that shape the metabolic identity of a cell or a tissue. Integration with EpiATLAS-derived cell type-specific gene-level chromatin states and enhancer-gene interactions identified enhancers, transcription factors, and key nodes controlling core and unique metabolism. Transport and first reactions of pathways were enriched for high expression, active chromatin state, and Polycomb-mediated repression in cell types where pathways are inactive, suggesting that key nodes are targets of repression. This integrative analysis forms the basis for identifying regulation points that control metabolic identity in human cells.

Multi-omic analysis of bat versus human fibroblasts reveals altered central metabolism

Bats have unique characteristics compared to other mammals, including increased longevity and higher resistance to cancer and infectious disease. While previous studies have analyzed the metabolic requirements for flight, it is still unclear how bat metabolism supports these unique features, and no study has integrated metabolomics, transcriptomics, and proteomics to characterize bat metabolism. In this work, we performed a multi-omics data analysis using a computational model of metabolic fluxes to identify fundamental differences in central metabolism between primary lung fibroblast cell lines from the black flying fox fruit bat (Pteropus alecto) and human. Bat cells showed higher expression levels of Complex I components of electron transport chain (ETC), but, remarkably, a lower rate of oxygen consumption. Computational modeling interpreted these results as indicating that Complex II activity may be low or reversed, similar to an ischemic state. An ischemic-like state of bats was also supported by decreased levels of central metabolites and increased ratios of succinate to fumarate in bat cells. Ischemic states tend to produce reactive oxygen species (ROS), which would be incompatible with the longevity of bats. However, bat cells had higher antioxidant reservoirs (higher total glutathione and higher ratio of NADPH to NADP) despite higher mitochondrial ROS levels. In addition, bat cells were more resistant to glucose deprivation and had increased resistance to ferroptosis, one of the characteristics of which is oxidative stress. Thus, our studies revealed distinct differences in the ETC regulation and metabolic stress responses between human and bat cells.

Cell factory design with advanced metabolic modelling empowered by artificial intelligence

Advances in synthetic biology and artificial intelligence (AI) have provided new opportunities for modern biotechnology. High-performance cell factories, the backbone of industrial biotechnology, are ultimately responsible for determining whether a bio-based product succeeds or fails in the fierce competition with petroleum-based products. To date, one of the greatest challenges in synthetic biology is the creation of high-performance cell factories in a consistent and efficient manner. As so-called white-box models, numerous metabolic network models have been developed and used in computational strain design. Moreover, great progress has been made in AI-powered strain engineering in recent years. Both approaches have advantages and disadvantages. Therefore, the deep integration of AI with metabolic models is crucial for the construction of superior cell factories with higher titres, yields and production rates. The detailed applications of the latest advanced metabolic models and AI in computational strain design are summarized in this review. Additionally, approaches for the deep integration of AI and metabolic models are discussed. It is anticipated that advanced mechanistic metabolic models powered by AI will pave the way for the efficient construction of powerful industrial chassis strains in the coming years.

Panera: An innovative framework for surmounting uncertainty in microbial community modeling using pan-genera metabolic models

Utilization of 16S rRNA data in constraint-based modeling to characterize microbial communities confronts a major hurdle of lack of species-level resolution, impeding the construction of community models. We introduce "Panera," an innovative framework designed to model communities under this uncertainty and yet perform metabolic inferences using pan-genus metabolic models (PGMMs). We demonstrated PGMMs' utility for comprehending the metabolic capabilities of a genus and in characterizing community models using amplicon data. The unique, adaptable nature of PGMMs unlocks their potential in building hybrid communities, combining genome-scale metabolic models (GSMMs) and PGMMs. Notably, these models provide predictions comparable to the standard GSMM-based community models, while achieving a nearly 46% reduction in error compared to the genus model-based communities. In essence, "Panera" presents a potent and effective approach to aid in metabolic modeling by enabling robust predictions of community metabolic potential when dealing with amplicon data, and offers insights into genus-level metabolic landscapes.

Mixotrophic growth of a ubiquitous marine diatom

Diatoms are major players in the global carbon cycle, and their metabolism is affected by ocean conditions. Understanding the impact of changing inorganic nutrients in the oceans on diatoms is crucial, given the changes in global carbon dioxide levels. Here, we present a genome-scale metabolic model (iMK1961) for Cylindrotheca closterium, an in silico resource to understand uncharacterized metabolic functions in this ubiquitous diatom. iMK1961 represents the largest diatom metabolic model to date, comprising 1961 open reading frames and 6718 reactions. With iMK1961, we identified the metabolic response signature to cope with drastic changes in growth conditions. Comparing model predictions with Tara Oceans transcriptomics data unraveled C. closterium's metabolism in situ. Unexpectedly, the diatom only grows photoautotrophically in 21% of the sunlit ocean samples, while the majority of the samples indicate a mixotrophic (71%) or, in some cases, even a heterotrophic (8%) lifestyle in the light. Our findings highlight C. closterium's metabolic flexibility and its potential role in global carbon cycling.

Understanding pulmonary hypertension: the need for an integrative metabolomics and transcriptomics approach

Pulmonary hypertension (PH), characterised by mean pulmonary arterial pressure (mPAP) >20 mm Hg at rest, is a complex pathophysiological disorder associated with multiple clinical conditions. The high prevalence of the disease along with increased mortality and morbidity makes it a global health burden. Despite major advances in understanding the disease pathophysiology, much of the underlying complex molecular mechanism remains to be elucidated. Lack of a robust diagnostic test and specific therapeutic targets also poses major challenges. This review provides a comprehensive update on the dysregulated pathways and promising candidate markers identified in PH patients using the transcriptomics and metabolomics approach. The review also highlights the need of using an integrative multi-omics approach for obtaining insight into the disease at a molecular level. The integrative multi-omics/pan-omics approach envisaged to help in bridging the gap from genotype to phenotype is outlined. Finally, the challenges commonly encountered while conducting omics-driven studies are also discussed.

Revealing microbiota characteristics and predicting flavor-producing sub-communities in Nongxiangxing baijiu pit mud through metagenomic analysis and metabolic modeling

The microorganisms of the pit mud (PM) of Nongxiangxing baijiu (NXXB) have an important role in the synthesis of flavor substances, and they determine attributes and quality of baijiu. Herein, we utilize metagenomics and genome-scale metabolic models (GSMMs) to investigate the microbial composition, metabolic functions in PM microbiota, as well as to identify microorganisms and communities linked to flavor compounds. Metagenomic data revealed that the most prevalent assembly of bacteria and archaea was Proteiniphilum, Caproicibacterium, Petrimonas, Lactobacillus, Clostridium, Aminobacterium, Syntrophomonas, Methanobacterium, Methanoculleus, and Methanosarcina. The important enzymes ofPMwere in bothGH and GT familymetabolism. A total of 38 high-quality metagenome-assembled genomes (MAGs) were obtained, including those at the family level (n = 13), genus level (n = 17), and species level (n = 8). GSMMs of the 38 MAGs were then constructed. From the GSMMs, individual and community capabilities respectively were predicted to be able to produce 111 metabolites and 598 metabolites. Twenty-three predicted metabolites were consistent with the metabonomics detected flavors and served as targets. Twelve sub-community of were screened by cross-feeding of 38 GSMMs. Of them, Methanobacterium, Sphaerochaeta, Muricomes intestini, Methanobacteriaceae, Synergistaceae, and Caloramator were core microorganisms for targets in each sub-community. Overall, this study of metagenomic and target-community screening could help our understanding of the metabolite-microbiome association and further bioregulation of baijiu.

Genome-resolved metagenomics: a game changer for microbiome medicine

Recent substantial evidence implicating commensal bacteria in human diseases has given rise to a new domain in biomedical research: microbiome medicine. This emerging field aims to understand and leverage the human microbiota and derivative molecules for disease prevention and treatment. Despite the complex and hierarchical organization of this ecosystem, most research over the years has relied on 16S amplicon sequencing, a legacy of bacterial phylogeny and taxonomy. Although advanced sequencing technologies have enabled cost-effective analysis of entire microbiota, translating the relatively short nucleotide information into the functional and taxonomic organization of the microbiome has posed challenges until recently. In the last decade, genome-resolved metagenomics, which aims to reconstruct microbial genomes directly from whole-metagenome sequencing data, has made significant strides and continues to unveil the mysteries of various human-associated microbial communities. There has been a rapid increase in the volume of whole metagenome sequencing data and in the compilation of novel metagenome-assembled genomes and protein sequences in public depositories. This review provides an overview of the capabilities and methods of genome-resolved metagenomics for studying the human microbiome, with a focus on investigating the prokaryotic microbiota of the human gut. Just as decoding the human genome and its variations marked the beginning of the genomic medicine era, unraveling the genomes of commensal microbes and their sequence variations is ushering us into the era of microbiome medicine. Genome-resolved metagenomics stands as a pivotal tool in this transition and can accelerate our journey toward achieving these scientific and medical milestones.

Biobased short chain fatty acid production - Exploring microbial community dynamics and metabolic networks through kinetic and microbial modeling approaches

In recent years, there has been growing interest in harnessing anaerobic digestion technology for resource recovery from waste streams. This approach has evolved beyond its traditional role in energy generation to encompass the production of valuable carboxylic acids, especially volatile fatty acids (VFAs) like acetic acid, propionic acid, and butyric acid. VFAs hold great potential for various industries and biobased applications due to their versatile properties. Despite increasing global demand, over 90% of VFAs are currently produced synthetically from petrochemicals. Realizing the potential of large-scale biobased VFA production from waste streams offers significant eco-friendly opportunities but comes with several key challenges. These include low VFA production yields, unstable acid compositions, complex and expensive purification methods, and post-processing needs. Among these, production yield and acid composition stand out as the most critical obstacles impacting economic viability and competitiveness. This paper seeks to offer a comprehensive view of combining complementary modeling approaches, including kinetic and microbial modeling, to understand the workings of microbial communities and metabolic pathways in VFA production, enhance production efficiency, and regulate acid profiles through the integration of omics and bioreactor data.

A long-term growth stable Halomonas sp. deleted with multiple transposases guided by its metabolic network model Halo-ecGEM

Microbial instability is a common problem during bio-production based on microbial hosts. Halomonas bluephagenesis has been developed as a chassis for next generation industrial biotechnology (NGIB) under open and unsterile conditions. However, the hidden genomic information and peculiar metabolism have significantly hampered its deep exploitation for cell-factory engineering. Based on the freshly completed genome sequence of H. bluephagenesis TD01, which reveals 1889 biological process-associated genes grouped into 84 GO-slim terms. An enzyme constrained genome-scale metabolic model Halo-ecGEM was constructed, which showed strong ability to simulate fed-batch fermentations. A visible salt-stress responsive landscape was achieved by combining GO-slim term enrichment and CVT-based omics profiling, demonstrating that cells deploy most of the protein resources by force to support the essential activity of translation and protein metabolism when exposed to salt stress. Under the guidance of Halo-ecGEM, eight transposases were deleted, leading to a significantly enhanced stability for its growth and bioproduction of various polyhydroxyalkanoates (PHA) including 3-hydroxybutyrate (3HB) homopolymer PHB, 3HB and 3-hydroxyvalerate (3HV) copolymer PHBV, as well as 3HB and 4-hydroxyvalerate (4HB) copolymer P34HB. This study sheds new light on the metabolic characteristics and stress-response landscape of H. bluephagenesis, achieving for the first time to construct a long-term growth stable chassis for industrial applications. For the first time, it was demonstrated that genome encoded transposons are the reason for microbial instability during growth in flasks and fermentors.

Emerging tools and best practices for studying gut microbial community metabolism

The human gut microbiome vastly extends the set of metabolic reactions catalysed by our own cells, with far-reaching consequences for host health and disease. However, our knowledge of gut microbial metabolism relies on a handful of model organisms, limiting our ability to interpret and predict the metabolism of complex microbial communities. In this Perspective, we discuss emerging tools for analysing and modelling the metabolism of gut microorganisms and for linking microorganisms, pathways and metabolites at the ecosystem level, highlighting promising best practices for researchers. Continued progress in this area will also require infrastructure development to facilitate cross-disciplinary synthesis of scientific findings. Collectively, these efforts can enable a broader and deeper understanding of the workings of the gut ecosystem and open new possibilities for microbiome manipulation and therapy.

hopsy - a methods marketplace for convex polytope sampling in Python

SUMMARY

Effective collaboration between developers of Bayesian inference methods and users is key to advance our quantitative understanding of biosystems. We here present hopsy, a versatile open-source platform designed to provide convenient access to powerful Markov chain Monte Carlo sampling algorithms tailored to models defined on convex polytopes (CP). Based on the high-performance C++ sampling library HOPS, hopsy inherits its strengths and extends its functionalities with the accessibility of the Python programming language. A versatile plugin-mechanism enables seamless integration with domain-specific models, providing method developers with a framework for testing, benchmarking, and distributing CP samplers to approach real-world inference tasks. We showcase hopsy by solving common and newly composed domain-specific sampling problems, highlighting important design choices. By likening hopsy to a marketplace, we emphasize its role in bringing together users and developers, where users get access to state-of-the-art methods, and developers contribute their own innovative solutions for challenging domain-specific inference problems.

AVAILABILITY AND IMPLEMENTATION

Sources, documentation and a continuously updated list of sampling algorithms are available at https://jugit.fz-juelich.de/IBG-1/ModSim/hopsy, with Linux, Windows and MacOS binaries at https://pypi.org/project/hopsy/.

Identification of metabolites reproducibly associated with Parkinson's Disease via meta-analysis and computational modelling

Many studies have reported metabolomic analysis of different bio-specimens from Parkinson's disease (PD) patients. However, inconsistencies in reported metabolite concentration changes make it difficult to draw conclusions as to the role of metabolism in the occurrence or development of Parkinson's disease. We reviewed the literature on metabolomic analysis of PD patients. From 74 studies that passed quality control metrics, 928 metabolites were identified with significant changes in PD patients, but only 190 were replicated with the same changes in more than one study. Of these metabolites, 60 exclusively increased, such as 3-methoxytyrosine and glycine, 54 exclusively decreased, such as pantothenic acid and caffeine, and 76 inconsistently changed in concentration in PD versus control subjects, such as ornithine and tyrosine. A genome-scale metabolic model of PD and corresponding metabolic map linking most of the replicated metabolites enabled a better understanding of the dysfunctional pathways of PD and the prediction of additional potential metabolic markers from pathways with consistent metabolite changes to target in future studies.

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective

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.

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Integrating Experimental and Computational Analyses of Yeast Protein Profiles for Optimizing the Production of High-Quality Microbial Proteins

Microbial proteins represent a promising solution to address the escalating global demand for protein, particularly in regions with limited arable land. Yeasts, such as Saccharomyces cerevisiae, are robust and safe protein-producing strains. However, the utilization of non-conventional yeast strains for microbial protein production has been hindered, partly due to a lack of comprehensive understanding of protein production traits. In this study, we conducted experimental analyses focusing on the growth, protein content, and amino acid composition of nine yeast strains, including one S. cerevisiae strain, three Yarrowia lipolytica strains, and five Pichia spp. strains. We identified that, though Y. lipolytica and Pichia spp. strains consumed glucose at a slower rate compared to S. cerevisiae, Pichia spp. strains showed a higher cellular protein content, and Y. lipolytica strains showed a higher glucose-to-biomass/protein yield and methionine content. We further applied computational approaches to explain that metabolism economy was the main underlying factor for the limited amount of scarce/carbon-inefficient amino acids (such as methionine) within yeast cell proteins. We additionally verified that the specialized metabolism was a key reason for the high methionine content in Y. lipolytica strains, and proposed Y. lipolytica strain as a potential producer of high-quality single-cell protein rich in scarce amino acids. Through experimental evaluation, we identified Pichia jadinii CICC 1258 as a potential strain for high-quality protein production under unfavorable pH/temperature conditions. Our work suggests a promising avenue for optimizing microbial protein production, identifying the factors influencing amino acid composition, and paving the way for the use of unconventional yeast strains to meet the growing protein demands.

Accelerating Genetic Sensor Development, Scale-up, and Deployment Using Synthetic Biology

Living cells are exquisitely tuned to sense and respond to changes in their environment. Repurposing these systems to create engineered biosensors has seen growing interest in the field of synthetic biology and provides a foundation for many innovative applications spanning environmental monitoring to improved biobased production. In this review, we present a detailed overview of currently available biosensors and the methods that have supported their development, scale-up, and deployment. We focus on genetic sensors in living cells whose outputs affect gene expression. We find that emerging high-throughput experimental assays and evolutionary approaches combined with advanced bioinformatics and machine learning are establishing pipelines to produce genetic sensors for virtually any small molecule, protein, or nucleic acid. However, more complex sensing tasks based on classifying compositions of many stimuli and the reliable deployment of these systems into real-world settings remain challenges. We suggest that recent advances in our ability to precisely modify nonmodel organisms and the integration of proven control engineering principles (e.g., feedback) into the broader design of genetic sensing systems will be necessary to overcome these hurdles and realize the immense potential of the field.

Machine-learning-based integrative -'omics analyses reveal immunologic and metabolic dysregulation in environmental enteric dysfunction

Environmental enteric dysfunction (EED) is a subclinical enteropathy challenging to diagnose due to an overlap of tissue features with other inflammatory enteropathies. EED subjects (n = 52) from Pakistan, controls (n = 25), and a validation EED cohort (n = 30) from Zambia were used to develop a machine-learning-based image analysis classification model. We extracted histologic feature representations from the Pakistan EED model and correlated them to transcriptomics and clinical biomarkers. In-silico metabolic network modeling was used to characterize alterations in metabolic flux between EED and controls and validated using untargeted lipidomics. Genes encoding beta-ureidopropionase, CYP4F3, and epoxide hydrolase 1 correlated to numerous tissue feature representations. Fatty acid and glycerophospholipid metabolism-related reactions showed altered flux. Increased phosphatidylcholine, lysophosphatidylcholine (LPC), and ether-linked LPCs, and decreased ester-linked LPCs were observed in the duodenal lipidome of Pakistan EED subjects, while plasma levels of glycine-conjugated bile acids were significantly increased. Together, these findings elucidate a multi-omic signature of EED.

Pathogenic Bacteroides fragilis strains can emerge from gut-resident commensals

Bacteroides fragilis is a prominent member of the human gut microbiota, playing crucial roles in maintaining gut homeostasis and host health. Although it primarily functions as a beneficial commensal, B. fragilis can become pathogenic. To determine the genetic basis of its duality, we conducted a comparative genomic analysis of 813 B. fragilis strains, representing both commensal and pathogenic origins. Our findings reveal that pathogenic strains emerge across diverse phylogenetic lineages, due in part to rapid gene exchange and the adaptability of the accessory genome. We identified 16 phylogenetic groups, differentiated by genes associated with capsule composition, interspecies competition, and host interactions. A microbial genome-wide association study identified 44 genes linked to extra-intestinal survival and pathogenicity. These findings reveal how genomic diversity within commensal species can lead to the emergence of pathogenic traits, broadening our understanding of microbial evolution in the gut.

Microbiome modeling: a beginner's guide

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.

gMCSpy: efficient and accurate computation of genetic minimal cut sets in Python

MOTIVATION

The identification of minimal genetic interventions that modulate metabolic processes constitutes one of the most relevant applications of genome-scale metabolic models (GEMs). The concept of Minimal Cut Sets (MCSs) and its extension at the gene level, genetic Minimal Cut Sets (gMCSs), have attracted increasing interest in the field of Systems Biology to address this task. Different computational tools have been developed to calculate MCSs and gMCSs using both commercial and open-source software.

RESULTS

Here, we present gMCSpy, an efficient Python package to calculate gMCSs in GEMs using both commercial and non-commercial optimization solvers. We show that gMCSpy substantially overperforms our previous computational tool GMCS, which exclusively relied on commercial software. Moreover, we compared gMCSpy with recently published competing algorithms in the literature, finding significant improvements in both accuracy and computation time. All these advances make gMCSpy an attractive tool for researchers in the field of Systems Biology for different applications in health and biotechnology.

AVAILABILITY AND IMPLEMENTATION

The Python package gMCSpy and the data underlying this manuscript can be accessed at: https://github.com/PlanesLab/gMCSpy.

Engineering natural microbiomes toward enhanced bioremediation by microbiome modeling

Engineering natural microbiomes for biotechnological applications remains challenging, as metabolic interactions within microbiomes are largely unknown, and practical principles and tools for microbiome engineering are still lacking. Here, we present a combinatory top-down and bottom-up framework to engineer natural microbiomes for the construction of function-enhanced synthetic microbiomes. We show that application of herbicide and herbicide-degrader inoculation drives a convergent succession of different natural microbiomes toward functional microbiomes (e.g., enhanced bioremediation of herbicide-contaminated soils). We develop a metabolic modeling pipeline, SuperCC, that can be used to document metabolic interactions within microbiomes and to simulate the performances of different microbiomes. Using SuperCC, we construct bioremediation-enhanced synthetic microbiomes based on 18 keystone species identified from natural microbiomes. Our results highlight the importance of metabolic interactions in shaping microbiome functions and provide practical guidance for engineering natural microbiomes.

Bridging systems biology and tissue engineering: Unleashing the full potential of complex 3D in vitro tissue models of disease

Rapid advances in tissue engineering have resulted in more complex and physiologically relevant 3D in vitro tissue models with applications in fundamental biology and therapeutic development. However, the complexity provided by these models is often not leveraged fully due to the reductionist methods used to analyze them. Computational and mathematical models developed in the field of systems biology can address this issue. Yet, traditional systems biology has been mostly applied to simpler in vitro models with little physiological relevance and limited cellular complexity. Therefore, integrating these two inherently interdisciplinary fields can result in new insights and move both disciplines forward. In this review, we provide a systematic overview of how systems biology has been integrated with 3D in vitro tissue models and discuss key application areas where the synergies between both fields have led to important advances with potential translational impact. We then outline key directions for future research and discuss a framework for further integration between fields.

Rewiring Saccharomyces cerevisiae metabolism for optimised Taxol® precursors production

Saccharomyces cerevisiae has been conveniently used to produce Taxol® anticancer drug early precursors. However, the harmful impact of oxidative stress by the first cytochrome P450-reductase enzymes (CYP725A4-POR) of Taxol® pathway has hampered sufficient progress in yeast. Here, we evolved an oxidative stress-resistant yeast strain with three-fold higher titre of their substrate, taxadiene. The performance of the evolved and parent strains were then evaluated in galactose-limited chemostats before and under the oxidative stress by an oxidising agent. The interaction of evolution and oxidative stress was comprehensively evaluated through transcriptomics and metabolite profiles integration in yeast enzyme-constrained genome scale model. Overall, the evolved strain showed improved respiration, reduced overflow metabolites production and oxidative stress re-induction tolerance. The cross-protection mechanism also potentially contributed to better heme, flavin and NADPH availability, essential for CYP725A4 and POR optimal activity in yeast. The results imply that the evolved strain is a robust cell factory for future efforts towards Taxol© production.

Emerging methods for genome-scale metabolic modeling of microbial communities

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.

GEM-based computational modeling for exploring metabolic interactions in a microbial community

Microbial communities play fundamental roles in every complex ecosystem, such as soil, sea and the human body. The stability and diversity of the microbial community depend precisely on the composition of the microbiota. Any change in the composition of these communities affects microbial functions. An important goal of studying the interactions between species is to understand the behavior of microbes and their responses to perturbations. These interactions among species are mediated by the exchange of metabolites within microbial communities. We developed a computational model for the microbial community that has a separate compartment for exchanging metabolites. This model can predict possible metabolites that cause competition, commensalism, and mutual interactions between species within a microbial community. Our constraint-based community metabolic modeling approach provides insights to elucidate the pattern of metabolic interactions for each common metabolite between two microbes. To validate our approach, we used a toy model and a syntrophic co-culture of Desulfovibrio vulgaris and Methanococcus maripaludis, as well as another in co-culture between Geobacter sulfurreducens and Rhodoferax ferrireducens. For a more general evaluation, we applied our algorithm to the honeybee gut microbiome, composed of seven species, and the epiphyte strain Pantoea eucalypti 299R. The epiphyte strain Pe299R has been previously studied and cultured with six different phyllosphere bacteria. Our algorithm successfully predicts metabolites, which imply mutualistic, competitive, or commensal interactions. In contrast to OptCom, MRO, and MICOM algorithms, our COMMA algorithm shows that the potential for competitive interactions between an epiphytic species and Pe299R is not significant. These results are consistent with the experimental measurements of population density and reproductive success of the Pe299R strain.

Integrating Genome-Scale Metabolic Models with Patient Plasma Metabolome to Study Endothelial Metabolism In Situ

Patient blood samples are invaluable in clinical omics databases, yet current methodologies often fail to fully uncover the molecular mechanisms driving patient pathology. While genome-scale metabolic models (GEMs) show promise in systems medicine by integrating various omics data, having only exometabolomic data remains a limiting factor. To address this gap, we introduce a comprehensive pipeline integrating GEMs with patient plasma metabolome. This pipeline constructs case-specific GEMs using literature-based and patient-specific metabolomic data. Novel computational methods, including adaptive sampling and an in-house developed algorithm for the rational exploration of the sampled space of solutions, enhance integration accuracy while improving computational performance. Model characterization involves task analysis in combination with clustering methods to identify critical cellular functions. The new pipeline was applied to a cohort of trauma patients to investigate shock-induced endotheliopathy using patient plasma metabolome data. By analyzing endothelial cell metabolism comprehensively, the pipeline identified critical therapeutic targets and biomarkers that can potentially contribute to the development of therapeutic strategies. Our study demonstrates the efficacy of integrating patient plasma metabolome data into computational models to analyze endothelial cell metabolism in disease contexts. This approach offers a deeper understanding of metabolic dysregulations and provides insights into diseases with metabolic components and potential treatments.

Microbial Abundances Retrieved from Sequencing data-automated NCBI Taxonomy (MARS): a pipeline to create relative microbial abundance data for the Microbiome Modelling Toolbox and utilizing homosynonyms for efficient mapping to resources

Motivation

Computational approaches to the functional characterization of the microbiome, such as the Microbiome Modelling Toolbox, require precise information on microbial composition and relative abundances. However, challenges arise from homosynonyms-different names referring to the same taxon, which can hinder the mapping process and lead to missed species mapping when using microbial metabolic reconstruction resources, such as AGORA and APOLLO.

Results

We introduce the integrated MARS pipeline, a user-friendly Python-based solution that addresses these challenges. MARS automates the extraction of relative abundances from metagenomic reads, maps species and genera onto microbial metabolic reconstructions, and accounts for alternative taxonomic names. It normalizes microbial reads, provides an optional cut-off for low-abundance taxa, and produces relative abundance tables apt for integration with the Microbiome Modelling Toolbox. A sub-component of the pipeline automates the task of identifying homosynonyms, leveraging web scraping to find taxonomic IDs of given species, searching NCBI for alternative names, and cross-reference them with microbial reconstruction resources. Taken together, MARS streamlines the entire process from processed metagenomic reads to relative abundance, thereby significantly reducing time and effort when working with microbiome data.

Availability and implementation

MARS is implemented in Python. It can be found as an interactive application here: https://mars-pipeline.streamlit.app/along with a detailed documentation here: https://github.com/ThieleLab/mars-pipeline.

Designing function-specific minimal microbiomes from large microbial communities

Microorganisms exist in large communities of diverse species, exhibiting various functionalities. The mammalian gut microbiome, for instance, has the functionality of digesting dietary fibre and producing different short-chain fatty acids. Not all microbes present in a community contribute to a given functionality; it is possible to find a minimal microbiome, which is a subset of the large microbiome, that is capable of performing the functionality while maintaining other community properties such as growth rate and metabolite production. Such a minimal microbiome will also contain keystone species for SCFA production in that community. In this work, we present a systematic constraint-based approach to identify a minimal microbiome from a large community for a user-proposed function. We employ a top-down approach with sequential deletion followed by solving a mixed-integer linear programming problem with the objective of minimising the L1-norm of the membership vector. Notably, we consider quantitative measures of community growth rate and metabolite production rates. We demonstrate the utility of our algorithm by identifying the minimal microbiomes corresponding to three model communities of the gut, and discuss their validity based on the presence of the keystone species in the community. Our approach is generic, flexible and finds application in studying a variety of microbial communities. The algorithm is available from https://github.com/RamanLab/minMicrobiome .

BN-BacArena: Bayesian network extension of BacArena for the dynamic simulation of microbial communities

MOTIVATION

Simulating gut microbial dynamics is extremely challenging. Several computational tools, notably the widely used BacArena, enable modeling of dynamic changes in the microbial environment. These methods, however, do not comprehensively account for microbe-microbe stimulant or inhibitory effects or for nutrient-microbe inhibitory effects, typically observed in different compounds present in the daily diet.

RESULTS

Here, we present BN-BacArena, an extension of BacArena consisting on the incorporation within the native computational framework of a Bayesian network model that accounts for microbe-microbe and nutrient-microbe interactions. Using in vitro experiments, 16S rRNA gene sequencing data and nutritional composition of 55 foods, the output Bayesian network showed 23 significant nutrient-bacteria interactions, suggesting the importance of compounds such as polyols, ascorbic acid, polyphenols and other phytochemicals, and 40 bacteria-bacteria significant relationships. With test data, BN-BacArena demonstrates a statistically significant improvement over BacArena to predict the time-dependent relative abundance of bacterial species involved in the gut microbiota upon different nutritional interventions. As a result, BN-BacArena opens new avenues for the dynamic modeling and simulation of the human gut microbiota metabolism.

AVAILABILITY AND IMPLEMENTATION

MATLAB and R code are available in https://github.com/PlanesLab/BN-BacArena.

Hepatic glucose metabolism in the steatotic liver

The liver is central in regulating glucose homeostasis, being the major contributor to endogenous glucose production and the greatest reserve of glucose as glycogen. It is both a target and regulator of the action of glucoregulatory hormones. Hepatic metabolic functions are altered in and contribute to the highly prevalent steatotic liver disease (SLD), including metabolic dysfunction-associated SLD (MASLD) and metabolic dysfunction-associated steatohepatitis (MASH). In this Review, we describe the dysregulation of hepatic glucose metabolism in MASLD and MASH and associated metabolic comorbidities, and how advances in techniques and models for the assessment of hepatic glucose fluxes in vivo have led to the identification of the mechanisms related to the alterations in glucose metabolism in MASLD and comorbidities. These fluxes can ultimately increase hepatic glucose production concomitantly with fat accumulation and alterations in the secretion and action of glucoregulatory hormones. No pharmacological treatment has yet been approved for MASLD or MASH, but some antihyperglycaemic drugs approved for treating type 2 diabetes have shown positive effects on hepatic glucose metabolism and hepatosteatosis. A deep understanding of how MASLD affects glucose metabolic fluxes and glucoregulatory hormones might assist in the early identification of at-risk individuals and the use or development of targeted therapies.

Genome-scale metabolic network models for industrial microorganisms metabolic engineering: Current advances and future prospects

The construction of high-performance microbial cell factories (MCFs) is the centerpiece of biomanufacturing. However, the complex metabolic regulatory network of microorganisms poses great challenges for the efficient design and construction of MCFs. The genome-scale metabolic network models (GSMs) can systematically simulate the metabolic regulation process of microorganisms in silico, providing effective guidance for the rapid design and construction of MCFs. In this review, we summarized the development status of 16 important industrial microbial GSMs, and further outline the technologies or methods that continuously promote high-quality GSMs construction from five aspects: I) Databases and modeling tools facilitate GSMs reconstruction; II) evolving gap-filling technologies; III) constraint-based model reconstruction; IV) advances in algorithms; and V) developed visualization tools. In addition, we also summarized the applications of GSMs in guiding metabolic engineering from four aspects: I) exploring and explaining metabolic features; II) predicting the effects of genetic perturbations on metabolism; III) predicting the optimal phenotype; IV) guiding cell factories construction in practical experiment. Finally, we discussed the development of GSMs, aiming to provide a reference for efficiently reconstructing GSMs and guiding metabolic engineering.

A model-driven approach to upcycling recalcitrant feedstocks in Pseudomonas putida by decoupling PHA production from nutrient limitation

Bacterial polyhydroxyalkanoates (PHAs) have emerged as promising eco-friendly alternatives to petroleum-based plastics since they are synthesized from renewable resources and offer exceptional properties. However, their production is limited to the stationary growth phase under nutrient-limited conditions, requiring customized strategies and costly two-phase bioprocesses. In this study, we tackle these challenges by employing a model-driven approach to reroute carbon flux and remove regulatory constraints using synthetic biology. We construct a collection of Pseudomonas putida-overproducing strains at the expense of plastics and lignin-related compounds using growth-coupling approaches. PHA production was successfully achieved during growth phase, resulting in the production of up to 46% PHA/cell dry weight while maintaining a balanced carbon-to-nitrogen ratio. Our strains are additionally validated under an upcycling scenario using enzymatically hydrolyzed polyethylene terephthalate as a feedstock. These findings have the potential to revolutionize PHA production and address the global plastic crisis by overcoming the complexities of traditional PHA production bioprocesses.

Design of a proteolytic module for improved metabolic modeling of Bacteroides caccae

The gut microbiota plays a crucial role in health and is significantly modulated by human diets. In addition to Western diets which are rich in proteins, high-protein diets are used for specific populations or indications, mainly weight loss. In this study, we investigated the effect of protein supplementation on Bacteroides caccae, a Gram-negative gut symbiont. The supplementation with whey proteins led to a significant increase in growth rate, final biomass, and short-chain fatty acids production. A comprehensive genomic analysis revealed that B. caccae possesses a set of 156 proteases with putative intracellular and extracellular localization and allowed to identify amino acid transporters and metabolic pathways. We developed a fully curated genome-scale metabolic model of B. caccae that incorporated its proteolytic activity and simulated its growth and production of fermentation-related metabolites in response to the different growth media. We validated the model by comparing the predicted phenotype to experimental data. The model accurately predicted B. caccae's growth and metabolite production (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.