Corticosteroid Pharmacodynamic Modeling: Osteocalcin Suppression by Prednisolone

Pharmaceutical Research -     

Pharmacokinetic/Pharmacodynamic Evaluation of Deflazacort in Comparison to Methylprednisolone and Prednisolone

Purpose. The pharmacokinetics and pharmacodynamics of deflazacort after oral administration (30 mg) to healthy volunteers were determined and compared with those of 20 mg of methylprednisolone and 25 mg of prednisolone. Methods. Methylprednisolone, prednisolone and the active metabolite of deflazacort, 21-desacetyldeflazacort, were measured in plasma using HPLC. For the assessment of pharmacodynamics, differential white blood cell counts were obtained over 24 hours. An integrated pharmacokinetic-pharmacodynamic (PK-PD) model was applied to link corticosteroid concentrations to the effect on lymphocytes and granulocytes. Results. Deflazacort is an inactive prodrug which is converted rapidly to the active metabolite 21-desacetyldeflazacort. Maximum concentrations of 21-desacetyldeflazacort averaged 116 ng/ml and were observed after 1.3 h. The average area under the curve was 280 ng/ml · h, and the terminal half-life was 1.3 h. 21-Desacetyldeflazacort was cleared significantly faster than both methylprednisolone and prednisolone. The PK-PD-model was suitable to describe time course and magnitude of the observed effects. The results were consistent with reported values for glucocorticoid receptor binding affinities for the investigated compounds. Conclusions. Due to the short pharmacokinetic half-life of its active metabolite, pharmacodynamic effects of deflazacort are of shorter duration than those of methylprednisolone and prednisolone. The PK-PD model allows good prediction of pharmacodynamic effects based on pharmacokinetic and receptor binding data.     

Pharmacokinetic/Pharmacodynamic Modeling of Total Lymphocytes and Selected Subtypes After Oral Budesonide

In the present pharmacokinetic/pharmacodynamic (PK/PD) evaluation, cortisol, total lymphocytes, and lymphocyte subpopulations were monitored following single oral doses of two oral formulations of 3 mg budesonide (BUD) (Dosage Forms A and B) in order to assess the differential effects that BUD may have on cortisol suppression and the modulation of blood lymphocyte subtypes. On a single occasion, five subjects received one 3 mg capsule of Dosage Form A, four received three capsules of Dosage Form A (single dose of 9 mg), and five received three capsules of Dosage Form B (single dose of 9 mg). Placebo capsules were administered to six subjects in the study. Cortisol concentrations, total lymphocyte counts, and lymphocyte subpopulation counts for the CD3, CD4, CD8, CD19, and CD56/16 were fitted to an in direct PK/PD response model that described the effects of BUD on serum cortisol concentrations as well as the combined effects of BUD and cortisol on total lymphocytes and the CD3, CD4, CD8, CD19, and CD56/16 subtypes. Data were also analyzed using noncompartmental methods. The PK/PD model fitted the data with the exception of data for CD56/16. The IC₅₀ value for unbound BUD acting on total lymphocytes was 0.276 ng/ml while the IC₅₀ values for unbound BUD acting on lymphocyte subtypes ranged from 0.150 ng/ml for CD4 to 0.364 ng/ml for CD8. The IC₅₀ values for the effects of BUD on serum cortisol were lower (0.079 ng/ml). The results of PK/PD modeling and noncompartmental analysis indicate that BUD has a smaller effect on the CD8 subtype and larger effects on the CD4 and CD19 subtypes, relative to the effect on total lymphocytes, and that cortisol suppression, although not a direct immunological biomarker, may be a more sensitive marker for the systemic effect of corticosteroids.     

Pharmacodynamic Assessment of the Benztropine Analogues AHN-1055 and AHN-2005 Using Intracerebral Microdialysis to Evaluate Brain Dopamine Levels and Pharmacokinetic/Pharmacodynamic Modeling

