Deep Learning for Alzheimer’s Disease Drug Repurposing using Knowledge Graph and Multi-level Evidence

Developing drugs for treating Alzheimer’s disease (AD) has been extremely challenging and costly due to limited knowledge on underlying biological mechanisms and therapeutic targets. Repurposing drugs or their combination has shown potential in accelerating drug development due to the reduced drug toxicity while targeting multiple pathologies. To address the challenge in AD drug development, we developed a multi-task machine learning pipeline to integrate a comprehensive knowledge graph on biological/pharmacological interactions and multi-level evidence on drug efficacy, to identify repurposable drugs and their combination candidates. We developed and computationally validated a heterogeneous graph representation model with transfer learning from universal biomedical databases and with joint optimization with AD risk genes. Using the drug embedding from the heterogeneous graph representation model, we ranked drug candidates based on evidence from post-treatment transcriptomic patterns, mechanistic efficacy in preclinical models, population-based treatment effect, and Phase II/III clinical trials. We experimentally validated the top-ranked candidates in neuronal cells, identifying drug combinations with efficacy in reducing oxidative stress and safety in maintaining neuronal viability and morphology. Our neuronal response experiments confirmed several biologically efficacious drug combinations (e.g., Galantamine + Mifepristone). This pipeline showed that harmonizing heterogeneous and complementary data/knowledge, including human interactome, transcriptome patterns, experimental efficacy, and real-world patient data shed light on the drug development of complex diseases.


INTRODUCTION
Developing drugs for treating Alzheimer's disease (AD) has been extremely challenging and costly. The total cost of developing new AD drugs, including failures, is estimated at $5.7 billion -seven times more than the cost of developing cancer medicines. (1) The most recent FDA approval of Aducanumab is an amyloid beta-directed monoclonal antibody, which aims to treat patients by reducing the buildup of βamyloid.
(2) Despite showing some promising efficacy in removing β-amyloid, only patients with mild stage and amyloid burden can be benefited from the expensive treatment (estimated to cost $56,000 per year).(3) To develop potential treatments in a timely and cost-effective manner, drug repurposing has shown good potential due to the reduced risk of drug toxicity of previously approved drugs. What is even more promising is to study combinatorial drug repositioning, which has demonstrated their potential for synergistically treating complicated diseases, including cancer,(4) diabetes,(5) metabolic syndrome, (6) cardiovascular disease, (7) thanks to their ability to target multiple pathologies.
Most current AD drug repurposing investigates various perspectives, including drug perturbation transcriptome profiles, network pharmacology, and treatment effects in real-world patient data.
Transcriptomic-based strategy compares drug-induced gene expression with gene expression in AD specimen, (8)(9)(10) which captures integrated molecular changes in different AD pathology. Network pharmacology is another approach that represents a drug's multi-target capacity in a human interaction network. (11)(12)(13) This method aims to identify hidden interaction among drugs and proteins (or disease targets) by estimating the proximity between entities in the network. Alternatively, the real-world data approach leverages large sets of heterogeneous patients' data to obtain off-label indication via counterfactual treatment effect estimation. (14,15) Among these drug repurposing strategies, each strategy captures different aspects in multimodality and multiscale of the AD treatment landscape. These different approaches provide disparate but complementary evidence into potential drugs. (16) Success in a single aspect, however, cannot guarantee clinical effectiveness due to AD's heterogeneous pathogenesis and associated comorbid conditions. Consequently, there is a critical need to integrate diverse evidence of transcriptomic observations, heterogeneous biological interactions, preclinical experimental results, and real-world observation to target multiple perspectives in AD drug development.
Thus, our objective is to develop a multi-task learning pipeline for AD drug repurposing and drug combination prediction. We first built an AD knowledge graph that integrates the heterogeneous human interactome of drugs, genes, pathways, and functional ontology directly or indirectly related to AD. A knowledge graph representation learning was applied to encapsulate the entities and heterogeneous interactions into embedding. Next, repurposable drug candidates were identified based on multi-level evidence from drug perturbation transcriptome profiles, clinical trial, and preclinical trial history, as well as population-based treatment effect estimation. We then mechanistically validated the efficacy and safety of our computationally predicted drug candidates using mouse brain cell lines. Our cell-based validation showed that our proposed framework can predict highly promising repurposable drug combinations.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

