Background Stroke now represents the condition with the highest need for physical rehabilitation worldwide, with only low or moderate-level evidence testing telerehabilitation compared to in-person care. We compared functional ambulation in subacute patients with stroke following telerehabilitation and matched in-person controls with no biopsychosocial differences at baseline.Methods We conducted a matched case-control study to compare functional ambulation between individuals with stroke following telerehabilitation and in-person rehabilitation, assessed using the Functional Ambulation Categories (FAC) and the Functional Independence Measure (TM) (FIM).Results The telerehabilitation group (n = 38) achieved significantly higher FAC gains (1.5 (1.3) vs 1.0 (1.0)) than the in-person rehabilitation group, with no differences in ambulation efficiency, in individuals: admitted to rehabilitation within 60 days after stroke onset; aged 49.8 (+/- 11.4) years at admission; 55.3% female sex; moderate stroke severity; 42.1% with 'good' motor FIM at baseline; mostly living with sentimental partner (73.7%); with 21.1% holding an university education degree.Conclusions The groups showed no significant differences in ambulation efficiency, though the telerehabilitation group achieved higher FAC gains. Our results suggest that home telerehabilitation can be considered a good alternative to in-person rehabilitation when addressing ambulation in patients with moderate stroke severity and whose home situation mostly includes a cohabiting partner.

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 chi 2 -test (187.0 < T-chi 2 = 226.4) and independent validation data (70.8 < T-chi 2 = 93.5). The model can also be personalized using anthropometric data from a specific individual and can thus be used as a physiologically-based digital twin. This twin is also able to connect short-term consumption of alcohol to the long-term dynamics of PEth levels in the blood, a clinical biomarker of alcohol consumption. Here we illustrate how connecting short-term consumption to long-term markers allows for a new way to determine patient alcohol consumption from measured PEth levels. An additional use case of the twin could include the combined evaluation of patient-reported AUDIT forms and measured PEth levels. Finally, we integrated the new model into an eHealth application, which could help guide individual users or clinicians to help reduce dangerous drinking.

Background & aims: Non-alcoholic fatty liver disease (NAFLD) is a common metabolic disorder, characterized by the accumulation of excess fat in the liver, and is a driving factor for various severe liver diseases. These multi -factorial and multi-timescale changes are observed in different clinical studies, but these studies have not been integrated into a uni fied framework. In this study, we aim to present such a uni fied framework in the form of a dynamic mathematical model. Methods: For model training and validation, we collected data for dietary or drug -induced interventions aimed at reducing or increasing liver fat. The model was formulated using ordinary differential equations (ODEs) and the mathematical analysis, model simulation, model formulation and the model parameter estimation were all performed in MATLAB. Results: Our mathematical model describes accumulation of fat in the liver and predicts changes in lipid fluxes induced by both dietary and drug interventions. The model is validated using data from a wide range of drug and dietary intervention studies and can predict both short-term (days) and long-term (weeks) changes in liver fat. Importantly, the model computes the contribution of each individual lipid flux to the total liver fat dynamics. Furthermore, the model can be combined with an established bodyweight model, to simulate even longer scenarios (years), also including the effects of insulin resistance and body weight. To help prepare for corresponding eHealth applications, we also present a way to visualize the simulated changes, using dynamically changing lipid droplets, seen in images of liver biopsies. Conclusion: In conclusion, we believe that the minimal model presented herein might be a useful tool for future applications, and to further integrate and understand data regarding changes in dietary and drug induced changes in ectopic TAG in the liver. With further development and validation, the minimal model could be used as a disease progression model for steatosis. (c) 2024 The Author(s). Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

Liver metabolism is central to human physiology and influences the pathogenesis of common metabolic diseases. Yet, our understanding of human liver metabolism remains incomplete, with much of current knowledge based on animal or cell culture models that do not fully recapitulate human physiology. Here, we perform in-depth measurement of metabolism in intact human liver tissue ex vivo using global C-13 tracing, non-targeted mass spectrometry and model-based metabolic flux analysis. Isotope tracing allowed qualitative assessment of a wide range of metabolic pathways within a single experiment, confirming well-known features of liver metabolism but also revealing unexpected metabolic activities such as de novo creatine synthesis and branched-chain amino acid transamination, where human liver appears to differ from rodent models. Glucose production ex vivo correlated with donor plasma glucose, suggesting that cultured liver tissue retains individual metabolic phenotypes, and could be suppressed by postprandial levels of nutrients and insulin, and also by pharmacological inhibition of glycogen utilization. Isotope tracing ex vivo allows measuring human liver metabolism with great depth and resolution in an experimentally tractable system.

