Pharmacokinetic-pharmacodynamic modelling and simulation of methotrexate dosing in patients with rheumatoid arthritis.

AIMS: To develop a non-linear mixed-effects population pharmacokinetic and pharmacodynamic (PK-PD) model describing the change in the concentration of methotrexate polyglutamates in erythrocytes (ery-MTX-PGn with "n" number of glutamate, representing PK component) and how this relates to modified 28-joint Disease Activity Score incorporating erythrocyte sedimentation rate (DAS-28-3) for rheumatoid arthritis (RA), representing PD component. METHODS: An existing PK model was fitted to data from a study consisting of 117 RA patients. The estimation of population PK-PD parameters was performed using stochastic approximation expectation maximisation algorithm in Monolix 2021R2. The model was used to perform Monte Carlo simulations of a loading dose regimen (50mg subcutaneous methotrexate as loading doses, then 20mg weekly oral methotrexate) compared to a standard dosing regimen (10mg weekly oral methotrexate for 2 weeks, then 20mg weekly oral methotrexate). RESULTS: Every 40 nmol/L increase in ery-MTX-PG3-5 total concentration correlated with 1-unit reduction in DAS-28-3. Significant covariate effects on the therapeutic response of methotrexate included the use of prednisolone in the first 4 weeks (positive use correlated with 25% reduction in DAS-28-3 when other variables were constant) and patient age (every 10-year increase in age correlated with 3.4% increase in DAS-28-3 when other variables were constant). 4 methotrexate loading doses led to a higher percentage of patients achieving a good/moderate response compared to the standard regimen (Week 4: 87.6% vs. 39.8%; Week 10: 64.7% vs. 57.0%). CONCLUSIONS: A loading dose regimen was more likely to achieve higher ery-MTX-PG concentration and better therapeutic response after 4 weeks of methotrexate treatment.

Spatio-temporal trends in livestock exposure to per- and polyfluoroalkyl substances (PFAS) inform risk assessment and management measures.

The migration of per- and polyfluoroalkyl substances (PFAS) onto agricultural properties has resulted in the accumulation of PFAS in livestock. The environmental determinants of PFAS accumulation in livestock from the grazing environment are poorly understood, resulting in limited capacity to manage livestock exposure and subsequent transfer of PFAS through the food chain. Analytical- (n = 978 samples of soil, water, pasture, and serum matrices), farm management/practice- and livestock physiology data were collated and interrogated from environmental PFAS investigations across ten farms, from four agro-ecological regions of Victoria (Australia). Statistical analysis identified perfluorooctane sulfonate (PFOS) and perfluorohexane sulfonate (PFHxS) as key analytes of concern for livestock bioaccumulation. PFOS and PFHxS concentrations in livestock drinking water were positively correlated with serum concentrations while other intake pathways (pasture and soil) had weaker correlations. Seasonal trends in PFAS body burden (serum concentrations) were identified and suggested to be linked to seasonal grazing behaviours and physiological water requirements. The data showed for the first time that livestock exposure to PFAS is dynamic and with relatively short elimination half-lives, there is opportunity for exposure management. Meat from cattle, grazed on PFAS impacted sites, may exceed health-based guideline values for PFAS, especially for markets with low limits (like the European Commission Maximum Limits or EC MLs). This study found that sites with mean livestock drinking water concentrations as low as 0.003 µg PFOS/L may exceed the EC ML for PFOS in cattle meat. Risk assessment can be used to prioritise site cleanup and development of management plans to reduce PFAS body burden by considering timing of stock rotation and/or supplementation of primary exposure sources.

Population pharmacokinetics of lenalidomide in patients with B-cell malignancies.

