Bayesian adaptive model selection design for optimal biological dose finding in phase I/II clinical trials.

Identification of the optimal dose presents a major challenge in drug development with molecularly targeted agents, immunotherapy, as well as chimeric antigen receptor T-cell treatments. By casting dose finding as a Bayesian model selection problem, we propose an adaptive design by simultaneously incorporating the toxicity and efficacy outcomes to select the optimal biological dose (OBD) in phase I/II clinical trials. Without imposing any parametric assumption or shape constraint on the underlying dose-response curves, we specify curve-free models for both the toxicity and efficacy endpoints to determine the OBD. By integrating the observed data across all dose levels, the proposed design is coherent in dose assignment and thus greatly enhances efficiency and accuracy in pinning down the right dose. Not only does our design possess a completely new yet flexible dose-finding framework, but it also has satisfactory and robust performance as demonstrated by extensive simulation studies. In addition, we show that our design enjoys desirable coherence properties, while most of existing phase I/II designs do not. We further extend the design to accommodate late-onset outcomes which are common in immunotherapy. The proposed design is exemplified with a phase I/II clinical trial in chronic lymphocytic leukemia.

A Bayesian latent-subgroup platform design for dose optimization.

The US Food and Drug Administration launched Project Optimus to reform the dose optimization and dose selection paradigm in oncology drug development, calling for the paradigm shift from finding the maximum tolerated dose to the identification of optimal biological dose (OBD). Motivated by a real-world drug development program, we propose a master-protocol-based platform trial design to simultaneously identify OBDs of a new drug, combined with standards of care or other novel agents, in multiple indications. We propose a Bayesian latent subgroup model to accommodate the treatment heterogeneity across indications, and employ Bayesian hierarchical models to borrow information within subgroups. At each interim analysis, we update the subgroup membership and dose-toxicity and -efficacy estimates, as well as the estimate of the utility for risk-benefit tradeoff, based on the observed data across treatment arms to inform the arm-specific decision of dose escalation and de-escalation and identify the OBD for each arm of a combination partner and an indication. The simulation study shows that the proposed design has desirable operating characteristics, providing a highly flexible and efficient way for dose optimization. The design has great potential to shorten the drug development timeline, save costs by reducing overlapping infrastructure, and speed up regulatory approval.

Adaptive designs in public health: Vaccine and cluster randomized trials go Bayesian.

Clinical trials in public health-particularly those conducted in low- and middle-income countries-often involve communicable and non-communicable diseases with high disease burden and unmet needs. Trials conducted in these regions often are faced with resource limitations, so improving the efficiencies of these trials is critical. Adaptive trial designs have the potential to save trial time and resources and reduce the number of patients receiving ineffective interventions. In this paper, we provide a detailed account of the implementation of vaccine and cluster randomized trials within the framework of Bayesian adaptive trials, with emphasis on computational efficiency and flexibility with regard to stopping rules and allocation ratios. We offer an educated approach to selecting prior distributions and a data-driven empirical Bayes method for plug-in estimates for nuisance parameters.

A Bayesian adaptive design approach for stepped-wedge cluster randomized trials.

BACKGROUND: The Bayesian group sequential design has been applied widely in clinical studies, especially in Phase II and III studies. It allows early termination based on accumulating interim data. However, to date, there lacks development in its application to stepped-wedge cluster randomized trials, which are gaining popularity in pragmatic trials conducted by clinical and health care delivery researchers. METHODS: We propose a Bayesian adaptive design approach for stepped-wedge cluster randomized trials, which makes adaptive decisions based on the predictive probability of declaring the intervention effective at the end of study given interim data. The Bayesian models and the algorithms for posterior inference and trial conduct are presented. RESULTS: We present how to determine design parameters through extensive simulations to achieve desired operational characteristics. We further evaluate how various design factors, such as the number of steps, cluster size, random variability in cluster size, and correlation structures, impact trial properties, including power, type I error, and the probability of early stopping. An application example is presented. CONCLUSION: This study presents the incorporation of Bayesian adaptive strategies into stepped-wedge cluster randomized trials design. The proposed approach provides the flexibility to stop the trial early if substantial evidence of efficacy or futility is observed, improving the flexibility and efficiency of stepped-wedge cluster randomized trials.

Adaptive phase I-II clinical trial designs identifying optimal biological doses for targeted agents and immunotherapies.

