Putting the pieces together in diabetes research: towards a hierarchical model of whole-body glucose homeostasis

A multi-scale digital twin for adiposity-driven insulin resistance in humans: diet and drug effects

BACKGROUND

The increased prevalence of insulin resistance is one of the major health risks in society today. Insulin resistance involves both short-term dynamics, such as altered meal responses, and long-term dynamics, such as the development of type 2 diabetes. Insulin resistance also occurs on different physiological levels, ranging from disease phenotypes to organ-organ communication and intracellular signaling. To better understand the progression of insulin resistance, an analysis method is needed that can combine different timescales and physiological levels. One such method is digital twins, consisting of combined mechanistic mathematical models. We have previously developed a model for short-term glucose homeostasis and intracellular insulin signaling, and there exist long-term weight regulation models. Herein, we combine these models into a first interconnected digital twin for the progression of insulin resistance in humans.

METHODS

The model is based on ordinary differential equations representing biochemical and physiological processes, in which unknown parameters were fitted to data using a MATLAB toolbox.

RESULTS

The interconnected twin correctly predicts independent data from a weight increase study, both for weight-changes, fasting plasma insulin and glucose levels, and intracellular insulin signaling. Similarly, the model can predict independent weight-change data in a weight loss study with the weight loss drug topiramate. The model can also predict non-measured variables.

CONCLUSIONS

The model presented herein constitutes the basis for a new digital twin technology, which in the future could be used to aid medical pedagogy and increase motivation and compliance and thus aid in the prevention and treatment of insulin resistance.

Deciphering signal transduction networks in the liver by mechanistic mathematical modelling

In health and disease, liver cells are continuously exposed to cytokines and growth factors. While individual signal transduction pathways induced by these factors were studied in great detail, the cellular responses induced by repeated or combined stimulations are complex and less understood. Growth factor receptors on the cell surface of hepatocytes were shown to be regulated by receptor interactions, receptor trafficking and feedback regulation. Here, we exemplify how mechanistic mathematical modelling based on quantitative data can be employed to disentangle these interactions at the molecular level. Crucial is the analysis at a mechanistic level based on quantitative longitudinal data within a mathematical framework. In such multi-layered information, step-wise mathematical modelling using submodules is of advantage, which is fostered by sharing of standardized experimental data and mathematical models. Integration of signal transduction with metabolic regulation in the liver and mechanistic links to translational approaches promise to provide predictive tools for biology and personalized medicine.

A systems biology analysis of lipolysis and fatty acid release from adipocytes in vitro and from adipose tissue in vivo

Lipolysis and the release of fatty acids to supply energy fuel to other organs, such as between meals, during exercise, and starvation, are fundamental functions of the adipose tissue. The intracellular lipolytic pathway in adipocytes is activated by adrenaline and noradrenaline, and inhibited by insulin. Circulating fatty acids are elevated in type 2 diabetic individuals. The mechanisms behind this elevation are not fully known, and to increase the knowledge a link between the systemic circulation and intracellular lipolysis is key. However, data on lipolysis and knowledge from in vitro systems have not been linked to corresponding in vivo data and knowledge in vivo. Here, we use mathematical modelling to provide such a link. We examine mechanisms of insulin action by combining in vivo and in vitro data into an integrated mathematical model that can explain all data. Furthermore, the model can describe independent data not used for training the model. We show the usefulness of the model by simulating new and more challenging experimental setups in silico, e.g. the extracellular concentration of fatty acids during an insulin clamp, and the difference in such simulations between individuals with and without type 2 diabetes. Our work provides a new platform for model-based analysis of adipose tissue lipolysis, under both non-diabetic and type 2 diabetic conditions.

CaliPro: A Calibration Protocol That Utilizes Parameter Density Estimation to Explore Parameter Space and Calibrate Complex Biological Models

Introduction

