Virtual metabolic human dynamic model for pathological analysis and therapy design for diabetes

Longitudinal Metabolomics Data Analysis Informed by Mechanistic Models

Background: Metabolomics measurements are noisy, often characterized by a small sample size and missing entries. While data-driven methods have shown promise in terms of analyzing metabolomics data, e.g., revealing biomarkers of various phenotypes, metabolomics data analysis can significantly benefit from incorporating prior information about metabolic mechanisms. This paper introduces a novel data analysis approach to incorporate mechanistic models in metabolomics data analysis. Methods: We arranged time-resolved metabolomics measurements of plasma samples collected during a meal challenge test from the COPSAC2000 cohort as a third-order tensor: subjects by metabolites by time samples. Simulated challenge test data generated using a human whole-body metabolic model were also arranged as a third-order tensor: virtual subjects by metabolites by time samples. Real and simulated data sets were coupled in the metabolites mode and jointly analyzed using coupled tensor factorizations to reveal the underlying patterns. Results: Our experiments demonstrated that the joint analysis of simulated and real data had better performance in terms of pattern discovery, achieving higher correlations with a BMI (body mass index)-related phenotype compared to the analysis of only real data in males, while in females, the performance was comparable. We also demonstrated the advantages of such a joint analysis approach in the presence of incomplete measurements and its limitations in the presence of wrong prior information. Conclusions: The joint analysis of real measurements and simulated data (generated using a mechanistic model) through coupled tensor factorizations guides real data analysis with prior information encapsulated in mechanistic models and reveals interpretable patterns.

Modeling sex-specific whole-body metabolic responses to feeding and fasting

Men generally favor carbohydrate metabolism, while women lean towards lipid metabolism, resulting in significant sex-based differences in energy oxidation across various metabolic states such as fasting and feeding. These differences are influenced by body composition and inherent metabolic fluxes, including increased lipolysis rates in women. However, understanding how sex influences organ-specific metabolism and systemic manifestations remains incomplete. To address these gaps, we developed a sex-specific, whole-body metabolic model for feeding and fasting scenarios in healthy young adults. Our model integrates organ metabolism with whole-body responses to mixed meals, particularly high-carbohydrate and high-fat meals. Our predictions suggest that differences in liver and adipose tissue nutrient storage and oxidation patterns drive systemic metabolic disparities. We propose that sex differences in fasting hepatic glucose output may result from the different handling of free fatty acids, glycerol, and glycogen. We identified a metabolic pathway, possibly more prevalent in female livers, redirecting lipids towards carbohydrate metabolism to support hepatic glucose production. This mechanism is facilitated by the TG-FFA cycle between adipose tissue and the liver. Incorporating sex-specific data into multi-scale frameworks offers insights into how sex modulates human metabolism.

Integrated modeling of labile and glycated hemoglobin with glucose for enhanced diabetes detection and short-term monitoring

Metabolic biomarkers, particularly glycated hemoglobin and fasting plasma glucose, are pivotal in the diagnosis and control of diabetes mellitus. Despite their importance, they exhibit limitations in assessing short-term glucose variations. In this study, we propose labile hemoglobin as an additional biomarker, providing insightful perspectives into these fluctuations. By utilizing datasets from 40,652 retrospective general participants and conducting glucose tolerance tests on 60 prospective pediatric subjects, we explored the relationship between plasma glucose and labile hemoglobin. A mathematical model was developed to encapsulate short-term glucose kinetics in the pediatric group. Applying dimensionality reduction techniques, we successfully identified participant subclusters, facilitating the differentiation between diabetic and non-diabetic individuals. Intriguingly, by integrating labile hemoglobin measurements with plasma glucose values, we were able to predict the likelihood of diabetes in pediatric subjects, underscoring the potential of labile hemoglobin as a significant glycemic biomarker for diabetes research.

Analyzing postprandial metabolomics data using multiway models: a simulation study

BACKGROUND

Analysis of time-resolved postprandial metabolomics data can improve the understanding of metabolic mechanisms, potentially revealing biomarkers for early diagnosis of metabolic diseases and advancing precision nutrition and medicine. Postprandial metabolomics measurements at several time points from multiple subjects can be arranged as a subjects by metabolites by time points array. Traditional analysis methods are limited in terms of revealing subject groups, related metabolites, and temporal patterns simultaneously from such three-way data.

RESULTS