Targeted agents and immunotherapies have revolutionized cancer treatment, offering promising options for various cancer types. Unlike traditional therapies the principle of "more is better" is not always applicable to these new therapies due to their unique biomedical mechanisms. As a result, various phase I-II clinical trial designs have been proposed to identify the optimal biological dose that maximizes the therapeutic effect of targeted therapies and immunotherapies by jointly monitoring both efficacy and toxicity outcomes. This review article examines several innovative phase I-II clinical trial designs that utilize accumulated efficacy and toxicity outcomes to adaptively determine doses for subsequent patients and identify the optimal biological dose, maximizing the overall therapeutic effect. Specifically, we highlight three categories of phase I-II designs: efficacy-driven, utility-based, and designs incorporating multiple efficacy endpoints. For each design, we review the dose-outcome model, the definition of the optimal biological dose, the dose-finding algorithm, and the software for trial implementation. To illustrate the concepts, we also present two real phase I-II trial examples utilizing the EffTox and ISO designs. Finally, we provide a classification tree to summarize the designs discussed in this article.

The score-goldilocks design for phase 3 clinical trials.

In this paper, we propose a new Bayesian adaptive design, score-goldilocks design, which has the same algorithmic idea as goldilocks design. The score-goldilocks design leads to a uniform formula for calculating the probability of trial success for different endpoint trials by using the normal approximation. The simulation results show that the score-goldilocks design is not only very similar to the goldilocks design in terms of operating characteristics such as type 1 error, power, average sample size, probability of stop for futility, and probability of early stop for success, but also greatly saves the calculation time and improves the operation efficiency.

T3 + 3: 3 + 3 Design With Delayed Outcomes.

Delayed outcome is common in phase I oncology clinical trials. It causes logistic difficulty, wastes resources, and prolongs the trial duration. This article investigates this issue and proposes the time-to-event 3 + 3 (T3 + 3) design, which utilizes the actual follow-up time for at-risk patients with pending toxicity outcomes. The T3 + 3 design allows continuous accrual without unnecessary trial suspension and is costless and implementable with pretabulated dose decision rules. Besides, the T3 + 3 design uses the isotonic regression to estimate the toxicity rates across dose levels and therefore can accommodate for any targeted toxicity rate for maximum tolerated dose (MTD). It dramatically facilitates the trial preparation and conduct without intensive computation and statistical consultation. The extension to other algorithm-based phase I dose-finding designs (e.g., i3 + 3 design) is also studied. Comprehensive computer simulation studies are conducted to investigate the performance of the T3 + 3 design under various dose-toxicity scenarios. The results confirm that the T3 + 3 design substantially shortens the trial duration compared with the conventional 3 + 3 design and yields much higher accuracy in MTD identification than the rolling six design. In summary, the T3 + 3 design addresses the delayed outcome issue while keeping the desirable features of the 3 + 3 design, such as simplicity, transparency, and costless implementation. It has great potential to accelerate early-phase drug development.

An alternative metric for evaluating the potential patient benefit of response-adaptive randomization procedures.

When planning a two-arm group sequential clinical trial with a binary primary outcome that has severe implications for quality of life (e.g., mortality), investigators may strive to find the design that maximizes in-trial patient benefit. In such cases, Bayesian response-adaptive randomization (BRAR) is often considered because it can alter the allocation ratio throughout the trial in favor of the treatment that is currently performing better. Although previous studies have recommended using fixed randomization over BRAR based on patient benefit metrics calculated from the realized trial sample size, these previous comparisons have been limited by failures to hold type I and II error rates constant across designs or consider the impacts on all individuals directly affected by the design choice. In this paper, we propose a metric for comparing designs with the same type I and II error rates that reflects expected outcomes among individuals who would participate in the trial if enrollment is open when they become eligible. We demonstrate how to use the proposed metric to guide the choice of design in the context of two recent trials in persons suffering out of hospital cardiac arrest. Using computer simulation, we demonstrate that various implementations of group sequential BRAR offer modest improvements with respect to the proposed metric relative to conventional group sequential monitoring alone.

BASIC: A Bayesian adaptive synthetic-control design for phase II clinical trials.