AIMS: Lenalidomide is an immunomodulatory imide drug used broadly in the treatment of multiple myeloma and lymphoma. It continues to be evaluated in chronic lymphocytic leukaemia (CLL) at lower doses due to dose-related toxicities including tumour flare and tumour lysis syndrome. This study aimed to develop a population pharmacokinetic model for lenalidomide in multiple cancers, including CLL, to identify any disease-related differences in disposition. METHODS: Lenalidomide concentrations from 4 clinical trials were collated (1999 samples, 125 subjects), covering 4 cancers (multiple myeloma, CLL, acute myeloid leukaemia and acute lymphoblastic leukaemia) and a large dose range (2.5-75 mg). A population pharmacokinetic model was developed with NONMEM and patient demographics were tested as covariates. RESULTS: The data were best fitted by a 1-compartment kinetic model with absorption described by 7 transit compartments. Clearance and volume of distribution were allometrically scaled for fat-free mass. The population parameter estimates for apparent clearance, apparent volume of distribution and transit rate constant were 12 L/h (10.8-13.6), 68.8 L (61.8-76.3), and 13.5 h-1 (11.9-36.8) respectively. Patients with impaired renal function (creatinine clearance <30 mL/min) exhibited a 22% reduction in lenalidomide clearance compared to patients with creatinine clearance of 90 mL/min. Cancer type had no discernible effect on lenalidomide disposition. CONCLUSIONS: This is the first report of a lenalidomide population pharmacokinetic model to evaluate lenalidomide pharmacokinetics in patients with CLL and compare its pharmacokinetics with other B-cell malignancies. As no differences in pharmacokinetics were found between the observed cancer-types, the unique toxicities observed in CLL may be due to disease-specific pharmacodynamics.

Converting from Transdermal to Buccal Formulations of Buprenorphine: A Pharmacokinetic Meta-Model Simulation in Healthy Volunteers.

Objective: To develop a model to predict buprenorphine plasma concentrations during transition from transdermal to buccal administration. Design: Population pharmacokinetic model-based meta-analysis of published data. Methods: A model-based meta-analysis of available buprenorphine pharmacokinetic data in healthy adults, extracted as aggregate (mean) data from published literature, was performed to explore potential conversion from transdermal to buccal buprenorphine. The time course of mean buprenorphine plasma concentrations following application of transdermal patch or buccal film was digitized from available literature, and a meta-model was developed using specific pharmacokinetic parameters (e.g., absorption rate, apparent clearance, and volumes of distribution) derived from analysis of pharmacokinetic data for intravenously, transdermally, and buccally administered buprenorphine. Results: Data from six studies were included in this analysis. The final transdermal absorption model employed a zero-order input rate that was scaled to reflect a nominal patch delivery rate and time after patch application (with decline in rate over time). The transdermal absorption rate constant became zero following patch removal. Buccal absorption was a first-order process with a time lag and bioavailability term. Simulations of conversion from transdermal 20 mcg/h and 10 mcg/h to buccal administration suggest that transition can be made rapidly (beginning 12 hours after patch removal) using the recommended buccal formulation titration increments (75-150 mcg) and schedule (every four days) described in the product labeling. Conclusions: Computer modeling and simulations using a meta-model built from data extracted from publications suggest that rapid and straightforward conversion from transdermal to buccal buprenorphine is feasible.

Population in vitro-in vivo pharmacokinetic model of first-pass metabolism: itraconazole and hydroxy-itraconazole.

The aim of this study was to develop a population in vitro-in vivo pharmacokinetic model that simultaneously describe the absorption and accumulation kinetics of itraconazole (ICZ) and hydroxy-itraconazole (HICZ) in healthy subjects. The model integrated meta-models of gastrointestinal pH and gastrointestinal transit time and in vitro dissolution models of ICZ with the absorption and disposition kinetics of ICZ and HICZ. Mean concentration intravenous data, and single- and multi-dose oral data were used for model development. Model development was conducted in NONMEM in a stepwise manner. First, a model of intravenous data (systemic kinetics) was established and then extended to include the oral data. The latter was then extended to establish the in vitro-in vivo pharmacokinetic model. The systemic disposition of ICZ was best described by a 3-compartment model with oral absorption described by 4-transit compartments and HICZ distribution by a 1-compartment model. ICZ clearance was best described using a mixed inhibition model that allowed HICZ concentrations to inhibit the clearance of parent drug. HICZ clearance was described by Michaelis-Menten elimination kinetics. An in vitro-in vivo model was successfully established for both formulations. The presented model was able to describe ICZ and HICZ plasma concentrations over a wide range of oral and intravenous doses and allowed the exploration of complexities associated with the non-linear ICZ and HICZ kinetics. The model may provide insight into the variability in exposure of ICZ with respect to relating in vivo dissolution characteristics with in vivo disposition kinetics.

