Reporting a Population Pharmacokinetic–Pharmacodynamic Study: A Journal’s Perspective

Population pharmacokinetics of olanzapine in pediatric patients with psychiatric disorders

OBJECTIVE

To develop and validate a population pharmacokinetic (PPK) model of oral olanzapine in pediatric Chinese patients in order to individualize therapy in this population.

METHODS

A total of 897 serum concentrations from 269 pediatric patients taking oral olanzapine (ages 8-17 years) were collected. Demographic parameters, biological characteristics and concomitant medications were investigated as covariates. The data were analyzed using a nonlinear mixed-effects modeling approach. Bootstrapping (1000 runs), normalized prediction distribution error (NPDE), and external validation of 62 patients were employed. Simulations were performed to explore the individualized dosing regimens in various situations.

RESULTS

The one-compartment model with first-order absorption and elimination had an apparent clearance (CL/F) of 10.38 L/h, a distribution volume (V/F) of 9.41 L/kg and an absorption rate constant (Ka) fixed at 0.3 h-1. The equation was CL∕F (L∕h) = 10.38 × (body weight∕60)0.25 ×1.33 (if male) × 0.71 (if co-occurrence of infection) × 0.51 (if co-therapy with fluvoxamine) × 1.27 (if co-therapy with sertraline) × 1.43 (if co-therapy with valproate). The final model had satisfactory stability, robustness, and predictive ability. The results from a simulation suggested the oral olanzapine doses required for male and female pediatric patients weighing between 40 and 60 kg without co-medication were 10-15 mg/day and 7.5-10 mg/day, respectively, and dosage adjustments should be based on sex and body weight; and co-administrated with valproate, sertraline, or fluvoxamine.

CONCLUSION

This model may help individualize optimum dosing of oral olanzapine for pediatric patients.

Population pharmacokinetic analyses of methotrexate in pediatric patients: a systematic review

BACKGROUND AND OBJECTIVES

Methotrexate is widely utilized in the chemotherapy of malignant tumors and autoimmune diseases in the pediatric population, but dosing can be challenging. Several population pharmacokinetic models were developed to characterize factors influencing variability and improve individualization of dosing regimens. However, significant covariates included varied across studies. The primary objective of this review was to summarize and discuss population pharmacokinetic models of methotrexate and covariates that influence pharmacokinetic variability in pediatric patients.

METHODS

Systematic searches were conducted in the PubMed and EMBASE databases from inception to 7 July 2023. Reporting Quality was evaluated based on a checklist with 31 items. The characteristics of studies and information for model construction and validation were extracted, summarized, and discussed.

RESULTS

Eighteen studies (four prospective studies and fourteen retrospective studies with sample sizes of 14 to 772 patients and 2.7 to 93.1 samples per patient) were included in this study. Two-compartment models were the commonly used structural models for methotrexate, and the clearance range of methotrexate ranged from 2.32 to 19.03 L/h (median: 6.86 L/h). Body size and renal function were found to significantly affect the clearance of methotrexate for pediatric patients. There were limited reports on the role of other covariates, such as gene polymorphisms and co-medications, in the pharmacokinetic parameters of methotrexate pediatric patients. Internal and external evaluations were used to assess the performance of the population pharmacokinetic models.

CONCLUSION

A more rigorous external evaluation needs to be performed before routine clinical use to select the appropriate PopPK model. Further research is necessary to incorporate larger cohorts or pool analyses in specific susceptible pediatric populations to improve the understanding of predicted exposure profiles and covariate identification.

Development and Quality Control of a Population Pharmacokinetic Model Library for Caspofungin

BACKGROUND

Caspofungin is an echinocandin antifungal agent commonly used as the first-line therapy for invasive candidiasis, salvage therapy for invasive aspergillosis, and empirical therapy for presumed fungal infections. Pharmacokinetic (PK) variabilities and suboptimal exposure have been reported for caspofungin, increasing the risk of insufficient efficacy.

OBJECTIVE

This work aimed to develop a caspofungin population pharmacokinetic (popPK) library and demonstrate its utility by assessing the probability of target attainment across diverse settings.

METHODS

We established a caspofungin popPK model library following a rigorous literature review, re-implementing selected models in R with rxode2. Quality control procedures included a comparison of different studies and assessing covariate impacts. Model libraries were primarily used to perform Monte Carlo simulations to estimate target attainment and guide personalized dosing in Candida infections.

RESULTS

A total of 13 models, one- or two-compartment models, were included. The most significant covariates were body size (weight and body surface area), liver function, and albumin level. The results show that children and adults showed considerable differences in pharmacokinetics. For C. albicans and C. parapsilosis, none of the populations achieved a PTA of ≥90% at their respective susceptible MIC values. In contrast, for C. glabrata, 70% of the adult studies reached a PTA of ≥90%, while all pediatric studies achieved the same PTA level.

CONCLUSION

At the recommended dosage, adult patients showed notably lower exposure to caspofungin compared to pediatric patients. Considering body size, liver function, and serum albumin is crucial when determining caspofungin dosage regimens. Furthermore, further research is required to comprehensively understand the pharmacokinetics of caspofungin in pediatric patients.

