Abstract
Introduction:
Achieving optimal adherence to zoledronic acid (ZOL) among osteoporosis (OP) patients is a challenging task. Here, we aimed to develop and validate a precise and efficient prediction tool for ZOL nonadherence risk in OP patients.
Methods:
We prospectively collected and analyzed survey data from a clinical registry. A total of 1010 OP patients treated for the first time with ZOL in two separate hospitals were selected for nonadherence analysis. The evaluation included a 16-item ZOL Nonadherence Questionnaire and potential risk factors for ZOL nonadherence were assessed via univariate and multivariate analyses. We next developed and validated two distinct-stage nomograms. Discrimination, calibration, and clinical usefulness of the predicting models were assessed using the area under the curve (AUC), calibration curves, and decision curve analysis (DCA).
Results:
The total nonadherence rate was 20.30% after the first ZOL infusion. To generate a model predicting ZOL nonadherence risk, six predictors of 16 items were retained. Model 2 (AUC, 0.8486; 95% confidence interval [CI], 0.8171–0.8801) exhibited considerably more discrimination in desirable functional outcomes, relative to Model 1 (AUC, 0.7644; 95% CI, 0.7265–0.8024). The calibration curves displayed good calibration. DCA revealed that a cutoff probability of 5–54% (Model 1) and 1–85% (Model 2) indicated that the models were clinically useful. External validation also exhibited good discrimination and calibration.
Conclusions:
This study developed and validated two novel, distinct-stage prediction nomograms that precisely estimate nonadherence risk among OP patients receiving the first infusion of ZOL. However, additional evaluation and external validation are necessary prior to widespread implementation.
Introduction
Osteoporosis (OP) is manifested by low bone mass and disintegrated bone microarchitecture that severely enhance fracture risk. 1 In mainland China, OP is diagnosed in 30–40% women and 10–20% men over 50 years of age. 2 Based on the large population in China, OP and its associated fractures signify a massive societal expense that rises each year with a growing aged population. 3 In addition, OP and related fracture management also produces a considerable economic burden that is on a steep rise all over the world. 4
Unfortunately, low adherence is far too common in authorized anti-OP therapies. In a prior report, only ~50% of patients complied with the anti-OP therapeutic regimen after 1 year and about a third followed the regimen after 2 years. 2 As expected, nonadherence to anti-OP therapy is strongly related to elevated fracture risk, morbidity, mortality, and expenses. 5 Fortunately, altering dosing interval of the anti-OP medication to weekly or monthly applications can improve patient compliance, but it still remains suboptimal.6–8 Zoledronic acid (ZOL), the most commonly used intravenous bisphosphonates for diagnosed osteoporosis, suppresses osteoclast-driven bone resorption and lowers vertebral fracture risk. 9 But, achieving optimal clinical efficacy requires the patient to return for the next infusion. However, a number of studies reported that 30–60% OP patients discontinued ZOL treatment by the end of the first year, following initiation.10–13
The reasons for ZOL nonadherence are multifactorial. It not only involves patients, but is also dependent on the specified therapy or healthcare system, which may influence a patient’s ability to adhere to the treatment regimen. 14 The World Health Organization (WHO) highlights the multifaceted causes of medication nonadherence, 15 which include socio-economy, patient, therapy, condition, and healthcare system. 15 Often times, nonadherence is the end result of multiple causes. Moreover, a patient’s perception and motivation for compliance depends on a self-perceived risk–benefit ratio. 15 Here, we generated a nonadherence prediction model, using readily available and relevant data, to aid clinicians in personalizing treatment strategies that encourage compliance.
Methods
Study patients
This is a retrospective evaluation of ZOL nonadherence, based on medical data from a clinical registry. This work is registered in the Chinese Clinical Trial Registry (ChiCTR2000036375). Among 1595 consecutive OP patients treated with ZOL in the Affiliated Kunshan Hospital of Jiangsu University (AKSHJU), between March 2018 and June 2020, we chose, for analysis, patients who met our strict inclusion criteria, the details of which are as follows: (1) the first infusion of ZOL (n = 1328) and (2) able to monitor and record data, as required by this study (n = 1221). In all, 211 patients were eliminated from analysis. The reasons were as follows: (1) allergy to any components of ZOL (n = 12), (2) serum total calcium level of < 2.13 mmol/L (8.5 mg/dL) (n = 30), (3) prior usage of other types of bisphosphonates (n = 91), (4) atrial fibrillation and other serious cardiovascular diseases (n = 34), and (5) loss to follow-up (n = 44). Ultimately, our analysis included 1010 patients. Among them, 756 patients from the main hospital of AKSHJU were assigned to the development set (DS), which included 157 nonadherent patients and 599 adherent patients. In addition, 254 patients from the You-yi branch hospital of AKSHJU were recruited as the validation set (VS), and it included 48 nonadherent and 206 adherent patients (Figure 1). We obtained written informed consent from all study subjects or their caregivers and our work received approval from the ethical committee at the Affiliated Kunshan Hospital of Jiangsu University (approval No. 2021-06-015-K01). Patient information was initially documented for hospital’s quality improvement purposes. Data analyzers were blinded to the identity of patients. Finally, we followed the guidelines of the Declaration of Helsinki and its recent amendments.
Flowchart.
OP management
The OP treatment regimen followed institutional protocols and were guided by doctors directly involved with patient cases. All patients received routine OP management services or OP liaison services (OLS), which entailed extensive assessment and care encompassing two primary services: medication management (MMS) and fracture liaison (FLS). Enrollment in the FLS program at AKSHJU was based on the following criteria: (1) age ⩾ 50 years; (2) newly identified fragility fracture; (3) local resident and expected to complete the 2-year outpatient follow-up evaluations. Patients who experienced fracture from major trauma, cancer, femoral shaft atypical fracture, or those who participated in a different fracture or OP therapy study were eliminated from this study. Patients were recruited for the MMS program if they met the following requirement: (1) aged ⩾ 50 years; (2) newly initiated OP medication or poor compliance to current OP medication; (3) local resident expected to complete the 2-year outpatient follow-up evaluations. The exclusion criteria for this program resembled that of the FLS program.
Our hospital specially developed an electronic system for the OLS program, which consists of three independent modules: coordinators, medical staffs, as well as patients and their families (Software Copyright ID Number: 2021SR0760693, 2021SR0802762, and 2021SR0791538, CN). The FLS program of the electronic management system was designed according to the 13 Best Practice Framework (BPF) standards established by the International Osteoporosis Foundation (IOF).16,17 These actions allowed coordinators to readily coordinate the electronic systems for standardized, personalized, and efficient OP management.
Survey structured questionnaire
A multidisciplinary team of experts developed this structured, self-administered questionnaire (16-item ZOL Nonadherence Questionnaire [ZOLNQ-16], see Supplementary File 1) related to the potential predictors of medication nonadherence. The potential predictors were stratified according to the WHO’s five categories of medication adherence. 15 The effectiveness of ZOLNQ-16 was examined. Our aim was to test this questionnaire under different conditions to assess its precision and feasibility. The ZOLNQ-16 questionnaire consists of 16 items encompassing socio-economic (4 items), patient-associated (5 items), therapy-associated (3 items), healthcare system-associated (2 items), and condition-associated (2 items) areas. ZOLNQ-16 can be divided into two stages, based on the patients’ follow-up time in relation to the initial ZOL treatment. The stage 1 questionnaire can be completed as early as 1 week after initial ZOL infusion, and it includes items 1 to 13. The stage 2 questionnaire can be completed between 1 month and 12 months after the initial ZOL infusion, and it includes items 1 to 16. The ZOLNQ-16 is designed to be easy to complete and attain optimal adherence and results. Hence, it consists of a small number of items that takes approximately 8 min to answer.
Definition of items
The detailed definitions of items in ZOLNQ-16 are summarized in Table 1. Musculoskeletal pain evaluation was based on a numeric rating scale (NRS; 0 = no pain, 10 = the worst possible pain) 18 before (NRS1, item 10) and after (NRS2) ZOL treatment. Since the baseline NRS scores were vastly different among patients, the difference value between the two NRS (NRSD, NRS2 − NRS1; item 16) was employed to assess the extent of pain relief. Acute-phase response (APR, item 11) was defined as fever, musculoskeletal pain, and other flu-like symptoms occurring within 72 h of ZOL infusion.19,20
Sixteen-item Zoledronic Acid Nonadherence Questionnaire (ZOLNQ-16).
APR, acute-phase response; BMI, body mass index; NSAID, nonsteroidal anti-inflammatory drug; NRS, numeric rating scale; OLS, osteoporosis liaison service; ZOL, zoledronic acid.
Data collection
All patients completed ZOLNQ-16 at 12 months after the initial ZOL infusion by themselves or with the help of others. Data were prospectively collected over the phone or electronically by highly skilled OP coordinators. At 14 months after initial ZOL infusion, all patients were contacted to confirm whether the next expected ZOL dose was administered or not. Nonadherence was defined as a delay of >56 days in receiving the next expected dose.21,22 All participating coordinators were specially trained to maintain standardization of study procedures to minimize systematic error.
Statistical analysis
Data analyses were done with the R Studio software (R version 3.6.1) and a p value < 0.05 was set as the significance threshold. The 16-item characteristics between adherence and nonadherence cohorts were analyzed via χ2 test or independent sample t test. Risk factor analyses employed both univariate and multivariate logistic regression analyses. Statistically significant variables from univariate analysis were entered into multivariate logistic regression analysis (p value < 0.1). Multicollinearity among the covariates was tested by examining whether the variance inflation factor (VIF) is ⩾4, and the multivariate logistic regression analysis was employed to build the models. We developed both early- (Model 1) and late-stage models (Model 2), depending on the stage of predictors included. Along with each model, we also expressed odds ratios with 95% confidence intervals (CI) to signify the impact of each risk factor.
Using independent variables regression coefficients from the DS, we established two nomogram prediction models associated with nonadherence.23,24 The prediction models were assessed for discrimination and calibration. The area under the curve (AUC) of the receiver operating characteristic (ROC) curve was computed to establish discrimination of both models of nonadherence prediction. In addition, we also established the optimal threshold for nonadherence prediction. Using ROC analysis, we also generated values representing optimal sensitivity and specificity for prediction of a given threshold. We established that an AUC of 0.5 ± 0.75 is acceptable, whereas AUC > 0.75 represents outstanding discrimination. 25 We next determined the integrated discrimination index (IDI) and net reclassification index (NRI) between the two models.26,27 Significant differences between ROC curves were further confirmed with the DeLong’s test for two associated ROC curves (R package pROC). 28 Calibration was assessed with calibration curves.
Furthermore, the decision curve analysis (DCA) was performed to assess the clinical feasibility of these models. This was done by quantifying the net benefits at various threshold probabilities in the DS. 29 The net benefit was calculated by subtracting the proportion of all patients who are false positive from the proportion of the patients who are true positive and by weighing the relative harm of forgoing interventions compared with the negative consequences of an unnecessary intervention. 30
Then, internal validation was completed using bootstrapping (BS) techniques (1000 BS resampling). Finally, external validation was performed in the VS using ROC and calibration curves using the above methods.
Sample size calculation
We calculated a minimum requirement of 88–107 participants to achieve a 99% power in our ROC investigation, carrying a two-sided α of 0.05, considering that the AUC is typically 0.75, and 30–60% OP patients discontinued ZOL treatment by the end of the first year, after initiation.10–13
Results
Sample characteristics
Please refer to Table 2 for all patient information, including the ZOL nonadherence risk-related 16 item characteristics from DS and VS. The ages of all participants were between 50 and 95 years, with an average of 68.12 years (SD 9.08). Among them, 82.87% (n = 837) were female. In addition, the total nonadherence rate was 20.30% after the first infusion of ZOL.
The 16-item characteristics of the development and validation sets.
NRSD = NRS2 − NRS1 (NRS1, NRS score before ZOL infusion; NRS2, NRS score after ZOL infusion).
APR, acute-phase response; BMI, body mass index; NRS, numeric rating scale; OLS, osteoporosis liaison service; ZOL, zoledronic acid.
Nomogram development
Based on univariate analysis, the ZOL nonadherence risk factors included insurance, APR, prophylaxis for APR, OLS, inconvenience, fracture, and NRSD. Variables reaching significance in univariate analysis were then entered into the non-conditional binary multivariate logistic regression. Subsequently, we revealed that insurance, APR, OLS, inconvenience, fracture, and NRSD were independent risk factors for medication nonadherence (Table 3). We next conducted collinearity analysis on the aforementioned independent risk factors (VIFs = 1) and showed no multicollinearity among these factors.
Univariate and multivariate logistic regression models of ZOL nonadherence in the development set.
NRSD = NRS2 − NRS1 (NRS1, NRS score before ZOL infusion; NRS2, NRS score after ZOL infusion).
APR, acute-phase response; BMI, body mass index; CI, confidence interval; NRS, numeric rating scale; OLS, osteoporosis liaison service; OR, odds ratio; ZOL, zoledronic acid.
Based on the logistic multivariate regression analysis, the six independent risk factors were finally included (Table 4). We next constructed an early-stage nomogram prediction model (Model 1) and a late-stage nomogram prediction model (Model 2) for the prediction of patient ZOL nonadherence risk. Schemes of two developed nomograms are shown in Figure 2. Model 1 can be used as early as 1 week after initial ZOL infusion and it includes early-stage predictors (insurance, APR, and OLS). In contrast, Model 2 can be used between 1 month and 12 months after the initial ZOL infusion, and it was constructed with Model 1 adding late-stage predictors (inconvenience, fracture, and NRSD). Both nomograms are presented in Figure 3. The nomogram can be used to retrieve values for each predictor, which can be summed up to receive a total score. The estimated risk corresponding to the total score reveals the likelihood of ZOL nonadherence.
Multivariable logistic regression models in the development set.
NRSD = NRS2 − NRS1 (NRS1, NRS score before ZOL infusion; NRS2, NRS score after ZOL infusion).
APR, acute-phase response; CI, confidence interval; NRS, numeric rating scale; OLS, osteoporosis liaison service; OR, odds ratio; ZOL, zoledronic acid.
Schemes of two developed nomograms for the prediction of ZOL nonadherence risk. Model 1 (black) can be used as early as 1 week after initial ZOL infusion and it includes three predictors. Model 2 (red) can be used from 1 month to 12 months after the initial ZOL infusion. Its construction is similar to Model 1, with the addition of three late-stage predictors. NRSD = NRS2 − NRS1 (NRS1, NRS score before ZOL infusion; NRS2, NRS score after ZOL infusion; NRS, numeric rating scale).
Two distinct-stage individualized nomograms for the prediction of ZOL nonadherence risk in osteoporosis patients. Model 1 (early-stage) includes insurance, APR, and OLS. Model 2 (late-stage) includes insurance, APR, OLS, inconvenience, fracture, and NRSD. The nomogram is used as follows: First, each predictor is located on the corresponding axis, then a line is drawn from the predictor to the points axis to obtain the point value; subsequently, all points from all predictors are added together and another line is drawn from the total point axis to the lower line of the nomogram to establish ZOL nonadherence risk. NRSD = NRS2 − NRS1 (NRS1, NRS score before ZOL infusion; NRS2, NRS score after ZOL infusion; NRS, numeric rating scale).
For example, in Model 1, the scores were provided as follows: an OP patient who received OLS (0 point), had no medical insurance for ZOL purchase (100 points), and experienced APR after ZOL infusion (45 points). The cumulative score of the various predictive indicators was 0 + 100 + 45 = 145, and the corresponding predicted risk of nonadherence was 0.37 (37%) (Figure 4).
An example prediction nomogram for the estimation of ZOL nonadherence risk in osteoporosis patients.
Model performance
The Model 1 threshold was 0.1865, carrying a sensitivity of 77.07% and specificity of 66.78%. The Model 2 threshold was 0.1672, carrying a sensitivity of 84.08% and specificity of 69.45%. According to the AUC differences between the two models, Model 2 (AUC, 0.8486; 95% CI, 0.8171–0.8801) displayed remarkably high discrimination of desirable functional outcome, compared with Model 1 (AUC, 0.7644; 95% CI, 0.7265–0.8024; Table 5; Figure 5(a)). In addition, to estimate the discrimination performance by late-stage predictors more sensitively, we compared the probabilities of events and nonevents of models using the relative IDI and category-free NRI. Model 2 showed significant improvement of both IDI and NRI as compared with Model 1.
Comparisons of the two models by AUC, IDI, and NRI in the development set.
AUC, area under the curve; CI, confidence interval; IDI, integrated discrimination index; NRI, net reclassification index.
Evaluation of two models for the prediction of ZOL nonadherence risk in the development set. The ROC curve AUC of (a) calculated to assess the discriminative nature of both models. The AUC of Model 1 (black) and Model 2 (red) were 0.7644 (95% CI, 0.7265–0.8024) and 0.8486 (95% CI, 0.8171–0.8801), respectively. There was statistical significance between two models (p < 0.0001). The calibration curves assessed calibrations of Model 1 (b) and Model 2 (c).
Based on the internal bootstrap validation method, at a probability of 0–0.45, the Model 1 bias-corrected curve was comparable to the ideal curve. But, at a high probability (>0.45), Model 1 may overestimate the probability of ZOL nonadherence [Figure 5(b)]. Upon addition of late-stage predictors, the Model 2 bias-corrected curve was close to the ideal curve, and good agreement was observed between predictions and observations (Figure 5(c)), demonstrating a good fit.
DCA
The DCA comparing clinical feasibility of Models 1 and 2 are presented in Figure 6. The probability cutoff for ZOL nonadherence was placed in the x-axis, whereas the standard net benefit of employing the model was placed in the y-axis. Based on the DCA, threshold probabilities of 5–54% for Model 1 and 1–85% for Model 2 were more beneficial than the intervention-all or intervention-none schemes. Besides, Models 2 was better than 1 within the limits of this range.
Decision curve analysis for the nonadherence models. The y-axis represents standardized net benefit. The blue (Model 1) and red (Model 2) line denotes ZOL nonadherence risk models. The thin oblique solid line assumes that all patients are noncompliant to ZOL. The thick horizontal solid line assumes that all patients are compliant to ZOL. Based on the decision curve, a cutoff of 5–54% for Model 1 and 1–85% for Model 2 are far more beneficial to patients than the intervention-all or intervention-none schemes.
Model validation
We performed BS-corrected ROC analysis on both models, and the results are presented in Supplementary File 2. High AUC value was still reached in the interval validation. Furthermore, the two models were next externally validated using ROC and calibration curves in VS (Figure 7). The Model 1 AUC for VS was 0.7714 (95% CI 0.6954–0.8473), carrying a sensitivity of 43.75% and a specificity of 97.09% at the optimal threshold. Alternately, the Model 2 AUC for VS was 0.8480 (95% CI: 0.7807–0.9154), carrying a sensitivity of 64.58% and a specificity of 93.20% at the optimal threshold. Finally, the calibration curves of both nomograms exhibited satisfactory agreement between the estimated and observed values.
External validation for ZOL nonadherence nomograms using ROC curves (a) and calibration curves (b and c) corresponding to the validation set.
APR risk findings
To further predict APR risk, all potential predictors in stage 1 (item 1 to 13 except 12) were employed to establish a prediction model in the DS by using the above methods. Unfortunately, the performance was poor (AUC, 0.6281; 95% CI, 0.5883–0.6679). Table 6 summarizes the results of univariate and multivariate analyses of APR risk. Based on multivariate analysis of APR risk, age (>79 years old [OR, 0.32; 95% CI, 0.18–0.59; p = 0.0003]), prophylaxis for APR (OR, 0.64; 95% CI, 0.47–0.87; p = 0.0042), and OLS (OR, 0.0212; 95% CI, 0.51–0.05; p = 0.0212) offer independent protection against APR.
Univariate and multivariate analyses of the APR risk in the development set.
NRS1 indicates NRS score before ZOL infusion.
APR, acute-phase response; BMI, body mass index; CI, confidence interval; NRS, numeric rating scale; OLS, osteoporosis liaison service; OR, odds ratio; ZOL, zoledronic acid.
Discussion
Nomograms are increasing popular as prognostic tools in oncology and medicine. They utilize user-friendly digital interfaces, enhanced precision, and comprehensible prognoses to build enhanced treatment plans. 31 This study offers a novel look into the application of nomograms in OP and ZOL nonadherence. We established and validated two new, stage-specific prediction models for nonadherence risk in OP patients receiving the first infusion of ZOL, using merely three or six easily available variables. External validation confirmed satisfactory discrimination and calibration power. In particular, the high Model 2 AUC in VS verified that the nomograms can be extensively and precisely used.
Based on our analysis, about 20% of OP patients failed to receive a second ZOL infusion. This result is consistent with the Curtis et al. 10 study, where it was reported that nonadherence reached 18% among U.S. Medicare enrollees. In other studies, nonadherence was relatively high, compared with our study. A retrospective observational study, using data from Humana, a US administrative healthcare claims database, revealed that among the 3128 OP patients receiving the first ZOL infusion, nonadherence was 59.2% at 12 months and 79.8% at 24 months. 11 A study of a German database revealed that the proportion of patients that remained on ZOL treatment at 12 months was 65.6% (nonadherence = 34.5%). 13 In a study examining 259 Korean patients, only 36% (nonadherence = 64%) received a second ZOL infusion. 12 Although, compared with other treatment option, an annual intravenous administration of 5 mg ZOL is a highly efficient and favorable, 32 it is not uncommon to discontinue ZOL administration from the second year, based on the above data.
A systematic review identified 24 factors with 139 sub-factors that influence medication compliance in OP patients, based on the WHO’s five-domain framework. 33 Among the condition-associated factors regulating medication noncompliance was polypharmacy, and enhanced medication compliance usually resulted from a history of falls. Patient-associated factors resulting in poor medication compliance was old age and OP-related misconception. Therapy-associated factors promoting medication noncompliance were frequent dosing and medication side effects. Healthcare system–associated factors driving medication noncompliance were care under different medical specialties and lack of patient education. Finally, among the socio-economic–associated factors influencing medication noncompliance was smoking history and absence of medical insurance coverage. 33
Similar to the above results, the predictors of ZOL nonadherence identified in our study were insurance, APR, OLS, inconvenience, fracture, and NRSD. The nomograms suggested that lack of medical insurance coverage, APR occurrence, lack of OLS management, perception of medication regime being complicated, new fractures, and no pain relief after treatment are the main factors regulating ZOL nonadherence risk in OP patients.
Lack of medical insurance was the major socio-economic predictor of ZOL nonadherence. This concurs with an additional four studies.34–37 It is essential to assess patient financial situation before the prescription of expensive anti-OP medications. In case of patients in need of financial aid, they must be referred to appropriate service providers, such as hospital social services for financial assistance. In addition, regional governments must cover cost-effective anti-OP medications in patient medical insurance.
The occurrence of APR was found, unsurprisingly, as a predictor of ZOL nonadherence. This is in accordance with other publications38–41 that indicated that patients often forgo medications for fear of undesirable side effects. Although an adult may be familiar with the advantages of long-term medical compliance in managing OP, its benefits are not immediately apparent. Therefore, presence of side effects like APR dramatically reduces their motivation to follow through with the treatment regimen. So, regular evaluation and management of APR, with the help of a prediction model, are essential. Despite our analysis demonstrating poor performance of APR risk prediction, several independent protective factors associated with APR were identified. Older age (>79 years old) is among the first independent protective factors of APR, which is in accordance with the Ding et al. 20 study. Our prior assumption was that the relatively low proinflammatory cytokines present in the older population lead to the difference in APR. 42 The second independent protective factor for APR was nonsteroidal anti-inflammatory drug (NSAID) usage immediately after ZOL infusion to alleviate pain and fever associated with prophylactics. This is consistent with a clinical study. 43 OLS management was the third independent protective factor of APR. Once patients were enrolled into our hospital OLS, they were advised to drink more water. Adequate hydration may play an important role in reducing APR occurrence. Hence, clinicians can take special steps to mitigate APR prior to therapy or design alternative routes of drug administration under an elevated APR risk. In future, laboratory indicators can also be introduced as predictors, such as serum 25(OH)D levels, to better develop the APR risk prediction model.
Among healthcare team and system-related factors, lack of OLS management was found to be a predictor of ZOL nonadherence. Previous studies demonstrated that early interventions like medication reminders, drug monitoring, and family support can benefit patients with elevated noncompliance risk.44,45 When it comes to adhering to an infrequent dose intravenous therapy, the most important factor is to ensure the patient is scheduled and shows up to the appointment for an infusion treatment. 46 Therefore, employing low-cost reminding tools can augment compliance in high-risk patients. 45 Therefore, our healthcare team developed an electronic management system for the OLS program. Prior to the commencement of the first treatment, the electronic OLS system automatically generates a paper liaison leaflet containing the patient’s fracture/BMD/lab data and long-term treatment plan. In addition, coordinators, clinicians, OP patients, and their families can easily log into the electronic system using mobile phones and computers. This electronic OLS system provides intervention, as well as education, exercise training, medication adjustment, and referral to physiatrists. Meanwhile, automated medication reminders can schedule the following infusion appointment and send reminders to the patient to ensure adherence. If a patient fails to respond, the system can remind the coordinator to call by phone. The OLS system also provides a green reservation channel for patients to come to the hospital for regular follow-ups. The final results of our study demonstrated that a dedicated healthcare team and properly designed system are key to enhancing medication compliance in OP patients.
Similar to other studies,47–49 our analysis also confirmed that a patient’s false perception of a complicated treatment regimen can give rise to medication noncompliance. This may be due to the fact that during the 3–5 years of ZOL treatment, patients need to come to the hospital every year for an intravenous ZOL infusion. In addition, supplements like calcium and vitamin D must be taken throughout the treatment period. Polypharmacy is a well-recognized issue in the older population as they are often afflicted with multiple chronic diseases and are under numerous medications at the same time. 15 The need to remember and consume large number of medications on a daily basis, over a prolonged period of time, can be overwhelming to the elderly, thus giving the perception of a complicated treatment regimen. This can negatively affect medication compliance. Though ZOL has longer dosing intervals, it is still necessary to strengthen patient education and provide necessary support.
Experiencing fracture and no pain relief after ZOL infusion are two of the main therapy-related indicators of ZOL nonadherence. It may be attributed to a patient’s lack of understanding of therapeutic benefit, which is a well-established factor associated with poor medication adherence.50–52 Typically, nonadherence increases when treatment does not appear as effective as patients’ expectations. In this case, healthcare providers need to constantly guide patients to the realization that OP is a chronic disease and anti-OP treatment is a long and arduous task. In addition, patients need to be aware that a single fracture while on therapy does not indicate treatment failure. 53 However, it is also necessary for clinicians to consider alternative therapy or reassessment for causes of secondary OP in patients experiencing two or more refractures or significant bone loss during therapy. 53
As expected, OP patients who adhere to the treatment plan experience better outcomes than patients who were nonadherent.54,55 Therefore, establishing patient nonadherence prediction models can guide personalized risk assessment and interventions. Thus, we established a precise ZOL noncompliance risk prediction tool that can aid clinicians in evaluating the ZOL noncompliance risk, design personalized treatment plans, prevent additional interventions in low-risk patients, and avoid delays or discontinuity in intervention. Supporting interventions like patient education, monitoring and supervision, alteration in medication regimen, providing patient support, and interdisciplinary collaborations, 56 or an OLS system (as was done in our hospital) can highly benefit patients at elevated noncompliance risk who are starting their ZOL infusions. Moreover, this can also aid researchers in selecting optimal OP patients (exhibiting good adherence to anti-OP therapy) for clinical research. Likewise, models can also be used to exclude patients with poor compliance in retrospective studies, thereby securing more reliable analysis.
Our work encountered certain limitations. First, the data were mostly self-reported, which introduces certain social desirability and recall bias, and the values may not be accurate. 57 However, inaccuracies in reported nonadherence are believed to be minimal since independent, experienced, and trained OP coordinators, not patients’ healthcare providers, were responsible for the follow-up evaluations. Furthermore, in order to minimize potential recall bias, two distinct-stage prediction models were generated, depending on the phase of treatment the patients were in. In addition, from the patient’s perspective, when reporting ZOL nonadherence information, it was believed that ZOL was easily identified during the telephone or network follow-ups, due to its unique characteristics (intravenous administration, once a year). Second, the patients or clinicians were not blinded to the intervention, which may contribute to both selection and observational biases. Third, the follow-up duration was relatively short (14 months). This study did not include assessments of long-term ZOL nonadherence risk. Fourth, risk factor analysis of the current study did not include all potential factors that can regulate ZOL nonadherence. Possible predictors, such as living status and education level, need to be identified in the future. Finally, despite extensive examination of our models with VS, the feasibility of employing our nomograms on OP populations from other regions and countries remains to be examined.
Conclusion
The total nonadherence rate was 20.30% after the first infusion of ZOL. This study developed and validated two novel, distinct-stage prediction nomograms, with satisfactory precision to aid clinicians in accessing nonadherence risk among OP patients receiving the first infusion of ZOL. Lack of medical insurance coverage, APR occurrence, lack of OLS management, perception of medication regime being complicated, new fractures, and no pain relief after treatment are among the main factors determining ZOL noncompliance. Corresponding interventions, such as OLS management, can enhance ZOL compliance in high-risk patients. Additional prospective validation, involving a large patient population, is warranted to confirm our results.
Supplemental Material
sj-docx-1-taj-10.1177_20406223221114214 – Supplemental material for Predicting the individualized risk of nonadherence to zoledronic acid among osteoporosis patients receiving the first infusion of zoledronic acid: development and validation of new predictive nomograms
Supplemental material, sj-docx-1-taj-10.1177_20406223221114214 for Predicting the individualized risk of nonadherence to zoledronic acid among osteoporosis patients receiving the first infusion of zoledronic acid: development and validation of new predictive nomograms by Chong Li, Ke Lu, Qin Shi and Ya-qin Gong in Therapeutic Advances in Chronic Disease
Supplemental Material
sj-docx-2-taj-10.1177_20406223221114214 – Supplemental material for Predicting the individualized risk of nonadherence to zoledronic acid among osteoporosis patients receiving the first infusion of zoledronic acid: development and validation of new predictive nomograms
Supplemental material, sj-docx-2-taj-10.1177_20406223221114214 for Predicting the individualized risk of nonadherence to zoledronic acid among osteoporosis patients receiving the first infusion of zoledronic acid: development and validation of new predictive nomograms by Chong Li, Ke Lu, Qin Shi and Ya-qin Gong in Therapeutic Advances in Chronic Disease
Footnotes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