AD knowledge graph representation
We integrated various sources of information from several human-curated databases to build a comprehensive AD knowledge graph using drugs, genes, pathways, and gene ontology, together with edges for drug-target interaction, gene-gene interactions, drug-drug structural similarity, pathway-gene, GO-gene, and GO-drug associations (Fig. 1a, Method 1). Using the AD knowledge graph, we encapsulated the molecular interactions into representation (embedding) using the knowledge graph representation method we designed (Method 2). The AD knowledge graph we built is specific to AD-related molecular interactions. Underlying AD etiology is an open research area and current knowledge of AD-related interactions is likely incomplete. To address the incomplete knowledge of AD, we utilized transfer learning from the universal biological database representation (Method 3). Additionally, we encouraged the gene embedding in the knowledge graph to distinguish high-risk AD-related genes, by jointly optimizing the target classification task and the graph representation (Method 4). After training the model to learn heterogeneous relationships in the AD knowledge graph, we evaluated whether the embeddings faithfully condense the AD knowledge. We investigated whether the node embedding can restore the local interactions (Method 5, Table 1) and visually examined the global position of the node embedding via UMAP plot (Fig. 2). (20) Overall, the interaction (edge) prediction accuracy of embeddings with transfer learning and fine-tuning achieved AUROC=0.991 and AUPRC=0.988, which is higher than DRKG's universal embedding (AUROC=0.565, AUPRC=0.648) and AD knowledge graph representation without transfer learning (AUROC=0.959, AUPRC=0.958) ( Table 1). This higher accuracy with transfer learning implies that the universal pharmacological interactions encompassing the ADspecific molecular interactions can complement the AD knowledge graph representation. By ablating one edge type and measuring the edge prediction accuracy (Method 5), we found that the drug-target or drugpathway interactions are the most unique interactions that are hard to infer if they are excluded (Table S1).
For the AD-related gene prediction task (Method 4), the gene embedding was well trained to distinguish AD-related genes from other remaining genes (Table S2). Visualizing the embedding via UMAP, we found . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint that the overall node distribution is generally dispersed throughout the plot rather than clustered with the same node type. Undergoing clinical trial drugs are aggregated in three main groups, each with similar therapeutic or pharmacological categories, implying the node embedding well reflects the contextual information. In all, the high accuracy in missing interaction prediction and AD-related gene prediction, together with global distribution in UMAP plot, implies that the node embedding faithfully captures the local and global network topology of the AD knowledge graph regardless of edge types, with the help of universal embedding that provides complementary information to restore the masked AD knowledge. This implied the potential applicability of our derived embeddings for drug repurposing tasks.

Identifying drug candidates using transcriptomic, mechanistic, and epidemiological evidence
After we synthesized the heterogeneous AD knowledge into the embedding, we identified repurposable drugs using the drug embedding that encapsulates molecular interactions and pharmacological interactions.
Our strategy is to prioritize drugs with multi-level evidence on the drug's efficacy, such as transcriptomic reversed patterns, mechanistic efficacy, and population-based treatment effects (Method 7, Table S3), using multi-task learning (Method 6). As a result, the ranking model accuracy was AUROC=0.910 and AUPRC=0.472 to predict drugs with at least one of the drug efficacy labels. We presented the top-ranked drugs in Table S4 (precision at top 600=0.42) and highlighted the top 10 in Table 2.

