Skip to main content

Abstract

Objective:

Pimavanserin, a serotonin 5-HT2 antagonist, is indicated for treatment of hallucinations and delusions associated with Parkinson’s disease psychosis. In premarketing trials in patients with Parkinson’s disease psychosis, 11% of patients died during open-label pimavanserin treatment. Antipsychotics, which are used off-label in Parkinson’s disease psychosis, increase mortality in dementia patients. The authors compared mortality with pimavanserin and atypical antipsychotics in a large database.

Methods:

This was a retrospective new-user cohort study of Medicare beneficiaries with Parkinson’s disease initiating pimavanserin (N=3,227) or atypical antipsychotics (N=18,442) from April 2016 to March 2019. All-cause mortality hazard ratios and 95% confidence intervals were estimated for pimavanserin compared with atypical antipsychotics, using segmented proportional hazards regression over 1–180 and 181+ days of treatment. Potential confounding was addressed through inverse probability of treatment weighting (IPTW).

Results:

Pimavanserin users had a mean age of approximately 78 years, and 45% were female. Before IPTW, some comorbidities were more prevalent in atypical antipsychotic users; after IPTW, comorbidities were well balanced between groups. In the first 180 days of treatment, mortality was approximately 35% lower with pimavanserin than with atypical antipsychotics (hazard ratio=0.65, 95% CI=0.53, 0.79), with approximately one excess death per 30 atypical antipsychotic–treated patients; however, during treatment beyond 180 days, there was no additional mortality advantage with pimavanserin (hazard ratio=1.05, 95% CI=0.82, 1.33). Pimavanserin showed no mortality advantage in nursing home patients.

Conclusions:

Pimavanserin use was associated with lower mortality than atypical antipsychotic use during the first 180 days of treatment, but only in community-dwelling patients, not nursing home residents.
Parkinson’s disease is the second most common neurodegenerative disorder after Alzheimer’s disease (1), and its prevalence is expected to increase to more than 1 million patients in North America by 2030 (2). Parkinson’s disease is commonly associated with psychosis (3), and in patients with Parkinson’s disease, lifetime prevalence of visual hallucinations may approach 50% (4).
Pimavanserin, a serotonin 5-HT2 receptor antagonist, is the only drug approved in the United States for treatment of hallucinations and delusions associated with Parkinson’s disease psychosis. In premarketing trials of pimavanserin in patients with Parkinson’s disease psychosis, serious adverse events were more frequent during treatment with pimavanserin than with placebo (7.9% taking pimavanserin, compared with 3.5% taking the placebo in the 6-week trials), and during open-label pimavanserin treatment, 11% of patients died (5). Postmarketing pharmacovigilance data on pimavanserin includes multiple reports of deaths (6), but the data do not point to any specific, unexpected types of serious adverse events (7).
A meta-analysis of placebo-controlled trials in patients with dementia-related psychosis that analyzed data on 3,353 patients treated with atypical antipsychotics and 1,757 treated with placebo found that the risk of mortality was roughly 50% higher with atypical antipsychotics (118 deaths [3.5%] compared with 40 [2.3%]) (8). With respect to the Parkinson’s disease population, Parkinson’s disease psychosis has been treated off-label with atypical antipsychotics, which have dual dopamine and serotonin antagonism (9). A prospective cohort study of 230 patients with Parkinson’s disease found that presence of psychotic symptoms was associated with roughly a 45% increase in mortality, although antipsychotic use itself was not a significant independent predictor of mortality (10). However, a large retrospective cohort study of Veterans Health Administration patients with Parkinson’s disease (most without dementia) found higher mortality associated with antipsychotic use; as the authors noted, the comparison they made to antipsychotic nonusers may have been biased if the need for treatment with an antipsychotic was a risk factor for mortality (a situation known as confounding by indication) (11). More recently, a retrospective cohort study of Medicare beneficiaries with Parkinson’s disease who were residing in nursing homes found that pimavanserin users had a 24% higher 30-day risk of hospitalization than nonusers. Also, mortality was higher for long-term pimavanserin users than nonusers; the adjusted hazard ratio for death at 1 year was 1.56 (95% CI=1.42, 1.72). However, as with the Veterans Health Administration study, the authors noted that the comparison between users and nonusers could be confounded by indication (12). A single-center observational study of pimavanserin (N=113) and quetiapine (N=505) in Parkinson’s disease found, relative to untreated control subjects (N=784), an odds ratio for death of 1.74 (95% CI=1.15, 2.62) with quetiapine and 1.23 (95% CI=0.57, 2.68) with pimavanserin, adjusted for age and gender (13). Another single-center retrospective cohort study included 92 patients with either Parkinson’s disease or Lewy body dementia who were treated with quetiapine or pimavanserin; death during follow-up occurred in seven (15%) of the quetiapine-treated patients and four (9%) of the pimavanserin-treated patients (p=0.88) (14). Of note, mortality among patients with Parkinson’s disease using pimavanserin or atypical antipsychotics has not been compared in a large, nationally representative study.
To address knowledge gaps regarding mortality risks in patients with Parkinson’s disease treated with pimavanserin, we undertook a retrospective cohort study using U.S. Medicare data to evaluate all-cause mortality in patients with Parkinson’s disease using pimavanserin and atypical antipsychotics.