Population pharmacokinetics of immunoglobulin G after intravenous, subcutaneous, or hyaluronidase-facilitated subcutaneous administration in immunoglobulin-naive patients with primary immunodeficiencies

Immunoglobulin G (IgG) replacement therapy is the standard of care for patients with primary immunodeficiencies with antibody deficiencies. Intravenous (IVIG), subcutaneous (SCIG), and hyaluronidase-facilitated subcutaneous immunoglobulin (fSCIG) therapies differ in their pharmacokinetic (PK) profiles, administration routes, and dosing regimens. Information on use of subcutaneous therapy in IgG treatment-naive patients is limited. This study used population pharmacokinetic (popPK) model-based simulations to characterize IgG PKs in IgG-naive patients with varying disease severity across several IVIG, SCIG, and fSCIG dosing regimens. An integrated popPK model, developed and validated using data from eight clinical trials, was utilized to simulate scenarios that varied by therapy, loading regimen, maintenance dose (equivalent to 400, 600, or 800 mg/kg every 4 weeks [Q4W]), and baseline endogenous total IgG concentration (1.5 or 4.0 g/L). Simulations were performed for age groups of 2-<6, 6-<12, 12-<18, and ≥18 years. Steady-state serum trough IgG concentrations (Cmin,ss), proportion of patients achieving Cmin,ss ≥ 7 g/L, and days taken to reach this threshold were summarized. SCIG provided greater mean Cmin,ss values than IVIG and fSCIG for any scenario. Across all therapies, Cmin,ss tended to increase with age, dose, and endogenous concentration. Although the findings are model-based and not a summarization of real-world observations, doses ≥ 800 mg/kg Q4W with corresponding loading regimens are likely to be clinically appropriate for achieving target IgG concentrations in treatment-naive patients in a timely manner, especially at low endogenous starting concentrations. Therapy-specific dose adjustment based on baseline endogenous IgG concentration, clinical status, and patient characteristics may be warranted.

Insights into the population pharmacokinetics and pharmacodynamics of quetiapine: a systematic review

INTRODUCTION

Quetiapine exhibits notable pharmacokinetic and pharmacodynamic (PK/PD) variability, the origins of which are poorly understood. This systematic review summarizes published population PK/PD studies and identifies significant covariates accounting for this variability to inform precision dosing.

METHODS

We systematically searched the PubMed, Web of Science, and Embase databases and compared study characteristics, model parameters, and covariate effects. Visual predictive distributions were used to compare different models. Forest plots and Monte Carlo simulations were used to assess the influence of covariates.

RESULTS

Six population PK and three population PK/PD studies were included. The median apparent clearance in adults was 87.7 L/h. Strong and moderate cytochrome P450 3A4 inducers increased the apparent clearance approximately fourfold, while strong cytochrome P450 3A4 inhibitors reduced it by 93%. The half-maximum effect concentrations were 82.8 ng/mL for the Brief Psychiatric Rating Scale and 583 ng/mL for dopamine D2 receptor occupancy. Both treatment duration and quetiapine exposure were associated with weight gain.

CONCLUSIONS

Concurrent administration of potent or moderate CYP3A4 inducers and inhibitors need to be avoided in quetiapine-treated patients. When co-medication is required, it is recommended to adjust the dosage based on therapeutic drug monitoring. Additional research is warranted to delineate the dose-exposure-response relationships of quetiapine and active metabolite norquetiapine in pediatrics, geriatrics, hepatically-impaired patients, and women using contraceptives or are pregnant or menopausal.

PROSPERO REGISTRATION

CRD42023446654.

External evaluation and systematic review of population pharmacokinetic models for high-dose methotrexate in cancer patients

Several population pharmacokinetic (PPK) models have been established to optimize the therapeutic regimen and reduce the toxicity of high-dose methotrexate (HDMTX) in patients with cancer. However, their predictive performance when extrapolated to different clinical centers was unknown. In this study, we aimed to externally evaluate the predictive ability of HDMTX PPK models and determine the potential influencing factors. We searched the literature and determined the predictive performance of the selected models using methotrexate concentrations in 721 samples from 60 patients in the First Affiliated Hospital of the Navy Medical University. Prediction-based diagnostics and simulation-based normalized prediction distribution errors (NPDE) were used to evaluate the predictive performance of the models. The influence of prior information was also assessed using Bayesian forecasting, and the potential factors affecting model predictability were investigated. Thirty models extracted from published PPK studies were assessed. Prediction-based diagnostics showed that the number of compartments potentially influenced model transferability, and simulation-based NPDE indicated model misspecification. Bayesian forecasting significantly improved the predictive performance of the models. Various factors, including bioassays, covariates, and population diagnosis, influence model extrapolation. The published models were unsatisfactory for all prediction-based diagnostics, except for the 24 h methotrexate concentration monitoring and simulation-based diagnostics, making them inappropriate for direct extrapolation. Moreover, Bayesian forecasting combined therapeutic drug monitoring could improve the predictive performance of the models.

Parametric population pharmacokinetics of isoniazid: a systematic review

INTRODUCTION

