Steroid-tacrolimus drug-drug interaction and the effect of CYP3A genotypes.

AIMS: Tacrolimus, metabolized by CYP3A4 and CYP3A5 enzymes, is susceptible to drug-drug interactions (DDI). Steroids induce CYP3A genes to increase tacrolimus clearance, but the effect is variable. We hypothesized that the extent of the steroid-tacrolimus DDI differs by CYP3A4/5 genotypes. METHODS: Kidney transplant recipients (n = 2462) were classified by the number of loss of function alleles (LOF) (CYP3A5*3, *6 and *7 and CYP3A4*22) and steroid use at each tacrolimus trough in the first 6 months post-transplant. A population pharmacokinetic analysis was performed by nonlinear mixed-effect modelling (NONMEM) and stepwise covariate modelling to define significant covariates affecting tacrolimus clearance. A stochastic simulation was performed and translated into a Shiny application with the mrgsolve and Shiny packages in R. RESULTS: Steroids were associated with modestly higher (3%-11.8%) tacrolimus clearance. Patients with 0-LOF alleles receiving steroids showed the greatest increase (11.8%) in clearance compared to no steroids, whereas those with 2-LOFs had a negligible increase (2.6%) in the presence of steroids. Steroid use increased tacrolimus clearance by 5% and 10.3% in patients with 1-LOF and 3/4-LOFs, respectively. CONCLUSIONS: Steroids increase the clearance of tacrolimus but vary slightly by CYP3A genotype. This is important in individuals of African ancestry who are more likely to carry no LOF alleles, may more commonly receive steroid treatment, and will need higher tacrolimus doses.

Limited Sampling Strategies Fail to Accurately Predict Mycophenolic Acid Area Under the Curve in Kidney Transplant Recipients and the Impact of Enterohepatic Recirculation.

BACKGROUND: Therapeutic drug monitoring for mycophenolic acid (MPA) is challenging due to difficulties in measuring the area under the curve (AUC). Limited sampling strategies (LSSs) have been developed for MPA therapeutic drug monitoring but come with risk of unacceptable performance. The authors hypothesized that the poor predictive performance of LSSs were due to the variability in MPA enterohepatic recirculation (EHR). This study is the first to evaluate LSSs models performance in the context of EHR. METHODS: Adult kidney transplant recipients (n = 84) receiving oral mycophenolate mofetil underwent intensive MPA pharmacokinetic sampling. MPA AUC0-12hr and EHR were determined. Published MPA LSSs in kidney transplant recipients receiving tacrolimus were evaluated for their predictive performance in estimating AUC0-12hr in our full cohort and separately in individuals with high and low EHR. RESULTS: None of the evaluated LSS models (n = 12) showed good precision or accuracy in predicting MPA AUC0-12hr in the full cohort. In the high EHR group, models with late timepoints had better accuracy but low precision, except for 1 model with late timepoints at 6 and 10 hours postdose, which had marginally acceptable precision. For all models, the good guess of predicted AUC0-12hr (±15% of observed AUC0-12hr) was highly variable (range, full cohort = 19%-61.9%; high EHR = 4.5%-65.9%; low EHR = 27.5%-62.5%). CONCLUSIONS: The predictive performance of the LSS models varied according to EHR status. Timepoints ≥5 hours postdose in LSS models are essential to capture EHR. Models and strategies that incorporate EHR during development are required to accurately ascertain MPA exposure.

Population Pharmacokinetic Analysis of Atomoxetine and its Metabolites in Children and Adolescents with Attention-Deficit/Hyperactivity Disorder.