Methods

Data Source

Medicare is the U.S. government’s health insurance program for individuals age 65 years and older and for persons with end-stage renal disease or disability. We studied Medicare beneficiaries enrolled in fee-for-service Medicare Part A (hospitalization), Part B (outpatient medical care), and Part D (prescription drug coverage). We used anonymized Medicare administrative claims records to obtain data on inpatient or outpatient diagnoses and procedures, outpatient drug prescriptions, and deaths. As this retrospective analysis was considered a public health surveillance activity, ethics committee review and signed patient consent were not required by the U.S. Food and Drug Administration.

Study Population

We employed a new user retrospective cohort design to compare patients initiating treatment with pimavanserin or oral atypical antipsychotics. Eligible atypical antipsychotics were those available in the United States aside from clozapine, as it is available only through a special program because of the risk of severe neutropenia. From April 29, 2016, when pimavanserin was approved, through March 2019, patients were eligible for study inclusion if, on the date of qualifying prescription fill (index date), they had at least 1 year of continuous Medicare Parts A, B, and D enrollment, were age 65 or older, and during the year before their index date had at least one Parkinson’s disease diagnosis based on ICD-9 and ICD-10 coding (ICD-9 code 332, ICD-10 code G20), had at least one levodopa prescription, and had no prescriptions for pimavanserin or any antipsychotic. Patients were excluded if they received prescriptions for multiple atypical antipsychotics or received both pimavanserin and antipsychotics on the index date. Patients were also excluded if they were hospitalized on the index date or received hospice care in the prior year. Additionally, patients with a stay in a skilled nursing facility in the 30 days prior to the index date were excluded, as their prescriptions may not be reliably captured while in a skilled nursing facility. Patients with a chronic dialysis claim during the baseline period or who entered Medicare with end-stage renal disease were also excluded, as these patients may have greater mortality risk. We excluded patients with diagnoses of schizophrenia, schizoaffective disorder, or schizophreniform disorder during the baseline period, as these diagnoses may indicate a potential alternative indication for chronic atypical antipsychotic use.

Covariates

For the 1-year baseline period, Medicare claims data on chronic medical conditions, health care utilization, and dispensed medications were collected. We also used these data to produce values for the Charlson comorbidity index (15) and the frailty index (16, 17).

Follow-Up

Follow-up began the day after cohort entry and continued until whichever of the following came first: disenrollment from Medicare Parts A, B, or D; stopping treatment, defined by a gap in days’ supply between successive study drug prescriptions exceeding 14 days; dispensing of a nonstudy antipsychotic; switching from pimavanserin to an atypical antipsychotic or vice versa; end of the study period; or death.

Outcome

The primary outcome was death from any cause. Previously, we found high concordance between deaths identified in Medicare and the National Death Index (kappa=0.98) (18). As a sensitivity analysis, we analyzed a composite outcome of death or admission to hospice.

Statistical Analysis

