Pharmacolibrary - Modelica 2025

Created by Tomas Kulhanek

pharmacolibrary

#modelica, #pharmacolibrary

Pharmacolibrary - Free Library to Model Pharmacology

Tomáš Kulhánek¹²³, Filip Ježek²⁴, Jiří Kofránek²³,

Marek Mateják²⁵, Stef Romes¹

¹Flemish Institute of Technology and Research (VITO), Mol, BE

²1st Faculty of Medicine, Charles University, Prague, CZ

³Creative Connections s.r.o, Prague, CZ

⁴University of Michigan, Ann Arbor, USA

⁵Institute of Clinical end Experimental Medicine, Prague, CZ

VITO

Flemish Institute for Technological Research

Mol, Belgium

Environmental unit

Digital Precision Health group

RAIDO

PRIBE - personal reference intervals

brave ZR8tgl5Xml

DigiTwin4PH - PGx digital twin

ONCOSCREEN

Pharmacology

Pharmacokinetics (PK) body → drug

Pharmacodynamics (PD) drug → body

Pharmacogenomics (PGx) genetics→drug&body

pharmacokinetics(PK): body → drug

Administration(absorption), Distribution, Metabolism, Elimination(excretion)

connectors

Administration (absorption)

parenteral (e.g. intravenous, main parameters: adminMass,adminDuration):

enteral (e.g. oral, main parameters: adminMass, Tlag, F absorption fraciton, ka absorption rate)

Distribution

simplified - ideal mixing, no volume flow rate, compartment (Vd volume of distribution)

$q_m =\frac{dM}{dt}$

$C=\frac{M}{V_d}$

and 1st order kinetics

q_m=Cl_a*C_a-Cl_b*C_b

Distribution - physiology based

compartment - with concentration mixing based on volumetric flow through organs and tissues:

$\frac{dM}{dt} = q_{v_{a}}*actualStream(c_a)+q_{v_{b}}*actualStream(c_b)+q_m$

fixed volumetric flow component:

Metabolism

simplified - using 1st order rate

e.g. paracetamol

acetaminophen → sulfate_conjugate

for more complex - use Chemical library

Elimination

via kidney or liver using clearance rate Cl

q_m=Cl*C

1-compartment model

parameters for intravenous administration of gentamicin, type of antibiotic against bacterial infections

equivalent 1-compartment model (in equations)

VirtualBoxVM 4UK0GmE2mR

gentamicin - 1 day simulation

gentamicin - 1 day simulation logscale

pk 1c simlog

2-compartment model

parameters for intravenous administration of midazolam, a short-acting benzodiazepine used for sedation, anesthesia, ...

midazolam - 12 h simulation

logscale - 12h simulation

logscale - 12h with peripheral

pk2csimlog

3-compartment model

fentanyl - 10h simulation

fentanyl - 10h simulation logscale

pk3csim

physiology based

compartment takes physiological features into account

cardiac output
lymphatic system flow
detail organ tissue

PBPK model

PBPK simulation

5496 models of drugs by Anatomical Therapeutical Chemical (ATC) classification - knowledge base extracted and combined from Pubmed, Drugbank using LLM

ATC level 1

ATC level 3

ATC level 5

Drug model variants

caffeine

1-compartment model - caffeine

1 cup = 100 mg of caffeine in 8 hours

caffeine in 24 hours

pharmacodynamics (PD): drug → body

What effect has a drug to body

linear effect

effectlin

Emax effect

effectemax

sigmoid Emax effect

effectsigmoid

PK PD model Midazolam on EEG activity

midazolam

PK PD simulation

PK PD model digoxin on hemodynamics of cardiovascular system

PKPD HemodynamicsMeurs flatNorm

PK PD simulation

pkpd sim1

pkpd sim2

pkpd sim4

pkpd sim3

pressure volume (PV) diagram

PV shift from right to left → less maximal volume, slight increase in systolic pressure.

