There are numerous non-biologic and biologic disease-modifying anti-rheumatic drugs (bDMARDs) for rheumatoid arthritis (RA). Typical sequences of bDMARDs are not clear. Future treatment policies and trials should be informed by quantitative estimates of current treatment practice.
We used data from Corrona, a large real-world RA registry, to develop a method for quantifying sequential patterns in treatment with bDMARDs. As a proof of concept, we study patients who eventually use tocilizumab monotherapy (TCZm), an IL-6 antagonist with similar benefits used as monotherapy or in combination. Patients starting a bDMARD were included and were followed using a discrete-state Markov model, observing changes in treatments every 6 months and determining whether they used TCZm. A supervised machine learning algorithm was then employed to determine longitudinal patient factors associated with TCZm use.
7300 patients starting a bDMARD were followed for up to 5 years. Their median age was 58 years, 78% were female, median disease duration was 5 years, and 57% were seropositive. During follow-up, 287 (3.9%) reported use of TCZm with median time until use of 25.6 (11.5, 56.0) months. Eighty-two percent of TCZm use began within 3 years of starting any bDMARD. Ninety-three percent of TCZm users switched from TCZ combination, a TNF inhibitor, or another bDMARD. Very few patients are given TCZm as their first DMARD (0.6%). Variables associated with the use of TCZm included prior use of TCZ combination therapy, older age, longer disease duration, seronegative, higher disease activity, and no prior use of a TNF inhibitor.
Improved understanding of treatment sequences in RA may help personalize care. These methods may help optimize treatment decisions using large-scale real-world data.