BACKGROUND: Randomized controlled trials are considered the gold standard for evaluating experimental treatments but often require large sample sizes. Single-arm trials require smaller sample sizes but are subject to bias when using historical control data for comparative inferences. This article presents a Bayesian adaptive synthetic-control design that exploits historical control data to create a hybrid of a single-arm trial and a randomized controlled trial. METHODS: The Bayesian adaptive synthetic control design has two stages. In stage 1, a prespecified number of patients are enrolled in a single arm given the experimental treatment. Based on the stage 1 data, applying propensity score matching and Bayesian posterior prediction methods, the usefulness of the historical control data for identifying a pseudo sample of matched synthetic-control patients for making comparative inferences is evaluated. If a sufficient number of synthetic controls can be identified, the single-arm trial is continued. If not, the trial is switched to a randomized controlled trial. The performance of The Bayesian adaptive synthetic control design is evaluated by computer simulation. RESULTS: The Bayesian adaptive synthetic control design achieves power and unbiasedness similar to a randomized controlled trial but on average requires a much smaller sample size, provided that the historical control data patients are sufficiently comparable to the trial patients so that a good number of matched controls can be identified in the historical control data. Compared to a single-arm trial, The Bayesian adaptive synthetic control design yields much higher power and much smaller bias. CONCLUSION: The Bayesian adaptive synthetic-control design provides a useful tool for exploiting historical control data to improve the efficiency of single-arm phase II clinical trials, while addressing the problem of bias when comparing trial results to historical control data. The proposed design achieves power similar to a randomized controlled trial but may require a substantially smaller sample size.

Enhancement of Bayesian optimal interval design by accounting for overdose and underdose errors trade-offs.

Model-assisted designs, a new class of dose-finding designs for determining the maximum tolerated dose (MTD), model only the dose-limiting toxicity (DLT) data observed at the current dose based on a simple binomial model and offer the boundaries of DLT for the determination of dose escalation, retention, or de-escalation before beginning the trials. The boundaries for dose-escalation and de-escalation decisions are relevant to the operating characteristics of the design. The well-known model-assisted design, Bayesian Optimal Interval (BOIN), selects these boundaries to minimize the probability of incorrect decisions at each dose allocation but does not distinguish between overdose and underdose allocations caused by incorrect decisions when calculating the probability of incorrect decisions. Distinguishing between overdose and underdose based on the decision error in the BOIN design is expected to increase the accuracy of MTD determination. In this study, we extended the BOIN design to account for the decision probabilities of incorrect overdose and underdose allocations separately. To minimize the two probabilities simultaneously, we propose utilizing multiple objective optimizations and formulating an approach for determining the boundaries for dose escalation and de-escalation. Comprehensive simulation studies using fixed and randomly generated scenarios of DLT probability demonstrated that the proposed method is superior or comparable to existing interval designs, along with notably better operating characteristics of the proposed method.

Bayesian optimal phase II designs with dual-criterion decision making.

The conventional phase II trial design paradigm is to make the go/no-go decision based on the hypothesis testing framework. Statistical significance itself alone, however, may not be sufficient to establish that the drug is clinically effective enough to warrant confirmatory phase III trials. We propose the Bayesian optimal phase II trial design with dual-criterion decision making (BOP2-DC), which incorporates both statistical significance and clinical relevance into decision making. Based on the posterior probability that the treatment effect reaches the lower reference value (statistical significance) and the clinically meaningful value (clinical significance), BOP2-DC allows for go/consider/no-go decisions, rather than a binary go/no-go decision. BOP2-DC is highly flexible and accommodates various types of endpoints, including binary, continuous, time-to-event, multiple, and coprimary endpoints, in single-arm and randomized trials. The decision rule of BOP2-DC is optimized to maximize the probability of a go decision when the treatment is effective or minimize the expected sample size when the treatment is futile. Simulation studies show that the BOP2-DC design yields desirable operating characteristics. The software to implement BOP2-DC is freely available at www.trialdesign.org.

A Bayesian piecewise exponential phase II design for monitoring a time-to-event endpoint.

A robust Bayesian design is presented for a single-arm phase II trial with an early stopping rule to monitor a time to event endpoint. The assumed model is a piecewise exponential distribution with non-informative gamma priors on the hazard parameters in subintervals of a fixed follow up interval. As an additional comparator, we also define and evaluate a version of the design based on an assumed Weibull distribution. Except for the assumed models, the piecewise exponential and Weibull model based designs are identical to an established design that assumes an exponential event time distribution with an inverse gamma prior on the mean event time. The three designs are compared by simulation under several log-logistic and Weibull distributions having different shape parameters, and for different monitoring schedules. The simulations show that, compared to the exponential inverse gamma model based design, the piecewise exponential design has substantially better performance, with much higher probabilities of correctly stopping the trial early, and shorter and less variable trial duration, when the assumed median event time is unacceptably low. Compared to the Weibull model based design, the piecewise exponential design does a much better job of maintaining small incorrect stopping probabilities in cases where the true median survival time is desirably large.