Identifying drug combination
As indicated by the complexity of the AD interactome network, using single drugs to treat AD might result in limited effects. To improve treatment efficacy, we identified potential drug combinations from the top 600 ranked drugs and existing FDA-approved AD drugs. The FDA has approved four drugs for the treatment of symptoms and memory loss in AD -rivastigmine, galantamine, donepezil, memantine, and memantine combined with donepezil. These drugs have two main mechanisms of action, as cholinesterase (AChE) inhibitors or N-methyl D-aspartate (NMDA) antagonists. Hippocampal cholinergic hypofunction is closely related to the cognitive deficits of Alzheimer's disease. These pharmacological treatments are . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint intended to maintain or stabilize neuronal function, management of behavioral symptoms, and slow down the rate of memory loss by regulating neurotransmitters. In general, AChE inhibitors are recommended for the treatment of mild-to-moderate AD, whereas NMDA antagonist is prescribed to treat moderate-to-severe AD.(32) Unfortunately, after an expected period of daily administration of these AChE inhibitors and NMDA antagonist drugs, some side effects are experienced by AD patients, including nausea, vomiting, diarrhea, indigestion, muscle cramps, fatigue, loss of appetite, weight loss, dizziness, headache, and confusion. Here, we proposed to discover alternative and effective treatments for AD by combining the FDA-approved AD drugs with the top-ranked drugs identified in our computational model.
To identify drug combinations with synergistic interactions without degradation in safety, we leveraged the "Complementary Exposure pattern" (Method 8) indicating that a drug combination is therapeutically effective if the targets of the drugs hit the disease module without overlap. (33) We presented the identified drug combinations in Table S5.

In vitro validation on safety, tolerability, and neuronal responses of drug combinations
We mechanistically validated the top-ranked drug combinations by measuring cytotoxicity and oxidative stress on hippocampal neurons (Method 9, Fig. 1f). We identified significant efficacy in multiple drug combinations involving existing and novel drug candidates for AD.
Memantine functions as an antagonist of glutamatergic NMDA receptors, type 3 serotonin (5HT3) receptor, and nicotinic acetylcholine receptor (nAChRs), allowing neuronal availability of neurotransmitters and improving cognition and behavior in moderate-to-severe AD patients. In combination with Tolcapone, an FDA-approved catechol-O-methyltransferase inhibitor used in the management of Parkinson's disease to increase levels of peripheral dopamine, we observed safety and tolerability in cholinergic neurons during the dose-and time-course treatments. The production of reactive oxygen species (ROS) was reduced by 5.6±2.2% with Tolca, and in the presence of Aβ, only the combined dosage of Mem+Tolca was able to reestablish oxidative stress to basal levels, whereas the single Mem treatment presented a 5.5±2.7% increase (Fig. 3). Memantine is usually administered in combination with donepezil, but in combination with . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Tolcapone, the drugs could have significant benefits on global and functional outcomes in AD subjects. (34,35) Donepezil, an AChE inhibitor, allows an increase in acetylcholine in cholinergic neurons. This compound has been used for 25 years to delay cognitive decline and symptoms of AD. High-rank drugs with potential activity tested included octyl-methoxycinnamate, nebivolol, desogestrel, and salicylic acid. The combination of Don+octyl-methoxycinnamate showed high tolerability and safety. Functionally, donepezil effectively reduced the amount of ROS by 6.82-7.98%, but its combination with octyl-met showed a slight increased by 4.3-9.1% (Fig. 3), this could be related to the limited solubility of this active compound and that it is usually employed in topical formations.
Galantamine is a weak AChE inhibitor and allosteric potentiator of acetylcholine receptors, intended to treat the progression of symptoms, memory loss, and thinking ability in mild-to-moderate AD. Gal prevents acetylcholine processing and stimulation of its release by nicotinic receptors. When galantamine was combined with mifepristone, caffeine, or diclofenac, it caused a reduction in oxidative stress responses by 5.8-6.3%. These neuroprotective effects were more evident in the neurons dosed with Aβ, observing a reduction of ROS produced of 17.7% with gal+mifep, 13.5% with gal+caffeine, and 2% with gal+diclofenac, whereas single galantamine treatment presented a 5.5% increase in neuronal oxidative stress (Fig. 3). The neuroprotective effect of mifepristone is related to preventing stress-induced responses, protecting neurons from undergoing programmed cell death (apoptosis) and increasing AMPA receptor expression,(36) also mifepristone may promote rapid repair of the synaptic alterations in the hippocampus. (37) In contrast, caffeine activates noradrenaline neurons and the releasing of dopamine, and counteracting the development of cognitive impairment as an antioxidant. (38) Diclofenac is an antiinflammatory drug targeting voltage-dependent K + channels, with neuroprotective effects in neurons. (39) In pilot studies, diclofenac showed positive effects to reduce the risk of AD in veterans. (40) Rivastigmine is an AChE and butyrylcholinesterase inhibitor, used to treat memory loss and to improve thinking abilities in AD patients. The tolerability and safety of rivastigmine or in combination with lithium, vitamin E, and diclofenac was tested, observing that these drugs were well-tolerated by cholinergic neurons.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Neuronal responses indicated that these active compounds effectively reduced oxidative stress in neurons by 4.8-6.2%. In the presence of Aβ, the rivastigmine, riv+lithium, and riv+vit E were able to regulate ROS with similar levels for all these treatments, whereas riv+diclof reduced by 6.8% the amount of intracellular ROS (Fig. 3). At the neuronal level, lithium can reduce excitatory neurotransmitters (dopamine and glutamate) and increase GABA.