Introduction Cerebrovascular 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. Methods We 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 discussion Our 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.

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. A human-cell-based pancreas-liver microphysiological system serves as a preclinical platform for studying glucose-insulin homeostasis and disease mechanisms of glucose dysregulation, offering a tool for identifying targets and testing drugs.

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.

BackgroundThe 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.MethodsThe model is based on ordinary differential equations representing biochemical and physiological processes, in which unknown parameters were fitted to data using a MATLAB toolbox.ResultsThe 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.ConclusionsThe 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.

Insulin resistance causes compensatory insulin production, which in humans can eventually progress to beta-cell failure and type 2 diabetes (T2D). This disease progression involves multi-scale processes, ranging from intracellular signaling to organ and whole-body level regulations, on timescales from minutes to years. T2D progression is commonly studied using overfed and genetically modified rodents. Available multi-scale data from rodents is too complex to fully comprehend using traditional analysis, not based on mathematical modelling. To help resolve these issues, we here present an in silico mouse model, featuring 38 ordinary differential equations and 78 parameters. This is the first mathematical model that simultaneously explains (chi-square cost=28.1 <51 =cut-off, p = 0.05) multi-scale mouse insulin resistance data on all three levels - cells, organs, body - ranging from minutes to months. The model predicts new independent multi-scale simulations, on e.g., weight and meal response changes, which are corroborated by our own new experimental data. The thus validated model provides insights and non-trivial predictions regarding complex non-measured processes, such as the relation between insulin resistance and insulin-dependent glucose uptake for adipose tissue. Finally, we add a beta-cell failure module to the in silico mouse model to simulate different human-like scenarios of progression towards T2D. In summary, our in silico mouse model is an extendable and interactive knowledge-base for the study of T2D, which could help simulate treatment scenarios in rodents and translate results to the human situation.

Author summaryThe neurovascular coupling (NVC) is the basis for functional magnetic resonance imaging (fMRI), since the NVC connects neural activity with the observed hemodynamic changes. This connection is highly complex, which warrants a model-based analysis. However, even though NVC-data from several species and many relevant variables are available, a mathematical model for all these data is still missing. Herein, we combine experimental data from mice, monkeys, and humans, to develop a comprehensive model for NVC. Importantly, our new approach to modelling propagates the qualitative insights from each species to the subsequent analysis of data from other species. In mice, we unravel the role of different neuronal sub-populations when producing a biphasic response to prolonged sensory stimulations. The qualitative role of these sub-populations is preserved when analysing primate data. These primate data add knowledge on the interplay between local field potential (LFP) and vascular changes. Similarly, these pre-clinical qualitative insights are propagated to analysis of human data, which contain additional insights regarding blood flow and volume in arterioles and venules, during both positive and negative responses. This work illustrates how data with complementary information from different species can be combined, so that qualitative insights from animals are preserved in the quantitative analysis of human data. 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.

BackgroundLiver 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.MethodsIn 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.DiscussionThe 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 registrationClinicaltrials.gov, registration identifier NCT05502198. Registered on 16 August 2022. Link: https://classic.clinicaltrials.gov/ct2/show/NCT05502198.

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 trans-port, 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.

Type 2 diabetes (T2D) and hypertension increase the risk of cardiovascular diseases mediated by whole-body changes to metabolism, cardiovascular structure and haemodynamics. The haemodynamic changes related to hypertension and T2D are complex and subject-specific, however, and not fully understood. We aimed to investigate the haemodynamic mechanisms in T2D and hypertension by comparing the haemodynamics between healthy controls and subjects with T2D, hypertension, or both. For all subjects, we combined 4D flow magnetic resonance imaging data, brachial blood pressure and a cardiovascular mathematical model to create a comprehensive subject-specific analysis of central haemodynamics. When comparing the subject-specific haemodynamic parameters between the four groups, the predominant haemodynamic difference is impaired left ventricular relaxation in subjects with both T2D and hypertension compared to subjects with only T2D, only hypertension and controls. The impaired relaxation indicates that, in this cohort, the long-term changes in haemodynamic load of co-existing T2D and hypertension cause diastolic dysfunction demonstrable at rest, whereas either disease on its own does not. However, through subject-specific predictions of impaired relaxation, we show that altered relaxation alone is not enough to explain the subject-specific and group-related differences; instead, a combination of parameters is affected in T2D and hypertension. These results confirm previous studies that reported more adverse effects from the combination of T2D and hypertension compared to either disease on its own. Furthermore, this shows the potential of personalized cardiovascular models in providing haemodynamic mechanistic insights and subject-specific predictions that could aid in the understanding and treatment planning of patients with T2D and hypertension.

