Holford NHG and Sheiner LB "Understanding the Dose-Effect Relationship-Clinical Application of Pharmacokinetic-Pharmacodynamic Models", Clin Pharmacokin 6:429-453 (1981)-The Backstory

Safety, tolerability, pharmacokinetics and pharmacodynamics of a novel farnesoid X receptor (FXR) agonist-TQA3526 in healthy Chinese volunteers: a double-blind, randomized, placebo-controlled, dose-escalation, food effect phase I study

Background: TQA3526 is a novel farnesoid X receptor agonist developed to treat non-alcoholic steatohepatitis (NASH) or primary biliary cholangitis (PBC). This study aimed to evaluate the safety, tolerability, pharmacokinetics (PK), and pharmacodynamics (PD) of TQA3526 in healthy Chinese patients.Methods: Healthy subjects aged 18-55 years were enrolled in this double-blinded, first-in-human, placebo-controlled single ascending dose (1, 2, 5, and 10 mg) comprising food effect investigation (10 mg) and multiple dose study (2 mg and 0.2 + 0.5 + 1 mg). Safety was assessed on the basis of adverse events. The TQA3526 concentrations were analysed in the PK study. Alkaline phosphatase (ALP), fibroblast growth factor-19 (FGF19), bile acid precursor C4 (7α-hydroxy-cholest-4-ene-3-one), cholesterol, and bile acid were selected for PD analysis.Results: TQA3526 was well tolerated, and the primary adverse drug reaction was pruritus, as expected. The exposure to TQA3526 increased in a dose-dependent manner after a single dose of 1-10 mg. The exposure was higher after food intake. A steady state was reached around 5 days, and obvious plasma accumulation of TQA3526 was observed in the multiple dose study. TQA3526 increased circulating FGF-19 and decreased C4 levels in a dose-dependent manner. ALP increased only mildly in the 2 mg multiple dose cohort.Conclusions: TQA3526 (<10 mg/day) was safe and tolerable in healthy Chinese subjects. The safety profile and PK/PD characteristics of TQA3526 support further evaluation of patients with NASH or PBC. This study was registered at https://www.chictr.org.cn/ under the identifier ChiCTR1800019570.

Riluzole in Spinal Cord Injury Study (RISCIS)-Pharmacokinetic (PK) Sub-Study: An Analysis of Pharmacokinetics, Pharmacodynamics, and Impact on Axonal Degradation of Riluzole in Patients With Traumatic Cervical Spinal Cord Injury Enrolled in the RISCIS Phase III Randomized Controlled Trial

