Xpose--an S-PLUS based population pharmacokinetic/pharmacodynamic model building aid for NONMEM

Daptomycin Population Pharmacokinetics in Patients Affected by Severe Gram-Positive Infections: An Update

Daptomycin pharmacokinetics may not depend on renal function only and it significantly differs between healthy volunteers and severely ill patients. Herein, we propose a population pharmacokinetics model based on 424 plasma daptomycin concentrations collected from 156 patients affected by severe Gram-positive infections during a routine therapeutic drug monitoring protocol. Model building and validation were performed using NONMEM 7.2 (ICON plc), Xpose4 and Perl-speaks-to-NONMEM. The final pop-PK model was a one-compartment first-order elimination model, with a 2.7% IIV for drug clearance (Cl), influence of creatinine clearance on drug clearance and of sex on distribution volume. After model validation, we simulated 10,000 patients with the Monte-Carlo method to predict the efficacy and tolerability of different daptomycin daily dosages. For the most common 6 mg/kg daily dose, the simulated probability of overcoming the toxic minimum concentration (24.3 mg/L) was 14.8% and the efficacy (expressed as a cumulative fraction of response) against methicillin-resistant S. aureus, S. pneumoniae and E. faecium was 95.77%, 99.99% and 68%, respectively. According to the model-informed precision dosing paradigm, pharmacokinetic models such as ours could help clinicians to perform patient-tailored antimicrobial dosing and maximize the odds of therapy success without neglecting toxicity risks.

Development of a population pharmacokinetic/pharmacodynamic model for various oral paclitaxel formulations co-administered with ritonavir and thrombospondin-1 based on data from early phase clinical studies

Purpose

Orally administered paclitaxel offers increased patient convenience while providing a method to prolong exposure without long continuous, or repeated, intravenous infusions. The oral bioavailability of paclitaxel is improved through co-administration with ritonavir and application of a suitable pharmaceutical formulation, which addresses the dissolution-limited absorption of paclitaxel. We aimed to characterize the pharmacokinetics of different paclitaxel formulations, co-administered with ritonavir, and to investigate a pharmacodynamic relationship between low-dose metronomic (LDM) treatment with oral paclitaxel and the anti-angiogenic marker thrombospondin-1 (TSP-1).

Methods

Fifty-eight patients treated with different oral paclitaxel formulations were included for pharmacokinetic analysis. Pharmacodynamic data was available for 36 patients. All population pharmacokinetic/pharmacodynamic modelling was performed using non-linear mixed-effects modelling.

Results

A pharmacokinetic model consisting of gut, liver, central, and peripheral compartments was developed for paclitaxel. The gastrointestinal absorption rate was modelled with a Weibull function. Relative gut bioavailabilities of the tablet and capsule formulations, as fractions of the gut bioavailability of the drinking solution, were estimated to be 0.97 (95%CI: 0.67–1.33) and 0.46 (95%CI: 0.34–0.61), respectively. The pharmacokinetic/pharmacodynamic relationship between paclitaxel and TSP-1 was modelled using a turnover model with paclitaxel plasma concentrations driving an increase in TSP-1 formation rate following an E_max relationship with an EC₅₀ of 284 ng/mL (95%CI: 122–724).

Conclusion

The developed pharmacokinetic model adequately described the paclitaxel plasma concentrations for the different oral formulations co-administered with ritonavir. This model, and the established pharmacokinetic/pharmacodynamic relationship with TSP-1, may facilitate future development of oral paclitaxel.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-022-04445-z.

Population Pharmacokinetic Modelling of Intravenous Immunoglobulin Treatment in Patients with Guillain-Barré Syndrome

Background and Objective

Intravenous immunoglobulin (IVIg) at a standard dosage is the treatment of choice for Guillain–Barré syndrome. The pharmacokinetics, however, is highly variable between patients, and a rapid clearance of IVIg is associated with poor recovery. We aimed to develop a model to predict the pharmacokinetics of a standard 5-day IVIg course (0.4 g/kg/day) in patients with Guillain–Barré syndrome.

Methods

Non-linear mixed-effects modelling software (NONMEM^®) was used to construct a pharmacokinetic model based on a model-building cohort of 177 patients with Guillain–Barré syndrome, with a total of 589 sequential serum samples tested for total immunoglobulin G (IgG) levels, and evaluated on an independent validation cohort that consisted of 177 patients with Guillain–Barré syndrome with 689 sequential serum samples.

Results

The final two-compartment model accurately described the daily increment in serum IgG levels during a standard IVIg course; the initial rapid fall and then a gradual decline to steady-state levels thereafter. The covariates that increased IgG clearance were a more severe disease (as indicated by the Guillain–Barré syndrome disability score) and concomitant methylprednisolone treatment. When the current dosing regimen was simulated, the percentage of patients who reached a target ∆IgG > 7.3 g/L at 2 weeks decreased from 74% in mildly affected patients to only 33% in the most severely affected and mechanically ventilated patients (Guillain–Barré syndrome disability score of 5).

Conclusions

This is the first population-pharmacokinetic model for standard IVIg treatment in Guillain–Barré syndrome. The model provides a new tool to predict the pharmacokinetics of alternative regimens of IVIg in Guillain–Barré syndrome to design future trials and personalise treatment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-022-01136-z.

Impact of Obesity on Postprandial Triglyceride Contribution to Glucose Homeostasis, Assessed with a Semimechanistic Model

The integrated glucose-insulin model is a semimechanistic model describing glucose and insulin after a glucose challenge. Similarly, a semiphysiologic model of the postprandial triglyceride (TG) response in chylomicrons and VLDL-V6 was recently published. We have developed the triglyceride-insulin-glucose-GLP-1 (TIGG) model by integrating these models and active GLP-1. The aim was to characterize, using the TIGG model, the postprandial response over 13 hours following a high-fat meal in 3 study populations based on body mass index categories: lean, obese, and very obese. Differential glucose and lipid regulation were observed between the lean population and obese or very obese populations. A population comparison revealed further that fasting glucose and insulin were elevated in obese and very obese when compared with lean; and euglycemia was achieved at different times postmeal between the obese and very obese populations. Postprandial insulin was incrementally elevated in the obese and very obese populations compared with lean. Postprandial chylomicrons TGs were similar across populations, whereas the postprandial TGs in VLDL-V6 were increased in the obese and very obese populations compared with lean. Postprandial active GLP-1 was diminished in the very obese population compared with lean or obese. The TIGG model described the response following a high-fat meal in individuals who are lean, obese, and very obese and provided insight into the possible regulation of glucose homeostasis in the extended period after the meal by utilizing lipids. The TIGG-model is the first model to integrate glucose and insulin regulation, incretin effect, and postprandial TGs response in chylomicrons and VLDL-V6.

Pharmacokinetic modelling of caspofungin to develop an extended dosing regimen in paediatric patients

BACKGROUND

Echinocandins are commonly used in treatment and prophylaxis of invasive fungal diseases. Intravenous daily dosing for prophylaxis in the outpatient setting can however become a hurdle for adequate compliance in the paediatric population.

OBJECTIVES

Simulations were performed to assess extended twice-weekly dosing for antifungal prophylaxis using caspofungin.

METHODS

A population pharmacokinetic model was developed based on previously published data from children aged 3 months to 17 years. Using the final model, Monte Carlo simulations were performed to assess the dose needed for adequate exposure in a twice-weekly setting. Mean weekly AUC0-24 h/MIC together with reported AUC0-24 h from previously reported paediatric trials were used to guide adequate exposure.

RESULTS AND CONCLUSIONS

A two-compartment model with linear elimination and allometric scaling using fixed exponents was found most adequate to describe the given paediatric populations. Simulations showed that a 200 mg/m2 twice-weekly regimen with maximal 200 mg total dose should result in exposures matching registered daily dosing as well as commonly used pharmacokinetic/pharmacodynamic targets.

Improving the tolerability of osimertinib by identifying its toxic limit

Background:

Osimertinib is the cornerstone in the treatment of epidermal growth factor receptor-mutated non-small cell lung cancer (NSCLC). Nonetheless, ±25% of patients experience severe treatment-related toxicities. Currently, it is impossible to identify patients at risk of severe toxicity beforehand. Therefore, we aimed to study the relationship between osimertinib exposure and severe toxicity and to identify a safe toxic limit for a preventive dose reduction.

Methods:

In this real-life prospective cohort study, patients with NSCLC treated with osimertinib were followed for severe toxicity (grade ⩾3 toxicity, dose reduction or discontinuation, hospital admission, or treatment termination). Blood for pharmacokinetic analyses was withdrawn during every out-patient visit. Primary endpoint was the correlation between osimertinib clearance (exposure) and severe toxicity. Secondary endpoint was the exposure–efficacy relationship, defined as progression-free survival (PFS) and overall survival (OS).

Results:

In total, 819 samples from 159 patients were included in the analysis. Multivariate competing risk analysis showed osimertinib clearance (c.q. exposure) to be significantly correlated with severe toxicity (hazard ratio 0.93, 95% CI: 0.88–0.99). An relative operating characteristic curve showed the optimal toxic limit to be 259 ng/mL osimertinib. A 50% dose reduction in the high-exposure group, that is 25.8% of the total cohort, would reduce the risk of severe toxicity by 53%. Osimertinib exposure was not associated with PFS nor OS.

Conclusion:

Osimertinib exposure is highly correlated with the occurrence of severe toxicity. To optimize tolerability, patients above the toxic limit concentration of 259 ng/mL could benefit from a preventive dose reduction, without fear for diminished effectiveness.

Swine as the Animal Model for Testing New Formulations of Anti-Inflammatory Drugs: Carprofen Pharmacokinetics and Bioavailability of the Intramuscular Route

