Survival analyses from the CENTAUR trial in amyotrophic lateral sclerosis: Evaluating the impact of treatment crossover on outcomes.

INTRODUCTION/AIMS: Trials incorporating placebo-to-active treatment crossover are encouraged in fatal conditions like amyotrophic lateral sclerosis (ALS) but may underestimate active treatment survival benefit. Here, we apply methods for modeling survival without crossover, including the rank-preserving structural failure time model (RPSFTM), to data from the CENTAUR trial of sodium phenylbutyrate and taurursodiol (PB and TURSO) in ALS incorporating both randomized placebo-controlled and open-label extension (OLE) phases.
METHODS: Intent-to-treat (ITT) and RPSFTM survival analyses were performed with final data at a July 2020 cutoff date. Analyses of subgroups based on randomized treatment and OLE phase participation were also performed.
RESULTS: Hazard ratios (95% confidence intervals) of death for PB and TURSO versus participants initially on placebo were 0.57 (0.35-0.92) on ITT analysis and 0.39 (0.17-0.88) in the primary on-treatment RPSFTM analysis (p = .023). Median ITT survival duration for PB and TURSO (25.8 mo) was 6.9 mo longer than placebo (18.9 mo) on ITT analysis and 10.6 mo longer than the median RPSFTM-adjusted survival duration for placebo (15.2 mo). Median survival duration was 18.8 mo longer in the PB and TURSO-randomized subgroup who continued into the OLE phase versus the placebo-randomized subgroup who did not continue into the OLE phase (p < .0001), although OLE phase selection bias may have potentially confounded these results. DISCUSSION: Similar to the prespecified ITT analysis, post hoc analyses adjusting for treatment crossover in CENTAUR showed a significant survival benefit for PB and TURSO. Such methods may provide clinical context for observed survival outcomes in future ALS crossover trials.

INTRODUCTION

Of the two approved disease‐modifying therapies for amyotrophic lateral sclerosis (ALS), only riluzole has been shown to improve survival (specifically, tracheostomy‐free survival) in clinical trials. , Recently, an orally administered, fixed‐dose coformulation of sodium phenylbutyrate and taurursodiol (PB and TURSO) was shown to significantly slow functional decline and to prolong overall survival in people with ALS in the phase 2 CENTAUR trial, which encompassed both randomized placebo‐controlled (NCT03127514) and open‐label extension (OLE) (NCT03488524) phases. Median survival duration was 6.5 mo longer in those originally randomized to PB and TURSO versus placebo in an interim intent‐to‐treat (ITT) analysis spanning both the randomized and OLE phases, but the crossover design of trials like CENTAUR is expected to underestimate the effect of active treatment on overall survival compared with a true placebo. , , , While incorporating an option for receiving active treatment in an OLE phase is of critical importance in the design of trials for a rapidly progressive and fatal disease like ALS, such a design may lead to underestimation of the clinical effect of investigational therapies when comparing earlier initiation of treatment to later initiation, with subsequent effects on both clinical decision‐making and health technology assessments. , ,

Analyses evaluating crossover trial subgroups based on both originally randomized treatment and OLE phase participation provide a simple approach that may provide insight into the potential effect of treatment switching in the control group; however, the results of such analyses are likely biased by differences in prognosis between the OLE phase and non–OLE phase subgroups. , More complex statistical methods for modeling survival benefit in the absence of treatment crossover exist that theoretically eliminate such bias, including rank‐preserving structural failure time models (RPSFTMs). These models are used in the analysis of oncology trials, which frequently incorporate a crossover design based on ethical considerations, and are generally regarded as appropriate methods for assessing the potential effect of treatment crossover on clinical trial survival estimates. ,

In this article, we present the final ITT analyses of survival data from CENTAUR and apply the aforementioned analytic methods to these data post hoc to provide insight into the potential effect of treatment crossover in the placebo group on the estimated survival effect of PB and TURSO.

METHODS

CENTAUR trial

CENTAUR was conducted at 25 Northeast ALS Consortium centers. Protocol approval for both the randomized and OLE phases was provided by a central institutional review board, the Partners Human Research Committee, for all trial sites. Participants provided written informed consent before entering each trial phase. Detailed methods of both trial phases have been published. , Briefly, adults with definite ALS (revised El Escorial criteria) who were ≤18 mo from symptom onset were randomized 2:1 to receive daily PB and TURSO or placebo by mouth or feeding tube. Participants completing the randomized phase were eligible to enroll in the OLE phase and receive daily PB and TURSO for up to 40 mo total. Continuation of a stable dose of riluzole at baseline was permitted, as was initiation or continuation of edaravone during both phases of the trial.