No HeadingPurpose. The benztropine (BZT) analogues bind with high affinity to the dopamine transporter (DAT) and demonstrate a behavioral and pharmacokinetic profile unlike that of cocaine. The development of a predictive pharmacokinetic/pharmacodynamic (PK/PD) model to characterize the concentration-effect relationship between the BZT analogues and brain dopamine (DA) levels is an important step in the evaluation of these compounds as potential cocaine abuse pharmacotherapies. Hence, the objective of this study was to mathematically characterize the PD of BZT analogues and cocaine, using appropriate PK/PD models.Methods. Dialysis probes were stereotaxically implanted into the nucleus accumbens of Sprague-Dawley rats (275–300 g). Extracellular fluid (ECF) DA levels were measured after intravenous administration of the BZT analogues AHN-1055 and AHN-2005, as well as cocaine using high performance liquid chromatography-electrochemical detection (HPLC-ECD). PD models were used to describe the relationship between the BZT analogues or cocaine and brain microdialysate DA, and suitability was based on standard goodness-of-fit criteria.Results. The BZT analogues produced a sustained increase in brain microdialysate DA levels in comparison to cocaine. The time of maximum concentration (T_max) for brain microdialysate DA was 2 h for AHN-1055 and 1 h for AHN-2005 compared to a T_max of 10 min for cocaine. The duration of brain microdialysate DA elevation was 12–24 h for the BZTs in comparison to 1 h for cocaine. An indirect model with inhibition of loss of response and a sigmoid E_max model best described the PK/PD for the BZT analogues and cocaine, respectively. The 50% of maximum inhibition (IC₅₀) of the loss of DA was lower for AHN-2005 (226 ± 27.5 ng/ml) compared to AHN-1055 (321 ± 19.7 ng/ml). In addition, the EC₅₀ for cocaine was 215 ± 11.2 ng/ml.Conclusions. The slow onset and long duration of BZT analogue–induced DA elevation may avoid the reinforcing effects and craving of cocaine. Further, the developed models will be useful in characterizing the PK/PD of other analogues and aid in the assessment of the therapeutic efficacy of the BZT analogues as substitute medications for cocaine abuse.     

Pharmacokinetic and Pharmacodynamic Evaluation for Tissue-Selective Inhibition of Cholesterol Synthesis by Pravastatin

The tissue-selective inhibition of cholesterol synthesis by pravastatin was evaluated pharmacokinetically and pharmacodynamically. Plasma, tissue, urine, and bile concentrations were measured after iv bolus injection of pravastatin to rats at various doses. The total body clearance and steady state volume of distribution decreased with increasing dose. A saturable biliary excretion was also observed. The time course of plasma and liver concentrations was described by a three-compartment model, consisting of a central compartment, a deep compartment with an nonsaturable uptake process, and a shallow compartment with saturable uptake and nonsaturable elimination processes. It suggests that a mechanism for the decrease in the total body clearance and distribution volume might be explained by a saturation of pravastatin uptake into the liver. Plasma concentration data after oral administration was also fitted to the same model by connecting an absorption compartment to the shallow compartment. The inhibitory activity of pravastatin against cholesterol synthesis in liver could be related to the concentration in the shallow compartment via a sigmoidal Emax model and the obtained pharmacodynamic parameters were comparable to those in vitro. Results suggest that the carrier-mediated hepatic uptake of pravastatin is actually responsible for the hepatoselective inhibition of cholesterol synthesis under physiological conditions.     

Pharmacodynamic Modeling of Prednisolone Effects on Natural Killer Cell Trafficking

Pharmaceutical Research -     

Sex Specificity in Methadone Analgesia in the Rat: A Population Pharmacokinetic and Pharmacodynamic Approach

