Serological biomarkers of joint tissue turnover predict tocilizumab response at baseline.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by polyarticular joint inflammation resulting in massive tissue turnover. The turnover is partly mediated by an up-regulation of proteolytic enzymes, such as matrix metalloproteinases (MMPs). Matrix metalloproteinase 3 is 1 of the MMPs responsible for the degradation of the extracellular matrix (ECM). The MMP-mediated degradation of the main joint ECM proteins (eg, types I and III collagen) results in the release of specific biomarkers such as the connective tissue biomarkers C1M and C3M,[4,5] known as protein fingerprints. These biomarkers are direct measures of changes to the tissue affected by the disease, in contrast to measurement of acute reactants such as C-reactive protein (CRP) or interleukin 6 (IL-6), which are upstream of tissue changes. Thus, protein fingerprint biomarkers may be more sensitive tools for measuring disease changes or changes caused by intervention. Protein fingerprint biomarkers have been associated with disease progression and response to therapy.[7,8] CRPM is a protein fingerprint formed through degradation of CRP. In response to IL-6, CRP is secreted by the liver as an acute phase reactant. C-reactive proteinaccumulates in inflamed tissue, where it is degraded by MMPs, resulting in the release of CRPM. The ratio of C3M to CRPM may depict MMP3 is 1 of the MMPs responsible for the degradation of the ECM; its expression is highly elevated in the affected joint and may therefore be a relevant marker of proteolytic activity.

Tocilizumab (TCZ) is approved in 2 doses for intravenous infusion:4 and 8 mg/kg. Although both doses provide structural progression and symptomatic relief, 8 mg/kg generally affords a higher level of response. Composite quantifiable measures depending on CRP were more reduced in 8 mg/kg compared with 4 mg/kg.[12,13] As there are more adverse events in the higher dose, identification of those patients who respond most optimally to 4 mg/kg would significantly improve the benefit-to-risk assessment. The aim of present study was to identify responders to 4 mg/kg TCZ by measuring protein fingerprints at baseline.

METHODS

Study Design and Serum Samples

Fasting RA patients’ serum samples (n = 200) from the 4 mg/kg TCZ treatment arm of the LITHE phase III study were analyzed for the following protein fingerprint biomarkers: C1M, C3M, and CRPM were measured by manual competitive enzyme-linked immunosorbent assays (Nordic Bioscience, Herlev, Denmark), and serum total MMP-3 was measured by a 2-site enzyme-linked immunosorbent assay (Quantikine; R&D Systems, Lille, France). Five additional markers were measured; however, these showed no discriminative power when plotted (area under the receiver operating characteristic curve). Treatment response was recorded according to the American College of Rheumatology criteria for 50% improvement (ACR50) at week 52. Baseline patient characteristics are summarized in the Table. The study was approved by the ethics committee at each recruiting institution and was conducted according to the Principles of Good Clinical Practice and according to the Declaration of Helsinki.

Statistical Analysis

Biomarker data were log transformed to reach normal distribution. The biomarker data were plotted separating ACR50 responders and nonresponders, and cutoffs were determined by areas under the receiver operating characteristic curves. Primary cutoffs were selected at sensitivity of a minimum of 70% (Sen70% bootstrapping). Secondary cutoffs were selected for MMP-3 and C1M at sensitivity of a minimum of 60% (Sen60% bootstrapping). The same was done for the ratio between C3M and CRPM; however, these were also plotted as a scatterplot (responder vs nonresponder, not shown), where distribution patterns were investigated to identify a threshold range including must responders (minimum 70%). The odds ratio (OR) for likelihood of being an ACR50 responder with a biomarker level at baseline above/below the set cutoff levels was determined by 2 × 2 tables. A decision tree was used to segregate ACR50 responders and nonresponders.

RESULTS

Determination of Cutoffs to Be Used in the Decision Tree for Predicting ACR50 Response

The MMP3 and C1M scatterplot showed that the variances were similar for the responder and nonresponder groups, whereas the variance of C3M/CRPM was markedly lower for responder group compared with the nonresponder group (Fig. 1). High level of MMP3 was significantly associated with ACR50 response with ORs of 3.1 for both Sen70% and Sen60% (P < 0.001), and cutoffs were set as 1.57 and 1.64 (Fig. 1A). There were trends toward a lower level of C1M in the responder group; however, these were not significant (Fig. 1B). Cutoffs were set at 2.04 and 2.00 for Sen70% and Sen60%, respectively (Fig. 1B). The scatterplot of the C3M/CRPM showed that 70% of responders felt in the range between 0.30 and 0.51 (Fig. 1C) with an OR for response of 2.4 (P = 0.019).