Digoxin is a cardiac glycoside derived from the foxglove plant Digitalis lanata. It is primarily used in the treatment of various heart conditions, notably atrial fibrillation, atrial flutter, and sometimes heart failure that cannot be controlled by other medications. Digoxin is approved and widely used in clinical practice today.

Pharmacogenomics (PGx)

Howto systematically approach PGx in Modelica and Pharmacolibrary?

genotype
phenotype

genotype

phenotype

SLC22A2

OTC2_phenotype

alleles

*1, *1

808T, *1

clearance

normal

intermediate

CYP2C19

CYP2C19_phenotype

alleles

*1, *1

*2, *1

*17, *1

metabolizer

normal

poor

rapid

PK PGx model

gentamicin

PK PGx simulation (gentamicin)

population PGx model

pgmodelexperiment

population PGx simulation

pkpgx experiment sim

optimal dose of lithium carbonate

brave OZ24t5RrJ2

https://egolem.online/pharma

iron metabolism visualisation

brave 5rh0hlzHNL

eGolem powered by bodylight.js

https://egolem.online/irm/ https://bodylight.physiome.cz

howto find a drug

VirtualBoxVM quF8UzRISk

workflow - 1. export model to FMI unit 2. simulate in Python & fmpy & 3. get data from external data (PK-DB.com) 4. compare/fit/...

VirtualBoxVM TIsoOjFLH2

Conclusion

Pharmacolibrary: https://github.com/creative-connections/Pharmacolibrary

pharmacokinetic: 5496 drugs models - most generated

pharmacodynamics: various effects

pharmacogenomics: various genotype/phenotype influence

current and future research and development: contribute to democratize pharmacology knowledge, use PBPK for quantitative in vitro to in vivo extrapolation (QIVIVE), validate PK models of drugs, add PK of dietary supplements and toxicokinetics of chemicals, cover published PD and common pharmacogenes,

Pharmacolibrary / GITHUB

https://github.com/creative-connections/Pharmacolibrary

acknowledgment

ONCOSCREEN project, funded by the European Union’s Horizon Europe grant No. 101097036.

subaward to R01 HL173346. funded by NIH

Readability & Plausability

small models - equation based

VirtualBoxVM 4UK0GmE2mR

standard Modelica library 4.0.0

respiration

Dymola sdXGmJf24o

medium/complex models- composition, inheritance, polymorphism (OOM) of smaller models in diagrams

Dymola 9CIA5zGnUy

standard Modelica language

Commercial tools:

Dymola,

, …

Open-source:

Physiomodel

firefox SHzAoizS4q

Complex model of physiology in Modelica and Physiolibrary

https://physiomodel.org

blood-gas exchange

Dymola qmDZxxtrfj

hemodynamics of cardiovascular system

Dymola 9CIA5zGnUy

comprehensibility

why Modelica?

co-simulation of multiple models via FMI in static web app with bodylight.js

https://egolem.online/demo/

Guyton's original scheme part of circulatory dynamics

equivalent implementation in component and acausal modeling language

custom libraries for specific domains

Physiolibrary 3.0.0

www.physiolibrary.org

Dymola qmDZxxtrfj

Frank-Starling law is much more visible in component model - "force of the heart muscle depends on pressure of incoming blood".

model of complex system may lead to implementation errors

A C Guyton, T G Coleman, H J Granger, Circulation: overall regulation, Annu Rev Physiol.1972;34:13-46.doi: 10.1146/annurev.ph.34.030172.000305.

Correction of schema after reimplementation in 2010

J.Kofránek, et al., Restoration of Guyton‘s Diagram for Regulation of the Circulation ...Physiol Res. 2010;59(6):897-908.doi: 10.33549/physiolres.931838.

Pharmacolibrary 25.06 (beta)

connectors and components:

effects

PGx

40-50 pharmacogenes commercially available in testing panels

majority cytochrome P450 enzymes ( CYP2D6,CYP2C19,...)