DISCUSSION
The objective of this study was to provide proof-of-concept hypotheses on drug combinations to treat AD by knowledge graph embedding of heterogeneous interactome and multi-task learning of complementary evidence on drug efficacy. Our AD knowledge graph embedding boosted by pretrained universal pharmacological embedding integrates not only the human interactomes but also millions of interactions in public biomedical literature. Our multi-task learning of AD-related genes and multi-level drug efficacy closes the gap between the scattered knowledge in AD research.
Our methodology tackles important challenges to harmonize fragmented multimodal data. Recently, several prior works have used the graph neural network approach for drug development, either for general purpose, (44,45) for rare diseases, (46) or for specific diseases including COVID-19 (47)(48)(49) and AD. (50) Our methodological strength over the prior approach is i) transfer learning of universal embedding to integrate broader modality to address the incomplete knowledge, ii) multi-task learning for multimodal interactions, AD-related genes, and multi-level drug efficacy to mitigate the lack of established therapeutic targets, and iii) in vitro validation of the computational model via hippocampal neuron cells.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Drug combinations that we validated are less toxic (FDA-approved drugs) and economically viable (generic drugs). AD is multifactorial and the use of multi-target therapeutic interventions addressing several molecular targets seems to be an alternative approach to modify the course of AD progression.(51) Different clinical studies demonstrate that a combination therapy is more efficient than monotherapy, especially in mild-to-moderate AD to slow down the rate of cognitive impairment. Our integrative approach has identified promising drug combinations adjuvant to existing FDA-approved AD drugs. All drug combinations we identified are available in generic formulations, making them economically viable and available as an alternative to patented formulations with a high monthly cost (>500 USD).
Our study has several limitations. Instead of a mechanism-based approach that identifies therapeutic targets first and identifies drugs interacting with them, our approach relies on historical or data-driven evidence, such as estimated drug efficacy (transcriptomic reversed patterns, population-based treatment effect) and existing efficacy (preclinical or clinical trials). If this data-driven drug efficacy is biased, our model's drug ranking can be also biased. For example, drugs targeting amyloid beta are much more studied than drugs targeting other pathologies, thus our model has a potential risk to trap in the circularity of failed drugs targeting amyloid beta. Population-based treatment effects in real-world data can partly mitigate the bias in historical research data toward amyloid-beta. Another limitation is that the identified drug combinations are based on neither pharmacokinetics/pharmacodynamics nor biological experimental data on synergistic drug combinations in treating AD. Instead, our approach was to use the complementary exposure pattern, which oversimplifies the synergistic mechanism of actions.
Despite the limitations, our drug repurposing strategy has a broad impact not only on AD but also on complex diseases, in which complex mechanisms of actions and therapeutic targets are an open question.
The well-known deficiency of knowledge and insufficiency of data in AD research, when studied separately, prevent us from capturing the complex AD's heterogeneous pathogenesis. Our pipeline integrates complementary information at multiple stages (such as transcriptomic patterns, preclinical efficacy, clinical partial effectiveness, and population-based treatment effect) by jointly optimization to . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. better identify promising drug combinations. Hypotheses generated from this integrative approach may shed light on complex disease drug development.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Build AD knowledge graph
We built the AD knowledge graph from various information sources. We queried "Alzheimer's Disease" and downloaded AD-related multimodal interactions from the Comparative Toxicogenomics Database is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint drug-AD interactions, or a gene-GO annotation in AD, or both. These expanded our AD knowledge graph with 6,817 drug-GO interactions, including 2,377 unique GO terms and 108 unique drugs.