We introduce an unsupervised multiway analysis approach based on the CANDECOMP/PARAFAC (CP) model for improved analysis of postprandial metabolomics data guided by a simulation study. Because of the lack of ground truth in real data, we generate simulated data using a comprehensive human metabolic model. This allows us to assess the performance of CP models in terms of revealing subject groups and underlying metabolic processes. We study three analysis approaches: analysis of fasting-state data using principal component analysis, T0-corrected data (i.e., data corrected by subtracting fasting-state data) using a CP model and full-dynamic (i.e., full postprandial) data using CP. Through extensive simulations, we demonstrate that CP models capture meaningful and stable patterns from simulated meal challenge data, revealing underlying mechanisms and differences between diseased versus healthy groups.

CONCLUSIONS

Our experiments show that it is crucial to analyze both fasting-state and T0-corrected data for understanding metabolic differences among subject groups. Depending on the nature of the subject group structure, the best group separation may be achieved by CP models of T0-corrected or full-dynamic data. This study introduces an improved analysis approach for postprandial metabolomics data while also shedding light on the debate about correcting baseline values in longitudinal data analysis.

Data-driven analysis and prediction of dynamic postprandial metabolic response to multiple dietary challenges using dynamic mode decomposition

Motivation

In the field of precision nutrition, predicting metabolic response to diet and identifying groups of differential responders are two highly desirable steps toward developing tailored dietary strategies. However, data analysis tools are currently lacking, especially for complex settings such as crossover studies with repeated measures.Current methods of analysis often rely on matrix or tensor decompositions, which are well suited for identifying differential responders but lacking in predictive power, or on dynamical systems modeling, which may be used for prediction but typically requires detailed mechanistic knowledge of the system under study. To remedy these shortcomings, we explored dynamic mode decomposition (DMD), which is a recent, data-driven method for deriving low-rank linear dynamical systems from high dimensional data.Combining the two recent developments "parametric DMD" (pDMD) and "DMD with control" (DMDc) enabled us to (i) integrate multiple dietary challenges, (ii) predict the dynamic response in all measured metabolites to new diets from only the metabolite baseline and dietary input, and (iii) identify inter-individual metabolic differences, i.e., metabotypes. To our knowledge, this is the first time DMD has been applied to analyze time-resolved metabolomics data.

Results

We demonstrate the potential of pDMDc in a crossover study setting. We could predict the metabolite response to unseen dietary exposures on both measured (R2 = 0.40) and simulated data of increasing size (R max 2 = 0.65), as well as recover clusters of dynamic metabolite responses. We conclude that this method has potential for applications in personalized nutrition and could be useful in guiding metabolite response to target levels.

Availability and implementation

The measured data analyzed in this study can be provided upon reasonable request. The simulated data along with a MATLAB implementation of pDMDc is available at https://github.com/FraunhoferChalmersCentre/pDMDc.

Technologies for whole-cell modeling: Genome-wide reconstruction of a cell in silico

With advances in high-throughput, large-scale in vivo measurement and genome modification techniques at the single-nucleotide level, there is an increasing demand for the development of new technologies for the flexible design and control of cellular systems. Computer-aided design is a powerful tool to design new cells. Whole-cell modeling aims to integrate various cellular subsystems, determine their interactions and cooperative mechanisms, and predict comprehensive cellular behaviors by computational simulations on a genome-wide scale. It has been applied to prokaryotes, yeasts, and higher eukaryotic cells, and utilized in a wide range of applications, including production of valuable substances, drug discovery, and controlled differentiation. Whole-cell modeling, consisting of several thousand elements with diverse scales and properties, requires innovative model construction, simulation, and analysis techniques. Furthermore, whole-cell modeling has been extended to multiple scales, including high-resolution modeling at the single-nucleotide and single-amino acid levels and multicellular modeling of tissues and organs. This review presents an overview of the current state of whole-cell modeling, discusses the novel computational and experimental technologies driving it, and introduces further developments toward multihierarchical modeling on a whole-genome scale.

A computational model of the effects of macronutrients absorption and physical exercise on hormonal regulation and metabolic homeostasis

Regular physical exercise and appropriate nutrition affect metabolic and hormonal responses and may reduce the risk of developing chronic non-communicable diseases such as high blood pressure, ischemic stroke, coronary heart disease, some types of cancer, and type 2 diabetes mellitus. Computational models describing the metabolic and hormonal changes due to the synergistic action of exercise and meal intake are, to date, scarce and mostly focussed on glucose absorption, ignoring the contribution of the other macronutrients. We here describe a model of nutrient intake, stomach emptying, and absorption of macronutrients in the gastrointestinal tract during and after the ingestion of a mixed meal, including the contribution of proteins and fats. We integrated this effort to our previous work in which we modeled the effects of a bout of physical exercise on metabolic homeostasis. We validated the computational model with reliable data from the literature. The simulations are overall physiologically consistent and helpful in describing the metabolic changes due to everyday life stimuli such as multiple mixed meals and variable periods of physical exercise over prolonged periods of time. This computational model may be used to design virtual cohorts of subjects differing in sex, age, height, weight, and fitness status, for specialized in silico challenge studies aimed at designing exercise and nutrition schemes to support health.

