Modelling the progression of illicit substance use patterns from real-world evidence.

AIMS: To investigate an innovative pharmacometrics approach that addresses the challenges of using real-world evidence to model the progression of illicit substance use. METHODS: The modelling strategy analysed real-world data from the National Longitudinal Study of Adolescent to Adult Health (AddHealth) survey using survival analyses and differential equations. Respondents were categorized into drug-naïve, active users and nonusers. The transitions between categories were modelled using interval-censored parametric survival analysis. The resulting hazard rate functions were used as time-dependent rate constants in a differential equation system. Covariate models for sex and depression status were assessed. RESULTS: AddHealth enrolled 6504 American teenagers (median age 16 years, range 11-21 years); this cohort was followed with five interviews over a 22-year period; the median age at the last interview was 38 years (range 34-45 years). The percentages of illicit drug users at Interviews 1-5 were 7.7%, 5.9%, 15.8%, 21.4% and 0.98%, respectively. The generalized gamma distribution emerged as the preferred model for the survival functions for transitions between categories. Age-dependent prevalence was obtained from the differential equation system. Active drug use was more prevalent in males, increased in adolescence and college years, peaked at 24 years, and decreased to low levels by 35 years. Depression, which was more frequent in females, increased the drug-naïve-active user transition rates but not the active user-nonuser and nonuser-active user transition rates. The evidence did not support an interaction between sex and depression. CONCLUSIONS: The model provided a satisfactory approximation for the age-dependent progression of illicit substance use from preadolescence to early middle age.

Unraveling the Influence of Age, IQ, Education, and Negative Symptoms on Neurocognitive Performance in Schizophrenia: A Conditional Inference Tree Analysis.

INTRODUCTION: The complex nature of neurocognitive impairment in schizophrenia has been discussed in light of the mixed effects of antipsychotic drugs, psychotic symptoms, dopamine D2 receptor blockade, and intelligence quotient (IQ). These factors have not been thoroughly examined before. METHODS: This study conducted a comprehensive re-analysis of the CATIE data using machine learning techniques, in particular Conditional Inference Tree (CTREE) analysis, to investigate associations between neurocognitive functions and moderating factors such as estimated trough dopamine D2 receptor blockade with risperidone, olanzapine, or ziprasidone, Positive and Negative Syndrome Scale (PANSS), and baseline IQ in 573 patients with schizophrenia. RESULTS: The study reveals that IQ, age, and education consistently emerge as significant predictors across all neurocognitive domains. Furthermore, higher severity of PANSS-negative symptoms was associated with lower cognitive performance scores in several domains. CTREE analysis, in combination with a genetic algorithm approach, has been identified as particularly insightful for illustrating complex interactions between variables. Lower neurocognitive function was associated with factors such as age>52 years, IQ<94/95,<12/13 education years, and more pronounced negative symptoms (score<26). CONCLUSIONS: These findings emphasize the multifaceted nature of neurocognitive functioning in patients with schizophrenia, with the PANSS-negative score being an important predictor. This gives rise to a role in addressing negative symptoms as a therapeutic objective for enhancing cognitive impairments in these patients. Further research must examine nonlinear relationships among various moderating factors identified in this work, especially the role of D2 occupancy.

Evaluation of ChatGPT and Gemini large language models for pharmacometrics with NONMEM.