Isoniazid (INH) plays an important role in prevention and treatment of tuberculosis (TB). However, large pharmacokinetic (PK) variations are observed in patients receiving standard INH dosages. Considering the influence of PK variations on INH efficacy or adverse reactions, we reviewed the population PK studies of INH and explored significant covariates that influence INH PK.

METHODS

The PubMed and Embase databases were systematically searched from their inception to 30 January 2023. PPK studies on INH using a parametric nonlinear mixed-effect approach were included in this review. The characteristics and identified significant covariates of the included studies were summarized.

RESULTS

Twenty-one studies conducted in adults, and seven in pediatrics were included. A two-compartment model with first-order absorption and elimination was the frequently used structural model for INH. NAT2 genotype, body size, and age were identified as significant covariates affecting INH PK variation. The median clearance (CL) value in the fast metabolizers was 2.55-fold higher than that in the slow metabolizers. Infants and children had higher CL per weight values than adults with the same metabolic phenotype. In pediatric patients, CL value increased with postnatal age.

CONCLUSIONS

Compared with slow metabolizers, the daily dose of INH should be increased by 200-600 mg in fast metabolizers. To achieve effective treatment, pediatric patients need a higher dose per kilogram than adults. Further PPK studies of anti-tuberculosis drugs are needed to comprehensively understand the covariates that affect their PK characteristics and to achieve accurate dose adjustments.

Is a Lower Dose of Rivaroxaban Required for Asians? A Systematic Review of a Population Pharmacokinetics and Pharmacodynamics Analysis of Rivaroxaban

Rivaroxaban has been widely used to prevent and treat various thromboembolic diseases for more than a decade. However, whether a lower dose of rivaroxaban is required for Asians is still debatable. This review aimed to explore the potential ethnic difference in pharmacokinetic/pharmacodynamic (PK/PD) characteristics between Asians and Caucasians. A systematic search was conducted and twenty-four studies were identified, of which 10 were conducted on Asian adults, 11 on predominantly Caucasian adults, and 3 on Caucasian pediatrics. The apparent clearance (CL/F) of rivaroxaban in Caucasian adults with non-valvular atrial fibrillation (6.45-7.64 L/h) was about 31-43% higher than that in Asians (4.46-5.98 L/h) taking 10~20 mg rivaroxaban every 24 h. Moreover, there was no obvious difference in CL/F among Japanese, Chinese, Thai, and Irani people. Regarding PK/PD relationship, prothrombin time was linked to rivaroxaban concentration in a linear or near-linear manner, and Factor Xa activity was linked with the E_max model. The exposure-response relationship was comparable between Asians and Caucasians. Renal function has a significant influence on CL/F, and no covariate was recognized for exposure-response relationship. In conclusion, a lower dose of rivaroxaban might be required for Asians, and further studies are warranted to verify this ethnic difference to facilitate optimal dosing regimens.

EA, Udupa KS, Gota V, Mallayasamy S. Factors Influencing Pharmacokinetics of Tamoxifen in Breast Cancer Patients: A Systematic Review of Population Pharmacokinetic Models

BACKGROUND

Tamoxifen is useful in managing breast cancer and it is reported to have significant variability in its pharmacokinetics. This review aimed to summarize reported population pharmacokinetics studies of tamoxifen and to identify the factors affecting the pharmacokinetics of tamoxifen in adult breast cancer patients.

METHOD

A systematic search was undertaken in Scopus, Web of Science, and PubMed for papers published in the English language from inception to 20 August 2022. Studies were included in the review if the population pharmacokinetic modeling was based on non-linear mixed-effects modeling with a parametric approach for tamoxifen in breast cancer patients.

RESULTS

After initial selection, 671 records were taken for screening. A total of five studies were selected from Scopus, Web of Science, PubMed, and by manual searching. The majority of the studies were two-compartment models with first-order absorption and elimination to describe tamoxifen and its metabolites' disposition. The CYP2D6 phenotype and CYP3A4 genotype were the main covariates that affected the metabolism of tamoxifen and its metabolites. Other factors influencing the drug's pharmacokinetics included age, co-medication, BMI, medication adherence, CYP2B6, and CYP2C19 genotype.

CONCLUSION

The disposition of tamoxifen and its metabolites varies primarily due to the CYP2D6 phenotype and CYP3A4 genotype. However, other factors, such as anthropometric characteristics and menopausal status, should also be addressed when accounting for this variability. All these studies should be externally evaluated to assess their applicability in different populations and to use model-informed dosing in the clinical setting.

Ulcerative Colitis and Acute Severe Ulcerative Colitis Patients Are Overlooked in Infliximab Population Pharmacokinetic Models: Results from a Comprehensive Review

