Covariate modeling in pharmacometrics: General points for consideration.

Modeling the relationships between covariates and pharmacometric model parameters is a central feature of pharmacometric analyses. The information obtained from covariate modeling may be used for dose selection, dose individualization, or the planning of clinical studies in different population subgroups. The pharmacometric literature has amassed a diverse, complex, and evolving collection of methodologies and interpretive guidance related to covariate modeling. With the number and complexity of technologies increasing, a need for an overview of the state of the art has emerged. In this article the International Society of Pharmacometrics (ISoP) Standards and Best Practices Committee presents perspectives on best practices for planning, executing, reporting, and interpreting covariate analyses to guide pharmacometrics decision making in academic, industry, and regulatory settings.

Lanadelumab in Patients 2 to Less Than 12 Years Old With Hereditary Angioedema: Results From the Phase 3 SPRING Study.

BACKGROUND: Symptoms of hereditary angioedema (HAE) often first occur during childhood, and HAE attacks in children can be severe and substantially affect health-related quality of life (HRQoL). However, there are no approved long-term prophylaxis treatments for children aged less than 6 years. OBJECTIVE: The SPRING Study (NCT04070326) evaluated the safety, pharmacokinetics, and efficacy of lanadelumab and HRQoL in patients aged 2 to less than 12 years. METHODS: Over 52 weeks of treatment, patients aged 2 to less than 6 years received lanadelumab 150 mg every 4 weeks (Q4W) and patients aged 6 to less than 12 years received 150 mg every 2 weeks (Q2W) but could switch to Q4W if they were attack-free for 26 weeks. RESULTS: We enrolled 21 patients (aged 2 to less than 6 years: n = 4; aged 6 to less than 12 years: n = 17), 20 of whom completed the study. There were no reported serious treatment-emergent adverse events or discontinuations resulting from such events. Treatment-emergent adverse events were reported for 17 patients (81.0%). The most common TEAE was injection site pain. Overall systemic exposure was comparable for both age groups. The mean (SD) attack rate during treatment decreased by 94.8% from baseline (1.84 [1.53] to 0.08 [0.17] attacks/mo), and 16 (76.2%) patients were attack-free. The attack rate reduction in both age groups was similar during the first 26-week fixed-dosing treatment. Seven patients switched from Q2W to Q4W and remained attack-free. A large, clinically meaningful increase in the Pediatric Quality of Life Inventory Generic Core Scale Total Score and a large increase in the Pediatric Quality of Life Inventory Generic Core Scale-Family Impact Module Total Score from baseline to end of study (better HRQoL) were observed. CONCLUSIONS: Findings support safety, efficacy, and improved HRQoL with lanadelumab 150 mg Q2W and Q4W regimens for the prevention of HAE attacks in patients aged 2 to less than 12 years.

Model-Informed Approaches and Innovative Clinical Trial Design for Adeno-Associated Viral Vector-Based Gene Therapy Product Development: A White Paper.

The promise of viral vector-based gene therapy (GT) as a transformative paradigm for treating severely debilitating and life-threatening diseases is slowly coming to fruition with the recent approval of several drug products. However, they have a unique mechanism of action often necessitating a tortuous clinical development plan. Expertise in such complex therapeutic modality is still fairly limited in this emerging class of adeno-associated virus (AAV) vector-based gene therapies. Because of the irreversible mode of action and incomplete understanding of genotype-phenotype relationship and disease progression in rare diseases careful considerations should be given to GT product's benefit-risk profile. In particular, special attention needs to be paid to safe dose selection, reliable dose exposure response (including clinically relevant endpoints), or creative approaches in study design targeting small patient populations during clinical development. We believe that quantitative tools encompassed within model-informed drug development (MIDD) framework fits quite well in the development of such novel therapies, as they enable us to benefit from the totality of data approach in order to support dose selection as well as optimize clinical trial designs, end point selection, and patient enrichment. In this thought leadership paper, we provide our collective experiences, identify challenges, and suggest areas of improvement in applications of modeling and innovative trial design in development of AAV-based GT products and reflect on the challenges and opportunities for incorporating MIDD tools and more in rational development of these products.

Achieving big with small: quantitative clinical pharmacology tools for drug development in pediatric rare diseases.