To assess ChatGPT 4.0 (ChatGPT) and Gemini Ultra 1.0 (Gemini) large language models on NONMEM coding tasks relevant to pharmacometrics and clinical pharmacology. ChatGPT and Gemini were assessed on tasks mimicking real-world applications of NONMEM. The tasks ranged from providing a curriculum for learning NONMEM, an overview of NONMEM code structure to generating code. Prompts in lay language to elicit NONMEM code for a linear pharmacokinetic (PK) model with oral administration and a more complex model with two parallel first-order absorption mechanisms were investigated. Reproducibility and the impact of "temperature" hyperparameter settings were assessed. The code was reviewed by two NONMEM experts. ChatGPT and Gemini provided NONMEM curriculum structures combining foundational knowledge with advanced concepts (e.g., covariate modeling and Bayesian approaches) and practical skills including NONMEM code structure and syntax. ChatGPT provided an informative summary of the NONMEM control stream structure and outlined the key NONMEM Translator (NM-TRAN) records needed. ChatGPT and Gemini were able to generate code blocks for the NONMEM control stream from the lay language prompts for the two coding tasks. The control streams contained focal structural and syntax errors that required revision before they could be executed without errors and warnings. The code output from ChatGPT and Gemini was not reproducible, and varying the temperature hyperparameter did not reduce the errors and omissions substantively. Large language models may be useful in pharmacometrics for efficiently generating an initial coding template for modeling projects. However, the output can contain errors and omissions that require correction.

Authors' response to letter to editor.

Authors' Response to Letter to Editor from Hinpetch Daungsupawong and Viroj Wiwanitkit.

Integration of DNA sequencing with population pharmacokinetics to improve the prediction of irinotecan exposure in cancer patients.

BACKGROUND: Irinotecan (CPT-11) is an anticancer agent widely used to treat adult solid tumours. Large interindividual variability in the clearance of irinotecan and SN-38, its active and toxic metabolite, results in highly unpredictable toxicity. METHODS: In 217 cancer patients treated with intravenous irinotecan single agent or in combination, germline DNA was used to interrogate the variation in 84 genes by next-generation sequencing. A stepwise analytical framework including a population pharmacokinetic model with SNP- and gene-based testing was used to identify demographic/clinical/genetic factors that influence the clearance of irinotecan and SN-38. RESULTS: Irinotecan clearance was influenced by rs4149057 in SLCO1B1, body surface area, and co-administration of 5-fluorouracil/leucovorin/bevacizumab. SN-38 clearance was influenced by rs887829 in UGT1A1, pre-treatment total bilirubin, and EGFR rare variant burden. Within each UGT1A1 genotype group, elevated pre-treatment total bilirubin and/or presence of at least one rare variant in EGFR resulted in significantly lower SN-38 clearance. The model reduced the interindividual variability in irinotecan clearance from 38 to 34% and SN-38 clearance from 49 to 32%. CONCLUSIONS: This new model significantly reduced the interindividual variability in the clearance of irinotecan and SN-38. New genetic factors of variability in clearance have been identified.

A Physiologically-Based Pharmacokinetic Model of the Brain Considering Regional Lipid Variance.

Modeling and simulation of the central nervous system provides a tool for understanding and predicting the distribution of small molecules throughout the brain tissue and cerebral spinal fluid (CSF), and these efforts often rely on empirical data to make predictions of distributions to move toward a better mechanistic understanding. A physiologically based pharmacokinetic model presented here incorporates multiple means of drug distribution to assemble a model for understanding potential factors that may determine the distribution of drugs across various regions of the brain, including both intra- and extracellular regions. Two classes of parameters are presented. The first concerns regional gross anatomic variability of the brain; the second concerns estimation of unbound fractions of drugs using know membrane phospholipid heterogeneity derived from regional lipid content. The model was then tested by comparing its outcomes to data from published human pharmacokinetic studies of acetaminophen, morphine, and phenytoin. The alignment of model predictions in the plasma, CSF, and tissue concentrations with the published data from studies of those three drugs suggests that the model can be a template for identifying drug localization in the brain. Clearly, knowledge of differentiated drug distribution in the brain is a requisite step in postulating pharmacodynamic and certain disease mechanisms. SIGNIFICANCE STATEMENT: This study concerns the application of heterogenous lipid distribution in brain tissue to predict regional variations in drug distribution in the brain via a mathematical model, thus expanding upon the current understanding of mechanisms of drug distribution in the central nervous system.