To date, no drug therapy has shown significant efficacy in improving functional outcomes in patients with acute spinal cord injury (SCI). Riluzole is an approved benzothiazole sodium channel blocker to attenuate neurodegeneration in amyotrophic lateral sclerosis (ALS) and is of interest for neuroprotection in SCI. In a Phase I clinical trial (ClinicalTrials.gov Identifier: NCT00876889), riluzole was well tolerated with a 2-week treatment at the dose level approved for ALS and exhibited potential efficacy in patients with SCI. The acute and progressive nature of traumatic SCI and the complexity of secondary injury processes alter the pharmacokinetics (PK) of therapeutics. In the PK sub-study of the multi-center, randomized, placebo-controlled, double-blinded Riluzole in Spinal Cord Injury Study (RISCIS) Phase II/III trial (ClinicalTrials.gov Identifier: NCT01597518), a total of 32 SCI patients were enrolled, and most of our patients were middle-age Caucasian males with head and neck injuries. We studied the PK and pharmacodynamics (PD) of riluzole on motor recovery, measured by International Standards for Neurological Classification of SCI (ISNCSCI) Motor Score at injury and at 3-month and 6-month follow-ups, along with levels of the axonal injury biomarker phosphorylated neurofilament heavy chain (pNF-H), during the 2-week treatment. PK modeling, PK/PD correlations were developed to identify the potential effective exposure of riluzole for intended PD outcomes. The longitudinal impacts of SCI on the PK of riluzole are characterized. A time-varying population PK model of riluzole is established, incorporating time-varying clearance and volume of distribution from combined data of Phase I and Phase II/III trials. With the developed model, a rational, optimal dosing scheme can be designed with time-dependent modification to preserve the required therapeutic exposure of riluzole. The PD of riluzole and the relationship between PK and neurological outcomes of the treatment were established. The time course of efficacy in total motor score improvement (ΔTMS) and pNF-H were monitored. A three-dimensional (3D) PK/PD correlation was established for ΔTMS at 6 months with overall riluzole exposure area under the curve for Day 0-Day14 (AUCD0-D14) and baseline TMS for individual patients. Patients with baseline TMS between 1 and 36 benefited from the optimal exposure range of 16-48 mg*h/mL. The PD models of pNF-H revealed the riluzole efficacy, as treated subjects exhibited a diminished increase in progression of pNF-H, indicative of reduced axonal breakdown. The independent parameter of area between effective curves (ABEC) between the time profiles of pNF-H in placebo and treatment groups was statistically identified as a significant predictor for the treatment effect on the biomarker. A mechanistic clinical outcomes (CO)/PD (pNF-H) model was established, and the proposed structure demonstrated the feasibility of PK/PD/CO correlation model. No appreciable hepatic toxicity was observed with the current riluzole treatment regimen. The development of effective treatment for SCI is challenging. However, the future model-informed and PK-guided drug development and regimen modification can be rationally executed with the optimal dosing regimen design based on the developed 3D PK/PD model. The PK/PD/CO model can serve as a rational guide for future drug development, PKPD model refinement, and extension to other studies in SCI settings.

Modelling the time course of antimalarial parasite killing: a tour of animal and human models, translation and challenges

Malaria remains a global public health concern and current treatment options are suboptimal in some clinical settings. For effective chemotherapy, antimalarial drug concentrations must be sufficient to remove completely all of the parasites in the infected host. Optimized dosing therefore requires a detailed understanding of the time course of antimalarial response, whilst simultaneously considering the parasite life cycle and host immune elimination. Recently, the World Health Organization (WHO) has recommended the development of mathematical models for understanding better antimalarial drug resistance and management. Other international groups have also suggested that mechanistic pharmacokinetic (PK) and pharmacodynamic (PD) models can support the rationalization of antimalarial dosing strategies. At present, artemisinin-based combination therapy (ACT) is recommended as first line treatment of falciparum malaria for all patient groups. This review summarizes the PK-PD characterization of artemisinin derivatives and other partner drugs from both preclinical studies and human clinical trials. We outline the continuous and discrete time models that have been proposed to describe antimalarial activity on specific stages of the parasite life cycle. The translation of PK-PD predictions from animals to humans is considered, because preclinical studies can provide rich data for detailed mechanism-based modelling. While similar sampling techniques are limited in clinical studies, PK-PD models can be used to optimize the design of experiments to improve estimation of the parameters of interest. Ultimately, we propose that fully developed mechanistic models can simulate and rationalize ACT or other treatment strategies in antimalarial chemotherapy.

Understanding the hysteresis loop conundrum in pharmacokinetic/pharmacodynamic relationships

Hysteresis loops are phenomena that sometimes are encountered in the analysis of pharmacokinetic and pharmacodynamic relationships spanning from pre-clinical to clinical studies. When hysteresis occurs it provides insight into the complexity of drug action and disposition that can be encountered. Hysteresis loops suggest that the relationship between drug concentration and the effect being measured is not a simple direct relationship, but may have an inherent time delay and disequilibrium, which may be the result of metabolites, the consequence of changes in pharmacodynamics or the use of a non-specific assay or may involve an indirect relationship. Counter-clockwise hysteresis has been generally defined as the process in which effect can increase with time for a given drug concentration, while in the case of clockwise hysteresis the measured effect decreases with time for a given drug concentration. Hysteresis loops can occur as a consequence of a number of different pharmacokinetic and pharmacodynamic mechanisms including tolerance, distributional delay, feedback regulation, input and output rate changes, agonistic or antagonistic active metabolites, uptake into active site, slow receptor kinetics, delayed or modified activity, time-dependent protein binding and the use of racemic drugs among other factors. In this review, each of these various causes of hysteresis loops are discussed, with incorporation of relevant examples of drugs demonstrating these relationships for illustrative purposes. Furthermore, the effect that pharmaceutical formulation has on the occurrence and potential change in direction of the hysteresis loop, and the major pharmacokinetic / pharmacodynamic modeling approaches utilized to collapse and model hysteresis are detailed.