If users of pimavanserin differed in important characteristics from users of atypical antipsychotics, bias (confounding) could occur when comparing mortality risk. To address potential confounding, we used inverse probability of treatment weighting (IPTW) to estimate the average treatment effect in pimavanserin-treated patients (see the Supplementary Statistical Methods section in the online supplement). We assessed balance on baseline characteristics (see Table S1 in the online supplement) using standardized mean differences calculated before and after weighting, considering a standardized mean difference of ≤0.10 as a negligible imbalance in the baseline covariate (19). Weighted Kaplan-Meier plots were used to characterize risk of death over time (20). Weighted Cox proportional hazards regression with robust variance estimation was used to estimate the hazard ratio and 95% confidence interval for comparison of duration of survival among the atypical antipsychotic and pimavanserin groups. Our initial comparison revealed that the hazard ratio was not constant over the duration of follow-up, indicating that the proportional hazard assumption was violated. Consequently, we performed segmented proportional hazards modeling with time intervals determined primarily from inspection of Schoenfeld residual plots (21). In addition, we calculated the cumulative number needed to harm for one additional death, based on time-to-event data (22). To assess the potential for bias in the primary analysis, we calculated the E-value, which is the minimum strength of residual confounding required to fully explain away the observed treatment effect (23).
In a secondary analysis, we limited our comparison to patients who received their index prescription for pimavanserin or atypical antipsychotics from neurologists, to enhance comparability and further reduce the possibility that atypical antipsychotic treatment was unrelated to Parkinson’s disease psychosis. Additionally, to assess whether the hazard ratio for mortality was different between beneficiaries who had been admitted to hospice and beneficiaries who had not, we estimated different hazard ratios during the interval before a beneficiary entered hospice and during the interval after a beneficiary was observed to have a hospice admission.
We conducted a priori subgroup analyses by age, gender, comorbid dementia (defined by a diagnosis or use of an antidementia drug during the year prior to index date), comorbid depression (defined by a diagnosis in the year prior to index date), nursing home residency on index date, Charlson comorbidity index, and frailty score. Post hoc, we analyzed dementia patients identified with antidementia drug use alone.
A priori sensitivity analyses included increasing the gap between prescription refills from 14 to 30 days; restricting the sample to patients with exposure defined by at least two prescriptions; restricting the sample to patients receiving exactly 34 mg/day of pimavanserin, the recommended daily dose; restricting the sample to patients receiving care from neurologists in the 90 days prior to the index date, to account for possible differences in management of Parkinson’s disease psychosis for these patients; and censoring for entry to a skilled nursing facility, to ensure capture of dispensed prescriptions. We also performed an analysis restricted to patients who received quetiapine as their initial atypical antipsychotic, and in a post hoc analysis, we compared pimavanserin users to users of non-quetiapine atypical antipsychotics. To assess whether there was a difference in risk of clinical deterioration preceding mortality, we also assessed a composite outcome of death or admission to hospice. Because pimavanserin prolongs the QT interval, it should be avoided by patients with congenital QT prolongation or other risk factors for QT prolongation, and it should not be combined with other QT-prolonging drugs (24). Accordingly, we conducted an analysis accounting for concomitant use of potent QT-prolonging drugs (25) (see Table S2 in the online supplement) as a time-varying covariate. Finally, we repeated the analysis without any adjustment and without weighting, using multivariable Cox regression with all covariates, and covariates selected with LASSO regression (26). We performed all secondary, subgroup, and sensitivity analyses using segmented Cox regression.
We used R, version 3.6.0 (R Foundation for Statistical Computing, Vienna) and SAS, version 9.4 (SAS Institute, Inc., Cary, N.C.) for statistical analysis.

Results

We identified 3,227 eligible pimavanserin-treated patients and 18,442 eligible atypical antipsychotic–treated patients (see Table S3 in the online supplement). The median exposed follow-up time was 78 days (interquartile range=44–202) for pimavanserin and 74 days (interquartile range=44–170) for atypical antipsychotics. The reasons for censoring during follow-up were similar by cohort (see Table S4 in the online supplement). Among atypical antipsychotic users, the most frequent initial atypical antipsychotic prescription was for quetiapine (in 78% of atypical antipsychotic users), followed by risperidone (9%), olanzapine (6%), aripiprazole (5%), and other atypical antipsychotics (1%). The mean age of pimavanserin-treated patients was 78 years, and 45% were female. Prior to weighting, compared to pimavanserin users, atypical antipsychotic users were more likely to visit a primary care provider and to have had a hospital admission in the 30 days before study entry. Certain comorbidities, including atrial fibrillation, chronic obstructive pulmonary disease, and diabetes, were more common in atypical antipsychotic users, and these patients’ mean Charlson comorbidity index and frailty scores were higher. Additionally, pimavanserin-treated patients tended to use higher doses and more classes of anti-Parkinson medications. After weighting, the cohorts were closely balanced on all covariates (Table 1; see also Table S1 in the online supplement).
TABLE 1. Baseline characteristics of the sample of patients with Parkinson’s disease treated with pimavanserin or atypical antipsychoticsa
CharacteristicPimavanserinAtypical Antipsychotics
Before WeightingAfter Weightingb
 NN N 
Eligible population3,22718,442 3,251 
 MeanMeanStd. mean diff.cMeanStd. mean diff.c
Age (years)77.978.90.1577.80.01
 %%Std. mean diff.c%Std. mean diff.c
