Nephrol Dial Transplant (2003) 18: 1201-1208
© 2003 European Renal Association-European Dialysis and Transplant Association
A compartmental pharmacokinetic model of cyclosporin and its predictive performance after Bayesian estimation in kidney and simultaneous pancreaskidney transplant recipients
1 Department of Clinical Pharmacy and Toxicology, 2 Department of Nephrology, 3 Centre for Human Drug Research and 4 Department of Clinical Chemistry, Leiden University Medical Center, Leiden, The Netherlands
| Abstract |
|---|
|
|
|---|
Background. Therapeutic drug monitoring of cyclosporin A (CsA) is an obvious necessity because of its unpredictable absorption and narrow therapeutic window. The use of limited sampling models (LSMs) has improved the estimation of the systemic exposure [area under curve (AUC)] compared with C0h monitoring, but these equations are rigid and not reliable in patients with an abnormal absorption profile. We developed and validated a limited sampling (t=0, 2 and 3 h) strategy, based on a compartmental population pharmacokinetic (PK) model for CsA after kidney transplantation alone (KTA) and simultaneous pancreaskidney transplant (SPKT) recipients, a group of patients with unpredictable absorption kinetics.
Methods. A two-compartment model with lag time and first-order absorption was calculated using a PK software package from data of 20 KTA and SPKT recipients and validated prospectively in 20 KTA and 20 SPKT recipients. Calculated population PK parameters were individualized for each of the remaining 40 patients based on their CsA dosing and on one or a combination of measured CsA blood concentrations using the Bayesian fitting method. AUCs were calculated from individualized PK parameters. AUCs were also calculated using previously published LSMs. Relationships between AUCs calculated by the models and the golden standard AUC (trapezoidal rule) were investigated by Pearson correlation test.
Results and conclusions. A population two-compartment model is presented to reliably estimate the CsA AUC in KTA and SPKT recipients. The performance of the model to estimate the AUC is comparable to the performance of two published LSMs in KTA patients, but markedly better in SPKT patients. Combined with Bayesian fitting, the model is very flexible since sampling times are not rigid and can be varied as long as dosing and sampling times are recorded accurately. The model has already proven to be clinically useful and is currently used to further investigate CsA in an integrated pharmacokinetic/pharmacodynamic model.
Keywords: Bayesian; cyclosporin; kidney transplant; pharmacokinetics; simultaneous pancreaskidney transplant
| Introduction |
|---|
|
|
|---|
Since the introduction of cyclosporin A (CsA), it was apparent that defining the appropriate dose was problematic because of its poor and unpredictable absorption and narrow therapeutic window. Therapeutic drug monitoring (TDM) is therefore an obvious necessity and, by convention, dosing adjustments became targeted to trough levels (C0h), i.e. the CsA blood concentration immediately before the next scheduled dose. Although the correlation between C0h level and its efficacy to prevent acute rejection or chronic nephrotoxicity is poor, this parameter is still widely used to guide CsA dosing after transplantation [1,2]. Despite maintaining CsA trough levels within the therapeutic range, substantial groups of patients experience either acute rejection episodes or nephrotoxicity [3,4]. The CsA area under the drug concentration vs time curve (AUC) seems a better measure of systemic drug exposure and efficacy to prevent acute rejection episodes [5], but AUC monitoring has not gained popularity, largely because of the inconvenience of multiple blood samplings over a 12 h period and cost considerations.
Because AUC reflects systemic drug exposure, several groups have used sparse-sampling algorithms [limited sampling models (LSMs)] as a way of predicting AUCs without the need for large numbers of blood-level measurements [6,7]. Exposure during the first 4 h (AUC04h), the area of greatest inter- and intrapatient variability, is a reliable estimate of total drug exposure throughout the dosing interval (AUC012h) [8]. Recently, dosing based on a single drug level measurement, i.e. C2h, within this absorption phase has been associated with improved results after de novo heart, liver and kidney transplantation [9]. It should be stressed, however, that estimation of the true AUC will always be more reliable when more samples are taken. This may be particularly true, or even mandatory, in patients with an unpredictable absorption profile, such as diabetic recipients [10].
Compared with LSMs, a compartmental population pharmacokinetic model for CsA in renal transplant recipients combined with the maximum a posteriori Bayesian fitting method offers in clinical practice the important advantage of flexibility [11,12]. The exact sampling time is no longer an issue as long as the exact times of drug administration and blood drawing are recorded. In the present study we describe the development and prospective validation of such a model in non-diabetic patients who received a kidney transplant alone (KTA) and simultaneous pancreaskidney transplant (SPKT) recipients.
| Subjects and methods |
|---|
|
|
|---|
Patients
Sixty (32 KTA and 28 SPKT) transplant recipients (38 male, 22 female) treated with CsA-based immunosuppression were studied. Thirty patients were early (<3 months) post-transplantation, while the other 30 had been transplanted at least 3 months before the present study. All patients received the microemulsion formulation of CsA (Neoral; Novartis Pharmaceuticals, Basel, Switzerland) with prior dosing based on trough levels. Patients who used drugs known to interact with CsA pharmacokinetics were excluded. Patient characteristics are listed in Table 1
|
Pharmacokinetics
Using the Kinpop module of the pharmacokinetic software package MW/Pharm version 3.33 (Mediware, Groningen, the Netherlands), a population two-compartment model with a lag time and first-order absorption pharmacokinetics was calculated from the CsA dosing and the blood concentration values of the first 20 patients (12 KTA and 8 SPKT recipients). This program uses an iterative two-stage Bayesian procedure and calculates means, medians and SDs of the pharmacokinetic parameters [11]. During the iterative two-stage Bayesian procedure, pharmacokinetic parameters were set to be distributed log-normally and bioavailability was fixed at 0.5.
The calculated mean population pharmacokinetic parameters were individualized for each of the remaining 40 patients (20 non-diabetic KTA and 20 SPKT recipients) based on their CsA dosing and one or a combination of measured blood concentrations (0 h; 1 h; 2 h; 3 h; 4 h; 0+1 h; 0+2 h; 0+3 h; 0+4 h; 2+4 h; 0+1+2 h; 0+1+3 h; 0+2+3 h; 1+2+3 h; 0+1+2+3 h; 0+1+2+4 h; 0+1+2+3+4+6+8+12 h) according to the maximum a posteriori (MAP) Bayesian fitting method [12], using the MW/Pharm computer program. By means of MAP Bayesian fitting any available information, i.e. a priori population parameters, drug dosage regimen and measured blood concentrations can be used to estimate the a posteriori pharmacokinetic parameters of the individual patients. These a posteriori pharmacokinetic parameters of the individual patient are the maximum likelihood estimates obtained by MAP Bayesian fitting, minimizing the deviations of measured and predicted concentrations and of population pharmacokinetic parameters and pharmacokinetic parameters of the individual patient [12]. This approach is very flexible and it ensures an optimal use of available information, both from a population and from the individual patient. From the individualized pharmacokinetic parameters the area under the CsA blood concentration time curve (AUC012h) was calculated for each combination of measured blood concentrations.
The AUC012h of the 40 patients was also calculated using the two-point limited sampling strategy as described by Wacke et al. [7]:
|
| (001) |
|
| (002) |
Drug analysis
Whole blood concentrations of CsA in the first 20 patients were determined by radioimmunoassay (RIA) (Cyclotrac; IncStar, Stillwater, MN). Because of a change of equipment, whole blood concentrations of CsA in the remaining 40 patients were determined by fluorescence polarization immunoassay (FPIA) (Axsym; Abbott Diagnostics, Abbott Park, IL).
Statistical analysis
Statistical analysis was performed three times. First, the KTA recipient group was investigated. Secondly, the SPKT recipient group was investigated. Finally, these two groups were analysed together.
The AUCs calculated by the different methods were compared with the golden standard AUC by linear regression analysis and Pearson correlation coefficient. Calculations were carried out by means of the SPSS software (version 9.0). Predictive performance of the different methods was also investigated by calculating the prediction precision and bias according to Sheiner and Beal [14]. Prediction bias was calculated as the mean prediction error (MPE), i.e. the mean of differences between the AUC according to the different methods and the golden standard AUC. Prediction precision was calculated as the mean absolute prediction error (MAPE), i.e. the mean of the absolute differences between the AUC according to the several different methods and the golden standard AUC. Smaller values for MPE and MAPE indicate less bias and greater precision.
The actual description of blood concentration in time by the model was investigated prospectively by calculating residuals between measured and predicted CsA blood concentrations in the 40 patients. The concentrations were predicted for each patient using either the whole data set or only limited sampling at t=0, 2 and 3 h after drug administration.
| Results |
|---|
|
|
|---|
CsA pharmacokinetics was adequately described by a two-compartment model with a lag time. Pharmacokinetic parameters as calculated by Kinpop are listed in Table 2
|
|
In Figures 2
|
|
|
Bias and precision of the different models are listed in Tables 3
|
|
In KTA recipients (Table 3
In SPKT recipients (Table 4
), it appears that the first LSM (1, 3) is less useful to estimate the AUC012h. Bias, precision and correlation of the compartmental model are also worse in several combinations of blood concentration data compared with KTA recipients. The best relationships are found with the blood concentration data 0+3 h, 0+2+3 h and 0+1+2+3 h when only data within 3 h after drug administration are taken into consideration. Bias and precision of these combinations are good and comparable to the good performance of the second LSM (2, 6).
When the performance of the models is validated on both KTA and SPKT recipients together (data not shown), the combinations of blood concentration data that best describe the systemic exposure are 0+2 h, 0+2+3 h and 0+1+2+3 h. The performance of LSM (2, 6) is very good; the performance of LSM (1, 3) seems less useful.
Figure 5
displays the residuals of the model as a function of time after the MAP Bayesian fitting procedure based on all concentrations and on limited sampling at t=0, 2 and 3 h in the 20 KTA and 20 SPKT recipients. In both groups the concentrations are scattered around the x-axis, indicating a good description of the blood concentrations by the model. The largest deviations from the measured concentrations are around t=1 h. The peak concentration measured is usually reached at 1 or 2 h. Limited sampling at t=0, 2 and 3 h causes somewhat more deviation from the measured concentrations. However, the residuals remain more or less distributed around 0 µg/l, especially at t=1 and 2 h. The residuals at the other time points do not appear to change during limited sampling at 0, 2 and 3 h. The residuals in the SPKT recipients are also scattered around the x-axis. However, the pattern is somewhat different, indicating different pharmacokinetics in this group.
|
| Discussion |
|---|
|
|
|---|
At present, calcineurin inhibitors still constitute the cornerstone of immunosuppressive regimens for the prevention of allograft rejection in renal transplant recipients. Recent data have indicated that the long-term outcome of both cadaveric and live donor renal allografts has improved in the past decade, especially if acute rejection episodes could be prevented [15]. In healthy volunteers and stable renal transplant recipients, the CsA microemulsion formulation under the name Neoral has increased dose linearity with AUC and has resulted in a better correlation between C0h and AUC012h than cyclosporin gelcaps or Sandimmune [16]. However, in de novo transplant recipients, the interpatient variability of pharmacokinetic parameters is not different between these formulations during the first post-operative weeks and in clinical practice C0h has an equally poor correlation with AUC012h [17]. Conversely, AUC estimates correlate with early clinical events [18], but there is little information about the usefulness of AUC monitoring to improve long-term outcome [19].
Estimating the systemic exposure to CsA with LSMs, but also with C2h, is rigid. It is not allowed to deviate too much from the predefined time points. Moreover, once a time point is missed, a LSM formula becomes useless, since then one of the determinants in the formula is lacking. Deviations from predefined time points and missing of samples occur frequently in daily clinical practice. In contrast to the LSMs, these everyday problems are dealt with easily by a Bayesian fitting procedure. For example, the pharmacokinetic model can be used to calculate systemic exposure to CsA in a data set that consists of 0, 1.5 and 3.4 h time points instead of 0, 2 and 3 h or a data set that consists of only 2 and 3 h time points, illustrating the flexibility of the model and the clinical applicability.
The performance of the current population model to estimate the AUC012h was investigated using several sampling strategies. As expected, the estimation of the AUC012h always improves when the number of sampling points is increased, as more information is introduced. However, using the time points t=0, 2 and 3 h post-dose with the population model and MAP Bayesian fitting, the AUC was estimated with a mean accuracy of -2.3% and a mean precision of 8.3%. This is well within accepted and clinical relevant limits.
Using individualized population pharmacokinetic models, CsA exposure in renal transplant patients can be estimated well. However, intrapatient pharmacokinetics is known to change in the early period (68 weeks) after transplantation. The model was developed and validated in a population with varying time after transplantation. Consequently, the model is able to deal with the changes in pharmacokinetics after transplantation, mainly because this change has been incorporated into the range of the pharmacokinetic parameters. However, using the population pharmacokinetic model as described, the systemic exposure can be estimated, but actual prediction of the systemic exposure during a period of great changes in pharmacokinetics is not possible. In order to ensure an adequate systemic exposure, several approaches can be made. Either the estimation of the AUC of CsA is performed frequently in the early period after transplantation or a pharmacokinetic model is used in which factors specifically influencing this change in pharmacokinetics are incorporated. One of those factors is probably time after transplantation. However, clinically applicable pharmacokinetic models, in which time after transplantation has been incorporated, enabling the actual prediction of systemic exposure in a period of changing pharmacokinetics of CsA, are scarce. For our own approach, in which we thus estimate the AUC012h of CsA at a certain moment, we are currently investigating still the optimal frequency of drug monitoring, both short- and long-term after transplantation. In our prospective study we use every other week estimations of the AUC012h during the first 12 weeks after transplantation, which is followed by an estimation of the AUC012h every 36 months.
The study population consisted of both patients that received kidney transplant as well as SPKT recipients. Patients with insulin-dependent diabetes mellitus often suffer from (severe) gastrointestinal dysfunction and in these patients CsA pharmacokinetics is often variable and unpredictable [10]. Inclusion of these patients may have skewed the pharmacokinetic parameters of the population model, leading to deviation of predicted from measured CsA concentrations in the KTA patients. The advantage, however, is that the model is suitable to describe CsA pharmacokinetics in the population of interest, i.e. recipients with an abnormal absorption profile. The differences in CsA pharmacokinetics between KTA and SPKT recipients were adequately detected by the compartment model and by one of the previously described LSMs [6]. The performance of the other LSM was significantly worse [7], indicating that caution should be taken when monitoring CsA in various patient groups with defined sampling time-points. This is illustrated especially by the lack of correlation between C2h and AUC04h or AUC012h in SPKT recipients as compared with the patients without diabetes mellitus, who received a kidney transplant. In Figure 6
, two clinical examples are given illustrating the above-mentioned issues and how the model deals with them.
|
Although the drug is metabolized to at least 25 metabolites, blood concentrations of CsA after oral administration, as determined by the immunoassays, are adequately described by a two-compartment model with a lag time and first-order absorption pharmacokinetics. The model was superior to a one-compartment model, while an extra compartment did not add significantly to the description of the data. Fitting procedures by the calculation algorithm were improved by fixing the oral availability at 0.5. Absorption pharmacokinetics therefore is characterized mainly by the lag time and the absorption rate-constant. This may also explain why deviation of the CsA concentrations predicted by the model was largest at 1 h post-dose. A recent study by Debord et al. [20] has applied a gamma distribution model of absorption to the pharmacokinetics of cyclosporin in stable renal transplant recipients and indeed yielded a better fit during the absorption phase compared with a classical exponential model with lag-time. An underestimation of the Cmax using the latter model, however, was not observed in the current population with scattering of the residuals around zero at both 1 and 2 h after ingestion of the oral dose. Compared with the FPIA, the RIA method slightly overestimated the data in our laboratory. Recently, Keown et al. [9] pointed out that differences between assays may have consequences for trough-level-based monitoring, while these assays performed quite comparable during the absorption phase of CsA in renal transplant recipients. The presented pharmacokinetic model, however, is flexible enough to deal with the apparent differences. The model therefore performs adequately despite the immunoassay used.
Calcineurin-inhibitor trough levels are a poor indicator of drug exposure and drug exposure should be quantitated by means of more accurate methods. To prevent structural underimmunosuppression or graft dysfunction due to chronic calcineurin-inhibitor toxicity, a compartmental population pharmacokinetic model with maximum a posteriori Bayesian fitting method offers several advantages. The presented model is reliable and more flexible than LSMs and does not depend on exact blood sampling time-points. In addition, it can be used in transplant recipients with gastrointestinal dysfunction and even offers the possibility to integrate pharmacodynamic parameters [21]. These practical and theoretical advantages argue in favour of this approach for future therapeutic drug monitoring.
Conflict of interest statement. None declared.
| Notes |
|---|
Correspondence and offprint requests to: S. C. L. M. Cremers, Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, PO Box 9600, NL-2300 RC Leiden, The Netherlands. Email: S.C.L.M.Cremers{at}lumc.nl
| References |
|---|
|
|
|---|
- Lindholm A, Kahan BD. Influence of cyclosporine pharmacokinetics, trough concentrations, and AUC monitoring on outcome after kidney transplantation. Clin Pharmacol Ther 1993; 54:205218[Web of Science][Medline]
- Kahan BD, Welsh M, Schoenberg L et al. Variable oral absorption of cyclosporine. A biopharmaceutical risk factor for chronic renal allograft rejection. Transplantation 1996; 62:599606[CrossRef][Web of Science][Medline]
- Nankivell BJ, Hibbins M, Chapman JR. Diagnostic utility of whole blood cyclosporine measurements in renal transplantation using triple therapy. Transplantation 1994; 58:989996[Medline]
- Oellerich M, Armstrong VW, Kahan B et al. Lake Louise Consensus Conference on cyclosporine monitoring in organ transplantation: report of the consensus panel. Ther Drug Monit 1995; 17:642654[Web of Science][Medline]
- Meyer MM, Munar M, Udeaja J, Bennett W. Efficacy of area under the curve cyclosporine monitoring in renal transplantation. J Am Soc Nephrol 1993; 4:13061315[Abstract]
- Amante AJ, Kahan BD. Abbreviated AUC strategy for monitoring cyclosporine microemulsion therapy in the immediate posttransplant period. Transplant Proc 1996; 28:21622163[Medline]
- Wacke R, Rohde B, Engel G et al. Comparison of several approaches of therapeutic drug monitoring of cyclosporin A based on individual pharmacokinetics. Eur J Clin Pharmacol 2000; 56:4348[Medline]
- Johnston A, David OJ, Cooney GF. Pharmacokinetic validation of Neoral absorption profiling. Transplant Proc 2000; 32:53S56S[CrossRef][Web of Science][Medline]
- Keown P, Cole E, Landsberg D et al. Absorption profiling of cyclosporine microemulsion (Neoral) during the first 2 weeks after renal transplantation. Transplantation 2001; 72:10241032[CrossRef][Medline]
- van der Pijl JW, Srivastava N, Denouel J et al. Pharmacokinetics of the conventional and microemulsion formulations of cyclosporine in pancreaskidney transplant recipients with gastroparesis. Transplantation 1996; 62:456462[Medline]
- Jelliffe RW, Schumitzky A, Bayard D et al. Model-based, goal-oriented, individualised drug therapy. Linkage of population modelling, new multiple model dosage design, Bayesian feedback and individualised target goals. Clin Pharmacokinet 1998; 34:5777[Medline]
- Proost JH. Adaptive control of drug dosage regimens using maximum a posteriori probability Bayesian fitting. Int J Clin Pharmacol Ther 1995; 33:531536[Medline]
- Rowland M, Tozer TN. Clinical Pharmacokinetics, Concepts and Applications. Lea & Febiger, Philadelphia, PA, 1989;459463
- Sheiner LB, Beal SL. Some suggestions for measuring predictive performance. J Pharmacokinet Biopharm 1981; 9:503512[CrossRef][Web of Science][Medline]
- Hariharan S, Johnson CP, Bresnahan BA, Taranto SE, McIntosh MJ, Stablein D. Improved graft survival after renal transplantation in the United States, 1988 to 1996. N Engl J Med 2000; 342:605612
[Abstract/Free Full Text] - Belitsky P, Levy GA, Johnston A. Neoral absorption profiling: an evolution in effectiveness. Transplant Proc 2000; 32:45S52S[Medline]
- Kovarik JM, Mueller EA, Richard F et al. Evidence for earlier stabilization of cyclosporine pharmacokinetics in de novo renal transplant patients receiving a microemulsion formulation. Transplantation 1996; 62:759763[CrossRef][Medline]
- Mahalati K, Belitsky P, Sketris I, West K, Panek R. Neoral monitoring by simplified sparse sampling area under the concentrationtime curve: its relationship to acute rejection and cyclosporine nephrotoxicity early after kidney transplantation. Transplantation 1999; 68:5562[CrossRef][Web of Science][Medline]
- Kahan BD, Welsh M, Urbauer DL et al. Low intraindividual variability of cyclosporin A exposure reduces chronic rejection incidence and health care costs. J Am Soc Nephrol 2000; 11:11221131
[Abstract/Free Full Text] - Debord J, Risco E, Harel M et al. Application of a gamma model of absorption to oral cyclosporine. Clin Pharmacokinet 2001; 40:375382[Medline]
- Halloran PF, Helms LM, Kung L, Noujaim J. The temporal profile of calcineurin inhibition by cyclosporine in vivo. Transplantation 1999; 68:13561361[CrossRef][Web of Science][Medline]
Accepted in revised form: 20.12.02
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
H. H. van Rossum, F. P.H.T.M. Romijn, K. J. Sellar, N. P.M. Smit, P. J.M. van der Boog, J. W. de Fijter, and J. van Pelt Variation in Leukocyte Subset Concentrations Affects Calcineurin Activity Measurement: Implications for Pharmacodynamic Monitoring Strategies Clin. Chem., March 1, 2008; 54(3): 517 - 524. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. M. Scholten, A. T. Rowshani, S. Cremers, F. J. Bemelman, M. Eikmans, E. van Kan, M. J. Mallat, S. Florquin, J. Surachno, I. J. ten Berge, et al. Untreated Rejection in 6-Month Protocol Biopsies Is Not Associated with Fibrosis in Serial Biopsies or with Loss of Graft Function J. Am. Soc. Nephrol., September 1, 2006; 17(9): 2622 - 2632. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. C. Roos-van Groningen, E. M. Scholten, P. M. Lelieveld, A. T. Rowshani, H. J. Baelde, I. M. Bajema, S. Florquin, F. J. Bemelman, E. de Heer, J. W. de Fijter, et al. Molecular Comparison of Calcineurin Inhibitor-Induced Fibrogenic Responses in Protocol Renal Transplant Biopsies J. Am. Soc. Nephrol., March 1, 2006; 17(3): 881 - 888. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||