Inflammation is one of the vital mechanisms through which the immune system responds to harmful stimuli. During inflammation, proinflammatory and anti-inflammatory cytokines interplay to orchestrate fine-tuned and dynamic immune responses. The cytokine interplay governs switches in the inflammatory response and dictates the propagation and development of the inflammatory response. Molecular pathways underlying the interplay are complex, and time-resolved monitoring of mediators and cytokines is necessary as a basis to study them in detail. Our understanding can be advanced by mathematical models that enable to analyze the system of interactions and their dynamical interplay in detail. We, therefore, used a mathematical modeling approach to study the interplay between prominent proinflammatory and anti-inflammatory cytokines with a focus on tumor necrosis factor and interleukin 10 (IL-10) in lipopolysaccharide-primed primary human monocytes. Relevant time-resolved data were generated by experimentally adding or blocking IL-10 at different time points. The model was successfully trained and could predict independent validation data and was further used to perform simulations to disentangle the role of IL-10 feedbacks during an acute inflammatory event. We used the insight to obtain a reduced predictive model including only the necessary IL-10-mediated feedbacks. Finally, the validated reduced model was used to predict early IL-10–tumor necrosis factor switches in the inflammatory response. Overall, we gained detailed insights into fine-tuning of inflammatory responses in human monocytes and present a model for further use in studying the complex and dynamic process of cytokine-regulated acute inflammation. © 2023 The Authors

BACKGROUND: Drug 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.

AIMS: Our 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 RESULTS: We 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.

CONCLUSION: Our 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.