Purpose. To quantify the extent to which a sex-specific dichotomy in the temporal evolution of the analgesic effect, after intravenous (i.v.) methadone injection in the rat, relates to the pharmacokinetics (PK) and pharmacodynamics (PD) that mediate the dose-to-effect pathway. Methods. Tail-flick analgesia was measured after i.v. methadone injection (0.35 mg/kg) in female (n = 16) and male (n = 16) Sprague-Dawley rats. The PK were evaluated in separate female (n = 56) and male (n = 56) rats after they had received the same dose of methadone i.v. (0.35 mg/kg). A bicompartmental model described the kinetics and a sigmoid E_max model-related drug effect vs. simulated concentrations (pharmacodynamics) at the times of effect measurement. All model parameters as well as interanimal and assay variabilities were estimated with a mixed-effects population method using the program NONMEM. Results. The area under the effect-time curve (AUCE_0-120) was (mean ± interanimal SD) 1859 ± 346 min in the females, which was significantly lower than the 4871 ± 393 min in the males (P < 0.0001). On the contrary, the profiles of concentration vs. time were higher in females and, therefore, corresponded inversely to the effect vs. time-relative magnitudes. The central volume of distribution, V1, was 1.94 ± 0.37 l/kg for female rats and 3.01 ± 0.33 l/kg for male rats. Also, the central clearance was 0.077 ± 0.006 l/min/kg and 0.102 ± 0.005 l/min/kg, respectively, for female and male rats. Both parameters differed significantly between sexes (P < 0.0001). The pharmacodynamic maximum observed effect parameter (E_max) was 37% ± 29% in female rats and 85% ± 16% in male rats, and these values were significantly different (P < 0.0001). The parameter for the concentration eliciting half of E_max (EC₅₀) was 24.1 ± 7.5 g/l in female rats and 20.3 ± 2.9 g/l in male rats, and the Hill-related exponent, , was 6.3 ± 3.9 in female rats and 5.5 ± 4.1 in male rats. These parameters did not differ significantly (at the P < 0.05 level). Conclusions. A sex-specific dichotomy in the methadone antinociceptive effect, in the rat, was not proportionally related to plasma concentrations. Each sex corresponded to a distinct subpopulation of the PK parameters and one of the pharmacodynamic parameters (E_max). When the course of a drug involves PK or PD subpopulations, PK/PD modeling can afford the safest prediction of the effect-time evolution for a particular dose.     

Evaluating Pharmacokinetic/Pharmacodynamic Models Using the Posterior Predictive Check

The posterior predictive check (PPC) is a model evaluation tool. It assigns a value (p ^PPC ) to the probability that the value of a given statistic computed from data arising under an analysis model is as or more extreme than the value computed from the real data themselves. If this probability is too small, the analysis model is regarded as invalid for the given statistic. Properties of the PPC for pharmacokinetic (PK) and pharmacodynamic (PD) model evaluation are examined herein for a particularly simple simulation setting: extensive sampling of a single individual's data arising from simple PK/PD and error models. To test the performance characteristics of the PPC, repeatedly, real data are simulated and for a variety of statistics, the PPC is applied to an analysis model, which may (null hypothesis) or may not (alternative hypothesis) be identical to the simulation model. Five models are used here: (PK1) mono-exponential with proportional error, (PK2) biexponential with proportional error, (PK2) biexponential with additive error, (PD1) E _max model with additive error under the logit transform, and (PD2) sigmoid E _max model with additive error under the logit transform. Six simulation/analysis settings are studied. The first three, (PK1/PK1), (PK2/PK2), and (PD1/PD1) evaluate whether the PPC has appropriate type-I error level, whereas the second three (PK2/PK1), (PK2/PK2), and (PD2/PD1) evaluate whether the PPC has adequate power. For a set of 100 data sets simulated/analyzed under each model pair according to a stipulated extensive sampling design, the p ^PPC is computed for a number of statistics in three different ways (each way uses a different approximation to the posterior distribution on the model parameters). We find that in general; (i) The PPC is conservative under the null in the sense that for many statistics, prob(p ^PPC )< for small . With respect to such statistics, this means that useful models will rarely be regarded incorrectly as invalid. A high correlation of a statistic with the parameter estimates obtained from the same data used to compute the statistic (a measure of statistical sufficiency) tends to identify the most conservative statistics. (ii) Power is not very great, at least for the alternative models we tested, and it is especially poor with statistics that are in part a function of parameters as well as data. Although there is a tendency for nonsufficient statistics (as we have measured this) to have greater power, this is by no means an infallible diagnostic. (iii) No clear advantage for one or another method of approximating the posterior distribution on model parameters is found.     