Ulcerative colitis (UC) is part of the inflammatory bowels diseases, and moderate to severe UC patients can be treated with anti-tumour necrosis α monoclonal antibodies, including infliximab (IFX). Even though treatment of UC patients by IFX has been in place for over a decade, many gaps in modelling of IFX PK in this population remain. This is even more true for acute severe UC (ASUC) patients for which early prediction of IFX pharmacokinetic (PK) could highly improve treatment outcome. Thus, this review aims to compile and analyse published population PK models of IFX in UC and ASUC patients, and to assess the current knowledge on disease activity impact on IFX PK. For this, a semi-systematic literature search was conducted, from which 26 publications including a population PK model analysis of UC patients receiving IFX therapy were selected. Amongst those, only four developed a model specifically for UC patients, and only three populations included severe UC patients. Investigations of disease activity impact on PK were reported in only 4 of the 14 models selected. In addition, the lack of reported model codes and assessment of predictive performance make the use of published models in a clinical setting challenging. Thus, more comprehensive investigation of PK in UC and ASUC is needed as well as more adequate reports on developed models and their evaluation in order to apply them in a clinical setting.

Influence of N-acetyltransferase 2 (NAT2) genotype/single nucleotide polymorphisms on clearance of isoniazid in tuberculosis patients: a systematic review of population pharmacokinetic models

Purpose

Significant pharmacokinetic variabilities have been reported for isoniazid across various populations. We aimed to summarize population pharmacokinetic studies of isoniazid in tuberculosis (TB) patients with a specific focus on the influence of N-acetyltransferase 2 (NAT2) genotype/single-nucleotide polymorphism (SNP) on clearance of isoniazid.

Methods

A systematic search was conducted in PubMed and Embase for articles published in the English language from inception till February 2022 to identify population pharmacokinetic (PopPK) studies of isoniazid. Studies were included if patient population had TB and received isoniazid therapy, non-linear mixed effects modelling, and parametric approach was used for building isoniazid PopPK model and NAT2 genotype/SNP was tested as a covariate for model development.

Results

A total of 12 articles were identified from PubMed, Embase, and hand searching of articles. Isoniazid disposition was described using a two-compartment model with first-order absorption and linear elimination in most of the studies. Significant covariates influencing the pharmacokinetics of isoniazid were NAT2 genotype, body weight, lean body weight, body mass index, fat-free mass, efavirenz, formulation, CD4 cell count, and gender. Majority of studies conducted in adult TB population have reported a twofold or threefold increase in isoniazid clearance for NAT2 rapid acetylators compared to slow acetylators.

Conclusion

The variability in disposition of isoniazid can be majorly attributed to NAT2 genotype. This results in a trimodal clearance pattern with a multi-fold increase in clearance of NAT2 rapid acetylators compared to slow acetylators. Further studies exploring the generalizability/adaptability of developed PopPK models in different clinical settings are required.

Application of a Physiologically Based Pharmacokinetic Model to Predict Cefazolin and Cefuroxime Disposition in Obese Pregnant Women Undergoing Caesarean Section

Intravenous (IV) cefuroxime and cefazolin are used prophylactically in caesarean sections (CS). Currently, there are concerns regarding sub-optimal dosing in obese pregnant women compared to lean pregnant women prior to CS. The current study used a physiologically based pharmacokinetic (PBPK) approach to predict cefazolin and cefuroxime pharmacokinetics in obese pregnant women at the time of CS as well as the duration that these drug concentrations remain above a target concentration (2, 4 or 8 µg/mL or µg/g) in plasma or adipose tissue. Cefazolin and cefuroxime PBPK models were first built using clinical data in lean and in obese non–pregnant populations. Models were then used to predict cefazolin and cefuroxime pharmacokinetics data in lean and obese pregnant populations. Both cefazolin and cefuroxime models sufficiently described their total and free levels in the plasma and in the adipose interstitial fluid (ISF) in non–pregnant and pregnant populations. The obese pregnant cefazolin model predicted adipose exposure adequately at different reference time points and indicated that an IV dose of 2000 mg can maintain unbound plasma and adipose ISF concentration above 8 µg/mL for 3.5 h post dose. Predictions indicated that an IV 1500 mg cefuroxime dose can achieve unbound plasma and unbound ISF cefuroxime concentration of ≥8 µg/mL up to 2 h post dose in obese pregnant women. Re-dosing should be considered if CS was not completed within 2 h post cefuroxime administration for both lean or obese pregnant if cefuroxime concentrations of ≥8 µg/mL is required. A clinical study to measure cefuroxime adipose concentration in pregnant and obese pregnant women is warranted.

Published population pharmacokinetic models of valproic acid in adult patients: a systematic review and external validation in a Chinese sample of inpatients with bipolar disorder

BACKGROUND

This study reviewed all published valproic acid (VPA) population pharmacokinetic (PPK) models in adult patients and assessed them using external validation methods to determine predictive performance.

METHODS

Thirteen published PPK models (labeled with letters A to M) not restricted to children were identified in PubMed, Embase, and Web of Science databases. They were evaluated in a sample totaling 411 serum concentrations from 146 adult inpatients diagnosed with bipolar disorder in a Chinese hospital. Serum concentrations of VPA were analyzed by validated ultra-performance liquid chromatography-tandem mass spectrometry. Performance was assessed by 4 tests (prediction-based diagnostics, visual predictive checks, normalized prediction distribution error, and Bayesian forecasting).

RESULTS

Models K and L, developed in large samples of Chinese and Thai patients, showed good performance in our Chinese dataset. Models H and J demonstrated good performance in Tests 2 and 3 of the 4 tests, respectively. Another 7 models exhibited intermediate performance. The models with the worst performance, F and M, could not be improved by Bayesian forecasting.

