Simulating realistic patient profiles from pharmacokinetic models by a machine learning postprocessing correction of residual variability.

We address the problem of model misspecification in population pharmacokinetics (PopPK), by modeling residual unexplained variability (RUV) by machine learning (ML) methods in a postprocessing step after conventional model building. The practical purpose of the method is the generation of realistic virtual patient profiles and the quantification of the extent of model misspecification, by introducing an appropriate metric, to be used as an additional diagnostic of model quality. The proposed methodology consists of the following steps: After developing a PopPK model, the individual residual errors IRES = DV-IPRED, are computed, where DV are the observations and IPRED the individual predictions and are modeled by ML to obtain IRESML. Correction of the IPREDs can then be carried out as DVML = IPRED + IRESML. The methodology was tested in a PK study of ropinirole, for which a PopPK model was developed while a second deliberately misspecified model was also considered. Various supervised ML algorithms were tested, among which Random Forest gave the best results. The ML model was able to correct individual predictions as inspected in diagnostic plots and most importantly it simulated realistic profiles that resembled the real data and canceled out the artifacts introduced by the elevated RUV, even in the case of the heavily misspecified model. Furthermore, a metric to quantify the extent of model misspecification was introduced based on the R2 between IRES and IRESML, following the rationale that the greater the extent of variability explained by the ML model, the higher the degree of model misspecification present in the original model.

Authors' response to letter to editor.

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

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.

Model-Based Prediction of Irinotecan-Induced Grade 4 Neutropenia in Cancer Patients: Influence of Incorporating Germline Genetic Factors in the Model.

Neutropenia is the major dose-limiting toxicity of irinotecan-based therapy. The objective of this study was to assess whether inclusion of germline genetic variants into a population pharmacokinetic/pharmacodynamic model can improve prediction of irinotecan-induced grade 4 neutropenia and identify novel variants of clinical value. A semimechanistic population pharmacokinetic/pharmacodynamic model was used to predict neutrophil response over time in 197 patients receiving irinotecan. Covariate analysis was performed for demographic/clinical factors and 4,781 genetic variants in 84 drug response- and toxicity-related genes to identify covariates associated with neutrophil response. We evaluated the predictive value of the model for grade 4 neutropenia reflecting different clinical scenarios of available data on identified demographic/clinical covariates, baseline and post-treatment absolute neutrophil counts (ANCs), individual pharmacokinetics, and germline genetic variation. Adding 8 genetic identified covariates (rs10929302 (UGT1A1), rs1042482 (DPYD), rs2859101 (HLA-DQB3), rs61754806 (NR3C1), rs9266271 (HLA-B), rs7294 (VKORC1), rs1051713 (ALOX5), and ABCB1 rare variant burden) to a model using only baseline ANCs improved prediction of irinotecan-induced grade 4 neutropenia from area under the receiver operating characteristic curve (AUC-ROC) of 50-64% (95% confidence interval (CI), 54-74%). Individual pharmacokinetics further improved the prediction to 74% (95% CI, 64-84%). When weekly ANC was available, the identified covariates and individual pharmacokinetics yielded no additional contribution to the prediction. The model including only ANCs at baseline and at week 1 achieved an AUC-ROC of 78% (95% CI, 69-88%). Germline DNA genetic variants may contribute to the prediction of irinotecan-induced grade 4 neutropenia when incorporated into a population pharmacokinetic/pharmacodynamic model. This approach is generalizable to drugs that induce neutropenia and ultimately allows for personalized intervention to enhance patient safety.

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.

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.

Deep Learning Methods Applied to Drug Concentration Prediction of Olanzapine.

Pharmacometrics and the utilization of population pharmacokinetics play an integral role in model-informed drug discovery and development (MIDD). Recently, there has been a growth in the application of deep learning approaches to aid in areas within MIDD. In this study, a deep learning model, LSTM-ANN, was developed to predict olanzapine drug concentrations from the CATIE study. A total of 1527 olanzapine drug concentrations from 523 individuals along with 11 patient-specific covariates were used in model development. The hyperparameters of the LSTM-ANN model were optimized through a Bayesian optimization algorithm. A population pharmacokinetic model using the NONMEM model was constructed as a reference to compare to the performance of the LSTM-ANN model. The RMSE of the LSTM-ANN model was 29.566 in the validation set, while the RMSE of the NONMEM model was 31.129. Permutation importance revealed that age, sex, and smoking were highly influential covariates in the LSTM-ANN model. The LSTM-ANN model showed potential in the application of drug concentration predictions as it was able to capture the relationships within a sparsely sampled pharmacokinetic dataset and perform comparably to the NONMEM model.

The MPRINT Hub Data, Model, Knowledge and Research Coordination Center: Bridging the gap in maternal-pediatric therapeutics research through data integration and pharmacometrics.