Male55.455.30.0055.70.01
Race     
 White88.789.10.0188.90.01
 Black4.33.80.034.30.00
 Other7.07.10.006.80.01
Prescriber specialty     
 Neurology71.538.30.7171.90.01
 Psychiatry3.08.20.232.90.01
 Primary care24.651.00.5724.30.01
 Other0.92.50.120.90.00
In nursing home on index date15.213.40.0515.10.00
Medical comorbidities     
 Neuropsychiatric     
  Depressive disorder42.148.60.1341.60.01
  Anxiety32.940.40.1532.80.00
  Dementia56.560.20.0756.30.00
 Cardiovascular     
  Atrial fibrillation16.020.70.1215.90.00
  Congestive heart failure14.318.80.1214.40.00
  Other ischemic heart disease28.634.20.1228.70.00
  Peripheral vascular disease32.534.90.0532.40.00
  Hospitalized stroke0.71.50.070.70.00
  Hypertension76.180.10.1076.00.00
 Other cerebrovascular disease19.524.70.1219.30.00
 General medical conditions     
  Anemia35.237.90.0635.20.00
  Chronic obstructive pulmonary disease12.017.30.1512.00.00
  Diabetes25.631.80.1425.40.01
  Hypercholesterolemia36.038.10.0435.90.00
  Chronic kidney disease17.020.40.0916.80.00
  Malignancy16.017.10.0316.00.00
  Nicotine dependency21.926.70.1122.00.00
QT-prolonging drug prescription covering the index date33.738.10.0933.50.01
Parkinson’s disease treatment medications     
 Number of days covered by levodopa, past 0–365 days     
  1–180 days10.820.50.2710.70.00
  181–365 days89.279.50.2789.30.00
 Dopamine agonists35.025.50.2135.40.01
 Monoamine oxidase B inhibitors20.711.80.2421.40.02
 Amantadine12.99.60.1013.20.01
Medications for treatment of cardiovascular conditions     
 ACEI/ARBsd33.638.10.0933.50.00
 Beta-blockers31.537.80.1331.30.00
 Calcium channel blockers20.924.00.0820.70.00
 Diuretics (loop)18.320.90.0718.20.00
 Drugs for hypotension13.48.70.1513.60.00
 Statins45.349.70.0945.10.01
Medications for treatment of neuropsychiatric disorders     
 Antidementia drugs45.442.20.0645.40.00
 Antianxiety drugs22.631.70.2122.50.00
 Antidepressants59.467.00.1659.20.00
 MeanMeanStd. mean diff.cMeanStd. mean diff.c