A Bayesian phase I/II design to determine subgroup-specific optimal dose for immunotherapy sequentially combined with radiotherapy.

Sequential administration of immunotherapy following radiotherapy (immunoRT) has attracted much attention in cancer research. Due to its unique feature that radiotherapy upregulates the expression of a predictive biomarker for immunotherapy, novel clinical trial designs are needed for immunoRT to identify patient subgroups and the optimal dose for each subgroup. In this article, we propose a Bayesian phase I/II design for immunotherapy administered after standard-dose radiotherapy for this purpose. We construct a latent subgroup membership variable and model it as a function of the baseline and pre-post radiotherapy change in the predictive biomarker measurements. Conditional on the latent subgroup membership of each patient, we jointly model the continuous immune response and the binary efficacy outcome using plateau models, and model toxicity using the equivalent toxicity score approach to account for toxicity grades. During the trial, based on accumulating data, we continuously update model estimates and adaptively randomize patients to admissible doses. Simulation studies and an illustrative trial application show that our design has good operating characteristics in terms of identifying both patient subgroups and the optimal dose for each subgroup.

AIDE: Adaptive intrapatient dose escalation designs to accelerate Phase I clinical trials.

Designing Phase I clinical trials is challenging when accrual is slow or sample size is limited. The corresponding key question is: how to efficiently and reliably identify the maximum tolerated dose (MTD) using a sample size as small as possible? We propose model-assisted and model-based designs with adaptive intrapatient dose escalation (AIDE) to address this challenge. AIDE is adaptive in that the decision of conducting intrapatient dose escalation depends on both the patient's individual safety data, as well as other enrolled patient's safety data. When both data indicate reasonable safety, a patient may perform intrapatient dose escalation, generating toxicity data at more than one dose. This strategy not only provides patients the opportunity to receive higher potentially more effective doses, but also enables efficient statistical learning of the dose-toxicity profile of the treatment, which dramatically reduces the required sample size. Simulation studies show that the proposed designs are safe, robust, and efficient to identify the MTD with a sample size that is substantially smaller than conventional interpatient dose escalation designs. Practical considerations are provided and R code for implementing AIDE is available upon request.

Biomarker-based precision dose finding for immunotherapy combined with radiotherapy.

Recent success of sequential administration of immunotherapy following radiotherapy (RT), often referred to as immunoRT, has sparked the urgent need for novel clinical trial designs to accommodate the unique features of immunoRT. For this purpose, we propose a Bayesian phase I/II design for immunotherapy administered after standard-dose RT to identify the optimal dose that is personalized for each patient according to his/her measurements of PD-L1 expression at baseline and post-RT. We model the immune response, toxicity, and efficacy as functions of dose and patient's baseline and post-RT PD-L1 expression profile. We quantify the desirability of the dose using a utility function and propose a two-stage dose-finding algorithm to find the personalized optimal dose. Simulation studies show that our proposed design has good operating characteristics, with a high probability of identifying the personalized optimal dose.

Bayesian adaptive trial design for a continuous biomarker with possibly nonlinear or nonmonotone prognostic or predictive effects.

As diseases like cancer are increasingly understood on a molecular level, clinical trials are being designed to reveal or validate subpopulations in which an experimental therapy has enhanced benefit. Such biomarker-driven designs, particularly "adaptive enrichment" designs that initially enroll an unselected population and then allow for later restriction of accrual to "marker-positive" patients based on interim results, are increasingly popular. Many biomarkers of interest are naturally continuous, however, and most existing design approaches either require upfront dichotomization or force monotonicity through algorithmic searches for a single marker threshold, thereby excluding the possibility that the continuous biomarker has a nondisjoint and truly nonlinear or nonmonotone prognostic relationship with outcome or predictive relationship with treatment effect. To address this, we propose a novel trial design that leverages both the actual shapes of any continuous marker effects (both prognostic and predictive) and their corresponding posterior uncertainty in an adaptive decision-making framework. At interim analyses, this marker knowledge is updated and overall or marker-driven decisions are reached such as continuing enrollment to the next interim analysis or terminating early for efficacy or futility. Using simulations and patient-level data from a multi-center Children's Oncology Group trial in Acute Lymphoblastic Leukemia, we derive the operating characteristics of our design and compare its performance to a traditional approach that identifies and applies a dichotomizing marker threshold.