Adipocyte cellular signaling, normally and in type 2 diabetes, is far from fully studied. We have earlier developed detailed dynamic mathematical models for some well-studied, and 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 is key. There exists such data for the insulin response of adipocytes, as well as prior knowledge on possible protein-protein interactions associated with a confidence level. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. In our new method, we first establish a core model by connecting our partially overlapping models of adipocyte cellular signaling with focus on: 1) lipolysis and fatty acid release, 2) glucose uptake, and 3) the release of adiponectin. We use the phosphoproteome data and prior knowledge to identify phosphosites adjacent to the core model, and then try to add the adjacent phosphosites to the model. The additions of the adjacent phosphosites is tested in a parallel, pairwise approach with low computation time. We then iteratively collect the accepted additions into a layer, and use the newly added layer to find new adjacent phosphosites. We find that the first 15 layers (60 added phosphosites) with the highest confidence can correctly predict independent inhibitor-data (70-90 % correct), and that this ability decrease when we add layers of decreasing confidence. In total, 60 layers (3926 phosphosites) can be added to the model and still keep predictive ability. Finally, we use the comprehensive adipocyte model to simulate systems-wide alterations in adipocytes in type 2 diabetes. This new method provide a tool to create large models that keeps track of varying confidence.Competing Interest StatementThe authors have declared no competing interest.

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 individuals sex, weight, height, as well as to the individuals historical data on metabolite dynamics. This tool enables an offline digital twin technology.

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.

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. Author summary The neurovascular and neurometabolic couplings are highly central for several clinical imaging techniques since these frequently use blood oxygenation (the BOLD signal) as a proxy for neuronal activity. This relationship is described by the highly complex neurovascular and neurometabolic couplings, which describe the balancing between increased metabolic demand and blood flow, and which involve several cell types and regulatory systems, which all change dynamically over time. While there are previous works that describe the neurovascular coupling in detail, neither we nor others have developed connections to corresponding mechanistic models for the third aspect, the metabolic aspect. Furthermore, magnetic resonance spectroscopy (MRS) data for such modelling readily is available. In this paper we present a minimal mechanistic model that can describe the metabolic response to visual stimuli. The model is trained to describe experimental data for the relative change in metabolic concentrations of several metabolites in the visual cortex during stimulation. The model is also validated against independent validation data, that was not used for model training. Finally, we also connect this metabolic model to a detailed mechanistic model of the neurovascular coupling. Showing that the model can describe both the metabolic response and a neurovascular response simultaneously.

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 chi(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.

Circulating levels of the adipocyte hormone adiponectin are typically reduced in obesity, and this deficiency has been linked to metabolic diseases. It is thus important to understand the mechanisms controlling adiponectin exocytosis. This understanding is hindered by the high complexity of both the available data and the underlying signaling network. To deal with this complexity, we have previously investigated how different intracellular concentrations of Ca2+, cAMP, and ATP affect adiponectin exocytosis, using both patch-clamp recordings and systems biology mathematical modeling. Recent work has shown that adiponectin exocytosis is physiologically triggered via signaling pathways involving adrenergic beta(3) receptors (beta(3)ARs). Therefore, we developed a mathematical model that also includes adiponectin exocytosis stimulated by extracellular epinephrine or the beta(3)AR agonist CL 316243. Our new model is consistent with all previous patch-clamp data as well as new data (collected from stimulations with a combination of the intracellular mediators and extracellular adrenergic stimuli) and can predict independent validation data. We used this model to perform new in silico experiments where corresponding wet lab experiments would be difficult to perform. We simulated adiponectin exocytosis in single cells in response to the reduction of beta(3)ARs that is observed in adipocytes from animals with obesity-induced diabetes. Finally, we used our model to investigate intracellular dynamics and to predict both cAMP levels and adiponectin release by scaling the model from single-cell to a population of cells-predictions corroborated by experimental data. Our work brings us one step closer to understanding the intricate regulation of adiponectin exocytosis.

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.

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.

BackgroundNon-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.MethodsWe 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.DiscussionWe 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.

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.

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.

BACKGROUND AND AIMS: Accurate biomarkers for quantifying liver fibrosis are important for clinical practice and trial end-points. We compared the diagnostic performance of magnetic resonance imaging (MRI), including gadoxetate-enhanced MRI and ³¹P-MR spectroscopy, with fibrosis stage and serum fibrosis algorithms in a clinical setting. Also, in a subset of patients, MR- and transient elastography (MRE and TE) was evaluated when available.

METHODS: Patients were recruited prospectively if they were scheduled to undergo liver biopsy on a clinical indication due to elevated liver enzyme levels without decompensated cirrhosis. Within a month of the clinical work-up, an MR-examination and liver needle biopsy were performed on the same day. Based on late-phase gadoxetate-enhanced MRI, a mathematical model calculated hepatobiliary function (relating to OATP1 and MRP2). The hepatocyte gadoxetate uptake rate (K_Hep) and the normalised liver-to-spleen contrast ratio (LSC_N10) were also calculated. Nine serum fibrosis algorithms were investigated (GUCI, King's Score, APRI, FIB-4, Lok-Index, NIKEI, NASH-CRN regression score, Forns' score, and NAFLD-fibrosis score).

RESULTS: The diagnostic performance (AUROC) for identification of significant fibrosis (F2-4) was 0.78, 0.80, 0.69, and 0.78 for MRE, TE, LSC_N10, and GUCI, respectively. For the identification of advanced fibrosis (F3-4), the AUROCs were 0.93, 0.84, 0.81, and 0.82 respectively.

CONCLUSION: MRE and TE were superior for non-invasive identification of significant fibrosis. Serum fibrosis algorithms developed for specific liver diseases are applicable in this cohort of diverse liver diseases aetiologies. Gadoxetate-MRI was sufficiently sensitive to detect the low function losses associated with fibrosis. None was able to efficiently distinguish between stages within the low fibrosis stages.Lay summaryExcessive accumulation of scar tissue, fibrosis, in the liver is an important aspect in chronic liver disease. To replace the invasive needle biopsy, we have explored non-invasive methods to assess liver fibrosis. In our study we found that elastographic methods, which assess the mechanical properties of the liver, are superior in assessing fibrosis in a clinical setting. Of interest from a clinical trial point-of-view, none of the tested methods was sufficiently accurate to distinguish between adjacent moderate fibrosis stages.

Both initiation and suppression of inflammation are hallmarks of the immune response. If not balanced, the inflammation may cause extensive tissue damage, which is associated with common diseases, e.g., asthma and atherosclerosis. Anti-inflammatory drugs come with side effects that may be aggravated by high and fluctuating drug concentrations. To remedy this, an anti-inflammatory drug should have an appropriate pharmacokinetic half-life or better still, a sustained anti-inflammatory drug response. However, we still lack a quantitative mechanistic understanding of such sustained effects. Here, we study the anti-inflammatory response to a common glucocorticoid drug, dexamethasone. We find a sustained response 22 hours after drug removal. With hypothesis testing using mathematical modeling, we unravel the underlying mechanism-a slow release of dexamethasone from the receptor-drug complex. The developed model is in agreement with time-resolved training and testing data and is used to simulate hypothetical treatment schemes. This work opens up for a more knowledge-driven drug development to find sustained anti-inflammatory responses and fewer side effects.

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.

Author summary

Being able to accurately and reliably estimate liver function is important when monitoring the progression of patients with liver disease, as well as when identifying drug-induced liver injury during drug development. A promising method for quantifying liver function is to use magnetic resonance imaging combined with gadoxetate. Gadoxetate is a liver-specific contrast agent, which is taken up by the hepatocytes and excreted into the bile. We have previously developed a mechanistic model for gadoxetate dynamics using averaged data from healthy volunteers. In this work, we extended our model with a non-linear mixed-effects modeling framework to give patient-specific estimates of the gadoxetate transport-rates. We validated the model by recruiting 100 patients with liver disease, covering a range of severity and etiologies. All patients underwent an MRI-examination and provided both blood and liver biopsies. Our validated model provides a new and deeper look into how the mechanisms of liver function varies across a wide variety of liver diseases. The basic mechanisms remain the same, but increasing fibrosis reduces uptake and increases excretion of gadoxetate.

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

Type 2 diabetes is characterized by insulin resistance, which arises from malfunctions in the intracellular insulin signaling network. Knowledge of the insulin signaling network is fragmented, and because of the complexity of this network, little consensus has emerged for the structure and importance of the different branches of the network. To help overcome this complexity, systems biology mathematical models have been generated for predicting both the activation of the insulin receptor (IR) and the redistribution of glucose transporter 4 (GLUT4) to the plasma membrane. Although the insulin signal transduction between IR and GLUT4 has been thoroughly studied with modeling and time-resolved data in human cells, comparable analyses in cells from commonly used model organisms such as rats and mice are lacking. Here, we combined existing data and models for rat adipocytes with new data collected for the signaling network between IR and GLUT4 to create a model also for their interconnections. To describe all data (amp;gt;140 data points), the model needed three distinct pathways from IR to GLUT4: (i) via protein kinase B (PKB) and Akt substrate of 160 kDa (AS160), (ii) via an AS160-independent pathway from PKB, and (iii) via an additional pathway from IR, e.g. affecting the membrane constitution. The developed combined model could describe data not used for training the model and was used to generate predictions of the relative contributions of the pathways from IR to translocation of GLUT4. The combined model provides a systems-level understanding of insulin signaling in rat adipocytes, which, when combined with corresponding models for human adipocytes, may contribute to model-based drug development for diabetes.

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.

The molecular mechanisms of insulin resistance in Type 2 diabetes have been extensively studied in primary human adipocytes, and mathematical modelling has clarified the central role of attenuation of mammalian target of rapamycin (mTOR) complex 1 (mTORC1) activity in the diabetic state. Attenuation of mTORC1 in diabetes quells insulin-signalling network-wide, except for the mTOR in complex 2 (mTORC2)-catalysed phosphorylation of protein kinase B (PKB) at Ser(473) (PKB-S473P), which is increased. This unique increase could potentially be explained by feedback and interbranch cross-talk signals. To examine if such mechanisms operate in adipocytes, we herein analysed data from an unbiased phosphoproteomic screen in 3T3-L1 adipocytes. Using a mathematical modelling approach, we showed that a negative signal from mTORC1-p70 S6 kinase (S6K) to rictor-mTORC2 in combination with a positive signal from PKB to SIN1-mTORC2 are compatible with the experimental data. This combined cross-branch signalling predicted an increased PKB-S473P in response to attenuation of mTORC1 - a distinguishing feature of the insulin resistant state in human adipocytes. This aspect of insulin signalling was then verified for our comprehensive model of insulin signalling in human adipocytes. Introduction of the cross-branch signals was compatible with all data for insulin signalling in human adipocytes, and the resulting model can explain all data network-wide, including the increased PKB-S473P in the diabetic state. Our approach was to first identify potential mechanisms in data from a phosphoproteomic screen in a cell line, and then verify such mechanisms in primary human cells, which demonstrates how an unbiased approach can support a direct knowledge-based study.

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.

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 naive 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) of hemodynamic changes captured in the blood oxygen level-dependent (BOLD) response contains information of brain activity. The BOLD response is the result of a complex neurovascular coupling and comes in at least two fundamentally different forms: a positive and a negative deflection. Because of the complexity of the signaling, mathematical modelling can provide vital help in the data analysis. For the positive BOLD response, there are plenty of mathematical models, both physiological and phenomenological. However, for the negative BOLD response, no physiologically based model exists. Here, we expand our previously developed physiological model with the most prominent mechanistic hypothesis for the negative BOLD response: the neural inhibition hypothesis. The model was trained and tested on experimental data containing both negative and positive BOLD responses from two studies: 1) a visual-motor task and 2) a workin-gmemory task in conjunction with administration of the tranquilizer diazepam. Our model was able to predict independent validation data not used for training and provides a mechanistic underpinning for previously observed effects of diazepam. The new model moves our understanding of the negative BOLD response from qualitative reasoning to a quantitative systems-biology level, which can be useful both in basic research and in clinical use.