Atomoxetine (ATX) is a non-stimulant used to treat attention-deficit/hyperactivity disorder (ADHD) and systemic exposure is highly variable due to polymorphic cytochrome P450 2D6 (CYP2D6) activity. The objective of this study was to characterize the time course of ATX and metabolites (4-hydroxyatomoxetine (4-OH); N-desmethylatomoxetine (NDA); and 2-carboxymethylatomoxetine (2-COOH)) exposure following oral ATX dosing in children with ADHD to support individualized dosing. A nonlinear mixed-effect modeling approach was used to analyze ATX, 4-OH, and NDA plasma and urine, and 2-COOH urine profiles obtained over 24-72 hours from children with ADHD (n = 23) following a single oral ATX dose. Demographics and CYP2D6 activity score (AS) were evaluated as covariates. Simulations were performed to explore the ATX dosing in subjects with various CYP2D6 AS. A simultaneous pharmacokinetic modeling approach was used in which a model for ATX, 4-OH, and NDA in plasma and urine, and 2-COOH in urine was developed. Plasma ATX, 4-OH, and NDA were modeled using two-compartment models with first-order elimination. CYP2D6 AS was a significant determinant of ATX apparent oral clearance (CL/F), fraction metabolized to 4-OH, and systemic exposure of NDA. CL/F of ATX varied almost 7-fold across the CYP2D6 AS groups: AS 2: 20.02 L/hour; AS 1: 19.00 L/hour; AS 0.5: 7.47 L/hour; and AS 0: 3.10 L/hour. The developed model closely captures observed ATX, 4-OH, and NDA plasma and urine, and 2-COOH urine profiles. Application of the model shows the potential for AS-based dosing recommendations for improved individualized dosing.

Population Pharmacokinetic Model of Intravenous Busulfan in Hematopoietic Cell Transplantation: Systematic Review and Comparative Simulations.

BACKGROUND: Busulfan is commonly used in the chemotherapy prior to hematopoietic cell transplantation (HCT). Busulfan has a narrow therapeutic window and a well-established exposure-response relationship with important clinical outcomes. Model-informed precision dosing (MIPD) based on population pharmacokinetic (popPK) models has been implemented in the clinical settings. We aimed to systematically review existing literature on popPK models of intravenous busulfan. METHODS: We systematically searched Ovid MEDLINE, EMBASE, Cochrane Library, Scopus, and Web of Science databases from inception to December 2022 to identify original popPK models (nonlinear mixed-effect modeling) of intravenous busulfan in HCT population. Model-predicted busulfan clearance (CL) was compared using US population data. RESULTS: Of the 44 eligible popPK studies published since 2002, 68% were developed predominantly in children, 20% in adults, and 11% in both children and adults. The majority of the models were described using first-order elimination or time-varying CL (69% and 26%, respectively). All but three included a body-size descriptor (e.g., body weight, body surface area). Other commonly included covariates were age (30%) and GSTA1 variant (15%). Median between-subject and between-occasion variabilities of CL were 20% and 11%, respectively. Between-model variabilities in predicted median CL were < 20% in all of the weight tiers (10-110 kg) in the simulation based on US population data. CONCLUSION: Busulfan PK is commonly described using a first-order elimination or time-varying CL. A simple model with limited covariates were generally sufficient to attain relatively small unexplained variabilities. However, therapeutic drug monitoring may still be necessary to attain a narrow target exposure.

Higher abatacept exposure after transplant decreases acute GVHD risk without increasing adverse events.

