Whole-body metabolic modelling reveals microbiome and genomic interactions on reduced urine formate levels in Alzheimer's disease

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimer's disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilised whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalised whole-body metabolic models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

Whole-body metabolic modelling predicts isoleucine dependency of SARS-CoV-2 replication

We aimed at investigating host-virus co-metabolism during SARS-CoV-2 infection. Therefore, we extended comprehensive sex-specific, whole-body organ resolved models of human metabolism with the necessary reactions to replicate SARS-CoV-2 in the lung as well as selected peripheral organs. Using this comprehensive host-virus model, we obtained the following key results: 1. The predicted maximal possible virus shedding rate was limited by isoleucine availability. 2. The supported initial viral load depended on the increase in CD4+ T-cells, consistent with the literature. 3. During viral infection, the whole-body metabolism changed including the blood metabolome, which agreed well with metabolomic studies from COVID-19 patients and healthy controls. 4. The virus shedding rate could be reduced by either inhibition of the guanylate kinase 1 or availability of amino acids, e.g., in the diet. 5. The virus variants differed in their maximal possible virus shedding rates, which could be inversely linked to isoleucine occurrences in the sequences. Taken together, this study presents the metabolic crosstalk between host and virus and emphasises the role of amino acid metabolism during SARS-CoV-2 infection, in particular of isoleucine. As such, it provides an example of how computational modelling can complement more canonical approaches to gain insight into host-virus crosstalk and to identify potential therapeutic strategies.

Two-species community design of lactic acid bacteria for optimal production of lactate