Pharmaceutical induction of metabolically active beige adipocytes in the normally energy storing white adipose tissue has potential to reduce obesity. Mitochondrial uncoupling in beige adipocytes, as in brown adipocytes, has been reported to occur via the uncoupling protein 1 (UCP1). However, several previous in vitro characterizations of human beige adipocytes have only measured UCP1 mRNA fold increase, and assumed a direct correlation with metabolic activity. Here, we provide an example of pharmaceutical induction of beige adipocytes, where increased mRNA levels of UCP1 are not translated into increased protein levels, and perform a thorough analysis of this example. We incorporate mRNA and protein levels of UCP1, time-resolved mitochondrial characterizations, and numerous perturbations, and analyze all data with a new fit-for-purpose mathematical model. The systematic analysis challenges the seemingly obvious experimental conclusion, i.e., that UCP1 is not active in the induced cells, and shows that hypothesis testing with iterative modeling and experimental work is needed to sort out the role of UCP1. The analyses demonstrate, for the first time, that the uncoupling capability of human beige adipocytes can be obtained without UCP1 activity. This finding thus opens the door to a new direction in drug discovery that targets obesity and its associated comorbidities. Furthermore, the analysis advances our understanding of how to evaluate UCP1-independent thermogenesis in human beige adipocytes.

