Target-mediated Drug Disposition Model Research Articles

Application of target-mediated drug disposition (TMDD) models in the development of monoclonal antibody

Application of target-mediated drug disposition (TMDD) models in the development of monoclonal antibody

Translational Pharmacokinetic/Pharmacodynamic Characterization and Target-Mediated Drug Disposition Modeling of an Anti–Tissue Factor Pathway Inhibitor Antibody, PF-06741086

Tissue factor pathway inhibitor (TFPI) exhibits multiple isoforms, which are known to present in multiple locations such as plasma, endothelium, and platelets. TFPI is an endogenous negative modulator of the coagulation pathway, and therefore, neutralization of TFPI function can potentially increase coagulation activity. A human monoclonal antibody, PF-06741086, which interacts with all isoforms of TFPI is currently being tested in clinic for treating hemophilia patients with and without inhibitors. To support clinical development of PF-06741086, pharmacokinetics (PK) and pharmacodynamics of PF-06741086 were characterized in monkeys. In addition, a mechanistic model approach was used to estimate PK parameters in monkeys and simulate PK profiles in human. The results show that PF-06741086 exhibited target-mediated drug disposition and had specific effects on various hemostatic markers including diluted prothrombin time, thrombin generation, and thrombin-antithrombin complex in monkeys after administration. The model-predicted and observed human exposures were compared retrospectively, and the result indicates that the exposure prediction was reasonable within less than 2-fold deviation. This study demonstrated in vivo efficacy of PF-06741086 in monkeys and the utility of a rational mechanistic approach to describe PK for a monoclonal antibody with complex target binding.

Mathematical description of drug-target interactions: application to biologics that bind to targets with two binding sites.

The emerging discipline of mathematical pharmacology occupies the space between advanced pharmacometrics and systems biology. A characteristic feature of the approach is application of advance mathematical methods to study the behavior of biological systems as described by mathematical (most often differential) equations. One of the early application of mathematical pharmacology (that was not called this name at the time) was formulation and investigation of the target-mediated drug disposition (TMDD) model and its approximations. The model was shown to be remarkably successful, not only in describing the observed data for drug-target interactions, but also in advancing the qualitative and quantitative understanding of those interactions and their role in pharmacokinetic and pharmacodynamic properties of biologics. The TMDD model in its original formulation describes the interaction of the drug that has one binding site with the target that also has only one binding site. Following the framework developed earlier for drugs with one-to-one binding, this work aims to describe a rigorous approach for working with similar systems and to apply it to drugs that bind to targets with two binding sites. The quasi-steady-state, quasi-equilibrium, irreversible binding, and Michaelis-Menten approximations of the model are also derived. These equations can be used, in particular, to predict concentrations of the partially bound target (RC). This could be clinically important if RC remains active and has slow internalization rate. In this case, introduction of the drug aimed to suppress target activity may lead to the opposite effect due to RC accumulation.

Target-mediated drug disposition model for drugs with two binding sites that bind to a target with one binding site.

The paper extended the TMDD model to drugs with two identical binding sites (2-1 TMDD). The quasi-steady-state (2-1 QSS), quasi-equilibrium (2-1 QE), irreversible binding (2-1 IB), and Michaelis-Menten (2-1 MM) approximations of the model were derived. Using simulations, the 2-1 QSS approximation was compared with the full 2-1 TMDD model. As expected and similarly to the standard TMDD for monoclonal antibodies (mAb), 2-1 QSS predictions were nearly identical to 2-1 TMDD predictions, except for times of fast changes following initiation of dosing, when equilibrium has not yet been reached. To illustrate properties of new equations and approximations, several variations of population PK data for mAbs with soluble (slow elimination of the complex) or membrane-bound (fast elimination of the complex) targets were simulated from a full 2-1 TMDD model and fitted to 2-1 TMDD models, to its approximations, and to the standard (1-1) QSS model. For a mAb with a soluble target, it was demonstrated that the 2-1 QSS model provided nearly identical description of the observed (simulated) free drug and total target concentrations, although there was some minor bias in predictions of unobserved free target concentrations. The standard QSS approximation also provided a good description of the observed data, but was not able to distinguish between free drug concentrations (with no target attached and both binding site free) and partially bound drug concentrations (with one of the binding sites occupied by the target). For a mAb with a membrane-bound target, the 2-1 MM approximation adequately described the data. The 2-1 QSS approximation converged 10 times faster than the full 2-1 TMDD, and its run time was comparable with the standard QSS model.

