General relationship between transit compartments and lifespan models

Delay differential equations based models in NONMEM

Delay differential equations (DDEs) are commonly used in pharmacometric models to describe delays present in pharmacokinetic and pharmacodynamic data analysis. Several DDE solvers have been implemented in NONMEM 7.5 for the first time. Two of them are based on algorithms already applied elsewhere, while others are extensions of existing ordinary differential equations (ODEs) solvers. The purpose of this tutorial is to introduce basic concepts underlying DDE based models and to show how they can be developed using NONMEM. The examples include previously published DDE models such as logistic growth, tumor growth inhibition, indirect response with precursor pool, rheumatoid arthritis, and erythropoiesis-stimulating agents. We evaluated the accuracy of NONMEM DDE solvers, their ability to handle stiff problems, and their performance in parameter estimation using both first-order conditional estimation (FOCE) and the expectation-maximization (EM) method. NONMEM control streams and excerpts from datasets are provided for all discussed examples. All DDE solvers provide accurate and precise solutions with the number of significant digits controlled by the error tolerance parameters. For estimation of population parameters, the EM method is more stable than FOCE regardless of the DDE solver.

Exact solutions and equi-dosing regimen regions for multi-dose pharmacokinetics models with transit compartments

Compartmental models which yield linear ordinary differential equations (ODEs) provide common tools for pharmacokinetics (PK) analysis, with exact solutions for drug levels or concentrations readily obtainable for low-dimensional compartment models. Exact solutions enable valuable insights and further analysis of these systems. Transit compartment models are a popular semi-mechanistic approach for generalising simple PK models to allow for delayed kinetics, but computing exact solutions for multi-dosing inputs to transit compartment systems leading to different final compartments is nontrivial. Here, we find exact solutions for drug levels as functions of time throughout a linear transit compartment cascade followed by an absorption compartment and a central blood compartment, for the general case of n transit compartments and M equi-bolus doses to the first compartment. We further show the utility of exact solutions to PK ODE models in finding constraints on equi-dosing regimen parameters imposed by a prescribed therapeutic range. This leads to the construction of equi-dosing regimen regions (EDRRs), providing new, novel visualisations which summarise the safe and effective dosing parameter space. EDRRs are computed for classical and transit compartment models with two- and three-dimensional parameter spaces, and are proposed as useful graphical tools for informing drug dosing regimen design.

Distributed transit compartments for arbitrary lifespan distributions in aging populations

Transit compartment models (TCM) are often used to describe aging populations where every individual has its own lifespan. However, in the TCM approach these lifespans are gamma-distributed which is a serious limitation because often the Weibull or more complex distributions are realistic. Therefore, we extend the TCM concept to approximately describe any lifespan distribution and call this generalized concept distributed transit compartment models (DTCMs). The validity of DTCMs is obtained by convergence investigations. From the mechanistic perspective the transit rates are directly controlled by the lifespan distribution. Further, DTCMs could be used to approximate the convolution of a signal with a probability density function. As example a stimulatory effect of a drug in an aging population with a Weibull-distributed lifespan is presented where distribution and model parameters are estimated based on simulated data.

Lifespan based pharmacokinetic-pharmacodynamic model of tumor growth inhibition by anticancer therapeutics

Accurate prediction of tumor growth is critical in modeling the effects of anti-tumor agents. Popular models of tumor growth inhibition (TGI) generally offer empirical description of tumor growth. We propose a lifespan-based tumor growth inhibition (LS TGI) model that describes tumor growth in a xenograft mouse model, on the basis of cellular lifespan T. At the end of the lifespan, cells divide, and to account for tumor burden on growth, we introduce a cell division efficiency function that is negatively affected by tumor size. The LS TGI model capability to describe dynamic growth characteristics is similar to many empirical TGI models. Our model describes anti-cancer drug effect as a dose-dependent shift of proliferating tumor cells into a non-proliferating population that die after an altered lifespan T_A. Sensitivity analysis indicated that all model parameters are identifiable. The model was validated through case studies of xenograft mouse tumor growth. Data from paclitaxel mediated tumor inhibition was well described by the LS TGI model, and model parameters were estimated with high precision. A study involving a protein casein kinase 2 inhibitor, AZ968, contained tumor growth data that only exhibited linear growth kinetics. The LS TGI model accurately described the linear growth data and estimated the potency of AZ968 that was very similar to the estimate from an established TGI model. In the case study of AZD1208, a pan-Pim inhibitor, the doubling time was not estimable from the control data. By fixing the parameter to the reported in vitro value of the tumor cell doubling time, the model was still able to fit the data well and estimated the remaining parameters with high precision. We have developed a mechanistic model that describes tumor growth based on cell division and has the flexibility to describe tumor data with diverse growth kinetics.

Solution and implementation of distributed lifespan models

We consider a population where every individual has a unique lifespan. After expiring of its lifespan the individual has to leave the population. A realistic approach to describe these lifespans is by a continuous distribution. Such distributed lifespan models (DLSMs) were introduced earlier in the indirect response context and consist of the mathematical convolution operator to describe the rate of change. Therefore, DLSMs could not directly be implemented in standard PKPD software. In this work we present the solution representation of DLSMs with and without a precursor population and an implementation strategy for DLSMs in ADAPT , NONMEM and MATLAB . We fit hemoglobin measurements from literature and investigate computational properties.