Mathematical and computational modeling have a long history of uncovering mechanisms and making predictions for biological systems. However, to create a model that can provide relevant quantitative predictions, models must first be calibrated by recapitulating existing biological datasets from that system. Current calibration approaches may not be appropriate for complex biological models because: 1) many attempt to recapitulate only a single aspect of the experimental data (such as a median trend) or 2) Bayesian techniques require specification of parameter priors and likelihoods to experimental data that cannot always be confidently assigned. A new calibration protocol is needed to calibrate complex models when current approaches fall short.

Methods

Herein, we develop CaliPro, an iterative, model-agnostic calibration protocol that utilizes parameter density estimation to refine parameter space and calibrate to temporal biological datasets. An important aspect of CaliPro is the user-defined pass set definition, which specifies how the model might successfully recapitulate experimental data. We define the appropriate settings to use CaliPro.

Results

We illustrate the usefulness of CaliPro through four examples including predator-prey, infectious disease transmission, and immune response models. We show that CaliPro works well for both deterministic, continuous model structures as well as stochastic, discrete models and illustrate that CaliPro can work across diverse calibration goals.

Conclusions

We present CaliPro, a new method for calibrating complex biological models to a range of experimental outcomes. In addition to expediting calibration, CaliPro may be useful in already calibrated parameter spaces to target and isolate specific model behavior for further analysis.

Type 2 diabetes progression: A regulatory network approach

In order to elucidate central elements underlying type 2 diabetes, we constructed a regulatory network model involving 37 components (molecules, receptors, processes, etc.) associated to signaling pathways of pancreatic beta-cells. In a first approximation, the network topology was described by Boolean rules whose interacting dynamics predicted stationary patterns broadly classified as health, metabolic syndrome, and diabetes stages. A subsequent approximation based on a continuous logic analysis allowed us to characterize the progression of the disease as transitions between these states associated to alterations of cell homeostasis due to exhaustion or exacerbation of specific regulatory signals. The method allowed the identification of key transcription factors involved in metabolic stress as essential for the progression of the disease. Integration of the present analysis with existent mathematical models designed to yield accurate account of experimental data in human or animal essays leads to reliable predictions for beta-cell mass, insulinemia, glycemia, and glycosylated hemoglobin in diabetic fatty rats.

Mathematical model of blood glucose dynamics by emulating the pathophysiology of glucose metabolism in type 2 diabetes mellitus

Mathematical modelling has established itself as a theoretical tool to understand fundamental aspects of a variety of medical-biological phenomena. The predictive power of mathematical models on some chronic conditions has been helpful in its proper prevention, diagnosis, and treatment. Such is the case of the modelling of glycaemic dynamics in type 2 diabetes mellitus (T2DM), whose physiology-based mathematical models have captured the metabolic abnormalities of this disease. Through a physiology-based pharmacokinetic-pharmacodynamic approach, this work addresses a mathematical model whose structure starts from a model of blood glucose dynamics in healthy humans. This proposal is capable of emulating the pathophysiology of T2DM metabolism, including the effect of gastric emptying and insulin enhancing effect due to incretin hormones. The incorporation of these effects lies in the implemented methodology since the mathematical functions that represent metabolic rates, with a relevant contribution to hyperglycaemia, are adjusting individually to the clinical data of patients with T2DM. Numerically, the resulting model successfully simulates a scheduled graded intravenous glucose test and oral glucose tolerance tests at different doses. The comparison between simulations and clinical data shows an acceptable description of the blood glucose dynamics in T2DM. It opens the possibility of using this model to develop model-based controllers for the regulation of blood glucose in T2DM.

Computational Modeling in Liver Surgery

The need for extended liver resection is increasing due to the growing incidence of liver tumors in aging societies. Individualized surgical planning is the key for identifying the optimal resection strategy and to minimize the risk of postoperative liver failure and tumor recurrence. Current computational tools provide virtual planning of liver resection by taking into account the spatial relationship between the tumor and the hepatic vascular trees, as well as the size of the future liver remnant. However, size and function of the liver are not necessarily equivalent. Hence, determining the future liver volume might misestimate the future liver function, especially in cases of hepatic comorbidities such as hepatic steatosis. A systems medicine approach could be applied, including biological, medical, and surgical aspects, by integrating all available anatomical and functional information of the individual patient. Such an approach holds promise for better prediction of postoperative liver function and hence improved risk assessment. This review provides an overview of mathematical models related to the liver and its function and explores their potential relevance for computational liver surgery. We first summarize key facts of hepatic anatomy, physiology, and pathology relevant for hepatic surgery, followed by a description of the computational tools currently used in liver surgical planning. Then we present selected state-of-the-art computational liver models potentially useful to support liver surgery. Finally, we discuss the main challenges that will need to be addressed when developing advanced computational planning tools in the context of liver surgery.

Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions

