Path to improving the life cycle and quality of genome-scale models of metabolism

Scaling metabolic model reconstruction up to the pan-genome level: A systematic review and prospective applications to photosynthetic organisms

Advances in genomics technologies have generated large data sets that provide tremendous insights into the genetic diversity of taxonomic groups. However, it remains challenging to pinpoint the effect of genetic diversity on different traits without performing resource-intensive phenotyping experiments. Pan-genome-scale metabolic models (panGEMs) extend traditional genome-scale metabolic models by considering the entire reaction repertoire that enables the prediction and comparison of metabolic capabilities within a taxonomic group. Here, we systematically review the state-of-the-art methodologies for constructing panGEMs, focusing on used tools, databases, experimental datasets, and orthology relationships. We highlight the unique advantages of panGEMs compared to single-species GEMs in predicting metabolic phenotypes and in guiding the experimental validation of genome annotations. In addition, we emphasize the disparity between the available (pan-)genomic data on photosynthetic organisms and their under-representation in current (pan)GEMs. Finally, we propose a perspective for tackling the reconstruction of panGEMs for photosynthetic eukaryotes that can help advance our understanding of the metabolic diversity in this taxonomic group.

Influence of Amino Acids on Quorum Sensing-Related Pathways in Pseudomonas aeruginosa PAO1: Insights from the GEM iJD1249

BACKGROUND/OBJECTIVES

Amino acids (AAs) play a critical role in diseases such as cystic fibrosis where Pseudomonas aeruginosa PAO1 adapts its metabolism in response to host-derived nutrients. The adaptation influences virulence and complicates antibiotic treatment mainly for the antimicrobial resistance context. D- and L-AAs have been analyzed for their impact on quorum sensing (QS), a mechanism that regulates virulence factors. This research aimed to reconstruct the genome-scale metabolic model (GEM) of P. aeruginosa PAO1 to investigate the metabolic roles of D- and L-AAs in QS-related pathways.

METHODS

The updated GEM, iJD1249, was reconstructed by using protocols to integrate data from previous models and refined with well-standardized in silico media (LB, M9, and SCFM) to improve flux balance analysis accuracy. The model was used to explore the metabolic impact of D-Met, D-Ala, D-Glu, D-Ser, L-His, L-Glu, L-Arg, and L-Ornithine (L-Orn) at 5 and 50 mM in QS-related pathways, focusing on the effects on bacterial growth and carbon flux distributions.

RESULTS

Among the tested AAs, D-Met was the only one that did not enhance the growth rate of P. aeruginosa PAO1, while L-Arg and L-Orn increased fluxes in the L-methionine biosynthesis pathway, influencing the metH gene. These findings suggest a differential metabolic role for D-and L-AAs in QS-related pathways.

CONCLUSIONS

Our results shed some light on the metabolic impact of AAs on QS-related pathways and their potential role in P. aeruginosa virulence. Future studies should assess D-Met as a potential adjuvant in antimicrobial strategies, optimizing the concentration in combination with antibiotics to maximize its therapeutic effectiveness.

Metabolic Flux Modeling in Marine Ecosystems

Ocean metabolism constitutes a complex, multiscale ensemble of biochemical reaction networks harbored within and between the boundaries of a myriad of organisms. Gaining a quantitative understanding of how these networks operate requires mathematical tools capable of solving in silico the resource allocation problem each cell faces in real life. Toward this goal, stoichiometric modeling of metabolism, such as flux balance analysis, has emerged as a powerful computational tool for unraveling the intricacies of metabolic processes in microbes, microbial communities, and multicellular organisms. Here, we provide an overview of this approach and its applications, future prospects, and practical considerations in the context of marine sciences. We explore how flux balance analysis has been employed to study marine organisms, help elucidate nutrient cycling, and predict metabolic capabilities within diverse marine environments, and highlight future prospects for this field in advancing our knowledge of marine ecosystems and their sustainability.

Auxotrophy-based curation improves the consensus genome-scale metabolic model of yeast

Saccharomyces cerevisiae, a widely utilized model organism, has seen continuous updates to its genome-scale metabolic model (GEM) to enhance the prediction performance for metabolic engineering and systems biology. This study presents an auxotrophy-based curation of the yeast GEM, enabling facile upgrades to yeast GEMs in future endeavors. We illustrated that the curation bolstered the predictive capability of the yeast GEM particularly in predicting auxotrophs without compromising accuracy in other simulations, and thus could be an effective manner for GEM refinement. Last, we leveraged the curated yeast GEM to systematically predict auxotrophs, thereby furnishing a valuable reference for the design of nutrient-dependent cell factories and synthetic yeast consortia.

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. This versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