CONCLUSION

In our validation study the most important factors contributing to good performance were absence of children, Asian ethnicity, one-compartment models and inclusion of body weight and VPA dose in previously published models.

Population Pharmacokinetics of Vancomycin in Adult Patients with Long Bones’ Fractures

Vancomycin is a tricyclic glycopeptide antibiotic, mostly used in the treatment of severe staphylococcal and enterococcal infections, especially in orthopedic surgery.

The purpose of this analysis was to develop a population pharmacokinetic (PPK) model of vancomycine in hospitalized patients with bone fractures and identify important factors which influence its clearance (CL).

A total of ninety-nine measurements of vancomycin serum concentrations were used in our population modeling. A two-compartment model was applied to describe the pharmacokinetics of vancomycin using subroutines ADVAN3 and TRANS4.

The study population included patients of both sexes, with the mean age of 62.12±14.69 years and body weight of 80.32±12.44kg. Vancomycin was administered as intravenous infusion with average daily dose of 1772.73±521.34mg. Out of twenty different factors evaluated in the study (including demographic, clinical and laboratory data), only daily dose of vancomycin (DD) and co-medication with piperacillin/tazobactam (PT) showed significant effect on clearance of vancomycin. The final model was described by the following equation: CL (l/h) = 0.03 + 0.000468 x DD + 0.675 x PT. Bootstrapping was used for validation of the final model.

In conclusion, the main causes of variability in the clearance of vancomycin among adult patients with bone fractures are daily dose of vancomycin and co-medication with piperacillin/tazobactam.

Population pharmacokinetics of oxcarbazepine: a systematic review

INTRODUCTION

Oxcarbazepine is commonly used as first-line treatment for partial and generalized tonic-clonic seizures. Owing to the high pharmacokinetic variability, several population pharmacokinetic models have been developed for oxcarbazepine to explore potential covariates that affect its pharmacokinetic variation.

AREAS COVERED

This review summarizes the published population pharmacokinetic studies of oxcarbazepine in children and adults available in PubMed and Embase databases. The quality of the retrieved studies was evaluated, and significant covariates that may have an impact on the dosage regimen of oxcarbazepine were explored.

EXPERT OPINION

The pharmacokinetics of oxcarbazepine was founded to be affected by body weight and co-administration with enzyme inducers. Pediatric patients require a higher dose per kilogram than adults because children generally have a higher clearance than adults. Moreover, to maintain the target concentration, patients co-administrate with enzyme inducers need a higher dose than monotherapy due to higher clearance in those patients. Because limited information is available for exposure-response relationship, additional pharmacokinetic/pharmacodynamics investigations of oxcarbazepine need to be conducted to optimize the dosage regimen in clinical practice.

Can Population Pharmacokinetics of Antibiotics be Extrapolated? Implications of External Evaluations

BACKGROUND AND OBJECTIVE

External evaluation is an important issue in the population pharmacokinetic analysis of antibiotics. The purpose of this review was to summarize the current approaches and status of external evaluations and discuss the implications of external evaluation results for the future individualization of dosing regimens.

METHODS

We systematically searched the PubMed and EMBASE databases for external evaluation studies of population analysis and extracted the relevant information from these articles. A total of 32 studies were included in this review.

RESULTS

Vancomycin was investigated in 17 (53.1%) articles and was the most studied drug. Other studied drugs included gentamicin, tobramycin, amikacin, amoxicillin, ceftaroline, meropenem, fluconazole, voriconazole, and rifampicin. Nine (28.1%) studies were prospective, and the sample size varied widely between studies. Thirteen (40.6%) studies evaluated the population pharmacokinetic models by systematically searching for previous studies. Seven (21.9%) studies were multicenter studies, and 27 (84.4%) adopted the sparse sampling strategy. Almost all external evaluation studies of antibiotics (93.8%) used metrics for prediction-based diagnostics, while relatively fewer studies were based on simulations (46.9%) and Bayesian forecasting (25.0%).

CONCLUSION

The results of external evaluations in previous studies revealed the poor extrapolation performance of existing models of prediction- and simulation-based diagnostics, whereas the posterior Bayesian method could improve predictive performance. There is an urgent need for the development of standards and guidelines for external evaluation studies.

What "Impact" Do NLME Publications Have Outside Our Community?

The number of publications applying nonlinear mixed‐effect (NLME) modeling has increased yearly since its first appearance in 1979. Here, we evaluated articles that have used NLME modeling, were published in journals that attract a broader audience, and we discussed the standard of presentation of these to stimulate target audience‐specific improvements for increased impact in the future.

Modelling pharmacokinetics and pharmacodynamics of the selective S1P1 receptor modulator cenerimod in healthy subjects and systemic lupus erythematosus patients

AIMS

Assessment of time to attain steady state as well as drug accumulation following long-term treatment with the selective sphingosine-1-phosphate 1 receptor modulator cenerimod and prediction of the incidence of low total lymphocyte (LY) counts. Differences in pharmacokinetics and pharmacodynamics based on demographic characteristics and between healthy subjects and systemic lupus erythematosus (SLE) patients were to be identified.