Carprofen (CP) is a non-steroidal anti-inflammatory drug (NSAID) frequently used to treat respiratory diseases in numerous small animals, but also in large species. CP is a formidable candidate for further therapeutic research of human inflammatory diseases using the pig as an animal model. However, CP administration in swine is very uncommon and respective pharmacokinetics/bioavailability studies are scarce. A simultaneous population pharmacokinetic analysis after CP intravenous and intramuscular administrations in pigs has shown high extent and rate of absorption and a similar distribution profile with respect to man and other mammals. However, clearance and half-life values found in swine suggest a slower elimination process than that observed in man and some other animal species. Although not reported in other species, liver and kidney concentrations achieved at 48 h post-intramuscular administration in pigs were ten times lower than those found in plasma. Simulations pointed to 4 mg/kg every 24 h as the best dosage regimen to achieve similar therapeutic levels to those observed in other animal species. All these findings support the use of pig as an animal model to study the anti-inflammatory effects of CP in humans.

Pharmacokinetics and Dose Optimization Strategies of Para-Aminosalicylic Acid in Children with Rifampicin-Resistant Tuberculosis

Treatment options for children with Rifampicin-resistant tuberculosis (RR-TB) remain limited, and para-aminosalicylic acid (PAS) is still a relevant component of treatment regimens. Prevention of resistance to companion drugs by PAS is dose related, and at higher concentrations, PAS may exhibit significant bactericidal activity in addition to its bacteriostatic properties. The optimal dosing of PAS in children is uncertain, specifically for delayed-release granule preparations, which are the most used. A population pharmacokinetic model was developed describing PAS pharmacokinetics in children receiving routine RR-TB treatment. Model-based simulations evaluated current World Health Organization (WHO) weight-band doses against the adult pharmacokinetic target of 50 to 100 mg/liter for peak concentrations. Of 27 children included, the median (range) age and weight were 3.87 (0.58 to 13.7) years and 13.3 (7.15 to 30.5) kg, respectively; 4 (14.8%) were HIV positive. PAS followed one-compartment kinetics with first-order elimination and transit compartment absorption. The typical clearance in a 13-kg child was 9.79 liters/h. Increased PAS clearance was observed in both pharmacokinetic profiles from the only patient receiving efavirenz. No effect of renal function, sex, ethnicity, nutritional status, HIV status, antiretrovirals (lamivudine, abacavir, and lopinavir-ritonavir), or RR-TB drugs was detected. In simulations, target concentrations were achieved only using the higher WHO dose range of 300 mg/kg once daily. A transit compartment adequately describes absorption for the slow-release PAS formulation. Children should be dosed at the higher range of current WHO-recommended PAS doses and in a once-daily dose to optimize treatment.

A semimechanistic population pharmacokinetic and pharmacodynamic model incorporating autoinduction for the dose justification of TAS-114

TAS-114 is a dual deoxyuridine triphosphatase (dUTPase) and dihydropyrimidine dehydrogenase (DPD) inhibitor expected to widen the therapeutic index of capecitabine. Its maximum tolerated dose (MTD) was determined from a safety perspective in a combination study with capecitabine; however, its inhibitory effects on DPD activity were not assessed in the study. The dose justification to select its MTD as the recommended dose in terms of DPD inhibition has been required, but the autoinduction profile of TAS-114 made it difficult. To this end, an approach using a population pharmacokinetic (PPK)/pharmacodynamic (PD) model incorporating autoinduction was planned; however, the utility of this approach in the dose justification has not been reported. Thus, the aim of this study was to demonstrate the utility of a PPK/PD model incorporating autoinduction in the dose justification via a case study of TAS-114. Plasma concentrations of TAS-114 from 185 subjects and those of the endogenous DPD substrate uracil from 24 subjects were used. A two-compartment model with first-order absorption with lag time and an enzyme turnover model were selected for the pharmacokinetic (PK) model. Moreover, an indirect response model was selected for the PD model to capture the changes in plasma uracil concentrations. Model-based simulations provided the dose justification that DPD inhibition by TAS-114 reached a plateau level at the MTD, whereas exposures of TAS-114 increased dose dependently. Thus, the utility of a PPK/PD model incorporating autoinduction in the dose justification was demonstrated via this case study of TAS-114.

Population Pharmacokinetics of Intravenous Acyclovir in Oncologic Pediatric Patients

Background: Acyclovir represents the first-line prophylaxis and therapy for herpes virus infections. However, its pharmacokinetics in children exposes them to the risk of ineffective or toxic concentrations. The study was aimed at investigating the population pharmacokinetics (POP/PK) of intravenous (IV) acyclovir in oncologic children.

Methods: Patients (age, 8.6 ± 5.0 years, 73 males and 47 females) received IV acyclovir for prophylaxis (n = 94) and therapy (n = 26) under a therapeutic drug monitoring (i.e., minimum and maximal plasma concentrations, >0.5 and <25 mg/L, respectively). Plasma concentrations were fitted by nonlinear mixed effect modeling and a simulation of dosing regimens was performed. Findings were stratified according to an estimated glomerular filtration rate (eGFR) threshold of 250 ml/min/1.73 m².

Results: The final 1-compartment POP/PK model showed that eGFR had a significant effect on drug clearance, while allometric body weight influenced both clearance and volume of distribution. The population clearance (14.0 ± 5.5 L/h) was consistent across occasions. Simulation of standard 1-h IV infusion showed that a 10-mg/kg dose every 6 h achieved target concentrations in children with normal eGFR (i.e., ≤250 ml/min/1.73 m²). Increased eGFR values required higher doses that led to an augmented risk of toxic peak concentrations. On the contrary, simulated prolonged (i.e., 2 and 3-h) or continuous IV infusions at lower doses increased the probability of target attainment while reducing the risk of toxicities.

Conclusion: Due to the variable pharmacokinetics of acyclovir, standard dosing regimens may not be effective in some patients. Prospective trials should confirm the therapeutic advantage of prolonged and continuous IV infusions

Time-Dissociated Pharmacokinetic Pharmacodynamic Model of Cyclosporine Among Malaysian Renal Transplant Recipients

BACKGROUND

Cyclosporine is an essential component of many immunosuppressive regimens. However, its pharmacokinetic and pharmacodynamic (PKPD) modeling has not been widely investigated. This study aims to develop a time-dissociated PKPD model of cyclosporine in renal transplant patients.

METHODS

Medical records of renal transplant patients at Penang General Hospital were retrospectively analyzed. A time-dissociated PKPD model with covariate effects was developed using NONMEM to evaluate renal graft function response, quantified as estimated glomerular filtration rate (eGFR), toward the cyclosporine cumulative exposure (area under the concentration-time curve). The final model was integrated into a tool to predict the potential outcome. Individual eGFR predictions were evaluated based on the clinical response recorded as acute rejection/nephrotoxicity events.

RESULTS

A total of 1256 eGFR readings with 2473 drug concentrations were obtained from 107 renal transplant patients receiving cyclosporine. An Emax drug effect with a linear drug toxicity model best described the data. The baseline renal graft level (E0), maximum effect (Emax), area under the concentration-time curve achieving 50% of the maximum effect, and nephrotoxicity slope were estimated as 12.9 mL·min-1·1.73 m-2, 50.7 mL·min-1·1.73 m-2, 1740 ng·h·mL-1, and 0.00033, respectively. The hemoglobin level was identified as a significant covariate affecting the E0. The model discerned acute rejection from nephrotoxicity in 19/24 cases.

CONCLUSIONS

A time-dissociated PKPD model successfully described a large number of observations and was used to develop an online tool to predict renal graft response. This may help discern early rejection from nephrotoxicity, especially for patients unwilling to undergo a biopsy or those waiting for biopsy results.

Population Pharmacodynamic Modeling of Eflornithine-Based Treatments Against Late-Stage Gambiense Human African Trypanosomiasis and Efficacy Predictions of L-eflornithine-Based Therapy

Eflornithine is a recommended treatment against late-stage gambiense human African trypanosomiasis, a neglected tropical disease. Standard dosing of eflornithine consists of repeated intravenous infusions of a racemic mixture of L- and D-eflornithine. Data from three clinical studies, (i) eflornithine intravenous monotherapy, (ii) nifurtimox-eflornithine combination therapy, and (iii) eflornithine oral monotherapy, were pooled and analyzed using a time-to-event pharmacodynamic modeling approach, supported by in vitro activity data of the individual enantiomers. Our aim was to assess (i) the efficacy of the eflornithine regimens in a time-to-event analysis and (ii) the feasibility of an L-eflornithine-based therapy integrating clinical and preclinical data. A pharmacodynamic time-to-event model was used to estimate the total dose of eflornithine, associated with 50% reduction in baseline hazard, when administered as monotherapy or in the nifurtimox-eflornithine combination therapy. The estimated total doses were 159, 60 and 291 g for intravenous eflornithine monotherapy, nifurtimox-eflornithine combination therapy and oral eflornithine monotherapy, respectively. Simulations suggested that L-eflornithine achieves a higher predicted median survival, compared to when racemate is administered, as treatment against late-stage gambiense human African trypanosomiasis. Our findings showed that oral L-eflornithine-based monotherapy would not result in adequate efficacy, even at high dose, and warrants further investigations to assess the potential of oral L-eflornithine-based treatment in combination with other treatments such as nifurtimox. An all-oral eflornithine-based regimen would provide easier access to treatment and reduce burden on patients and healthcare systems in gambiense human African trypanosomiasis endemic areas. Graphical abstract.

Population Pharmacokinetic Modeling and Simulation of TQ-B3101 to Inform Dosing in Pediatric Patients With Solid Tumors

Background: TQ-B3101 is a novel kinase inhibitor currently in development for the treatment of advanced malignant solid tumor and relapsed or refractory ALK-positive anaplastic large cell lymphoma.

Methods: A population pharmacokinetic model was developed using data collected from a Phase 1 study and a Phase 2 study to characterize the pharmacokinetic of TQ-B3101 and its active metabolite (TQ-B3101M). The final model was used to optimize dosing of TQ-B3101 for pediatric patients (6-<18 years) with anaplastic large cell lymphoma.