Population Pharmacokinetic Model of Doxycycline Plasma Concentrations Using Pooled Study Data.

The literature presently lacks a population pharmacokinetic analysis of doxycycline. This study aimed to develop a population pharmacokinetic model of doxycycline plasma concentrations that could be used to assess the power of bioequivalence between Doryx delayed-release tablets and Doryx MPC. Doxycycline pharmacokinetic data were available from eight phase 1 clinical trials following single/multiple doses of conventional-release doxycycline capsules, Doryx delayed-release tablets, and Doryx MPC under fed and fasted conditions. A population pharmacokinetic model was developed in a stepwise manner using NONMEM, version 7.3. The final covariate model was developed according to a forward inclusion (P < 0.01) and then backward deletion (P < 0.001) procedure. The final model was a two-compartment model with two-transit absorption compartments. Structural covariates in the base model included formulation effects on relative bioavailability (F), absorption lag (ALAG), and the transit absorption rate (KTR) under the fed status. An absorption delay (lag) for the fed status (FTLAG2 = 0.203 h) was also included in the model as a structural covariate. The fed status was observed to decrease F by 10.5%, and the effect of female sex was a 14.4% increase in clearance. The manuscript presents the first population pharmacokinetic model of doxycycline plasma concentrations following oral doxycycline administration. The model was used to assess the power of bioequivalence between Doryx delayed-release tablets and Doryx MPC, and it could potentially be used to critically examine and optimize doxycycline dose regimens.

Intracellular CD3+ T Lymphocyte Teriflunomide Concentration Is Poorly Correlated with and Has Greater Variability Than Unbound Plasma Teriflunomide Concentration.

Leflunomide's active metabolite teriflunomide inhibits dihydro-oroate dehydrogenase, an enzyme essential to proliferation of T lymphocytes. As teriflunomide must reach the target site to have this effect, this study assessed the distribution of teriflunomide into T lymphocytes, as intracellular concentrations may be a superior response biomarker to plasma concentrations. CD3 MicroBeads (Miltenyi Biotec, Bergisch Gladbach, Germany) were used to extract CD3+ T cells from the peripheral blood of patients with rheumatoid arthritis who were taking a stable dose of leflunomide. Unbound plasma and intra-CD3+ T cell teriflunomide concentrations were quantified using liquid chromatography-mass spectrometry. Concentration (log transformed) and partition differences were assessed through paired Student t tests. Sixteen patients provided plasma steady-state teriflunomide samples, and eight provided a sample 6-12 weeks later. At time-point one, the geometric mean teriflunomide concentration (range) in CD3+ T cells was 18.12 µg/L (6.15-42.26 µg/L) compared with 69.75 µg/L (32.89-263.1 µg/L) unbound in plasma (P < 0.001). The mean partition coefficient (range) for unbound plasma teriflunomide into CD3+ T cells was 0.295 (0.092-0.632), which was significantly different from unity (P < 0.001). The median (range) change in teriflunomide concentration between the two time points was 14% (-10% to 40%) in unbound plasma and -29% (-69 to 138%) for CD3+ T cells. Because teriflunomide concentrations in CD3+ T cells were lower and displayed a higher intraindividual variability than the unbound plasma concentrations, its applicability as a therapeutic drug-monitoring marker may be limited.

Comparison of non-compartmental and mixed effect modelling methods for establishing bioequivalence for the case of two compartment kinetics and censored concentrations.