FIGURE 1

Determination of cutoffs to be used in the decision tree for predicting ACR50 response. The cutoffs were set for the biomarkers at sensitivity levels 70% and 60% (Sen70% and Sen60%) for MMP3 (A) and C1M (B) and at Sen70% for C3M/CRPM (C), shown by the dotted lines. Odds ratios for being a responder for the sensitivity levels are shown in the top right corner of each graph.

Segregation ACR50 (Week 52) Responders and Nonresponders

The cohort’s overall ACR50 response rate was 27% (Fig. 2). As serum MMP3 was the strongest predictor of ACR50 response by logistic regression, it was entered into the first level of CART followed by C3M/CRPM and C1M. In the first decision tree (Fig. 2A), 22% of the population was selected; response rate increased to 54%. In the tree, 86% of the nonresponders were identified as nonresponders, whereas 43% of the responders were identified as responders, giving an OR of 4.7 for prediction of response.

FIGURE 2

Segregation of ACR50 (week 52) responders and nonresponders by measurement of baseline biomarkers. A positive test was set for selection of those patients who were most likely to respond to TCZ and a negative test for identification of those patients who are most likely to be nonresponders. Dx pop indicates the percentage of patients with a positive test. OR, the OR of being a responder; true negative, the rate of nonresponders with a negative test.

In the second tree (Fig. 2B), the secondary cutoff value for MMP3 was used; 15% of the population was selected; response rate increased to 63%. In addition, 93% of the nonresponders had a negative test, whereas 35% of the responders had a positive test, giving an OR of 6.7 for prediction of response.

In the third tree (Fig. 2C), the secondary cutoff for C1M was used; 10% of the population was selected with a response rate of 70%. Ninety-six percent of the nonresponders were deselected, and 26% of the responders were selected, with an OR of 8.2 for prediction of response.

DISCUSSION

Biologic RA therapies provide on average 20% to 45% ACR50 response rates in phase III clinical studies, demonstrating that a significant number of patients derive insufficient benefit from therapeutic intervention. In this small cohort of a phase III clinical study, significantly improved response rates were achieved through analysis of 4 protein fingerprint biomarkers measured at baseline. In the first decision tree, 86% of the ACR50 nonresponders were positively identified by the negative test (low MMP, high/low C3M/CRPM, and high C1M at baseline). Thus, the biomarkers may provide means for deselection of patients who may not respond sufficiently. In addition, 22% of the patients had a positive test; patients with a positive test had 4.7 time chance of benefit from treatment with TCZ. It seemed that by measuring the level of biomarkers at baseline, it may enable selection of a treatment population where the response rate can be increased from 27% to 54%. These percentages could be further refined to increase response rates to 70%, albeit on the expense of selection of a smaller subpopulation.

The model was also tested for predictability for Disease Activity Score in 28 Joints (DAS28) remission rate at week 52. Of the 200 patients, only 112 had their DAS28 recorded at week 52, resulting in low power; thus, data were not shown. However, the model could indeed somewhat predict who experienced remission after treatment with 4 mg/kg TCZ with an OR of 2.5 (P = 0.057). This needs to be validated in a larger cohort.

The biomarkers measured in current cohort are protein fingerprints, which are direct measures of connective tissue degradation and inflammation that are downstream of the proinflammatory pathways, for example, the IL-6 or tumor necrosis factor α signaling pathways. As these are downstream biomarkers, it may be speculated that this class of biomarkers will show less fluctuation than cytokine markers or acute phase proteins, such as CRP. This hypothesis needs to be tested further. Current work presents an example of how biomarkers potentially can be used for personalized medicine, which should be validated in another TCZ trial before finally concluding its validity. The limitation of the cohort is that it is from a phase III clinical trial, thus not representing the general heterogeneous RA population.

This exploratory analysis was undertaken on a relatively small, but well-characterized, patient cohort and thus requires validation in an independent second cohort, which will also allow us to determine if this response profile is specific for anti–IL-6 receptor intervention or may be generalized to other populations. However, the perspective is that an intervention-predictive model would potentially be of value for both the industry and physicians and hopefully in the end for the patients. The RA field is in need of personalized medicine, and as the choice of treatment is becoming more complex and there is a demand for higher response rates, there is a home for easily assessable tools.

TABLE 1

Baseline Description of the Study Biomarker Study Population