Pediatric populations represent a major fraction of rare diseases and compound the intrinsic challenges of pediatric drug development and drug development for rare diseases. The intertwined complexities of pediatric and rare disease populations impose unique challenges to clinical pharmacologists and require integration of novel clinical pharmacology and quantitative tools to overcome multiple hurdles during the discovery and development of new therapies. Drug development strategies for pediatric rare diseases continue to evolve to meet the inherent challenges and produce new medicines. Advances in quantitative clinical pharmacology research have been a key component in advancing pediatric rare disease research to accelerate drug development and inform regulatory decisions. This article will discuss the evolution of the regulatory landscape in pediatric rare diseases, the challenges encountered during the design of rare disease drug development programs and will highlight the use of innovative tools and potential solutions for future development programs.

Clinical Pharmacology in Drug Development for Rare Diseases in Neurology: Contributions and Opportunities.

Several challenges are associated with rare disease drug development in neurology. In this article, we summarize the US Food and Drug Administration's experience with clinical drug development for rare neurological diseases and discuss clinical pharmacology's critical contributions to drug development for rare diseases. We used publicly available information to identify and screen drug products approved for rare neurological indications between 1983 and 2019. We highlighted cases in which clinical pharmacology contributed to the evidence of drug efficacy, dose selection for pivotal clinical trials, dose optimization based on intrinsic and extrinsic factors, pharmacokinetic bridging for formulations, and efficacy bridging across different racial groups. Fifty-one approved drug products were identified since the introduction of the Orphan Drug Act in 1983. Interestingly, the number of approvals in the last few years increased significantly, probably due to advances in genomic research and targeted drug modalities. Evaluation of dose selection in patient populations showed that in 52% of cases, the sponsors did not evaluate efficacy for more than one or two dose levels throughout the development program. Clinical pharmacology studies to evaluate the effect of intrinsic or extrinsic factors were adequately characterized in most of the applications. With the expansion of model informed drug development applications, (e.g., quantitative systems pharmacology and deep learning neural network models), the role and impact of clinical pharmacology is expected to grow exponentially in the next decade and enhance the development of novel treatment modalities for neurological rare diseases.

Concordance of Exposure Changes Because of Renal Impairment Between Results of Dedicated Renal Impairment Studies and Population Pharmacokinetic Predictions.

This analysis compared the results from noncompartmental analysis and population pharmacokinetic (PopPK) predictions of exposure changes in patients with renal impairment (RI) for 27 new molecular entities (NMEs) approved between 2000 and 2015. Renal function was identified as a covariate in the final PopPK model for 17 NMEs. The final PopPK model was used to simulate (n = 1000 replicates/individual) the results of a dedicated PK study in subjects with renal impairment. For the majority of NMEs, concordance between observed, and predicted area under the curve (AUC) geometric mean ratio (GMR) was observed (ie, in 17, 11, and 11 NMEs for mild, moderate, and severe renal impairment groups, respectively, the observed and predicted AUC GMR were within the same fold of change). Inclusion of colinear covariates in the PopPK model appeared to be the major driver for the NMEs for which there was discordance. PopPK, when done properly, is a valuable tool for supporting labeling recommendations for subjects with renal impairment.

An integrated PK-PD model for cortisol and the 17-hydroxyprogesterone and androstenedione biomarkers in children with congenital adrenal hyperplasia.

AIMS: The aim of this study was to characterize the pharmacokinetic/pharmacodynamic relationships of cortisol and the adrenal biomarkers 17-hydroxyprogesterone and androstenedione in children with congenital adrenal hyperplasia (CAH). METHODS: A nonlinear mixed-effect modelling approach was used to analyse cortisol, 17-hydroxyprogesterone and androstenedione concentrations obtained over 6 hours from children with CAH (n = 50). A circadian rhythm was evident and the model leveraged literature information on circadian rhythm in untreated children with CAH. Indirect response models were applied in which cortisol inhibited the production rate of all three compounds using an Imax model. RESULTS: Cortisol was characterized by a one-compartment model with apparent clearance and volume of distribution estimated at 22.9 L/h/70 kg and 41.1 L/70 kg, respectively. The IC50 values of cortisol concentrations for cortisol, 17-hydroxyprogesterone and androstenedione were estimated to be 1.36, 0.45 and 0.75 µg/dL, respectively. The inhibitory effect was found to be more potent on 17OHP than D4A, and the IC50 values were higher in salt-wasting subjects than simple virilizers. Production rates of cortisol, 17-hydroxyprogesterone and androstenedione were higher in simple-virilizer subjects. Half-lives of cortisol, 17-hydroxyprogesterone and androstenedione were 60, 47 and 77 minutes, respectively. CONCLUSION: Rapidly changing biomarker responses to cortisol concentrations highlight that single measurements provide volatile information about a child's disease control. Our model closely captured observed cortisol, 17-hydroxyprogesterone and androstenedione concentrations. It can be used to predict concentrations over 24 hours and allows many novel exposure metrics to be calculated, e.g., AUC, AUC-above-threshold, time-within-range, etc. Our long-range goal is to uncover dose-exposure-outcome relationships that clinicians can use in adjusting hydrocortisone dose and timing.