Non-compartmental analysis (NCA) is regarded as the standard for establishing bioequivalence, despite its limitations and the existence of alternative methods such as non-linear mixed effects modelling (NLMEM). Comparisons of NCA and NLMEM in bioequivalence testing have been limited to drugs with one-compartment kinetics and have included a large number of different approaches. A simulation tool was developed with the ability to rapidly compare NCA and NLMEM methods in determining bioequivalence using both R and NONMEM and applied to a drug with two-compartment pharmacokinetics. Concentration-time profiles were simulated where relative bioavailability, random unexplained variability (RUV) at the lower limit of quantification (LLOQ) differed between simulations. NLMEM analyses employed either the M1 or M3 methods for dealing with values below the LLOQ. It was used to elucidate the impact of changes in (i) RUV at the LLOQ, (ii) the extent of censoring data below the LLOQ and (iii) the concentration sampling times. The simulations showed NLMEM having a consistent 20-40% higher accuracy and sensitivity in identifying bioequivalent studies when compared to NCA, while NCA was found to have a 1-10% higher specificity than NLMEM. Increasing data censoring by increasing the LLOQ resulted in decreases of ~10% to the accuracy and sensitivity of NCA, with minimal effects on NLMEM. The tool provides a platform for comparing NCA and NLMEM methods and its use can be extended beyond the scenarios reported here. In the situations examined it is seen that NLMEM is more accurate than NCA and may offer some advantages in the determination of bioequivalence.

A model-based evaluation of single metrics for discriminating changes in rheumatoid arthritis disease activity.

AIMS: Composite indices for quantifying rheumatoid arthritis (RA) disease activity such as the 28-joint disease activity score (DAS28) are comprised of single parameters ('metrics') in various combinations. Population modelling methods were used to evaluate single metrics for their ability to reflect changes in disease activity with a view to understanding and improving composite indices. METHODS: A total of 11 single metrics of RA disease activity (tender and swollen joint counts, acute phase reactants and global health, pain and physical function assessments) were obtained from 203 patients with recent onset RA. Participants received combination disease-modifying anti-rheumatic drugs (DMARDs) according to a treat-to-target approach with a pre-defined protocol for treatment intensification. Models describing each metric's magnitude and variability of change from baseline to a single 'treated' state in the population were developed using nonmem(®) . Measures that displayed uniformly large changes between states across the population were ranked higher in terms of discriminatory capacity. RESULTS: Joint counts demonstrated a greater ability to discriminate changes in RA disease activity than others. Correlations between metrics demonstrated that erythrocyte sedimentation rate (ESR) had limited relationships with others for baseline scores and changes in RA disease activity (r generally < 0.2). However it appeared to be important in describing changes for those individuals where ESR levels were initially elevated. CONCLUSION: It appears unlikely that a single group of metrics may be suitable to capture disease activity changes across all RA patients and defining the most appropriate metric(s) for individual patients will be an important area of future research.

Genetic polymorphism of CYP1A2 but not total or free teriflunomide concentrations is associated with leflunomide cessation in rheumatoid arthritis.

AIM: Leflunomide, via its active metabolite teriflunomide, is used in rheumatoid arthritis (RA) treatment, yet approximately 20 to 40% of patients cease due to toxicity. The aim was to develop a time-to-event model describing leflunomide cessation due to toxicity within a clinical cohort and to investigate potential predictors of cessation such as total and free teriflunomide exposure and pharmacogenetic influences. METHODS: This study included individuals enrolled in the Early Arthritis inception cohort at the Royal Adelaide Hospital between 2000 and 2013 who received leflunomide. A time-to-event model in nonmem was used to describe the time until leflunomide cessation and the influence of teriflunomide exposure and pharmacogenetic variants. Random censoring of individuals was simultaneously described. The clinical relevance of significant covariates was visualized via simulation. RESULTS: Data from 105 patients were analyzed, with 34 ceasing due to toxicity. The baseline dropout hazard and baseline random censoring hazard were best described by step functions changing over discrete time intervals. No statistically significant associations with teriflunomide exposure metrics were identified. Of the screened covariates, carriers of the C allele of CYP1A2 rs762551 had a 2.29 fold increase in cessation hazard compared with non-carriers (95% CI 2.24, 2.34, P = 0.016). CONCLUSIONS: A time-to-event model described the time between leflunomide initiation and cessation due to side effects. The C allele of CYP1A2 rs762551 was linked to increased leflunomide toxicity, while no association with teriflunomide exposure was identified. Future research should continue to investigate exposure-toxicity relationships, as well as potentially toxic metabolites.

Population In Vitro-In Vivo Correlation Model Linking Gastrointestinal Transit Time, pH, and Pharmacokinetics: Itraconazole as a Model Drug.