A Bayesian phase I/II platform design for co-developing drug combination therapies for multiple indications.

There is a growing trend to combine a new targeted or immunotherapy agent with the cancer-specific standard of care to treat different types of cancers. We propose a master-protocol-based, Bayesian phase I/II platform design to co-develop combination (BPCC) therapies in multiple indications. Under the BPCC design, only a single master protocol is needed, and the combined drug is evaluated in different indications in a concurrent or staggered fashion. For each indication, we jointly model dose-toxicity and -efficacy relationships and employ Bayesian hierarchical models to borrow information across them for more efficient indication-specific decision-making. To account for the characteristic of targeted or immunotherapy agents that their efficacy may not monotonically increase with the dose, and often plateau at high doses, we use the utility to quantify the risk-benefit tradeoff of the treatment. At each interim, we update the toxicity and efficacy model, as well as the estimate of the utility, based on the observed data across indications to inform the indication-specific decision of dose escalation and de-escalation and identify the optimal biological dose for each indication. Simulation study shows that the BPCC design has desirable operating characteristics, and that it provides an efficient approach to accelerate the development of combination therapies.

TITE-BOIN12: A Bayesian phase I/II trial design to find the optimal biological dose with late-onset toxicity and efficacy.

In the era of immunotherapies and targeted therapies, the focus of early phase clinical trials has shifted from finding the maximum tolerated dose to identifying the optimal biological dose (OBD), which maximizes the toxicity-efficacy trade-off. One major impediment to using adaptive designs to find OBD is that efficacy or/and toxicity are often late-onset, hampering the designs' real-time decision rules for treating new patients. To address this issue, we propose the model-assisted TITE-BOIN12 design to find OBD with late-onset toxicity and efficacy. As an extension of the BOIN12 design, the TITE-BOIN12 design also uses utility to quantify the toxicity-efficacy trade-off. We consider two approaches, Bayesian data augmentation and an approximated likelihood method, to enable real-time decision making when some patients' toxicity and efficacy outcomes are pending. Extensive simulations show that, compared to some existing designs, TITE-BOIN12 significantly shortens the trial duration while having comparable or higher accuracy to identify OBD and a lower risk of overdosing patients. To facilitate the use of the TITE-BOIN12 design, we develop a user-friendly software freely available at http://www.trialdesign.org.

Optimal adaptive promising zone designs.

We develop optimal decision rules for sample size re-estimation in two-stage adaptive group sequential clinical trials. It is usual for the initial sample size specification of such trials to be adequate to detect a realistic treatment effect δa with good power, but not sufficient to detect the smallest clinically meaningful treatment effect δmin . Moreover it is difficult for the sponsors of such trials to make the up-front commitment needed to adequately power a study to detect δmin . It is easier to justify increasing the sample size if the interim data enter a so-called "promising zone" that ensures with high probability that the trial will succeed. We have considered promising zone designs that optimize unconditional power and promising zone designs that optimize conditional power and have discussed the tension that exists between these two objectives. Where there is reluctance to base the sample size re-estimation rule on the parameter δmin we propose a Bayesian option whereby a prior distribution is assigned to the unknown treatment effect δ , which is then integrated out of the objective function with respect to its posterior distribution at the interim analysis.

BIPSE: A biomarker-based phase I/II design for immunotherapy trials with progression-free survival endpoint.

A Bayesian biomarker-based phase I/II design (BIPSE) is presented for immunotherapy trials with a progression-free survival (PFS) endpoint. The objective is to identify the subgroup-specific optimal dose, defined as the dose with the best risk-benefit tradeoff in each biomarker subgroup. We jointly model the immune response, toxicity outcome, and PFS with information borrowing across subgroups. A plateau model is used to describe the marginal distribution of the immune response. Conditional on the immune response, we model toxicity using probit regression and model PFS using the mixture cure rate model. During the trial, based on the accumulating data, we continuously update model estimates and adaptively randomize patients to doses with high desirability within each subgroup. Simulation studies show that the BIPSE design has desirable operating characteristics in selecting the subgroup-specific optimal doses and allocating patients to those optimal doses, and outperforms conventional designs.