ITT survival analyses

Results from an interim ITT survival analysis performed at a cutoff date of July 20, 2020 (longest follow‐up, 35 mo after randomization) have been previously reported. Updated analyses of the finalized dataset from CENTAUR incorporating updated participant vital status information as of the prior interim ITT analysis, at both the prior July 2020 cutoff date and the date of the final OLE phase participant visit (March 1, 2021), are reported here (Supporting Information Figure S1).

The prespecified ITT survival analyses compared time to death (all‐cause mortality) between participants originally randomized to PB and TURSO versus placebo (Supporting Information Figure S1). Methods for participant vital status determination have been described. For the post hoc subgroup analysis, the study population was broken into four subgroups based on a) whether participants were originally randomized to PB and TURSO versus placebo and b) whether participants enrolled in the OLE phase or not (Supporting Information Figure S1). For both the prespecified ITT survival analysis and post hoc subgroup survival analysis, the hazard ratio (HR) of death in the group originally randomized to PB and TURSO versus the group originally randomized to placebo was estimated using a Cox proportional hazards model with covariates of age at randomization, pre‐baseline Amyotrophic Lateral Sclerosis Functional Rating Scale–Revised (ALSFRS‐R) slope, and baseline ALSFRS‐R total score. Median survival durations and 95% confidence intervals (CIs) were estimated from Kaplan–Meier plots; tests were declared significant if two‐tailed p values were ≤ .05.

RPSFTM analyses

For consistency with the primary analyses, the entire randomized population was used for the RPSFTM analyses. The RPSFTM end point was defined in the same way as for the primary analyses, using all‐cause mortality (Supporting Information Figure S1).

In RPSFTM analyses, the crossover treatment effect duration can be assumed to last from the first dose of active drug until death or censoring (“treatment group” approach) or encompass only the days the participant received active drug (“on‐treatment” approach). The on‐treatment approach was used as the primary analysis herein, with the treatment group approach used as a sensitivity analysis. For the on‐treatment approach, the number of days on PB and TURSO was defined as the date of last dose (in the randomized or OLE phase) minus the date of first dose (in the randomized phase or OLE phase) +1.

An acceleration factor (AF) corresponding with the extent by which active treatment either extended or decreased survival time was determined using G‐estimation and was then used to adjust the survival estimate for participants in the placebo arm who switched to PB and TURSO; additional details are provided in Supporting Information Appendix S1. Given the potential for informative censoring bias, multiple recensoring approaches were applied as sensitivity analyses (see Supporting Information Appendix S1 for details).

Kaplan–Meier plots were produced for the adjusted survival in the arm originally randomized to placebo using both on‐treatment and treatment group approaches. Adjusted survival estimates in the group originally randomized to placebo were compared with observed survival in the group originally randomized to PB and TURSO from the prespecified ITT analysis (Supporting Information Figure S1). As for the primary analyses, HRs were estimated using a Cox proportional hazards model with covariates of age at randomization, pre‐baseline ALSFRS‐R slope, and baseline ALSFRS‐R total score. Confidence intervals and p values were based on the ITT p value.

RESULTS

Participants

A total of 137 participants were randomized in CENTAUR (PB and TURSO, n = 89; placebo, n = 48); 90 of 98 eligible participants continued into the OLE phase, including 71% (34/48) of participants originally randomized to placebo. In this final analysis, vital status was obtainable for all but one participant randomized in CENTAUR; this participant was censored at the date of last follow‐up visit. Additional participant disposition as well as baseline data have been previously published.

ITT analyses

In the final overall survival analysis encompassing all randomized participants at the July 2020 cutoff date, median (95% CI) survival duration was 25.8 mo (19.0 mo–not reached [NR]) in the group originally randomized to PB and TURSO and 18.9 (13.5–28.7) mo in the group originally randomized to placebo (6.9‐mo difference), with an HR of 0.57 (Figure 1); mean PB and TURSO exposure durations were 10.2 mo in the group originally randomized to PB and TURSO and 4.6 mo in the group originally randomized to placebo (all in the OLE phase). Results of the survival analysis at the March 2021 cutoff date were concordant, showing a significantly lower hazard of death and longer median survival duration in the group originally randomized to PB and TURSO (Supporting Information Figure S2).