It is possible to estimate hepatic triglyceride content by calculating the proton density fat fraction (PDFF), using proton magnetic resonance spectroscopy (less thansuperscriptgreater than1less than/superscriptgreater thanH-MRS), instead of collecting and analyzing liver biopsies to detect steatosis. However, the current PDFF cut-off value (5%) used to define steatosis by magnetic resonance was derived from studies that did not use histopathology as the reference standard. We performed a prospective study to determine the accuracy of less thansuperscriptgreater than1less than/superscriptgreater thanH-MRS PDFF in measurement of steatosis using histopathology analysis as the standard. We collected clinical, serologic, less thansuperscriptgreater than1less than/superscriptgreater thanH-MRS PDFF, and liver biopsy data from 94 adult patients with increased levels of liver enzymes (6 months or more) referred to the Department of Gastroenterology and Hepatology at Linköping University Hospital in Sweden from 2007 through 2014. Steatosis was graded using the conventional histopathology method and fat content was quantified in biopsy samples using stereological point counts (SPCs). We correlated less thansuperscriptgreater than1less than/superscriptgreater thanH-MRS PDFF findings with SPCs (r = 0.92; P less than.001). less thansuperscriptgreater than1less than/superscriptgreater thanH-MRS PDFF results correlated with histopathology results (ρ = 0.87; P less than.001), and SPCs correlated with histopathology results (ρ = 0.88; P less than.001). All 25 subjects with PDFF values of 5.0% or more had steatosis based on histopathology findings (100% specificity for PDFF). However, of 69 subjects with PDFF values below 5.0% (negative result), 22 were determined to have steatosis based on histopathology findings (53% sensitivity for PDFF). Reducing the PDFF cut-off value to 3.0% identified patients with steatosis with 100% specificity and 79% sensitivity; a PDFF cut-off value of 2.0% identified patients with steatosis with 94% specificity and 87% sensitivity. These findings might be used to improve non-invasive detection of steatosis.

