Prediction of brain clozapine and norclozapine concentrations in humans from a scaled pharmacokinetic model for rat brain and plasma pharmacokinetics
© Li et al.; licensee BioMed Central Ltd. 2014
Received: 14 April 2014
Accepted: 8 July 2014
Published: 20 August 2014
Clozapine is highly effective in treatment-resistant schizophrenia, although, there remains significant variability in the response to this drug. To better understand this variability, the objective of this study was to predict brain extracellular fluid (ECF) concentrations and receptor occupancy of clozapine and norclozapine in human central nervous system by translating plasma and brain ECF pharmacokinetic (PK) relationships in the rat and coupling these with known human disposition of clozapine in the plasma.
Unbound concentrations of clozapine and norclozapine were measured in rat brain ECF using quantitative microdialysis after subcutaneous administration of a 10 mg/kg single dose of clozapine or norclozapine. These data were linked with plasma concentrations obtained in the same rats to develop a plasma–brain ECF compartmental model. Parameters describing brain ECF disposition were then allometrically scaled and linked with published human plasma PK to predict human ECF concentrations. Subsequently, prediction of human receptor occupancy at several CNS receptors was based on an effect model that related the predicted ECF concentrations to published concentration-driven receptor occupancy parameters.
A one compartment model with first order absorption and elimination best described clozapine and norclozapine plasma concentrations in rats. A delay in the transfer of clozapine and norclozapine from plasma to the brain ECF compartment was captured using a transit compartment model approach. Human clozapine and norclozapine concentrations in brain ECF were simulated, and from these the median percentage of receptor occupancy of dopamine-2, serotonin-2A, muscarinic-1, alpha-1 adrenergic, alpha-2 adrenergic and histamine-1 for clozapine, and dopamine-2 for norclozapine were consistent with values reported in the literature.
A PK model that relates clozapine and norclozapine disposition in rat plasma and brain, including blood–brain barrier transport, was developed. Using allometry and published human plasma PK, the model was successfully translated to predict clozapine and norclozapine concentrations and accordant receptor occupancy of both agents in human brain. These predicted exposure and occupancy measures at several receptors that bind clozapine may be employed to extend our understanding of clozapine’s complex behavioral effects in humans.
KeywordsClozapine Norclozapine Receptor occupancy Population pharmacokinetics Translational modeling NONMEM
Schizophrenia is a debilitating disorder that affects approximately 1% of the global population without regard to race, sex or socioeconomic status . Its onset occurs typically in the late-teen years or early twenties and is characterized by a high rate of morbidity and mortality. Given these high personal and societal costs, investment in research aiming at understanding the biology of the disease, its genetic components and their interplay with environmental factors, continues on many levels. Over the past 50 years, pharmacotherapeutic support has been instrumental in managing primarily the positive symptoms of schizophrenia. It hinges on suppression of a central circuitry dysfunction that can be normalized by antagonism of dopamine D2 receptors in the striatum . Introduction of clozapine, the first so-called atypical antipsychotic approximately 25 years ago, represented a significant advance in our understanding of schizophrenia from a systems biology perspective in that this drug did not have the typical side effects of the first generation neuroleptics. This reduction in side effects was attributed to higher 5HT2A than D2 receptor binding .
However, clozapine pharmacology is not limited to D2 and 5HT2A antagonism. Albeit unintentionally, the drug binds to several other dopamine and serotonin receptor subtypes, muscarinic M1/M4 receptors, and alpha-1 adrenergic receptors with pharmacologically relevant affinity . From a clinical perspective, this broad receptor coverage may account for clozapine’s unique superiority in treatment resistant schizophrenia (TRS), even amongst other atypical antipsychotics. From a research perspective, the broad receptor coverage of clozapine conceivably makes the drug a useful tool to advance our understanding of complex pharmacotherapy that incorporates multiple interacting receptor systems.
The use of positron emission tomography (PET) imaging to measure receptor occupancy of clozapine and other atypical antipsychotics in humans has been invaluable in demonstrating the importance of D2 and 5HT2A receptor antagonism contributing to the efficacy of these drugs [5–7]. However, broader application of this non-invasive technique has been limited by the lack of ligands specific for other receptors to which clozapine has affinity. In this regard, availability of other approaches that are complementary to PET imaging would be useful. One possibility is to link non-clinical measurements of clozapine disposition in the brain with clinical studies of clozapine systemic exposure using a translational PK modeling approach. Prediction of clozapine CNS exposure could then be related to its receptor binding kinetics at multiple receptors to impart a virtual predicted pharmacodynamic component to a model. This approach has been used recently to predict CNS concentrations of atomoxetine and duloxetine that were in the range of receptor affinities associated with therapeutic doses . In a related manner, a population pharmacokinetic-pharmacodynamic (PK-PD) modeling approach was used to predict D2 receptor occupancy of olanzapine in humans , and the D2 and 5HT2A receptor occupancy of risperidone and its active metabolite paliperidone (9-OH risperidone) [5–7, 10]. These studies, as well as earlier PK-PD models applied to other CNS drugs [11, 12], provide confidence in the ability of this approach to deepen our understanding of drug action in human brain.
A recent study measured clozapine and its N-desmethyl metabolite, norclozapine, in extracellular fluid (ECF) of rat medial prefrontal cortex using quantitative microdialysis, and these results provided evidence of net efflux from brain across the blood–brain barrier (BBB) . This suggests that plasma concentrations may not be a good predictor of brain concentration for clozapine or norclozapine. Therefore, prediction of clozapine exposure in the ECF of human brain using a translational PK modeling approach could be cross-validated against PET results at D2 and 5HT2A receptor occupancy in humans, and subsequently used to estimate clozapine receptor occupancy at the drug’s other receptor targets for which PET tracers do not exist. Such comprehensive PK-PD model could potentially support individualized dosing of clozapine to improve its efficacy and CNS tolerability. It would also support research aimed at discovering new approaches for the treatment of schizophrenia in its different forms.
The purpose of this study was; (1) to build a PK model that accounted for both plasma and brain concentrations measured in rats; (2) to utilize this model to predict concentrations of clozapine and norclozapine in human brain. This would allow for the prediction of expected receptor occupancy in humans.
A single dose of clozapine (10 mg/kg) was administered subcutaneously to four male Wistar rats with an average weight of 0.35 kg purchased from Harlan (Zeist, The Netherlands). Three days prior to administration a microdialysis guide cannula was surgically implanted in the medial prefrontal cortex; at the same time, a catheter for blood sample collection was placed in the right jugular vein and was exteriorized through an incision at the top of the head. This vascular cannulation enabled an equivalent volume of saline replacement for each blood sample. A MetaQuant probe (6 mm, cellulose membrane, BrainLink, The Netherlands) was inserted into the guide cannula 24 hours prior to drug administration to enable sampling of brain extracellular fluid (ECF). Concentrations of clozapine and its N-desmethylated metabolite, norclozapine, were measured in plasma and brain ECF by HPLC with tandem mass spectrometry in the positive ion mode as previously described . For each rat, the unbound concentrations in each compartment were measured at 9 time points (0, 15, 30, 60, 90, 120, 240, 360 and 480 minutes) in plasma and 18 time points (-30, 0, 30, 60, 90, 120, 150, 180, 210, 240, 270, 300, 330, 360, 390, 420, 450 and 480 minutes) in brain. A single dose of norclozapine (10 mg/kg) was also administered subcutaneously to another five male Wistar rats with an average weight of 0.34 kg bought from the same Harlan laboratories. Concentrations were measured in plasma and brain ECF, and the same time points were used as those specified for clozapine.
Different model structures were initially evaluated using the system dynamics software VENSIM (Ventana Systems, Inc., MA, US). Thereafter a population approach was used to describe the pharmacokinetics of clozapine and norclozapine. Population PK parameters were estimated using a nonlinear mixed effect modeling approach, as implemented in NONMEM version 7.2 (Icon Development Solutions, Hanover, Maryland) using Wings for NONMEM version 7 . The first-order conditional estimation method (FOCE) with interaction was used to estimate the structural PK parameters and the random effects parameters.
Model development was started with an assessment of clozapine PK in plasma. One and two compartment models with first order absorption for clozapine in plasma were tested. A peripheral compartment structure was subsequently implemented to represent the brain extracellular fluid concentrations. The transfer characteristics of clozapine between the plasma and the brain compartment were evaluated using an intercompartmental clearance, CLin/CLout, as well as incorporating delay functions . These delay functions included a lag time and transit compartment approaches. Once the structural model for clozapine was established, the plasma compartment of norclozapine was integrated and then connected to the brain compartment. The same strategy was utilized in building the structural model for norclozapine concentrations that were measured following norclozapine administration. Clozapine and norclozapine concentration measurements were then combined from the 9 rats and modeled simultaneously in the final structural model. The volume of distribution of clozapine and norclozapine in brain were tested with and without fixing this parameter to a literature reported value . A parallel metabolic pathway from the extravascular space was also explored.
Between-animal variability (BAV) for PK parameters was assumed to be log-normally distributed and evaluated using an exponential model P i = P TV x e ηp where P i is the parameter estimate for the ith animal, and P TV is the typical parameter value at the population level. The variability between ith individual and population parameter values was described by η p , which was identically distributed with mean equal to 0 and variance, ω η 2. A combined additive and proportional model was first used to describe the intra-animal variability. If one of the elements of the model was found to be negligible and not significant, it was then removed from the residual error model. Residual error parameters were assumed to be normally distributed with mean equal to 0 and variance, σ2.
Model selection and evaluation
Model evaluation was based on a likelihood ratio test using the objective function value (OFV) provided by NONMEM. The minimum OFV returned by NONMEM is approximately equal to -2 × log likelihood (-2LL) and served as a guide during model design. A decrease in -2LL of 6.63 points for 1 degree of freedom was regarded as a significant model improvement, corresponding to a p value of 0.01 for nested models. The final model was further examined using goodness-of-fit plots generated using R version 2.13 based on the conditional weighted residuals distribution and the predicted versus observed concentrations at both the population and individual levels. Furthermore, the final pharmacokinetic model was also evaluated using a visual predictive check (VPC), and the uncertainty on each parameter was determined using a non-parametric bootstrap sampling with replacement 1000 times from the original dataset.
Human translation and expected receptor occupancy
After the pharmacokinetic model of clozapine in rat was finalized, the PK model framework was adapted by scaling PK parameters with allometric principles to predict human concentrations in brain. The following exponents were utilized scaling body weight to: clearance 0.75; volume of distribution 1; and first order rate constants 0.25 . A 50% of conversion from clozapine to norclozapine in humans was assumed in the model based on prior reports [19–21], and this was implemented in the simulated model assuming CLclo/F is equal to CLclo-p/F. Model performance was evaluated by comparing model simulated plasma concentrations to published human clozapine plasma concentrations  at steady state following 200, 300 and 400 mg daily doses. The published human clozapine data were reported as total concentrations, and these concentrations were converted to free concentration using 3% unbound fraction  prior to the comparison. After model validation, the simulated human clozapine and norclozapine concentrations were used to calculate the expected human receptor occupancy for the following receptors: dopamine 2 (D2); serotonin 2A (5-HT2A); muscarinic-1 (M1); alpha-1 adrenergic (α1); alpha-2 adrenergic (α2); and histamine-1 (H1) using published equilibrium dissociation constants (Kd) for clozapine [24–26] and norclozapine (D2 only) .
Rat population pharmacokinetics
Parameter estimates of final population pharmacokinetic model
Bootstrap analysis Median [5-95th percentiles]
Residual error (proportional)
The model presented represents a unique PK model developed from directly measured concentrations of clozapine and norclozapine in rat plasma and brain ECF. A multiple transit compartment model was used to account for a delay in the transport of clozapine and norclozapine from plasma to brain across the rat blood brain barrier. Some evidence suggests that Pgp may be involved in the process of clozapine transport [13, 28] across the blood brain barrier. The need to incorporate transit compartments in the present model is consistent with a Pgp role in clozapine transport across this barrier. Using an animal model, drug exposure can be measured by microdialysis at the target site. Based on a previously published non-compartmental analysis, the ratio of AUC between parent and metabolite in the rat indicated that only about 10% of parent drug was eliminated through metabolism . This is consistent with the ratio of norclozapine to clozapine clearance (Clclo-met = 9.91% of Clclo). As the results revealed, even with relatively rich sampling profiles, the uncertainty of some parameters, in particular of between-animal variability, was large likely because of the small number of animals in this study.
Simulated human plasma concentrations were based on previously published human plasma concentration data . The unbound plasma concentrations at steady state after a range of doses overlapped with published data corrected for the unbound fraction of clozapine (3%) . Subsequently, plasma exposures were linked to the plasma – brain structural PK parameters, using allometric scaling, that described clozapine and norclozapine transport between plasma and brain in the rat to ultimately predict human brain ECF exposure. This population pharmacokinetic approach, based on a transit compartmental approach as opposed to explicit assumption of a Pgp role and its associated interspecies scaling, enabled translational representation of the system across species to predict human brain ECF concentrations. As an atypical antipsychotic drug, clozapine targets D2 receptors as well as acts as an agonist or antagonist at several other receptors found in the CNS. In order to get a more complete profile of PK-PD linkage, percentage receptor occupancy of each receptor was calculated from 6 to 24 hour after three dosage levels. Our results show that the median percentage D2 receptor occupancy was in a range of 42% to 59% 6 hours after administration of a daily dose of 200–400 mg. This range is congruent with the 33% to 67% range reported by Nordström et al. , and agrees with the widely recognized understanding of low D2 receptor occupancy of therapeutic doses of clozapine relative to those obtained with therapeutic doses of other antipsychotics (first and second generation). In addition to D2, percent 5-HT2a receptor occupancy also overlapped with the results of Nordström et al.
The proposed PK model thus demonstrated the ability to extrapolate human systemic exposure to predict clozapine brain concentrations and associated receptor occupancy profiles in humans at clinically relevant doses. In addition, the model simultaneously captured parent and metabolite in the system, which is relevant since norclozapine also has activity at multiple receptors . However, the model can be improved in the precision of the PK parameter estimates by increasing the sample size. With this limitation taken into consideration, the model framework reported shows promise in predicting clozapine receptor occupancy at multiple receptors in human CNS, which can then be probed as a correlate to response and/or toxicity.
We demonstrated that a translational PK modeling approach was able to predict clozapine and norclozapine CNS exposures in humans, and these CNS exposures could then be related to their receptor binding kinetics at multiple receptors. Such a modeling approach could be foundational to the design of comprehensive PK-PD models and extend our understanding of clozapine’s complex behavioral effects in humans.
Central nervous system
Positron emission tomography
First order conditional estimation method
Objective function value
Treatment resistant schizophrenia.
Sources of funding
Brains On-Line, Inc. supported the rat microdialysis study live phase and sample analysis.
Eli Lilly and Company through the Indiana Clinical and Translational Sciences Institute (CTSI) supported the effort of Claire Li, Dr. Velez de Mendizabal and Dr. Bies.
- Mathers CD, Lopez AD, Murray CJL: The burden of disease and mortality by condition: data, methods and results for 2001. Global Burden of Disease and Risk Factors. Edited by: Lopez AD, Mathers CD, Ezzati M, Jamison DT, Murray CJL. 2006, New York: The World Bank and Oxford University PressGoogle Scholar
- Murray RM, Lappin J, Di FM: Schizophrenia: from developmental deviance to dopamine dysregulation. Eur Neuropsychopharmacol. 2008, 18 (Suppl 3): S129-S134.View ArticlePubMedGoogle Scholar
- Meltzer HY, Matsubara S, Lee JC: Classification of typical and atypical antipsychotic drugs on the basis of Dopamine D-1, D-2 and Serotonin2 Pki Values. J Pharmacol Exp Ther. 1989, 251 (1): 238-246.PubMedGoogle Scholar
- Horacek J, Bubenikova-Valesova V, Kopecek M, Palenicek T, Dockery C, Mohr P, Höschl C: Mechanism of action of atypical antipsychotic drugs and the neurobiology of schizophrenia. CNS Drugs. 2006, 20: 389-409.View ArticlePubMedGoogle Scholar
- Takeuchi H, Suzuki T, Bies RR, Remington G, Mamo DC, Pollock BG, Mimura M, Uchida H: Estimated dopamine D2 receptor occupancy from plasma concentrations of atypical antipsychotics and subjective experience/drug attitude in schizophrenia: an analysis of the CATIE data. Schizophr Res. 2013, 150 (2–3): 373-379.View ArticlePubMedGoogle Scholar
- Moriguchi S, Bies RR, Remington G, Suzuki T, Mamo DC, Watanabe K, Mimura M, Pollock BG, Uchida H: Estimated dopamine D2 receptor occupancy and remission in schizophrenia: analysis of the CATIE data. J Clin Psychopharmacol. 2013, 33 (5): 682-685.View ArticlePubMedGoogle Scholar
- Tsuboi T, Bies RR, Suzuki T, Mamo DC, Pollock BG, Graff-Guerrero A, Mimura M, Uchida H: Hyperprolactinemia and estimated dopamine D2 receptor occupancy in patients with schizophrenia: analysis of the CATIE data. Prog Neuropsychopharmacol Biol Psychiatry. 2013, 45: 178-182.View ArticlePubMedGoogle Scholar
- Kielbasa W, Stratford RE: Exploratory translational modeling approach in drug development to predict human brain pharmacokinetics and pharmacologically relevant clinical doses. Drug Metab Dispos. 2012, 40: 1-12.View ArticleGoogle Scholar
- Johnson M, Kozielska M, Pilla Reddy V, Vermeulen A, Li C, Grimwood S, de Greef R, Groothuis GM, Danhof M, Proost JH: Mechanism-based pharmacokinetic-pharmacodynamic modeling of the dopamine D2 receptor occupancy of olanzapine in rats. Pharm Res. 2011, 28: 2490-2504.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozlelska M, Johnson M, Pilla Reddy V, Vermeulen A, Li C, Grimwood S, de Greef R, Groothuis GM, Danhof M, Proost JH: Phamacokinetic-Pharmacodynamic modeling of the D2 and 5-HT(2A) receptor occupancy of risperidone and paliperidone in rats. Pharm Res. 2012, 29 (7): 1932-1948.View ArticleGoogle Scholar
- de Lange ECM, Ravenstijn PGM, Groenendaal D, van Steeg TS: Towards the prediction of CNS drug effect profiles in physiological and pathological conditions using microdialysis and mechanism-based pharmacokinetic-pharmacodynamic modeling. AAPS J. 2005, 7 (3): E532-E543.PubMed CentralView ArticlePubMedGoogle Scholar
- de Lange ECM: The mastermind approach to CNS drug therapy: Translational prediction of human brain distribution, target site kinetics and therapeutic effects. Fluids and Barriers of the CNS. 2013, 10: 12-28.PubMed CentralView ArticlePubMedGoogle Scholar
- Cremers TIFH, Flik G, Hofland C, Stratford RE: Microdialysis evaluation of clozapine and n-desmethylclozapine pharmacokinetics in rat brain. Drug Metab Dispos. 2012, 40: 1909-1916.View ArticlePubMedGoogle Scholar
- Holford NHG: Wings for NONMEM. [http://wfn.sourceforge.net/]
- Savic RM, Jonker DM, Kerbusch T, Karlsson MO: Implementation of a transit compartment model for describing drug absorption in pharmacokinetic studies. J Pharmacokinet Pharmacodyn. 2007, 34: 711-726.View ArticlePubMedGoogle Scholar
- Fridén M, Gupta A, Antonsson M, Bredberg U, Hammarlund-Udenaes M: In vitro methods for estimating unbound drug concentrations in brain interstitial and intracellular fluids. Drug Metab Dispos. 2007, 35: 1711-1719.View ArticlePubMedGoogle Scholar
- Feng Y, Pollock BG, Ferrell RE, Kimak MA, Reynold CF, Bies RR: Paroxetine: population pharmacokinetic analysis in late-life depression using sparse concentration sampling. Br J Clin Pharmacol. 2006, 61: 558-569.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharma V, McNeill JH: To scale or not to scale: the principles of dose extrapolation. Br J Clin Pharmacol. 2009, 157: 907-921.View ArticleGoogle Scholar
- Centorrino F, Baldessarrrini RJ, Kando JC, Frankenburg FR, Volpicelli SA, Flood JG: Clozapine and metabolites: concentrations in serum and clinical findings during treatment of chronically psychotic patients. J Clin Psychopharmacol. 1994, 14: 119-125.View ArticlePubMedGoogle Scholar
- Raedler TJ, Hinkelmann K, Wiedemann K: Variability of the in vitro metabolism of clozapine. Clin Neuropharmacol. 2008, 31: 347-352.View ArticlePubMedGoogle Scholar
- Couchman L, Morgan PE, Spencer EP, Flanagan RJ: Plasma clozapine, norclozapine, and the clozapine:norclozapine ratio in relation to prescribed dose and other factors: Data from a therapeutic drug monitoring service, 1993–2007. Ther Drug Monit. 2010, 32: 438-447.View ArticlePubMedGoogle Scholar
- Ismail Z, Wessels A, Ucbida H, Ng W, Mamo DC, Rajji TK, Pollock BG, Mulsant BH, Bies RR: Age and sex impact clozapine plasma concentrations in inpatients and outpatients with schizophrenia. Am J Geriatr Psychiatry. 2012, 20 (1): 53-60.View ArticlePubMedGoogle Scholar
- Gold Standard, Inc. Clozapine: Clinical Pharmacology [database online]. 2013, Available at: http://www.clinicalpharmacology.com. Accessed May 2, 2013Google Scholar
- Seeman PT: Atypical antipsychotics:mechanism of action. Can J Psychiatry. 2002, 47: 27-38.PubMedGoogle Scholar
- Bymaster FP, Calligaro DO, Falcone JF, Marsh RD, Moore NA, Tye NC, Seeman P, Wong DT: Radioreceptor binding profile of the atypical antipsychotic olanzapine. Neuropsychopharmacology. 1996, 14: 87-96.View ArticlePubMedGoogle Scholar
- Kroeze WK, Hufeisen SJ, Popadak BA, Renock SM, Steinberg S, Ernsberger P, Jayathilake K, Meltzer HY, Roth BL: H1-histamine receptor affinity predicts short-term weight gain for typical and atypical antipsychotic drugs. Neuropsychopharmacology. 2003, 28: 519-526.View ArticlePubMedGoogle Scholar
- Lidow MS: Neurotransmitter Receptors in Actions of Antipsychotic Medications. 2000, Boca Raton, FL: CRC PressView ArticleGoogle Scholar
- Doran A, Obach RS, Smith BJ, Hosea NA, Becker S, Callegari E, Chen C, Chen X, Choo E, Cuabfrigba J, Cox LM, Gibbs JP, Gibbs MA, Hatch H, Hop CE, Kasman IN, Laperle J, Liu J, Liu X, Logman M, Maclin D, Nedza FM, Nelson F, Olson E, Rahematpura S, Raunig D, Rogers S, Schmidt K, Spracklin DK, Szewc M:The impact of P-glycoprotein on the disposition of drugs targeted for indications of the central nervous system: evaluation using the MDR1A/1B knockout mouse model. Drug Metab Dispos. 2005, 33: 165-174.View ArticlePubMedGoogle Scholar
- Olsen CK, Brennum LT, Kreilgaard M: Using pharmacokinetic-pharmacodynamic modelling as a tool forprediction of therapeutic effective plasma levels of antipsychotics. Eur J Pharmacol. 2008, 584: 318-327.View ArticlePubMedGoogle Scholar
- Nordström AL, Farde L, Nyberg S, Karlsson P, Halldin C, Sedvall G: D1, D2, and 5–HT2 receptor occupancy in relation to clozapine serum concentration: a PET study of schizophrenia. Am J Psychiatry. 1995, 152: 1444-1449.View ArticlePubMedGoogle Scholar
- Bishara D, Taylor D: Upcoming agents for the treatment of schizophrenia: mechanism of action, efficacy and tolerability. Drugs. 2008, 68 (16): 2269-2292.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.