In the ABA2 study, the T-cell costimulation blockade agent, abatacept, was safe and effective in preventing acute graft-versus-host disease (aGVHD) after unrelated-donor hematopoietic cell transplant (HCT), leading to US Food and Drug Administration approval. Here, we performed a determination of abatacept pharmacokinetics (PK), which enabled an examination of how abatacept exposure-response relationships affected clinical outcomes. We performed a population PK analysis of IV abatacept using nonlinear mixed-effect modeling and assessed the association between abatacept exposure and key transplant outcomes. We tested the association between the trough after dose 1 (Ctrough_1) and grade (GR) 2 or 4 aGVHD (GR2-4 aGVHD) through day +100. An optimal Ctrough_1 threshold was identified via recursive partitioning and classification tree analysis. This demonstrated that abatacept PK was characterized by a 2-compartment model with first-order elimination. The ABA2 dosing regimen was based on previous work targeting a steady-state abatacept trough of 10 µg/mL. However, a higher Ctrough_1 (≥39 µg/mL, attained in ∼60% of patients on ABA2) was associated with a favorable GR2-4 aGVHD risk (hazard ratio, 0.35; 95% confidence interval, 0.19-0.65; P < .001), with a Ctrough_1 <39 µg/mL associated with GR2-4 aGVHD risk indistinguishable from placebo (P = .37). Importantly, no significant association was found between Ctrough_1 and key safety indicators, including relapse, and cytomegalovirus or Epstein-Barr virus viremia. These data demonstrate that a higher abatacept Ctrough_1 (≥39 µg/mL) was associated with a favorable GR2-4 aGVHD risk, without any observed exposure-toxicity relationships. This trial was registered at www.clinicaltrials.gov as #NCT01743131.

Influence of Renal Function on Phosphoramide Mustard Exposure: A Nonlinear Mixed-Effects Analysis.

Phosphoramide mustard (PM) is the final cytotoxic metabolite formed from the parent compound cyclophosphamide through a complex metabolic pathway, primarily through hepatic metabolism. Little is known about the effect of renal elimination on the disposition of PM. We evaluated the effect of renal function on PM exposure after single doses of cyclophosphamide in 85 patients undergoing allogeneic hematopoietic cell transplantation using nonlinear mixed-effects modeling. Mixed linear and nonlinear elimination pathways were required to adequately describe the disposition of PM. Creatinine clearance (CrCL) was incorporated as a covariate associated with first-order elimination, representing renal clearance (ClR ) of PM. For a 70-kg patient, ClR was 14.9 L/h, Volume of distribution was 525 L, maximum rate was 81.2 mg/h, and the concentration to achieve 50% of maximum rate was 0.51 mg/L. We conducted simulations to explore the impact of CrCL as a measure of renal function and observed that when CrCL decreases from 120 to 40 mL/min, PM area under the plasma concentration-time curve (AUC) from time 0 to 8 hours and AUC increases by 9.2% and 80.9% on average after a single dose, respectively. Our data suggest that renal function has limited influence on PM exposure during the first 8 hours after dosing but has a large impact on the total exposure. Dose adjustment of cyclophosphamide may not be necessary in hematopoietic cell transplant recipients with moderate to severe kidney dysfunction to attain targeted exposures based on AUC from time 0 to 8 hours. However, dose reduction may be necessary if demonstrated at some future time that total AUC is a better surrogate for safety or toxicity.

Comparison of Dose Adjustment Strategies for Obesity in High-dose Cyclophosphamide Among Adult Hematopoietic Cell Transplantation Recipients: Pharmacokinetic Analysis.

