Design considerations and analysis planning of a phase 2a proof of concept study in rheumatoid arthritis in the presence of possible non-monotonicity

BackgroundIt is important to quantify the dose response for a drug in phase 2a clinical trials so the optimal doses can then be selected for subsequent late phase trials. In a phase 2a clinical trial of new lead drug being developed for the treatment of rheumatoid arthritis (RA), a U-shaped dose response curve was observed. In the light of this result further research was undertaken to design an efficient phase 2a proof of concept (PoC) trial for a follow-on compound using the lessons learnt from the lead compound.MethodsThe planned analysis for the Phase 2a trial for GSK123456 was a Bayesian Emax model which assumes the dose-response relationship follows a monotonic sigmoid “S” shaped curve. This model was found to be suboptimal to model the U-shaped dose response observed in the data from this trial and alternatives approaches were needed to be considered for the next compound for which a Normal dynamic linear model (NDLM) is proposed. This paper compares the statistical properties of the Bayesian Emax model and NDLM model and both models are evaluated using simulation in the context of adaptive Phase 2a PoC design under a variety of assumed dose response curves: linear, Emax model, U-shaped model, and flat response.ResultsIt is shown that the NDLM method is flexible and can handle a wide variety of dose-responses, including monotonic and non-monotonic relationships. In comparison to the NDLM model the Emax model excelled with higher probability of selecting ED90 and smaller average sample size, when the true dose response followed Emax like curve. In addition, the type I error, probability of incorrectly concluding a drug may work when it does not, is inflated with the Bayesian NDLM model in all scenarios which would represent a development risk to pharmaceutical company.The bias, which is the difference between the estimated effect from the Emax and NDLM models and the simulated value, is comparable if the true dose response follows a placebo like curve, an Emax like curve, or log linear shape curve under fixed dose allocation, no adaptive allocation, half adaptive and adaptive scenarios. The bias though is significantly increased for the Emax model if the true dose response follows a U-shaped curve.ConclusionsIn most cases the Bayesian Emax model works effectively and efficiently, with low bias and good probability of success in case of monotonic dose response. However, if there is a belief that the dose response could be non-monotonic then the NDLM is the superior model to assess the dose response.

[1]  Z. Q. John Lu,et al.  Bayesian methods for data analysis, third edition , 2010 .

[2]  Ute Beyer,et al.  Bayesian Forecasting And Dynamic Models , 2016 .

[3]  Brenda Gaydos,et al.  Application of Adaptive Design Methodology in Development of a Long-Acting Glucagon-like Peptide-1 Analog (Dulaglutide): Statistical Design and Simulations , 2012, Journal of diabetes science and technology.

[4]  G. E. Dille Guidance , 1930 .

[5]  Stefano Zamuner,et al.  Safety, tolerability, pharmacokinetics and pharmacodynamics of an anti- oncostatin M monoclonal antibody in rheumatoid arthritis: results from phase II randomized, placebo-controlled trials , 2013, Arthritis Research & Therapy.

[6]  Frank Bretz,et al.  Analysis of Dose–Response Studies—Modeling Approaches , 2006 .

[7]  Adrian Raftery,et al.  The Number of Iterations, Convergence Diagnostics and Generic Metropolis Algorithms , 1995 .

[8]  Brian K. Shoichet,et al.  Colloidal Drug Formulations Can Explain “Bell-Shaped” Concentration–Response Curves , 2014, ACS chemical biology.

[9]  Andrew P Grieve,et al.  ASTIN: a Bayesian adaptive dose–response trial in acute stroke , 2005, Clinical trials.

[10]  G. Stucki,et al.  Rheumatoid arthritis measures: Disease activity score (DAS), disease activity score-28 (DAS28), rapid assessment of disease activity in rheumatology (RADAR), and rheumatoid arthritis disease activity index (RADAI) , 2003 .

[11]  D. Spiegelhalter,et al.  Bayesian Approaches to Clinical Trials and Health-Care Evaluation: Spiegelhalter/Clinical Trials and Health-Care Evaluation , 2004 .

[12]  Naitee Ting Dose Finding in Drug Development , 2006 .

[13]  Byron J. T. Morgan,et al.  Modelling Population Dynamics , 2014 .

[14]  A. Reynolds,et al.  Potential Relevance of Bell-Shaped and U-Shaped Dose-Responses for the Therapeutic Targeting of Angiogenesis in Cancer , 2010, Dose-response : a publication of International Hormesis Society.

[15]  P. Tak,et al.  MOR103, a human monoclonal antibody to granulocyte–macrophage colony-stimulating factor, in the treatment of patients with moderate rheumatoid arthritis: results of a phase Ib/IIa randomised, double-blind, placebo-controlled, dose-escalation trial , 2014, Annals of the rheumatic diseases.

[16]  N. Thomas,et al.  Meta-Analysis of Clinical Dose–Response in a Large Drug Development Portfolio , 2014 .

[17]  Dose Response Information to Support Drug Registration Step 5 NOTE FOR GUIDANCE ON DOSE RESPONSE INFORMATION , 2006 .

[18]  Michael D. Smith,et al.  Adaptive Bayesian Designs for Dose-Ranging Drug Trials , 2002 .

[19]  M. A. Best Bayesian Approaches to Clinical Trials and Health‐Care Evaluation , 2005 .

[20]  Darren J. Wilkinson,et al.  Bayesian hierarchical modelling for inferring genetic interactions in yeast , 2014, Journal of the Royal Statistical Society. Series C, Applied statistics.

[21]  M. Genovese,et al.  Efficacy and Safety of Belimumab in Patients with Rheumatoid Arthritis: A Phase II, Randomized, Double-blind, Placebo-controlled, Dose-ranging Study , 2013, The Journal of Rheumatology.

[22]  Andrew P Grieve,et al.  Implementation of a Bayesian adaptive design in a proof of concept study , 2006, Pharmaceutical statistics.

[23]  J. Arrowsmith,et al.  Trial Watch: Phase II and Phase III attrition rates 2011–2012 , 2013, Nature Reviews Drug Discovery.

[24]  Byron J. T. Morgan,et al.  Modelling population dynamics : model formulation, fitting and assessment using state-space methods , 2014 .

[25]  S. Kirby,et al.  Fitting Emax models to clinical trial dose-response data when the high dose asymptote is ill defined. , 2014, Pharmaceutical statistics.

[26]  E. Calabrese,et al.  U-shaped dose-responses in biology, toxicology, and public health. , 2001, Annual review of public health.

[27]  Jane Ruth Temple Adaptive designs for dose-finding trials , 2012 .

[28]  H. Leffers,et al.  Dual effects of phytoestrogens result in u-shaped dose-response curves. , 2002, Environmental health perspectives.