Maternal and pediatric populations have historically been considered "therapeutic orphans" due to their limited inclusion in clinical trials. Physiologic changes during pregnancy and lactation and growth and maturation of children alter pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. Precision therapy in these populations requires knowledge of these effects. Efforts to enhance maternal and pediatric participation in clinical studies have increased over the past few decades. However, studies supporting precision therapeutics in these populations are often small and, in isolation, may have limited impact. Integration of data from various studies, for example through physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling or bioinformatics approaches, can augment the value of data from these studies, and help identify gaps in understanding. To catalyze research in maternal and pediatric precision therapeutics, the Obstetric and Pediatric Pharmacology and Therapeutics Branch of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) established the Maternal and Pediatric Precision in Therapeutics (MPRINT) Hub. Herein, we provide an overview of the status of maternal-pediatric therapeutics research and introduce the Indiana University-Ohio State University MPRINT Hub Data, Model, Knowledge and Research Coordination Center (DMKRCC), which aims to facilitate research in maternal and pediatric precision therapeutics through the integration and assessment of existing knowledge, supporting pharmacometrics and clinical trials design, development of new real-world evidence resources, educational initiatives, and building collaborations among public and private partners, including other NICHD-funded networks. By fostering use of existing data and resources, the DMKRCC will identify critical gaps in knowledge and support efforts to overcome these gaps to enhance maternal-pediatric precision therapeutics.

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.

Model-Based Prediction of Irinotecan-Induced Grade 4 Neutropenia in Advanced Cancer Patients: Influence of Demographic and Clinical Factors.

Severe neutropenia is the major dose-liming toxicity of irinotecan-based chemotherapy. The objective was to assess to what extent a population pharmacokinetic/pharmacodynamic model including patient-specific demographic/clinical characteristics, individual pharmacokinetics, and absolute neutrophil counts (ANCs) can predict irinotecan-induced grade 4 neutropenia. A semimechanistic population pharmacokinetic/pharmacodynamic model was developed to describe neutrophil response over time in 197 patients with cancer receiving irinotecan. For covariate analysis, sex, race, age, pretreatment total bilirubin, and body surface area were evaluated to identify significant covariates on system-related parameters (mean transit time (MTT) and É£) and sensitivity to neutropenia effects of irinotecan and SN-38 (SLOPE). The model-based simulation was performed to assess the contribution of the identified covariates, individual pharmacokinetics, and baseline ANC alone or with incremental addition of weekly ANC up to 3 weeks on predicting irinotecan-induced grade 4 neutropenia. The time course of neutrophil response was described using the model assuming that irinotecan and SN-38 have toxic effects on bone marrow proliferating cells. Sex and pretreatment total bilirubin explained 10.5% of interindividual variability in MTT. No covariates were identified for SLOPE and γ. Incorporating sex and pretreatment total bilirubin (area under the receiver operating characteristic curve (AUC-ROC): 50%, 95% CI 50-50%) or with the addition of individual pharmacokinetics (AUC-ROC: 62%, 95% CI 53-71%) in the model did not result in accurate prediction of grade 4 neutropenia. However, incorporating ANC only at baseline and week 1 in the model achieved a good prediction (AUC-ROC: 78%, 95% CI 69-88%). These results demonstrate the potential applicability of a model-based approach to predict irinotecan-induced neutropenia, which ultimately allows for personalized intervention to maximize treatment outcomes.

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.

Application of machine learning to predict reduction in total PANSS score and enrich enrollment in schizophrenia clinical trials.

Clinical trial efficiency, defined as facilitating patient enrollment, and reducing the time to reach safety and efficacy decision points, is a critical driving factor for making improvements in therapeutic development. The present work evaluated a machine learning (ML) approach to improve phase II or proof-of-concept trials designed to address unmet medical needs in treating schizophrenia. Diagnostic data from the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) trial were used to develop a binary classification ML model predicting individual patient response as either "improvement," defined as greater than 20% reduction in total Positive and Negative Syndrome Scale (PANSS) score, or "no improvement," defined as an inadequate treatment response (<20% reduction in total PANSS). A random forest algorithm performed best relative to other tree-based approaches in model ability to classify patients after 6 months of treatment. Although model ability to identify true positives, a measure of model sensitivity, was poor (<0.2), its specificity, true negative rate, was high (0.948). A second model, adapted from the first, was subsequently applied as a proof-of-concept for the ML approach to supplement trial enrollment by identifying patients not expected to improve based on their baseline diagnostic scores. In three virtual trials applying this screening approach, the percentage of patients predicted to improve ranged from 46% to 48%, consistently approximately double the CATIE response rate of 22%. These results show the promising application of ML to improve clinical trial efficiency and, as such, ML models merit further consideration and development.

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.

Uridine Glucuronosyltransferase 2B7 Polymorphism-Based Pharmacogenetic Dosing of Epirubicin in FEC Chemotherapy for Early-Stage Breast Cancer.