On the use of multivariable piecewise-linear models for predicting human response to anesthesia

The standard modeling paradigm used to describe the relationship between input anesthetic agents and output patient endpoint variables are single-input single-output pharmacokinetic-pharmacodynamic (PK-PD) compartment models. In this paper, we propose the use of multivariable piecewise-linear models to describe the relations between inputs that include anesthesia, surgical stimuli and disturbances to a variety of patient output variables. Subspace identification methods are applied to clinical data to construct the models. A comparison of predicted and measured responses is completed, which includes predictions from PK-PD models, and piecewise-linear time-invariant models.

Pharmacodynamic models for the cardiovascular effects of moxonidine in patients with congestive heart failure

AIMS

To assess the pharmacodynamics of moxonidine in patients with functional NYHA Class II-III congestive heart failure (CHF).

METHODS

A parallel population pharmacokinetic/pharmacodynamic (PK/PD) analysis was performed to assess the effect of moxonidine (0.1, 0.2, 0.3 mg twice daily) and placebo treatment on plasma noradrenaline (NA) levels, standing systolic blood pressure (SBP), and heart rate (HR) over 12 weeks in 97 patients with CHF using a parallel group design with dose escalation. A sequential analysis was also developed, where the relative changes in NA concentration were related to both SBP and HR.

RESULTS

In the parallel PD analysis, an effect delay was shown for all three end points (NA, SBP, and HR). An inhibitory Emax model was used to characterize the concentration-effect relationships. For SBP and HR, the EC50 value increased over time. For NA, there was a positive baseline drift over the 12 weeks; this was interpreted as disease progression. Moxonidine delayed this increase by 9.8 weeks. For SBP, there was a circadian pattern at baseline. In the sequential PD analysis, the relationship between the drug response (NA) and SBP or HR was best described by an inhibitory Emax model. No effect delays between the response and effects were found.

CONCLUSIONS

Effects of moxonidine on NA, SBP, and HR could be quantified by an effect compartment model in the presence of disease progression and circadian variations. Disease progression, as judged by increasing NA levels with time, was delayed by moxonidine. A direct relationship was found between NA and SBP/HR.

Target concentration intervention: beyond Y2K

Target concentration intervention (TCI) is proposed as an alternative conceptual strategy to therapeutic drug monitoring (TDM). It is argued that the idea of a therapeutic range has limited the interpretation of measured drug concentrations and diminished the anticipated clinical benefit to patients by use of an oversimplified pharmacodynamic model. TCI on the other hand embraces pharmacokinetic and pharmacodynamic concepts and uses the idea of a target effect and associated target concentration to make rational individual dose decisions.

Pharmacokinetic-pharmacodynamic model for educational simulations

Pharmacokinetic-pharmacodynamic (PK-PD) models play an important role in educational simulations. The parameters of PK-PD models described in the scientific literature are obtained from studies in which the drug concentrations and the drug-effect data are measured simultaneously. Simultaneous PK-PD studies cannot be expected to incorporate all possible combinations of drugs and patient physiology that are desired for educational simulations. To solve this problem, we elaborate on the traditional simultaneous PK-PD model, creating a new model that accepts parameter data from different, more readily available, nonsimultaneous pharmacologic studies. These data are incorporated in the model using a novel estimation procedure for the parameters kc0 and EC50. A sensitivity analysis of the parameter estimation procedure confirms that the time of peak effect following a bolus and the dose-response curve are accurately reflected by the new model. It also demonstrates how inconsistencies among the different parameter sets affect simulation of the recovery phase. The model is extended to incorporate any monotonic parametric or nonparametric dose-response curve. For the neuromuscular relaxant vecuronium, we demonstrate that data from different pharmacologic studies are available, and that the described estimation procedure leads to parameter estimates that are within the standard deviations of the parameters determined in a simultaneous PK-PD study.