Cyclophosphamide (CY) is an alkylating agent widely used in the field of oncology and hematopoietic cell transplantation (HCT). It is recommended to use an adjusted body weight with an adjustment factor of 0.25 (ABW25) for dosing of CY in obese patients undergoing HCT. However, evidence based on the pharmacokinetics (PK) of CY to support this recommendation is lacking. We aimed to identify a dosing strategy of CY that achieves equivalent exposures among obese and nonobese patients. The present study is a secondary analysis of a previously conducted observational PK study of phosphoramide mustard (PM), the final cytotoxic metabolite of CY. Data were collected from 85 adults with hematologic malignancy who received a single infusion of CY 50 mg/kg, fludarabine, ± anti-thymocyte globulin, and a single fraction of total body irradiation as HCT conditioning therapy. A previously developed population PK model in these patients was used for simulations. Using individualized PK parameters from that analysis, simulations were performed to assess cumulative exposures of PM (i.e., area-under-the-curve [AUC]) resulting from 8 different dosing strategies according to various measures of body size: (1) "mg/kg" by total body weight (TBW); (2) "mg/kg" by ideal body weight (IBW); (3) "mg/kg" by fat free mass; (4) "mg/m2" by body surface area (BSA); (5) "mg/kg" by TBW combined with ABW25 (TBW-ABW25); (6) "mg/kg" by IBW combined with ABW25 (IBW-ABW25); (7) "mg/kg" by TBW combined with ABW by adjustment factor of 0.50 (TBW-ABW50); and (8) "mg" by fixed-dose. We defined equivalent exposure as the effect of obesity on PM AUC within ±20% from the PM AUC in the nonobese group, where obesity is defined based on TBW/IBW ratio (i.e., nonobese, <1.2; mildly obese, 1.2-1.5; and moderately/severely obese, >1.5). Primary and secondary outcomes were PM AUC0-8hours and PM AUC0-infinity, respectively. In the 85 patients, with the median age of 63 years (range 21-75), 46% were classified as mildly and 25% were moderately/severely obese based on the TBW/IBW ratio. Negative correlations (i.e., higher the extent of obesity, lower the PM AUC) were shown when dosing simulations were based on IBW, TBW-ABW25, and fixed dosing (P < .05). Positive correlations were shown when dosing was simulated by TBW (P < .05). None of the 8 dosing strategies attained equivalent PM AUC0-8hours between patients with versus without obesity, whereas dosing by BSA and TBW-ABW50 attained equivalent PM AUC0-infinity (P < .05). Our study predicted that the recommended ABW25 dose adjustment may result in lower exposure of CY therapy in obese patients than in nonobese. A CY dosing strategy that would result in similar PM concentrations between obese and nonobese was not identified for early exposure (i.e., PM AUC0-8hours). The data suggest though that CY dosing based on "mg/m2" by BSA or "mg/kg" by TBW-ABW50 would result in similar total exposure (i.e., PM AUC0-infinity) and may minimize exposure differences in obese and nonobese patients.

Identification of a novel fusion Iduronidase with improved activity in the cardiovascular system.

Background: Lysosomal diseases are a group of over 70 rare genetic conditions in which a protein deficiency (most often an enzyme deficiency) leads to multi-system disease. Current therapies for lysosomal diseases are limited in their ability to treat certain tissues that are major contributors to morbidity and mortality, such as the central nervous system (CNS) and cardiac valves. For this study, the lysosomal disease mucopolysaccharidosis type I (MPS I) was selected as the disease model. In MPS I, mutations in the IDUA gene cause a deficiency of the α-L-iduronidase (IDUA) enzyme activity, leading to disease pathology in tissues throughout the body, including the CNS and cardiac valves. Current therapies have been unable to prevent neurodevelopmental deficits and cardiac valvular disease in patients with MPS I. This study aimed to evaluate the delivery of IDUA enzyme, via a novel gene therapy construct, to target tissues. Methods: MPS I mice were hydrodynamically injected through the tail vein with plasmids containing either a codon-optimized cDNA encoding the wild-type IDUA protein or one of four modified IDUAs under the control of the liver-specific human α1-antitrypsin (hAAT) promoter. Two modified IDUAs contained a ligand for the CB1 receptor, which is a highly expressed receptor in the CNS. Iduronidase activity levels were measured in the tissues and plasma using an enzyme activity assay. Results: The modified IDUAs did not appear to have improved activity levels in the brain compared with the unmodified IDUA. However, one modified IDUA exhibited higher activity levels than the unmodified IDUA in the heart (p = 0.0211). This modified iduronidase (LT-IDUA) contained a sequence for a six amino acid peptide termed LT. LT-IDUA was further characterized using a noncompartmental pharmacokinetic approach that directly analyzed enzyme activity levels after gene delivery. LT-IDUA had a 2-fold higher area under the curve (AUC) than the unmodified IDUA (p = 0.0034) when AUC was estimated using enzyme activity levels in the plasma. Conclusion: The addition of a six amino acid peptide improved iduronidase's activity levels in the heart and plasma. The short length of this LT peptide facilitates its use as fusion enzymes encoded as gene therapy or administered as enzyme replacement therapy. More broadly, the LT peptide may aid in developing therapies for numerous lysosomal diseases.