PURPOSE: To establish an in vitro-in vivo correlation (IVIVC) model for Sporanox and SUBA-itraconazole formulations and to understand the impact of gastrointestinal (GI) pH and transit times on itraconazole dissolution and absorption. METHODS: IVIVC was developed based on fed/fasted pharmacokinetic data from randomized cross-over trials, in vitro dissolution studies, and prior information about typical and between subject variability of GI pH and transit times. Data were analysed using the population modelling approach as implemented in NONMEM. RESULTS: Dissolution kinetics were described using first order models. The in vivo pharmacokinetics of itraconazole was described with a 2-compartment model with 4-transit absorption compartments. Pharmacokinetic profiles for fasted itraconazole periods were described based on the in vitro dissolution model, in vivo disposition model, and the prior information on GI pH and transit times. The IVIVC model indicated that drug dissolution in the fed state required an additional pH-independent dissolution pathway. The IVIVC models were presented in a 'Shiny' application. CONCLUSION: An IVIVC model was established and internally evaluated for the two itraconazole formulations. The IVIVC model provides more insight into the observed variability of itraconazole pharmacokinetics and indicated that GI pH and transit times influence in vivo dissolution and exposure.

An introduction to physiologically-based pharmacokinetic models.

Physiologically-based pharmacokinetic (PBPK) models represent drug kinetics in one or more 'real' organs (and hence require submodels of organs/tissues) and they describe 'whole-body' kinetics by joining together submodels with drug transport by blood flow as dictated by anatomy. They attempt to reproduce 'measureable' physiological and/or pharmacokinetic processes rather than more abstract rate constants and volumes. PBPK models may be built using a 'bottom-up' approach, where parameters are chosen from first principles, literature, or in vitro data as opposed to a 'top-down' approach, where all parameters are estimated from data. The basic principles of PBPK models are described, focusing on the equations for three individual organs: a single flow-limited compartment describing distribution only, a membrane-limited compartment describing distribution, and a single flow-limited compartment with elimination. These organ models are linked to make a basic three-compartment physiological model of the whole body. PBPK models are particularly suited to scaling kinetics across body size (e.g., adult to neonate) and species (e.g., animal to first-in-man) as physiology and pharmacology can be represented by independent parameters. Maturation models can be incorporated as for compartmental models. PBPK models are now available in commercial software packages, and are perhaps now more accessible than ever. Alternatively, even complex PBPK models can be represented in generic differential equation-solving software using the simple principles described here. The relative ease of constructing the code for PBPK models belies the most difficult aspect of their implementation-collecting, collating, and justifying the data used to parameterize the model.

Population pharmacokinetic modeling of itraconazole and hydroxyitraconazole for oral SUBA-itraconazole and sporanox capsule formulations in healthy subjects in fed and fasted states.

Itraconazole is an orally active antifungal agent that has complex and highly variable absorption kinetics that is highly affected by food. This study aimed to develop a population pharmacokinetic model for itraconazole and the active metabolite hydroxyitraconazole, in particular, quantifying the effects of food and formulation on oral absorption. Plasma pharmacokinetic data were collected from seven phase I crossover trials comparing the SUBA-itraconazole and Sporanox formulations of itraconazole. First, a model of single-dose itraconazole data was developed, which was then extended to the multidose data. Covariate effects on itraconazole were then examined before extending the model to describe hydroxyitraconazole. The final itraconazole model was a 2-compartment model with oral absorption described by 4-transit compartments. Multidose kinetics was described by total effective daily dose- and time-dependent changes in clearance and bioavailability. Hydroxyitraconazole was best described by a 1-compartment model with mixed first-order and Michaelis-Menten elimination for the single-dose data and a time-dependent clearance for the multidose data. The relative bioavailability of SUBA-itraconazole compared to that of Sporanox was 173% and was 21% less variable between subjects. Food resulted in a 27% reduction in bioavailability and 58% reduction in the transit absorption rate constant compared to that with the fasted state, irrespective of the formulation. This analysis presents the most extensive population pharmacokinetic model of itraconazole and hydroxyitraconazole in the literature performed in healthy subjects. The presented model can be used for simulating food effects on itraconazole exposure and for performing prestudy power analysis and sample size estimation, which are important aspects of clinical trial design of bioequivalence studies.