Application of Clinical Pharmacology Principles in Drug Development of Modified-Release Products: Leveraging Exposure-Response Information to Support Approval.

The development of modified-release (MR) drug products aims to address a clinical need such as improving patient compliance. There are multiple pathways and development strategies for the registration and approval of MR products. The development strategy of an MR product is usually dependent on the availability and pharmacokinetic/pharmacodynamics (PK/PD) characteristics of the reference drug product, that is, an immediate-release (IR) product or a reference MR. Compared with a reference IR product, an MR product is likely to have a different PK profile over the least common dosing time due to unequal dosing intervals. In case of differences in PK profiles between the MR product and the reference product, confirmatory efficacy and safety studies may be needed to support registration. In some cases, however, a thorough clinical PK/PD characterization may provide sufficient basis to support the approval of the proposed MR product without the need for additional safety and efficacy studies. This article summarizes the US Food and Drug Administration experience and the regulatory considerations supporting the approval of MR products in the past 6 years and discusses cases in which clinical pharmacology and PK/PD information were leveraged to support approval without the need for additional clinical studies. Details of all these cases are available in the public domain. In 2 cases a well-characterized exposure-response relationship provided sufficient justification that differences in the shape of the PK profiles were not clinically relevant. In the remaining 3 cases a thorough characterization of the PK profile along with a risk-based approach provided bases for approval.

Pivotal Considerations for Optimal Deployment of Healthy Volunteers in Oncology Drug Development.

Oncology drug development is among the most challenging of any therapeutic area, with first-in-human trials expected to deliver information on both safety and activity. Until recently, therapeutic approaches in oncology focused on cytotoxic chemotherapy agents, ruling out even the possibility of enrolling normal healthy volunteers (NHVs) in clinical trials due to safety considerations. The emergence of noncytotoxic modalities, including molecularly targeted agents with more favorable safety profiles, however, has led to increasing numbers of clinical pharmacology studies of these agents being conducted in NHVs. Beyond rapid enrollment and cost savings, there are other advantages of conducting specific types of studies in NHVs with the goal of more appropriate dosing decisions in certain subsets of the intended patient populations, allowing for enrollment of such patients in therapeutic trials from which they might otherwise have been excluded. Nevertheless, the decision must be carefully weighed against potential disadvantages, and although the considerations surrounding conduct of clinical trials using NHVs are generally well-defined in most other therapeutic areas, they are less well-defined in oncology.

Orphan drug development: the increasing role of clinical pharmacology.

Over the last few decades there has been a paradigm shift in orphan drug research and development. The development of the regulatory framework, establishment of rare disease global networks that support drug developments, and advances in technology, has resulted in tremendous growth in orphan drug development. Nevertheless, several challenges during orphan drug development such as economic constraints; insufficient clinical information; fewer patients and thus inadequate power; etc. still exist. While the standard regulatory requirements for drug approval stays the same, applications of scientific judgment and regulatory flexibility is significantly important to help meeting some of the immense unmet medical need in rare diseases. Clinical pharmacology presents a vital role in accelerating orphan drug development and overcoming some of these challenges. This review highlights the critical contributions of clinical pharmacology in orphan drug development; for example, dose finding, optimizing clinical trial design, indication expansion, and population extrapolation. Examples of such applications are reviewed in this article.

Mycophenolic Acid and Its Metabolites in Kidney Transplant Recipients: A Semimechanistic Enterohepatic Circulation Model to Improve Estimating Exposure.