Physiologically-Based Pharmacokinetic Model Development, Validation, and Application for Prediction of Eliglustat Drug-Drug Interactions.

Eliglustat is a glucosylceramide synthase inhibitor indicated as a long-term substrate reduction therapy for adults with type 1 Gaucher disease, a lysosomal rare disease. It is primarily metabolized by cytochrome P450 2D6 (CYP2D6), and variants in the gene encoding this enzyme are important determinants of eliglustat pharmacokinetics (PK) and drug-drug interactions (DDIs). The existing drug label addresses the DDIs to some extent but has omitted scenarios where both metabolizing CYPs (2D6 and 3A4) are mildly or moderately inhibited. The objectives of this study were (i) to develop and validate an eliglustat physiologically-based pharmacokinetic (PBPK) model with and without drug interactions, (ii) to simulate untested DDI scenarios, and (iii) to explore potential dosing flexibility using lower dose strength of eliglustat (commercially not available). PK data from healthy adults receiving eliglustat with or without interacting drugs were obtained from literature and used for the PBPK model development and validation. The model-predicted single-dose and steady-state maximum concentration (Cmax ) and area under the concentration-time curve (AUC) of eliglustat were within 50-150% of the observed values when eliglustat was administered alone or coadministered with ketoconazole or paroxetine. Then as model-based simulations, we illustrated eliglustat exposure as a victim of interaction when coadministered with fluvoxamine following the US Food and Drug Administration (FDA) dosing recommendations. Second, we showed that with lower eliglustat doses (21 mg, 42 mg once daily) the exposure in participants of intermediate and poor metabolizer phenotypes was within the outlined safety margin (Cmax <250 ng/mL) when eliglustat was administered with ketoconazole, where the current recommendation is a contraindication of coadministration (84 mg). The present study demonstrated that patients with CYP2D6 deficiency may benefit from lower doses of eliglustat.

Reduced Enterohepatic Recirculation of Mycophenolate and Lower Blood Concentrations Are Associated with the Stool Bacterial Microbiome after Hematopoietic Cell Transplantation.