Digital twin predicting diet response before and after long-term fasting

Today, there is great interest in diets proposing new combinations of macronutrient compositions and fasting schedules. Unfortunately, there is little consensus regarding the impact of these different diets, since available studies measure different sets of variables in different populations, thus only providing partial, non-connected insights. We lack an approach for integrating all such partial insights into a useful and interconnected big picture. Herein, we present such an integrating tool. The tool uses a novel mathematical model that describes mechanisms regulating diet response and fasting metabolic fluxes, both for organ-organ crosstalk, and inside the liver. The tool can mechanistically explain and integrate data from several clinical studies, and correctly predict new independent data, including data from a new study. Using this model, we can predict non-measured variables, e.g. hepatic glycogen and gluconeogenesis, in response to fasting and different diets. Furthermore, we exemplify how such metabolic responses can be successfully adapted to a specific individual’s sex, weight, height, as well as to the individual’s historical data on metabolite dynamics. This tool enables an offline digital twin technology.

Fasting and diet are central components of prevention against cardiovascular disease. Unfortunately, there is little consensus regarding which diet schemes are optimal. This is partially because different clinical studies contribute with different non-connected pieces of knowledge, which have not been fully integrated into a useful and interconnected big picture. In principle, mathematical models describing meal responses could be used for such an integration. However, today’s models still lack critical mechanisms, such as protein metabolism and a dynamic glycogen regulation. Herein, we present a) a new expanded model structure including these mechanisms; b) a set of parameters which can simultaneously describe a wide array of complementary estimation data, in both healthy and diabetic populations; c) a personalisation-script, which allows these generic parameters to be tuned to an individual/sub-population, using demographics (age, weight, height, diabetes status) and historic metabolic data. We exemplify how this personalisation can be used to predict new independent data, including a new clinical study, where a qualitatively new prediction is validated: that an oral protein tolerance test gives a clear response in plasma glucose, after, but not before, a 48h fasting period. Our combined model, parameters, and fitting script lay the foundation for an offline digital twin.

The transition between acute and chronic infections in light of energy control: a mathematical model of energy flow in response to infection

BACKGROUND

Different parts of an organism like the gut, endocrine, nervous and immune systems constantly exchange information. Understanding the pathogenesis of various systemic chronic diseases increasingly relies on understanding how these subsystems orchestrate their activities.

METHODS

We started from the working hypothesis that energy is a fundamental quantity that governs activity levels of all subsystems and that interactions between subsystems control the distribution of energy according to acute needs. Based on physiological knowledge, we constructed a mathematical model for the energy flow between subsystems and analysed the resulting organismal responses to in silico infections.

RESULTS

The model reproduces common behaviour in acute infections and suggests several host parameters that modulate infection duration and therapeutic responsiveness. Moreover, the model allows the formulation of conditions for the induction of chronic infections and predicts that alterations in energy released from fat can lead to the transition from clearance of acute infections to a chronic inflammatory state.

IMPACT

These results suggest a fundamental role for brain and fat in controlling immune response through systemic energy control. In particular, it suggests that lipolysis resistance, which is known to be involved in obesity and ageing, might be a survival programme for coping with chronic infections.

New Approaches to Profile the Microbiome for Treatment of Neurodegenerative Disease

Progressive neurodegenerative diseases represent some of the largest growing treatment challenges for public health in modern society. These diseases mainly progress due to aging and are driven by microglial surveillance and activation in response to changes occurring in the aging brain. The lack of efficacious treatment options for Alzheimer's disease (AD), as the focus of this review, and other neurodegenerative disorders has encouraged new approaches to address neuroinflammation for potential treatments. Here we will focus on the increasing evidence that dysbiosis of the gut microbiome is characterized by inflammation that may carry over to the central nervous system and into the brain. Neuroinflammation is the common thread associated with neurodegenerative diseases, but it is yet unknown at what point and how innate immune function turns pathogenic for an individual. This review will address extensive efforts to identify constituents of the gut microbiome and their neuroactive metabolites as a peripheral path to treatment. This approach is still in its infancy in substantive clinical trials and requires thorough human studies to elucidate the metabolic microbiome profile to design appropriate treatment strategies for early stages of neurodegenerative disease. We view that in order to address neurodegenerative mechanisms of the gut, microbiome and metabolite profiles must be determined to pre-screen AD subjects prior to the design of specific, chronic titrations of gut microbiota with low-dose antibiotics. This represents an exciting treatment strategy designed to balance inflammatory microglial involvement in disease progression with an individual's manifestation of AD as influenced by a coercive inflammatory gut.