A population pharmacokinetic model based on HPTN 077 of long-acting injectable cabotegravir for HIV PrEP.

AIMS: Cabotegravir delivered as a long-acting intramuscular injection has shown superior efficacy to oral tenofovir-emtricitabine as pre-exposure prophylaxis (PrEP) for HIV. Cabotegravir pharmacokinetics (PK), like those of other long-acting depot preparations, exhibit variability between individuals and between injection occasions. The aim of this study is to describe the population pharmacokinetics of long-acting cabotegravir (CAB-LA). METHODS: Using available PK measurements from 133 participants in the HIV Prevention Trials Network (HPTN) 077 trial, we analysed CAB-LA PK data using nonlinear mixed-effects modelling to develop a population PK model. RESULTS: A two-compartment model with first order absorption best described the CAB-LA PK. The analysis identified between-occasion variability (BOV, i.e., differences in PK within one individual from one injection to the next) as a significant covariate affecting the absorption rate, with an estimated contribution of BOV to PK variability on the absorption rate (ka ) of 38.5%. Sex and body weight were identified as significant covariates influencing the absorption rate and apparent clearance of CAB-LA after intramuscular injection at various doses and frequencies. Male participants had 67% higher ka than female participants. Serially adding to the model body weight on clearance, sex on ka , and BOV on ka led to a decrease in the objective function value (OFV) of 24.4, 36 and 321.4, respectively. CONCLUSION: The public availability of this model will facilitate and enable a wide variety of future clinically relevant simulations to inform the optimal use of CAB-LA.

Population pharmacokinetic modeling of intramuscular and oral dexamethasone and betamethasone in Indian women.