Microbial communities that metabolise pentose and hexose sugars are useful in producing high-value chemicals, resulting in the effective conversion of raw materials to the product, a reduction in the production cost, and increased yield. Here, we present a computational analysis approach called CAMP (Co-culture/Community Analyses for Metabolite Production) that simulates and identifies appropriate communities to produce a metabolite of interest. To demonstrate this approach, we focus on the optimal production of lactate from various Lactic Acid Bacteria. We used genome-scale metabolic models (GSMMs) belonging to Lactobacillus, Leuconostoc, and Pediococcus species from the Virtual Metabolic Human (VMH; https://vmh.life/) resource and well-curated GSMMs of L. plantarum WCSF1 and L. reuteri JCM 1112. We analysed 1176 two-species communities using a constraint-based modelling method for steady-state flux-balance analysis of communities. Flux variability analysis was used to detect the maximum lactate flux in the communities. Using glucose or xylose as substrates separately or in combination resulted in either parasitism, amensalism, or mutualism being the dominant interaction behaviour in the communities. Interaction behaviour between members of the community was deduced based on variations in the predicted growth rates of monocultures and co-cultures. Acetaldehyde, ethanol, acetate, among other metabolites, were found to be cross-fed between community members. L. plantarum WCSF1 was found to be a member of communities with high lactate yields. In silico community optimisation strategies to predict reaction knock-outs for improving lactate flux were implemented. Reaction knock-outs of acetate kinase, phosphate acetyltransferase, and fumarate reductase in the communities were found to enhance lactate production.

Modelling microbial communities: Harnessing consortia for biotechnological applications

Microbes propagate and thrive in complex communities, and there are many benefits to studying and engineering microbial communities instead of single strains. Microbial communities are being increasingly leveraged in biotechnological applications, as they present significant advantages such as the division of labour and improved substrate utilisation. Nevertheless, they also present some interesting challenges to surmount for the design of efficient biotechnological processes. In this review, we discuss key principles of microbial interactions, followed by a deep dive into genome-scale metabolic models, focussing on a vast repertoire of constraint-based modelling methods that enable us to characterise and understand the metabolic capabilities of microbial communities. Complementary approaches to model microbial communities, such as those based on graph theory, are also briefly discussed. Taken together, these methods provide rich insights into the interactions between microbes and how they influence microbial community productivity. We finally overview approaches that allow us to generate and test numerous synthetic community compositions, followed by tools and methodologies that can predict effective genetic interventions to further improve the productivity of communities. With impending advancements in high-throughput omics of microbial communities, the stage is set for the rapid expansion of microbial community engineering, with a significant impact on biotechnological processes.

Single-cell Digital Twins for Cancer Preclinical Investigation

Laboratory models derived from clinical samples represent a solid platform in preclinical research for drug testing and investigation of disease mechanisms. The integration of these laboratory models with their digital counterparts (i.e., predictive mathematical models) allows to set up digital twins essential to fully exploit their potential to face the enormous molecular complexity of human organisms. In particular, due to the close integration of cell metabolism with all other cellular processes, any perturbation in cellular physiology typically reflect on altered cells metabolic profiling. In this regard, changes in metabolism have been shown, also in our laboratory, to drive a causal role in the emergence of cancer disease. Nevertheless, a unique metabolic program does not describe the altered metabolic profile of all tumour cells due to many causes from genetic variability to intratumour heterogeneous dependency on nutrients consumption and metabolism by multiple co-existing subclones. Currently, fluxomics approaches just match with the necessity of characterizing the overall flux distribution of cells within given samples, by disregarding possible heterogeneous behaviors. For the purpose of stratifying cancer heterogeneous subpopulations, quantification of fluxes at the single-cell level is needed. To this aim, we here present a new computational framework called single-cell Flux Balance Analysis (scFBA) that aims to set up digital metabolic twins in the perspective of being better exploited within a framework that makes also use of laboratory patient cell models. In particular, scFBA aims at integrating single-cell RNA-seq data within computational population models in order to depict a snapshot of the corresponding single-cell metabolic phenotypes at a given moment, together with an unsupervised identification of metabolic subpopulations.

Cofactors revisited - Predicting the impact of flavoprotein-related diseases on a genome scale

Flavin adenine dinucleotide (FAD) and its precursor flavin mononucleotide (FMN) are redox cofactors that are required for the activity of more than hundred human enzymes. Mutations in the genes encoding these proteins cause severe phenotypes, including a lack of energy supply and accumulation of toxic intermediates. Ideally, patients should be diagnosed before they show symptoms so that treatment and/or preventive care can start immediately. This can be achieved by standardized newborn screening tests. However, many of the flavin-related diseases lack appropriate biomarker profiles. Genome-scale metabolic models can aid in biomarker research by predicting altered profiles of potential biomarkers. Unfortunately, current models, including the most recent human metabolic reconstructions Recon and HMR, typically treat enzyme-bound flavins incorrectly as free metabolites. This in turn leads to artificial degrees of freedom in pathways that are strictly coupled. Here, we present a reconstruction of human metabolism with a curated and extended flavoproteome. To illustrate the functional consequences, we show that simulations with the curated model - unlike simulations with earlier Recon versions - correctly predict the metabolic impact of multiple-acyl-CoA-dehydrogenase deficiency as well as of systemic flavin-depletion. Moreover, simulations with the new model allowed us to identify a larger number of biomarkers in flavoproteome-related diseases, without loss of accuracy. We conclude that adequate inclusion of cofactors in constraint-based modelling contributes to higher precision in computational predictions.

PFA toolbox: a MATLAB tool for Metabolic Flux Analysis

BACKGROUND

Metabolic Flux Analysis (MFA) is a methodology that has been successfully applied to estimate metabolic fluxes in living cells. However, traditional frameworks based on this approach have some limitations, particularly when measurements are scarce and imprecise. This is very common in industrial environments. The PFA Toolbox can be used to face those scenarios.

RESULTS

Here we present the PFA (Possibilistic Flux Analysis) Toolbox for MATLAB, which simplifies the use of Interval and Possibilistic Metabolic Flux Analysis. The main features of the PFA Toolbox are the following: (a) It provides reliable MFA estimations in scenarios where only a few fluxes can be measured or those available are imprecise. (b) It provides tools to easily plot the results as interval estimates or flux distributions. (c) It is composed of simple functions that MATLAB users can apply in flexible ways. (d) It includes a Graphical User Interface (GUI), which provides a visual representation of the measurements and their uncertainty. (e) It can use stoichiometric models in COBRA format. In addition, the PFA Toolbox includes a User's Guide with a thorough description of its functions and several examples.

CONCLUSIONS

The PFA Toolbox for MATLAB is a freely available Toolbox that is able to perform Interval and Possibilistic MFA estimations.

The steady-state assumption in oscillating and growing systems

The steady-state assumption, which states that the production and consumption of metabolites inside the cell are balanced, is one of the key aspects that makes an efficient analysis of genome-scale metabolic networks possible. It can be motivated from two different perspectives. In the time-scales perspective, we use the fact that metabolism is much faster than other cellular processes such as gene expression. Hence, the steady-state assumption is derived as a quasi-steady-state approximation of the metabolism that adapts to the changing cellular conditions. In this article we focus on the second perspective, stating that on the long run no metabolite can accumulate or deplete. In contrast to the first perspective it is not immediately clear how this perspective can be captured mathematically and what assumptions are required to obtain the steady-state condition. By presenting a mathematical framework based on the second perspective we demonstrate that the assumption of steady-state also applies to oscillating and growing systems without requiring quasi-steady-state at any time point. However, we also show that the average concentrations may not be compatible with the average fluxes. In summary, we establish a mathematical foundation for the steady-state assumption for long time periods that justifies its successful use in many applications. Furthermore, this mathematical foundation also pinpoints unintuitive effects in the integration of metabolite concentrations using nonlinear constraints into steady-state models for long time periods.

A constraint-based modelling approach to metabolic dysfunction in Parkinson's disease

One of the hallmarks of sporadic Parkinson's disease is degeneration of dopaminergic neurons in the pars compacta of the substantia nigra. The aetiopathogenesis of this degeneration is still not fully understood, with dysfunction of many biochemical pathways in different subsystems suggested to be involved. Recent advances in constraint-based modelling approaches hold great potential to systematically examine the relative contribution of dysfunction in disparate pathways to dopaminergic neuronal degeneration, but few studies have employed these methods in Parkinson's disease research. Therefore, this review outlines a framework for future constraint-based modelling of dopaminergic neuronal metabolism to decipher the multi-factorial mechanisms underlying the neuronal pathology of Parkinson's disease.