Population pharmacokinetics of orally administered mefloquine in healthy volunteers and patients with uncomplicated Plasmodium falciparum malaria.

BACKGROUND: The determination of dosing regimens for the treatment of malaria is largely empirical and thus a better understanding of the pharmacokinetic/pharmacodynamic properties of antimalarial agents is required to assess the adequacy of current treatment regimens and identify sources of suboptimal dosing that could select for drug-resistant parasites. Mefloquine is a widely used antimalarial, commonly given in combination with artesunate. PATIENTS AND METHODS: Mefloquine pharmacokinetics was assessed in 24 healthy adults and 43 patients with Plasmodium falciparum malaria administered mefloquine in combination with artesunate. Population pharmacokinetic modelling was conducted using NONMEM. RESULTS: A two-compartment model with a single transit compartment and first-order elimination from the central compartment most adequately described mefloquine concentration-time data. The model incorporated population parameter variability for clearance (CL/F), central volume of distribution (VC/F) and absorption rate constant (KA) and identified, in addition to body weight, malaria infection as a covariate for VC/F (but not CL/F). Monte Carlo simulations predict that falciparum malaria infection is associated with a shorter elimination half-life (407 versus 566 h) and T>MIC (766 versus 893 h). CONCLUSIONS: This is the first known population pharmacokinetic study to show falciparum malaria to influence mefloquine disposition. Protein binding, anaemia and other factors may contribute to differences between healthy individuals and patients. As VC/F is related to the earlier portion of the concentration-time profiles, which occurs during acute malaria, and CL/F is more related to the terminal phase during convalescence after treatment, this may explain why malaria was found to be a covariate for VC/F but not CL/F.

A population model of early rheumatoid arthritis disease activity during treatment with methotrexate, sulfasalazine and hydroxychloroquine.

AIMS: To develop a population model describing the disease activity (DAS28) time course in patients with early rheumatoid arthritis (RA) treated with triple disease-modifying anti-rheumatic drug (DMARD) therapy (methotrexate, sulfasalazine and hydroxychloroquine). METHODS: DAS28 was obtained in 263 patients with early RA from initiation of therapy until 60 weeks. Using NONMEM(®), base models (DAS28 vs. time) and covariate influences were investigated for the population. RESULTS: The best model was an exponential model of DAS28 vs. time that was additive to baseline DAS28, with covariance between parameters, and a combined residual error model. Age and patient smoking history were covariates significantly affecting response to therapy. Population estimates were baseline DAS28 (5.7), extent of change in DAS28 (-2.8) and the half-life of disease activity (6.2 weeks; time to steady disease state achieved within approximately 30 weeks). Older individuals exhibited more severe baseline DAS28, described by a power function centred around 57 years (baseline DAS28 for 40- and 70-year-old patients were 5.4 vs. 5.8, respectively) and current smokers took longer to achieve a steady disease state (approximately 50 weeks). There was considerable within-patient random variability in DAS28 over time (empirical 90% CI for DAS28 in a population typical patient at 60 weeks: 1.8, 4.2 with median value of 2.8). CONCLUSIONS: This is the first report of a disease activity model for early RA treated with triple DMARD therapy. Smoking and age were identified as covariates.

Pharmacokinetics of tramadol after subcutaneous administration in a critically ill population and in a healthy cohort.