Recent and rapidly evolving progress on high-throughput measurement techniques and computational performance has led to the emergence of new disciplines, such as systems medicine and translational systems biology. At the core of these disciplines lies the desire to produce multiscale models: mathematical models that integrate multiple scales of biological organization, ranging from molecular, cellular and tissue models to organ, whole-organism and population scale models. Using such models, hypotheses can systematically be tested. In this review, we present state-of-the-art multiscale modelling of bacterial and fungal infections, considering both the pathogen and host as well as their interaction. Multiscale modelling of the interactions of bacteria, especially Mycobacterium tuberculosis, with the human host is quite advanced. In contrast, models for fungal infections are still in their infancy, in particular regarding infections with the most important human pathogenic fungi, Candida albicans and Aspergillus fumigatus. We reflect on the current availability of computational approaches for multiscale modelling of host-pathogen interactions and point out current challenges. Finally, we provide an outlook for future requirements of multiscale modelling.

Requirements for multi-level systems pharmacology models to reach end-usage: the case of type 2 diabetes

We are currently in the middle of a major shift in biomedical research: unprecedented and rapidly growing amounts of data may be obtained today, from in vitro, in vivo and clinical studies, at molecular, physiological and clinical levels. To make use of these large-scale, multi-level datasets, corresponding multi-level mathematical models are needed, i.e. models that simultaneously capture multiple layers of the biological, physiological and disease-level organization (also referred to as quantitative systems pharmacology-QSP-models). However, today's multi-level models are not yet embedded in end-usage applications, neither in drug research and development nor in the clinic. Given the expectations and claims made historically, this seemingly slow adoption may seem surprising. Therefore, we herein consider a specific example-type 2 diabetes-and critically review the current status and identify key remaining steps for these models to become mainstream in the future. This overview reveals how, today, we may use models to ask scientific questions concerning, e.g., the cellular origin of insulin resistance, and how this translates to the whole-body level and short-term meal responses. However, before these multi-level models can become truly useful, they need to be linked with the capabilities of other important existing models, in order to make them 'personalized' (e.g. specific to certain patient phenotypes) and capable of describing long-term disease progression. To be useful in drug development, it is also critical that the developed models and their underlying data and assumptions are easily accessible. For clinical end-usage, in addition, model links to decision-support systems combined with the engagement of other disciplines are needed to create user-friendly and cost-efficient software packages.

A comprehensive compartmental model of blood glucose regulation for healthy and type 2 diabetic subjects

We have expanded a former compartmental model of blood glucose regulation for healthy and type 2 diabetic subjects. The former model was a detailed physiological model which considered the interactions of three substances, glucose, insulin and glucagon on regulating the blood sugar. The main drawback of the former model was its restriction on the route of glucose entrance to the body which was limited to the intravenous glucose injection. To handle the oral glucose intake, we have added a model of glucose absorption in the gastrointestinal tract to the former model to address the resultant variations of blood glucose concentrations following an oral glucose intake. Another model representing the incretins production in the gastrointestinal tract along with their hormonal effects on boosting pancreatic insulin production is also added to the former model. We have used two sets of clinical data obtained during oral glucose tolerance test and isoglycemic intravenous glucose infusion test from both type 2 diabetic and healthy subjects to estimate the model parameters and to validate the model results. The estimation of model parameters is accomplished through solving a nonlinear optimization problem. The results show acceptable precision of the estimated model parameters and demonstrate the capability of the model in accurate prediction of the body response during the clinical studies.