Knowledge graph representation
Graph neural network (GNN) is one field of deep neural networks that derive a vectorized representation of nodes, edges, or whole graphs. Adopting GNN into the biomedical network facilitates the integration of multimodal and complex relationships. The graph node embedding can preserve the node's local role and global position in the graph via iterative and nonlinear message passing and aggregation. It learns the structural properties of the neighborhood and the graph's overall topological structure (64), thus allowing us not only to restore the known interactions but also infer unknown interactions. The knowledge graph representation that contains such putative interactions can help us find new indications of drugs. (65) Recently GNN has demonstrated a great advance in predicting hidden interactions (e.g., PPIs, drug-drug adverse interactions, and drug-target interactions) and the discovery of new molecules (45,66). GNN is also used to derive representation from a graph with multiple types of relations and nodes (i.e., In this study, we developed the knowledge graph representation model by customizing the deep variational graph neural autoencoders (VGAE) to incorporate multiple types of relations or edges (Fig. 1b). (17) The self-supervised graph autoencoder method encodes the nodes into a latent vector (embedding) and reconstructs the given graph structure (i.e., graph adjacency matrix) with the encoded latent vector. The variational graph approach can encapsulate the complex network into a probabilistic distribution (rather . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint than deterministic vector representation) for nodes considering the uncertainties of our knowledge for AD; therefore, it alleviates the overfitting/underfitting issues due to partial knowledge. We customized the VGAE to incorporate the different types of edge by setting different weight matrices for each type of edge.
This multi-relational approach respects the different propagation for each interaction type. Since our objective is to derive node embedding that reflects the overall graph topology, the multi-relational VGAE model was trained to reconstruct the missing interaction using the node embeddings as an autoencoding manner. For the optimization, in the message-passing step, each node (entity)'s embedding is iteratively updated by aggregating the neighbors embedding, in which the aggregation function is a mean of the neighbor's features, concatenation with current embedding, and a single layer of a neural network on the concatenated one. We set the embedding size as 128 after several trials. We used PyTorch Geometric for implementation (72). The model structure was (1×400) → Graph convolution to (1×256) → RELU → Dropout → Summation of multiple edge types → Batch norm → Graph convolution to 1×128 (mean) and 1×128 (variance). We randomly split the knowledge graph edge (i.e., positive edges) into 90% for training and 10 % for testing. We randomly generated the same number of fake or negative edges from which the positive edges will be predicted.

Transfer learning from a large drug knowledge graph
We boosted the node embedding's biological relevance in the AD knowledge graph with millions of interactions from public biomedical literature by transfer learning. Transfer learning is to transfer knowledge from previously learned universal models to domain-specific models. By transferring and finetuning the universal knowledge, the AD network embedding can respect universal pharmacological interactions while prioritizing AD-related interactions. Injecting this universal knowledge into the AD domain is critically important in addressing the uncertainty in AD knowledge. We used the DRKG's pretrained node embedding as initial node features in the AD knowledge graph. A recent attempt (e.g., DRKG network(49)) has been made to obtain the universal embedding of 15 million pharmacological entities considering 39 different types of interactions among genes and compounds from seven large . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint biomedical databases including DrugBank, HetioNet, GNBR, STRING, IntAct, and DGLdb. (44) We retrieved the pretrained embedding for drugs, genes, pathways, and functional ontology from DRKG. We set the initial node feature of the AD knowledge graph as the pre-trained embedding of DRKG and finetuned them to the AD knowledge graph. When nodes in the knowledge graph were not matched to entities in the universal embedding, we learned the embedding as parameters (e.g., one hot to embedding vectors).

Multi-task learning to distinguish high-confidence AD-related genes from other genes
One of the most important tasks in AD drug development is identifying driver genes related to AD. We encouraged the knowledge graph representation on genes to discriminate the AD-related genes out of the 16,997 genes in the graph. We obtained 743 high-confidence AD-related genes from Agora(54) and CTD.
We hypothesized that one can predict AD-related genes using gene interaction with other genes, GO, and pathways. We created a node classification task to predict whether a gene is one of the AD-related genes using the gene representation in the AD knowledge graph. The node classification consists of two graph convolution layers (GraphSAGE) that sample a node's adjacent neighborhoods and aggregate their representation. The dimensions of those two convolution layers were Graph convolution (128, 64) → RELU → Dropout → Graph convolution (64,1). We split the gene nodes into 90% for training and 10% for testing.
As the labels are not balanced (only 743 AD-related genes out of 16,997 genes), we randomly split the gene nodes into 90% for training and 10% for testing while maintaining the ratio of AD gene and non-AD gene.
This AD-related gene classification task was jointly optimized with the relational VGAE so that the AD knowledge graph representation can distinguish AD-related genes from all other genes. We used a multitask optimization technique that resolves conflict among multiple tasks' gradients by normal vector projections, (73) to achieve a better local optimum.

Validating the quality of knowledge graph representation
We evaluated the quality of AD knowledge graph representation by checking whether the learned node embedding can restore the known edges. The edge prediction task is to predict whether an edge exists . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint between two nodes using the two nodes' embedding. The edge prediction task's accuracy will be high if the node embedding encapsulates the edge information faithfully. For training and test, we randomly set aside 10% of the edges for evaluation. We randomly sampled negative or fake edges that do not exist in the knowledge graph (within the same node types) and predicted whether a given edge is real or fake. (17) We measured the edge prediction accuracy by each edge type, to check whether any edge types were neglected or underrepresented in the node embedding. The model was trained using all edge types and evaluated separately for each edge type. We also ablated one edge type and observed how one edge type contributed to predicting all the edge types (Table S1). This experiment allowed us to determine which interactomes are unique and hard to be replaced or implicitly inferred by other edge types. The edge prediction accuracy was measured as the area under the receiver operating curve (AUROC) and the area under the precisionrecall curve (AUPRC). For the node classification task to distinguish between AD-related genes and remaining genes, we also measured the classification accuracy using the AUROC and AUPRC.

Transcriptomic, mechanistic, and epidemiological evidence for drug efficacy
We built a pipeline to prioritize drugs that are similar to drugs with historical evidence. Particularly, we focused on drugs with multiple historical evidence. AD is a complex disease and its effective therapeutic targets are largely unknown, thus drugs interacting with a few targets can have limited efficacy.(74) A holistic evidence may alternatively capture the efficacy in the systemic and complex disease. We believe the most desirable drug candidates would have multiple levels of complementary evidence to support their efficacy, such as post-treatment transcriptomic patterns, mechanistic efficacy in preclinical animal models, and epidemiological treatment effects in real-world patient data. We identified each evidence as follow: Transcriptomic reversed patterns: We examined reversed patterns of the AD model's genetic signature and drug-induced gene signature using gene set enrichment analysis. We leveraged Gareth et al. 's iPSCderived cortical neurons studies (9) and identified 149 drugs with significant transcriptomic reversed patterns. Among the 149 drugs, 76 drugs were matched to our 6,543 drug candidates.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Mechanistic efficacy in preclinical models: Alzheimer's Disease Preclinical Efficacy Database (AlzPED) accumulated pre-clinical evidence for effective AD drugs. We identified 338 drugs with significant efficacy from 791 animal-based studies targeting various pathologies including amyloid-beta, tau, and neuroinflammation. Among the 338 drugs, 189 drugs were matched to our 6,543 drug candidates.
Population-based treatment effect: We estimated drugs' population-based treatment effect as real-world evidence of drug efficacy. Healthcare claim data is a unique resource to investigate treatment effects of prescribed drugs on reducing AD onset risk during the preclinical stage of AD. We used Optum's deidentified Clinformatics ® Data Mart subscribed by UTHealth. This billing-purpose claim data is inherently noisy and sparse, thus requires careful data preparation (e.g., setting observation windows to exclude subjects who are not old enough to have AD onset, identifying AD onset via billing-purpose diagnosis and medication codes, grouping high-dimensional diagnosis codes into clinical comorbidities as confounding variables). We followed the data preparation process in our previous work(75) based on target trials. (76) We calculated the average treatment effect among treated (ATT) using inverse probability of treatment weighting (IPTW). (77,78) See details at Supplementary Method 1. We identified 126 drugs with positive treatment effects (i.e., ATT<0 to reduce AD onset). Among the 126 drugs, 101 drugs were matched to our 6,543 drug candidates.
Clinical trials failed at Phase II/III: Many drugs in the failed clinical trials passed several phases and showed non-trivial response rates (e.g., good response on certain subpopulations). We hypothesize individual drugs (focusing on different targets) from previously failed clinical trials respond differently on sub-populations, and therefore, cocktails of complementary drugs might cover a larger population. Also, because drug testing in clinical trials is partially efficacious in treating AD, a drug that is similar to these trial drugs (and their combinations) can have potential efficacy too. We collected 294 drugs from 2,675 interventional clinical trials for AD from ClinicalTrial.gov. Among the 294 drugs, 210 drugs were matched to our 6,543 drug candidates.

Prioritizing repurposable drugs with multiple evidence and multi-task learning
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint We extracted four positive labels from extensive calculation and search (Method 7). We applied the multitask learning (73) to jointly build a drug ranking model that ranks drugs high if the drugs satisfy multiple evidence. We designed a ranking model with neural networks and Bayesian pairwise ranking loss. (79,80) The positive samples to prioritize were the drugs with positive evidence for each task; the unlabeled or unknown samples were the remaining drugs in the AD knowledge graph. The architecture was two fully connected layers (with the size of 128→128→1) with residual connection, nonlinear activation (ReLU), dropout, batch norm in the middle, and the optimization loss (Bayesian pairwise ranking loss). The model was trained to minimize four loss functions for each task. We validated the ranking model by evaluating whether the positive drugs (under clinical trials) are ranked high. We measured the accuracy of the drug ranking model using the area under the receiver operating curve (AUROC) and area under the precisionrecall curve (AUPRC) of top k drugs with 50% training and 50% test cross-validation. We purposely set the portion of the training set lower because the clinical trials are not our sole "gold standard" to prioritize drugs.

Drug combination search
We identified efficacious drug combinations from top-ranked drugs. Our approach is to leverage the network of high-confidence drugs' targets and AD risk genes. Our hypothesis was that "a drug combination is therapeutically effective only if the targets of the drugs both hit the disease module, but they target a separate neighborhood (Complementary Exposure pattern)",(33) which has been successfully applied to hypertension (33) and COVID-19. (48) We searched the drug combinations within the top-ranked drugs. We counted the number of genes in the AD module that a drug combination hits, where the drug combination's targets are disjoint. Imaging. Neurons were analyzed in fluorescence, bright-field, phase contrast, live-cell, and Z-stack imaging with Digital Imaging System Celena-S (Logos). Micrographs were analyzed with Fiji-Image J.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Statistical analysis. Data acquired from bioassays were normalized to control for each data set, the normalized average is presented as mean ± standard error of the mean (SEM). SigmaPlot was used for the analysis.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021.

Data Availability
All study data are publicly available from CTD, STRING, Agora, ClinicalTrial.gov, and AlzPED.
Experimental results are available at supplementary materials. The code is available at https://github.com/freshnemo/AD-KG. Table S1. Ablation study on edge prediction. Table S2. Node classification task to predict therapeutic targets out of all gene nodes. Table S3. Drugs with statistically significant treatment effects in reducing AD onset in Optum claim data.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint text. Four main AD risk genes (green triangle) and undergoing clinical trial drug (red circle) clusters were highlighted in subplots. Genes, which are in gray triangles, were mixed with drugs, which are in the black rounds. Gene ontologies (pink diamond) and pathways (light grey circle) are adjacent to genes and drugs.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint Table 1. Edge prediction accuracy using AD knowledge graph's node embedding. To examine which edge type is difficult/easy to contain in the node embedding, we measured the edge prediction accuracy for each edge type separately. The accuracy values for all edge types were consistent without any edge types being under-fitted. For AD knowledge graph representation without transfer learning, we initialized the node embedding with random initialization. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted December 5, 2021. ; https://doi.org/10.1101/2021.12.03.21267235 doi: medRxiv preprint