Inotropic Effect of Digoxin in Humans: Mechanistic Pharmacokinetic/Pharmacodynamic Model Based on Slow Receptor Binding

Purpose. The purpose of this study was to construct a mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model for digoxin that describes the relationship between plasma concentration and inotropic response. Methods. On the basis of results obtained in the isolated perfused rat heart, a PK/PD model for digoxin in humans was developed. In fitting the model to previously published bolus dose and concentration clamp data (shortening of electromechanical systole), the plasma concentration-time curves were used as forcing functions in the computer program ADAPT II. Results. The mechanistic approach allowed a modeling of digoxin pharmacodynamics which is consistent with available inotropic response data. The estimates of the receptor binding parameters were in the same order of magnitude as those measured in vitro for ouabain. The mechanistic model explained the parameters of the empirical link model (EC₅₀, E _max and delay time ) in terms of the underlying processes, suggesting that the long equilibration half-time of 13 h is due to slow receptor binding. The empirical link model, in contrast, is not compatible with a noninstantaneous receptor binding process and led to estimates of the delay time that were dependent on the digoxin administration schedule. Conclusions. The new, mechanistic model may provide a rationale for better understanding of digoxin pharmacodynamics and could become a tool to bridge the gap between in vitro and in vivo studies.     

A Pharmacokinetic/Pharmacodynamic Approach to Predict Total Prednisolone Concentrations in Human Plasma

Prednisolone and prednisone are two widely used corticosteroids for various inflammatory and immune diseases. Prednisolone is the active form of prednisone in vivo. Total prednisolone in plasma exhibits nonlinear pharmacokinetics mainly due to its nonlinear protein binding. Other factors such as reversible metabolism (or interconversion between prednisolone and prednisone), competitive protein binding from endogenous cortisol, cortisol circadian rhythm, and prednisolone mediated cortisol suppression complicate prednisolone pharmacokinetics. This study was aimed to develop a new approach to describe the nonlinear pharmacokinetics of total prednisolone and predict total prednisolone concentrations in plasma. Based on literature datasets, a linear two-compartment pharmacokinetic model was developed to adequately describe the reversible metabolism between free prednisone and prednisolone. Cortisol and prednisolone protein binding were described via the sum of a Langmuir and linear type binding. The endogenous cortisol circadian rhythm and cortisol suppression during prednisone or prednisolone exposure were described with a previously reported linear release rate pharmacokinetic/pharmacodynamic (PK/PD) model. By combining the pharmacokinetic models for free prednisone and prednisolone, the linear release rate model for cortisol suppression, and competitive protein binding between cortisol and prednisolone, we were able to predict total prednisolone concentrations in plasma. The predicted total prednisolone concentrations in plasma were in good agreement with the literature reported data. Thus, this PK/PD approach shows that the combination of nonlinear protein binding, cortisol circadian rhythm, and cortisol suppression could account for the nonlinearity of total prednisolone. In addition, it also allows a valid prediction of total prednisolone in plasma after either prednisone or prednisolone administration.     

A Pharmacokinetic/Pharmacodynamic Approach to Predict the Cumulative Cortisol Suppression of Inhaled Corticosteroids

