Abstract
Background Psychotic symptoms are reported as one of the most common complications of Alzheimer’s disease (AD), affecting approximately half of AD patients, in whom they are associated with more rapid deterioration and increased mortality. Empiric treatments, namely first and second-generation antipsychotics, confer modest efficacy in AD patients with psychosis (AD+P) and themselves increase mortality. A recent genome-wide meta-analysis and early clinical trials suggest the use and beneficial effects of antidepressants among AD+P patients. This motivates our rationale for exploring their potential as a novel combination therapy option amongst these patients.
Methods We included University of Pittsburgh Medical Center (UPMC) electronic medical records (EMRs) of 10,260 AD patients from January 2004 and October 2019 in our study. Survival analysis was performed to assess the effects of the combination of antipsychotics and antidepressants on the mortality of these patients. To provide more valuable insights on the hidden mechanisms of the combinatorial therapy, a protein-protein interaction (PPI) network representing AD+P was built, and network analysis methods were used to quantify the efficacy of these drugs on AD+P. An indicator score combining the measurements on the separation between drugs and the proximity between the drugs and AD+P was used to measure the effect of an antipsychotic-antidepressant drug pair against AD+P.
Results Our survival analyses replicated that antipsychotic usage is strongly associated with increased mortality in AD patients while the co-administration of antidepressants with antipsychotics showed a significant beneficial effect in reducing mortality. Our network analysis showed that the targets of antipsychotics and antidepressants are well-separated, and antipsychotics and antidepressants have similar proximity scores to AD+P. Eight drug pairs, including some popular recommendations like Aripiprazole/Sertraline and other pairs not reported previously like Iloperidone/Maprotiline showed higher than average indicator scores which suggest their potential in treating AD+P via strong synergetic effects as seen in our study.
Conclusion Our proposed combinations of antipsychotics and antidepressants therapy showed a strong superiority over current antipsychotics treatment for AD+P. The observed beneficial effects can be further strengthened by optimizing drug-pair selection based on our systems pharmacology analysis.
1. Introduction
Approximately more than 50% of Alzheimer’s disease (AD) patients experience psychotic symptoms and other neuropsychiatric complications 1. AD patients with psychosis (AD+P) are considered a subgroup of patients who experience more severe symptoms including greater cognitive impairment and a quicker cognitive decline 2. AD+P is also associated with higher rates of co-occurring agitation 3, aggression 3, 4, depression 5, 6, mortality7, functional impairment 8, and increased caregiver burden 9 than AD patients without psychosis (AD-P).
Currently, no specific medication has been approved by FDA for AD+P. Second-generation antipsychotics (SGAs) have been widely used and recommended by geriatric experts in the management of psychosis in AD 10-12. However, their usage is greatly limited by increased risk of adverse events and the likelihood of co-morbid health problems 13, 14. This prompted the FDA to issue a “black-box” warning in 2005 to highlight the increased risk of mortality for patients with dementia who are treated with SGAs 15. Meanwhile, antipsychotics have demonstrated modest efficacy in treating psychosis, aggression, and agitation in individuals with dementia 16-18. This emphasizes the need to develop safer and more effective treatment options for AD+P.
Based on our previous studies, antipsychotics do not engage the underlying biology of AD+P, and therefore their modest effectiveness is unsuprising19. In order to identify safe and effective treatments for AD+P, it is essential to have a comprehensive understanding of the underlying biology of AD+P. A recent study reported that the heritability of psychosis in AD is estimated to be 61%, thus suggesting a strong association between AD+P and genetic variations20. Another study performed a large genome-wide meta-analysis on 12,317 AD subjects with or without psychosis 21. The authors reported that AD⍰+ ⍰P was not significantly genetically correlated with schizophrenia, but it was negatively correlated with bipolar disorder and positively correlated with depression. These associations provide a biologic rationale for repurposing antidepressant agents as novel treatment options for AD+P in our current study.
To answer if antidepressants are effective in managing neuropsychiatric symptoms in AD patients, nine clinical trials involving 692 patients were conducted. Five of them compared antidepressants with placebo and 4 compared with antipsychotics 22-32. However, only two selective serotonin reuptake inhibitors (SSRIs) sertraline (Zoloft) and citalopram (Celexa) were studied and the antipsychotics in the study were typical antipsychotics (haloperidol, perphenazine) while only 1 trial studied an SGA, risperidone. Among the five studies comparing SSRIs with placebo, two of them reported a significant benefit for Citalopram against AD+P 22, 24. Meanwhile, no significant difference was reported between the efficacy of SSRIs and risperidone. Therefore, testing more antidepressants, especially other classes of antidepressants such as serotonin-noradrenaline reuptake inhibitors, tricyclic antidepressants, and monoamine oxidase inhibitors, may be worthwhile to provide a better understanding of the impact of antidepressants on AD+P.
While antipsychotics and antidepressants both have shown beneficial effects against AD+P, they have different targets which in turn can modulate different biological pathways. Thus, the combination of these drugs can potentially provide multiple advantages like enhanced efficacy, decreased dosage with an equal or increased level of efficacy, and delayed development of drug resistance 33. Due to the excessive time and cost it takes to clinically test the drug combination effects, exhaustive computational methods can be used to predict drug synergy. By integrating information from drugs and diseases we can obtain a comprehensive picture of the potential synergetic effects of these drug combinations. The goal of our study is to further identify key combinations of antipsychotics and antidepressants that possess potential synergetic effects against AD+P with the help of our state-of-the-art quantitative systems pharmacology approaches.
2. Methods and Material
2.1 Dataset collection
To explore the beneficial effect of combination therapy, we examined the data from January 2004 to October 2019 from the Neptune system at the University of Pittsburgh Medical Center (UPMC), which manages the use of patient EMRs from the UPMC health system for research purposes (rio.pitt.edu/services)34. The database includes demographic information, diagnoses, encounters, medication prescriptions, prescription fill history, and laboratory tests. AD patients were identified using ICD9/10 codes (331.0, G30.0, G30.1, G30.9) and the onset of psychosis were defined by ICD9/10 codes (780.1, F06.0, R44.2, R44.1, R44.3, R44.0, 298.8, F22, F23, F28, F29, 293.82, 298.9, 290.11, 293, 290.3) based on the suggestions from UPMC clinicians.
2.2 Survival analysis
We included patients who met the following inclusion criteria: 1) Patient had an AD diagnosis; 2) Patients did not take antidepressants nor antipsychotics one year prior to the diagnosis of AD. Nine comorbidities, including MDD (major depressive disorder), stroke, COPD (chronic obstructive pulmonary disease), ASCVD (atherosclerotic cardiovascular disease), T2DM (type 2 diabetes), HTN (hypertension), CKD (chronic kidney disease), HF (heart failure) and cancer, were considered as confounders in our survival analysis (ICD codes listed in supplementary material) 35, 36. The time origin for each patient in the survival analysis is the first AD diagnosis date and time to all-cause death is the outcome. Patients are marked with the above comorbidities if they were diagnosed before the AD diagnosis. Only medications that are prescribed to the same patients more than 2 times with more than 30 days apart were considered to eliminate short-term usage during hospitalization. The records of patients up to 5 years after the AD diagnosis were used in the analysis.
Time to all-cause death was constructed as the time between the first date of AD diagnosis and death. Patients who were alive by the end of 5 years since AD diagnosis were censored. Survival analysis was performed to evaluate the association between medications and mortality. To accommodate the change of drug usage during the follow-up, we fitted a time-dependent Cox’s proportional hazards model37, 38 with antipsychotic drug effect (yes or no) and antidepressant drug effect (yes or no) as time dependent covariates. Specifically, the 5-year follow-up period was divided into 60 months, and we assumed the drug effect from one prescription will last 2 months which covers 2 intervals in our study. As shown in Figure 1, if a patient had an antipsychotics prescription in the first month, we consider the patient under drug effect for that month and the month after. If a patient was prescribed both antipsychotics and antidepressants (boxed in Figure 1), we consider the patient under combinational therapy. The drug effect on one patient may change over time between four statuses: no drug, antipsychotics only, antidepressants only, and the combination. Baseline demographics and comorbidities were also included in the model. Contrasts between different drug groups were performed with hazard ratios and p values reported. The data were analyzed using both R (version 4.1.0) and python (version 3.7.12) packages.
Schematic diagram for identifying drug usage status of subjects. In the upper two rows of Figure 1, the sections with 1 indicate that there are antipsychotics/antidepressants prescriptions in that month while 0 means there are no prescriptions for the two kinds of medications. In the lower rows, the sections under 1 are extended for 1 month to reflect the drug effect thus these markers show the time that the subject was under drug effect.
2.3 Prediction of synergetic effect among antipsychotic-antidepressant pairs
Since the analysis conducted are a mixed effect of any combination of antipsychotics and antidepressants, we further examined if certain antipsychotic-antidepressant pairs may possess higher synergetic effect in managing psychotic symptoms in AD. We first attempted this evaluation by using EMR data. However, because of the large number of antipsychotics and antidepressants used clinically, the appearances of specific drug pairs are limited (most co-administrated pair: Quetiapine and Sertraline, 408 times. More data in supplementary material ST1), and their overlap times are too short to produce enough statistical power. Therefore, we turned to systems pharmacology approaches to identify antipsychotic-antidepressant pairs that are most suited in treating AD+P.
For the systems pharmacology study, the postsynaptic density (PSD) proteome was used to build the AD+P network. This proteomic signature was generated by Dr. Sweet’s team39. Information about antipsychotics, antidepressants and their targets were extracted from DrugBank (https://www.drugbank.ca/)40. The pharmacological action label of a drug provides information about whether binding to a target contributes to the pharmacological effects. PPI data was collected from STRING (https://string-db.org/)41. The PPI networks were constructed and analyzed using the python package networkx (https://networkx.github.io/)42. The interaction network was shown in the molecular action view with medium confidence level (> 0.4) 43. AD+P-related proteins were joined with targets of antipsychotics and antidepressants to construct the disease-target network. In addition, we also included the proteins bridging proteins from disease module and proteins from the target module in our disease-targets networks. Gene signature data were used to calculate the proximity score for drugs and AD+P. The post-treatment gene signature data were obtained from the LINCS L1000 database44.
To predict potential drug combinations for AD+P, we adopted the methods from Chen, S., et al. 45 and modified them by incorporating differentially expressed genes (DEGs) after drug treatment to minimize the bias caused by module sizes. Based on previous studies, for a drug pair to have a therapeutic effect on a disease, both target modules (green and yellow circles in figure 2) of the two drugs must overlap with the disease module (pink circle in figure 2) 45. In addition, the two target modules need to be overlapped with the disease module independently to form a complementary exposure to have synergetic effects with each other as shown in Figure 2. To be specific, the targets of the two drugs both need to be overlapped with the disease module in the PPI network, but these two target modules can’t overlap 45. Therefore, two network approaches are applied to predict the possible drug combinations for AD+P: (1) network-based separation between targets of two drugs 46; (2) gene signature-based proximity between the disease (AD+P) module and the target modules of the two drugs 47.
Schematic diagram for the network-based complementary exposure relationship between two drug–target modules and one disease module on a drug–drug–disease combination (adopted from Guney, E.; Menche, J.; Vidal, M.; Barábasi, A.-L., Network-based in silico drug efficacy screening. Nature communications 2016, 7 (1), 1-13. 47). Big red circle: AD+P modules composed of AD+P related proteins/genes (small red rounds). Green and yellow circle: Drug modules composed of drug targets.
Separation evaluation
The separation score (SAB) of drug modules A and B are calculated for all possible combinations between antipsychotics and antidepressants. The separation score (SAB) of drug modules A and B can be calculated as:
where <dAA>, <dBB> and <dAB>, are the mean shortest distance for genes within each module. It compares the mean shortest distance between the targets of each drug. For better understanding, if SAB < 0, it means that the targets in the two drug modules are in the same network neighborhood which is not separated; if SAB ≥ 0, it means that the two drug modules are topologically separated from each other. We filter the combinations based on their ability to achieve complementary exposure with the AD+P module.
Proximity evaluation
Level 5 LINCS L1000 data, a collection of gene expression profiles for thousands of perturbagens at a variety of time points, doses, and cell lines, were downloaded from the GEO database (accession numbers: GSE70138 and GSE92742). Gene expression profiles were included only if they are tested on a cell line of central nervous systems and their dose should be beyond 1 uM. To identify genes that are significantly differentially expressed in the data, their Z scores from multiple tests were averaged and if their |Z|>1, the genes are considered as significant for a drug48.
The association between the drug and AD+P was quantitatively evaluated with Signed Jaccard Score49. The index ranges from +1 to −1, where +1 and −1 indicate a same pattern and an inverse pattern of two identical gene sets, respectively. Zero indicates that the two sets have no overlap, or the positive and negative correlations cancel out.
3. Results
After applying all inclusion criteria to our dataset, 10,206 unique AD patients were included, and their baseline characteristics are shown in supplementary material (ST2).
3.1 Use of antipsychotics in AD patients is associated with increased mortality
cPrior literature has reported that the use of antipsychotics is associated with increased mortality in AD patients50-52; we therefore first sought to replicate this finding to validate the integrity of our methodological approach. The time varying Cox model was conducted with our data to test the significance.
Table 1 Multivariate Cox regression analyses of association between antipsychotics and all-cause mortality in AD patients As indicated in Table 1, antipsychotic usage is significantly associated with increased mortality in AD patients (HR=2.47, p<0.001). The results showed in this table are in accordance with current literature which suggest that our data and methods are valid for further analysis.
3.2 Survival analysis revealed significant beneficial effect combining antidepressants and antipsychotics in AD patients
After knowing that antipsychotics may increase mortality in AD patients while antidepressants may decrease mortality53, the effect of a combination therapy comes into play. We would like to examine the protective effects of adding antidepressants to the existing antipsychotics therapy. We performed another survival analysis to examine three mutually exclusive medication use groups: antipsychotics only, antidepressants only and combination.
The results are shown in Table 2 and the combination group is the reference group in this model. Based on the results from Table 2, the combination group showed a significant beneficial effect relative to antipsychotics only group (HR=0.654, p=0.012), which means that combining antidepressants with antipsychotic treatment was associated with significantly protective effects in AD patients, reducing mortality. In addition, marginal significant difference was observed between no drug group and combination group (HR=1.294, p=0.056), which means that by using combination therapy, the increase in mortality due to using antipsychotics was mitigated to some extent in these patients.
Table 2 Multivariate Cox regression analyses of association among treatments and all-cause mortality in AD patients
Based on our findings shown in Table 1 and Table 2, we can conclude that by combining antipsychotics and antidepressants, we can significantly mitigate the increase in mortality associated with antipsychotics. For more direct comparison between different treatment groups, a table with pair-wise comparison among groups is included in supplementary table 2 (ST3).
In addition to these results, we were interested to see if the effect of the treatments will change over time. Therefore, we conducted 6 analyses with follow-up times ranging from 1 to 6 years. We compared the effects of three different treatments (antipsychotics only, antidepressants only, no drug) to the drug combination group, their effects showed moderate fluctuation within the first 3 years and stabilized after 4 years (SF1). In comparison to co-administration of antidepressants and antipsychotics, the antidepressants only group showed a consistent lower mortality throughout the 6 years. While the antipsychotics only group had marginally increased mortality compared to the combination group at the first 3 years, subsequent worse outcomes were clearly evident in 4, 5 and 6 years of follow-up. Finally, the combination group showed comparable effects with the patients with no drug treatment throughout the 6 years period and demonstrated almost significant beneficial effects in 4, 5 and 6 years of the follow-up. Our results are attached in supplementary materials (ST4).
3.2 Systems pharmacology studies provided possible explanations for the synergetic effects and proposed drug pairs may achieve optimized efficacy
In total, 21 antipsychotics and 17 antidepressants commonly used in the clinic are included in our study along with 75 targets for antipsychotics and 32 targets for antidepressants. The PPI network was built with 240 AD+P proteins, targets for antipsychotics and antipsychotics. A PPI network with 321 nodes and 1,363 edges was generated. A total of 357 pairs of antipsychotics and antidepressants are evaluated in the network and their separation scores are calculated as shown in supplementary figure (SF2).
We found that some antidepressants showed great separation (Vortioxetine, Vilazodone, Mirtazapine, Maprotiline), and most drugs pairs showed a separation score above 0 (blue). This suggests an existing difference in the mechanism which can be the key condition to the synergetic effect in the combinational therapy.
To evaluate the proximity between AD+P and medications, 148 and 78 eligible expression profiles for antipsychotics and antidepressants were collected based on the inclusion criteria described earlier. Signed Jaccard scores were calculated between them and AD+P protein expressions.
Post-treatment gene expression data for 16 antipsychotics and 13 antidepressants were exacted, and their Signed Jaccard scores were calculated accordingly and are shown in the tables below.
Table 3 Signed Jaccard scores of antipsychotics and antidepressants with AD+P
To comprehensively evaluate the potential for these drug combinations, separation scores and proximity scores were normalized to a [-1, 1] interval and combined. A combined score was calculated for every drug pair by subtracting two proximity scores of the drugs from their separation score. The combined scores for drug pairs are shown in the heatmap below (Figure 5).
As shown in Figure 3, most drug pairs showed a combined score around 0 which indicate that the synergistic effect between antipsychotics and antidepressants may not be easily achieved within the two drug categories. However, some drugs showed promising results through the panel like Aripiprazole and Thioridazine. For drug combinations that may possess beneficial synergistic effect, they must meet the following criteria: 1) The separation score between two drugs > 0; 2) The proximity scores of the two drugs <0. Table 4 summarized the drug combinations that met the criteria.
Combined scores for antipsychotics and antidepressants combinations. Though most drug pairs showed a combined score around 0, some drugs do show promising scores across the board including Aripiprazole and Thioridazine.
Table 4 Antidepressants and antipsychotics combinations with highest combined scores
Among the nominated drug combinations, Aripiprazole is the only SGAs and is also the most recommended treatment options for AD+P. As a matter of fact, the combination of Aripiprazole and Sertraline has already been tested in clinical trials and showed superior efficacy in treating major depressive disorders (MDD) with no reports of safety concerns54. However, none of these drug pairs have been tested for AD+P in clinical trials. As a matter of fact, no combination drug therapy has ever been tested for AD+P at present.
4. Discussion and conclusion
In this study, we applied a combination of different state-of-the-art systems pharmacology techniques to explore the potential synergetic effect of combining antipsychotics and antidepressants in treating AD+P. This study incorporated different data types including clinical EMR data, protein expressions, post-treatment gene expressions and protein-protein interaction networks. Our results indicate that the combination of antipsychotics and antidepressants may be a safer and more efficient treatment option for AD+P and provided informative insights for drug pairing choices for future therapies.
When dealing with real-world data, like EMRs, there is always a challenge that the compliance of patients presented in the EMR will not be as ideal as we get from a carefully performed clinical trial. In this study, by analyzing real-world EMR data through Cox model with time-dependent covariates, we were able to accommodate the complex usage patterns and allowed the maximum utilization of the data. The beneficial effect of antidepressants in AD patients were reported by multiple studies55-58, though its mechanisms remain unclear. We also found a strong signal in our results (Table 3 and Figure 3) which further substantiated our claim that antidepressants may aid in reducing mortality in AD patients. Our finding provided the fundamental support necessary for our hypothesis that by combining antipsychotics and antidepressants, we can decrease the severe side effect that is constraining the use of antipsychotics in AD therapy.
In order to provide mechanistic support for our observations and identify the most potent drug combinations for treating AD+P, we took advantage of multiple categories of data including PPI network, post-treatment gene expression profiles to quantitatively evaluate the potential synergistic effect of antipsychotics and antidepressants. Our analysis yielded several pairs of drugs that may possess better synergistic effects in treating AD+P. As shown in Table 6, four antipsychotics: Aripiprazole, Thioridazine, were reported, and seven antidepressants: Sertraline, Maprotiline, Nefazodone, Mirtazapine, Trazodone, Paroxetine and Bupropion were mentioned. Between the 2 antipsychotics, Thioridazine is first-generation antipsychotic and it was withdrawn worldwide in 2005 due to its association with cardiac arrythmias59. For Aripiprazole, it has been reported with significant better efficacy in AD patients against psychological symptoms60 which further consolidate our conclusions. For the 7 antidepressants, Nefazodone is discontinued on 2004 because of its association with drug-induced hepatic injuries61. For the rest 6 antidepressants, they belong to 4 classes. Sertraline, Trazodone, Paroxetine are selective serotonin reuptake inhibitors (SSRI) which is the most used antidepressant class62. Maprotiline is a tetracyclic antidepressant with similar pharmacological properties to tricyclic antidepressants (TCAs), it can inhibit neuronal norepinephrine reuptake, possesses some anticholinergic activity, and does not affect monoamine oxidase activity63. Mirtazapine is a tetracyclic piperazino-azepine antidepressant with its effect can be observed as early as 1 week after beginning therapy64, it has also been reported to be efficacious in the off-label management of various other conditions. It may improve the symptoms of neurological disorders, reverse weight loss caused by medical conditions, improve sleep, and prevent nausea and vomiting after surgery65. Bupropion is a norepinephrine/dopamine-reuptake inhibitor (NDRI) antidepressant, and it is a unique option for the treatment of MDD as it lacks any clinically relevant serotonergic effects, typical of other mood medications, or any effects on histamine or adrenaline receptors66, 67. As for combination therapy consisting of antipsychotics and antidepressants, though they were never tested specifically for AD or AD+P, multiple studies have tested their safety and efficacy profile against other disorders. For example, a meta-analysis consisting of eight randomized, placebo-controlled studies reported that antidepressant-antipsychotic cotreatment was superior to monotherapy with either drug class in the acute treatment of psychotic depression68 and another study reported that adding SGAs to antidepressants yielded highly significant superiority in treating MDD, too69. In that study, Aripiprazole, Olanzapine, Risperidone, and Ziprasidone were found to be more effective than other SGAs69.
The results of our study aims to provide a comprehensive and quantitative overview of the underlying relationship among antipsychotics, antidepressants, and AD+P. Our results supported the efficacy of antipsychotics and suggested the most promising antidepressants such as Sertraline and Maprotiline, can be added as supplementary treatment. In addition, since they are all marketed drugs and some of their combinations are already tested by clinical trials for other indications, safety profile will not be a major concern when proposing their long-term usage as an alternative treatment option for AD patients.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
5. Author Contribution
P.F. wrote the manuscript; P.F., Z.L., Y.D. and L.W. designed the research; P.F. and Z.L. performed the research; P.F. and Z.L. analyzed the data; J.K., J.S., J.K. and R.A.S. contributed new reagents/analytical tools.
6. Funding
This study is supported by the following National Institutes of Health grants:
S10 OD028483/OD/NIH HHS/United States
P50 AG005133/AG/NIA NIH HHS/United States
R01 AG027224/AG/NIA NIH HHS/United States
R01 MH116046/MH/NIMH NIH HHS/United States
P30 AG066468/AG/NIA NIH HHS/United States
UL1 TR001857/TR/NCATS NIH HHS/United States
7. Study Highlights
What is the current knowledge on the topic?
Current treatments to AD+P are restricted by their modest efficacy and poor safety profiles. Better treatment options are in urgent need.
What question did this study address?
Will the combination of antipsychotics and antidepressants be a better treatment option? Which antipsychotics and antidepressants are the best pair?
What does this study add to our knowledge?
This study proved that the combination of antipsychotics and antidepressants showed superiority over antipsychotics monotherapy and proposed several drug pairs that may possess higher efficacy.
How might this change drug discovery, development, and/or therapeutics?
Our results supported the efficacy of antipsychotics and suggested the most promising antidepressants such as Sertraline and Maprotiline, can be added as supplementary treatment. In addition, since they are all marketed drugs and some of their combinations are already tested by clinical trials for other indications, safety profile will not be a major concern when proposing their long-term usage as an alternative treatment option for AD patients.
Footnotes
Peihao Fan has no conflict of interest to report.
Lang Zeng has no conflict of interest to report.
Ying Ding has no conflict of interest to report
Jonathan C. Silverstein has no conflict of interest to report.
Joshua Krivinko has no conflict of interest to report.
Robert A Sweet has no conflict of interest to report.
LiRong Wang has no conflict of interest to report.
Email: sweetra{at}upmc.edu
References