Pharmacokinetic-Pharmacodynamic Relationship of Erenumab (AMG 334) and Capsaicin-Induced Dermal Blood Flow in Healthy and Migraine Subjects

PurposeCapsaicin-induced dermal blood flow (CIDBF) is a validated biomarker used to evaluate the target engagement of potential calcitonin gene-related peptide-blocking therapeutics for migraine. To characterize the pharmacokinetics (PK) and quantify the inhibitory effects of erenumab (AMG 334) on CIDBF, CIDBF data were pooled from a single- and a multiple-dose study in healthy and migraine subjects.MethodsRepeated capsaicin challenges and DBF measurements were performed and serum erenumab concentrations determined. A population analysis was conducted using a nonlinear mixed-effects modeling approach. Effects of body weight, gender, and age on model parameters were evaluated.ResultsTwo-compartment target-mediated drug disposition (TMDD) model assuming binding of erenumab in the central compartment best described the nonlinear PK of erenumab. Subcutaneous absorption half-life was 1.6 days and bioavailability was 74%. Erenumab produced a maximum inhibition of 89% (95% confidence interval: 87–91%). Erenumab concentrations required for 50% and 99% of maximum inhibition were 255 ng/mL and 1134 ng/mL, respectively. Increased body weight was associated with increased erenumab clearance but had no effect on the inhibitory effect on CIDBF.ConclusionsOur results show that erenumab pharmacokinetics was best characterized by a TMDD model and resulted in potent inhibition of CIDBF.

Mechanistic Modeling of the Pharmacodynamic and Pharmacokinetic Relationship of Tissue Factor Pathway Inhibitor-Neutralizing Antibody (BAY 1093884) in Cynomolgus Monkeys

BAY 1093884 is a fully human monoclonal antibody against the tissue factor pathway inhibitor (TFPI) in development as prophylaxis in patients with hemophilia with or without inhibitors. In vitro, BAY 1093884 binds to human, mouse, and monkey TFPI. The objective of this study was to find a pharmacodynamic (PD) biomarker after administration of BAY 1093884 to normal monkeys. In monkey plasma, BAY 1093884 exhibited an IC50 (concentration that inhibits 50%) of 4.65 and 6.19 nM for free TFPI and diluted prothrombin time (dPT), respectively. The BAY 1093884 pharmacokinetic (PK) profile and its PD effects on dPT and free TFPI levels were assessed after intravenous and subcutaneous administration of BAY 1093884 (5 and 20 mg/kg) to female cynomolgus monkeys. Free TFPI concentrations in plasma decreased rapidly and increased to baseline in a dose-dependent manner. dPT clotting time was shortened and correlated with free TFPI levels and drug concentration in plasma, demonstrating the relationship between PD activities (dPT clotting time and free TFPI levels) and drug concentration. BAY 1093884 exhibited nonlinear PK, and a target-mediated drug disposition model was used to characterize the BAY 1093884 versus TFPI concentration–response relationship. We concluded that a mechanism-based PK/PD binding model could be useful for predicting human response to BAY 1093884. For the first-in-human study, measurement of free TFPI will be included as part of the dose-escalation design.

Target-Mediated Drug Disposition Pharmacokinetic–Pharmacodynamic Model of Bosentan and Endothelin-1

Bosentan is a competitive antagonist on endothelin receptor A and B (ETA and ETB), displacing the endogenous binding partner endothelin-1 (ET-1) from its binding sites. After administration of escalating single doses of 10-750mg as an intravenous (i.v.) infusion, bosentan showed dose-dependent pharmacokinetics (PK). The aim of this analysis was to develop a PK model of bosentan after i.v. administration including competitive antagonism with ET-1 and to analyze its influence on blood pressure and heart rate with a combined pharmacokinetic/pharmacodynamic (PK/PD) model. PK/PD data from 70 young male Caucasian subjects were analyzed after single i.v. administration of 10, 50, 250, 500, and 750mg of bosentan. Population analyses, simulations, and evaluation were performed using a non-linear mixed-effects modeling approach. The PK of bosentan was best described by a two-compartment, target-mediated drug disposition (TMDD) model. ET-1 plasma and urine profiles were successfully integrated into the bosentan two-compartment, TMDD model encompassing competition for the same receptor. A multiple-peak phenomenon of bosentan plasma concentrations after i.v. administration was best described by a diurnal expression or reappearance of ET receptors on the cell surface. Blood pressure was best described by an E max model; heart rate was modeled as a compensatory effect of changes in blood pressure. The developed competitive PK/PD model of bosentan and ET-1 after i.v. administration provides a first step towards understanding the complex PK properties of bosentan and offers a valuable tool for future PK/PD research.