Leveraging genome-scale metabolic models to understand aerobic methanotrophs

Genome-scale metabolic models (GEMs) are valuable tools serving systems biology and metabolic engineering. However, GEMs are still an underestimated tool in informing microbial ecology. Since their first application for aerobic gammaproteobacterial methane oxidizers less than a decade ago, GEMs have substantially increased our understanding of the metabolism of methanotrophs, a microbial guild of high relevance for the natural and biotechnological mitigation of methane efflux to the atmosphere. Particularly, GEMs helped to elucidate critical metabolic and regulatory pathways of several methanotrophic strains, predicted microbial responses to environmental perturbations, and were used to model metabolic interactions in cocultures. Here, we conducted a systematic review of GEMs exploring aerobic methanotrophy, summarizing recent advances, pointing out weaknesses, and drawing out probable future uses of GEMs to improve our understanding of the ecology of methane oxidizers. We also focus on their potential to unravel causes and consequences when studying interactions of methane-oxidizing bacteria with other methanotrophs or members of microbial communities in general. This review aims to bridge the gap between applied sciences and microbial ecology research on methane oxidizers as model organisms and to provide an outlook for future studies.

Control Theory and Systems Biology: Potential Applications in Neurodegeneration and Search for Therapeutic Targets

Control theory, a well-established discipline in engineering and mathematics, has found novel applications in systems biology. This interdisciplinary approach leverages the principles of feedback control and regulation to gain insights into the complex dynamics of cellular and molecular networks underlying chronic diseases, including neurodegeneration. By modeling and analyzing these intricate systems, control theory provides a framework to understand the pathophysiology and identify potential therapeutic targets. Therefore, this review examines the most widely used control methods in conjunction with genomic-scale metabolic models in the steady state of the multi-omics type. According to our research, this approach involves integrating experimental data, mathematical modeling, and computational analyses to simulate and control complex biological systems. In this review, we find that the most significant application of this methodology is associated with cancer, leaving a lack of knowledge in neurodegenerative models. However, this methodology, mainly associated with the Minimal Dominant Set (MDS), has provided a starting point for identifying therapeutic targets for drug development and personalized treatment strategies, paving the way for more effective therapies.

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Laboratory evolution, transcriptomics, and modeling reveal mechanisms of paraquat tolerance

Relationships between the genome, transcriptome, and metabolome underlie all evolved phenotypes. However, it has proved difficult to elucidate these relationships because of the high number of variables measured. A recently developed data analytic method for characterizing the transcriptome can simplify interpretation by grouping genes into independently modulated sets (iModulons). Here, we demonstrate how iModulons reveal deep understanding of the effects of causal mutations and metabolic rewiring. We use adaptive laboratory evolution to generate E. coli strains that tolerate high levels of the redox cycling compound paraquat, which produces reactive oxygen species (ROS). We combine resequencing, iModulons, and metabolic models to elucidate six interacting stress-tolerance mechanisms: (1) modification of transport, (2) activation of ROS stress responses, (3) use of ROS-sensitive iron regulation, (4) motility, (5) broad transcriptional reallocation toward growth, and (6) metabolic rewiring to decrease NADH production. This work thus demonstrates the power of iModulon knowledge mapping for evolution analysis.

Lifespan-extending interventions induce consistent patterns of fatty acid oxidation in mouse livers

Aging manifests as progressive deteriorations in homeostasis, requiring systems-level perspectives to investigate the gradual molecular dysregulation of underlying biological processes. Here, we report systemic changes in the molecular regulation of biological processes under multiple lifespan-extending interventions. Differential Rank Conservation (DIRAC) analyses of mouse liver proteomics and transcriptomics data show that mechanistically distinct lifespan-extending interventions (acarbose, 17α-estradiol, rapamycin, and calorie restriction) generally tighten the regulation of biological modules. These tightening patterns are similar across the interventions, particularly in processes such as fatty acid oxidation, immune response, and stress response. Differences in DIRAC patterns between proteins and transcripts highlight specific modules which may be tightened via augmented cap-independent translation. Moreover, the systemic shifts in fatty acid metabolism are supported through integrated analysis of liver transcriptomics data with a mouse genome-scale metabolic model. Our findings highlight the power of systems-level approaches for identifying and characterizing the biological processes involved in aging and longevity.