Results: The pharmacokinetic of TQ-B3101 and TQ-B3101M was adequately described by a 1-compartment model with first-order absorption and elimination for parent drug coupled with a 2-compartment model with time-dependent clearance for the metabolite. The clearance of TQ-B3101M decreased over time with a maximum fractional reduction of 0.41. The estimated apparent clearance and apparent volume of distribution of TQ-B3101 were 2850 L/h and 4200 L, respectively. The elimination half-life of TQ-B3101 was 1.0 h. The distribution and elimination half-lives of TQ-B3101M at steady state were 4.9 and 39.4 h, respectively. The projected exposure of TQ-B3101M in virtual pediatric population following the body surface area tiered dosing regimen was similar to that in children pediatric patients after the recommended pediatric dose of crizotinib (280 mg/m2 twice daily), an analog of TQ-B3101M.

Conclusion: A population pharmacokinetic model was developed to provide optimal dose of regimen for further development of TQ-B3101 in pediatric patients with anaplastic large cell lymphoma.

Explorative Analysis of Treatment Outcomes of Levofloxacin- and Moxifloxacin-Based Regimens and Outcome Predictors in Ethiopian MDR-TB Patients: A Prospective Observational Cohort Study

Purpose/Background

Although Ethiopia is among the thirty high multi-drug resistant tuberculosis (MDR-TB) burden countries in the world, comparative therapeutic efficacy of moxifloxacin and levofloxacin has not been explored, particularly in MDR-TB patients. We therefore aimed to prospectively compare clinical outcomes and determine potential predictors of the outcomes among patients on moxifloxacin or levofloxacin-based MDR-TB drug regimens.

Methods

We analyzed clinical parameters and laboratory data of eighty MDR-TB patients on moxifloxacin- or levofloxacin-based regimens. The clinical outcomes were compared using the Kaplan–Meier survival functions and the outcome definitions of the 2013 World Health Organization. Monthly sputum culture conversions and a molecular line probe assay results were also assessed. Observed outcomes and patient-related variables between the two groups were compared using chi-square, Wilcoxon Rank and Fisher exact tests. We also determined the potential predictors influencing treatment outcomes of moxifloxacin and levofloxacin using Cox proportional hazard model.

Results

The levofloxacin-based treatment group had a lower failure rate and adverse drug events as well as better treatment success than the moxifloxacin-based group. Overall treatment success was 65%. Disaggregating the data revealed that 53.8% were cured, 11.2% completed treatment, 10.0% died, 11.2% failed, and 13.8% were lost-to-follow-up. The line probe assay result showed that 11.3% of the clinical isolates were resistant to fluoroquinolones and 3.8% were resistant to both fluoroquinolones and injectable anti-TB agents. Treatment regimen type, culture conversion rate, alcohol use, cavity lesion, serum levels of creatinine and alanine aminotransferase were independent predictors of treatment outcome.

Conclusion

The levofloxacin-based regimen group has a better overall treatment success than the moxifloxacin-based group among MDR-TB patients. Clinical parameters and substance use history of the patients influenced treatment outcomes. We recommend further broader clinical studies to substantiate our findings as an input to review MDR-TB treatment guidelines.

Characterizing Absorption Properties of Dispersible Pretomanid Tablets Using Population Pharmacokinetic Modelling

BACKGROUND AND INTRODUCTION

The dispersible tablet formulation (DTF) of pretomanid has been developed to facilitate future use in children. This work aimed to assess the pharmacokinetics (PK) and relative bioavailability of the DTF compared to the marketed formulation (MF) and the potential influence of dose.

METHODS

Pretomanid DTF was investigated in a single-dose, randomized, four-period, cross-over study, with 7 days of washout between doses. Forty-eight healthy volunteers were enrolled and randomized into one of two panels to receive doses either in the fasted state or after a high-fat meal. Each volunteer received doses of 10, 50, and 200 mg DTF, and 200 mg MF pretomanid. Blood samples for pharmacokinetic assessment were drawn following a rich schedule up to 96 h after each single dose. The study data from the panel receiving the high-fat meal were analyzed using a nonlinear mixed-effects modeling approach, and all data were characterized with noncompartmental methods.

RESULTS

A one-compartment model with first-order elimination and absorption through a transit compartment captured the mean and variability of the observed pretomanid concentrations with acceptable precision. No significant difference in bioavailability was found between formulations. The mean absorption time for the DTF was typically 137% (86-171%) of that for the MF. The bioavailability was found to be dose dependent with a small positive and larger negative correlation under fed and fasted conditions, respectively.

CONCLUSION

Using data from a relative bioavailability study in healthy adult volunteers, a mathematical model has been developed to inform dose selection for the investigation of pretomanid in children using the new dispersible tablet formulation. Under fed conditions and at the currently marketed adult dose of 200 mg, the formulation type was found to influence the absorption rate, but not the bioavailability. The bioavailability of the DTF was slightly positively correlated with doses when administered with food.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT04309656, first posted on 16 March 2020.

Population Pharmacokinetic Properties of Omecamtiv Mecarbil in Healthy Subjects and Patients with Heart Failure with Reduced Ejection Fraction

ABSTRACT

Omecamtiv mecarbil is a small molecule that has been shown to improve cardiac function in patients with heart failure with reduced ejection fraction and is currently being developed as an oral modified release (MR) tablet for subjects with chronic HF. The objectives of this study were to analyze the pharmacokinetic (PK) properties of omecamtiv mecarbil and to investigate the effects of potential covariates on pertinent PK parameters using population PK modeling of data from 3 clinical trials in healthy subjects (N=85) and 3 clinical trials in patients with heart failure (N=4261). The population PK analysis was performed using a non-linear mixed effects modeling approach. Omecamtiv mecarbil has a clearance of 11.7 L/hr (0.701 % RSE [relative standard error]) and a central volume of distribution of 275 L (2.12% RSE). The estimated half-life of omecamtiv mecarbil was 33 hours. Body weight and estimated glomerular filtration rate (eGFR) were significant covariates, but their effect on exposure was modest and lacked clinical relevance. Additional covariates including sex, race, bilirubin, albumin, concomitant medications, New York Heart Association Functional Classification, N-terminal-pro hormone B-type natriuretic peptide (NT-proBNP), troponin I, creatine kinase MB, serum hemoglobin, tablet formulation, aspartate aminotransferase, and serum urea was tested and found to have no impact on omecamtiv mecarbil exposures. The results of this integrated evaluation of the impact of covariates on the systemic exposure of omecamtiv mecarbil suggest dose adjustment is not required for the studied subpopulations of patients with heart failure.

Imeglimin population pharmacokinetics and dose adjustment predictions for renal impairment in Japanese and Western patients with type 2 diabetes

Imeglimin is an orally administered first-in-class drug to treat type 2 diabetes mellitus (T2DM) and is mainly excreted unchanged by the kidneys. The present study aimed to define the pharmacokinetic (PK) characteristics of imeglimin using population PK analysis and to determine the optimal dosing regimen for Japanese patients with T2DM and chronic kidney disease (CKD). Imeglimin plasma concentrations in Japanese and Western healthy volunteers, and patients with T2DM, including patients with mild to severe CKD with an estimated glomerular filtration rate (eGFR) > 14 mL/min/1.73 m² were included in a population PK analysis. Pharmacokinetic simulations were conducted using a population PK model, and the area under concentration-time curve (AUC) was extrapolated with power regression analysis to lower eGFR. The influence of eGFR, weight, and age on apparent clearance and of dose on relative bioavailability were quantified by population PK analysis. Simulations and extrapolation revealed that the recommended dosing regimen based on the AUC was 500 mg twice daily (b.i.d.) for patients with eGFR 15-45 mL/min/1.73 m² , and 500 mg with a longer dosing interval was suggested for those with eGFR < 15. Simulations revealed that differences in plasma AUCs between Japanese and Western patients at the same dose were mainly driven by a difference in the eGFR and that the plasma AUC after 1,000 and 1,500 mg b.i.d. in Japanese and Western patients, respectively, was comparable in the phase IIb studies. These results indicate suitable dosages of imeglimin in the clinical setting of T2DM with renal impairment.

Pharmacology of enflicoxib, a new coxib drug: Efficacy and dose determination by clinical and pharmacokinetic-guided approach for the treatment of osteoarthritis in dogs based on an acute arthritis induction model

Enflicoxib is a newly developed NSAID of the coxib class. The optimal therapeutic dose to be confirmed in the field studies was established using a combination of pharmacokinetic (PK) modelling and pharmacodynamic (PD) studies. First, a PK study was performed to determine the plasmatic profile of enflicoxib and its active pyrazol metabolite in dogs. Thereafter, two studies using a urate crystal‐induced acute arthritis model allowed to correlate efficacy with plasmatic concentrations. Finally, a population PK model was developed to establish the Minimum Effective Concentration (MEC) and the Maximum Tolerated Concentration (MTC). Enflicoxib plasma concentrations were highest for the first 48 h. Thereafter, pyrazol metabolite concentrations were higher and persisted up to the end of the study. No reduction on the lameness (CLS) or pain scores (PS) was observed in the first hours after enflicoxib administration so no MEC could be established for the parent compound. Both CLS and PS were greatly reduced when the pyrazol metabolite achieved concentrations of 411 ng/ml or higher, so this concentration was established as the MEC for the pyrazol metabolite. Enflicoxib MTC was established at 6723 ng/ml whereas for the pyrazol metabolite it was 4258 ng/ml. The population PK model showed that a loading enflicoxib dose of 8 mg/kg followed by weekly maintenance doses of 4 mg/kg would achieve stable concentrations of the pyrazol metabolite within the therapeutic window (between the MEC and the MTC), and it was considered the most adequate posology to be confirmed in the field clinical studies for the treatment of canine osteoarthritis.

Pharmacogenetics-based population pharmacokinetic analysis for dose optimization of ritonavir-boosted atazanavir in Thai adult HIV-infected patients

BACKGROUND