FIGURE 1

Kaplan–Meier analyses: Prespecified ITT survival analysis (red and blue lines) and rank‐preserving structural failure time model–adjusted survival (yellow line, recensoring acceleration factor only; on‐treatment approach). All analyses incorporate final data from the July 20, 2020, cutoff date, including updated participant vital status information as of the prior published interim ITT analysis at this cutoff date

On‐treatment RPSFTM analyses of data were performed through the July 2020 cutoff date and applying recensoring of the AF. Compared with the observed estimate of 25.8 mo for the group originally randomized to PB and TURSO in the prespecified ITT analysis, RPSFTM‐adjusted median survival duration was 15.2 mo in the group originally randomized to placebo (difference of 10.6 mo); HR was 0.39 (Figure 1). Similar results were obtained with full recensoring, with an HR of 0.44 (Supporting Information Figure S3), consistent with a beneficial effect (AF <1) of PB and TURSO treatment on overall survival. AF could not be estimated in assessments of on‐treatment RPSFTM without applying recensoring. Sensitivity analyses using the treatment group RPSFTM yielded similar results to the on‐treatment analyses (Supporting Information Figure S4). The same on‐treatment and treatment group RPSFTM analyses were performed on data through the March 2021 cutoff date and also yielded consistent results (Supporting Information Figures S2 and S5, respectively).

Post hoc subgroup analysis

Results of the post hoc survival analysis of subgroups based on originally randomized treatment and OLE phase participation through the July 2020 cutoff date are shown in Figure 2. Median survival duration decreased and hazard of death increased with decreasing PB and TURSO exposure duration (all p ≤ .01 compared with the subgroup of participants who were originally randomized to PB and TURSO and continued into the OLE phase). Median survival duration was 18.8 mo longer in the group who were originally randomized to PB and TURSO and continued into the OLE phase (33.6 mo) than in the group who were originally randomized to placebo and did not cross over to active treatment in the OLE phase (14.8 mo; p < .0001). Among the other subgroups, median survival duration was longer in the group who were originally randomized to placebo and continued into the OLE phase, who had a longer mean PB and TURSO exposure duration, compared with the group who were originally randomized to PB and TURSO but did not continue into the OLE phase. Subgroup analysis at the March 2021 cutoff date yielded similar results (Supporting Information Figure S6).

FIGURE 2

Cox proportional hazards analysis of time to death: Subgroups based on originally randomized treatment and OLE phase participation (July 2020 cutoff date). aCompared with subgroup originally randomized to sodium phenylbutyrate and taurursodiol that entered the OLE phase

DISCUSSION

In our RPSFTM analyses of data from CENTAUR, statistical methods that adjust for treatment crossover in the group originally randomized to placebo suggested a greater survival benefit with PB and TURSO use than seen in the ITT analysis.

Post hoc assessment of subgroups from CENTAUR based on randomization group as well as enrollment in the OLE phase demonstrated that earlier and longer exposure to PB and TURSO was associated with longer survival estimates. The finding of longer median survival duration in the subgroup of participants who were originally randomized to placebo and crossed over to active treatment in the OLE phase versus the subgroup who were originally randomized to PB and TURSO but did not continue into the OLE phase further supports an association between PB and TURSO exposure duration and survival outcome in CENTAUR.

It is important to note that the post hoc subgroup analysis applied herein has limitations as the subgroups are small, and potential confounding differences among the groups were not controlled for, although covariates of age at randomization, pre‐baseline ALSFRS‐R slope, and baseline ALSFRS‐R total score attempted to control for some of this bias. In the RPSFTM analysis, it was not possible to estimate the AF for the RPSFTM without recensoring using the on‐treatment duration of effect assumption. No factor could be selected that balanced the counterfactual survival (ie, the model‐estimated survival time in the absence of PB and TURSO) without recensoring between the treatment arms. It may be due to a prognostic imbalance at baseline between randomized arms, although this seems unlikely as baseline characteristics were generally well balanced and the AF could be estimated from other models. Alternatively, it may be that the on‐treatment duration of effect assumption is not suitable for this study, but this was judged to be a more plausible assumption than the treatment group approach. Similar results were seen for the analyses using both duration of effect assumptions, suggesting that varying this assumption does not have a large impact.

