Metabolic trade-offs constrain the cell size ratio in a nitrogen-fixing symbiosis.

Biological dinitrogen (N2) fixation is a key metabolic process exclusively performed by prokaryotes, some of which are symbiotic with eukaryotes. Species of the marine haptophyte algae Braarudosphaera bigelowii harbor the N2-fixing endosymbiotic cyanobacteria UCYN-A, which might be evolving organelle-like characteristics. We found that the size ratio between UCYN-A and their hosts is strikingly conserved across sublineages/species, which is consistent with the size relationships of organelles in this symbiosis and other species. Metabolic modeling showed that this size relationship maximizes the coordinated growth rate based on trade-offs between resource acquisition and exchange. Our findings show that the size relationships of N2-fixing endosymbionts and organelles in unicellular eukaryotes are constrained by predictable metabolic underpinnings and that UCYN-A is, in many regards, functioning like a hypothetical N2-fixing organelle (or nitroplast).

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

Host-Microbe-Drug-Nutrient Screen Identifies Bacterial Effectors of Metformin Therapy.

Metformin is the first-line therapy for treating type 2 diabetes and a promising anti-aging drug. We set out to address the fundamental question of how gut microbes and nutrition, key regulators of host physiology, affect the effects of metformin. Combining two tractable genetic models, the bacterium E. coli and the nematode C. elegans, we developed a high-throughput four-way screen to define the underlying host-microbe-drug-nutrient interactions. We show that microbes integrate cues from metformin and the diet through the phosphotransferase signaling pathway that converges on the transcriptional regulator Crp. A detailed experimental characterization of metformin effects downstream of Crp in combination with metabolic modeling of the microbiota in metformin-treated type 2 diabetic patients predicts the production of microbial agmatine, a regulator of metformin effects on host lipid metabolism and lifespan. Our high-throughput screening platform paves the way for identifying exploitable drug-nutrient-microbiome interactions to improve host health and longevity through targeted microbiome therapies. VIDEO ABSTRACT.

Interpreting metabolic complexity via isotope-assisted metabolic flux analysis.

Isotope-assisted metabolic flux analysis (iMFA) is a powerful method to mathematically determine the metabolic fluxome from experimental isotope labeling data and a metabolic network model. While iMFA was originally developed for industrial biotechnological applications, it is increasingly used to analyze eukaryotic cell metabolism in physiological and pathological states. In this review, we explain how iMFA estimates the intracellular fluxome, including data and network model (inputs), the optimization-based data fitting (process), and the flux map (output). We then describe how iMFA enables analysis of metabolic complexities and discovery of metabolic pathways. Our goal is to expand the use of iMFA in metabolism research, which is essential to maximizing the impact of metabolic experiments and continuing to advance iMFA and biocomputational techniques.

Modeling blood metabolite homeostatic levels reduces sample heterogeneity across cohorts.

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

Metabolic models predict fotemustine and the combination of eflornithine/rifamycin and adapalene/cannabidiol for the treatment of gliomas.

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

Genome-Scale Metabolic Modeling of the Human Microbiome in the Era of Personalized Medicine.

The human microbiome plays an important role in human health and disease. Meta-omics analyses provide indispensable data for linking changes in microbiome composition and function to disease etiology. Yet, the lack of a mechanistic understanding of, e.g., microbiome-metabolome links hampers the translation of these findings into effective, novel therapeutics. Here, we propose metabolic modeling of microbial communities through constraint-based reconstruction and analysis (COBRA) as a complementary approach to meta-omics analyses. First, we highlight the importance of microbial metabolism in cardiometabolic diseases, inflammatory bowel disease, colorectal cancer, Alzheimer disease, and Parkinson disease. Next, we demonstrate that microbial community modeling can stratify patients and controls, mechanistically link microbes with fecal metabolites altered in disease, and identify host pathways affected by the microbiome. Finally, we outline our vision for COBRA modeling combined with meta-omics analyses and multivariate statistical analyses to inform and guide clinical trials, yield testable hypotheses, and ultimately propose novel dietary and therapeutic interventions.