This population pharmacokinetic-pharmacogenetic study aimed to investigate the optimal dose of RTV-boosted ATV (ATV/RTV) for Thai adult HIV-infected patients.

METHODS

A total of 1460 concentrations of ATV and RTV from 544 patients receiving an ATV/RTV-based regimen were analyzed. The CYP3A5 6986 A > G, ABCB1 3435 C > T, ABCB1 2677 G > T, SLCO1B1 521 T > C, and NR1I2 63396 C > T were genotyped. A population pharmacokinetic model was performed using a nonlinear mixed-effect model (NONMEM®). Monte Carlo simulations were conducted to compare the percentages of patients achieving the therapeutic range of ATV through concentrations (C_trough).

RESULTS

The apparent oral clearance of ATV (CL/F_ATV) without RTV was 7.69 L/h with interindividual variability (IIV) of 28.7%. Patients with CYP3A5 6986 GG had a 7.1% lower CL/F_ATV than those with AA or AG genotype. The CL/F_ATV decreased by 10.8% for females compared with males. Simulation results showed higher percentages (~70%) of patient receiving doses of 200/100 or 200/50 mg achieved the target ATV C_trough, while more patients (~40%) receiving a standard dose (300/100 mg) had ATV C_trough above this target.

CONCLUSIONS

Both CYP3A5 6986 A > G and female decreased CL/F_ATV in Thai HIV-infected patients. Simulations supported that the reduced dose of ATV/RTV was sufficient to achieve the target concentration for Thai population.

Dosing strategies of imipenem in neonates based on pharmacometric modelling and simulation

OBJECTIVES

Imipenem is a broad-spectrum antibacterial agent used in critically ill neonates after failure of first-line treatments. Few studies have described imipenem disposition in this population. The objectives of our study were: (i) to characterize imipenem population pharmacokinetics (PK) in a cohort of neonates; and (ii) to conduct model-based simulations to evaluate the performance of six different dosing regimens aiming at optimizing PK target attainment.

METHODS

A total of 173 plasma samples from 82 neonates were collected over 15 years at the Lausanne University Hospital, Switzerland. The majority of study subjects were preterm neonates with a median gestational age (GA) of 27 weeks (range: 24-41), a postnatal age (PNA) of 21 days (2-153) and a body weight (BW) of 1.16 kg (0.5-4.1). PK data were analysed using non-linear mixed-effect modelling (NONMEM).

RESULTS

A one-compartment model best characterized imipenem disposition. Population PK parameters estimates of CL and volume of distribution were 0.21 L/h and 0.73 L, with an interpatient variability (CV%) of 20.1% on CL in a representative neonate (GA 27 weeks, PNA 21 days, BW 1.16 kg, serum creatinine, SCr 46.6 μmol/L). GA and PNA exhibited the greatest impact on PK parameters, followed by SCr. These covariates explained 36% and 15% of interindividual variability in CL, respectively.Simulated regimens using a dose of 20-25 mg/kg every 6-12 h according to postnatal age led to the highest PTA (T>MIC over 100% of time).

CONCLUSIONS

Dosing adjustment according to BW, GA and PNA optimizes imipenem exposure in neonates.

Pharmacokinetic-pharmacodynamic modelling of the anti-FcRn monoclonal antibody rozanolixizumab: Translation from preclinical stages to the clinic

Rozanolixizumab is a fully humanized high‐affinity anti‐human neonatal Fc receptor (FcRn) monoclonal antibody (mAb) that accelerates the removal of circulating immunoglobulin G (IgG), including pathogenic IgG autoantibodies, via the natural lysosomal degradation pathway. The aim of this study was to develop a pharmacokinetic/pharmacodynamic (PK/PD) model characterizing the effect of rozanolixizumab on IgG levels in cynomolgus monkeys, translate it into humans to support the first‐in‐human (FIH) rozanolixizumab clinical trial study design, and, ultimately, develop a PK/PD model in humans. Simulations from the preclinical model were performed to predict IgG responses in humans and select clinically relevant doses in the FIH study. Good alignment was observed between predicted and observed reductions in IgG, which increased with increasing dose in the FIH study. The model successfully described the PK of the 4 and 7 mg/kg intravenous (i.v.) dose groups, although the PKs were underpredicted for the 1 mg/kg i.v. dose group. Updating the model with subsequent human data identified parameters that deviated from preclinical assumptions. The updated PK/PD model was able to effectively characterize the PK FcRn‐IgG nonlinear system in response to rozanolixizumab in the FIH data.

Integrated exposure-response analysis of efficacy and safety of lurbinectedin to support the dose regimen in small-cell lung cancer

Purpose

These exposure–response (E–R) analyses integrated lurbinectedin effects on key efficacy and safety variables in relapsed SCLC to determine the adequacy of the dose regimen of 3.2 mg/m² 1-h intravenous infusion every 3 weeks (q3wk).

Methods

Logistic models and Cox regression analyses were applied to correlate lurbinectedin exposure metrics (AUC_tot and AUC_u) with efficacy and safety endpoints: objective response rate (ORR) and overall survival (OS) in SCLC patients (n = 99) treated in study B-005 with 3.2 mg/m² q3wk, and incidence of grade 4 (G4) neutropenia and grade 3–4 (G ≥ 3) thrombocytopenia in a pool of cancer patients from single-agent phase I to III studies (n = 692) treated at a wide range of doses. A clinical utility index was used to assess the appropriateness of the selected dose.

Results

Effect of lurbinectedin AUC_u on ORR best fitted to a sigmoid-maximal response (E_max) logistic model, where E_max was dependent on chemotherapy-free interval (CTFI). Cox regression analysis with OS found relationships with both CTFI and AUC_u. An E_max logistic model for G4 neutropenia and a linear logistic model for G ≥ 3 thrombocytopenia, which retained platelets and albumin at baseline and body surface area, best fitted to AUC_tot and AUC_u. AUC_u between approximately 1000 and 1700 ng·h/L provided the best benefit/risk ratio, and the dose of 3.2 mg/m² provided median AUC_u of 1400 ng·h/L, thus maximizing the proportion of patients within that lurbinectedin target exposure range.

Conclusions

The relationships evidenced in this integrated E–R analysis support a favorable benefit-risk profile for lurbinectedin 3.2 mg/m² q3wk.

Trial registration

Clinicaltrials.gov: NCT02454972; registered May 27, 2015.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-021-04366-3.

Population Pharmacokinetic Modeling of Inebilizumab in Subjects with Neuromyelitis Optica Spectrum Disorders, Systemic Sclerosis, or Relapsing Multiple Sclerosis

Background and Objective

Inebilizumab is a humanized, affinity-optimized, afucosylated immunoglobulin (Ig)-G1κ monoclonal antibody that binds to CD19, resulting in effective depletion of peripheral B cells. It is being developed to treat various autoimmune diseases, including neuromyelitis optica spectrum disorders (NMOSD), systemic sclerosis (SSc), and relapsing multiple sclerosis (MS).

Methods

Pharmacokinetic data from a pivotal study in adult subjects with NMOSD and two early-stage studies in subjects with SSc or relapsing MS were pooled and simultaneously analyzed using a population approach.

Results

Upon intravenous administration, the pharmacokinetics of inebilizumab were adequately described by a two-compartment model with parallel first-order and time-dependent nonlinear elimination pathways. An asymptotic nonlinear elimination suggests that inebilizumab undergoes receptor (CD19)-mediated clearance. The estimated systemic clearance (CL) of the first-order elimination pathway (0.188 L/day) and the volume of distribution (V_d) (5.52 L) were typical for therapeutic immunoglobulins. The elimination half-life was approximately 18 days. The maximum velocity (V_max) of the nonlinear elimination pathway decreased with time, presumably due to the depletion of B cells upon inebilizumab administration. As for other therapeutic monoclonal antibodies, the CL and V_d of inebilizumab increased with body weight.

Conclusions

The presence of antidrug antibodies, status of hepatic or renal function, and use of small-molecule drugs commonly used by subjects with NMOSD had no clinically relevant impact on the pharmacokinetics of inebilizumab. The nonlinear elimination pathway at the 300 mg therapeutic dose level is not considered clinically relevant.

Semimechanistic models to relate noxious stimulation, movement, and pupillary dilation responses in the presence of opioids

Intraoperative targeting of the analgesic effect still lacks an optimal solution. Opioids are currently the main drug used to achieve antinociception, and although underdosing can lead to an increased stress response, overdose can also lead to undesirable adverse effects. To better understand how to achieve the optimal analgesic effect of opioids, we studied the influence of remifentanil on the pupillary reflex dilation (PRD) and its relationship with the reflex movement response to a standardized noxious stimulus. The main objective was to generate population pharmacodynamic models relating remifentanil predicted concentrations to movement and to pupillary dilation during general anesthesia. A total of 78 patients undergoing gynecological surgery under general anesthesia were recruited for the study. PRD and movement response to a tetanic stimulus were measured multiple times before and after surgery. We used nonlinear mixed effects modeling to generate a population pharmacodynamic model to describe both the time profiles of PRD and movement responses to noxious stimulation. Our model demonstrated that movement and PRD are equally depressed by remifentanil. Using the developed model, we changed the intensity of stimulation and simulated remifentanil predicted concentrations maximizing the probability of absence of movement response. An estimated effect site concentration of 2 ng/ml of remifentanil was found to inhibit movement to a tetanic stimulation with a probability of 81%.

Effects of Letermovir and/or Methylprednisolone Coadministration on Voriconazole Pharmacokinetics in Hematopoietic Stem Cell Transplantation: A Population Pharmacokinetic Study

Objective

The aim of this study was to identify factors affecting blood concentrations of voriconazole following letermovir coadministration using population pharmacokinetic (PPK) analysis in allogeneic hematopoietic stem cell transplant (allo-HSCT) recipients.

Methods