BACKGROUND: Tramadol is an atypical centrally acting analgesic agent available as both oral and parenteral preparations. For patients who are unable to take tramadol orally, the subcutaneous route of administration offers an easy alternative to intravenous or intramuscular routes. This study aimed to characterise the absorption pharmacokinetics of a single subcutaneous dose of tramadol in severely ill patients and in healthy subjects. METHODS/DESIGN: Blood samples (5 ml) taken at intervals from 2 minutes to 24 hours after a subcutaneous dose of tramadol (50 mg) in 15 patients (13 male, two female) and eight healthy male subjects were assayed using high performance liquid chromatography. Pharmacokinetic parameters were derived using a non-compartmental approach. RESULTS: There were no statistically significant differences between the two groups in the following parameters (mean ± SD): maximum venous concentration 0.44 ± 0.18 (patients) vs. 0.47 ± 0.13 (healthy volunteers) mcg/ml (p = 0.67); area under the plasma concentration-time curve 177 ± 109 (patients) vs. 175 ± 75 (healthy volunteers) mcg/ml*min (p = 0.96); time to maximum venous concentration 23.3 ± 2 (patients) vs. 20.6 ± 18.8 (healthy volunteers) minutes (p = 0.73) and mean residence time 463 ± 233 (patients) vs. 466 ± 224 (healthy volunteers) minutes (p = 0.97). CONCLUSIONS: The similar time to maximum venous concentration and mean residence time suggest similar absorption rates between the two groups. These results indicate that the same dosing regimens for subcutaneous tramadol administration may therefore be used in both healthy subjects and severely ill patients. TRIAL REGISTRATION: ACTRN12611001018909.

"Getting the Dose Right"-Revisiting the Topic With Focus on Biologic Agents.

Nearly two decades after the Peck and Cross article '"Getting the dose right: facts, a blueprint, and encouragements" was published, a review of dose recommendations for biologics shows that the success in getting the dose right appears to have improved given the relatively low incidence of drug withdrawals and dosing/label changes. However, the clinical experience with monoclonal antibodies (MAbs) following approval has been less than perfect. In inflammatory diseases, the disease burden changes with time and high treatment failure rates have been reported. In addition, the use of concomitant steroids and immunosuppressant drugs with MAbs is common. These concomitant agents have their own safety issues and many immunosuppressant agents are not well-tolerated although they have been shown to reduce the incidence of anti-drug antibodies (ADA). This same complexity is seen in MAbs used in oncology as well, although with these agents the doses appear to be higher than needed, which results in high treatment costs and incidence of adverse events. Given the complexity of MAb pharmacokinetics, which makes providing a detailed description of dose options difficult, product labeling should include the options for alternative dose strategies and potentially include the use of therapeutic drug monitoring with dose individualization which have been shown to improve clinical response and reduce the incidence of ADA. So, while the recommended dosing for biologics seems improved over the issues noted 17 years ago, we still have some work to do.

Application of a Precision-Dosing Model to a Real-World Cohort of Patients on Infliximab Maintenance Therapy: Drug Usage and Cost Analysis.

Precision-dosing models forecast infliximab doses to achieve targeted trough concentrations in patients with inflammatory bowel disease (IBD). These models have shown to reduce nonresponse and improve patient outcomes. We compared infliximab doses determined by iDOSE precision dosing with standard dosing, and the associated drug costs, in patients with IBD. In this retrospective study, patients with IBD treated with infliximab every 8 weeks at 5 mg/kg were included. An infliximab dose was named dose X if 3 previous infliximab doses, laboratory values including trough infliximab concentrations, and the patient's weight were recorded. The actual dose X was compared to an iDOSE-predicted dose X. Net drug use and costs were evaluated. A total of 174 patients-56% men; median age, 36 (interquartile range, 29-47) years; 135 with Crohn disease; and 31 with ulcerative colitis-were included, with 417 dose X recordings. Median prior infliximab therapy was 2 (0-4) years. Comparing actual dose X with predicted dose X, 52% and 32% of doses were subtherapeutic when aiming for trough concentrations of 5-10 and 3-7 µg/mL, respectively. Treatment costs increased by 102% and 29% for the 2 trough ranges, respectively. On multivariate regression analysis, subtherapeutic infliximab concentrations were associated with ulcerative colitis compared with Crohn disease (odds ratio, 9.81; 95% confidence interval, 1.28-75.40; P = .028) and predose X infliximab trough concentration [odds ratio, 0.07; 95% confidence interval, 0.03-0.15; P < .001]. Over half of maintenance infliximab drug doses were too low to achieve infliximab blood concentrations of 5 µg/mL or greater. While applying precision dosing may improve patient outcomes, drug costs could be considerably greater.

Treatment Targets Should Influence Choice of Infliximab Dose Intensification Strategy in Inflammatory Bowel Disease: A Pharmacokinetic Simulation Study.