Kinetic modeling identifies targets for engineering improved photosynthetic efficiency in potato (Solanum tuberosum cv. Solara).

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

DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates.

The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.

Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Engineering of Saccharomyces cerevisiae for enhanced metabolic robustness and L-lactic acid production from lignocellulosic biomass.

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

Revealing the dynamics and mechanisms of bacterial interactions in cheese production with metabolic modelling.

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

Soybean-mediated suppression of BjaI/BjaR1 quorum sensing in Bradyrhizobium diazoefficiens impacts symbiotic nitrogen fixation.

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

Vibrio natriegens genome-scale modeling reveals insights into halophilic adaptations and resource allocation.

Vibrio natriegens is a Gram-negative bacterium with an exceptional growth rate that has the potential to become a standard biotechnological host for laboratory and industrial bioproduction. Despite this burgeoning interest, the current lack of organism-specific qualitative and quantitative computational tools has hampered the community's ability to rationally engineer this bacterium. In this study, we present the first genome-scale metabolic model (GSMM) of V. natriegens. The GSMM (iLC858) was developed using an automated draft assembly and extensive manual curation and was validated by comparing predicted yields, central metabolic fluxes, viable carbon substrates, and essential genes with empirical data. Mass spectrometry-based proteomics data confirmed the translation of at least 76% of the enzyme-encoding genes predicted to be expressed by the model during aerobic growth in a minimal medium. iLC858 was subsequently used to carry out a metabolic comparison between the model organism Escherichia coli and V. natriegens, leading to an analysis of the model architecture of V. natriegens' respiratory and ATP-generating system and the discovery of a role for a sodium-dependent oxaloacetate decarboxylase pump. The proteomics data were further used to investigate additional halophilic adaptations of V. natriegens. Finally, iLC858 was utilized to create a Resource Balance Analysis model to study the allocation of carbon resources. Taken together, the models presented provide useful computational tools to guide metabolic engineering efforts in V. natriegens.

Isotopically nonstationary metabolic flux analysis of plants: recent progress and future opportunities.

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

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities.

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

Computer-assisted biosynthetic capacities elucidation accelerates our ability to interpret microbial interactions, allowing us to understand better and establish a balance within ecosystems.

Metabolic and mitochondria alterations induced by SARS-CoV-2 accessory proteins ORF3a, ORF9b, ORF9c and ORF10.

Antiviral signaling, immune response and cell metabolism are dysregulated by SARS-CoV-2, the causative agent of COVID-19. Here, we show that SARS-CoV-2 accessory proteins ORF3a, ORF9b, ORF9c and ORF10 induce a significant mitochondrial and metabolic reprogramming in A549 lung epithelial cells. While ORF9b, ORF9c and ORF10 induced largely overlapping transcriptomes, ORF3a induced a distinct transcriptome, including the downregulation of numerous genes with critical roles in mitochondrial function and morphology. On the other hand, all four ORFs altered mitochondrial dynamics and function, but only ORF3a and ORF9c induced a marked alteration in mitochondrial cristae structure. Genome-Scale Metabolic Models identified both metabolic flux reprogramming features both shared across all accessory proteins and specific for each accessory protein. Notably, a downregulated amino acid metabolism was observed in ORF9b, ORF9c and ORF10, while an upregulated lipid metabolism was distinctly induced by ORF3a. These findings reveal metabolic dependencies and vulnerabilities prompted by SARS-CoV-2 accessory proteins that may be exploited to identify new targets for intervention.

Development of generic metabolic Raman calibration models using solution titration in aqueous phase and data augmentation for in-line cell culture analysis.

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

Itaconic acid production by co-feeding of Ustilago maydis: A combined approach of experimental data, design of experiments, and metabolic modeling.

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

Guidelines for extracting biologically relevant context-specific metabolic models using gene expression data.

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