The following data were retrospectively collected: voriconazole trough levels, patient characteristics, concomitant drugs, and laboratory information. PPK analysis was performed with NONMEM^® version 7.4.3, using the first-order conditional estimation method with interaction. We collected data on plasma voriconazole steady-state trough concentrations at 216 timepoints for 47 patients. A nonlinear pharmacokinetic model with the Michaelis–Menten equation was applied to describe the relationship between steady-state trough concentration and daily maintenance dose of voriconazole. After stepwise covariate modeling, the final model was evaluated using a goodness-of-fit plot, case deletion diagnostics, and bootstrap methods.

Results

The maximum elimination rate (V_max) of voriconazole in patients coadministered letermovir and methylprednisolone was 1.72 and 1.30 times larger than that in patients not coadministered these drugs, respectively, resulting in decreased voriconazole trough concentrations. The developed PPK model adequately described the voriconazole trough concentration profiles in allo-HSCT recipients. Simulations clearly showed that increased daily doses of voriconazole were required to achieve an optimal trough voriconazole concentration (1–5 mg/L) when patients received voriconazole with letermovir and/or methylprednisolone.

Conclusions

The development of individualized dose adjustment is critical to achieve optimal voriconazole concentration, especially among allo-HSCT recipients receiving concomitant letermovir and/or methylprednisolone.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40268-021-00365-0.

Towards Fixed Dosing of Tocilizumab in ICU-Admitted COVID-19 Patients: Results of an Observational Population Pharmacokinetic and Descriptive Pharmacodynamic Study

BACKGROUND AND OBJECTIVE

In the randomized controlled trial REMAP-CAP, it was shown that next to dexamethasone, the interleukin (IL)-6 receptor antagonist tocilizumab improves outcome, including survival in intensive care unit (ICU)-admitted coronavirus disease 2019 (COVID)-19 patients. Therefore tocilizumab has been added to many COVID-19 treatment guidelines. Because obesity is a risk factor for the development of severe COVID-19, concerns have been raised about overtreatment, as well as undertreatment, through weight-based dosing of tocilizumab. The currently applied dose of 8 mg/kg is based on the use of this drug for other indications, however it has not formally been investigated for COVID-19. In this study, the pharmacokinetics and pharmacodynamics of tocilizumab were investigated in ICU-admitted COVID-19 patients.

METHODS

This was an open-label, single-centre, observational population pharmacokinetic and descriptive pharmacodynamic evaluation study. Enrolled patients, with polymerase chain reaction-confirmed COVID-19 were admitted to the ICU for mechanical ventilation or high flow nasal canula oxygen support. All patients were 18 years of age or older and received intravenous tocilizumab 8 mg/kg (maximum 800 mg) within 24 h after admission to the ICU and received dexamethasone 6 mg daily as concomitant therapy. For evaluation of the pharmacokinetics and pharmacodynamics of tocilizumab, all time points from day 0 to 20 days after dose administration were eligible for collection. A nonlinear mixed-effects model was developed to characterize the population pharmacokinetic parameters of tocilizumab in ICU-admitted COVID-19 patients. Covariate analysis was performed to identify potential covariates for dose individualization. For the development of alternative dosing schedules, Monte Carlo simulations using the final model were performed.

RESULTS

Overall, 29 patients were enrolled between 15 December 2020 and 15 March 2021. A total of 139 tocilizumab plasma samples were obtained covering the pharmacokinetic curve of day 0 to day 20 after tocilizumab initiation. A population pharmacokinetic model with parallel linear and nonlinear clearance (CL) was developed and validated. Average CL was estimated to be 0.725 L/day, average volume of distribution (V_d) was 4.34 L, maximum elimination rate (V_max) was 4.19 μg/day, and concentration at which the elimination pathway is half saturated (K_m) was 0.22 μg/mL. Interindividual variability was identified for CL (18.9%) and V_d (21%). Average area under the concentration versus time curve from time zero to infinity of the first dose (AUC_{inf 1st DOSE}) was 938 [±190] μg/mL*days. All patients had tocilizumab exposure above 1 μg/mL for at least 15 days. Bodyweight-based dosing increases variability in exposure compared with fixed dosing.

CONCLUSIONS

This study provides evidence to support a fixed dose of tocilizumab 600 mg in COVID-19 patients. Fixed dosing is a safe, logistically attractive, and drug expenses saving alternative compared with the current 8 mg/kg recommendation.

Pharmacokinetics and Drug-Drug Interactions of Abacavir and Lamuvudine Co-administered With Antituberculosis Drugs in HIV-Positive Children Treated for Multidrug-Resistant Tuberculosis

Given the high prevalence of multidrug-resistant (MDR)-TB in high HIV burden settings, it is important to identify potential drug-drug interactions between MDR-TB treatment and widely used nucleoside reverse transcriptase inhibitors (NRTIs) in HIV-positive children. Population pharmacokinetic models were developed for lamivudine (n = 54) and abacavir (n = 50) in 54 HIV-positive children established on NRTIs; 27 with MDR-TB (combinations of high-dose isoniazid, pyrazinamide, ethambutol, ethionamide, terizidone, fluoroquinolones, and amikacin), and 27 controls without TB. Two-compartment models with first-order elimination and transit compartment absorption described both lamivudine and abacavir pharmacokinetics, respectively. Allometric scaling with body weight adjusted for the effect of body size. Clearance was predicted to reach half its mature value ∼ 2 (lamivudine) and ∼ 3 (abacavir) months after birth, with completion of maturation for both drugs at ∼ 2 years. No significant difference was found in key pharmacokinetic parameters of lamivudine and abacavir when co-administered with routine drugs used for MDR-TB in HIV-positive children.

Development of a genetic algorithm and NONMEM workbench for automating and improving population pharmacokinetic/pharmacodynamic model selection

The current approach to selection of a population PK/PD model is inherently flawed as it fails to account for interactions between structural, covariate, and statistical parameters. Further, the current approach requires significant manual and redundant model modifications that heavily lend themselves to automation. Within the discipline of numerical optimization it falls into the "local search" category. Genetic algorithms are a class of algorithms inspired by the mathematics of evolution. GAs are general, powerful, robust algorithms and can be used to find global optimal solutions for difficult problems even in the presence of non-differentiable functions, as is the case in the discrete nature of including/excluding model components in search of the best performing mixed-effects PK/PD model. A genetic algorithm implemented in an R-based NONMEM workbench for identification of near optimal models is presented. In addition to the GA capabilities, the workbench supports modeling efforts by: (1) Organizing and displaying models in tabular format, allowing the user to sort, filter, edit, create, and delete models seamlessly, (2) displaying run results, parameter estimates and precisions, (3) integrating xpose4 and PsN to facilitate generation of model diagnostic plots and run PsN scripts, (4) running regression models between post-hoc parameter estimates and covariates. This approach will further facilitate the scientist to shift efforts to focus on model evaluation, hypotheses generation, and interpretation and applications of resulting models.

Population repeated time-to-event analysis of exacerbations in asthma patients: A novel approach for predicting asthma exacerbations based on biomarkers, spirometry, and diaries/questionnaires

Identification of covariates, including biomarkers, spirometry, and diaries/questionnaires, that predict asthma exacerbations would allow better clinical predictions, shorter phase II trials and inform decisions on phase III design, and/or initiation (go/no-go). The objective of this work was to characterize asthma-exacerbation hazard as a function of baseline and time-varying covariates. A repeated time-to-event (RTTE) model for exacerbations was developed using data from a 52-week phase IIb trial, including 502 patients with asthma randomized to placebo or 70 mg, 210 mg, or 490 mg astegolimab every 4 weeks. Covariate analysis was performed for 20 baseline covariates using the full random effects modeling approach, followed by time-varying covariate analysis of nine covariates using the stepwise covariate model (SCM) building procedure. Following the SCM, an astegolimab treatment effect was explored. Diary-based symptom score (difference in objective function value [dOFV] of -83.7) and rescue medication use (dOFV = -33.5), and forced expiratory volume in 1 s (dOFV = -14.9) were identified as significant time-varying covariates. Of note, time-varying covariates become more useful with more frequent measurements, which should favor the daily diary scores over others. The most influential baseline covariates were exacerbation history and diary-based symptom score (i.e., symptom score was important as both time-varying and baseline covariate). A (nonsignificant) astegolimab treatment effect was included in the final model because the limited data set did not allow concluding the remaining effect size as irrelevant. Without time-varying covariates, the treatment effect was statistically significant (p < 0.01). This work demonstrated the utility of a population RTTE approach to characterize exacerbation hazard in patients with severe asthma.

Pharmacokinetics of Ceftazidime in Children and Adolescents with Obesity

PURPOSE

The aim of this study was to evaluate ceftazidime pharmacokinetics (PK) in a cohort that includes a predominate number of children and adolescents with obesity and assess the efficacy of competing dosing strategies.

METHODS

A population PK model was developed using opportunistically collected plasma samples. For each dosing strategy, model-based probability of target attainment (PTA) estimates were computed for study participants using empirical Bayes estimates. In addition, the effects of body size and renal function on PTA were evaluated using stochastic model simulations with virtually generated subjects.

RESULTS

Twenty-nine participants, 24 of whom were obese, contributed data towards the analysis. The median (range) age, body weight, and body mass index of participants were 12.2 years (2.3-20.6), 59.2 kg (8.4-121), and 25.2 kg/m2 (13.8-42.9), respectively. Administration of 50 mg/kg intravenously (IV) every 8 hours (q8h; max 6 g/day) or 40 mg/kg IV q6h (max 6 g/day) resulted in PTA values of ≥ 90% (minimum inhibitory concentration 8 mg/L) for the subset of obese participants with estimated glomerular filtration rates (GFR) ≥ ~ 80 mL/min/1.73 m2. However, for both regimens, stochastic model simulations denoted lower PTA values (< 90%) with increasing body weight for virtual subjects with GFR ≥ 120 mL/min/1.73 m2. Alternatively, permitting for a maximum daily dose of 8 g/day using a 40 mg/kg IV q6h regimen provided PTA values that were near or above target (90%) for virtual subjects between 10 to 120 kg with GFR ≥ 80 mL/min/1.73 m2.