Another important point about the RPSFTM method is that it assumes a common treatment effect (ie, exposure‐response is the same, no matter when the treatment is received). This assumption may be unreliable in degenerative conditions such as ALS, although the time between randomization and crossover in CENTAUR was relatively short (only 6 mo). In addition, the subgroup of participants who were originally randomized to placebo and continued into the OLE phase had a lower risk of death and longer median survival than the subgroup who were originally randomized to active treatment and discontinued after the randomized phase, indicating that exposure to drug and not current severity of disease was more important to survival. Finally, PB and TURSO targets neuronal death, which is expected to be relevant at all stages of disease represented in the trial, although again, it cannot be known for certain whether treatment effects at earlier versus later time points might differ. Other methods such as inverse probability of censoring weighting and two‐stage models are often used in oncology to adjust overall survival for switching and were considered here; however, these methods are not suitable when most participants switch, as was the case in the CENTAUR trial, in which only three of the placebo‐randomized participants who were eligible to enroll in the OLE phase did not do so.

In conclusion, ALS clinical trials can gather robust survival data while incorporating study design elements, such as OLE phases, that are critical for facilitating access to investigational therapies for people with ALS. However, the inherent inability to sustain a true placebo group long enough to assess survival in such studies may lead to underestimation of the clinical effect of the therapies under investigation when having to resort to comparing earlier initiation of treatment to later initiation. By adjusting survival estimates in the presence of treatment switching, methods such as RPSFTM and subgroup analyses may provide additional clinical context beyond the observed survival outcomes in ALS trials incorporating this critical crossover design that may be informative for patients and other stakeholders. In post hoc analyses of the CENTAUR trial, these two methods yielded a 10.6‐ and 18.8‐mo adjusted survival benefit of PB and TURSO, respectively. Additional data from the ongoing phase 3 PHOENIX trial (NCT05021536) and clinical experience will provide further information regarding the effect of PB and TURSO on survival in people with ALS.

AUTHOR CONTRIBUTIONS

Dr Paganoni had full access to the trial data and takes full responsibility for the integrity of the data and the accuracy of the data analysis. Dr Paganoni, Ms Watkins, Mr Cawson, and Drs Timmons, Manuel, and Cudkowicz were involved in the conceptualization of the research described in the article. Ms Watkins, Mr Cawson, Drs Hendrix and Dickson, and Mr Knowlton curated and analyzed the data; oversaw the methodology and provided software for the analyses described in the article; and contributed to the validation and visualization of the analysis results. Drs Paganoni and Cudkowicz contributed to data acquisition in CENTAUR and provided oversight and leadership for the planning and execution of the research described in the article. Drs Timmons and Manuel provided project administration oversight. Ms Watkins, Mr Cawson, Drs Hendrix and Dickson, Mr Knowlton, and Drs Timmons and Manuel contributed to provision of resources, including study materials and analysis tools. Ms Watkins, Mr Cawson, and Dr Timmons contributed to the first‐draft development of the article, and all authors critically reviewed and revised the article at each stage of development and approved the final draft for submission.

CONFLICTS OF INTEREST

Sabrina Paganoni reports research grants from Amylyx Pharmaceuticals, Revalesio Corporation, Alector Therapeutics, UCB, Biohaven Pharmaceuticals, Clene Nanomedicine, Prilenia Therapeutics, Seelos Therapeutics, The ALS Association, the American Academy of Neurology, the Centers for Disease Control, ALS Finding a Cure®, the Salah Foundation, the Spastic Paraplegia Foundation, the Muscular Dystrophy Association, I AM ALS, Tambourine, Target ALS, Columbia University, and the Cullen Education and Research Fund, and is a site principle investigator for studies funded by Alector Therapeutics, Cytokinetics, and Novus (formerly Anelixis) Therapeutics. Dr Paganoni also reports institutional consulting agreements with Amylyx Pharmaceuticals, Frequency Therapeutics, and SOLA Biosciences; personal consulting agreements with Cytokinetics, Arrowhead Pharmaceuticals, and Orthogonal Neuroscience; and honoraria from Medscape. Claire Watkins and Matthew Cawson report fees via the Maple Health Group from Amylyx Pharmaceuticals for conducting some of the analyses described in the submitted work. Suzanne Hendrix, Samuel P. Dixon, and Newman Knowlton are employees of Pentara Corporation. Jamie Timmons and Machelle Manuel are full‐time employees of Amylyx Pharmaceuticals. Merit Cudkowicz reports consulting fees from Faze, Regeneron, AB Sciences, Avexis, MT Pharma, Revalasio, Cytokinetics, Disarm, ALS Pharma, Immunity Pharma, Wave, Sunovian, Transposon, Quralis, Helixsmith, Locust Walk, and RRD and is a board member for Praxis, all outside the submitted work.