Charlson comorbidity index2.803.300.162.810.00
Frailty score (past 0–183 days)0.310.350.210.310.01
a
The denominator for percentages is the number of eligible patients in the pimavanserin and atypical antipsychotic groups.
b
Denotes the weighted pseudo-population created by applying inverse probability of treatment weighting–average treatment effect on the treated (IPTW-ATT) weights for the cohort treated with atypical antipsychotic.
c
Standardized mean difference between the pimavanserin cohort and the atypical antipsychotic cohort.
d
ACEI/ARBs=Angiotensin-converting enzyme inhibitors or angiotensin receptor blockers.
During follow-up, 207 pimavanserin users and 1,752 atypical antipsychotic users died. Overall, pimavanserin users exhibited a higher cumulative survival probability through 360 days of follow-up (Figure 1). The risk of death was lower among pimavanserin users (hazard ratio=0.77, 95% CI=0.66, 0.90). However, statistical evaluation showed that the proportional hazard assumption was violated. Guided by inspection of the Schoenfeld residuals (see Figure S1 in the online supplement), we analyzed days 1–180 and days 181+ separately in a segmented Cox regression (Figure 2). These segmented analyses satisfied the proportional hazard assumption and showed that mortality risk with pimavanserin was lower than with atypical antipsychotics during the first 180 days of use (hazard ratio=0.65, 95% CI=0.53, 0.79), while mortality risk was similar thereafter (hazard ratio=1.05, 95% CI=0.82, 1.33). Table 2 enumerates the deaths during follow-up by time segment (1–180 days and 181+ days). During the first 180 days of follow-up, the adjusted death rates were 13 and 21 per 100 person-years with pimavanserin and atypical antipsychotic use, respectively. Corresponding rates of death after 180 days were 18 and 17 per 100 person-years for pimavanserin and atypical antipsychotics, respectively. At 75 days of follow-up, which was near the median duration of use for pimavanserin and atypical antipsychotics, 2.4% of pimavanserin initiators and 3.6% of atypical antipsychotic initiators had died, and at 180 days of follow-up, 3.7% of pimavanserin initiators and 5.4% of atypical antipsychotic initiators had died (see Table 2 for weighted results). After 75 days, the number needed to harm for one excess death with atypical antipsychotic treatment versus pimavanserin treatment was 64 (95% CI=38, 220), and at 180 days of use, the number needed to harm was 30 (95% CI=19, 73). The secondary analysis limited to patients whose study drug was prescribed by a neurologist showed results similar to those of the main analysis (days 1–180: hazard ratio=0.61, 95% CI=0.47, 0.80; days 181+: hazard ratio=1.05, 95% CI=0.76–1.45).
FIGURE 1. Kaplan-Meier plot of weighted cumulative survival probability among patients with Parkinson’s disease treated with pimavanserin or atypical antipsychotics
FIGURE 2. Hazard ratio and 95% confidence interval for death among patients with Parkinson’s disease treated with pimavanserin or atypical antipsychotics
TABLE 2. Counts and incidence rates for deaths during follow-up among patients with Parkinson’s disease treated with pimavanserin or atypical antipsychotics
Follow-Up Time PeriodPimavanserinAtypical Antipsychotics
UnweightedWeighteda
Overall   
Number of patients3,22718,4423,251
Number of deaths2071,752248
Follow-up time (person-years)1,3707,2231,267
Incidence rate (per person-year)0.150.240.20
1–75 days   
Number of patients3,22718,4423,251
Number of deaths78841117
Follow-up time (person-years)5252,974525
Incidence rate (per person-year)0.150.280.22
1–180 days   
Number of patients3,22718,4423,251
Number of deaths1181,262176
Follow-up time (person-years)8804,817849
Incidence rate (per person-year)0.130.260.21
181+ days   
Number of patients8974,510788
Number of deaths8949072
Follow-up time (person-years)4912,405418
Incidence rate (per person-year)0.180.200.17
a
Denotes the values for the atypical-antipsychotic weighted pseudo-population created by applying inverse probability of treatment weighting–average treatment effect on the treated (IPTW-ATT) weights for the cohort treated with atypical antipsychotics.
Adjusting for hospice admission as a time-varying covariate eliminated the difference in mortality risk between pimavanserin and atypical antipsychotics (hazard ratio for death prior to hospice care, 0.91, 95% CI=0.72, 1.16; hazard ratio for death after hospice entry, 1.06, 95% CI=0.79, 1.42), consistent with hospice admission being on the pathway to the outcome of death. A total of 311 pimavanserin patients and 2,758 atypical antipsychotic patients died or were admitted to hospice care. The results for this composite outcome agreed with those for death alone, with a lower risk for pimavanserin users compared with atypical antipsychotic users in the first 180 days (hazard ratio=0.60, 95% CI=0.51, 0.70), but no difference in risk after 180 days (hazard ratio=0.97, 95% CI=0.79, 1.20).
A priori subgroup (Figure 3) and sensitivity analyses (Figure 4) were largely consistent with the main analyses, indicating lower mortality risk with pimavanserin use compared with atypical antipsychotic use during the first 180 days and little difference in mortality risk beyond 180 days. However, in the roughly 15% of patients residing in nursing homes, mortality was not lower with pimavanserin use compared with atypical antipsychotic use in the first 180 days (hazard ratio=1.05, 95% CI=0.73, 1.52) and in fact was numerically higher with pimavanserin in the 181+ days period (hazard ratio=1.54, 95% CI=0.91, 2.63) (Figure 3). A statistical test for effect modification by nursing home residence was significant (p=0.008). In a post hoc analysis, we categorized patients by use of antidementia drugs in the baseline period. The difference in mortality between pimavanserin and atypical antipsychotic users was more pronounced in the subgroup treated with antidementia medications. In the first 180 days, hazard ratios for mortality with pimavanserin compared with atypical antipsychotics were 0.49 (95% CI=0.36, 0.67) and 0.82 (95% CI=0.63, 1.06) among users and nonusers of antidementia drugs, respectively. The hazard ratios for mortality during the 181+ days period were 0.92 (95% CI=0.66, 1.28) and 1.22 (95% CI=0.84, 1.75) for antidementia drug users and nonusers, respectively.
FIGURE 3. Hazard ratios and 95% confidence intervals for death among patients with Parkinson’s disease treated with pimavanserin or atypical antipsychotics: subgroup analyses
FIGURE 4. Hazard ratios and 95% confidence intervals for death among patients with Parkinson’s disease treated with pimavanserin or atypical antipsychotics: sensitivity analyses
In one sensitivity analysis, we limited the atypical antipsychotic reference group to users of quetiapine, which was used by 78% of our atypical antipsychotic cohort (N=14,463). Daily quetiapine doses were generally low (79% were 50 mg or less). For pimavanserin versus quetiapine, over the first 180 days the weighted hazard ratio for death was 0.69 (95% CI=0.57, 0.85), and the hazard ratio for 180+ days was 1.09 (95% CI=0.85, 1.40). For pimavanserin versus non-quetiapine atypical antipsychotics, the hazard ratios for death were 0.54 (95% CI=0.37, 0.78) for the first 180 days and 1.24 (95% CI=0.75, 2.07) for 180+ days. The E-value for the observed hazard ratio during the first 180 days of treatment was 2.45, indicating that an unmeasured confounder would have to be more than twice as prevalent among patients using one medication as it was among patients using the other, and would have to alter the risk of death more than twofold, to explain the observed association.