BACKGROUND: The optimal infliximab dose intensification strategy to address loss of response associated with subtherapeutic infliximab trough levels remains uncertain, as does whether post-intensification trough and treatment targets should influence this decision. OBJECTIVES: This pharmacokinetic simulation study aimed to identify infliximab dose intensification strategies capable of achieving post-intensification infliximab trough thresholds associated with clinical and objective treatment targets in Crohn's disease and ulcerative colitis. METHODS: A validated pharmacokinetic infliximab model, applied to 200 simulated patients, identified those with subtherapeutic (< 3.00 mg/L) trough levels after 30 weeks of standard (5 mg/kg 8-weekly) dosing, and subsequently applied 10 dose intensification strategies over a further 32 weeks. The proportion of simulations achieving 32-week post-intensification infliximab trough levels associated with endoscopic remission (ulcerative colitis > 7.50 mg/L, Crohn's disease > 9.70 mg/L) was the primary outcome, with perianal fistula healing (Crohn's disease > 10.10 mg/L) and clinical improvement (ulcerative colitis > 3.70 mg/L, Crohn's disease > 7.00mg/L) evaluated as secondary outcomes. All outcomes were stratified by intensity of dose intensification, with standard (≤ 10 mg/kg 8-weekly or 5 mg/kg 4-weekly; n = 5) and intensive (> 10 mg/kg 8-weekly or 5 mg/kg 4-weekly; n = 5) dosing strategies defined, respectively. RESULTS: The median pre-intensification infliximab trough level was 0.91 mg/L (interquartile range 1.37). Intensive dosing strategies were more likely to achieve infliximab trough concentrations associated with endoscopic remission (ulcerative colitis 36.48% vs. 10.80%, Crohn's disease 25.98 vs. 4.68%), perianal fistula healing (24.52% vs. 4.36%) and clinical improvement (ulcerative colitis 61.90% vs. 34.86%, Crohn's disease 40.32 vs. 12.08%) than standard intensification strategies (all p < 0.01). When controlling for cumulative (mg/kg) infliximab dose over 32 weeks, strategies that concurrently dose increased and interval shortened achieved the highest infliximab trough levels (all p < 0.01). CONCLUSION: This simulation-based analysis highlights the potential of using post-intensification infliximab trough thresholds associated with aspirational treatment targets in Crohn's disease and ulcerative colitis to guide choice of infliximab dose intensification strategy. Intensive dose intensification strategies, particularly those that concurrently dose increase and interval shorten, appear to achieve higher infliximab levels than standard dose intensification strategies. This may be particularly important in the pursuit of stringent endpoints, such as endoscopic remission and fistula healing, which have been consistently associated with higher infliximab trough levels. These findings require validation across real-world cohorts.

Dynamic exposure and body burden models for per- and polyfluoroalkyl substances (PFAS) enable management of food safety risks in cattle.

With increasing global focus on planetary boundaries, food safety and quality, the presence of per- and polyfluoroalkyl substances (PFAS) in the food chain presents a challenge for the sustainable production and supply of quality assured food. Consumption of food is the primary PFAS exposure route for the general population. At contaminated sites, PFAS have been reported in a range of agricultural commodities including cattle. Consumer exposure assessments are complicated by the lack of validated modelling approaches to estimate PFAS bioaccumulation in cattle. Previous studies have shown that PFAS bioaccumulation in livestock is influenced by environmental, spatial and temporal factors that necessitate a dynamic modelling approach. This work presents an integrated exposure and population toxicokinetic (PopTK) model for cattle that estimates serum and tissue concentrations of PFAS over time. Daily exposures were estimated from intakes of water, pasture, and soil, and considered animal growth, seasonal variability (pasture moisture content and temperature) and variable PFAS concentrations across paddocks. Modelled serum and tissue estimates were validated against monitoring data from Australian and Swedish cattle farms. The models were also used to develop and test practical management options for reducing PFAS exposure and to prioritise remediation for farms. Model outputs for exposure management scenarios (testing cattle rotation and targeted supplementation of feed and water) showed potential for marked reductions in consumer exposures from cattle produce.