CONCLUSION

Our analysis suggests administration of 40 mg/kg IV q6h (max 8 g/day) maximizes PTA in children and adolescents with obesity and GFR ≥ 80 mL/min/1.73 m2.

TRIAL REGISTRATION

Clinicaltrials.gov Identifier: NCT01431326.

Evaluation of Pharmacokinetics and Pharmacodynamics of Deferasirox in Pediatric Patients

Background: Deferasirox (DFX) is commonly used to reduce the chronic iron overload (IO) in pediatric patients. However, the drug is characterized by a large pharmacokinetic variability and approximately 10% of patients may discontinue the treatment due to toxicities. Therefore, the present retrospective study investigated possible correlations between DFX pharmacokinetics and drug-associated toxicities in 39 children (26 males), aged 2–17 years, who underwent an allogeneic hematopoietic stem cell transplantation. Methods: IO was diagnosed by an abdominal magnetic resonance imaging and DFX was started at a median dose of 500 mg/day. DFX plasma concentrations were measured by a high performance liquid chromatographic method with UV detection and they were analysed by nonlinear mixed-effects modeling. Results: The pharmacometric analysis demonstrated that DFX pharmacokinetics were significantly influenced by lean body mass (bioavailability and absorption constant), body weight (volume of distribution), alanine and aspartate transaminases, direct bilirubin, and serum creatinine (clearance). Predicted DFX minimum plasma concentrations (C_trough) accounted for 32.4 ± 23.2 mg/L (mean ± SD), and they were significantly correlated with hepatic/renal and hematological toxicities (p-value < 0.0001, T-test and Fisher’s exact tests) when C_trough threshold values of 7.0 and 11.5 mg/L were chosen, respectively. Conclusions: The population pharmacokinetic model described the interindividual variability and identified C_trough threshold values that were predictive of hepatic/renal and hematological toxicities associated with DFX.

Impact of a vancomycin loading dose on the achievement of target vancomycin exposure in the first 24 h and on the accompanying risk of nephrotoxicity in critically ill patients

Background

The advocated pharmacokinetic/pharmacodynamic (PK/PD) target for vancomycin, AUC/MIC ≥ 400 mg·h/L, may not be reached with a conventional fixed starting dose of 1000 mg in critically ill patients, but increasing the dose may cause nephrotoxicity.

Objectives

To evaluate the effect of a weight-based loading dose of 25 mg/kg vancomycin on PK/PD target attainment in the first 24 h (AUC_0–24) in critically ill patients and to evaluate whether this increases the risk of acute kidney injury (AKI).

Patients and methods

A prospective observational before/after study was performed in ICU patients, comparing the percentage of vancomycin courses with AUC_0–24 ≥ 400 mg·h/L and the incidence of AKI, defined as worsening of the risk, injury, failure, loss of kidney function and end-stage kidney disease (RIFLE) score. The conventional dose group received 1000 mg of vancomycin as initial dose; the loading dose group received a weight-based loading dose of 25 mg/kg. A population PK model developed using non-linear mixed-effects modelling was used to estimate AUC_0–24 in all patients.

Results

One hundred and four courses from 82 patients were included. With a loading dose, the percentage of courses achieving AUC_0–24 ≥ 400 mg·h/L increased significantly from 53.8% to 88.0% (P = 0.0006). The percentage of patients with new-onset AKI was not significantly higher when receiving a 25 mg/kg loading dose (28.6% versus 37.8%; P = 0.48). However, the risk of AKI was significantly higher in patients achieving AUC_0–24 > 400 mg·h/L compared with patients achieving AUC < 400 mg·h/L (39.0% versus 14.8%; P = 0.031).

Conclusions

A weight-based loading dose of 25 mg/kg vancomycin led to significantly more patients achieving AUC_0–24 ≥ 400 mg·h/L without increased risk of AKI. However, some harm cannot be ruled out since higher exposure was associated with increased risk of AKI.

The Course of AαVal541 as a Proteinase 3 Specific Neo-Epitope after Alpha-1-Antitrypsin Augmentation in Severe Deficient Patients

In alpha-1-antitrypsin deficiency (AATD), neutrophil serine proteases such as elastase and proteinase 3 (PR3) are insufficiently inhibited. A previous study in AATD patients showed a higher plasma level of the specific PR3-generated fibrinogen-derived peptide AαVal541, compared with healthy controls. Here, we analyzed the course of AαVal541 plasma levels during 4 weeks after a single iv dose of 240 mg/kg AAT in ten patients with genotype Z/Rare or Rare/Rare. To this end, we developed an immunoassay to measure AαVal541 in plasma and applied population pharmacokinetic modeling for AAT. The median AαVal541 plasma level before treatment was 140.2 nM (IQR 51.5–234.8 nM)). In five patients who received AAT for the first time, AαVal541 levels decreased to 20.6 nM (IQR 5.8–88.9 nM), and in five patients who already had received multiple infusions before, it decreased to 26.2 nM (IQR 22.31–35.0 nM). In 9 of 10 patients, AαVal541 levels were reduced to the median level of healthy controls (21.4 nM; IQR 16.7–30.1 nM). At 7–14 days after treatment, AαVal541 levels started to increase again in all patients. Our results show that fibrinopeptide AαVal541 may serve as a biochemical footprint to assess the efficacy of in vivo inhibition of PR3 activity in patients receiving intravenous AAT augmentation therapy.

A Repeated Time-to-Positive Symptoms Improvement among Malaysian Patients with Schizophrenia Spectrum Disorders Treated with Clozapine

Clozapine remains the drug of choice for resistant schizophrenia. However, its dose-response relationship is still controversial. The current investigation aimed to develop a repeated time-to-positive symptoms improvement following the onset of clozapine treatment in Malaysian schizophrenia spectrum disorder patients. Data from patients’ medical records in the Psychiatric Clinic, Penang General Hospital, were retrospectively analyzed. Several parametric survival models were evaluated using nonlinear mixed-effect modeling software (NONMEM 7.3.0). Kaplan–Meier-visual predictive check (KM-VPC) and sampling-importance resampling (SIR) methods were used to validate the final model. A total of 116 patients were included in the study, with a mean follow-up of 306 weeks. Weibull hazard function best fitted the data. The hazard of positive symptoms improvement decreased 4% for every one-year increase in age over the median of 41 years (adjusted hazard ratio (aHR), 0.96; 95% confidence intervals (95% CI), (0.93–0.98)). However, patients receiving a second atypical antipsychotic agent had four-folds higher hazard (aHR, 4.01; 95% CI, (1.97–7.17)). The hazard increased 2% (aHR, 1.02; 95% CI, (1.01–1.03)) for every 1 g increase in the clozapine six months cumulative dose over the median of 34 g. The developed model provides essential information on the hazard of positive symptoms improvement after the first clozapine dose administration, including modifiable predictors of high clinical importance.

Population Pharmacokinetics Modeling of Lamotrigine in Jordanian Epileptic Patients Using Dried Blood Spot Sampling

AIMS

To characterize the population pharmacokinetics of lamotrigine in Jordanian epileptic patients and to identify factors affecting therapeutic parameters.

PATIENTS AND METHODS

A population pharmacokinetics model for lamotrigine was established based on a prospectively collected data of 52 steady-state concentrations from 38 adult and pediatric patients with epilepsy. Lamotrigine concentrations were determined by a dried blood spot liquid chromatography method. Data were analyzed according to a one-compartment model with first-order absorption and elimination using the nonlinear mixed effect modeling program. The covariates effect of total body weight, gender, age, and co-medication with topiramate, carbamazepine, phenytoin, phenobarbital, and valproic acid on lamotrigine clearance were investigated using a stepwise forward addition followed by a stepwise backward elimination.

RESULTS

The final population pharmacokinetics model for lamotrigine clearance was as follows: CL/F_pop=θ₁*exp ^(θ3*age)*exp ^{(θ5*carbamazepine)}*exp ^{(θ6*valproic acid)} , where θ₁ is the relative clearance (L/hr) estimated, and θ₃, θ₅, and θ₆ are the fixed parameters relating to age and co-medication with carbamazepine and valproic acid, respectively.The population mean value of lamotrigine total clearance generated in the final model (with covariates) was 2.12 L/hr. Inter-individual variability and residual unexplained variability expressed as the coefficient of variation was 37.1 and 26.1%, respectively.

CONCLUSION

Lamotrigine total clearance in the Jordanian patients is comparable to that reported by others for Caucasian patients. Age and concomitant therapy with carbamazepine and valproic acid significantly affected lamotrigine clearance, and accounted for 48% of its inter-individual variability.

Linking categorical models for prediction of pleasantness score using individual predictions of sweetness and creaminess: An advancement of categorical modeling

The aim of this work was to develop and evaluate approaches of linked categorical models using individual predictions of probability. A model was developed using data from a study which assessed the perception of sweetness, creaminess, and pleasantness in dairy solutions containing variable concentrations of sugar and fat. Ordered categorical models were used to predict the individual sweetness and creaminess scores and these individual predictions were used as covariates in the model of pleasantness response. The model using individual predictions was compared to a previously developed model using the amount of fat and sugar as covariates driving pleasantness score. The model using the individual prediction of odds of sweetness and creaminess had a lower variability of pleasantness than the model using the content of sugar and fat in the test solutions, which indicates that the individual odds explain part of the variability in pleasantness. Additionally, simultaneous and sequential modeling approaches were compared for the linked categorical model. Parameter estimation was similar, but precision was better with sequential modeling approaches compared to the simultaneous modeling approach. The previous model characterizing the pleasantness response was improved by using individual predictions of sweetness and creaminess rather than the amount of fat and sugar in the solution. The application of this approach provides an advancement within categorical modeling showing how categorical models can be linked to enable the utilization of individual prediction. This approach is aligned with biology of taste sensory which is reflective of the individual perception of sweetness and creaminess, rather than the amount of fat and sugar in the solution.