Mycophenolate mofetil (MMF) is an important immunosuppressant used after allogeneic hematopoietic cell transplantation (HCT). MMF has a narrow therapeutic index, and blood concentrations of mycophenolic acid (MPA), the active component of MMF, are highly variable. Low MPA concentrations are associated with the risk of graft-versus-host disease (GVHD), whereas high concentrations are associated with toxicity. Reasons for variability are not well known and may include the presence of ß-glucuronidase-producing bacteria in the gastrointestinal tract, which enhance MPA enterohepatic recirculation (EHR) by transforming MPA metabolites formed in the liver back to MPA. This study was conducted to determine whether individuals with high MPA EHR have a greater abundance of ß-glucuronidase-producing bacteria in their stool and higher MPA concentrations compared with those with low EHR. We conducted a pharmacomicrobiomics study in 20 adult HCT recipients receiving a myeloablative or reduced-intensity preparative regimen. Participants received MMF 1 g i.v. every 8 hours with tacrolimus. Intensive pharmacokinetic sampling of MMF was conducted before hospital discharge; total MPA, MPA glucuronide (MPAG), and acyl-glucuronide metabolite (acylMPAG) were measured. EHR was defined as the ratio of MPA area under the concentration-versus-time curve (AUC)4-8 to MPA AUC0-8. Differences in stool microbiome diversity and composition, determined by shotgun metagenomic sequencing, were compared above and below the median EHR (22%; range, 5% to 44%). The median EHR was 12% in the low EHR group and 29% in the high EHR group. MPA troughs, MPA AUC4-8, and acyl-glucuronide metabolite (acylMPAG) AUC4-8/AUC0-8 ratio were greater in the high EHR group compared with the low EHR group (1.53 µg/mL versus .28 µg/mL [P = .0001], 7.33 hour·µg/mL versus 1.79 hour·µg/mL [P = .0003], and .33 hour·µg/mL versus .24 hour·µg/mL [P = .0007], respectively). MPA AUC0-8 was greater in the high EHR group than in the low EHR group, and the difference trended toward significance (22.8 hour·µg/mL versus 15.3 hour·µg/mL; P = .06). Bacteroides vulgatus, Bacteroides stercoris, and Bacteroides thetaiotaomicron were 1.2- to 2.4-fold more abundant (P = .039, .024, and .046, respectively) in the high EHR group. MPA EHR was positively correlated with B. vulgatus (⍴ = .58; P ≤ .01) and B. thetaiotaomicron (⍴ = .46; P < .05) and negatively correlated with Blautia hydrogenotrophica (⍴ = -.53; P < .05). Therapeutic MPA troughs were achieved in 80% of patients in the high EHR group but in no patients in the low EHR group. There was a trend toward differences in MPA AUC0-8 and MPA concentration at steady-state (µg/mL) between the high EHR group versus the low EHR group (P = .06). MPA EHR was variable. Patients with high MPA EHR had greater abundance of Bacteroides species in stool and higher MPA exposure compared with patients with low MPA EHR. Therefore, Bacteroides may be protective against poor outcomes, such as graft-versus-host disease, in some patients but may increase the risk of MPA adverse effects in others. These data need to be confirmed and studied after oral MMF therapy.

Precision Dosing for Tacrolimus Using Genotypes and Clinical Factors in Kidney Transplant Recipients of European Ancestry.

Genetic variation in the CYP3A4 and CYP3A5 (CYP3A4/5) genes, which encode the key enzymes in tacrolimus metabolism, is associated with tacrolimus clearance and dose requirements. Tacrolimus has a narrow therapeutic index with high intra- and intersubject variability, in part because of genetic variation. High tacrolimus clearance and low trough concentration are associated with a greater risk for rejection, whereas high troughs are associated with calcineurin-induced toxicity. The objective of this study was to develop a model of tacrolimus clearance with a dosing equation accounting for genotypes and clinical factors in adult kidney transplant recipients of European ancestry that could preemptively guide dosing. Recipients receiving immediate-release tacrolimus for maintenance immunosuppression from 2 multicenter studies were included. Participants in the GEN03 study were used for tacrolimus model development (n = 608 recipients) and was validated by prediction performance in the DeKAF Genomics study (n = 1361 recipients). Nonlinear mixed-effects modeling was used to develop the apparent oral tacrolimus clearance (CL/F) model. CYP3A4/5 genotypes and clinical covariates were tested for their influence on CL/F. The predictive performance of the model was determined by assessing the bias (median prediction error [ME] and median percentage error [MPE]) and the precision (root median squared error [RMSE]) of the model. CYP3A5*3, CYP3A4*22, corticosteroids, calcium channel blocker and antiviral drug use, age, and diabetes significantly contributed to the interindividual variability of oral tacrolimus apparent clearance. The bias (ME, MPE) and precision (RMSE) of the final model was good, 0.49 ng/mL, 6.5%, and 3.09 ng/mL, respectively. Prospective testing of this equation is warranted.

Pharmacokinetic-Pharmacodynamic Modeling of Tumor Targeted Drug Delivery Using Nano-Engineered Mesenchymal Stem Cells.

