Zake2021 - PBPK model of metformin in mice: single dose intavenous
View the 2021-12 Model of the Month entry for this model
Model Identifier
BIOMD0000001039
Short description
This model is supplementary material of publication "Physiologically based metformin pharmacokinetics model of mice and scale-up to humans for the estimation of concentrations in various tissues"
by Darta Maija Zake, Linda Zaharenko, JanisKurlovics, Vitalijs Komasilovs, Egils Stalidzans and Janis Klovins.
This is a whole-body model representing the pharmacokinetics of metformin in the mouse body. The model is in the form of ordinary differential equations and describes metformin concentration in 20 compartments.
The model consists of 20 compartments (“Compartments” in COPASI model) describing various tissues or tissue sub-compartments and body fluids of metformin action (venous and arterial plasma, intestine, kidney, heart, fat, muscle, brain, lungs, stomach, liver, portal vein, remainder urine and feces). Body weight and the weight of all compartments is expressed as a volume in mL and for the calculations it is assumed that 1mL = 1g. The volumes of most compartments are calculated as a fraction of the body weight/volume, and the fractions are determined from literature data, the volumes of the stomach lumen and intestine lumen are fixed and do not change depending on the body weight. Similarly, the volume of external urine and feces is set to 1mL, but those are “volumeless” compartments as they are only necessary for the calculation of metformin amount, not concentration.
The model consists of 20 species (“Species” in COPASI model) that correspond to the metformin concentrations in the 20 compartments. The initial concentrations for all the species are 0 nmol/mL as metformin is not produced in the body and can only be detected after dose administration.
The model consists of 33 reactions – they describe the transport processes of metformin in the body. The reactions include local parameters that are involved only in that particular reaction and global parameters – parameters that are used in multiple reactions or are calculated depending on another parameter e.g. scale-up coefficients.
The model consists of 52 global quantities – parameters involved in multiple reactions or necessary for another parameter calculation:
1.Parameters describing metformin dose – either in peroral (Metformin Dose in Lumen in mg) or intravenous (Metformin Dose in Plasma in mg).
2.Parameter describing mice physiology – body weight (in mL), cardiac output, blood flow to different compartments described as Q”compartment_name” (for example Qliver describes blood flow to the liver compartment). Qgfr refers to the glomerular filtration rate.
3.Tissue:plasma partition coefficients (Ktp) that are necessary for the scale-up to humans.
4.Parameters involved in the calculation of metformin amount in mg, these parameters are named mg”Compartment_name” (for example mgLiver describes the metformin amount in mg in the liver tissues).
The time points of dose release are defined as “events” in COPASI and can be changed as necessary.
Time course simulations can be accessed through the section “Time Course” in this section the time duration and intervals can be changed. When time-course simulations are run three plots are created – Metformin amount in the 20 compartments, metformin concentrations in the compartments and reaction fluxes of all the reactions (see “Output Specifications” -> “Plots” to activate or deactivate plots). Also plotting the species result after 0.5 hours will reproduce the literature results.
Format
SBML
(L3V1)
Related Publication
-
Physiologically based metformin pharmacokinetics model of mice and scale-up to humans for the estimation of concentrations in various tissues
- Darta Maija Zake, Janis Kurlovics, Linda Zaharenko, Vitalijs Komasilovs, Janis Klovins, Egils Stalidzans
- PLOS ONE , 4/ 2021 , Volume 16 , Issue 4 , pages: e0249594 , DOI: 10.1371/journal.pone.0249594
- Latvian Biomedical Research and Study Centre: Riga, LV; University of Latvia: Riga, LV
- Metformin is the primary drug for type 2 diabetes treatment and a promising candidate for other disease treatment. It has significant deviations between individuals in therapy efficiency and pharmacokinetics, leading to the administration of an unnecessary overdose or an insufficient dose. There is a lack of data regarding the concentration-time profiles in various human tissues that limits the understanding of pharmacokinetics and hinders the development of precision therapies for individual patients. The physiologically based pharmacokinetic (PBPK) model developed in this study is based on humans’ known physiological parameters (blood flow, tissue volume, and others). The missing tissue-specific pharmacokinetics parameters are estimated by developing a PBPK model of metformin in mice where the concentration time series in various tissues have been measured. Some parameters are adapted from human intestine cell culture experiments. The resulting PBPK model for metformin in humans includes 21 tissues and body fluids compartments and can simulate metformin concentration in the stomach, small intestine, liver, kidney, heart, skeletal muscle adipose, and brain depending on the body weight, dose, and administration regimen. Simulations for humans with a bodyweight of 70kg have been analyzed for doses in the range of 500-1500mg. Most tissues have a half-life (T1/2) similar to plasma (3.7h) except for the liver and intestine with shorter T1/2 and muscle, kidney, and red blood cells that have longer T1/2. The highest maximal concentrations (Cmax) turned out to be in the intestine (absorption process) and kidney (excretion process), followed by the liver. The developed metformin PBPK model for mice does not have a compartment for red blood cells and consists of 20 compartments. The developed human model can be personalized by adapting measurable values (tissue volumes, blood flow) and measuring metformin concentration time-course in blood and urine after a single dose of metformin. The personalized model can be used as a decision support tool for precision therapy development for individuals.
Contributors
Submitter of the first revision: Egils Stalidzans
Submitter of this revision: Krishna Kumar Tiwari
Modellers: Egils Stalidzans, Krishna Kumar Tiwari
Submitter of this revision: Krishna Kumar Tiwari
Modellers: Egils Stalidzans, Krishna Kumar Tiwari
Metadata information
is
isDescribedBy
hasTaxon
hasProperty
Mathematical Modelling Ontology
Ordinary differential equation model
ChEBI metformin
NCIt Intravenous Route of Administration
ChEBI metformin
NCIt Intravenous Route of Administration
unknownQualifier
Mathematical Modelling Ontology
Ordinary differential equation model
Mathematical Modelling Ontology Physiologically based pharmacokinetic model
BioModels Database MODEL2103020002
Mathematical Modelling Ontology Physiologically based pharmacokinetic model
BioModels Database MODEL2103020002
occursIn
Curation status
Curated
Modelling approach(es)
Tags
Connected external resources
Name | Description | Size | Actions |
---|---|---|---|
Model files |
|||
Zake2021_Metformin+Mice+IV.xml | SBML L3V1 file of the curated model | 221.33 KB | Preview | Download |
Additional files |
|||
Zake2021 - PBPK Metformin Mice IV single dose.cps | COPASI version: pharmacokinetics of metformin in mice, intra-venous, single dose | 281.39 KB | Preview | Download |
Zake2021 - PBPK Metformin Mice IV single dose.xml | SBML L2V4 pharmacokinetics of metformin in mice, intra-venous, single dose | 162.32 KB | Preview | Download |
Zake2021_Metformin+Mice+IV.cps | COPASI 4.34(Build251) file of the curated model | 300.39 KB | Preview | Download |
Zake2021_Metformin+Mice+IV.sedml | SEDML file of the curated model | 109.41 KB | Preview | Download |
- Model originally submitted by : Egils Stalidzans
- Submitted: Aug 27, 2021 4:39:14 PM
- Last Modified: Aug 27, 2021 4:39:14 PM
Revisions
Legends
: Variable used inside SBML models
: Variable used inside SBML models
Species
Species | Initial Concentration/Amount |
---|---|
mIntestineLumen UBERON:0018543 |
0.0 nmol |
mIntestineVascular C13311 ; intestine |
0.0 nmol |
mPlasmaArterial UBERON:0004572 |
0.0 nmol |
mPortalVein UBERON:0002017 |
0.0 nmol |
mStomach UBERON:0000945 |
0.0 nmol |
mIntestineEnterocytes BTO:0000398 ; intestine |
0.0 nmol |
mAdipose UBERON:0001013 |
0.0 nmol |
Reactions
Reactions | Rate | Parameters |
---|---|---|
mIntestineLumen => mFeces | k1*mIntestineLumen/IntestineLumen | k1=0.177684 |
mIntestineVascular => mPortalVein | QIntestineVascular*mIntestineVascular/IntestineVascular/Ktp_IntestineVascular | Ktp_IntestineVascular = 4.6; QIntestineVascular = 113.238 ml/h |
mPlasmaArterial => mKidneyPlasma | QKidney*mPlasmaArterial/PlasmaArterial | QKidney = 76.3308 ml/h |
mStomach => mPortalVein | QStomach*mStomach/Stomach/Ktp_Stomach | QStomach = 16.776 ml/h; Ktp_Stomach = 3.20000000001417 |
mIntestineLumen => mIntestineEnterocytes | k1*mIntestineLumen/IntestineLumen-k2*mIntestineEnterocytes/IntestineEnterocytes | k1=0.14; k2=0.14 |
mIntestineLumen => mIntestineVascular | k1*mIntestineLumen/IntestineLumen-k2*mIntestineVascular/IntestineVascular | k1=0.26; k2=0.26 |
mPlasmaArterial => mAdipose | QAdipose*mPlasmaArterial/PlasmaArterial | QAdipose = 4.194 ml/h |
mStomachLumen => mIntestineLumen | k1*mStomachLumen/StomachLumen | k1=3.16214 |