PK/PD model library

The PK/PD model library is combining our library of standard pharmacokinetic models described on the PK model library page, with a library of standard pharmacodynamic models.

Most PD models included in the PK/PD library can also be used alone (without a PK model) in the PD library, and some models from the PD library are only available alone.

Here we provide general guidelines to guide you towards the standard PD model that is the most suited to your dataset. The guidelines are also summarized in this guidelines pdf: summary_pkpd_library_v2.pdf
The complete list and full equations for PD models in the PD library is available in this document: PKPDlibrary.pdf Here we focus on the PD models of the PK/PD library because they are the most widely used.

PD models can be categorized by:

the type of response linking the concentration of the drug in the central compartment CC
to the effect E or to the response R. The effect can be either direct, or with an effect compartment, or an action on a turnover rate;
the type of drug action on the response or on the rate. It can be a stimulation or an inhibition;
the presence or absence of sigmoidicity in this drug action.

Type of response

The response links the drug concentration C_C to the effect E or the response R via a function A modeling the action of the drug. The type of response specifies whether this action impacts directly the effect, or via an effect compartment or turnover rates. The content of the action function will be specified in the next section: type of drug action. To explore the differences in the types of response, we will see how they impact the PD output for a 1 compartment PK linear infusion model, with a stimulation action of the drug and without sigmoidicity.

Direct

In case of a direct response, the effect E directly relates to the concentration in the central compartment C_Cvia

Therefore the dynamics of the effect follow the dynamics of C_C without any delay, meaning that as soon as C_C starts decreasing, E decreases as well, and the trajectory in the phase plane C_C vs E is a single line because it follows the same route during the increase and the decrease of the signals. Below is an example simulation of a direct model with a stimulation action of the drug on the PD.

Effect compartment

The site of action of the drug is often not be the same as the central compartment, and therefore a delay often appears between the PK and the PD measurements. To account for this delay it may be enough to add a single theoretical effect compartment with a drug concentration Ce equal to the concentration at the site of action. This theoretical compartment is assumed not to impact the concentration in the central compartment CC, therefore we write the transfer arrows between the central compartment and the effect compartment with dashed lines. A single transfer rate ke0 is used and the effect now directly relates to Ce:

Below is an example simulation of a model with an effect compartment and with a stimulation action of the drug on the PD. The higher the transfer parameter ke0, the higher the response and shorter the delay in the response.

Turnover

To model a PD response showing delay compared to the PK, it is also possible to consider that the PK concentration has an impact on the production or the degradation rate of a response variable R, as described by Sharma et al (1998). The response is then modeled with the following ODE:

if action on the production rate:
if action on the degradation rate:

In both cases, the model can be parametrized with the initial value of R R0 and the parameter kout, and kin is derived as .

We compare below the response in the case of production stimulation and degradation inhibition. kout plays a similar role as ke0 in the case of an effect compartment: the higher kout, the shorter the delay.

The difference between a turnover effect on the production or on the degradation is not obvious when checking the response to a single dose amount. A clear difference can be seen if two different dose amounts are given to the same individual. In the case of a turnover on production, the responses will show exactly the same delay, whereas in the case of a turnover on degradation, the delay in the response will be different for different dose amounts.

Type of drug action

The type of drug action will determine the function A relating the PK concentration to the PD effect. It can be a stimulation or an inhibition, and it is written differently depending on the type of response.

Stimulation

In the case of a direct model or an effect compartment, the action of stimulation is parametrized with Emax, EC50, γ if sigmoidicity is selected (otherwise γ=1), and E₀ if baseline is selected (otherwise E₀=0). The equation for A in the most general case is:

where C is the concentration in the central compartment C_C if the model is direct, and the concentration in the effect compartment C_Eif the model includes an effect compartment.

In the case of a turnover model, the action of stimulation can happen on the production or on the degradation rate. In any case, it is parametrized with Emax, EC50, and γ if sigmoidicity is selected, and the function always takes as an input the concentration in the central compartment CC. E₀ does not appear because turnover models already include a parameter R0 to model the initial value of the response R.

Below is the example response of a direct model with a stimulation action (called Emax in the drug action section of the library) and without sigmoidicity (γ=1). We show the influence of increasing E0, Emax, or the EC50. The influence of the sigmoidicity parameter γ is shown in the next section.

Inhibition

In the case of a direct model or an effect compartment, the action of inhibition is parametrized with E0, IC50, Imax if partial inhibition is selected (otherwise Imax=1), γ if sigmoidicity is selected (otherwise γ=1). The equation for A in the most general case is:

where C is the concentration in the central compartment C_C if the model is direct, and the concentration in the effect compartment C_Eif the model includes an effect compartment.

In the case of a turnover model, the action of inhibition can happen on the production or on the degradation rate. In any case, it is parametrized with IC50, Imax if partial inhibition is selected (otherwise Imax=1), and γ if sigmoidicity is selected, and the function always takes as an input the concentration in the central compartment CC. E0 does not appear because turnover models already include a parameter R0 to model the initial value of the response R.

Below is the example response of a direct model with an inhibition type of action (called Imax in the drug action section of the library) and without sigmoidicity (γ=1). We show the influence of increasing Imax, or the IC50. The influence of the sigmoidicity parameter γ is shown in the next section.

Sigmoidicity

Whether it is a stimulation or an inhibition type of action, the drug action function can always involve an exponent parameter γ. If this parameter is higher than 1, the response curve in the phase plane (C_C vs E) looks more like a sigmoid. Sigmoidicity can impact the dynamics in the PD response, making the response more switch-like due to threshold effects. If sigmoidicity is selected, the parameter γ will be estimated.

Below is the example response of a direct model with a stimulation type of action, and with sigmoidicity. The higher the gamma, the steeper the dose response, and the more square-like the time response.