Introduction

Over the past decades, several different non-invasive methods for staging hepatic fibrosis have been proposed. One such method is dynamic contrast enhanced MRI (DCE-MRI) using the contrast agent (CA) Gd-EOB-DTPA. Gd-EOB-DTPA is liver specific, which means that it is taken up specifically by the hepatocytes via the OATP3B1/B3 transporters and excreted into the bile via the MRP2 transporter. Several studies have shown that DCE-MRI and Gd-EOBDTPA can separate patients with advanced (F3-F4) from mild (F0-F2) hepatic fibrosis by measuring the signal intensity, where patients with advanced fibrosis have a lower signal intensity than the mild fibrosis cases.1 However, none of the studies up to date have been able to differentiate if the reduced signal intensity in the liver is because of an decreased uptake of CA or an increased excretion. Analyzing the DCE-MRI data with mechanistic mathematical modelling has the possibility of investigating such a differentiation.

Subjects and methods

88 patients with diffuse liver disease were examined using DCE-MRI (1.5 T Philips Achieva, two-point Dixon, TR=6.5 ms, TE=2.3/4.6 ms, FA=13) after a bolus injection of Gd-EOB-DTPA, followed by a liver biopsy. Regions of interest were placed within the liver, spleen and veins and a whole-body mechanistic pharmacokinetic model2 was fitted to the data. The fitted parameters in the model correspond to the rate of CA transport between different compartments, e.g. hepatocytes, blood plasma, and bile (Fig. 1).

Results

As can be seen in Fig. 2, the parameter corresponding to the transport of CA from the blood plasma to the hepatocytes, kph, is lower for patients with advanced fibrosis (p=0.01). Fig. 3 shows that the parameter corresponding to the CA excretion into the bile, khb, is higher for patients with advanced fibrosis (p<0.01).

Discussion/Conclusion

This work shows that the decreased signal intensity in DCE-MRI images in patients with advanced fibrosis depends on both a decreased uptake of CA in the hepatocytes and an increased excretion into the bile. Similar results have also been observed in a rat study3. In that study, rats with induced cirrhosis had a higher MRP2-activity than the healthy control rats.

References

1Norén et al: Eur. Radiol, 23(1), 174-181, 2013.

2Forsgren et al: PloS One, 9(4): e95700, 2014.

3Tsuda & Matsui: Radiol, 256(3): 767-773, 2010.

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.

Parameter estimation – the assignment of values to the parameters in a model – is an important and time-consuming task in computational biology. Recent computational and algorithmic developments have provided novel tools to improve this estimation step. One of these improvements concerns the optimization step, where the parameter space is explored to find interesting regions. In this chapter we review the parameter estimation problem, with a special emphasis on the associated optimization methods. In relation to this, we also provide concepts and tools to help you select the appropriate methodology for a specific scenario.

One of the most important properties of a mathematical model is the abilityto make predictions: to predict that which has not yet been measured. Suchpredictions can sometimes be obtained from a simple simulation, but that requiresthat the parameters in the model are known from before. In biology, theparameters are usually both not known from before and not identifiable, i.e.the parameter values cannot be determined uniquely from available data. Insuch cases of unidentifiability, the space of acceptable parameters is large, ofteninfinite in certain directions. For such large spaces, sampling-based approachesthat try to characterize the entire space have difficulties. Recently, a new type ofalternative approaches that circumvent this characterization problem has beenproposed: where one only searches those directions in the space of acceptable parametersthat are relevant for the uncertainty of a particular prediction. In thisreview chapter, these recently proposed methods are compared and contrasted,both regarding theoretical properties, and regarding user experience. The focusis on methods from the field of systems biology, but also methods from biostatistics,pharmacodynamics, and biochemometrics are discussed. The hope is thatthis review will increase the usefulness and understanding of already proposedmethods, and thereby help foster a tradition where predictions only are deemedinteresting if their uncertainties have been determined.

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.