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.

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.

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.

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.

A physiologically-based recirculatory meta-model for nasal fentanyl in man.

Pharmacokinetic (PK) and pharmacodynamic (PD) data were available from a study of a nasal delivery system for the opioid analgesic fentanyl, together with data on the kinetics of fentanyl in arterial blood in man, and in the lung and brain of sheep. Our aim was to reconcile these data using a physiologically-based population recirculatory PK-PD model, with emphasis on achieving a meta-model that could simultaneously account for the arterial and venous (arm) concentrations of fentanyl, could relate PD effects (pain scores) to the CNS concentrations of fentanyl, and could account for the effect of body size and age on fentanyl kinetics. Data on the concentration gradients of fentanyl across brain, lung and muscle were used to develop sub-models of fentanyl kinetics in these organs. The sub-models were incorporated into a "whole body" recirculatory model by adding additional sub-models for a venous mixing compartment, the liver and gut, the kidney and the "rest of the body" with blood flows and organ volumes based on values for a Standard Man. Inter-individual variability was achieved by allometric scaling of organ size and blood flows, evidence-based assumptions about the effect of weight and age on cardiac output, and inter-individual variability in the free fraction in plasma and hepatic extraction of fentanyl. Post-operative pain scores were found to be temporally related to the predicted brain concentrations of fentanyl. We conclude that a physiologically-based meta-modelling approach was able to describe clinical PK-PD studies of fentanyl whilst providing a mechanistic interpretation of key aspects of its disposition.

Population Pharmacokinetics and Pharmacodynamics of the Therapeutic and Adverse Effects of Ketamine in Patients With Treatment-Refractory Depression.

We aimed to develop population pharmacokinetic/pharmacodynamic (PK/PD) models that can effectively describe ketamine and norketamine PK/PD relationships for Montgomery-Åsberg Depression Rating Scale (MADRS) scores, blood pressure (BP), and heart rate (HR) following i.v., s.c., and i.m. ketamine administration in patients with treatment-refractory depression. Ketamine PK/PD data were collected from 21 treatment-refractory depressed participants who received ketamine (dose titration 0.1-0.5 mg/kg as single doses) by i.v., s.c., or i.m. administration. Model development used nonlinear mixed effect modeling. Ketamine and norketamine PK were best described using two-compartment models with first-order absorption after s.c. and i.m. administration. Estimated ketamine bioavailability after i.m. and s.c. was ~ 64% with indistinguishable first-order absorption rate constants. Allometric scaling of body weight on all clearance and volumes of distribution improved the model fit. The delay in the concentration-response relationship for MADRS scores was best described using a turnover model (turnover time ~ 42 hours), whereas for the BP and HR rates this was an immediate effect model. For all PD effects, ketamine alone was superior to models with norketamine concentration linked to an effect. No covariates were identified for PD effects. The estimated half-maximal effective concentration from the MADRS score, BP, and HR were 0.44, 468, and 7,580 ng/mL, respectively. The integrated population models were able to effectively describe the PK/PD relationships for MADRS scores, BP, and HR after i.v., s.c., and i.m. ketamine administration. These findings allow for a deeper understanding of the complex relationships between route of ketamine administration and clinical response and safety.

The association of depression and anxiety with treatment outcomes in patients with rheumatoid arthritis - a pooled analysis of five randomised controlled trials.

Background: Rheumatoid arthritis (RA) is an inflammatory autoimmune condition associated with an increased risk of developing depression and anxiety. Depression and anxiety are associated with worse outcomes in RA, but the magnitude of the effect of each condition on RA outcomes is unclear. It is also unknown how pharmacological treatment of depression affects RA outcomes. Objective: The primary aim of this study was to investigate the association of comorbid depression and anxiety with remission in patients with RA. Secondary aims were to determine the association between comorbid depression and anxiety on patient-reported outcomes and the relationship between concomitant use of antidepressants and remission in patients with depression. Design: Data from patients with moderate to severe RA were pooled from five randomised controlled trials investigating tocilizumab and conventional synthetic disease-modifying agents. Methods: Remission was defined as a clinical disease activity index (CDAI) of ⩽2.8 and simple disease activity index (SDAI) of ⩽3.3. The association between the time to reach remission and depression and anxiety was analysed using Cox proportional hazard analysis. Results: Individual patient data were available from 5502 subjects, of whom 511 had depression, 236 had anxiety and 387 were using antidepressants. Depression was significantly associated with reduced remission [adjusted HR (95% CI): 0.62 (0.48-0.80), p < 0.001 and adjusted HR (95% CI): 0.59 (0.44-0.79), p < 0.001] using CDAI and SDAI, respectively. Depression was associated with a lower likelihood of achieving more subjective outcomes (⩽1 physician global assessment, ⩽1 patient global assessment) and ⩽1 28-swollen joint count, but not ⩽1 28-tender joint count or C-reactive protein measurement. Treatment with antidepressants did not improve outcomes for patients with depression. Anxiety was not significantly associated with RA remission. Conclusion: Comorbid depression, but not anxiety, was associated with less frequent remission. Concomitant antidepressant use was not associated with improvements in RA outcomes in patients with depression.

