Computational Modeling of Dynamical Liver Function Tests: LiMAx
König M., Wünsch T., Lill D., Lock JF., Mohr J., Taheri P., Timmer J., Stockmann M.
Computational Modeling of Dynamical Liver Function Tests: LiMAx and Methacetin Breath Test (MBT)
Modeling LiMAx liver function test
Assessment of liver function is a key task in hepatology but accurate quantification of hepatic function has remained a clinical challenge. Dynamic liver function tests are a promising tool for the non-invasive evaluation of liver function in vivo. One class of such tests are breath tests based on the conversion of 13C-labeled substrates by the liver to 13CO2 subsequently measured in the breath. A commonly applied substrate is 13C-methacetin, converted to paracetamol and 13CO2 via cytochrome P450 1A2 (CYP1A2), used orally in the methacetin breath test (MBT) and intravenously in the LiMAx test. An important clinical question is which factors can affect MBT and LiMAx results. The aim of our study was to answer this question using computational modeling to derive basic information for a better understanding of the methacetin breath test and factors influencing its results.
A physiological based pharmacokinetics (PBPK) model for 13C-methacetin breath tests including absorption, distribution, metabolism and elimination of 13C-methacetin, paracetamol, and 13C-bicarbonate/13CO2 was developed. The model correctly predicts data from more than 20 clinical studies after oral and intravenous application under various dosing regimens of paracetamol, 13C-bicarbonate, and 13C-methacetin, based on retrospective analysis. Model predictions were validated based on data from multiple studies consisting of 13C-bicarbonate kinetics, LiMAx dosing study and retrospective analysis of LiMAx data from healthy smokers and non-smokers. Sensitivity analysis was performed to identify key factors influencing MBT and LiMAx.The model was applied to study the effect of clinically relevant parameters like CYP1A2 content, hepatic perfusion or lifestyle factors like smoking on LiMAx. In summary, we present a valuable tool for the evaluation of dynamical liver function tests based on 13C-methacetin.
Modeling Hepatic Central Metabolism
Berndt N., Bulik S., Wallach I., Wünsch T., König M., Stockmann M., Meierhofer M., Holzhütter HG.
HEPATOKIN1: A Biochemistry-Based Model of Liver Metabolism Suited for Applications in Medicine and Pharmacology
Nat Commun. 2018 Jun 19;9(1):2386. doi: 10.1038/s41467-018-04720-9.
Modeling hepatic metabolism in health and disease
The epidemic increase of non-alcoholic fatty liver diseases (NAFLD) requires a deeper understanding of the regulatory circuits controlling the response of liver metabolism to nutritional challenges, medical drugs, and genetic enzyme variants. As in vivo studies of human liver metabolism are encumbered with serious ethical and technical issues, we were involved in the development of a comprehensive biochemistry-based kinetic model of the central liver metabolism including the regulation of enzyme activities by their reactants, allosteric effectors, and hormone-dependent phosphorylation.
Stratification and personalization of computational models
Workflow for extraction of pharmacokinetics data
Physiologically based pharmacokinetics (PBPK) models are an important tool for understanding dynamical liver function tests. Such models allow to simulate the absorption, distribution, metabolization and elimination (ADME) of test substances. Stratification and personalization of PBPK models is an important task in systems biology and systems medicine, e.g., to simulate different subpopulations like smokers vs. non-smokers, obese vs. lean people, or children vs. adults vs. elderly, or to simulate the individual person based for instance on anthropometric information like age, gender or bodyweight. Such PBPK models can be described in standard exchange formats like the systems biology markup language (SBML) which allow a reproducible, exchangable and interoperable encoding of these models. On the other hand data sets for stratification and personalization of such models are available from the research literature or clinical cooperation partners. An open task to implement an workflow for the stratification and personalization of PBPK models from existing datasets. This includes the automatic extraction of information via literature mining and machine learning, storage of information in a pharmacokinetics database, and parametrization of models base on subsets of information.
Multi-scale Model of Hepatic Metabolism
Multiscale model of human galactose metabolism - from single hepatocytes to individual liver function
Understanding how liver function arises from the complex interaction of morphology, perfusion, and metabolism from single cells up to the entire organ requires systems-levels computational approaches. We report a multiscale mathematical model of the Human liver comprising the scales from single hepatocytes, over representation of ultra-structure and micro-circulation in the hepatic tissue, up to the entire organ integrated with perfusion. The model was validated against data on multiple spatial and temporal scales. Herein we describe the model construction and application to hepatic galactose metabolism demonstrating its utility via i) the personalization of liver function tests based on galactose elimination capacity (GEC), ii) the explanation of changes in liver function with aging, and iii) the prediction of population variability in liver function based on variability in liver volume and perfusion. We conclude that physiology- and morphology-based multiscale models can improve the evaluation of individual liver function.
Liver Function Tests
Personalized evaluation of liver function tests
Understanding how liver function results from the complex interaction of hepatic morphology, perfusion and metabolism across multiple spatial scales, from the cellular level up to the entire organ, requires systems levels computational approaches. Despite the importance of quantifying liver function in the clinics, currently no computational model of the human liver is able to utilize the morphological and physiological information from imaging methods for a personalized evaluation of quantitative function tests. Our objectives are
- providing an open, freely available and reproducible multiscale model of the human liver which can readily be employed to evaluate the hepatic clearance of various test substances used in quantitative liver function tests in the clinics
- providing the infrastructure and methods to integrate experimental and clinical data on all scales (imaging data, anthropomorphic information, perfusion data multipleindicator dilution data, microcirculation and ultrastructure parameters, protein amounts, ...)
- create metabolic modules for the elimination of various test substances and drugs (methacetin, galactose, BSP, ICG, midazolam, codeine, caffeine, torsemide, pravastatin, talinolol) by the liver. This will allow the analysis of a wide range of clinical studies and create an integrative platform for the clearance of these substances by the liver.
Biomarker signatures in liver disease
Kerstin A.*, König M.*, Hoppe A., Thomas M., Müller I., Ebert M., Weng H., Holzhütter HG., Zanger UM., Bode J., Vollmar B. and Dooley S.
* equal contribution
Pathobiochemical signatures of cholestatic liver disease in bile duct ligated mice
BMC Syst Biol. 2015 Nov 20;9(1):83. doi: 10.1186/s12918-015-0229-0 [PubMed]
Pathobiochemical signatures of cholestatic liver disease
Disrupted bile secretion leads to liver damage characterized by inflammation, fibrosis, eventually cirrhosis, and hepatocellular cancer. As obstructive cholestasis often progresses insidiously, markers for the diagnosis and staging of the disease are urgently needed. To this end, we compiled a comprehensive data set of serum markers, histological parameters and transcript profiles at 8 time points of disease progression after bile duct ligation (BDL) in mice, aiming at identifying a set of parameters that could be used as robust biomarkers for transition of different disease progression phases.
High Performance Simulators and Simulation Software
Medley K., Choi K., König M., Smith L., Gu S., Hellerstein J., Sealfon S., and Sauro HM.
Tellurium Notebooks — An Environment for Reproducible Dynamical Modeling in Systems Biology.
PLoS computational biology 14, no. 6 (2018): e1006220.; doi: 10.1371/journal.pcbi.1006220
Tellurium Notebooks - An Environment for Dynamical Model Development, Reproducibility, and Reuse
The considerable difficulty encountered in reproducing the results of published dynamical models limits validation, exploration and reuse of this increasingly large biomedical research resource. To address this problem, we have developed Tellurium Notebook, a software system that facilitates building reproducible dynamical models and reusing models by 1) supporting the COMBINE archive format during model development for capturing model information in an exchangeable format and 2) enabling users to easily simulate and edit public COMBINE-compliant models from public repositories to facilitate studying model dynamics, variants and test cases. Tellurium Notebook, a Python--based Jupyter--like environment, is designed to seamlessly inter-operate with these community standards by automating conversion between COMBINE standards formulations and corresponding in--line, human--readable representations. Thus, Tellurium brings to systems biology the strategy used by other literate notebook systems such as Mathematica. These capabilities allow users to edit every aspect of the standards--compliant models and simulations, run the simulations in--line, and re--export to standard formats. We provide several use cases illustrating the advantages of our approach and how it allows development and reuse of models without requiring technical knowledge of standards. Adoption of Tellurium should accelerate model development, reproducibility and reuse.
Reproducible simulation experiments
tellurium-web: An online database for reproducible simulation experiments in computational biology
We present tellurium-web, a free and open online database and simulation tool for reproducible simulation experiments in computational biology. Tellurium-web is build on top of existing standards in computational biology, i.e., the SBML for models, SED-ML for description of simulation experiments and the COMBINE archive format for exchange of simulation studies. Key features are 1) large set of example archives (e.g., simulation studies from biomodels, SED-ML specification examples, SBML testsuite, JWS online, tellurium examples); 2) COMBINE archive browser with access to metadata and content of the archives; 3) Execution of archives using a high-performance simulator (libroadrunner) integrated into modeling platform for computational biology (tellurium); 4) REST web services for getting, posting and modifying archives providing means for simple integration into computational workflows; and 4) filtering of archives based on tags and users.
FluType: Influenza Classification based on Machine Learning
Development of an influenza sub-classification platform based on peptides
The FluTypeDB project is a web application for the data management of binding assays for the classification of influenza subtypes. FluTypeDB consists of a database and web interface with focus on various binding assays for the classification of influenza viruses and contains experimental data based on
- ELISA antibody binding assays
- Peptid binding assays on microwell plates
- Peptid bindings assays on peptide microarrays
Visualization of kinetic data for modeling
cy3sabiork: A Cytoscape app for visualizing kinetic data from SABIO-RK
[version 1; referees: awaiting peer review]. F1000Research 2016, 5:1736 (doi: 10.12688/f1000research.9211.1)
Cytoscape app for visualizing kinetic data from SABIO-RK
Kinetic data of biochemical reactions are essential for the creation of kinetic models of biochemical networks. One of the main resources of such information is SABIO-RK, a curated database for kinetic data of biochemical reactions and their related information. Despite the importance for computational modelling there has been no simple solution to visualize the kinetic data from SABIO-RK. In this work, I present cy3sabiork, an app for querying and visualization of kinetic data from SABIO-RK in Cytoscape. The kinetic information is accessible via a combination of graph structure and annotations of nodes, with provided information consisting of: (I) reaction details, enzyme and organism; (II) kinetic law, formula, parameters; (III) experimental conditions; (IV) publication; (V) additional annotations. cy3sabiork creates an intuitive visualization of kinetic entries in form of a species-reaction-kinetics graph, which reflects the reaction-centered approach of SABIO-RK. Kinetic entries can be imported in SBML format from either the SABIO-RK web interface or via web service queries. The app allows for easy comparison of kinetic data, visual inspection of the elements involved in the kinetic record and simple access to the annotation information of the kinetic record. The cy3sabiork approach was applied in the computational modelling of galactose metabolism in the human liver.