Discussion

In a large, nationally representative population of Medicare beneficiaries with Parkinson’s disease, pimavanserin was associated with 35% lower all-cause mortality compared with atypical antipsychotics over the first 180 days of treatment, which translated to the avoidance of one death for every 30 patients treated with pimavanserin rather than atypical antipsychotics. The relative protective effect of pimavanserin was observed in 85% of the sample who were community-dwelling patients with Parkinson’s disease but not among the 15% of patients residing in nursing homes. Reductions in mortality with pimavanserin compared with atypical antipsychotics were more pronounced over the first 180 days among patients receiving antidementia medication. For the composite outcome of death or transfer to hospice care, there was a 41% lower risk for pimavanserin users compared with atypical antipsychotic users over the first 180 days of treatment.
Our data point to lower mortality with use of pimavanserin compared with atypical antipsychotics, but with two important exceptions. The first involves duration of use, as the relative mortality advantage for pimavanserin was absent beyond 180 days. Further research may be needed to clarify the reasons for that, as high attrition beyond 180 days was a limitation of our study.
The second exception to lower mortality with pimavanserin treatment was among nursing home patients, who comprised 15% of our study population. Mortality was not reduced with pimavanserin relative to atypical antipsychotics among these patients, which was a statistically significantly different finding than the result in community-dwelling patients. Furthermore, our data suggest that pimavanserin could be associated with higher mortality than use of atypical antipsychotics in nursing home patients with longer-term use. However, the relatively small sample size renders these estimates highly uncertain. Reasons for the discrepancy in results between nursing home and community patients are unclear. One possibility is that high baseline mortality among nursing home residents overwhelms the contributions of Parkinson’s disease psychosis treatments to mortality. Additionally, both initiation and continuation of medications for nursing home residents is subject to quality review under the Medicare program, while there is no corresponding third-party review of medication regimens for community patients, so these drugs may be used differently in the two settings.
As noted, a recent study of patients with Parkinson’s disease reported higher mortality among pimavanserin users than pimavanserin nonusers (12). Our results do not conflict with this finding, however, because the reference group was nonusers; as psychotic symptoms in patients with Parkinson’s disease are associated with mortality (10), pimavanserin users would be expected to have higher mortality than pimavanserin nonusers, three-quarters of whom in our sample did not use other antipsychotics.
Our study has several strengths. It was national in scope, which afforded us a large sample of pimavanserin users. Confounding by indication is an important limitation of studies involving comparisons with nonusers, who presumably have a lower prevalence of psychosis. We addressed this weakness by restricting our study to patients with Parkinson’s disease receiving treatment for probable Parkinson’s disease psychosis. We performed subgroup and sensitivity analyses to examine the potential for treatment-effect heterogeneity and to determine whether our results were robust to changes in the specifications of the analysis. Except for patients in nursing homes, our results were consistent across these analyses.
Several limitations of our study should be borne in mind. We were unable to compare mortality risk with pimavanserin to no medication, as we did not have a reliable way to identify patients with Parkinson’s disease with psychosis who were not receiving medication for it. Accordingly, our analyses only address the relative safety of pimavanserin in Parkinson’s disease compared to use of atypical antipsychotics. Importantly, because of the observational, nonrandomized nature of our study, our results may be subject to residual confounding, notwithstanding our analytic strategies to mitigate the effects of potential confounders. Before weighting, certain medical comorbidities were more frequent among atypical antipsychotic users at baseline. Weighting successfully ameliorated between-group disparities in all measured baseline covariates, including important predictors of mortality such as the Charlson comorbidity index and frailty score. The E-value estimate of 2.45 in the first 180 days of treatment suggests that unmeasured confounding was unlikely. Nonetheless, our study was not a randomized trial, and unmeasured confounding as an explanation for our results, although unlikely, cannot be excluded.
Heterogeneity in characteristics of pimavanserin users may have arisen over the study period if patients who received prescriptions for pimavanserin when it was newly marketed differed from those who initiated pimavanserin later. Heterogeneity may also have been present in indications for atypical antipsychotics, which have been used off-label for insomnia or agitation.
Additionally, we were unable to analyze specific causes of death. However, adjusting for hospice entry accounted for the difference in mortality risk, which would be consistent with a more gradual acceleration of the natural trajectory toward death among patients with Parkinson’s disease psychosis, rather than an acute event leading to increased mortality. Finally, with median follow-up times of 78 and 74 days for pimavanserin and atypical antipsychotic users, respectively, we had relatively fewer data on prolonged use of these drugs, although we still had reasonable statistical precision to analyze longer-term use.
In summary, our findings indicate a 35% lower rate of all-cause mortality among pimavanserin-treated patients with Parkinson’s disease in the first 180 days of treatment compared with patients treated with atypical antipsychotics. This effect was absent among patients residing in nursing homes.