The suppression of endogenous cortisol release is one of the major systemic side effects of inhaled corticosteroids in the treatment of asthma. The circadian rhythm of the endogenous cortisol release and the resulting plasma concentrations as well as the release suppression during corticosteroid therapy could previously be described with an integrated PK/PD model. Based on this model, a PK/PD approach was developed to quantify and predict the cumulative cortisol suppression (CCS) as a surrogate marker for the systemic activity of inhaled corticosteroid therapy. The presented method was applied to predict CCS after single doses and during short-term multiple dosing of the inhaled corticosteroids flunisolide (FLU), fluticasone propionate (FP), and triamcinolone acetonide (TCA), and after oral methylprednisolone as systemic reference therapy. Drug-specific PK and PD parameters were obtained from previous single-dose studies and extrapolated to the multiple-dose situation. For single dosing, a similar CCS within the range of 16–21% was predicted for FP 250 g, FLU 500 g, and TCA 1000 g. For multiple dosing, a respective CCS of 28–33% was calculated for FLU 500 g bid, FP 250 g, bid, and TCA 1000 g bid. Higher cortisol suppression compared to these single and multiple dosing regimens of the inhaled corticosteroids was predicted after oral doses of only 1 mg and 2 mg methylprednisolone, respectively. The predictive power of the approach was evaluated by comparing the PK/PD-based simulations with data reported previously in clinical studies. The predicted CCS values were in good correlation with the clinically observed results. Hence, the presented PK/PD approach allows valid predictions of CCS for single and short-term multiple dosing of inhaled corticosteroids and facilitates comparisons between different dosing regimens and steroids.     

Receptor-mediated pharmacokinetic/pharmacodynamic model of interferon-beta 1a in humans

Purpose. An integrated receptor-based pharmacokinetic/pharmaco- dynamic (PK/PD) model of interferon- 1a (IFN- 1a) previously developed for monkeys was used to capture the time-course of drug and induced neopterin concentrations after intravenous (IV) and subcutaneous (SC) dosing in humans. Methods. Data were extracted from the literature by digitalization. Single-dose (3 IV doses and 1 SC dose) PK/PD profiles were simultaneously fitted using the basic model and the ADAPT II computer program. Additional submodels incorporating neutralizing antibody formation and negative feedback inhibition were applied to account for drug accumulation and lower than expected neopterin concentrations encountered after multiple-dosing (1 SC dose every 48 hs). Results. The basic model jointly-captured the nonlinear PK behavior of the drug and induced neopterin concentrations after all single doses. Slow and incomplete absorption (F = 0.33) of the SC dose resulted in prolonged drug concentrations reflective of flip-flop kinetics. Despite lower drug concentrations, SC dosing produced a similar neopterin profile as compared with the IV doses; however, with a longer time to peak effect and slightly higher neopterin concentrations at later time points. The PD component of the model represents a modified precursor-dependent indirect response model driven by the amount of internalized drug-receptor complex. The latter stimulated a 6-fold increase in the production of the neopterin precursor (S_max = 5.89). Drug accumulation and lower than expected neopterin concentrations after multiple dosing were also captured after the inclusion of the submodels. Conclusions. The present integrated PK/PD model for IFN- 1a is mechanistic in nature with receptor-mediated disposition and dynamics and was successfully applied to human clinical data.     

Pharmacokinetic/pharmacodynamic evaluation of inhalation drugs: application to targeted pulmonary delivery systems

Inhaled therapy with either glucocorticoids and/or β₂-adrenergic drugs remains the mainstay of asthma treatment. In the last few years, a number of new products have been introduced into the market with the goal of improving efficacy and safety. This review article summarises the pharmacokinetic and pharmacodynamic properties of inhaled drugs for topical delivery necessary to achieve this goal. Pharmacokinetic properties include a high pulmonary deposition, low oral bioavailability, optimised pulmonary residence time and a very high systemic clearance. Optimisation of pharmacodynamic properties, such as receptor selectivity, may also yield drugs with improved pulmonary selectivity. As existing drugs also provide high efficacy and safety profiles, future developments will represent only slight improvements and quantum leap improvements are unlikely to occur.     

Population Pharmacokinetic/ Pharmacodynamic Modeling of Cetrorelix,a Novel LH-RH Antagonist,and Testosterone in Rats and Dogs