METHODS

Data from 4 Phase I studies and 1 Phase II study were pooled to develop a population pharmacokinetic/pharmacodynamic model describing cenerimod concentration and its effect on LY count. Simulations addressed the objectives.

RESULTS

Simulations of 365 days of treatment indicated a time to steady state of 49 days and changes in exposure of 15% beyond 35 days. For a dose of 2 mg, the predicted incidences of LY counts below 0.5 and 0.2 × 10⁹ cells/L were 18.2 and 0.6% for healthy subjects and 25.9 and 1.0% for SLE patients, respectively. Incidence increased with higher dose and lower baseline LY counts. For body weights of 50 and 100 kg compared to 75 kg, exposure was predicted to change by +37% and -20%, respectively.

CONCLUSION

Long-term cenerimod administration is not expected to result in exposure and LY count reduction substantially different from results of completed studies. Low LY counts are predicted to occur more frequently in SLE patients compared to healthy subjects. Dose individualization based on the model is not considered necessary. Model-based simulations enable benefit-risk evaluations, supporting planning of late-phase clinical studies and scientific exchange with health authorities.

Population Pharmacokinetic and Exposure-Response Analysis of Finerenone: Insights Based on Phase IIb Data and Simulations to Support Dose Selection for Pivotal Trials in Type 2 Diabetes with Chronic Kidney Disease

BACKGROUND

Finerenone (BAY 94-8862) is a potent non-steroidal, selective mineralocorticoid receptor antagonist being developed for the treatment of patients with type 2 diabetes and chronic kidney disease.

METHODS

We present the population pharmacokinetics and pharmacodynamics (PD) analysis for efficacy and safety markers based on data from two clinical phase IIb studies: ARTS-DN (NCT01874431) and ARTS-DN Japan (NCT01968668).

RESULTS

The pharmacokinetics of finerenone were adequately characterized, with estimated glomerular filtration rate (eGFR) and body weight as influencing covariates. The area under the plasma concentration-time curve in Japanese patients did not differ from that in the global population, and the investigated pharmacokinetics were dose- and time-linear. In addition, the pharmacokinetic model provided robust individual exposure estimates to study exposure-response. The concentration-effect relationship over time for the efficacy marker urinary albumin:creatinine ratio (UACR) was well-characterized by a maximum effect model indicating saturation at high exposures. For the safety markers, a log-linear model and a power model were identified for serum potassium concentration and eGFR, respectively, indicating attenuation of effect gains at high exposures. There was no apparent ethnic effect on the investigated pharmacokinetic-pharmacodynamic relationships. The model-predicted times to reach the full (99%) steady-state drug effect on UACR, serum potassium, and eGFR were 138, 20, and 85 days, respectively, while the pharmacokinetic half-life was 2-3 h and steady state was achieved after 2 days, indicating timescale separation.

CONCLUSION

Our dose-exposure-response modeling and simulation indicates effects were largely saturated at finerenone 20 mg and doses of both 10 and 20 mg once daily appear safe and efficacious at reducing albuminuria.

Evaluation of Assumptions Underpinning Pharmacometric Models

Assumptions inherent to pharmacometric model development and use are not routinely acknowledged, described, or evaluated. The aim of this work is to present a framework for systematic evaluation of assumptions. To aid identification of assumptions, we categorise assumptions into two types: implicit and explicit assumptions. Implicit assumptions are inherent in a method or model and underpin its derivation and use. Explicit assumptions arise from heuristic principles and are typically defined by the user to enable the application of a method or model. A flowchart was developed for systematic evaluation of assumptions. For each assumption, the impact of assumption violation ('significant', 'insignificant', 'unknown') and the probability of assumption violation ('likely', 'unlikely', 'unknown') will be evaluated based on prior knowledge or the result of an additional bespoke study to arrive at a decision ('go', 'no-go') for both model building and model use. A table of assumptions with standardised headings has been proposed to facilitate the documentation of assumptions and evaluation of results. The utility of the proposed framework was illustrated using four assumptions underpinning a top-down model describing the warfarin-coagulation proteins' relationship. The next step of this work is to apply the framework to a series of other settings to fully assess its practicality and its value in identifying and making inferences from assumptions.

Reproducible pharmacokinetics

Reproducibility is a highly desired feature of scientific investigation in general, and it has special connotations for research in pharmacokinetics, a vibrant field with over 500,000 publications to-date. It is important to be able to differentiate between genuine heterogeneity in pharmacokinetic parameters from heterogeneity that is due to errors and biases. This overview discusses efforts and opportunities to diminish the latter type of undesirable heterogeneity. Several reporting and research guidance documents and standards have been proposed for pharmacokinetic studies, but their adoption is still rather limited. Quality problems in the methods used and model evaluations have been examined in some empirical studies of the literature. Standardization of statistical and laboratory tools and procedures can be improved in the field. Only a small fraction of pharmacokinetic studies become pre-registered and only 9995 such studies have been registered in ClinicalTrials.gov as of August 2018. It is likely that most pharmacokinetic studies remain unpublished. Publication bias affecting the results and inferences has been documented in case studies, but its exact extent is unknown for the field at-large. The use of meta-analyses in the field is still limited. Availability of raw data, detailed protocols, software and codes is hopefully improving with multiple ongoing initiatives. Several research practices can contribute to greater transparency and reproducibility for pharmacokinetic investigations.