Teasing out missing reactions in genome-scale metabolic networks through hypergraph learning

GEnome-scale Metabolic models (GEMs) are powerful tools to predict cellular metabolism and physiological states in living organisms. However, due to our imperfect knowledge of metabolic processes, even highly curated GEMs have knowledge gaps (e.g., missing reactions). Existing gap-filling methods typically require phenotypic data as input to tease out missing reactions. We still lack a computational method for rapid and accurate gap-filling of metabolic networks before experimental data is available. Here we present a deep learning-based method - CHEbyshev Spectral HyperlInk pREdictor (CHESHIRE) - to predict missing reactions in GEMs purely from metabolic network topology. We demonstrate that CHESHIRE outperforms other topology-based methods in predicting artificially removed reactions over 926 high- and intermediate-quality GEMs. Furthermore, CHESHIRE is able to improve the phenotypic predictions of 49 draft GEMs for fermentation products and amino acids secretions. Both types of validation suggest that CHESHIRE is a powerful tool for GEM curation to reveal unknown links between reactions and observed metabolic phenotypes.

Quantitative systems pharmacology of the eye: Tools and data for ocular QSP

Good eyesight belongs to the most-valued attributes of health, and diseases of the eye are a significant healthcare burden. Case numbers are expected to further increase in the next decades due to an aging society. The development of drugs in ophthalmology, however, is difficult due to limited accessibility of the eye, in terms of drug administration and in terms of sampling of tissues for drug pharmacokinetics (PKs) and pharmacodynamics (PDs). Ocular quantitative systems pharmacology models provide the opportunity to describe the distribution of drugs in the eye as well as the resulting drug-response in specific segments of the eye. In particular, ocular physiologically-based PK (PBPK) models are necessary to describe drug concentration levels in different regions of the eye. Further, ocular effect models using molecular data from specific cellular systems are needed to develop dose-response correlations. We here describe the current status of PK/PBPK as well as PD models for the eyes and discuss cellular systems, data repositories, as well as animal models in ophthalmology. The application of the various concepts is highlighted for the development of new treatments for postoperative fibrosis after glaucoma surgery.

Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks

Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.

Making waves: Applying systems biology principles in water distribution systems engineering

The complexity of modeling water quality variations in water distribution systems (WDS), studied for decades, stems from multiple constraints and variables involved and the complexity of the system behavior. The conventional macroscale-based WDS water quality models are founded on continuum mechanics. In attempts to provide a broad picture of the multi-species interactions, these models overlook the stochasticity corresponding to the reaction mechanisms within the WDS domain. Furthermore, owing to the black-box type modeling adopted in simulating the multi-species interactions, the existing state-of-the-art models have limitations in representing intermediates and/or by-products formation. Accordingly, they remain ineffective in describing the water chemistry-stoichiometric interactions within the WDS domain. Only a radically new modeling approach could overcome the limitations of the macroscale-based approaches and enables analyzing the stochastic WDS mechanisms by keeping the true nature of the system behavior. Stimulated by the metabolic network modeling principles in systems biology, this article outlines the prospect of developing an innovative 'water'bolic network modeling approach to provide a new outlook to the existing WDS water quality modeling research.

Recent Progress on Systems and Synthetic Biology of Diatoms for Improving Algal Productivity

Microalgae have drawn much attention for their potential applications as a sustainable source for developing bioactive compounds, functional foods, feeds, and biofuels. Diatoms, as one major group of microalgae with high yields and strong adaptability to the environment, have shown advantages in developing photosynthetic cell factories to produce value-added compounds, including heterologous bioactive products. However, the commercialization of diatoms has encountered several obstacles that limit the potential mass production, such as the limitation of algal productivity and low photosynthetic efficiency. In recent years, systems and synthetic biology have dramatically improved the efficiency of diatom cell factories. In this review, we discussed first the genome sequencing and genome-scale metabolic models (GEMs) of diatoms. Then, approaches to optimizing photosynthetic efficiency are introduced with a focus on the enhancement of biomass productivity in diatoms. We also reviewed genome engineering technologies, including CRISPR (clustered regularly interspaced short palindromic repeats) gene-editing to produce bioactive compounds in diatoms. Finally, we summarized the recent progress on the diatom cell factory for producing heterologous compounds through genome engineering to introduce foreign genes into host diatoms. This review also pinpointed the bottlenecks in algal engineering development and provided critical insights into the future direction of algal production.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.