Characterisation of pharmacokinetics, safety and tolerability in a first-in-human study for AZD8154, a novel inhaled selective PI3Kγδ dual inhibitor targeting airway inflammatory disease

AIMS

This 3-part, randomised, phase 1 first-in-human study (NCT03436316) investigated the safety, tolerability and pharmacokinetics (PK) of AZD8154, a dual phosphoinositide 3-kinase (PI3K) γδ inhibitor developed as a novel inhaled anti-inflammatory treatment for respiratory disease.

METHODS

Healthy men, and women of nonchildbearing potential, were enrolled to receive single and multiple ascending inhaled doses of AZD8154 in parts 1 and 3 of the study, respectively, while part 2 characterised the systemic PK after a single intravenous (IV) dose. In part 1, participants received 0.1-7.7 mg AZD8154 in 6 cohorts. In part 2, participants were given 0.15 mg AZD8154 as an IV infusion. In part 3, AZD8154 was given in 3 cohorts of 0.6, 1.8 and 3.1 mg, with a single dose on Day 1 followed by repeated once-daily doses on Days 4-12.

RESULTS

In total, 78 volunteers were randomised. All single inhaled, single IV and multiple inhaled doses were shown to be well tolerated without any safety concerns. A population PK model, using nonlinear mixed-effect modelling, was developed to describe the PK of AZD8154. The terminal mean half-life of AZD8154 was 18.0-32.0 hours. The geometric mean of the absolute pulmonary bioavailability of AZD8154 via the inhaled route was 94.1%.

CONCLUSION

AZD8154 demonstrated an acceptable safety profile, with no reports of serious adverse events and no clinically significant drug-associated safety concerns reported in healthy volunteers. AZD8154 demonstrated prolonged lung retention and a half-life supporting once-daily dosing.

Model-informed precision dosing for alemtuzumab in paediatric and young adult patients undergoing allogeneic haematopoietic cell transplantation

Alemtuzumab is a lymphodepleting monoclonal antibody utilized in conditioning regimens for allogeneic haematopoietic cell transplantation (HCT). A recently proposed therapeutic range of 0.15-0.6 μg/mL on the day of transplantation is associated with better HCT outcomes. The purpose of this study was to characterize alemtuzumab population pharmacokinetic/pharmacodynamic (PK/PD) and to propose individualized subcutaneous dosing schemes to achieve this optimal level for paediatric patients.

METHODS

Alemtuzumab concentration and absolute lymphocyte count (ALC) profiles were obtained from 29 paediatric and young adult patients (median age 6.4 y; range 0.28-21.4 y) with nonmalignant disorders undergoing HCT. PK/PD analyses were performed using nonlinear mixed effects modelling. Monte Carlo simulation was conducted to evaluate different improved dosing approaches.

RESULTS

A one-compartment model with sequential zero- and first-order absorption adequately described subcutaneously administered alemtuzumab PK. Model fit was significantly improved by including allometrically scaled body weight on clearance (0.080 L/h/70 kg) and volume of distribution (17.4 L/70 kg). ALC reduction following subcutaneous alemtuzumab was swift. An inhibitory E_max model best characterized the relationship between alemtuzumab concentration and ALC. E_max and EC₅₀ were estimated as 1.18 × 10³ /μL and 0.045 μg/mL, respectively. The currently used per kg dosing was found to cause uneven alemtuzumab exposure across different age and weight cohorts. Simulations indicated optimal target achieving dose as allometry-based dose of 18 mg × (weight/70)^0.75 or body surface area-based dose of 10 mg/m² , divided over 3 days, with a potential individualized top-up dose; both of which yielded similar results.

CONCLUSION

An allometry- or body surface area-based starting dosing regimen in combination with individualized Bayesian PK estimation using concentration feedback is proposed for alemtuzumab precision dosing in children undergoing allogeneic HCT.

Population pharmacokinetics and pharmacogenetics of ritonavir-boosted darunavir in the presence of raltegravir or tenofovir disoproxil fumarate/emtricitabine in HIV-infected adults and the relationship with virological response: a sub-study of the NEAT001/ANRS143 randomized trial

OBJECTIVES

NEAT001/ANRS143 demonstrated non-inferiority of once-daily darunavir/ritonavir (800/100 mg) + twice-daily raltegravir (400 mg) versus darunavir/ritonavir + tenofovir disoproxil fumarate/emtricitabine (245/200 mg once daily) in treatment-naive patients. We investigated the population pharmacokinetics of darunavir, ritonavir, tenofovir and emtricitabine and relationships with demographics, genetic polymorphisms and virological failure.

METHODS

Non-linear mixed-effects models (NONMEM v. 7.3) were applied to determine pharmacokinetic parameters and assess demographic covariates and relationships with SNPs (SLCO3A1, SLCO1B1, NR1I2, NR1I3, CYP3A5*3, CYP3A4*22, ABCC2, ABCC10, ABCG2 and SCL47A1). The relationship between model-predicted darunavir AUC0-24 and C24 with time to virological failure was evaluated by Cox regression.

RESULTS

Of 805 enrolled, 716, 720, 347 and 361 were included in the darunavir, ritonavir, tenofovir and emtricitabine models, respectively (11% female, 83% Caucasian). No significant effect of patient demographics or SNPs was observed for darunavir or tenofovir apparent oral clearance (CL/F); coadministration of raltegravir did not influence darunavir or ritonavir CL/F. Ritonavir CL/F decreased by 23% in NR1I2 63396C>T carriers and emtricitabine CL/F was linearly associated with creatinine clearance (P<0.001). No significant relationship was demonstrated between darunavir AUC0-24 or C24 and time to virological failure [HR (95% CI): 2.28 (0.53-9.80), P=0.269; and 1.82 (0.61-5.41), P=0.279, respectively].

CONCLUSIONS

Darunavir concentrations were unaltered in the presence of raltegravir and not associated with virological failure. Polymorphisms investigated had little impact on study-drug pharmacokinetics. Darunavir/ritonavir + raltegravir may be an appropriate option for patients experiencing NRTI-associated toxicity.

Scipion PKPD: an Open-Source Platform for Biopharmaceutics, Pharmacokinetics and Pharmacodynamics Data Analysis

PURPOSE

Biopharmaceutics examines the interrelationship of the drug's physical/chemical properties, the dosage form (drug product) in which the drug is given, and the administration route on the rate and extent of sys- temic drug absorption. Pharmacokinetics is the study of the movement of drugs in the body. It uses mathematical models to evaluate the movement of absorption, distribution, metabolism, and excretion (ADME) within an organism. Finally, Pharmacodynamics is the analysis of how these drugs af- fect that organism. Pharmacokinetics data normally comes in samples over time of the drug concentration either in plasma or in some specific tissue. Similarly, pharmacodynamics data comes normally in samples over time of some quantity of interest (biophysical quantity like temperature, blood pres- sure, etc.). The data is submitted to a non-parametric analysis, in which a description of the observed data is reported (e.g., the Area Under the Curve), or to a parametric analysis by fitting a model (normally based on differential equations) so that prediction about future events can be made. This paper aims to introduce Scipion PKPD, an open-source platform for data analysis of this kind in the three domains (Biopharmaceutics, Pharmacokinetics, and Pharmacodynamics). The platform implements the most popular models and is open to new ones. The platform provides almost 100 different high-level operations that we call protocols.

METHODS

We have developed a Python module integrated into the work- flow engine Scipion. The plugin implements the numerical analysis and meta- data handling tools to address multiple problems (see Suppl. Material for a detailed list of the tasks solved).

RESULTS

We illustrate the use of this package with an integrative exam- ple that involves all these areas.

CONCLUSIONS

We show that the package successfully addresses these kinds of analyses. Scipion PKPD is freely available at https://github. com/cossorzano/scipion-pkpd .

Evaluation of body-surface-area adjusted dosing of high-dose methotrexate by population pharmacokinetics in a large cohort of cancer patients

Background

The aim of this study was to identify sources of variability including patient gender and body surface area (BSA) in pharmacokinetic (PK) exposure for high-dose methotrexate (MTX) continuous infusion in a large cohort of patients with hematological and solid malignancies.

Methods

We conducted a retrospective PK analysis of MTX plasma concentration data from hematological/oncological patients treated at the University Hospital of Cologne between 2005 and 2018. Nonlinear mixed effects modeling was performed. Covariate data on patient demographics and clinical chemistry parameters was incorporated to assess relationships with PK parameters. Simulations were conducted to compare exposure and probability of target attainment (PTA) under BSA adjusted, flat and stratified dosing regimens.

Results

Plasma concentration over time data (2182 measurements) from therapeutic drug monitoring from 229 patients was available. PK of MTX were best described by a three-compartment model. Values for clearance (CL) of 4.33 [2.95–5.92] L h^− 1 and central volume of distribution of 4.29 [1.81–7.33] L were estimated. An inter-occasion variability of 23.1% (coefficient of variation) and an inter-individual variability of 29.7% were associated to CL, which was 16 [7–25] % lower in women. Serum creatinine, patient age, sex and BSA were significantly related to CL of MTX. Simulations suggested that differences in PTA between flat and BSA-based dosing were marginal, with stratified dosing performing best overall.

Conclusion

A dosing scheme with doses stratified across BSA quartiles is suggested to optimize target exposure attainment. Influence of patient sex on CL of MTX is present but small in magnitude.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08443-x.

IV and oral fosfomycin pharmacokinetics in neonates with suspected clinical sepsis

BACKGROUND

Fosfomycin has the potential to be re-purposed as part of a combination therapy to treat neonatal sepsis where resistance to current standard of care (SOC) is common. Limited data exist on neonatal fosfomycin pharmacokinetics and estimates of bioavailability and CSF/plasma ratio in this vulnerable population are lacking.

OBJECTIVES