Population Pharmacokinetics of Tacrolimus in Transplant Recipients: What Did We Learn About Sources of Interindividual Variabilities?

Tacrolimus, a calcineurin inhibitor, is a common immunosuppressant prescribed after organ transplantation and has notable inter- and intrapatient pharmacokinetic variability. The sources of variability have been investigated using population pharmacokinetic modeling over the last 2 decades. This article provides an updated synopsis on published nonlinear mixed-effects analyses developed for tacrolimus in transplant recipients. The objectives were to establish a detailed overview of the current data and to investigate covariate relationships determined by the models. Sixty-three published analyses were reviewed, and data regarding the study design, modeling approach, and resulting findings were extracted and summarized. Most of the studies investigated tacrolimus pharmacokinetics in adult and pediatric renal and liver transplants after administration of the immediate-release formulation. Model structures largely depended on the study sampling strategy, with ∼50% of studies developing a 1-compartment model using trough concentrations and a 2-compartment model with delayed absorption from intensive sampling. The CYP3A5 genotype, as a covariate, consistently impacted tacrolimus clearance, and dosing adjustments were required to achieve similar drug exposure among patients. Numerous covariates were identified as sources of interindividual variability on tacrolimus pharmacokinetics with limited consistency across these studies, which may be the result of the study designs. Additional analyses are required to further evaluate the potential impact of these covariates and the clinical implementation of these models to guide tacrolimus dosing recommendations. This article may be useful for guiding the design of future population pharmacokinetic studies and provides recommendations for the selection of an existing optimal model to individualize tacrolimus therapy.

An Evaluation of Vancomycin Area Under the Curve Estimation Methods for Children Treated for Acute Pulmonary Exacerbations of Cystic Fibrosis Due to Methicillin-Resistant Staphylococcus aureus

The prevalence of pulmonary methicillin-resistant Staphylococcus aureus infections in patients with cystic fibrosis (CF) has increased over the last 2 decades. Two concentrations-a postdistributive and a trough-are currently used to estimate the area under the curve (AUC) of vancomycin, an antibiotic routinely used to treat these infections, to achieve the target AUC/minimum inhibitory concentration of ≥400 mg·h/L in ensuring optimal dosing of this drug. This study evaluated precision and bias in estimating vancomycin AUCs obtained either from a population pharmacokinetic (PK) model by using a single trough concentration or from standard PK equation-based 2-point monitoring approach. AUCs were either obtained from a single trough concentration-fitted model or derived from a model fitted by 2 concentration points. Children ≥2 years of age with CF received intravenous vancomycin at 2 centers from June 2012 to December 2014. A population PK model was developed in Pmetrics to quantify the between-subject variability in vancomycin PK parameters, define the sources of PK variability, and leverage information from the population to improve individual AUC estimates. Twenty-three children with CF received 27 courses of vancomycin. The median age was 12.3 (interquartile range [IQR] 8.5-16.6) years. From the individual vancomycin PK parameter estimates from the population PK model, median AUC was 622 (IQR 529-680) mg·h/L. Values were not significantly different from the AUC calculated using the standard PK equation-based approach (median 616 [IQR 540-663] mg·h/L) (P = .89). A standard PK equation-based approach using 2 concentrations and a population PK model-based approach using a single trough concentration yielded unbiased and precise AUC estimates. Findings suggest that options exist to implement AUC-based pediatric vancomycin dosing in patients with CF. The findings of this study reveal that several excellent options exist for centers to implement AUC-based pediatric vancomycin dosing for patients with CF.

Pharmacokinetics of sildenafil in children with pulmonary arterial hypertension

BACKGROUND

Recently, sildenafil was introduced to treat pulmonary arterial hypertension (PAH); however, there are currently few studies on the pharmacokinetics of sildenalfil in children. Therefore, we aimed to carry out a pharmacokinetic study of sildenafil in children with PAH using a single dose.

METHODS

Twelve children diagnosed with PAH, consisting of with ten males and two females, were recruited for the study after obtaining written consent from their parents or guardians. Blood samples were obtained predose and at 0.25, 0.5, 1, 2, 4, 8 and 12 hours after the oral administration of 1 mg/kg of sildenafil using an extemporal pediatric formulation developed in our laboratory. The samples were analyzed using a previously validated high performance liquid chromatography method.

RESULTS

A pharmacokinetic analysis using the WinNonlin 3.1 program that considered the Akaike information criterion (AIC) for selecting a more adjustable model was performed. The following pharmacokinetic parameters were obtained: maximal concentration (C_max): 366±179 ng/mL, time to maximal concentration: 0.92±0.30 hours, elimination half-life (t_1/2): 2.41±1.18 hours, total clearance (CL_tot/F): 5.85±2.81 L/hour, volume of distribution (Vd/F): 20.13±14.5 L, absorption rate constants (Ka): 0.343 hour^-1, elimination rate (Ke): 0.35 hour^-1, area under curve from zero to infinity: 2061±618 ng/mL/hour. The data of all patients adjusted to the model of one compartment were corroborated using AIC.