Population analysis of pharmacokinetic data for five differing dosage forms and routes for dexamethasone and betamethasone in 48 healthy nonpregnant Indian women was performed that accounted for a partial and complex cross-over design. Single doses of 6 mg dexamethasone phosphate (DEX-P), betamethasone phosphate (BET-P), or 1:1 mixture of betamethasone phosphate and acetate (BET-PA) were administered orally (PO) or intramuscularly (IM). Plasma concentrations collected for two periods over 96 h were described with a two-compartment model with differing PO and IM first-order absorption inputs. Clearances and volumes were divided by the IM bioavailability [Formula: see text]. The homogeneous ages, body weights, and ethnicity of the women obviated covariate analysis. Parameter estimates were obtained by the Laplace estimation method implemented in NONMEM 7.4. Typical values for dexamethasone were clearance ([Formula: see text] of 9.29 L/h, steady-state volume ([Formula: see text] of 56.4 L, IM absorption constant [Formula: see text] of 0.460 1/h and oral absorption constant ([Formula: see text] of 0.936 1/h. Betamethasone parameters were CL/FIM of 5.95 L/h, [Formula: see text] of 72.4 L, [Formula: see text] of 0.971 1/h, and [Formula: see text] of 1.21 1/h. The PO to IM F values were close to 1.0 for both drugs. The terminal half-lives averaged about 7.5 h for DEX, 17 h for BET, and 78 h for BET from BET-PA with the latter reflecting very slow release of BET from the acetate ester. Overall, BET exhibited slower clearance, larger volume of distribution, faster absorption, and longer persistence than DEX. These data may be useful in considering exposures when substituting one form of corticosteroid for another.

Population pharmacodynamic modeling of intramuscular and oral dexamethasone and betamethasone effects on six biomarkers with circadian complexities in Indian women.

Population pharmacokinetic/pharmacodynamic (PK/PD) analysis was performed for extensive data for differing dosage forms and routes for dexamethasone (DEX) and betamethasone (BET) in 48 healthy nonpregnant Indian women in a partial and complex cross-over design. Single doses of 6 mg dexamethasone phosphate (DEX-P), betamethasone phosphate (BET-P), or 1:1 mixture of betamethasone phosphate and acetate (BET-PA) were administered orally (PO) or intramuscularly (IM) where each woman enrolled in a two-period cross-over study. Plasma concentrations collected over 96 h were described with a two-compartment model with differing PO and IM first-order absorption inputs. Overall, BET exhibited slower clearance, similar volume of distribution, faster absorption, and longer persistence than DEX with BET acetate producing extremely slow absorption but full bioavailability of BET. Six biomarkers were assessed over a 24-h baseline period with four showing circadian rhythms with complex baselines. These baselines and the strong responses seen after drug dosing were fitted with various indirect response models using the Laplace estimation methods in NONMEM 7.4. Both the PK and six biomarker responses were well-described with modest variability likely due to the homogeneous ages, weights, and ethnicities of the women. The drugs either inhibited or stimulated the influx processes with some models requiring joint inclusion of drug effects on circadian cortisol suppression. The biomarkers and order of sensitivity (lowest IC50/SC50 to highest) were: cortisol, T-helper cells, basophils, glucose, neutrophils, and T-cytotoxic cells. DEX sensitivities were generally greater than BET with corresponding mean ratios for these biomarkers of 2.86, 1.27, 1.72, 1.27, 2.69, and 1.06. Overall, the longer PK (e.g. half-life) of BET, but lesser PD activity (e.g. higher IC50), produces single-dose response profiles that appear quite similar, except for the extended effects from BET-PA. This comprehensive population modeling effort provides the first detailed comparison of the PK profiles and six biomarker responses of five commonly used dosage forms of DEX and BET in healthy women.

Model-Guided Antipsychotic Dose Reduction in Schizophrenia: A Pilot, Single-Blind Randomized Controlled Trial.

PURPOSE/BACKGROUND: Patients with schizophrenia as well as their psychiatrists are hesitant to reduce the antipsychotic dose in fear of relapse. To overcome such dilemmas, we developed models to individually calculate an oral dose that corresponds to a given target dopamine D2 receptor occupancy. METHODS/PROCEDURES: In this pilot, 52-week single-blind randomized controlled trial, 35 clinically stable patients with schizophrenia receiving either risperidone or olanzapine monotherapy were randomly assigned to dose reduction (n = 17) or dose maintenance group (n = 18). In the former group, baseline doses were reduced to the doses corresponding to 65% D2 occupancy (the lower end of therapeutic window) at trough that were calculated from randomly collected plasma concentrations using our models. FINDINGS/RESULTS: In the dose reduction group, doses of risperidone and olanzapine were decreased from 4.2 ± 1.9 to 1.4 ± 0.4 and 12.8 ± 3.9 to 6.7 ± 1.8 mg/d, whereas the doses in the dose maintenance group were 4.3 ± 1.9 and 15.8 ± 4.6 mg/d, respectively. Twelve subjects (70.5%) and 13 subjects (72.2%) in the dose reduction and dose maintenance groups completed the study (P = 0.604), whereas 3 subjects (18.8%) and none dropped out because of clinical worsening in the dose reduction and dose maintenance groups, respectively. There were not significant differences in score changes in Positive and Negative Syndrome Scale between the 2 groups but in Positive subscale scores in the Clinical Global Impression-Schizophrenia (0.4 ± 0.7 in the dose reduction group vs -0.1 ± 0.7 in the dose maintenance group, P = 0.029). IMPLICATIONS/CONCLUSIONS: Although our model-guided dose reduction strategy was found to be comparable with no-dose change in terms of dropout rates, safety issues have to be further examined.

Panos Macheras: a pioneering scientist in pharmaceutical science.

Professor Panos Macheras is a pioneering scientist in pharmacokinetics, pharmacodynamics and biopharmaceutics. His many important contributions to pharmaceutical science are reviewed.

Mathematical modeling of intracrine androgen metabolism in prostate cancer: Methodological aspects.

BACKGROUND: Progression of castration-recurrent/resistant prostate cancer (CRPC) relies in part on dihydrotestosterone derived from intratumoral androgen metabolism. Mathematical modeling provides a valuable tool for studies of androgen metabolism in CRPC. This modeling approach integrates existing knowledge about complex biologic systems and provides a means of interrogating the effects of various interventions. We sought to model a single reaction in the androgen biosynthesis network, namely the oxidation of androsterone (AND) to androstanedione (5α-dione) by four 3α-oxidoreductase enzymes, as an initial effort to establish the feasibility of our modeling approach. METHODS: Models were constructed for two cell culture systems, a non-prostate cancer cell line (CV-1) and a prostate cancer cell line (LAPC-4), using the SimBiology app (version 5.3) in MATLAB (version 8.6). The models included components for substrate (AND), product (5α-dione), each of the four enzymes, and each of the four enzyme-substrate complexes. Each enzymatic reaction consisted of a reversible enzyme-substrate binding step and an irreversible catalysis step. Rates of change for each component were described using ordinary differential equations. RESULTS: Mathematical models were developed with model parameter values derived from literature sources or from existing experimental data, which included gene expression measurements and substrate and product concentrations determined using liquid chromatography-tandem mass spectrometry. The models for both cell lines adequately described substrate and product concentrations observed after 12 h treatment with AND. CONCLUSIONS: This modeling approach represents an adaptable, extensible and mechanistic framework that reflects androgen metabolism. The models can be expanded systematically to describe the complex androgen metabolic pathways important for study of novel therapies for CRPC.

Population Pharmacokinetic Modeling To Estimate the Contributions of Genetic and Nongenetic Factors to Efavirenz Disposition.

Efavirenz pharmacokinetics is characterized by large between-subject variability, which determines both therapeutic response and adverse effects. Some of the variability in efavirenz pharmacokinetics has been attributed to genetic variability in cytochrome P450 genes that alter efavirenz metabolism, such as CYP2B6 and CYP2A6 While the effects of additional patient factors have been studied, such as sex, weight, and body mass index, the extent to which they contribute to variability in efavirenz exposure is inconsistently reported. The aim of this analysis was to develop a pharmacometric model to quantify the contribution of genetic and nongenetic factors to efavirenz pharmacokinetics. A population-based pharmacokinetic model was developed using 1,132 plasma efavirenz concentrations obtained from 73 HIV-seronegative volunteers administered a single oral dose of 600 mg efavirenz. A two-compartment structural model with absorption occurring by zero- and first-order processes described the data. Allometric scaling adequately described the relationship between fat-free mass and apparent oral clearance, as well as fat mass and apparent peripheral volume of distribution. Inclusion of fat-free mass and fat mass in the model mechanistically accounted for correlation between these disposition parameters and sex, weight, and body mass index. Apparent oral clearance of efavirenz was reduced by 25% and 51% in subjects predicted to have intermediate and slow CYP2B6 metabolizer status, respectively. The final pharmacokinetic model accounting for fat-free mass, fat mass, and CYP2B6 metabolizer status was consistent with known mechanisms of efavirenz disposition, efavirenz physiochemical properties, and pharmacokinetic theory. (This study has been registered at ClinicalTrials.gov under identifier NCT00668395.).

SSRI-antipsychotic combination in psychotic depression: sertraline pharmacokinetics in the presence of olanzapine, a brief report from the STOP-PD study.

OBJECTIVE: We recently reported an unexpected interaction between olanzapine and sertraline in a population being treated for psychotic depression. Contrary to knowledge of cytochrome p450 interactions sertraline increased apparent clearance of olanzapine by 30%. Here we examined the pharmacokinetics of sertraline in the same population. Existing studies suggest that sertraline apparent clearance is significantly increased in male subjects and suggested an age/sex interaction. METHODS: We studied subjects undergoing combination of sertraline/olanzapine treatment for psychotic depression in the Study of the Pharmacotherapy of Psychotic Depression. Nonlinear mixed effect modelling software was used to examine the sertraline pharmacokinetics, evaluating age, sex, race, and olanzapine exposure as covariates. RESULTS: Eighty-seven subjects (median age 62 years, 28 male subjects, 11 African-Americans) provided 138 samples for sertraline concentration. Olanzapine exposure had a 14.8-fold range. A one compartment model with combined residual error described the sertraline concentration data adequately. Half-life and sex effect on sertraline apparent clearance (males averaging 50% higher (p < 0.005); 96.6 l/h vs 64.8 in female subjects) were similar to previous reports. No other covariate (age, race or olanzapine exposure) had a significant impact on apparent clearance, and no age/sex interaction emerged. CONCLUSION: Sertraline pharmacokinetics were similar to historical descriptions in populations not taking antipsychotics. Unlike our unexpected finding that sertraline increases olanzapine apparent clearance, olanzapine exposure had no impact on sertraline pharmacokinetics. Copyright © 2016 John Wiley & Sons, Ltd.

A population pharmacokinetic model for R- and S-citalopram and desmethylcitalopram in Alzheimer's disease patients with agitation.

The citalopram for Alzheimer's disease trial evaluated citalopram for the management for agitation in Alzheimer's disease patients. Sparse data was available from this elderly patient population. A nonlinear mixed effects population pharmacokinetic modeling approach was used to describe the pharmacokinetics of R- and S-citalopram and their primary metabolite (desmethylcitalopram). A structural model with 4 compartments (one compartment/compound) with linear oral absorption and elimination described the data adequately. Overall, the model showed that clearance of the R-enantiomer was slower than the clearance of the S-enantiomer. Without accounting for any patient-specific covariates, the population estimate of the metabolic clearance of citalopram was 8.6 (R-citalopram) and 14 L/h (S-citalopram). The population estimate of the clearance of desmethylcitalopram was 23.8 (R-Dcit) and 38.5 L/h (S-Dcit). Several patient-specific covariates were found to have a significant effect on the pharmacokinetics of R,S-citalopram and desmethylcitalopram. A significant difference in the metabolic clearance of R-citalopram between males and females (13 vs 9.05 L/h) was identified in this analysis. Both R- and S-citalopram metabolic clearance decreased with age. Additionally, consistent with literature reports S-citalopram metabolic clearance increased with increasing body weight and was significantly influenced by CYPC19 genotype, with a difference of 5.8 L/h between extensive/rapid and intermediate/poor metabolizers. R,S-desmethylcitalopram clearance increased with increasing body weight. This model may allow for the opportunity to delineate the effect of R- and S-citalopram on pharmacodynamics outcomes related to the management of agitation in Alzheimer's disease.

Biomarkers for Diagnosis and Prognosis of Sinusoidal Obstruction Syndrome after Hematopoietic Cell Transplantation.

Reliable, noninvasive methods for diagnosing and prognosing sinusoidal obstruction syndrome (SOS) early after hematopoietic cell transplantation (HCT) are needed. We used a quantitative mass spectrometry-based proteomics approach to identify candidate biomarkers of SOS by comparing plasma pooled from 20 patients with and 20 patients without SOS. Of 494 proteins quantified, we selected 6 proteins (L-Ficolin, vascular cell adhesion molecule-1 [VCAM1], tissue inhibitor of metalloproteinase-1, von Willebrand factor, intercellular adhesion molecule-1, and CD97) based on a differential heavy/light isotope ratio of at least 2 fold, information from the literature, and immunoassay availability. Next, we evaluated the diagnostic potential of these 6 proteins and 5 selected from the literature (suppression of tumorigenicity-2 [ST2], angiopoietin-2 (ANG2), hyaluronic acid [HA], thrombomodulin, and plasminogen activator inhibitor-1) in samples from 80 patients. The results demonstrate that together ST2, ANG2, L-Ficolin, HA, and VCAM1 compose a biomarker panel for diagnosis of SOS. L-Ficolin, HA, and VCAM1 also stratified patients at risk for SOS as early as the day of HCT. Prognostic Bayesian modeling for SOS onset based on L-Ficolin, HA, and VCAM1 levels on the day of HCT and clinical characteristics showed >80% correct prognosis of SOS onset. These biomarkers may provide opportunities for preemptive intervention to minimize SOS incidence and/or severity.

A population PK model for citalopram and its major metabolite, N-desmethyl citalopram, in rats.

A population PK model was developed in order to simultaneously describe citalopram and its major metabolite, n-desmethyl citalopram, plasma concentrations in two different strain of rats after intravenous (IV) and oral (PO) administration of citalopram. Citalopram was administered to Sprague-Dawley (SD) rats at doses: 0.3, 1, 3, and 10 mg/kg IV and 10 mg/kg PO. The compound was dosed orally to Wistar rats at doses: 0.3, 1, 3, 10, 30 and 60 mg/kg. Plasma samples were collected for citalopram and metabolite. Pharmacokinetic analyses were conducted using NONMEM 7.2. Values below the quantification limit (BLQ < 0.1 ng/mL) were included in the analyses and treated as censored information. The disposition of citalopram was best described by a 3-compartment model and its desmethyl metabolite by a 2-compartment model. Several models for the absorption rate were explored (e.g. first, zero order and combined first and zero order absorption, Michaelis-Menten, lag time) in combination with dose and/or time dependent covariate effects. Dose dependent oral bioavailability properties were also identified in this analysis. Citalopram IV clearance and metabolite formation rate were adequately described as linear processes. Metabolite clearance was adequately described using a Michaelis-Menten clearance with different parameters depending on the strain. This analysis demonstrates a very complex absorption/metabolism process explaining the highly non-linear pharmacokinetics observed across all the doses. This is the first combined parent/metabolite population PK analysis in both SD and Wistar rats over a wide range of IV and PO dosages for citalopram, a compound that exhibits highly nonlinear oral pharmacokinetics in rats.

Pharmacokinetic modeling of tranexamic acid for patients undergoing cardiac surgery with normal renal function and model simulations for patients with renal impairment.

Tranexamic acid (TXA), an effective anti-fibrinolytic agent that is cleared by glomerular filtration, is used widely for cardiopulmonary bypass (CPB) surgery. However, an effective dosing regimen has not been fully developed in patients with renal impairment. The aims of this study were to characterize the inter-patient variability associated with pharmacokinetic parameters and to recommend a new dosing adjustment based on the BART dosing regimen for CPB patients with chronic renal dysfunction (CRD). Recently published data on CPB patients with normal renal function (n = 15) were re-examined with a two-compartment model using the ADAPT5 and NONMEMVII to identify covariates that explain inter-patient variability and to ascertain whether sampling strategies might affect parameter estimation. A series of simulations was performed to adjust the BART dosing regimen for CPB patients with renal impairment. Based on the two-compartmental model, the number of samples obtained after discontinuation of TXA infusion was found not to be critical in parameter estimation (p > 0.05). Both body weight and creatinine clearance were identified as significant covariates (p < 0.005). Simulations showed significantly higher than normal TXA concentrations in CRD patients who received the standard dosing regimen in the BART trial. Adjustment of the maintenance infusion rate based on the percent reduction in renal clearance resulted in predicted plasma TXA concentrations that were safe and therapeutic (~100 mg·L(-1) ). Our proposed dosing regimen, with consideration of renal function, is predicted to maintain effective target plasma concentrations below those associated with toxicity for patients with renal failure for CPB.

Prediction of brain clozapine and norclozapine concentrations in humans from a scaled pharmacokinetic model for rat brain and plasma pharmacokinetics.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.