Mycophenolic acid (MPA) is an approved immunosuppressive agent widely prescribed to prevent rejection after kidney transplantation. Wide between-subject variability (BSV) in MPA exposure exists which in part may be due to variability in enterohepatic recirculation (EHC). Several modeling strategies were developed to evaluate EHC as part of MPA pharmacokinetics, however mechanistic representation of EHC is limited. These models have not provided a satisfactory representation of the physiology of EHC in their modeling assumptions. The aim of this study was i) to develop an integrated model of MPA (total and unbound) and its metabolites (MPAG and acyl-MPAG) in kidney recipients, where this model provides a more physiological representation of EHC process, and ii) to evaluate the effect of donor and recipient clinical covariates and genotypes on MPA disposition. A five-compartment model with first-order input into an unbound MPA compartment connected to the MPAG, acyl-MPAG, and gallbladder compartment best fit the data. To represent the EHC process, the model was built based on the physiological concepts related to the hepatobiliary system and the gallbladder filling and emptying processes. The effect of cyclosporine versus tacrolimus on clearance of unbound MPA was included in the base model. Covariate analysis showed creatinine clearance to be significant on oral clearance of unbound MPA. The hepatic nuclear factor 1 alpha (HNF1A) genetic single nucleotide polymorphism (SNP) (rs2393791) in the recipient significantly affected the fraction of enterohepatically-circulated drug. Oral clearance of MPAG was affected by recipient IMPDH1 SNP (rs2288553), diabetes at the time of transplant, and donor sex.

Amikacin Pharmacokinetic-Pharmacodynamic Analysis in Pediatric Cancer Patients.

We performed pharmacokinetic-pharmacodynamic (PK-PD) and simulation analyses to evaluate a standard amikacin dose of 15 mg/kg once daily in children with cancer and to determine an optimal dosing strategy. A population pharmacokinetic model was developed from clinical data collected in 34 pediatric patients and used in a simulation study to predict the population probability of various dosing regimens to achieve accepted safety (steady-state unbound trough plasma concentration [fCmin] of <10 mg/liter)- and efficacy (free, unbound plasma concentration-to-MIC ratio [fCmax/MIC] of ≥8)-linked targets. In addition, an adaptive resistance PD (ARPD) model of Pseudomonas aeruginosa was built based on literature time-kill curve data and linked to the PK model to perform PK-ARPD simulations and compare results with those of the probability approach. Using the probability approach, an amikacin dose of 60 mg/kg administered once daily is expected to achieve the target fCmax/MIC in 80% of pediatric patients weighing 8 to 70 kg with a 97.5% probability, and almost all patients were predicted to have fCmin of <10 mg/liter. However, PK-ARPD simulation predicted that 60 mg/kg/day is unlikely to suppress bacterial resistance with repeated dosing. Furthermore, PK-ARPD simulation suggested that amikacin at 90 mg/kg, given in two divided doses (45 mg/kg twice a day), is expected to hit safety and efficacy targets and is associated with a lower rate of bacterial resistance. The disagreement between the two methods is due to the inability of the probability approach to predict development of drug resistance with repeated dosing. This originates from the use of PK-PD indices based on the MIC that neglects measurement errors, ignores the time course dynamic nature of bacterial growth and killing, and incorrectly assumes the MIC to be constant during treatment.

Stable Isotope Pharmacokinetic Studies Provide Insight into Effects of Age, Sex, and Weight on Levothyroxine Metabolism.

BACKGROUND: This study sought to determine whether levothyroxine pharmacokinetics (PKs) are affected by age, weight, and sex. METHODS: A PK study was performed after administration of a tracer dose of carbon-13-labeled LT4 (13C-LT4). The study was conducted at an academic medical center. Adults of any age being treated with levothyroxine for hypothyroidism were enrolled in the study. A single dose of 13C-LT4 was administered. Eighteen serial plasma samples were collected. One sample was obtained before the 13C-LT4 dose, and the majority of the remaining samples were collected over the 120-hour period post dosing. 13C-LT4 concentration was quantified using liquid chromatography tandem mass spectrometry. PK analysis was conducted using a linear log trapezoidal non-compartmental analysis using Phoenix 6.4. RESULTS: Eight males and 33 females with a median age of 50 years (range 22-78 years) and median weight of 65.9 kg (range 50-150 kg) were enrolled in the study. The median 13C-LT4 dose administered was 100 µg (range 70-300 µg). The median oral clearance rate (CL/F), apparent volume of distribution (V/F), time to peak concentration (Tmax), and dose-normalized peak concentration (Cmax) of 13C-LT4 were estimated to be 0.712 L/h, 164.9 L, 4 h, and 7.5 ng/L/µg, respectively. The dose-normalized area under the concentration-time curve from time 0 to 120 hours and half-life of the terminal distribution phase were 0.931 ng.h/mL/µg and 172.2 h, respectively. There was no significant difference in any 13C-LT4 PK parameter between patients aged >60 years (n = 10) and patients aged ≤60 years (n = 31), nor was there a relationship between age as a continuous variable and 13C-LT4 PK parameters. Sex only affected CL/F, V/F, and dose-normalized Cmax in univariate analyses. However, after adjusting for weight, sex was no longer a significant covariate. Weight was a significant predictor for CL/F, V/F and dose-normalized Cmax of 13C-LT4 in multivariate analyses. CONCLUSION: Prior studies suggest that patient age affects levothyroxine dose requirement. This study did not identify an effect of age and suggests that age-related changes in levothyroxine pharmacokinetics may be mediated by age-related weight differences. Physicians should consider a patient's weight, rather than age, for estimating levothyroxine dosage requirement.