CYP2C19-guided design of a proton pump inhibitor dose regimen to avoid the need for pharmacogenetic individualization in H. pylori eradication.

OBJECTIVE: It has been demonstrated that genetic variation in CYP2C19 has a significant influence upon H. pylori eradication rates. We have determined a dosage regimen of lansoprazole that will provide EMs with exposure approximately equivalent to that obtained by PMs treated with standard doses and determined the exposure that a PM would experience if they were to be treated with this 'EM optimised' lansoprazole dose. METHODS: Non-compartmental pharmacokinetic parameters (AUC, C(max), t(max)) for CYP2C19 genotypes were obtained from the literature. Primary pharmacokinetic parameters (CL, Vd, ka) for 200 virtual patients were calculated from the weighted non-compartmental variables and used to simulate a 7 day treatment course of twice daily lansoprazole, at standard and optimised doses for 1,000 patients. RESULTS: The administration of 180 mg twice daily to CYP2C19 EMs results in approximately equivalent exposure to lansoprazole as the administration of standard 30 mg twice daily doses to PMs. Administration of this six-fold dose increase to EMs is predicted to result in only a 2.5-fold increase in C(max) when compared with PMs receiving the standard 30 mg dose. CONCLUSION: We present a potential lansoprazole dosing regimen that should result in improved H. pylori eradication within CYP2C19 EMs and may not require individualization. Whilst the optimised dose represents a significant increase, it is below that reported in the chronic management of Zollinger-Ellison syndrome. On the basis that treatment is of limited duration and lansoprazole is generally well tolerated, such an approach warrants further in vivo evaluation to confirm drug exposure, efficacy and tolerability.

Concomitant beta-blocker use is associated with a reduced rate of remission in patients with rheumatoid arthritis treated with disease-modifying anti-rheumatic drugs: a post hoc multicohort analysis.

BACKGROUND: Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease associated with increased risk of cardiovascular disease (CVD). Treatment for CVD may involve pharmacological agents that antagonise beta adrenergic receptors. These receptors may play an important role in immunology, and the effects of beta-blockers (BB) in RA is unknown. The aim of this study was to investigate the association between BB use and remission in patients with RA initiating tocilizumab +/- conventional synthetic (cs-) DMARD therapy. METHODS: Data was pooled from five randomised trials investigating tocilizumab and/or csDMARD treatment in RA (primarily methotrexate). The association between BB use and remission according to the Simplified Disease Activity Index (SDAI) and Clinical Disease Activity Index (CDAI) was assessed by Cox proportional hazard analysis. Sensitivity analysis in patients with pre-existing CVD and an exploratory analysis of the impact of other CVD drugs were conducted. RESULTS: Data were available from 5502 participants, 594 (10.8%) of whom were using systemic BB. BB use was associated with less frequent SDAI remission in the total [adjusted hazard ratio (HR) 0.70, 95% confidence interval (CI) 0.57-0.87, p = 0.001] and CVD cohort [adjusted HR 0.72 (0.57-0.90, p = 0.005)]. The association was consistent between trials (interaction p = 0.44) and treatment arms (interaction p = 0.06). No significant association between remission and ß1-receptor selectivity was identified (p = 0.16), and the association was independent from other cardiovascular drug use. Similar associations between BB use and CDAI remission were observed. CONCLUSION: In a large, pooled cohort of RA patients initiating csDMARDs and/or tocilizumab, BB use was independently associated with less frequent remission.

Effect of Roux-en-Y gastric bypass on the pharmacokinetic-pharmacodynamic relationships of liquid and controlled-release formulations of oxycodone.