Clinical Pharmacology and Translational Aspects of Bispecific Antibodies.

Clinical Pharmacology and Translational Aspects of Bispecific Antibodies.

Target-mediated drug disposition with drug-drug interaction, Part I: single drug case in alternative formulations.

Target-mediated drug disposition (TMDD) describes drug binding with high affinity to a target such as a receptor. In application TMDD models are often over-parameterized and quasi-equilibrium (QE) or quasi-steady state (QSS) approximations are essential to reduce the number of parameters. However, implementation of such approximations becomes difficult for TMDD models with drug-drug interaction (DDI) mechanisms. Hence, alternative but equivalent formulations are necessary for QE or QSS approximations. To introduce and develop such formulations, the single drug case is reanalyzed. This work opens the route for straightforward implementation of QE or QSS approximations of DDI TMDD models. The manuscript is the first part to introduce DDI TMDD models with QE or QSS approximations.

Target mediated drug disposition with drug-drug interaction, Part II: competitive and uncompetitive cases.

We present competitive and uncompetitive drug-drug interaction (DDI) with target mediated drug disposition (TMDD) equations and investigate their pharmacokinetic DDI properties. For application of TMDD models, quasi-equilibrium (QE) or quasi-steady state (QSS) approximations are necessary to reduce the number of parameters. To realize those approximations of DDI TMDD models, we derive an ordinary differential equation (ODE) representation formulated in free concentration and free receptor variables. This ODE formulation can be straightforward implemented in typical PKPD software without solving any non-linear equation system arising from the QE or QSS approximation of the rapid binding assumptions. This manuscript is the second in a series to introduce and investigate DDI TMDD models and to apply the QE or QSS approximation.

Application of target-mediated drug disposition model to predict human pharmacokinetics and target occupancy of GC1118, an anti-epidermal growth factor receptor antibody

Application of target-mediated drug disposition model to predict human pharmacokinetics and target occupancy of GC1118, an anti-epidermal growth factor receptor antibody

A mathematical analysis of rebound in a target-mediated drug disposition model: II. With feedback

We consider the possibility of free receptor (antigen/cytokine) levels rebounding to higher than the baseline level after the application of an antibody drug using a target-mediated drug disposition model. It is assumed that the receptor synthesis rate experiences homeostatic feedback from the receptor levels. It is shown for a very fast feedback response, that the occurrence of rebound is determined by the ratio of the elimination rates, in a very similar way as for no feedback. However, for a slow feedback response, there will always be rebound. This result is illustrated with an example involving the drug efalizumab for patients with psoriasis. It is shown that slow feedback can be a plausible explanation for the observed rebound in this example.

Use of a Target-Mediated Drug Disposition Model to Predict the Human Pharmacokinetics and Target Occupancy of GC1118, an Anti-epidermal Growth Factor Receptor Antibody.

GC1118 is an anti-epidermal growth factor receptor (EGFR) monoclonal antibody that is currently under clinical development. In this study, the pharmacokinetics (PK) of GC1118 were modelled in monkeys to predict human PK and receptor occupancy (RO) profiles. The serum concentrations of GC1118 and its comparator (cetuximab) were assessed in monkeys with a non-compartmental analysis and a target-mediated drug disposition (TMDD) model after intravenous infusion (3-25 mg/kg) of these drugs. The scaling exponent of the EGFR synthesis rate was determined using a sensitivity analysis. The human cetuximab exposures were simulated by applying different exponents (0.7-1.0) for the EGFR synthesis rate in the allometric monkey PK model. Simulated Cmax and area under the curve values therein were compared with those previously reported in the literature to find the best exponent for the EGFR synthesis rate in human beings. The TMDD model appropriately described the monkey PK profile, which showed a decrease in clearance (CL; 1.2-0.4 ml/hr/kg) as the dose increased. The exponents for CL (0.75) and volume of distribution (Vd; 1.0) were used for the allometric scaling to predict human PK. The allometric coefficient for the EGFR synthesis rate chosen by the sensitivity analysis was 0.85, and the RO profiles that could not be measured experimentally were estimated based on the predicted concentrations of the total target and the drug-target complex. Our monkey TMDD model successfully predicts human PK and RO profiles of GC1118 and can be used to determine the appropriate dose for a first-in-human study investigating this drug.