CONCLUSIONS

The parameters Ka, Ke and t_1/2 were found to be similar to those reported in adults; however, the values of C_max and Vd/F were significantly higher. Based on these findings, we propose that treatment regimen of sildenafil be adjusted in children with PAH.

Prolonged infusion of sedatives and analgesics in adult intensive care patients: A systematic review of pharmacokinetic data reporting and quality of evidence

Although pharmacokinetic (PK) data for prolonged sedative and analgesic agents in intensive care unit (ICU) has been described, the number of publications in this important area appear relatively few, and PK data presented is not comprehensive. Known pathophysiological changes in critically ill patients result in altered drug PK when compared with non-critically ill patients. ClinPK Statement was recently developed to promote consistent reporting in PK studies, however, its applicability to ICU specific PK studies is unclear. In this systematic review, we assessed the overall ClinPK Statement compliance rate, determined the factors affecting compliance rate, graded the level of PK evidence and assessed the applicability of the ClinPK Statement to future ICU PK studies. Of the 33 included studies (n=2016), 22 (67%) were low evidence quality descriptive studies (Level 4). Included studies had a median compliance rate of 80% (IQR 66% to 86%) against the ClinPK Statement. Overall pooled compliance rate (78%, 95% CI 73% to 83%) was stable across time (P=0.38), with higher compliance rates found in studies fitting three compartments models (88%, P<0.01), two compartments models (83%, P<0.01) and one compartment models (77%, P=0.17) than studies fitting noncompartmental or unspecified models (69%) (P<0.01). Data unique to the interpretation of PK data in critically ill patients, such as illness severity (48%), organ dysfunction (36%) and renal replacement therapy use (32%), were infrequently reported. Discrepancy between the general compliance rate with ClinPK Statement and the under-reporting of ICU specific parameters suggests that the applicability of the ClinPK Statement to ICU PK studies may be limited in its current form.

Reporting guidelines for population pharmacokinetic analyses

The purpose of this work was to develop a consolidated set of guiding principles for the reporting of population pharmacokinetic (PK) analyses based on input from a survey of practitioners as well as discussions between industry, consulting, and regulatory scientists. The survey found that identification of population covariate effects on drug exposure and support for dose selection (in which population PK frequently serves as preparatory analysis for exposure-response modeling) are the main areas of influence for population PK analysis. The proposed guidelines consider 2 main purposes of population PK reports: (1) to present key analysis findings and their impact on drug development decisions, and (2) as documentation of the analysis methods for the dual purpose of enabling review of the analysis and facilitating future use of the models. This work also identified 2 main audiences for the reports: (1) a technically competent group responsible for in-depth review of the data, methodology, and results; and (2) a scientifically literate but not technically adept group, whose main interest is in the implications of the analysis for the broader drug development program. We recommend a generalized question-based approach with 6 questions that need to be addressed throughout the report. We recommend 8 sections (Synopsis, Introduction, Data, Methods, Results, Discussion, Conclusions, Appendix) with suggestions for the target audience and level of detail for each section. A section providing general expectations regarding population PK reporting from a regulatory perspective is also included. We consider this an important step toward industrialization of the field of pharmacometrics such that a nontechnical audience also understands the role of pharmacometric analyses in decision making. Population PK reports were chosen as representative reports to derive these recommendations; however, the guiding principles presented here are applicable for all pharmacometric reports including pharmacokinetics/pharmacodynamics and simulation reports.

Reporting guidelines for population pharmacokinetic analyses

The purpose of this work was to develop a consolidated set of guiding principles for reporting of population pharmacokinetic (PK) analyses based on input from a survey of practitioners as well as discussions between industry, consulting and regulatory scientists. The survey found that identification of population covariate effects on drug exposure and support for dose selection (where population PK frequently serves as preparatory analysis to exposure–response modeling) are the main areas of influence for population PK analysis. The proposed guidelines consider two main purposes of population PK reports (1) to present key analysis findings and their impact on drug development decisions, and (2) as documentation of the analysis methods for the dual purpose of enabling review of the analysis and facilitating future use of the models. This work also identified two main audiences for the reports: (1) a technically competent group responsible for in-depth review of the data, methodology, and results, and (2) a scientifically literate, but not technically adept group, whose main interest is in the implications of the analysis for the broader drug development program. We recommend a generalized question-based approach with six questions that need to be addressed throughout the report. We recommend eight sections (Synopsis, Introduction, Data, Methods, Results, Discussion, Conclusions, Appendix) with suggestions for the target audience and level of detail for each section. A section providing general expectations regarding population PK reporting from a regulatory perspective is also included. We consider this an important step towards industrialization of the field of pharmacometrics such that non-technical audience also understands the role of pharmacometrics analyses in decision making. Population PK reports were chosen as representative reports to derive these recommendations; however, the guiding principles presented here are applicable for all pharmacometric reports including PKPD and simulation reports.