Search

Abstract Current research on metabolic disorders and diabetes relies on animal models because multi-organ diseases cannot be well studied with standard in vitro assays. Here, we have connected cell models of key metabolic organs, the pancreas and liver, on a microfluidic chip to enable diabetes research in a human-based in vitro system. Aided by mechanistic mathematical modeling, we demonstrate that hyperglycemia and high cortisone concentration induce glucose dysregulation in the pancreas-liver microphysiological system (MPS), mimicking a diabetic phenotype seen in patients with glucocorticoid-induced diabetes. In this diseased condition, the pancreas-liver MPS displays beta-cell dysfunction, steatosis, elevated ketone-body secretion, increased glycogen storage, and upregulated gluconeogenic gene expression. Conversely, a physiological culture condition maintains glucose tolerance and beta-cell function. This method was reproducible in two laboratories and was effective in multiple pancreatic islet donors. The model also provides a platform to identify new therapeutic proteins, as demonstrated with a combined transcriptome and proteome analysis.

Abstract Alcohol consumption is associated with a wide variety of preventable health complications and is a major risk factor for all-cause mortality in the age group 15-47 years. To reduce dangerous drinking behavior, eHealth applications have shown promise. A particularly interesting potential lies in the combination of eHealth apps with mathematical models. However, existing mathematical models do not consider real-life situations, such as combined intake of meals and beverages, and do not connect drinking to clinical markers, such as phosphatidylethanol (PEth). Herein, we present such a model which can simulate real-life situations and connect drinking to long-term markers. The new model can accurately describe both estimation data according to a χ2 -test (187.0

Abstract Background Liver cirrhosis, the advanced stage of many chronic liver diseases, is associated with escalated risks of liver-related complications like decompensation and hepatocellular carcinoma (HCC). Morbidity and mortality in cirrhosis patients are linked to portal hypertension, sarcopenia, and hepatocellular carcinoma. Although conventional cirrhosis management centered on treating complications, contemporary approaches prioritize preemptive measures. This study aims to formulate novel blood- and imaging-centric methodologies for monitoring liver cirrhosis patients. Methods In this prospective study, 150 liver cirrhosis patients will be enrolled from three Swedish liver clinics. Their conditions will be assessed through extensive blood-based markers and magnetic resonance imaging (MRI). The MRI protocol encompasses body composition profile with Muscle Assement Score, portal flow assessment, magnet resonance elastography, and a abbreviated MRI for HCC screening. Evaluation of lifestyle, muscular strength, physical performance, body composition, and quality of life will be conducted. Additionally, DNA, serum, and plasma biobanking will facilitate future investigations. Discussion The anticipated outcomes involve the identification and validation of non-invasive blood- and imaging-oriented biomarkers, enhancing the care paradigm for liver cirrhosis patients. Notably, the temporal evolution of these biomarkers will be crucial for understanding dynamic changes. Trial registration Clinicaltrials.gov, registration identifier NCT05502198. Registered on 16 August 2022. Link: https://classic.clinicaltrials.gov/ct2/show/NCT05502198 .

Abstract 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.

IntroductionCerebrovascular reactivity imaging (CVR) is a diagnostic method for assessment of alterations in cerebral blood flow in response to a controlled vascular stimulus. The principal utility is the capacity to evaluate the cerebrovascular reserve, thereby elucidating autoregulatory functioning. In CVR, CO2 gas challenge is the most prevalent method, which elicits a vascular response by alterations in inspired CO2 concentrations. While several systems have been proposed in the literature, only a limited number have been devised to operate in tandem with mechanical ventilation, thus constraining the majority CVR investigations to spontaneously breathing individuals.MethodsWe have developed a new method, denoted Additional CO2, designed to enable CO2 challenge in ventilators. The central idea is the introduction of an additional flow of highly concentrated CO2 into the respiratory circuit, as opposed to administration of the entire gas mixture from a reservoir. By monitoring the main respiratory gas flow emanating from the ventilator, the CO2 concentration in the inspired gas can be manipulated by adjusting the proportion of additional CO2. We evaluated the efficacy of this approach in (1) a ventilator coupled with a test lung and (2) in spontaneously breathing healthy subjects. The method was evaluated by assessment of the precision in attaining target inspired CO2 levels and examination of its performance within a magnetic resonance imaging environment.Results and discussionOur investigations revealed that the Additional CO2 method consistently achieved a high degree of accuracy in reaching target inspired CO2 levels in both mechanical ventilation and spontaneous breathing. We anticipate that these findings will lay the groundwork for a broader implementation of CVR assessments in mechanically ventilated patients.

Abstract Adipocyte signaling, normally and in type 2 diabetes, is far from fully understood. We have earlier developed detailed dynamic mathematical models for several well-studied, partially overlapping, signaling pathways in adipocytes. Still, these models only cover a fraction of the total cellular response. For a broader coverage of the response, large-scale phosphoproteomic data and systems level knowledge on protein interactions are key. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. We have developed a method to first establish a core model by connecting existing models of adipocyte cellular signaling for: (1) lipolysis and fatty acid release, (2) glucose uptake, and (3) the release of adiponectin. Next, we use publicly available phosphoproteome data for the insulin response in adipocytes together with prior knowledge on protein interactions, to identify phosphosites downstream of the core model. In a parallel pairwise approach with low computation time, we test whether identified phosphosites can be added to the model. We iteratively collect accepted additions into layers and continue the search for phosphosites downstream of these added layers. For the first 30 layers with the highest confidence (311 added phosphosites), the model predicts independent data well (70–90% correct), and the predictive capability gradually decreases when we add layers of decreasing confidence. In total, 57 layers (3059 phosphosites) can be added to the model with predictive ability kept. Finally, our large-scale, layered model enables dynamic simulations of systems-wide alterations in adipocytes in type 2 diabetes.

Summary: Nuclear transport is an essential part of eukaryotic cell function. Here, we present scFRAP, a model-assisted fluorescent recovery after photobleaching (FRAP)- based method to determine nuclear import and export rates independently in individual live cells. To overcome the inherent noise of single-cell measurements, we performed sequential FRAPs on the same cell. We found large cell-to-cell variation in transport rates within isogenic yeast populations. For passive transport, the variability in NPC number might explain most of the variability. Using this approach, we studied mother-daughter cell asymmetry in the active nuclear shuttling of the transcription factor Ace2, which is specifically concentrated in daughter cell nuclei in early G1. Rather than reduced export in the daughter cell, as previously hypothesized, we found that this asymmetry is mainly due to an increased import in daughters. These results shed light on cell-to-cell variation in cellular dynamics and its sources.

Neurons regulate the activity of blood vessels through the neurovascular coupling (NVC). A detailed understanding of the NVC is critical for understanding data from functional imaging techniques of the brain. Many aspects of the NVC have been studied both experimentally and using mathematical models; various combinations of blood volume and flow, local field potential (LFP), hemoglobin level, blood oxygenation level-dependent response (BOLD), and optogenetics have been measured and modeled in rodents, primates, or humans. However, these data have not been brought together into a unified quantitative model. We now present a mathematical model that describes all such data types and that preserves mechanistic behaviors between experiments. For instance, from modeling of optogenetics and microscopy data in mice, we learn cell-specific contributions; the first rapid dilation in the vascular response is caused by NO-interneurons, the main part of the dilation during longer stimuli is caused by pyramidal neurons, and the post-peak undershoot is caused by NPY-interneurons. These insights are translated and preserved in all subsequent analyses, together with other insights regarding hemoglobin dynamics and the LFP/BOLD-interplay, obtained from other experiments on rodents and primates. The model can predict independent validation-data not used for training. By bringing together data with complementary information from different species, we both understand each dataset better, and have a basis for a new type of integrative analysis of human data.

BackgroundDrug induced liver injury (DILI) is a major concern when developing new drugs. A promising biomarker for DILI is the hepatic uptake rate of the contrast agent gadoxetate. This rate can be estimated using a novel approach combining magnetic resonance imaging and mathematical modeling. However, previous work has used different mathematical models to describe liver function in humans or rats, and no comparative study has assessed which model is most optimal to use, or focused on possible translatability between the two species.AimsOur aim was therefore to do a comparison and assessment of models for DILI biomarker assessment, and to develop a conceptual basis for a translational framework between the species.Methods and resultsWe first established which of the available pharmacokinetic models to use by identifying the most simple and identifiable model that can describe data from both human and rats. We then developed an extension of this model for how to estimate the effects of a hepatotoxic drug in rats. Finally, we illustrated how such a framework could be useful for drug dosage selection, and how it potentially can be applied in personalized treatments designed to avoid DILI.ConclusionOur analysis provides clear guidelines of which mathematical model to use for model-based assessment of biomarkers for liver function, and it also suggests a hypothetical path to a translational framework for DILI.

The neurovascular and neurometabolic couplings (NVC and NMC) connect cerebral activity, blood flow, and metabolism. This interconnection is used in for instance functional imaging, which analyses the blood-oxygen-dependent (BOLD) signal. The mechanisms underlying the NVC are complex, which warrants a model-based analysis of data. We have previously developed a mechanistically detailed model for the NVC, and others have proposed detailed models for cerebral metabolism. However, existing metabolic models are still not fully utilizing available magnetic resonance spectroscopy (MRS) data and are not connected to detailed models for NVC. Therefore, we herein present a new model that integrates mechanistic modelling of both MRS and BOLD data. The metabolic model covers central metabolism, using a minimal set of interactions, and can describe time-series data for glucose, lactate, aspartate, and glutamate, measured after visual stimuli. Statistical tests confirm that the model can describe both estimation data and predict independent validation data, not used for model training. The interconnected NVC model can simultaneously describe BOLD data and can be used to predict expected metabolic responses in experiments where metabolism has not been measured. This model is a step towards a useful and mechanistically detailed model for cerebral blood flow and metabolism, with potential applications in both basic research and clinical applications.

Microphysiological systems (MPS) are powerful tools for emulating human physiology and replicating disease progression in vitro. MPS could be better predictors of human outcome than current animal models, but mechanistic interpretation and in vivo extrapolation of the experimental results remain significant challenges. Here, we address these challenges using an integrated experimental-computational approach. This approach allows for in silico representation and predictions of glucose metabolism in a previously reported MPS with two organ compartments (liver and pancreas) connected in a closed loop with circulating medium. We developed a computational model describing glucose metabolism over 15 days of culture in the MPS. The model was calibrated on an experiment-specific basis using data from seven experiments, where HepaRG single-liver or liver-islet cultures were exposed to both normal and hyperglycemic conditions resembling high blood glucose levels in diabetes. The calibrated models reproduced the fast (i.e. hourly) variations in glucose and insulin observed in the MPS experiments, as well as the long-term (i.e. over weeks) decline in both glucose tolerance and insulin secretion. We also investigated the behaviour of the system under hypoglycemia by simulating this condition in silico, and the model could correctly predict the glucose and insulin responses measured in new MPS experiments. Last, we used the computational model to translate the experimental results to humans, showing good agreement with published data of the glucose response to a meal in healthy subjects. The integrated experimental-computational framework opens new avenues for future investigations toward disease mechanisms and the development of new therapies for metabolic disorders.

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.

Abstract Background Non-alcoholic fatty liver disease (NAFLD) affects 20–30% of the general adult population. NAFLD patients with type 2 diabetes mellitus (T2DM) are at an increased risk of advanced fibrosis, which puts them at risk of cardiovascular complications, hepatocellular carcinoma, or liver failure. Liver biopsy is the gold standard for assessing hepatic fibrosis. However, its utility is inherently limited. Consequently, the prevalence and characteristics of T2DM patients with advanced fibrosis are unknown. Therefore, the purpose of the current study is to evaluate the prevalence and severity of NAFLD in patients with T2DM by recruiting participants from primary care, using the latest imaging modalities, to collect a cohort of well phenotyped patients. Methods We will prospectively recruit 400 patients with T2DM using biomarkers to assess their status. Specifically, we will evaluate liver fat content using magnetic resonance imaging (MRI); hepatic fibrosis using MR elastography and vibration-controlled transient elastography; muscle composition and body fat distribution using water-fat separated whole body MRI; and cardiac function, structure, and tissue characteristics, using cardiovascular MRI. Discussion We expect that the study will uncover potential mechanisms of advanced hepatic fibrosis in NAFLD and T2DM and equip the clinician with better diagnostic tools for the care of T2DM patients with NAFLD. Trial registration: Clinicaltrials.gov, identifier NCT03864510. Registered 6 March 2019, https://clinicaltrials.gov/ct2/show/NCT03864510 .

Accurate measurements of metabolic fluxes in living cells are central to metabolism research and metabolic engineering. The gold standard method is model-based metabolic flux analysis (MFA), where fluxes are estimated indirectly from mass isotopomer data with the use of a mathematical model of the metabolic network. A critical step in MFA is model selection: choosing what compartments, metabolites, and reactions to include in the metabolic network model. Model selection is often done informally during the modelling process, based on the same data that is used for model fitting (estimation data). This can lead to either overly complex models (overfitting) or too simple ones (underfitting), in both cases resulting in poor flux estimates. Here, we propose a method for model selection based on independent validation data. We demonstrate in simulation studies that this method consistently chooses the correct model in a way that is independent on errors in measurement uncertainty. This independence is beneficial, since estimating the true magnitude of these errors can be difficult. In contrast, commonly used model selection methods based on the χ2-test choose different model structures depending on the believed measurement uncertainty; this can lead to errors in flux estimates, especially when the magnitude of the error is substantially off. We present a new approach for quantification of prediction uncertainty of mass isotopomer distributions in other labelling experiments, to check for problems with too much or too little novelty in the validation data. Finally, in an isotope tracing study on human mammary epithelial cells, the validation-based model selection method identified pyruvate carboxylase as a key model component. Our results argue that validation-based model selection should be an integral part of MFA model development.

Glucose homeostasis is the tight control of glucose in the blood. This complex control is important, due to its malfunction in serious diseases like diabetes, and not yet sufficiently understood. Due to the involvement of numerous organs and sub-systems, each with their own intra-cellular control, we have developed a multi-level mathematical model, for glucose homeostasis, which integrates a variety of data. Over the last 10 years, this model has been used to insert new insights from the intra-cellular level into the larger whole-body perspective. However, the original cell-organ-body translation has during these years never been updated, despite several critical shortcomings, which also have not been resolved by other modeling efforts. For this reason, we here present an updated multi-level model. This model provides a more accurate sub-division of how much glucose is being taken up by the different organs. Unlike the original model, we now also account for the different dynamics seen in the different organs. The new model also incorporates the central impact of blood flow on insulin-stimulated glucose uptake. Each new improvement is clear upon visual inspection, and they are also supported by statistical tests. The final multi-level model describes >300 data points in >40 time-series and dose-response curves, resulting from a large variety of perturbations, describing both intra-cellular processes, organ fluxes, and whole-body meal responses. We hope that this model will serve as an improved basis for future data integration, useful for research and drug developments within diabetes.

Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.

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.

The neurovascular coupling (NVC) connects neuronal activity to hemodynamic responses in the brain. This connection is the basis for the interpretation of functional magnetic resonance imaging data. Despite the central role of this coupling, we lack detailed knowledge about cell-specific contributions and our knowledge about NVC is mainly based on animal experiments performed during anesthesia. Anesthetics are known to affect neuronal excitability, but how this affects the vessel diameters is not known. Due to the high complexity of NVC data, mathematical modeling is needed for a meaningful analysis. However, neither the relevant neuronal subtypes nor the effects of anesthetics are covered by current models. Here, we present a mathematical model including GABAergic interneurons and pyramidal neurons, as well as the effect of an anesthetic agent. The model is consistent with data from optogenetic experiments from both awake and anesthetized animals, and it correctly predicts data from experiments with different pharmacological modulators. The analysis suggests that no downstream anesthetic effects are necessary if one of the GABAergic interneuron signaling pathways include a Michaelis-Menten expression. This is the first example of a quantitative model that includes both the cell-specific contributions and the effect of an anesthetic agent on the NVC.

Abstract Lumped parameter models of the cardiovascular system have the potential to assist researchers and clinicians to better understand cardiovascular function. The value of such models increases when they are subject specific. However, most approaches to personalize lumped parameter models have thus far required invasive measurements or fall short of being subject specific due to a lack of the necessary clinical data. Here, we propose an approach to personalize parameters in a model of the heart and the systemic circulation using exclusively non-invasive measurements. The personalized model is created using flow data from four-dimensional magnetic resonance imaging and cuff pressure measurements in the brachial artery. We term this personalized model the cardiovascular avatar. In our proof-of-concept study, we evaluated the capability of the avatar to reproduce pressures and flows in a group of eight healthy subjects. Both quantitatively and qualitatively, the model-based results agreed well with the pressure and flow measurements obtained in vivo for each subject. This non-invasive and personalized approach can synthesize medical data into clinically relevant indicators of cardiovascular function, and estimate hemodynamic variables that cannot be assessed directly from clinical measurements.

Estimation of liver function is important to monitor progression of chronic liver disease (CLD). A promising method is magnetic resonance imaging (MRI) combined with gadoxetate, a liver-specific contrast agent. For this method, we have previously developed a model for an average healthy human. Herein, we extended this model, by combining it with a patient-specific non-linear mixed-effects modeling framework. We validated the model by recruiting 100 patients with CLD of varying severity and etiologies. The model explained all MRI data and adequately predicted both timepoints saved for validation and gadoxetate concentrations in both plasma and biopsies. The validated model provides a new and deeper look into how the mechanisms of liver function vary across a wide variety of liver diseases. The basic mechanisms remain the same, but increasing fibrosis reduces uptake and increases excretion of gadoxetate. These mechanisms are shared across many liver functions and can now be estimated from standard clinical images.

Background: The possibility of non-invasively assessing load-independent parameters characterizing cardiac function is of high clinical value. Typically, these parameters are assessed during resting conditions. However, for diagnostic purposes, the parameter behavior across a physiologically relevant range of heart rate and loads is more relevant than the isolated measurements performed at rest. This study sought to evaluate changes in non-invasive estimations of load-independent parameters of left-ventricular contraction and relaxation patterns at rest and during dobutamine stress.Methods: We applied a previously developed approach that combines non-invasive measurements with a physiologically-based, reduced-order model of the cardiovascular system to provide subject-specific estimates of parameters characterizing left ventricular function. In this model, the contractile state of the heart at each time point along the cardiac cycle is modeled using a time-varying elastance curve. Non-invasive data, including four-dimensional magnetic resonance imaging (4D Flow MRI) measurements, were acquired in nine subjects without a known heart disease at rest and during dobutamine stress. For each of the study subjects, we constructed two personalized models corresponding to the resting and the stress state.Results: Applying the modeling framework, we identified significant increases in the left ventricular contraction rate constant [from 1.5 ± 0.3 to 2 ± 0.5 (p = 0.038)] and relaxation constant [from 37.2 ± 6.9 to 46.1 ± 12 (p = 0.028)]. In addition, we found a significant decrease in the elastance diastolic time constant from 0.4 ± 0.04 s to 0.3 ± 0.03 s (p = 0.008).Conclusions: The integrated image-modeling approach allows the assessment of cardiovascular function given as model-based parameters. The agreement between the estimated parameter values and previously reported effects of dobutamine demonstrates the potential of the approach to assess advanced metrics of pathophysiology that are otherwise difficult to obtain non-invasively in clinical practice.

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.

Abstract Intracellular signalling systems are highly complex. This complexity makes handling, analysis and visualisation of available knowledge a major challenge in current signalling research. Here, we present a novel framework for mapping signal‐transduction networks that avoids the combinatorial explosion by breaking down the network in reaction and contingency information. It provides two new visualisation methods and automatic export to mathematical models. We use this framework to compile the presently most comprehensive map of the yeast MAP kinase network. Our method improves previous strategies by combining (I) more concise mapping adapted to empirical data, (II) individual referencing for each piece of information, (III) visualisation without simplifications or added uncertainty, (IV) automatic visualisation in multiple formats, (V) automatic export to mathematical models and (VI) compatibility with established formats. The framework is supported by an open source software tool that facilitates integration of the three levels of network analysis: definition, visualisation and mathematical modelling. The framework is species independent and we expect that it will have wider impact in signalling research on any system.

In metabolic diseases such as Type 2 Diabetes and Non-Alcoholic Fatty Liver Disease, the systemic regulation of postprandial metabolite concentrations is disturbed. To understand this dysregulation, a quantitative and temporal understanding of systemic postprandial metabolite handling is needed. Of particular interest is the intertwined regulation of glucose and non-esterified fatty acids (NEFA), due to the association between disturbed NEFA metabolism and insulin resistance. However, postprandial glucose metabolism is characterized by a dynamic interplay of simultaneously responding regulatory mechanisms, which have proven difficult to measure directly. Therefore, we propose a mathematical modelling approach to untangle the systemic interplay between glucose and NEFA in the postprandial period. The developed model integrates data of both the perturbation of glucose metabolism by NEFA as measured under clamp conditions, and postprandial time-series of glucose, insulin, and NEFA. The model can describe independent data not used for fitting, and perturbations of NEFA metabolism result in an increased insulin, but not glucose, response, demonstrating that glucose homeostasis is maintained. Finally, the model is used to show that NEFA may mediate up to 30-45% of the postprandial increase in insulin-dependent glucose uptake at two hours after a glucose meal. In conclusion, the presented model can quantify the systemic interactions of glucose and NEFA in the postprandial state, and may therefore provide a new method to evaluate the disturbance of this interplay in metabolic disease.

OBJECTIVES: Diffuse liver disease (DLD), such as non-alcoholic fatty liver disease (NASH) and cirrhosis, is a rapidly growing problem throughout the Westernized world. Magnetic resonance imaging (MRI), based on uptake of the hepatocyte-specific contrast agent (CA) Gd-EOB-DTPA, is a promising non-invasive approach for diagnosing DLD. However, to fully utilize the potential of such dynamic measurements for clinical or research purposes, more advanced methods for data analysis are required. METHODS: A mathematical model that can be used for such data-analysis was developed. Data was obtained from healthy human subjects using a clinical protocol with high spatial resolution. The model is based on ordinary differential equations and goes beyond local diffusion modeling, taking into account the complete system accessible to the CA. RESULTS: The presented model can describe the data accurately, which was confirmed using chi-square statistics. Furthermore, the model is minimal and identifiable, meaning that all parameters were determined with small degree of uncertainty. The model was also validated using independent data. CONCLUSIONS: We have developed a novel approach for determining previously undescribed physiological hepatic parameters in humans, associated with CA transport across the liver. The method has a potential for assessing regional liver function in clinical examinations of patients that are suffering of DLD and compromised hepatic function.

Type 2 diabetes is characterized by insulin resistance of target organs, which is due to impaired insulin signal transduction. The skeleton of signaling mediators that provide for normal insulin action has been established. However, the detailed kinetics, and their mechanistic generation, remain incompletely understood. We measured time-courses in primary human adipocytes for the short-term phosphorylation dynamics of the insulin receptor (IR) and the IR substrate-1 in response to a step increase in insulin concentration. Both proteins exhibited a rapid transient overshoot in tyrosine phosphorylation, reaching maximum within 1 min, followed by an intermediate steady-state level after approximately 10 min. We used model-based hypothesis testing to evaluate three mechanistic explanations for this behavior: (A) phosphorylation and dephosphorylation of IR at the plasma membrane only; (B) the additional possibility for IR endocytosis; (C) the alternative additional possibility of feedback signals to IR from downstream intermediates. We concluded that (A) is not a satisfactory explanation; that (B) may serve as an explanation only if both internalization, dephosphorylation, and subsequent recycling are permitted; and that (C) is acceptable. These mechanistic insights cannot be obtained by mere inspection of the datasets, and they are rejections and thus stronger and more final conclusions than ordinary model predictions.