Acknowledgments

The authors gratefully acknowledge Dr. Glenn Mannheim for proposing this study.

Footnote

The views expressed here are those of the authors and not necessarily those of the U.S. Food and Drug Administration, the Centers for Medicare and Medicaid Services, or the Department of Health and Human Services.

Supplementary Material

File (appi.ajp.2021.21090876.ds001.pdf)

References

1.
Nussbaum RL, Ellis CE: Alzheimer’s disease and Parkinson’s disease. N Engl J Med 2003; 348:1356–1364
2.
Marras C, Beck JC, Bower JH, et al: Prevalence of Parkinson’s disease across North America. NPJ Parkinsons Dis 2018; 4:21
3.
Factor SA, Feustel PJ, Friedman JH, et al: Longitudinal outcome of Parkinson’s disease patients with psychosis. Neurology 2003; 60:1756–1761
4.
Fénelon G, Alves G: Epidemiology of psychosis in Parkinson’s disease. J Neurol Sci 2010; 289:12–17
5.
Mathis MV, Muoio BM, Andreason P, et al: The US Food and Drug Administration’s perspective on the new antipsychotic pimavanserin. J Clin Psychiatry 2017; 78:e668–e673
6.
Institute for Safe Medication Practices QuarterWatch report. 2017. https://www.ismp.org/quarterwatch/safety-signals-two-novel-drugs
7.
US Food and Drug Administration: FDA analysis finds no new or unexpected safety risks associated with Nuplazid (pimavanserin), a medication to treat the hallucinations and delusions of Parkinson’s disease psychosis. 2018. https://www.fda.gov/Drugs/DrugSafety/ucm621160.htm
8.
Schneider LS, Dagerman KS, Insel P: Risk of death with atypical antipsychotic drug treatment for dementia: meta-analysis of randomized placebo-controlled trials. JAMA 2005; 294:1934–1943
9.
Yuan M, Sperry L, Malhado-Chang N, et al: Atypical antipsychotic therapy in Parkinson’s disease psychosis: a retrospective study. Brain Behav 2017; 7:e00639
10.
Forsaa EB, Larsen JP, Wentzel-Larsen T, et al: What predicts mortality in Parkinson disease? A prospective population-based long-term study. Neurology 2010; 75:1270–1276
11.
Weintraub D, Chiang C, Kim HM, et al: Association of antipsychotic use with mortality risk in patients with Parkinson disease. JAMA Neurol 2016; 73:535–541
12.
Hwang YJ, Alexander GC, An H, et al: Risk of hospitalization and death associated with pimavanserin use in older adults with Parkinson disease. Neurology 2021; 97;e1266–e1275
13.
Moreno GM, Ghandi R, Lessig SL, et al: Mortality in patients with Parkinson disease psychosis receiving pimavanserin and quetiapine. Neurology 2018; 91:797–799
14.
Horn S, Richardson H, Xie SX, et al: Pimavanserin versus quetiapine for the treatment of psychosis in Parkinson’s disease and dementia with Lewy bodies. Parkinsonism Relat Disord 2019; 69:119–124
15.
Quan H, Sundararajan V, Halfon P, et al: Coding algorithms for defining comorbidities in ICD-9 and ICD-10 administrative data. Med Care 2005; 43:1130–1139
16.
Fried LP, Tangen CM, Walston J, et al: Frailty in older adults: evidence for a phenotype. J Gerontol A Biol Sci Med Sci 2001; 56:M146–M156
17.
Segal JB, Chang HY, Du Y, et al: Development of a claims-based frailty indicator anchored to a well-established frailty phenotype. Med Care 2017; 55:716–722
18.
Pennap D, Swain S, Akhtar S, et al: Comparing Centers for Medicare and Medicaid Services (CMS) enrollment database death dates to the National Death Index (NDI). Pharmacoepidemiol Drug Saf 2020; 29(suppl 1):401–402
19.
Austin PC: Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples. Stat Med 2009; 28:3083–3107
20.
Xie J, Liu C: Adjusted Kaplan-Meier estimator and log-rank test with inverse probability of treatment weighting for survival data. Stat Med 2005; 24:3089–3110
21.
Schoenfeld DA: Partial residuals for the proportional hazards regression model. Biometrika 1982; 69:239–241
22.
Altman DG, Andersen K: Calculating the number needed to treat for trials where the outcome is time to an event. BMJ 1999; 319:1492–1495
23.
VanderWeele TJ, Ding P: Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med 2017; 167:268–274
25.
Woosley RL, Heise CW, Romero KA: QTdrugs List. Azcert, Inc, Oro Valley, Ariz. www.Crediblemeds.org
26.
Tibshirani R: The lasso method for variable selection in the Cox model. Stat Med 1997; 16:385–395