The physiological changes following Roux-en-Y gastric bypass (RYGB) surgery may impact drug release from mechanistically different controlled-release tablets, making generic substitution inappropriate. This study aimed to characterise the pharmacokinetic-pharmacodynamic relationships of oxycodone from a lipid-based and water-swellable controlled-release tablet in RYGB patients. Twenty RYGB patients received 10-mg oral solution oxycodone or 20-mg controlled-release (water-swellable or lipid-based) oxycodone in a three-way, randomised, semiblinded and cross-over study. Blood sampling and pupillary recordings were conducted over a 24-h period. A previously established pharmacokinetic-pharmacodynamic model of these three formulations in healthy volunteers was used in the analysis as a reference model. No differences in absorption kinetics were seen between controlled-release formulations in patients. However, the absorption lag time was 11.5 min in patients vs 14 min in healthy volunteers for controlled-release tablets (P < 0.001). Furthermore, oral bioavailability was 14.4% higher in patients compared to healthy volunteers regardless of formulation type (P < 0.001). Oxycodone pharmacodynamics were not significantly affected by formulation or patient status. However, baseline pupil diameter was inversely correlated with age (P < 0.001) and plasma concentrations of oxycodone at half-maximum effect were 31% lower in males compared to females (P < 0.05). Generic substitution of monophasic lipid-based and water-swellable controlled-release oxycodone tablets may be considered safe in RYGB patients.

(R)- and (S)-methadone and buprenorphine concentration ratios in maternal and umbilical cord plasma following chronic maintenance dosing in pregnancy.

AIMS: The aim of this study was to compare the transfer of buprenorphine and methadone between maternal and cord blood in women under chronic dosing conditions and to determine if differences exist in the transfer of the two methadone enantiomers. METHODS: Maternal and cord blood samples were collected at delivery from women maintained on methadone (35, 25-140 mg day⁻¹) (median; range) or buprenorphine (6.00, 2-20 mg day⁻¹) during pregnancy. Plasma concentration ratios are presented as an indicator of foetal exposure relative to the mother. RESULTS: Methadone was quantified in all samples, with cord : maternal plasma methadone concentration ratios (n= 15 mother-infant pairs) being significantly higher (P < 0.0001; mean difference (MD) 0.07; 95% confidence interval (CI) 0.048, 0.092) for the active (R)-methadone enantiomer (0.41; 0.19, 0.56) (median; range) compared with (S)-methadone (0.36; 0.15, 0.53). (R)- : (S)-methadone concentration ratios were also significantly higher (P < 0.0001; MD 0.24 95% CI 0.300, 0.180) for cord (1.40; 0.95, 1.67) compared with maternal plasma (1.16; 0.81, 1.38). Half the infant buprenorphine samples were below the assay lower limit of quantification (LLOQ) (0.125 ng ml⁻¹). The latter was four-fold lower than the LLOQ for methadone (0.50 ng ml⁻¹). The cord : maternal plasma buprenorphine concentration ratio (n= 9 mother-infant pairs) was 0.35; 0.14, 0.47 and for norbuprenorphine 0.49; 0.24, 0.91. CONCLUSIONS: The transfer of the individual methadone enantiomers to the foetal circulation is stereoselective. Infants born to buprenorphine maintained women are not exposed to a greater proportion of the maternal dose compared with methadone and may be exposed to relatively less of the maternal dose compared with infants born to women maintained on methadone during pregnancy.

Association between obesity and remission in rheumatoid arthritis patients treated with disease-modifying anti-rheumatic drugs.

The aim of this study was to investigate the association between body-mass index (BMI) and remission in RA patients receiving conventional synthetic (cs-) or the biological Disease-Modifying Antirheumatic Drug (DMARD), tocilizumab. Individual participant data (IPD) were pooled from five trials investigating tocilizumab and/or csDMARDs therapy (primarily methotrexate) for RA. Time to first remission was recorded according to the Simplified Disease Activity Index (SDAI) and Clinical Disease Activity Index (CDAI). BMI was classified according to WHO definitions. Associations between baseline BMI and remission were assessed by Cox-proportional hazard analysis. IPD were available from 5428 patients treated with tocilizumab ± csDMARDs (n = 4098) or csDMARDs alone (n = 1330). Of these, 1839 (33.9%) had normal BMI, 1780 (32.8%) overweight, 1652 (30.4%) obese and 157 (2.9%) were underweight. Obesity, compared to normal BMI, was associated with less frequent remission using SDAI (adjusted HR 0.80 [95% CI 0.70-0.92]) and CDAI (adjusted HR 0.77 [0.68-0.87]). As continuous variable, increased BMI was associated with less frequent SDAI (P = 0.001) and CDAI (P = 0.001) defined remission. No heterogeneity in identified associations was observed between studies (P = 0.08) or treatments (P = 0.22). Obesity was negatively associated with RA disease remission regardless of RA therapy, suggesting that baseline BMI should be considered as a stratification factor in future RA trials.