Evaluating systems pharmacology models is different from evaluating standard pharmacokinetic-pharmacodynamic models

Based on the author's recent experience, there appears to be some confusion regarding the steps required to qualify a systems pharmacology model as adequate for the intended purpose. This manuscript outlines the model evaluation approach used in the author's recent publication¹on the systems pharmacology of a 5-lipoxygenase inhibitor and is an attempt to generate discussion on this topic within the pharmacometrics and systems pharmacology community.

The impact of mathematical modeling on the understanding of diabetes and related complications

Diabetes is a chronic and complex multifactorial disease caused by persistent hyperglycemia and for which underlying pathogenesis is still not completely understood. The mathematical modeling of glucose homeostasis, diabetic condition, and its associated complications is rapidly growing and provides new insights into the underlying mechanisms involved. Here, we discuss contributions to the diabetes modeling field over the past five decades, highlighting the areas where more focused research is required.

Multilevel model based glucose control for type-1 diabetes patients

Diabetes is a disease that involves alterations at multiple biological levels, ranging from intracellular signalling to organ processes. Since glucose homeostasis is the consequence of complex interactions that involve a number of factors, the control of diabetes should be based on a multilevel analysis. In this paper, a novel approach to design of closed-loop glucose controllers based on multilevel models is presented. A control scheme is proposed based on combining a pharmacokinetic/pharmacodynamic model with an insulin signal transduction model for type 1 diabetes mellitus patients. Based on this, an insulin feedback control schemes is designed. Two main advantages of explicitly utilizing information at the intracellular level were obtained. First, significant reduction of hypoglycaemic risk by reducing the undershoot in glucose levels in response to added insulin. Second, robust performance for inter-patient changes, demonstrated through application of the multilevel control strategy to a well established in silico population of diabetic patients.

Insulin signaling - mathematical modeling comes of age

Signaling pathways that only a few years ago appeared simple and understandable, albeit far from complete, have evolved into very complex multi-layered networks of cellular control mechanisms, which in turn are integrated in a similarly complex whole-body level of control mechanisms. This complexity sets limits for classical biochemical reasoning, such that a correct and complete analysis of experimental data while taking the full complexity into account is not possible. In this Opinion we propose that mathematical modeling can be used as a tool in insulin signaling research, and we demonstrate how recent developments in modeling - and the integration of modeling in the experimental process - provide new possibilities to approach and decipher complex biological systems more efficiently.

Zooming of states and parameters using a lumping approach including back-translation

Background

Systems biology models tend to become large since biological systems often consist of complex networks of interacting components, and since the models usually are developed to reflect various mechanistic assumptions of those networks. Nevertheless, not all aspects of the model are equally interesting in a given setting, and normally there are parts that can be reduced without affecting the relevant model performance. There are many methods for model reduction, but few or none of them allow for a restoration of the details of the original model after the simplified model has been simulated.

Results

We present a reduction method that allows for such a back-translation from the reduced to the original model. The method is based on lumping of states, and includes a general and formal algorithm for both determining appropriate lumps, and for calculating the analytical back-translation formulas. The lumping makes use of efficient methods from graph-theory and ϵ-decomposition and is derived and exemplified on two published models for fluorescence emission in photosynthesis. The bigger of these models is reduced from 26 to 6 states, with a negligible deviation from the reduced model simulations, both when comparing simulations in the states of the reduced model and when comparing back-translated simulations in the states of the original model. The method is developed in a linear setting, but we exemplify how the same concepts and approaches can be applied to non-linear problems. Importantly, the method automatically provides a reduced model with back-translations. Also, the method is implemented as a part of the systems biology toolbox for matlab, and the matlab scripts for the examples in this paper are available in the supplementary material.

Conclusions

Our novel lumping methodology allows for both automatic reduction of states using lumping, and for analytical retrieval of the original states and parameters without performing a new simulation. The two models can thus be considered as two degrees of zooming of the same model. This is a conceptually new development of model reduction approaches, which we think will stimulate much further research and will prove to be very useful in future modelling projects.