Nano-engineered mesenchymal stem cells (nano-MSCs) are promising targeted drug delivery platforms for treating solid tumors. MSCs engineered with paclitaxel (PTX) loaded poly(lactide-co-glycolide) (PLGA) nanoparticles (NPs) are efficacious in treating lung and ovarian tumors in mouse models. The quantitative description of pharmacokinetics (PK) and pharmacodynamics (PD) of nano-MSCs is crucial for optimizing their therapeutic efficacy and clinical translatability. However, successful translation of nano-MSCs is challenging due to their complex composition and physiological mechanisms regulating their pharmacokinetic-pharmacodynamic relationship (PK-PD). Therefore, in this study, a mechanism-based preclinical PK-PD model was developed to characterize the PK-PD relationship of nano-MSCs in orthotopic A549 human lung tumors in SCID Beige mice. The developed model leveraged literature information on diffusivity and permeability of PTX and PLGA NPs, PTX release from PLGA NPs, exocytosis of NPs from MSCs as well as PK and PD profiles of nano-MSCs from previous in vitro and in vivo studies. The developed PK-PD model closely captured the reported tumor growth in animals receiving no treatment, PTX solution, PTX-PLGA NPs and nano-MSCs. Model simulations suggest that increasing the dosage of nano-MSCs and/or reducing the rate of PTX-PLGA NPs exocytosis from MSCs could result in improved anti-tumor efficacy in preclinical settings.

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.

Pharmacokinetic/pharmacodynamic data extrapolation models for improved pediatric efficacy and toxicity estimation, with application to secondary hyperparathyroidism.

In most drug development settings, the regulatory approval process is accompanied by extensive studies performed to understand the drug's pharmacokinetic (PK) and pharmacodynamic (PD) properties. In this article, we attempt to utilize the rich PK/PD data to inform the borrowing of information from adults during pediatric drug development. In pediatric settings, it is especially crucial that we are parsimonious with the patients recruited for experimentation. We illustrate our approaches in the context of clinical trials of cinacalcet for treating secondary hyperparathyroidism in pediatric and adult patients with chronic kidney disease, where we model both parathyroid hormone (efficacy endpoint) and corrected calcium levels (safety endpoint). We use population PK/PD modeling of the cinacalcet data to quantitatively assess the similarity between adults and children, and use this information in various hierarchical Bayesian adult borrowing rules whose statistical properties can then be evaluated. In particular, we simulate the bias and mean square error performance of our approaches in settings where borrowing is and is not warranted to inform guidelines for the future use of our methods.

Bupivacaine Pharmacokinetics and Breast Milk Excretion of Liposomal Bupivacaine Administered After Cesarean Birth.

OBJECTIVE: To evaluate bupivacaine concentrations in maternal plasma and transfer into breast milk in women undergoing liposomal bupivacaine infiltration in the transversus abdominis plane after cesarean birth. METHODS: Prospective cohort study of healthy pregnant women who underwent cesarean birth at term followed by a transversus abdominis plane block using 52 mg bupivacaine hydrochloride 0.25% (20 mL) and 266 mg liposomal bupivacaine 1.3% (20 mL). Simultaneous blood and milk samples were collected in a staggered fashion, three to four samples per patient at the following timepoints after block administration: 2, 6, 12, 24, 48, 72, and 96 hours. Quantification of bupivacaine was performed by liquid chromatography-tandem mass spectrometry. Neonatal drug exposure was modeled by calculating milk/plasma area under the curve (AUC) ratios, neonatal dosage, and relative neonatal dosage of bupivacaine at each sampling time. RESULTS: Thirty patients were enrolled. Concentrations in breast milk peaked at 6 hours (mean 58 ng/mL), followed by constant and steady decline to low levels at 96 hours (mean 5.2 ng/mL). Maternal plasma concentrations had two peaks, first at 6 hours (mean 155.9 ng/mL) and then at 48 hours (mean 225.8 ng/mL), followed by steady decline. Milk/plasma AUC0-t ratios ranged between AUC0-2 of 0.45 (80% CI 0.38-0.52) and AUC0-96 of 0.15 (80% CI 0.14-0.17). Neonatal dosage ranged between a mean of 355.9 ng/kg at 0-2 hours and a mean of 15,155.4 ng/kg at 0-96 hours. Relative neonatal dosage was less than 1% at all time intervals. No serious adverse reactions occurred in any neonate. CONCLUSION: Bupivacaine is excreted in breast milk after local infiltration of liposomal bupivacaine and bupivacaine hydrochloride mixture into transversus abdominis plane blocks after cesarean birth. Relative neonatal dosages of less than 1% (less than 10% is considered to be unlikely to be of clinical concern) suggest minimal risks for breastfeeding healthy, term neonates after the administration of this combination of local anesthetics to mothers. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov, NCT03526419.