To generate data informing the appropriate dosing of IV and oral fosfomycin in neonates using a population pharmacokinetic analysis of plasma and CSF data.

METHODS

The NeoFosfo study (NCT03453177) was a randomized trial that examined the safety and pharmacokinetics of fosfomycin comparing SOC versus SOC plus fosfomycin. Sixty-one neonates received fosfomycin (100 mg/kg IV q12h for 48 h) and then they converted to oral therapy at the same dose. Two plasma pharmacokinetic samples were taken following the first IV and oral doses, sample times were randomized to cover the whole pharmacokinetic profile and opportunistic CSF pharmacokinetic samples were collected. A population pharmacokinetic model was developed in NONMEM and simulations were performed.

RESULTS

In total, 238 plasma and 15 CSF concentrations were collected. A two-compartment disposition model, with an additional CSF compartment and first-order absorption, best described the data. Bioavailability was estimated as 0.48 (95% CI = 0.347-0.775) and the CSF/plasma ratio as 0.32 (95% CI = 0.272-0.409). Allometric weight and postmenstrual age (PMA) scaling was applied; additional covariates included postnatal age (PNA) on clearance and CSF protein on CSF/plasma ratio.

CONCLUSIONS

Through this analysis a population pharmacokinetic model has been developed that can be used alongside currently available pharmacodynamic targets to select a neonatal fosfomycin dose based on an infant's PMA, PNA and weight.

Population Pharmacokinetics of Oxycodone and Metabolites in Patients with Cancer-Related Pain

Simple Summary

Patients with moderate to severe cancer-related pain are frequently treated with oxycodone, a strong-acting opioid. However, treatment with oxycodone does not always lead to sufficient analgesic action. In order to determine which factors affect treatment outcomes, we performed an observational study and developed a population pharmacokinetic model. The model described oxycodone, nor-oxycodone and nor-oxymorphone pharmacokinetics. The association between oxycodone or oxycodone metabolites’ exposure with pain scores and adverse events was not significant. The combined oxycodone, nor-oxycodone and nor-oxymorphone model is a good starting point for further unravelling the factors that affect the pharmacokinetic/pharmacodynamic relation of oxycodone and its metabolites.

Abstract

Oxycodone is frequently used for treating cancer-related pain, while not much is known about the factors that influence treatment outcomes in these patients. We aim to unravel these factors by developing a population-pharmacokinetic model to assess the pharmacokinetics of oxycodone and its metabolites in cancer patients, and to associate this with pain scores, and adverse events. Hospitalized patients with cancer-related pain, who were treated with oral oxycodone, could participate. Pharmacokinetic samples and patient-reported pain scores and occurrence and severity of nine adverse events were taken every 12 h. In 28 patients, 302 pharmacokinetic samples were collected. A one-compartment model for oxycodone and each metabolite best described oxycodone, nor-oxycodone, and nor-oxymorphone pharmacokinetics. Furthermore, oxycodone exposure was not associated with average and maximal pain scores, and oxycodone, nor-oxycodone, and nor-oxymorphone exposure were not associated with adverse events (all p > 0.05). This is the first model to describe the pharmacokinetics of oxycodone including the metabolites nor-oxycodone and nor-oxymorphone in hospitalized patients with cancer pain. Additional research, including more patients and a more timely collection of pharmacodynamic data, is needed to further elucidate oxycodone (metabolite) pharmacokinetic/pharmacodynamic relationships. This model is an important starting point for further studies to optimize oxycodone dosing regiments in patients with cancer-related pain.

Improved Decision-Making Confidence Using Item-Based Pharmacometric Model: Illustration with a Phase II Placebo-Controlled Trial

This study aimed to illustrate how a new methodology to assess clinical trial outcome measures using a longitudinal item response theory–based model (IRM) could serve as an alternative to mixed model repeated measures (MMRM). Data from the EXACT (Exacerbation of chronic pulmonary disease tool) which is used to capture frequency, severity, and duration of exacerbations in COPD were analyzed using an IRM. The IRM included a graded response model characterizing item parameters and functions describing symptom-time course. Total scores were simulated (month 12) using uncertainty in parameter estimates. The 50th (2.5th, 97.5th) percentiles of the resulting simulated differences in average total score (drug minus placebo) represented the estimated drug effect (95%CI), which was compared with published MMRM results. Furthermore, differences in sample size, sensitivity, specificity, and type I and II errors between approaches were explored. Patients received either oral danirixin 75 mg twice daily (n = 45) or placebo (n = 48) on top of standard of care over 52 weeks. A step function best described the COPD symptoms-time course in both trial arms. The IRM improved precision of the estimated drug effect compared to MMRM, resulting in a sample size of 2.5 times larger for the MMRM analysis to achieve the IRM precision. The IRM showed a higher probability of a positive predictive value (34%) than MMRM (22%). An item model–based analysis data gave more precise estimates of drug effect than MMRM analysis for the same endpoint in this one case study.

Population Pharmacokinetic Modelling and Simulation to Determine the Optimal Dose of Nanoparticulated Sorafenib to the Reference Sorafenib

Sorafenib, an oral multikinase inhibitor, exhibits a highly variable absorption profile due to enterohepatic reabsorption and poor solubility. SYO-1644 improved the solubility of sorafenib by nanoparticulation technology leading to enhanced bioavailability. To evaluate the pharmacokinetically equivalent dose of SYO-1644 to the reference Nexavar^® 200 mg, a randomized, open-label, replicated two-period study was conducted in healthy volunteers. A total of 32 subjects orally received a single dose of the following assigned treatment under a fasted state in the first period and repeated once more in the second period with a two-week washout: SYO-1644 100, 150 and 200 mg and Nexavar^® 200 mg. Pharmacokinetic (PK) samples were collected up to 168 h post-dose. The PK profile was evaluated by both non-compartmental analysis and population PK method. With the final model, 2 × 2 crossover trial scenarios with Nexavar^® 200 mg and each dose of SYO-1644 ranging from 100 to 150 mg were repeated 500 times by Monte Carlo simulation, and the proportion of bioequivalence achievement was assessed. Transit absorption compartments, followed by a one-compartment model with first-order elimination and enterohepatic reabsorption components were selected as the final model. The simulation results demonstrated that the SYO-1644 dose between 120 and 125 mg could yielded the highest proportion of bioequivalence.

Postprandial triglyceride reduction following acute treatment of a selective 5-hydroxytryptamine-2c agonist and characterization using a semi-physiological model

AIM

To investigate the tolerability, pharmacokinetics (PK) and postprandial triglyceride (TG) response of single, escalating oral doses of a selective 5-hydroxytryptamine-2c (5-HT_2c ) agonist in subjects with overweight/obesity and apply mechanistic population pharmacokinetic-pharmacodynamic modelling to identify a plausible drug mechanism of action.

MATERIALS AND METHODS

This phase 1, single-centre, double-blind, randomized, placebo-controlled, four-period, two-alternating cohorts study evaluated single escalating oral doses ranging from 5 to 130 mg of LY2140112 (LY) in subjects with overweight/obesity (body mass index: 27-39 kg/m² ). Postprandial TG response (total TG, chylomicrons and very low-density lipoprotein particles [VLDL]-V6) following a high-fat meal were assessed for 11 h postmeal for each dose level. The PK profile was assessed for 96 h postdose. Drug exposure and TG concentrations in chylomicrons and VLDL-V6 were used to characterize the drug mechanism of action using non-linear mixed-effect modelling.

RESULTS

Seventeen subjects entered the study and 16 subjects received at least one dose of LY. LY2140112 was generally well tolerated up to 75 mg. The PK of LY were described by a two-compartment model with first-order elimination. The 100 and 130 mg dose levels of LY significantly reduced the postprandial TG of VLDL-V6 by approximately 50%, while total TG and chylomicrons were not significantly different from placebo. The application of a published lipokinetic model successfully described the postprandial TG response in this study and indicated that LY reduced the conversion of TGs from chylomicron to VLDL-V6.

CONCLUSIONS

LY significantly reduced the postprandial TG of VLDL-V6 following a single dose, when food consumption was controlled. The data indicate that a selective 5-HT_2c agonist alters lipid metabolism, beyond the reported reduction in satiety. The application of a semi-physiological lipokinetic model enabled identification of a plausible drug mechanism of action of LY.

Methotrexate Triglutamate as a Determinant of Clinical Response in Mexican Patients With Rheumatoid Arthritis: Pharmacokinetics and Dose Recommendation

Methotrexate is the gold standard treatment in rheumatoid arthritis. Once absorbed, it is internalized in cells, where glutamate residues are added to produce polyglutamated forms, which are responsible for the effect of methotrexate. The aim of the current study is to determine the relationship between methotrexate triglutamate concentrations and the clinical evolution in rheumatoid arthritis patients, as well as to characterize the variability in both features to propose strategies for low-dose methotrexate optimization. The quantification of methotrexate triglutamate concentration in red blood cells was performed through ultra-performance liquid chromatography coupled with mass spectrometry. Polymorphisms of genes involved in the formation of polyglutamates were determined by real-time polymerase chain reaction. A multivariate regression was performed to determine the covariates involved in the variability of methotrexate triglutamate concentrations and a population pharmacokinetics model was developed through nonlinear mixed-effects modeling. Disease activity score changed according to methotrexate triglutamate concentrations; patients with good response to treatment had higher concentrations than moderate or nonresponding patients. The methotrexate triglutamate concentrations were related to time under treatment, dose, red blood cells, and body mass index. A 1-compartment open model was selected to estimate the pharmacokinetic parameters; the typical total clearance (L/day) was determined as 1.45 * (body mass index/28 kg/m² ) * (red blood cells/4.6 × 10⁶ cells/μL) and the volume of distribution was 52.4 L, with an absorption rate of 0.0346/day and a fraction metabolized of 1.03%. Through the application of the model, the initial dose of methotrexate is proposed on the basis of stochastic simulations and considering methotrexate triglutamate concentrations found in responders patients.