A cell-level model of pharmacodynamics-mediated drug disposition.

We aimed to develop a cell-level pharmacodynamics-mediated drug disposition (PDMDD) model to analyze in vivo systems where the PD response to a drug has an appreciable effect on the pharmacokinetics (PK). An existing cellular level model of PD stimulation was combined with the standard target-mediated drug disposition (TMDD) model and the resulting model structure was parametrically identifiable from typical in vivo PK and PD data. The PD model of the cell population was controlled by the production rate k in and elimination rate k out which could be stimulated or inhibited by the number of bound receptors on a single cell. Simulations were performed to assess the impact of single and repeated dosing on the total drug clearance. The clinical utility of the cell-level PDMDD model was demonstrated by fitting published data on the stimulatory effects of filgrastim on absolute neutrophil counts in healthy subjects. We postulated repeated dosing as a means of detecting and quantifying PDMDD as a single dose might not be sufficient to elicit the cellular response capable of altering the receptor pool to visibly affect drug disposition. In the absence of any PD effect, the model reduces down to the standard TMDD model. The applications of this model can be readily extended to include chemotherapy-induced cytopenias affecting clearance of endogenous hematopoietic growth factors, different monoclonal antibodies and immunogenicity effects on PK.

Asymptotic Analysis of a Target-Mediated Drug Disposition Model: Algorithmic and Traditional Approaches.

A detailed analysis is reported on a multiscale pharmacokinetic model, simulating the interaction of a drug with its target, the binding of the compounds and the outcome of their interaction. The analysis is based on the algorithmic computational singular perturbation (CSP) methodology. Among others, the analysis concludes that the partial equilibrium approximation and the quasi-steady-state approximation (PEA and QSSA) are valid in two distinct stages in the evolution of the process. Similar conclusions are reached from the algorithmic criteria for the validity of the QSSA and PEA. The reactions in the pharmacokinetic model that (i) generate the fast time scales, (ii) generate the constraints in which the system evolves and (iii) drive the system at various phases are identified, with the use of algorithmic CSP tools. These identifications are very important for the improvement of the model and for the identification of ways to control the evolution of the process. Regarding the qualitative understanding of the process, the present analysis systematises the findings in the literature and provides some new insights.

First-in-human application of the novel hepatitis B and hepatitis D virus entry inhibitor myrcludex B

Myrcludex B is a first-in-class compound, which blocks entry of hepatitis B and D virus into hepatocytes in vitro and in animal models. Based on the required preclinical data we aimed to translate this compound into the first application in humans. Single ascending doses of myrcludex B, a 47 amino acid peptide, were administered up to 20mg intravenously and 10mg subcutaneously in a prospective open first-in-human, phase I clinical trial to 36 healthy volunteers. Safety, tolerability and plasma concentrations of myrcludex B were assessed and a pharmacokinetic model was derived. Myrcludex B was well tolerated and no serious or relevant AEs representing off-target effects, and no immunogenic effects were observed up to the highest applied dose of 20mg (intravenously). Myrcludex B showed dose-dependent pharmacokinetics, best described by a 2-compartment target-mediated drug disposition model. Bioavailability of the subcutaneous application was large (85%). Interindividual variability was moderate. The pharmacokinetic model suggested that subcutaneous doses of 10mg and above reach a target saturation of over 80% for at least 15h. Myrcludex B showed excellent tolerability up to high doses. Pharmacologic properties followed a 2-compartment target-mediated drug disposition model. These findings are vital for planning of further multiple dose efficacy trials in patients. After showing antiviral activity in cell culture and animal models, myrcludex B, a new drug intended for the treatment of hepatitis B and D, has been administered the first time in humans. Healthy volunteers received the drug intravenously and subcutaneously up to high doses (20mg). The drug was well tolerated and the characteristics of the drug determining its way in the human body could be described. These results will allow testing myrcludex B in hepatitis B and D patients.

Population Pharmacokinetics and Pharmacodynamics of Lampalizumab Administered Intravitreally to Patients With Geographic Atrophy.