acceleration factor

amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis Functional Rating Scale–Revised

confidence interval

hazard ratio

intent‐to‐treat

not applicable

not reached

open‐label extension

sodium phenylbutyrate and taurursodiol

rank‐preserving structural failure time model

ETHICAL PUBLICATION STATEMENT

The authors confirm that we have read the Journal's position on issues involved in ethical publication and affirm that this report is consistent with those guidelines.

Figure S1 Schematic of survival analyses conducted on CENTAUR dataset. Time to all‐cause death was assessed over the periods spanning randomization through each cutoff date and compared among the groups shown for each analysis (table). The intent‐to‐treat (ITT) analysis was prespecified, while the other analyses were performed post hoc. aThe ITT population consisted of all participants randomized in CENTAUR who received ≥1 dose of trial drug in the randomized phase. OLE, open‐label extension; PB and TURSO, sodium phenylbutyrate and taurursodiol; RPSFTM, rank‐preserving structural failure time model.

Click here for additional data file.

Figure S2 Kaplan–Meier analyses: prespecified intent‐to‐treat (ITT) survival analysis (red and blue lines) and rank‐preserving structural failure time model (RPSFTM)–adjusted survival in the originally randomized placebo group (yellow line, recensoring acceleration factor [AF] only; purple line, full recensoring; on‐treatment approach). All analyses incorporate matured data from the March 1, 2021, cutoff date coinciding with the final participant visit in CENTAUR. The AF used to adjust survival for switching could not be estimated in the on‐treatment RPSFTM without applying recensoring. HR, hazard ratio; PB and TURSO, sodium phenylbutyrate and taurursodiol.

Click here for additional data file.

Figure S3 Kaplan–Meier analyses: prespecified intent‐to‐treat (ITT) survival analysis (red and blue lines) and rank‐preserving structural failure time model (RPSFTM)–adjusted survival (purple line, full recensoring; on‐treatment approach). All analyses incorporate final data from the July 20, 2020, cutoff date, including updated participant vital status information as of the prior published interim ITT analysis at this cutoff date [Paganoni S et al. Muscle Nerve. 2021;63(1):31–39]. The acceleration factor used to adjust survival for switching could not be estimated in the on‐treatment RPSFTM without applying recensoring. HR, hazard ratio; PB and TURSO, sodium phenylbutyrate and taurursodiol; RPSFTM, rank‐preserving structural failure time model.

Click here for additional data file.

Figure S4 Kaplan–Meier analysis: prespecified intent‐to‐treat survival analysis (red and blue lines) and rank‐preserving structural failure time model–adjusted survival (green line, no recensoring; yellow line, recensoring acceleration factor only; purple line, full recensoring; treatment group approach). All analyses incorporate final data from the July 20, 2020, cutoff date. AF, acceleration factor; HR, hazard ratio; ITT, intent‐to‐treat; PB and TURSO, sodium phenylbutyrate and taurursodiol; RPSFTM, rank‐preserving structural failure time model.

Click here for additional data file.

Figure S5 Kaplan–Meier analysis: prespecified intent‐to‐treat survival analysis (red and blue lines) and rank‐preserving structural failure time model–adjusted survival (green line, no recensoring; yellow line, recensoring acceleration factor only; purple line, full recensoring; treatment group approach). All analyses incorporate final data from the March 1, 2021, cutoff date. AF, acceleration factor; HR, hazard ratio; ITT, intent‐to‐treat; PB and TURSO, sodium phenylbutyrate and taurursodiol; RPSFTM, rank‐preserving structural failure time model.

Click here for additional data file.

Figure S6 Cox proportional hazards analysis of time to death: subgroups based on originally randomized treatment and open‐label extension (OLE) phase participation (March 2021 cutoff date). aCompared with subgroup originally randomized to sodium phenylbutyrate and taurursodiol that entered the OLE phase. HR, hazard ratio; NA, not applicable; NR, not reached; PB and TURSO, sodium phenylbutyrate and taurursodiol.

Click here for additional data file.

Appendix S1 Supporting information

Click here for additional data file.