Characteristics of indirect pharmacodynamic models and applications to clinical drug responses

This review describes four basic physiologic indirect pharmacodynamic response (IDR) models which have been proposed to characterize the pharmacodynamics of drugs that act by indirect mechanisms such as inhibition or stimulation of the production or dissipation of factors controlling the measured effect. The principles underlying IDR models and their response patterns are described. The applicability of these basic IDR models to characterize pharmacodynamic responses of diverse drugs such as inhibition of gastric acid secretion by nizatidine and stimulation of MX protein synthesis by interferon alpha-2a is demonstrated. A list of other uses of these models is provided. These models can be readily extended to accommodate additional complexities such as nonstationary or circadian baselines, equilibration delay, depletion or accumulation of a precursor pool, sigmoidicity, or other mechanisms. Indirect response models which have a logical mechanistic basis account for time-delays in many responses and are widely applicable in clinical pharmacology.

Influence of input rates on (+/-)-isradipine haemodynamics and concentration-effect relationship in healthy volunteers

Since the magnitude of the response to a drug may depend upon the drug input rate, the concentration-effect relationship of the new dihydropyridine (+/-)-isradipine was investigated using different administration modalities. Ten normotensive healthy volunteers were given, double-blind and in a crossover fashion, isradipine as a 1 mg iv infusion, 5 mg oral solution, 5 mg standard tablet, 10 mg slow release formulation, and a placebo. Blood pressure, heart rate, and plasma isradipine concentrations were recorded for 24 h. The maximal fall in diastolic blood pressure was similar after the infusion (-11.40 mmHg), the oral solution (-15.20 mmHg), and the standard tablet (-12.50 mmHg). In healthy volunteers the slow release form had no significant effect on blood pressure. The concentration-effect plots showed an increasing slope in the order infusion, solution, and table, and anticlockwise hysteresis. This was partly due to marked heart rate counter-regulation, the corresponding mean maximal heart rate increases being 24, 19, and 17 beats.min-1. The pronounced counter-regulation of the heart rate implies that a slow isradipine input rate would be more effective in decreasing blood pressure.

Algorithm to control "effect compartment" drug concentrations in pharmacokinetic model-driven drug delivery

In most computer-controlled pharmacokinetic model-driven drug infusion pumps, simulation of a linear compartmental pharmacokinetic model is used to compute the rate of intravenous drug infusion required to achieve setpoint central compartment (plasma) drug concentrations. For many drugs, it has been suggested that it is the drug concentration in a hypothetical "effect" compartment, rather than in the plasma, that should be manipulated to achieve maximum control over pharmacologic action. Controlling the effect compartment drug concentration is algorithmically more difficult than controlling the central compartment drug concentration because of the time delay between administration of drug into the central compartment and its subsequent appearance in the effect compartment. Presented in this paper is a model-based dosing algorithm for use in pharmacokinetic model-driven drug infusion devices that target the theoretical effect compartment drug concentration.

Pharmacodynamics of methylprednisolone phosphate after single intravenous administration to healthy volunteers

The pharmacokinetics and pharmacodynamics of methylprednisolone were investigated after intravenous administration of methylprednisolone phosphate to healthy subjects at seven different doses (16 to 1000 mg). Forty different pharmacodynamic parameters were followed for 1 week. The pharmacodynamic data were analyzed as a function of time as well as cumulative effects in form of the areas under the effect-time curves. Statistically significant dose-dependent effects of methylprednisolone were observed for 15 pharmacodynamic parameters. Highly significant (P less than or equal to 0.0001) effects were increases in glucose levels, number of white blood cells, and segmented granulocytes as well as a decrease in the number of lymphocytes. For these pharmacodynamic effects an integrated pharmacokinetic/pharmacodynamic model was derived that translates the methylprednisolone plasma concentration-time profiles into effect-time profiles. This model allows prediction of pharmacodynamic effects for any single dose in the range studied at any time point.