BACKGROUND: Epirubicin is metabolized by uridine glucuronosyltransferase 2B7 (UGT2B7). Patients homozygous for the minor allele (CC) in the UGT2B7 -161 promoter polymorphism have lower clearance and significantly higher rates of leukopenia compared to wild-type homozygote (TT) or heterozygote (CT) patients. This study was designed to determine if TT and CT genotype patients could tolerate a higher epirubicin dose compared to CC genotype patients. PATIENTS AND METHODS: We studied women with histologically confirmed non-metastatic, invasive breast cancer who were scheduled to receive at least three cycles of FE100C in the (neo)adjuvant setting. Patients received standard-dose FE100C during the first 21-day cycle. Based on genotype, the epirubicin dose was escalated in the second and third cycles to 115 and 130 mg/m2 or to 120 and 140 mg/m2 for CT and TT genotype patients, respectively. The main outcome measurements were myelosuppression and dose-limiting toxicity. These were analyzed for relationships with the three genotypes. RESULTS: Forty-five patients were enrolled (10 CC, 21 CT, and 14 TT genotypes) and received 100 mg/m2 of epirubicin in the first cycle. Twelve and 10 TT patients were dose escalated at the second and third cycles, respectively; 16 CT patients were dose escalated at the second and third cycles. Leukopenia, but not febrile neutropenia, was genotype and dose dependent and increased in patients with CT and TT genotypes as their dose was increased. However, the third-cycle leukopenia rates were comparable to patients with the CC genotype receiving standard-dose epirubicin. CONCLUSION: Pharmacogenetically guided epirubicin dosing is well tolerated and allowed dose escalation without increased toxicity.

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.

Optimal Sampling Strategies for Irinotecan (CPT-11) and its Active Metabolite (SN-38) in Cancer Patients.

Irinotecan (CPT-11) is an anticancer agent widely used in the treatment of a variety of adult solid tumors. The objective of this study was to develop an optimal sampling strategy model that accurately estimates pharmacokinetic parameters of CPT-11 and its active metabolite, SN-38. This study included 221 patients with advanced solid tumors or lymphoma receiving CPT-11 single or combination therapy with 5-fluorouracil (5-FU)/leucovorin (LV) (FOLFIRI) plus bevacizumab from 4 separate clinical trials. Population pharmacokinetic analysis of CPT-11 and SN-38 was performed by non-linear mixed effects modeling. The optimal sampling strategy model was developed using D-optimality with expected distribution approach. The pharmacokinetic profiles of CPT-11 and SN-38 were best described by a 3- and 2-compartment model, respectively, with first-order elimination. Body surface area and co-administration with 5-FU/LV plus bevacizumab were significant covariates (p < 0.01) for volumes of the central compartment of CPT-11 and SN-38, and clearance of CPT-11. Pre-treatment total bilirubin and co-administration with 5-FU/LV and bevacizumab were significant covariates (p < 0.01) for clearance of SN-38. Accurate and precise predictive performance (r2 > 0.99, -2 < bias (%ME) < 0, precision (% RMSE) < 12) of both CPT-11 and SN-38 was achieved using: (i) 6 fixed sampling times collected at 1.5, 3.5, 4, 5.75, 22, 23.5 hours post-infusion; or (ii) 1 fixed time and 2 sampling windows collected at 1.5, [3-5.75], [22-23.5] hours post-infusion. The present study demonstrates that an optimal sampling design with three blood samples achieves accurate and precise pharmacokinetic parameter estimates for both CPT-11 and SN-38.

Circadian patterns of hallucinatory experiences in patients with schizophrenia: Potentials for chrono-pharmacology.

The objective of this study was to investigate possible circadian pattern of psychotic symptoms in patients with schizophrenia, which could be reflected on the dosing schedule/regimen, i.e. chrono-pharmacology. Patients with schizophrenia (ICD-10) who reported having auditory hallucination, receiving monotherapy with risperidone, olanzapine or paliperidone for at least two weeks were included. The subjects were provided a diary and asked to record the time and duration of auditory hallucinations during the eight time periods (i.e. 00:00-03:00, 03:00-06:00, 06:00-09:00, 09:00-12:00, 12:00-15:00, 15:00-18:00, 18:00-21:00, and 21:00-24:00). In the diary, times of medication doses and sleep were also recorded. Time and degree of peak and trough dopamine D2 receptor blockade with antipsychotics were estimated from 2 sparsely collected plasma drug concentrations. The prevalence and duration of auditory hallucinations were statistically examined among the eight time periods, respectively. Forty-nine patients participated in this study (mean ±â€¯SD age, 50.7 ±â€¯14.8 years; 36 men (73.5%); 34 inpatients (69.4%)). Auditory hallucinations occurred most frequently and lasted for the longest duration in the period of 18:00-21:00 (75.5% (37/49) and 1.37 ±â€¯1.67 h). This happened despite the fact that the difference in D2 receptor occupancy between the peak and trough was less than 2%, indicating a stable drug delivery. Since the dopamine D2 receptor blockade by antipsychotics was stable, the nocturnal circadian pattern found in this study may reflect intrinsic dopaminergic fluctuation or generally quieter environments at night. These circadian patterns may be considered to devise individualized treatment approach in the context of "chrono-pharmacology" for patients with schizophrenia.

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.