Purpose. Population models for thepharmacokinetic-pharmacodynamic relationship for cetrorelix (CET), a luteinising hormone-releasinghormone (LH-RH) antagonist, and the pharmacodynamic response ontestosterone production were investigated in rats and dogs. Methods. The plasma concentrations of CET and testosterone weredetermined after intravenous and subcutaneous injections. Thepopulation PK/PD-models were developed using P-PHARM software. Results. Absolute bioavailability of cetrorelix was 100% in rats and97% in dogs. In rats, the pharmacokinetics was explained by atwo-compartment model with saturable absorption, while athree-compartment model was used in dogs. Testosterone suppression in both specieswas described by a sigmoid E_max model with maximum effect (E_max)considered as total hormonal suppression. The duration of testosteronesuppression in rats was longer at higher doses. The populationelimination half-lifes after iv-dose were 3.0 h in rats and 9.3 h in dogs.Population mean estimates of IC₅₀ were 1.39 and 1.24 ng/ml in ratsand dogs, respectively. Conclusions. A population pharmacokinetic model was developed toexplain the dissolution rate limited absorption from the injection site.The suppression of testosterone could be described by an indirectinhibitory sigmoid E_max model. In both species 1-2 ng/ml CET inplasma was necessary to suppress testosterone production.     

Simultaneous vs. Sequential Analysis for Population PK/PD Data II: Robustness of Methods

A model can be fit to joint PK/PD data (concentration and effect) either simultaneously or sequentially. The results of a companion paper suggested that when the data-analytic and true models agree, a particular sequential approach is computationally faster than the simultaneous one, yet produces hardly less precise PD parameter estimates, and for suitable designs, about as accurate PD standard error estimates. In this paper, we compare the performance of various methods for the case that the data-analytic model is misspecified. We illustrate these methods by applying them to a set of real data. Using NONMEM, population PK/PD observations were simulated under various study designs according to a one- or two-compartment PK model and direct Emax or sigmoid Emax model. A one-compartment PK model and Emax PD model were fit to the simulated observations by simultaneous and sequential methods. Predictive performance (interpolation and extrapolation) of PD and the type-I error rate of a likelihood ratio test are compared. The real data set consists of PK and (more frequent) PD observations after administration of the muscle relaxant vecuronium. When only the PK data-analytic model is misspecified, the simultaneous method has greater precision than the sequential methods. However a sequential method that uses a non-parametric PK model performs better than both other methods when PK model misspecification is severe. When the PD data-analytic model is misspecified, sequential and simultaneous methods perform similarly. The analysis of the real data shows that the PK fitted with the simultaneous method can be quite sensitive to PD model misspecification, yielding a possible diagnostic for this type of misspecification.     

Simultaneous vs. Sequential Analysis for Population PK/PD Data I: Best-Case Performance

Dose [-concentration]-effect relationships can be obtained by fitting a predictive pharmacokinetic (PK)-pharmacodynamic (PD) model to both concentration and effect observations. Either a model can befit simultaneously to all the data ("simultaneous" method), or first a model can befit to the PK data and then a model can be fit to the PD data, conditioning in some way on the PK data or on the estimates of the PK parameters ("sequential" method). Using simulated data, we compare the performance of the simultaneous method with that of three sequential method variants with respect to computation time, estimation precision, and inference. Using NONMEM, under various study designs, observations of one type of PK and one type of PD response from different numbers of individuals were simulated according to a one-compartment PK model and direct Emax PD model, with parameters drawn from an appropriate population distribution. The same PK and PD models were fit to these observations using simultaneous and sequential methods. Performance measures include computation time,fraction of cases for which estimates are successfully obtained, precision of PD parameter estimates, precision of PD parameter standard error estimates, and type-I error rates of a likelihood ratio test. With the sequential method, computation time is less, and estimates are more likely to be obtained. Using the First Order Conditional Estimation (FOCE) method, a sequential approach that conditions on both population PK parameter estimates and PK data, estimates PD parameters and their standard errors about as well as the "gold standard" simultaneous method, and saves about 40% computation time. Type-I error rates of likelihood ratio test for both simultaneous and sequential approaches are close to the nominal rates.