A hierarchical Bayesian approach for combining pharmacokinetic/pharmacodynamic modeling and Phase IIa trial design in orphan drugs: Treating adrenoleukodystrophy with Lorenzo's oil.

X-linked adrenoleukodystrophy (X-ALD) is a rare, progressive, and typically fatal neurodegenerative disease. Lorenzo's oil (LO) is one of the few X-ALD treatments available, but little has been done to establish its clinical efficacy or indications for its use. In this article, we analyze data on 116 male asymptomatic pediatric patients who were administered LO. We offer a hierarchical Bayesian statistical approach to understand LO pharmacokinetics (PK) and pharmacodynamics (PD) resulting from an accumulation of very long-chain fatty acids. We experiment with individual- and observational-level errors and various choices of prior distributions and deal with the limitation of having just one observation per administration of the drug, as opposed to the more usual multiple observations per administration. We link LO dose to the plasma erucic acid concentrations by PK modeling, and then link this concentration to a biomarker (C26, a very long-chain fatty acid) by PD modeling. Next, we design a Bayesian Phase IIa study to estimate precisely what improvements in the biomarker can arise from various LO doses while simultaneously modeling a binary toxicity endpoint. Our Bayesian adaptive algorithm emerges as reasonably robust and efficient while still retaining good classical (frequentist) operating characteristics. Future work looks toward using the results of this trial to design a Phase III study linking LO dose to actual improvements in health status, as measured by the appearance of brain lesions observed via magnetic resonance imaging.

A model-based approach to assess the exposure-response relationship of Lorenzo's oil in adrenoleukodystrophy.

AIMS: X-linked adrenoleukodystrophy (X-ALD) is a peroxisomal disorder, most commonly affecting boys, associated with increased very long chain fatty acids (C26:0) in all tissues, causing cerebral demyelination and adrenocortical insufficiency. Certain monounsaturated long chain fatty acids including oleic and erucic acids, known as Lorenzo's oil (LO), lower plasma C26:0 levels. The aims of this study were to characterize the effect of LO administration on plasma C26:0 concentrations and to determine whether there is an association between plasma concentrations of erucic acid or C26:0 and the likelihood of developing brain MRI abnormalities in asymptomatic boys. METHODS: Non-linear mixed effects modelling was performed on 2384 samples collected during an open label single arm trial. The subjects (n = 104) were administered LO daily at ~2-3 mg kg(-1) with a mean follow-up of 4.88 ± 2.76 years. The effect of erucic acid exposure on plasma C26:0 concentrations was characterized by an inhibitory fractional Emax model. A Weibull model was used to characterize the time-to-developing MRI abnormality. RESULTS: The population estimate for the fractional maximum reduction of C26:0 plasma concentrations was 0.76 (bootstrap 95% CI 0.73, 0.793). Our time-to-event analyses showed that every mg l(-1) increase in time-weighted average of erucic acid and C26:0 plasma concentrations was, respectively, associated with a 3.7% reduction and a 753% increase in the hazard of developing MRI abnormality. However, the results were not significant (P = 0.5344, 0.1509, respectively). CONCLUSIONS: LO administration significantly reduces the abnormally high plasma C26:0 concentrations in X-ALD patients. Further studies to evaluate the effect of LO on the likelihood of developing brain MRI abnormality are warranted.