Information & Authors

Information

Published In

Go to American Journal of Psychiatry
Go to American Journal of Psychiatry
American Journal of Psychiatry
Pages: 553 - 561
PubMed: 35702829

History

Received: 2 September 2021
Revision received: 20 December 2021
Revision received: 16 February 2022
Accepted: 14 March 2022
Published online: 15 June 2022
Published in print: August 2022

Keywords

  1. Antipsychotics
  2. Epidemiology
  3. Neurology

Authors

Details

Andrew D. Mosholder, M.D., M.P.H. [email protected]
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Yong Ma, Ph.D.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Sandia Akhtar, B.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Gerald D. Podskalny, D.O., M.P.H.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Yuhui Feng, M.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Hai Lyu, M.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Jiemin Liao, M.A.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Yuqin Wei, M.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Michael Wernecke, B.A.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Kira Leishear, Ph.D., M.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Lorene M. Nelson, Ph.D., M.S.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Thomas E. MaCurdy, Ph.D.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
Jeffrey A. Kelman, M.D., M.M.Sc.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).
David J. Graham, M.D., M.P.H.
Division of Epidemiology 1 (Mosholder, Leishear), Division of Neurology 1 (Podskalny), Office of Pharmacovigilance and Epidemiology (Graham), and Division of Biometrics 7 (Ma), U.S. Food and Drug Administration Center for Drug Evaluation and Research, Silver Spring, Md.; Acumen LLC, Burlingame, Calif. (Akhtar, Feng, Lyu, Liao, Wei, Wernecke, Nelson, MaCurdy); Guardant Health, Redwood City, Calif. (Liao); Department of Epidemiology and Population Health (Nelson) and Department of Economics (MaCurdy), Stanford University, Stanford, Calif.; Centers for Medicare and Medicaid Services, Washington, DC (Kelman).

Notes

Send correspondence to Dr. Mosholder ([email protected]).

Funding Information

The authors report no financial relationships with commercial interests.

Metrics & Citations

Metrics

Citations

Export Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

For more information or tips please see 'Downloading to a citation manager' in the Help menu.

Format
Citation style
Style
Copy to clipboard

View Options

View options

PDF/EPUB

View PDF/EPUB

Login options

Already a subscriber? Access your subscription through your login credentials or your institution for full access to this article.

Personal login Institutional Login Open Athens login
Purchase Options

Purchase this article to access the full text.

PPV Articles - American Journal of Psychiatry

PPV Articles - American Journal of Psychiatry

Not a subscriber?

Subscribe Now / Learn More

PsychiatryOnline subscription options offer access to the DSM-5-TR® library, books, journals, CME, and patient resources. This all-in-one virtual library provides psychiatrists and mental health professionals with key resources for diagnosis, treatment, research, and professional development.

Need more help? PsychiatryOnline Customer Service may be reached by emailing [email protected] or by calling 800-368-5777 (in the U.S.) or 703-907-7322 (outside the U.S.).

Media

Figures

Other

Tables

Share

Share

Share article link

Share