Intravitreally administered lampalizumab is an investigational complement inhibitor directed against complement factor D (CFD) for the treatment of geographic atrophy (GA) secondary to age-related macular degeneration. We sought to develop an integrated ocular and systemic pharmacokinetic/pharmacodynamic model for lampalizumab in patients with GA using the data from the clinical phase I and II studies. The kinetics of lampalizumab and CFD disposition were well described by the combined ocular/serum target-mediated drug disposition model using a quasi-steady-state approximation. This model takes into account the drug, target, and drug–target complex clearance, their transfer rates between ocular and serum compartments, and turnover kinetics of CFD. The constructed model provided a prediction of target occupancy in ocular tissues and supported that the two dosing regimens (10 mg q4w and 10 mg q6w) selected for the phase III studies are expected to be efficacious and able to achieve near-complete target engagement in the vitreous humor.

A Priori Identifiability of Target-Mediated Drug Disposition Models and Approximations.

A priori identifiability of mathematical models assures that for a given input/output experiment, the parameter set has one unique solution within a defined space, independent of the experimental design. Many biologic therapeutics exhibit target-mediated drug disposition (TMDD), and use of the full compartmental model describing this system is well documented. In practice, estimation of the full parameter set for TMDD models, given real-world clinical data, is characterized by convergence difficulties and unstable solutions. Still, the formal assessment of the a priori identifiability of these systems has yet to be reported. The exact arithmetic rank (EAR) approach was used to test the a priori identifiability of a TMDD model as well as model approximations. The full TMDD and quasi-equilibrium/rapid binding (QE/RB), quasi-steady state (QSS), and Michaelis-Menten (MM) approximations were fully identifiable, a priori, regardless of whether observations were taken from a single or multiple compartments. The results of these identifiability analyses indicated that the difficulty with TMDD model convergence, a posteriori, lies in the experimental design, not in the mathematical identifiability in the lack of samples from several compartments. Experiments can be tailored to resolve these structurally non-identifiable parameters, notwithstanding practical implementation challenges. This work highlights the importance of identifiability analyses, specifically how they can influence experimental design and selection of the appropriate model structure to describe a dynamic biological system.

A Tutorial on Target-Mediated Drug Disposition (TMDD) Models.

Target-mediated drug disposition (TMDD) is the phenomenon in which a drug binds with high affinity to its pharmacological target site (such as a receptor) to such an extent that this affects its pharmacokinetic characteristics.1 The aim of this Tutorial is to provide an introductory guide to the mathematical aspects of TMDD models for pharmaceutical researchers. Examples of Berkeley Madonna2 code for some models discussed in this Tutorial are provided in the Supplementary Materials.

A pharmacokinetic binding model for bevacizumab and VEGF165 in colorectal cancer patients.

Purpose To characterize the population pharmacokinetics of bevacizumab, its binding properties to VEGF165 and the effect of demographic data and VEGF-A polymorphisms on the interplay between bevacizumab serum pharmacokinetics and VEGF165 serum concentrations in patients with colorectal cancer stage IV.MethodsBevacizumab and VEGF165 data were collected from 19 adult patients with metastatic colorectal cancer enrolled in an observational clinical study. Bevacizumab was administered with one of the following combinations: 5-FU/Leucovorin/Irinotecan, 5-FU/Leucovorin/Oxaliplatin, Capecitabine/Irinotecan at doses ranging from 5 to 10 mg/kg every 2 or 3 weeks. Data analysis was performed using nonlinear mixed-effects modeling implemented in NONMEM 7.3.ResultsA target-mediated drug disposition model adequately described bevacizumab concentration changes over time and its binding characteristics to VEGF165. The estimated clearance of bevacizumab was 0.18 L/day, the free VEGF165 levels at baseline were 212 ng/L, and the elimination rate constant of free VEGF165 was 0.401 day−1. Body weight was allometrically included in all PK parameters.ConclusionThe final model adequately described the pre- and post-dose concentrations of total bevacizumab and free VEGF165 in patients with colorectal cancer. Model parameters were consistent with those previously reported for patients with solid tumors. Correlations between the binding affinity of bevacizumab and the VEGF-2578C/A and VEGF-634G/C polymorphisms were noticed.Electronic supplementary materialThe online version of this article (doi:10.1007/s00280-015-2701-3) contains supplementary material, which is available to authorized users.