Finding the Dose for Hydroxychloroquine Prophylaxis for COVID-19: The Desperate Search for Effectiveness.

Hydroxychloroquine is an antimalarial drug being tested as a potential treatment for the novel coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2. Although the efficacy of hydroxychloroquine for COVID-19 remains uncertain, it may serve as a potential prophylactic agent especially in those at high risk, such as healthcare workers, household contacts of infected patients, and the immunocompromised. Our aim was to identify possible hydroxychloroquine dosing regimens through simulation in those at high risk of infections by optimizing exposures above the in vitro generated half maximal effective concentration (EC50 ) and to help guide researchers in dose-selection for COVID-19 prophylactic studies. To maintain weekly troughs above EC50 in > 50% of subjects at steady-state in a pre-exposure prophylaxis setting, an 800 mg loading dose followed by 400 mg twice or 3 times weekly is required. In an exposure driven, post-exposure prophylaxis setting, 800 mg loading dose followed in 6 hours by 600 mg, then 600 mg daily for 4 more days achieved daily troughs above EC50 in > 50% subjects. These doses are higher than recommended for malaria chemoprophylaxis, and clinical trials are needed to establish safety and efficacy.

IL-1ß reduces cardiac lymphatic muscle contraction via COX-2 and PGE2 induction: Potential role in myocarditis.

The role of lymphatic vessels in myocarditis is largely unknown, while it has been shown to play a key role in other inflammatory diseases. We aimed to investigate the role of lymphatic vessels in myocarditis using in vivo model induced with Theiler's murine encephalomyelitis virus (TMEV) and in vitro model with rat cardiac lymphatic muscle cells (RCLMC). In the TMEV model, we found that upregulation of a set of inflammatory mediator genes, including interleukin (IL)-1ß, tumor necrosis factor (TNF)-αand COX-2 were associated with disease activity. Thus, using in vitro collagen gel contraction assays, we decided to clarify the role(s) of these mediators by testing contractility of RCLMC in response to IL-1ß and TNF-α individually and in combination, in the presence or absence of: IL-1 receptor antagonist (Anakinra); cyclooxygenase (COX) inhibitors inhibitors (TFAP, diclofenac and DuP-697). IL-1ß impaired RCLMC contractility dose-dependently, while co-incubation with both IL-1ß and TNF-α exhibited synergistic effects in decreasing RCLMC contractility with increased COX-2 expression. Anakinra maintained RCLMC contractility; Anakinra blocked the mobilization of COX-2 induced by IL-1ß with or without TNF-α. COX-2 inhibition blocked the IL-1ß-mediated decrease in RCLMC contractility. Mechanistically, we found that IL-1ß increased prostaglandin (PG) E2 release dose-dependently, while Anakinra blocked IL-1ß mediated PGE2 release. Using prostaglandin E receptor 4 (EP4) receptor antagonist, we demonstrated that EP4 receptor blockade maintained RCLMC contractility following IL-1ß exposure. Our results indicate that IL-1ß reduces RCLMC contractility via COX-2/PGE2 signaling with synergistic cooperation by TNF-α. These pathways may help provoke inflammatory mediator accumulation within the heart, driving progression from acute myocarditis into dilated cardiomyopathy.