Next Article in Journal
Combination Strategies for Immune-Checkpoint Blockade and Response Prediction by Artificial Intelligence
Next Article in Special Issue
Roles of Autophagy in Oxidative Stress
Previous Article in Journal
The Effects of a Weight-Loss Herbal Formula RCM-107 and Its Eight Individual Ingredients on Glucagon-Like Peptide-1 Secretion—An In Vitro and In Silico Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis

1
Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15213, USA
2
School of Medicine, Tsinghua University, Beijing 100084, China
3
Department of Pediatrics, School of Medicine, Washington University in St. Louis, St. Louis, MO 63130, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(8), 2855; https://doi.org/10.3390/ijms21082855
Submission received: 25 March 2020 / Revised: 14 April 2020 / Accepted: 16 April 2020 / Published: 19 April 2020
(This article belongs to the Special Issue Autophagy in Health, Ageing and Disease 2.0)

Abstract

:
Autophagy plays an essential role in cell survival/death and functioning. Modulation of autophagy has been recognized as a promising therapeutic strategy against diseases/disorders associated with uncontrolled growth or accumulation of biomolecular aggregates, organelles, or cells including those caused by cancer, aging, neurodegeneration, and liver diseases such as α1-antitrypsin deficiency. Numerous pharmacological agents that enhance or suppress autophagy have been discovered. However, their molecular mechanisms of action are far from clear. Here, we collected a set of 225 autophagy modulators and carried out a comprehensive quantitative systems pharmacology (QSP) analysis of their targets using both existing databases and predictions made by our machine learning algorithm. Autophagy modulators include several highly promiscuous drugs (e.g., artenimol and olanzapine acting as activators, fostamatinib as an inhibitor, or melatonin as a dual-modulator) as well as selected drugs that uniquely target specific proteins (~30% of modulators). They are mediated by three layers of regulation: (i) pathways involving core autophagy-related (ATG) proteins such as mTOR, AKT, and AMPK; (ii) upstream signaling events that regulate the activity of ATG pathways such as calcium-, cAMP-, and MAPK-signaling pathways; and (iii) transcription factors regulating the expression of ATG proteins such as TFEB, TFE3, HIF-1, FoxO, and NF-κB. Our results suggest that PKA serves as a linker, bridging various signal transduction events and autophagy. These new insights contribute to a better assessment of the mechanism of action of autophagy modulators as well as their side effects, development of novel polypharmacological strategies, and identification of drug repurposing opportunities.

1. Introduction

Autophagy is the process that enables self-clearance of potentially damaging protein aggregates, dysfunctional organelles (e.g., mitochondria by mitophagy), or entire cells, and as such it is one of the most significant regulatory processes that ensures healthy functioning of cells [1,2]. During starvation, autophagy can serve as an essential source of energy for survival, and in general it plays a housekeeping role to maintain homeostasis and viability of cells [3]. In addition, autophagy has been shown to be instrumental in erasing pathological responses or maintaining physiological functions of the cells such as anti-aging, tumor suppression, and antigen presentation [4].
Once autophagy is activated, e.g., by stress stimuli such as starvation, a cascade of reactions are triggered that lead to the formation of a double-membrane-bound vesicle called an autophagosome [5], which encapsulates the cytoplasmic aggregates or dysfunctional organelles. The subsequent fusion of the autophagosome with the lysosome to form the autolysosome permits the lysosomal enzymes to digest and degrade the contents of the autophagosome [5,6].
According to our current understanding, autophagic responses are modulated by many events, in which the kinase mammalian target of rapamycin (mTOR) plays a major role as a master regulator of autophagy [7]. mTOR, a serine/threonine kinase, is a downstream regulator of the phosphoinositide 3-kinase (PI3K) pathway. It suppresses autophagy induction by inhibiting the activation of the UNC-51-like autophagy activating kinase (ULK) complex [8]. The complex of activated ULK1 with autophagy related 1 (Atg1) protein, Atg1-ULK, would otherwise receive signals on cellular nutrient status and recruit downstream autophagy-related (ATG) proteins to the autophagosome formation site, thus initiating an autophagic response [9]. mTOR further acts as a negative regulator of autophagy by preventing the general expression of lysosomal and autophagy genes through phosphorylating transcription factor EB (TFEB) [7,8].
Other important regulators of autophagy include cAMP-dependent protein kinase A (PKA) [10] and other proteins involved in Ca2+ signaling pathways. Ca2+ ions have complex effects [11,12]: Cytoplasmic Ca2+ activates Ca2+/calmodulin-dependent kinase kinase β (CaMKKβ), which activates AMP-activated protein kinase (AMPK) to promote autophagy [13]. However, CaMKKβ also activates calpain, which inhibits autophagy by cleaving autophagy related 5 (Atg5) protein [14].
In the last few decades, there has been a surge in the number of autophagy-modulating drugs targeting these well-known autophagy regulators. For example, rapamycin that inhibits mTOR is known as one of the most important autophagy activators [15], and wortmannin which inhibits PI3K is known as a significant autophagy inhibitor [16,17]. In the meantime, autophagy modulation has been receiving increased attention as a strategy for developing therapeutics [18]. Autophagy has been widely accepted as cytoprotective against neurodegenerative diseases, and many clinical interventions are moving forward to increase autophagy as a therapeutic intervention [19]. For example, rapamycin is regarded as an important pharmacological agent against Alzheimer’s disease (AD) [20]. Various autophagy-enhancing drugs (such as carbamazepine [21], fluphenazine [22], and glibenclamide [23]) have been identified to have therapeutic effects against α1 antitrypsin deficiency (ATD) caused by the aggregation of the misfolded α1-antitrypsin Z (ATZ) proteins in the liver. Besides, autophagy dysfunction is also highly associated with cancer, and many autophagy-targeting drugs also serve as anticancer drugs [24]. For instance, metformin, an autophagy enhancer, is a drug against myeloma [1].
In this study, we carried out a comprehensive analysis of a set of 225 known autophagy modulators selected from DrugBank v5.1.1 [25] (Table S1), their 993 target proteins, and 294 associated pathways in order to obtain a systems-level understanding of the common mechanisms of autophagy modulation, as well as specific mechanisms of selected modulators. The set of modulators consists of 174 activators, 31 inhibitors, and 20 compounds acting as dual-modulators. We identified 1831 known interactions between these drugs and their targets, and predicted 368 novel interactions using our probabilistic matrix factorization (PMF)-based application programming interface [26,27,28]. Of these predictions, 75 were consistent with recently published experimental data (not yet deposited in DrugBank [25]).
Our results show that autophagy modulators are structurally diverse and functionally promiscuous. Their target proteins are commonly expressed in the brain and liver tissues and mainly function as tyrosine kinases, calcium channels, aminergic GPCRs, or nuclear hormone receptors. The most frequently targeted proteins do not directly take part in the autophagy pathway, but indirectly regulate the activity/expression level of core ATG proteins through differential signaling pathways. Through gene ontology (GO) annotation and pathway enrichment analysis, we distinguished the most important signaling pathways that regulate autophagy. Our results strongly point to the crucial role of PKA as a hub conveying autophagy signals, and directly regulating ATG proteins and their gene expression.

2. Results

2.1. Autophagy Modulators Are Structurally Diverse and Have Promiscuous Functions

We manually curated the role of the selected 225 autophagy modulators based on data/reports in the literature, and we classified them into three functional categories: activators (n = 174), inhibitors (n = 31), and dual-modulators (n = 20). Here, dual-modulators refer to drugs that can both positively and negatively regulate autophagy depending on the biological context. For instance, melatonin is a dual-modulator because it reduces autophagic activity in tumor trophoblast cells while it enhances it in normal cells [29].
As a first step, we examined to what extent the autophagy modulators within each group share structural and functional similarities. Figure S1 shows that the selected autophagy modulators are highly diverse. Within all three categories, most of the drug pairs have low structural and functional similarities (<0.5) evaluated based on their 2D fingerprints and drug–target interaction patterns extracted from DrugBank, respectively. However, we can distinguish clusters of drugs that have similar structures, which also share similar interaction patterns with targets. For example, rapamycin and its derivatives everolimus, temsirolimus, and ridaforolimus (Figure S1, panels a and b, activators #28-30, enclosed in light yellow ellipse) activate autophagy by inhibiting mTOR [30], and retain the same structure, presumably targeting the same site on mTOR. In contrast, the upper right white region in panel b (drugs #141-166, enclosed in pink ellipse) indicates several autophagy activators (e.g., fluphenazine, pimozide, clonidine, paroxetine, triflupromazine, chlorpromazine, citalopram, nortriptyline, fluspirilene, doxazosin, amiodarone, flunarizine, verapamil, and dronedarone) that share similar interaction patterns. Most of them are therapeutic agents for mental disorders such as schizophrenia, depression, and anxiety disorders, or for cardiovascular diseases such as high blood pressure, angina, and arrhythmia, regulating autophagy through cAMP- and Ca2+-signaling pathways. However, panel a indicates that they are structurally heterogeneous, suggesting that they have different targets on their shared pathways.
Further examination of Figure S1 shows that the macrolide antibiotics erythromycin, clarithromycin, and azithromycin (inhibitors #5-7 in panels c and d) inhibit autophagy via endoplasmic reticulum (ER) stress-mediated C/EBP homologous protein (CHOP) induction [31], but despite their close structural similarity (see enclosed region in panel c), they exhibit distinctive interaction patterns with their targets (panel d), suggesting different mechanisms of action. The phenothiazine antipsychotics thioridazine and trifluoperazine (panels e and f, dual-modulators #2-3) regulate autophagy through dopamine receptors [32] and share close structural similarities.

2.2. Selected Autophagy Modulators Are Distinguished by Their High Promiscuity

The space of proteins targeted by autophagy modulators is quite broad. We identified 993 such proteins (Table S2) composed of 706 targets associated with 174 activators, 374 associated with the 31 inhibitors, and 94 targets associated with the 20 dual-modulators, which show partial overlaps, as presented by the Venn diagram in Figure 1a (right). The space of drugs and targets may be viewed as a bipartite network, with multiple connections (drug–target interactions). The number of targets connected to a given modulator will be referred to as the degree of the modulator. The higher the degree, the more promiscuous the modulator.
Based on DrugBank v5.1.1, we identified a total of 1831 types of known drug–target interactions between these autophagy modulators and their target proteins. Specifically, the activators, inhibitors, and dual-modulators were involved in 1339, 389, and 103 known interactions with their targets (Figure 1a).
In principle, if each autophagy modulator had one distinctive target, we would have 225 targets (or even less as some targets could overlap), but this is not the case: Several modulators are highly promiscuous (Figure S2) in the sense that they interact with multiple targets. The detailed distribution of the multiplicity of targets among modulators is summarized in Figure 1b. Only 25%–32% of the modulators interact with a single drug (red portions of the bars), depending on the type of modulator. It is interesting to note that almost 20% of activators have more than 10 targets (yellow portion at the top of the leftmost bar), pointing to the high promiscuity of ~30 activators. Specifically, 31 out of 174 (17.8%) activators, 3 out of 31 (9.7%) inhibitors, and 3 out of 20 (15.0%) dual-modulators have at least 10 known targets.
Figure 2 provides more details on the identity of specific modulators with high promiscuity. The ordinates show the number of targets for autophagy activators (panel a), and all autophagy inhibitors (panel b) and dual-modulators (panel c), and the abscissa shows the modulators reordered in decreasing promiscuity. See also Figure S3 for the target number distribution of additional (less promiscuous) activators. We observed an average of 7.7 interactions per activator, 2.9 per inhibitor (if we exclude one outlier, fostamatinib; see below), and 5.3 per dual-modulator. The insets in the three panels of Figure 2 show the histograms of modulators with different promiscuity.
Notably, of the 1339 interactions in which activators are involved, 145 are made by the most promiscuous activator, copper (Figure 2a, leftmost bar), followed by zinc (124 interactions/targets), artenimol (78 targets), and olanzapine (48 targets). Transition metals copper and zinc are non-structural intracellular signaling mediators and essential elements in many enzymes. They both can induce autophagy by activating kinases such as mitogen-activated protein kinase (MAPK) [33,34]. A recent study shows that copper can even directly bind to ULK1/2 to enhance its pro-autophagic activity [35]. The antimalarial drug artenimol is a derivative of artemisinin that effectively kills P. falciparum by producing reactive oxygen [36] and targeting multiple biological processes otherwise essential for parasite survival [37]. The oxidative stresses induced by artemisinin and its derivatives also stimulate autophagy in the host cells [38,39]. As an antipsychotic drug, olanzapine mainly targets receptors of neurotransmitters such as dopamine, serotonin (5-HT), and γ-aminobutyric acid (GABA), which have been reported to regulate autophagy [32,40,41]. Further, olanzapine also induces autophagy by activating AMP-activated protein kinase (AMPK) signaling [42,43] and upregulating the expression of the products of ATG genes such as Atg5 and Beclin-1 [44].
The most promiscuous autophagy inhibitor is fostamatinib (303 targets). Fostamatinib is an FDA-approved drug for the treatment of chronic immune thrombocytopenia. It has been reported that fostamatinib suppresses autophagy by inhibiting the spleen tyrosine kinase (STK) [45]. Fostamatinib also inhibits critical autophagy regulators including ULK1/2, serine/threonine-protein kinase (TBK1) [46], death-associated protein kinase 1 and 2 (DAPK1/2) [47,48], and leucine-rich repeat serine/threonine-protein kinase 2 (LRRK2) [49].
The most promiscuous dual-modulator of autophagy is calcium (20 targets), followed by melatonin (10 targets). Ca2+ is a cofactor of many enzymes and also acts as a second messenger mediating various signal transduction pathways. It targets calcium channels (e.g., Ca2+ voltage-dependent channel subunit α1c, CACNA1C), calcium-transporting ATPase (e.g., ATP2C1), and downstream calcium-binding proteins (e.g., calmodulin). Cytoplasmic Ca2+ inhibits autophagy via activation of calpain and the inositol trisphosphate receptor (IP3R), while it also enhances autophagy via activation of AMPK [11]. Our previous study suggested that Ca2+/calmodulin-dependent kinase kinase β (CaMKKβ) plays a key role in regulating the balance between these opposing actions and acts as a determinant of the effect of competing roles of cytoplasmic Ca2+ in autophagy regulation [11].
To further identify potential drug-target interactions in DrugBank, we used our PMF model (see Materials and Methods) for predicting new targets for autophagy modulators. This analysis led to 368 novel interactions with high confidence scores (>0.6), which included 12 additional targets. These pairs (not reported in the DrugBank) are listed in the Table S3. The number of predicted targets for each drug are shown in Figure 2 by the light-shaded portions of the bars. Notably, 75 (out of 368) predicted pairs are consistent with already published experimental data [50,51,52]. Table 1 lists 20 of them, along with the corresponding experimentally measured binding affinities, and others can be seen from the entries in Table S3 where experimental binding affinities are listed. For example, the promiscuous drugs olanzapine and fostamatinib mentioned earlier were found to bind to serotonin receptor 1F (HTR1F) and fibroblast growth factor receptor 4 (FGFR4), respectively. Fluphenazine, fluspirilene, thioridazine, sertindole, and trifluoperazine were identified to bind to dopamine receptor D3 (DRD3), and experiments indicate their high affinities (Ki < 5 nM). We also detected various drug–target interactions with modest-to-weak affinities (e.g., the interaction between imatinib and Fms related receptor tyrosine kinase 3 (FLT3)).
Taken together, the 225 autophagy modulators we selected are diverse in terms of both structure and function. Several autophagy modulators, and especially autophagy activators, are distinguished by their highly promiscuous effects. The additional drug–target interactions we identified using the PMF model are plausible, thus complementing existing knowledge and also helping generate new hypotheses that await further testing. To have a better understanding of the effects of autophagy modulators, we next analyzed their target proteins.

2.3. Frequent Targets of Autophagy Modulators Are Not ATG Proteins But Their Regulators

As described above, the set of autophagy modulators have 1005 targets (993 known/reported in DrugBank plus 12 identified here). While drugs were highly promiscuous, the targets usually had a smaller number of connections to drugs in the bipartite network, representing drug–target pairs. We observed an average of 1.9 interactions per target, and 66.9% of the targets interacted with one single drug.
We show in Figure 3a the most frequently targeted proteins (51 of them), which interact with five or more known autophagy modulators. The bars display the corresponding numbers of modulators broken down by color-coded categories, with dark and light shaded portions referring to known and predicted interactions, respectively. We refer to this set of proteins as frequent targets (FTs) of autophagy modulators. The distribution of autophagy targets over protein families, extracted from family memberships given in Uniprot [53] is presented in Figure 3b. Kinases (32.17%), ion channels (6.76%, including ligand-gated ion channels, Ca2+ channels, and Na+ channels), and G protein coupled receptors (GPCRs) (5.85%) are three major protein families targeted by our autophagy modulators. However, this presumably reflects the historical dominance of these families among drug targets, noted in earlier studies [54]. More meaningful are the identities of the specific targets belonging to these families, which are analyzed next in relation to ATG proteins.
In order to see whether FTs overlap with the targets involved in autophagy pathways, we compared them to the 128 proteins participating in the canonical autophagy pathway (KEGG id: hsa04140; in the following, we refer to these 128 proteins as ATG proteins). Thirty-two out of 1005 targets of autophagy modulators were found to be ATG proteins, including mTOR, ULK1/2, AMPK, CAMKKβ, insulin (INS), AKT, PI3K, PKA, PKC, phosphoinositide-dependent kinase 1 (PDK1), Raf1, ERK1/2, JNK2/3, Beclin-1, BCL2, DAPK1/2/3, cathepsin L (CTSL), inositol-requiring enzyme 1 (IRE1), eukaryotic translation inhibition factor 2α kinase 4 (EIF2AK4), and hypoxia inducible factor 1α (HIF-1α) (see the rows highlighted in green in Table S2). Among them, 15 are targeted by autophagy activators, 26 by inhibitors, and two by dual-modulators. Among the ATG proteins, only mTOR is an FT. This analysis therefore showed that the majority of ATG proteins are not FTs of autophagy modulators.
The question was then to assess what the FTs are, and why they are the targets of many autophagy modulators. We next conducted a more detailed analysis of the FTs.
Figure 3a shows that a nuclear hormone receptor (NHR) (subfamily 1 group I member 2 (NR1I2) is the most frequently targeted protein. NR1I2 is found here to interact with 14 activators, one inhibitor, and one dual-modulator. The NR1I2 gene encodes the pregnane X receptor that regulates the transcription of many genes involved in the metabolism and secretion of potentially harmful endogenous and xenobiotic compounds, and in the homeostasis of glucose, lipid, cholesterol, bile acid, and bilirubin. The current analysis strongly suggests that it also activates the genes involved in autophagic elimination. There are six more NHRs that are FTs: estrogen receptor 1 and 2 (ESR1, ESR2), androgen receptor (ANDR), and peroxisome proliferator activated receptors γ, α, and δ (PPARG, PPARA, PPARD). Numerous studies have reported that these transcription factors regulate the expression of core ATG proteins. Specifically, ESR1 and ESR2 regulate mTOR expression (see review [55]). ANDR regulates the transcription of ULK1/2, Atg4B/D, and transcription factor EB (TFEB), a master regulator of autophagy genes [56]. PPARs are known to regulate autophagy via AMPK, phosphatase, and tensin homolog (PTEN) expression [57].
Eleven FTs directly participate in calcium signaling, including calmodulin and 10 ion channel family members: Ca2+ voltage-gated channel subunit α1c (CACNA1C) and other members of the same family (CACNA1D, CACNB2, CACNA1S, CACNA1H, CACNA2D1, CACNB1, CACNA2D2, K+ voltage-gated channel subfamily H member 2 (KCNH2), and Na+ voltage-gated channel α subunit 7 (SCN6A). We note that calmodulin interacts with 9 activators and 3 dual-modulators, consistent with the complex role of cytoplasmic Ca2+ in regulating autophagy.
Eighteen FTs are GPCRs, consisting of α1A adrenergic receptor ADRA1A and its family members (ADRA1B, ADRA1D, ADRB1, ADRB2, ADRA2A, ADRA2B, ADRA2C), dopamine receptors (DRD1, DRD2), histamine receptor HRH1, serotonin receptors (HTR2A, HTR2C), and muscarinic receptor CHRM1-5. GPCRs regulate autophagy through diverse downstream signaling pathways [58]. For example, ADRA1 activation induces protein kinase B (AKT)-mediated autophagy. β-adrenergic receptors (ADRB1 and ADRB2) promote autophagy via the cAMP pathway [59]. CHRM3 regulates autophagy via the CaMKKβ-AMPK pathway [60]. DRD1 stimulation induces autophagy via a cAMP-dependent but EPAC-independent mechanism [61], while the D2-like family of dopamine receptors (DRD2-4) positively regulates autophagy through AKT-mTOR and AMPK pathways [51,62]. Furthermore, activation of HTR2 inhibits autophagy via mTOR-independent p70S6K and 4E-EP1 phosphorylation [63].
As shown in Figure 3b, 32.17% of the targets are protein kinases mediating intracellular signal transduction. Among them, seven are FTs, including mTOR and several receptor tyrosine kinases (RTKs): epidermal growth factor receptor (EGFR), vascular endothelial growth factor receptor 2 (KDR), Abelson murine leukemia viral oncogene homolog 1 (ABL1), proto-oncogene c-KIT (KIT), and platelet-derived growth factor receptors A and B (PDGFRA and PDGFRB). Accumulating evidence suggests an interplay between RTK signaling and autophagy [64]. In particular, EGFR activation has been shown to inhibit autophagy via Beclin-1 phosphorylation [65]; KDR activation suppresses autophagy by upregulating FoxO1 expression; and ABL1 regulates early and late stages of autophagy by phosphorylating Beclin-1 and promoting lysosomal trafficking, respectively [66]. Furthermore, mutant KIT in acute myeloid leukemia (AML) cells triggers autophagy via STAT3 activation. Inhibition of PDGFRs can induce autophagy through the PI3K-AKT pathway [67].
This detailed analysis points to several targets whose pro- or anti-autophagic roles are consistent with experimental data. It also shows that the ATG proteins, except for mTOR, are not necessarily the direct targets of autophagy modulators, and a broad number of activators, inhibitors, and dual-modulators instead target proteins involved in the transcription regulation of genes involved in autophagy, or in intracellular signal transduction mainly via protein kinases and Ca2+ and Na+ channels. Notably, six targets interact with all three categories of autophagy modulators (see the Venn diagram in Figure 1 (right)): two nuclear hormone receptors (NR1I2 and ESR1), two proteins involved in potassium channeling and transport (KCNH2 and ATP4A, respectively), a GPCR (ADORA2A), and the casein kinase II subunit α (CSNK2A1, also known as CK2α) (highlighted in yellow in Table S2).

2.4. Autophagy Modulation Targets Are Highly Expressed in the Liver and in the Brain

Tissue-based enrichment analysis using the ARCHS4 database [68] suggests that autophagy modulation targets are highly expressed in the liver and in the brain (including various brain regions: cingulate gyrus, superior frontal gyrus, prefrontal cortex, dorsal striatum, and midbrain) (Figure 3c). This could be attributed to the occurrence of aggregation-prone diseases in both organs, e.g., ATZ aggregation in the liver, and many neurodegenerative disorders in the brain. Compared to other tissues, 165 targets were found to have elevated expression levels in the brain using the Human Protein Atlas [69] (Table S2).
Overall, our results in the above two subsections show that target proteins of autophagy modulators are commonly expressed in the brain and liver tissues and mainly function as tyrosine kinases, calcium channels, aminergic GPCRs, or nuclear hormone receptors. Few frequent targets directly take part in the canonical autophagy pathway, but many of them regulate the expression and activity of core ATG proteins through transcriptional regulation and signaling pathways. We next carried out an in-depth analysis of these pathways to understand the pharmacological mechanisms that underly autophagy modulation.

2.5. Functional Analysis of the Targets Reveals Enriched Pathways Implicated in Autophagy Modulation

Figure 4a shows the gene ontology (GO) annotation terms enriched in targets of autophagy modulators. The top 20 GO cellular components include the nucleus, extrinsic components of the cytoplasmic side of the plasma membrane, and voltage-gated calcium/sodium channel complex. The top 20 GO molecular functions include various protein kinase activity-related and ion channel-related terms. The top 20 biological processes comprise a variety of protein phosphorylation-related and signal transduction-related processes (see Table S4 for an extended list). These results are consistent with our observations in Section 2.3 and highlight the importance of RTK signaling, Ca2+ signaling, and downstream gene transcription regulation. To quantitatively assess the major pathways affected by autophagy modulators, we carried our pathway enrichment analysis using the KEGG database [70].
Our analysis indicated that the known and predicted targets of autophagy modulators take part in 294 pathways (Table S5). Figure 4b highlights the top 20 enriched pathways in three distinct KEGG groups: signal transduction, cellular processes, and organismal systems. The enriched pathways in other groups (human diseases, metabolism, and genetic information processing) can be found in Table S5. We found that the signal transduction group contains three tiers of autophagy regulation: (1) pathways that directly involve core ATG proteins: mTOR-, PI3K-AKT-, and AMPK-pathways; (2) upstream signaling pathways that regulate the activity of ATG proteins: Ca2+-, cAMP-, and RTK/Rap1/Ras/MAPK-pathways; and (3) transcription pathways that regulate the expression of ATG proteins, mainly HIF-1-, NF-κB-, and FoxO-pathways and their regulators such as nuclear hormone receptors NR1I2 and ESR1.
Though the enrichment scores of pathways in the cellular processes group are relatively low compared to those in the signal transduction group (due to the larger number of proteins involved in the cellular processes pathways), three autophagy pathways in KEGG were captured by our analysis, including autophagy-animal, autophagy-others, and mitophagy. Various programmed cell death pathways also exhibit relatively high scores, including apoptosis, necroptosis, and ferroptosis. This is due to a complex crosstalk between autophagy and programmed cell death pathways [71,72,73], which has been shown to process cellular stress signals and confer cell-fate decisions [11]. Interestingly, five cell adhesion pathways are top-ranked among cellular processes: regulation of actin cytoskeleton, focal adhesion, gap junction, adherens junction, and tight junction. These pathways regulate the cytoskeletal dynamics, which play a vital role in the biogenesis of autophagosomes from the ER membrane as well as their translocation and fusion with the lysosomes [74]. Other enriched pathways in the organismal systems class relate to the physiological functions of autophagy, such as neuron survival/functioning [75], insulin production/sensitivity [76], circadian rhythm [77], and longevity [78].
We further analyzed the distribution of pathways for each category of autophagy modulators using QuartataWeb [28]. Table S5 presents the list of pathways with high enrichment scores (p value < 0.05), including 63 targeted by activators, 24 targeted by inhibitors, and 106 targeted by dual-modulators. These pathways were further classified into finer subclasses of KEGG pathways. The subclasses of pathways enriched in targets of autophagy activators, inhibitors, and dual-modulators are shown in Figure 4c. Signal transduction pathways dominate all three categories, along with the endocrine system, a subclass of the organismal systems. Among the top-ranked endocrine system pathways, we distinguish the estrogen signaling pathway targeted by 131 activators, 25 inhibitors, and 14 dual-modulators; whereas melanogenesis is exclusively targeted by autophagy activators (see Table S5). Recent studies show that rottlerin, an autophagy activator, regulates melanogenesis by targeting the cAMP-CREB signaling pathway [79].
Cancer pathways are also commonly targeted by autophagy modulators (Figure 4c). This is particularly clear in the case of inhibitors. Our results reflect the fact that many autophagy modulators are anti-cancer drugs, which trigger programed cell death and introduce cellular stresses. In response to such stresses, cancer cells may trigger autophagy at the same time to protect themselves. Autophagy is a first rescue mechanism for cells before proceeding to cell death in the presence of stronger insults [11], and its interference may thus delay, if not prevent, tumor cell death. Thus, autophagy inhibition has been identified as a strategy of cancer therapeutics [80]. In particular, the breast cancer pathway is highly enriched in targets of autophagy inhibitors (p value = 1.95 × 10−5) but not activators or dual-modulators. This is consistent with the observation that autophagy inhibitors such as SB02024 increase the sensitivity of breast cancer cells to chemotherapy [81].
We also noticed that autophagy inhibitors target more immune system pathways and infection diseases than activators and dual-modulators (respective green and cyan portions of the bars in Figure 4c). Immune pathways specific to autophagy inhibition include Fc receptor (FcR) signaling, Toll-like receptor (TLR) signaling, T cell receptor (TCR) signaling, c-type lectin receptor (CLR) signaling, and chemokine signaling pathways (Table S5). It has been reported that stimulation of TLRs, FcRs, TCRs, and CLRs can induce autophagy [82,83]. We thus hypothesize that the mechanism of action of selected autophagy inhibitors is through blockage of these pathways.
Taken together, our enrichment analysis reveals that autophagy activators, inhibitors, and dual-modulators usually affect the same group of pathways, although a few selected pathways seem to be specifically targeted by autophagy activators (e.g., melanogenesis pathway) or inhibitors (e.g., TLR pathway). Signal transduction mechanisms shared by the three categories of autophagy modulators include those regulating the activity of ATG proteins directly (e.g., AKT pathway) or indirectly (e.g., calcium pathway), in addition to those involved in their transcriptional regulation (e.g., NF-κB and NHR pathways). We next investigated how different signal transduction pathways interconnect to control/regulate autophagy.

2.6. PKA, PI3K, AKT, and mTOR Play a Central Signal Transduction Role in Mediating Autophagy

We identified seven signal transduction pathways highly enriched in targets of all three categories of autophagy modulators: cAMP-, cGMP-PKG-, calcium-, Rap1, Ras, MAPK-, and PI3K-AKT-signaling pathways, schematically shown in Figure 5a. These pathways each contain a large number (between 43 and 108) of targets of autophagy modulators and they were each targeted by a large number (39 to 73) of modulators (see the rows colored light blue in Table S5).
Next, we analyzed the known and predicted targets of autophagy modulators and identified those shared by one or more of these pathways as well as the canonical autophagy pathway from KEGG. This led to 15 ATG-signaling targets. These target proteins, distinguished by their enhanced propensity to be involved in ATG signal transduction, are PKA, PI3K, AKT, MEK1/2, Raf1, JNK1, ERK, Bad, INS, IGFR, mTOR, PP2A, Bcl-2, PDK1, AMPK, and TBK1. Among them, PKA and PI3K each take part in four of the seven signal transduction pathways (Figure 5b), and AKT in six pathways, and as such they play a key role in mediating autophagic responses. A detailed mapping of these targets to signal transduction pathways is presented in Figure 5b. Three color-coded entries therein indicate whether the ATG-signaling target is associated with an autophagy activator (green), inhibitor (red), and/or dual-modulator (blue). We note in Figure 5b (and in Table S2) that PKA, PI3K, and mTOR are targeted by all three categories of autophagy modulators, further supporting the key role of PKA and PI3K, and highlighting that of mTOR as well. The diagram also shows that AKT, MEK1/2, Raf1, and ERK take part in at least five pathways. Interestingly, the MAPK pathway components MEK1/2, Raf1 and JNK1 are targeted by both activators and dual-modulators but not inhibitors. In contrast, AKT and IGFR are targeted by activators and inhibitors, but not dual-modulators.
Overall, this analysis draws attention to the central role of PKA as a mediator of ATG signal transduction events. PKA regulates three key players, mTOR, MEK, and Raf1, which in turn regulate autophagy (Figure 5a). We also note the critical role of PI3K that regulates AKT, which in turn regulates mTOR. We next constructed a comprehensive model of an ATG signal transduction network that provides an integrated view of these complex interactions.

2.7. An Integrative Roadmap for Autophagy Modulation

We constructed a unified signaling network that integrates the above described pathways, which allows us to visualize how targets of autophagy modulators regulate autophagy (see Figure 6a). Core ATG proteins are colored yellow; the associated transcription factors are in light blue. The above described major ATG-signal transduction pathways are displayed in distinct colors as described in the caption, and the crosstalk between these pathways are schematically shown in Figure 6b. We also indicate in gray boxes the six proteins that are targeted by all three categories of autophagy modulators. This comprehensive network of ATG-signal transduction pathways further highlights the role of PKA (at the center) as a key mediator of the effects of autophagy modulators, bridging various signal transduction events to autophagy and regulating the activity and gene expression of ATG proteins. Yellow ellipses highlight key ATG-signal transduction proteins.
This model provides a comprehensive framework for quantitative analysis of ATG-signal transduction pathways and clears the way to building polypharmacological treatment hypotheses. Figure 6a indicates several autophagy modulators that target different components, including those binding to upstream channels/receptors such as calcium or potassium channels (top left), the GPCR ADORA2A (top right), the downstream key signal transduction elements such as PI3K, AKT, and mTORC1, or transcription factors such as HIF-1α.

3. Discussion

Pharmacological modulation of autophagy has great therapeutic potential and attracted a great deal of attention in the past two decades [18]. A variety of autophagy-focused databases have been developed to provide annotation information for autophagy-related proteins (such as the Autophagy Database [84]), the interactions of autophagy components (e.g., Autophagy Regulatory Network [85], THANATOS [86]), and autophagy-modulating small molecules (e.g., AutophagySMDB [87], HAMdb [88]). These resources contain extensive data on autophagy that enables systems-level analyses of its regulation. In the present study, we mined these databases, manually curated the data in the light of recent literature, and identified a structurally and functionally diverse set of 225 drugs that are known to activate, inhibit, or modulate autophagy, as supported by direct experimental evidence.
The experimentally verified (known) targets of these drugs are available in databases such as DrugBank [25]. However, we are still far from a complete understanding of their actions and side effects. Quantitative systems pharmacology (QSP) methods/tools have been developed with an aim to utilize existing knowledge on networks of interactions among drugs, proteins, and genes to discover new drug targets, reveal underlying mechanisms, and identify repurposed drugs [89,90]. For example, matrix factorization-based machine learning algorithms have been applied to the prediction of novel drug-target interactions [26,27,91]. Pathway/network analysis tools have been developed to study drug effects [28,92,93,94].
In this study, our QSP analysis led to the identification of 1837 experimentally verified (known) and 368 computationally predicted (new) interactions between these autophagy modulators and their 993 known and 12 new target proteins, involved in 294 pathways. Many new interactions predicted using our machine learning method through QuartataWeb interface [28] were found to be consistent with recent experiments. For example, an interaction between fluphenazine and dopamine receptor D3 that was not reported in DrugBank (and therefore not included in our input dataset) was predicted with high confidence. Notably, this drug was recently reported to have a subnanomolar binding affinity on DRD3 (Ki = 0.11 nM) [50]. Among targets we predicted, the translocator protein (TSPO) has been reported to inhibit mitophagy downstream of the PTEN-induced kinase 1 (PINK1)-Parkin pathway [95]. Overall, we found experimental evidence in support of 75 (out of 368) predicted pairs, strongly suggesting that the other predictions could as well indicate possible interactions (either repurposable actions or side effects) of existing drugs. Such machine-learning based predictions may thus help generate novel hypotheses worth further testing or interpret experimental observations. As an example, we predict that olanzapine may interact with calmodulin, implying that the reported olanzapine-induced AMPK activation may be partially mediated by the calmodulin-CaMKKβ-AMPK pathway [42,43].
The dominant paradigm in drug discovery based on Ehrlich’s philosophy of magic bullets (drugs selective for single targets) has been challenged by the new concept of magic shotguns (selectively non-selective drugs) [96], due to the recognition of the complexity of diseases such as cancer and neurological disorders and the prevalence of multitarget drugs approved in recent years [97]. In many cases, perturbing individual targets has little effect on disease networks, and thus therapies are required to modulate multiple targets to achieve efficacy. Therefore, though the promiscuity of drugs may lead to unintended side-effects, it could also be an advantage in terms of improving therapeutic efficacy and preventing drug resistance [98]. Several autophagy modulators are found here to be highly promiscuous, including enzyme cofactors copper (145 targets) and zinc (124 targets), and compounds such as fostamatinib (303 targets), artenimol (78 targets), and olanzapine (48 targets). Fostamatinib, the most promiscuous autophagy modulator, is a type I kinase inhibitor targeting the ATP binding pocket of kinases, which is highly conserved across the human kinome. Pharmacological profiling showed that fostamatinib inhibits 120 kinases with IC50 values ranging from 3 nM to 3.49 μM [99]. These kinases include crucial ATG proteins such as mTOR, ULK1/2, CaMKKβ, PI3K, AKT, TBK1, PKG, PKA, PDK1, JNK, RAF1, and RTKs (Figure 6a). The atypical antipsychotic olanzapine, on the other hand, binds to receptors HTR2A/2C3/6, DRD1-4, HRH1, ADRA1s, and CHRN1-5 with high affinities (Ki ranges from 4 to 132 nM) at conserved binding sites (e.g., residue C/S3.36) [100]. Olanzapine has been reported to modulate autophagy in both neuronal (e.g., SH-SY5H cells) and non-neuronal cells (e.g., LN229 and T98 glioma cells) by upregulating reactive oxygen species (ROS) and ATG gene expression [44].
The dominance of the kinase family members among targets of autophagy modulators directed us to focus more closely on interconnected signal transduction pathways implicated in autophagy modulation (Figure 5). The dominant role of kinases is consistent with a study of the human kinase-substrate post-translational modification network where phosphorylation plays a key role in regulating autophagy [86]. In addition to the kinases highlighted in Figure 5, as illustrated in Figure 6, core autophagy proteins such as ULK1/2 and AMPK are all kinases. MAPKs and PKA crosstalk to other pathways and thus their roles in regulating autophagy depend on the cellular context. Since kinase pathways are activated in several cancers, differential regulation of autophagy by targeting multiple kinase pathways has been suggested for the development of better cancer therapeutics [101].
Another important class of proteins broadly used as drug targets, GPCRs, is also targeted for autophagy modulation. In this case, GPCRs regulate second messengers Ca2+ and cAMP, and downstream effectors MAPKs, which initiate the major upstream signal transduction that leads to autophagy. It has been suggested that the regulation of autophagy by GPCRs is essential to mediating the effects of hormones/neurotransmitters secreted in response to systemic nutrient fluctuations [58]. Many GPCRs (e.g., ADRA1s) and calcium channels (e.g., CACNA1s) are frequent targets, showing that the corresponding autophagy modulators are mediated by upstream signaling pathways, rather than the canonical mTOR pathway.
Turning to the NHRs, our analysis highlighted two important transcription factors directedly targeted by autophagy modulators HIF-1α (targeted by PX-478 and 2-methoxyestradiol) and NF-κB (targeted by glucosamine and acetylsalicylic acid). Notably, both of these, and other transcription factors such as FoxO and TFEB, are regulated by estrogen receptor/NHR ESR1, hence the frequent observation of ESR1 as a major target. Another NHR distinguished by its high promiscuity was NR1I2, which directly regulates the kinase pathways.
The key autophagy proteins that these transcription factors and/or NHRs regulate (e.g., Beclin-1, LC3, Atg5) are illustrated in Figure 6a. Though transcription factors such as TFEB, TFE3, and FoxO are not directly targeted by autophagy modulators, their regulators such as AKT, mTOR, and AMPK appear to be targets. The presence of FoxO, HIF-1-, and NF-κB signaling pathways among the top 20 enriched pathways (Figure 4) further confirmed that the regulation of autophagic gene expression is an important mechanism exploited by many autophagy modulators.
Note that Figure 6a only illustrated representative crosstalks between the pathways mentioned above. The signaling network regulating autophagy may be rewired by specific crosstalk in different cell types. For example, in KRAS-driven cancers, calmodulin binds KRAS, and can directly activate PI3K and AKT [102,103]. Further, our classification of autophagy modulators as activators, inhibitors, and dual-modulators is limited by the existing experimental evidence. Though many targets only interact with either activators or inhibitors, given the complex pathway crosstalk shown in Figure 6a, it is possible that current autophagy activators and inhibitors are in fact dual-modulators. For example, the most promiscuous autophagy inhibitor fostamatinib is known to inhibit serine/threonine-protein kinase tousled-like 2 (TLK2), a negative regulator of amino acid starvation-induced autophagy [99]. Thus, fostamatinib may enhance autophagy via TLK2, although there is no experimental data in support of this. Furthermore, Rev-Erb agonists SR-9009 and SR-9011 inhibit autophagy through the nuclear receptor subfamily 1 group D members 1 and 2 (NR1D1/2) [104]. Since NR1D1 mediates the regulation of autophagy rhythms by the circadian clock [104], whether the inhibition of NR1D1 leads to the enhancement or suppression of autophagy might depend on the timing and cellular state. Melanogenesis and TLR signaling pathways are exclusively targeted by autophagy activators or inhibitors, respectively. Whether it implies specific mechanisms underlying autophagy activation and inhibition requires further investigations.
PKA, downstream of cAMP signaling, takes part in all three tiers of autophagy regulation: (i) it inhibits core autophagy proteins AMPK and mTOR; (ii) it regulates upstream MAPK signaling; and (iii) it inhibits NF-κB signaling. The drugs targeting PKA involve autophagy activators (e.g., fasudil), inhibitors (e.g., fostamatinib), and dual-modulators (e.g., metformin). PKA is a key regulator that governs various cellular processes. Interestingly, a recent high-throughput screening study in yeast suggested that many substrates or downstream autophagy proteins (e.g., Atg10, Atg14, Atg15, and Atg16) regulated by PKA are, in turn, also PKA regulators [105]. This implies the presence of multiple positive and negative feedback loops coupling PKA to autophagic responses. It has been reported that mammalian PKA activates mTORC1 while, in yeast, TORC1 can indirectly activate PKA through a homolog of a ribosomal protein S6 kinase B1 (RPS6KB1) [106]. All these observations point to the utility of a thorough study of the feedback regulation between PKA and Atg proteins in human cells. In this context, dynamical modeling approaches promise to offer helpful pointers [11,107,108,109,110].
As shown in Figure 4, the targets of autophagy modulators are involved in neurological disease, infectious disease, cancer, cardiovascular disease, and metabolic disease pathways; they are upregulated in the liver and brain tissues (Figure 3c). This observation may be partially due to the presence of antipsychotic drugs (e.g., olanzapine and pimozide) and liver disease drugs (e.g., carbamazepine and resveratrol) in our input list of autophagy modulators, both developed for eliminating toxic aggregates. However, it also reflects the crucial role of autophagy in liver and brain health and diseases [111,112,113] and justifies the attempts of repurposing autophagy modulators against liver and brain diseases [21,22,23].
The targets involved in cardiovascular diseases are mainly calcium channels. Intracellular Ca2+ can both up- and downregulate autophagy and plays a key role in the crosstalk between autophagy and apoptosis [11]. The interplay between autophagy and apoptosis renders the opposing roles of autophagy as a therapeutic target for cardiovascular diseases. For example, in ischemic heart disease, enhanced autophagy is desired for cardioprotective effects via ROS reduction; however, alleviating chemotherapy-induced cardiotoxicity requires suppression of autophagy [114].
Apart from cardiovascular diseases, dual-modulation of autophagy is also relevant to designing therapeutic strategies for other diseases such as cancer, infectious diseases, and neurodegenerative diseases. For example, basal autophagy suppresses tumor initiation but promotes cancer progression. During infection, enhancing autophagy facilitates intracellular pathogen clearance but also suppresses immune responses [18]. Thus, precise modulation and time-sensitive strategies are required for autophagy-based therapeutic options. We also note that several drugs proposed for treating neurodegenerative diseases via autophagy are in fact dual-modulators, such as vorinostat, nimodipine, rosiglitazone, sitaglipin, and dexmedetomidine for AD, trifluoperazine and melatonin for Parkinson’s disease (PD), lithium for amyotrophic lateral sclerosis (ALS), berberine for trauma-induced neurodegeneration, ginsenoside Rb for Huntington’s disease, and emodin for hyperhomocysteinemia-induced dementia [115]. Fifty-seven targets of dual-modulators are involved in the calcium signaling pathway and five are microtubule-related proteins: tubulin β and β1 chain (TUBB and TUBB1) and microtubule-associated proteins 2, 4, and τ (MAP2/4/T), which are known to regulate neurodegeneration. For example, mutant MAPT induces disorder and instability in microtubules, which are associated with AD, PD, and ALS [116,117]. Therefore, a deeper understanding of the precise role of dual-modulators in specific disease settings is needed, and caution should be exercised in the use of autophagy dual-modulators for therapeutic purposes.
Our enrichment analysis also shows that many targets of autophagy modulators are involved in apoptosis (e.g., p53, Bcl-2, and caspase 8), necroptosis (e.g., caspase 8, RIPK1, and TNF), and ferroptosis (e.g., p53, 15LOX, and ACSL4) pathways (Figure 4 and Table S5). In recent work, we have demonstrated the pivotal roles of these targets in regulating distinct cell death pathways [107,118,119,120]. We also found that due to complex crosstalks among pathways, targeting autophagy, necroptosis, and ferroptosis using polypharmacological strategies can lead to better therapeutic effects in terms of cell survival [121,122]. The present study shows that the key regulators of cell death are targeted by autophagy modulators, in accord with the crosstalk between autophagy and these cell death pathways. The integrated network constructed based on our analysis (Figure 6) may serve as an excellent framework for pursuing quantitative analyses shedding light on the biological and biomedical implications of such crosstalks and opening new directions for designing polypharmacological strategies for monitoring or altering cell survival/death decisions.

4. Materials and Methods

4.1. Data Collection

We selected 225 drugs from DrugBank v5.1.1, which have at least one target and are also known as autophagy modulators (Table S1). We classified the selected modulators into three groups: activators, inhibitors, and dual-modulators by manually searching the literature for direct experimental evidence in support of their specific roles. A dataset of 993 known targets were retrieved from DrugBank (Table S2).

4.2. Drug–Target Interaction Prediction

New drug–target interactions were predicted using our PMF-based machine learning algorithm [26,27,28]. We represented the known interactions between N drugs and M targets as a sparse matrix R . We then decomposed R into a drug matrix U and a target matrix V , by learning the optimal D latent variables to represent each drug and each target using the PMF algorithm. The product of U T and V reconstructed R by assigning each of the unknown entries in R a value, which represents the confidence score for a novel drug–target interaction prediction:
R N × M = U N × D T V D × M
Using this method, we trained a PMF model based on 14,983 drug–target interactions between 5494 drugs and 2807 targets from DrugBank 5.1.1. We evaluated the confidence scores in the range [0, 1] for each predicted drug–target interaction. We selected the interactions with confidence scores higher than 0.6 within the top 40 predicted targets for each input drug. Performance benchmarking showed that using 0.6 as the threshold confidence score for predictions achieved a good compromise between sensitivity vs. specificity [28]. This led to 368 novel interactions (Table S3).

4.3. Drug–Drug Similarity Analysis

We evaluated the structural similarity between each pair of autophagy modulators by calculating the Tanimoto distance between their 2D structure fingerprints using the RDKit software [123]. We evaluated the functional similarity of two drugs based on their interaction patterns with known targets. Specifically, we represented each drug i by a vector d i with values of 0 or 1 (depending on the absence or existence of an interaction between this drug and the corresponding target) in the matrix R . Interaction-pattern similarities between drugs i and j were then evaluated by calculating the cosine distance between their vectors d i and d j :
1 ( d i · d j ) / ( | d i |     | d j | ) .

4.4. Pathway Enrichment Analysis

We enriched the KEGG pathways (version March 2018, homo sapiens) [124] by the 225 drugs with 993 known and 12 predicted targets. Two hundred ninety-four pathways were identified (see Table S5). We calculated the hypergeometric p values to rank enriched pathways. Our enrichment scores were defined as follows: Given a list of targets, the enrichment p value for pathway i ( P i ) is the probability of randomly drawing k 0 or more targets that belong to P i :
P i = k 0 k m ( K k ) ( M K m k ) ( M m )
where m is the total number of targets we identified, M is the total number of human pathways in KEGG, K is the number of targets that take part in the examined pathway i, and k0 is the number of targets we identified that are involved in pathway i. In order to account for multiple testing, we used the Benjamini–Hochberg method [125] to adjust the calculated p values by introducing a false discovery rate (FDR) correction. The FDR is the fraction of false significant pathways maximally expected from the significant pathways we identified, which is bounded by the cutoff of the adjusted p values. Specifically, the adjusted p value, p i * , for the ith pathway is:
p i * = m i n k = i m { min ( p k m i , 1 ) }
where m is the total number of pathways. Table S5 lists these p values for pathway enrichments based on both known and predicted targets.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/8/2855/s1. The source code and supplementary data are available at http://github.com/liubing1020/autophagy. Figure S1: Distribution of the autophagy modulators based on their structure (left) and interaction pattern (right) similarities; Figure S2: Population of autophagy modulators with different degrees (different numbers of associated targets); Figure S3: Additional data on the number of known and predicted targets of activators; Table S1: 225 autophagy modulators; Table S2: 993 known targets and 12 predicted targets of autophagy modulators; Table S3: 368 Predicted ATG interactions and experimental validation; Table S4: GO annotation enrichments based on ATG targets; Table S5: Pathways enriched in targets of autophagy modulators, corresponding targets, and autophagy modulators.

Author Contributions

Conceptualization, Q.S., F.P., B.L., and I.B.; methodology, Q.S. and F.P.; investigation, Q.S., F.P., B.L., and I.B.; data curation, Q.S.; writing—original draft preparation, Q.S. and B.L.; writing—review and editing, F.P., G.A.S., S.C.P., D.H.P. and I.B.; visualization, Q.S. and B.L.; supervision, B.L. and I.B.; funding acquisition, D.H.P. and I.B. All authors have read and agreed on the manuscript.

Funding

This research was funded by National Institutes of Health awards P01DK096990, P30DA035778, P19AI068021, and P41GM103712 to I.B.

Acknowledgments

F.P. wishes to thank D. Lansing Taylor for his mentorship.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Wang, Y.; Xu, W.; Yan, Z.; Zhao, W.; Mi, J.; Li, J.; Yan, H. Metformin induces autophagy and G0/G1 phase cell cycle arrest in myeloma by targeting the AMPK/mTORC1 and mTORC2 pathways. J. Exp. Clin. Cancer Res. 2018, 37, 63. [Google Scholar] [CrossRef] [PubMed]
  2. Levine, B.; Kroemer, G. Biological functions of autophagy genes: A disease perspective. Cell 2019, 176, 11–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Rabinowitz, J.D.; White, E. Autophagy and metabolism. Science 2010, 330, 1344–1348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mizushima, N. Autophagy: Process and function. Genes Dev. 2007, 21, 2861–2873. [Google Scholar] [CrossRef] [Green Version]
  5. Levine, B. Cell biology: Autophagy and cancer. Nature 2007, 446, 745–747. [Google Scholar] [CrossRef]
  6. Lorin, S.; Hamai, A.; Mehrpour, M.; Codogno, P. Autophagy regulation and its role in cancer. Semin. Cancer Biol. 2013, 23, 361–379. [Google Scholar] [CrossRef]
  7. Mathew, R.; Karantza-Wadsworth, V.; White, E. Role of autophagy in cancer. Nat. Rev. Cancer 2007, 7, 961–967. [Google Scholar] [CrossRef]
  8. Kim, Y.C.; Guan, K.L. mTOR: A pharmacologic target for autophagy regulation. J. Clin. Investig. 2015, 125, 25–32. [Google Scholar] [CrossRef] [Green Version]
  9. Mizushima, N. The role of the Atg1/ULK1 complex in autophagy regulation. Curr. Opin. Cell Biol. 2010, 22, 132–139. [Google Scholar] [CrossRef]
  10. Graef, M.; Nunnari, J. A role for mitochondria in autophagy regulation. Autophagy 2011, 7, 1245–1246. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, B.; Oltvai, Z.N.; Bayir, H.; Silverman, G.A.; Pak, S.C.; Perlmutter, D.H.; Bahar, I. Quantitative assessment of cell fate decision between autophagy and apoptosis. Sci. Rep. 2017, 7, 17605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Decuypere, J.P.; Bultynck, G.; Parys, J.B. A dual role for Ca(2+) in autophagy regulation. Cell Calcium 2011, 50, 242–250. [Google Scholar] [CrossRef]
  13. Racioppi, L.; Means, A.R. Calcium/calmodulin-dependent protein kinase kinase 2: Roles in signaling and pathophysiology. J. Biol. Chem. 2012, 287, 31658–31665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Yousefi, S.; Perozzo, R.; Schmid, I.; Ziemiecki, A.; Schaffner, T.; Scapozza, L.; Brunner, T.; Simon, H.U. Calpain-mediated cleavage of Atg5 switches autophagy to apoptosis. Nat. Cell Biol. 2006, 8, 1124–1132. [Google Scholar] [CrossRef] [PubMed]
  15. Pan, T.; Rawal, P.; Wu, Y.; Xie, W.; Jankovic, J.; Le, W. Rapamycin protects against rotenone-induced apoptosis through autophagy induction. Neuroscience 2009, 164, 541–551. [Google Scholar] [CrossRef]
  16. Li, Y.; Wang, X.; Wei, Z.; Mao, H.; Gao, M.; Liu, Y.; Ma, Y.; Liu, X.; Guo, C.; Zhang, L.; et al. Pretreatment with wortmannin alleviates lipopolysaccharide/d-galactosamine-induced acute liver injury. Biochem. Biophys. Res. Commun. 2014, 455, 234–240. [Google Scholar] [CrossRef]
  17. O’Reilly, L.P.; Long, O.S.; Cobanoglu, M.C.; Benson, J.A.; Luke, C.J.; Miedel, M.T.; Hale, P.; Perlmutter, D.H.; Bahar, I.; Silverman, G.A.; et al. A genome-wide RNAi screen identifies potential drug targets in a C. elegans model of alpha1-antitrypsin deficiency. Hum. Mol. Genet. 2014, 23, 5123–5132. [Google Scholar] [CrossRef] [Green Version]
  18. Galluzzi, L.; Bravo-San Pedro, J.M.; Levine, B.; Green, D.R.; Kroemer, G. Pharmacological modulation of autophagy: Therapeutic potential and persisting obstacles. Nat. Rev. Drug Discov. 2017, 16, 487–511. [Google Scholar] [CrossRef] [Green Version]
  19. Towers, C.G.; Thorburn, A. Therapeutic targeting of autophagy. EBioMedicine 2016, 14, 15–23. [Google Scholar] [CrossRef] [Green Version]
  20. Cai, Z.; Yan, L.J. Rapamycin, autophagy, and alzheimer’s disease. J. Biochem. Pharmacol. Res. 2013, 1, 84–90. [Google Scholar]
  21. Hidvegi, T.; Ewing, M.; Hale, P.; Dippold, C.; Beckett, C.; Kemp, C.; Maurice, N.; Mukherjee, A.; Goldbach, C.; Watkins, S.; et al. An autophagy-enhancing drug promotes degradation of mutant alpha1-antitrypsin Z and reduces hepatic fibrosis. Science 2010, 329, 229–232. [Google Scholar] [CrossRef] [PubMed]
  22. Li, J.; Pak, S.C.; O’Reilly, L.P.; Benson, J.A.; Wang, Y.; Hidvegi, T.; Hale, P.; Dippold, C.; Ewing, M.; Silverman, G.A.; et al. Fluphenazine reduces proteotoxicity in C. elegans and mammalian models of alpha-1-antitrypsin deficiency. PLoS ONE 2014, 9, e87260. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, Y.; Cobanoglu, M.C.; Li, J.; Hidvegi, T.; Hale, P.; Ewing, M.; Chu, A.S.; Gong, Z.; Muzumdar, R.; Pak, S.C.; et al. An analog of glibenclamide selectively enhances autophagic degradation of misfolded alpha1-antitrypsin Z. PLoS ONE 2019, 14, e0209748. [Google Scholar] [CrossRef]
  24. Cuomo, F.; Altucci, L.; Cobellis, G. Autophagy function and dysfunction: Potential drugs as anti-cancer therapy. Cancers (Basel) 2019, 11, 1465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef] [PubMed]
  26. Cobanoglu, M.C.; Liu, C.; Hu, F.; Oltvai, Z.N.; Bahar, I. Predicting drug-target interactions using probabilistic matrix factorization. J. Chem. Inf. Model. 2013, 53, 3399–3409. [Google Scholar] [CrossRef]
  27. Cobanoglu, M.C.; Oltvai, Z.N.; Taylor, D.L.; Bahar, I. BalestraWeb: Efficient online evaluation of drug-target interactions. Bioinformatics 2015, 31, 131–133. [Google Scholar] [CrossRef] [Green Version]
  28. Li, H.; Pei, F.; Taylor, D.L.; Bahar, I. QuartataWeb: An integrated chemical-protein interaction prediction and pathway inference server for polyphamalogical and chemogenomics analysis. Bioinformatics 2020. [Google Scholar] [CrossRef]
  29. Sagrillo-Fagundes, L.; Bienvenue-Pariseault, J.; Vaillancourt, C. Melatonin: The smart molecule that differentially modulates autophagy in tumor and normal placental cells. PLoS ONE 2019, 14, e0202458. [Google Scholar] [CrossRef] [Green Version]
  30. Brachmann, S.; Fritsch, C.; Maira, S.M.; Garcia-Echeverria, C. PI3K and mTOR inhibitors: A new generation of targeted anticancer agents. Curr. Opin. Cell Biol. 2009, 21, 194–198. [Google Scholar] [CrossRef]
  31. Moriya, S.; Che, X.F.; Komatsu, S.; Abe, A.; Kawaguchi, T.; Gotoh, A.; Inazu, M.; Tomoda, A.; Miyazawa, K. Macrolide antibiotics block autophagy flux and sensitize to bortezomib via endoplasmic reticulum stress-mediated CHOP induction in myeloma cells. Int. J. Oncol. 2013, 42, 1541–1550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, D.; Ji, X.; Liu, J.; Li, Z.; Zhang, X. Dopamine receptor subtypes differentially regulate autophagy. Int. J. Mol. Sci. 2018, 19, 1540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Trejo-Solis, C.; Jimenez-Farfan, D.; Rodriguez-Enriquez, S.; Fernandez-Valverde, F.; Cruz-Salgado, A.; Ruiz-Azuara, L.; Sotelo, J. Copper compound induces autophagy and apoptosis of glioma cells by reactive oxygen species and JNK activation. BMC Cancer 2012, 12, 156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Liuzzi, J.P.; Guo, L.; Yoo, C.; Stewart, T.S. Zinc and autophagy. Biometals 2014, 27, 1087–1096. [Google Scholar] [CrossRef]
  35. Tsang, T.; Posimo, J.M.; Guidiel, A.A.; Cicchini, M.; Feldser, D.M.; Brady, D.C. Copper is an essential regulator of the autohagic kinases ULK1/2 to dirve lung adenocarcinoma. Nat. Cell Biol. 2020, in press. [Google Scholar] [CrossRef]
  36. Olliaro, P.L.; Haynes, R.K.; Meunier, B.; Yuthavong, Y. Possible modes of action of the artemisinin-type compounds. Trends Parasitol. 2001, 17, 122–126. [Google Scholar] [CrossRef]
  37. Wang, J.; Zhang, C.J.; Chia, W.N.; Loh, C.C.; Li, Z.; Lee, Y.M.; He, Y.; Yuan, L.X.; Lim, T.K.; Liu, M.; et al. Haem-activated promiscuous targeting of artemisinin in Plasmodium falciparum. Nat. Commun. 2015, 6, 10111. [Google Scholar] [CrossRef]
  38. Shi, X.; Wang, L.; Li, X.; Bai, J.; Li, J.; Li, S.; Wang, Z.; Zhou, M. Dihydroartemisinin induces autophagy-dependent death in human tongue squamous cell carcinoma cells through DNA double-strand break-mediated oxidative stress. Oncotarget 2017, 8, 45981–45993. [Google Scholar] [CrossRef] [Green Version]
  39. Konstat-Korzenny, E.; Ascencio-Aragon, J.A.; Niezen-Lugo, S.; Vazquez-Lopez, R. Artemisinin and its synthetic derivatives as a possible therapy for cancer. Med. Sci. (Basel) 2018, 6, 19. [Google Scholar] [CrossRef] [Green Version]
  40. Niture, S.; Gyamfi, M.A.; Kedir, H.; Arthur, E.; Ressom, H.; Deep, G.; Kumar, D. Serotonin induced hepatic steatosis is associated with modulation of autophagy and notch signaling pathway. Cell Commun. Signal. 2018, 16, 78. [Google Scholar] [CrossRef] [Green Version]
  41. Kim, J.K.; Kim, Y.S.; Lee, H.M.; Jin, H.S.; Neupane, C.; Kim, S.; Lee, S.H.; Min, J.J.; Sasai, M.; Jeong, J.H.; et al. GABAergic signaling linked to autophagy enhances host protection against intracellular bacterial infections. Nat. Commun. 2018, 9, 4184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ikegami, M.; Ikeda, H.; Ishikawa, Y.; Ohsawa, M.; Ohashi, T.; Kai, M.; Kamei, A.; Kamei, J. Olanzapine induces glucose intolerance through the activation of AMPK in the mouse hypothalamus. Eur. J. Pharmacol. 2013, 718, 376–382. [Google Scholar] [CrossRef] [PubMed]
  43. Schmidt, R.H.; Jokinen, J.D.; Massey, V.L.; Falkner, K.C.; Shi, X.; Yin, X.; Zhang, X.; Beier, J.I.; Arteel, G.E. Olanzapine activates hepatic mammalian target of rapamycin: New mechanistic insight into metabolic dysregulation with atypical antipsychotic drugs. J. Pharmacol. Exp. Ther. 2013, 347, 126–135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Vucicevic, L.; Misirkic-Marjanovic, M.; Paunovic, V.; Kravic-Stevovic, T.; Martinovic, T.; Ciric, D.; Maric, N.; Petricevic, S.; Harhaji-Trajkovic, L.; Bumbasirevic, V.; et al. Autophagy inhibition uncovers the neurotoxic action of the antipsychotic drug olanzapine. Autophagy 2014, 10, 2362–2378. [Google Scholar] [CrossRef]
  45. Shinde, A.; Hardy, S.D.; Kim, D.; Akhand, S.S.; Jolly, M.K.; Wang, W.H.; Anderson, J.C.; Khodadadi, R.B.; Brown, W.S.; George, J.T.; et al. Spleen tyrosine kinase-mediated autophagy is required for epithelial-mesenchymal plasticity and metastasis in breast cancer. Cancer Res. 2019, 79, 1831–1843. [Google Scholar] [CrossRef] [Green Version]
  46. Duan, W.; Guo, M.; Yi, L.; Zhang, J.; Bi, Y.; Liu, Y.; Li, Y.; Li, Z.; Ma, Y.; Zhang, G.; et al. Deletion of Tbk1 disrupts autophagy and reproduces behavioral and locomotor symptoms of FTD-ALS in mice. Aging (Albany N.Y.) 2019, 11, 2457–2476. [Google Scholar] [CrossRef]
  47. Singh, P.; Ravanan, P.; Talwar, P. Death associated protein kinase 1 (DAPK1): A regulator of apoptosis and autophagy. Front. Mol. Neurosci. 2016, 9, 46. [Google Scholar] [CrossRef] [Green Version]
  48. Ber, Y.; Shiloh, R.; Gilad, Y.; Degani, N.; Bialik, S.; Kimchi, A. DAPK2 is a novel regulator of mTORC1 activity and autophagy. Cell Death Differ. 2015, 22, 465–475. [Google Scholar] [CrossRef] [Green Version]
  49. Manzoni, C.; Lewis, P.A. LRRK2 and Autophagy. Adv. Neurobiol. 2017, 14, 89–105. [Google Scholar] [CrossRef]
  50. Davis, M.I.; Hunt, J.P.; Herrgard, S.; Ciceri, P.; Wodicka, L.M.; Pallares, G.; Hocker, M.; Treiber, D.K.; Zarrinkar, P.P. Comprehensive analysis of kinase inhibitor selectivity. Nat. Biotechnol. 2011, 29, 1046–1051. [Google Scholar] [CrossRef]
  51. Li, P.; Snyder, G.L.; Vanover, K.E. Dopamine targeting drugs for the treatment of schizophrenia: Past, present and future. Curr. Top. Med. Chem. 2016, 16, 3385–3403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Besnard, J.; Ruda, G.F.; Setola, V.; Abecassis, K.; Rodriguiz, R.M.; Huang, X.P.; Norval, S.; Sassano, M.F.; Shin, A.I.; Webster, L.A.; et al. Automated design of ligands to polypharmacological profiles. Nature 2012, 492, 215–220. [Google Scholar] [CrossRef] [PubMed]
  53. Pundir, S.; Martin, M.J.; O’Donovan, C. UniProt Protein Knowledgebase. Methods Mol. Biol. 2017, 1558, 41–55. [Google Scholar] [CrossRef] [Green Version]
  54. Santos, R.; Ursu, O.; Gaulton, A.; Bento, A.P.; Donadi, R.S.; Bologa, C.G.; Karlsson, A.; Al-Lazikani, B.; Hersey, A.; Oprea, T.I.; et al. A comprehensive map of molecular drug targets. Nat. Rev. Drug Discov. 2017, 16, 19–34. [Google Scholar] [CrossRef] [PubMed]
  55. Xiang, J.; Liu, X.; Ren, J.; Chen, K.; Wang, H.L.; Miao, Y.Y.; Qi, M.M. How does estrogen work on autophagy? Autophagy 2019, 15, 197–211. [Google Scholar] [CrossRef] [Green Version]
  56. Blessing, A.M.; Rajapakshe, K.; Reddy Bollu, L.; Shi, Y.; White, M.A.; Pham, A.H.; Lin, C.; Jonsson, P.; Cortes, C.J.; Cheung, E.; et al. Transcriptional regulation of core autophagy and lysosomal genes by the androgen receptor promotes prostate cancer progression. Autophagy 2017, 13, 506–521. [Google Scholar] [CrossRef]
  57. Yan, S.; Yang, X.; Chen, T.; Xi, Z.; Jiang, X. The PPARgamma agonist Troglitazone induces autophagy, apoptosis and necroptosis in bladder cancer cells. Cancer Gene Ther. 2014, 21, 188–193. [Google Scholar] [CrossRef] [Green Version]
  58. Wauson, E.M.; Dbouk, H.A.; Ghosh, A.B.; Cobb, M.H. G protein-coupled receptors and the regulation of autophagy. Trends Endocrinol. Metab. 2014, 25, 274–282. [Google Scholar] [CrossRef] [Green Version]
  59. Lizaso, A.; Tan, K.T.; Lee, Y.H. beta-adrenergic receptor-stimulated lipolysis requires the RAB7-mediated autolysosomal lipid degradation. Autophagy 2013, 9, 1228–1243. [Google Scholar] [CrossRef] [Green Version]
  60. Merlin, J.; Evans, B.A.; Csikasz, R.I.; Bengtsson, T.; Summers, R.J.; Hutchinson, D.S. The M3-muscarinic acetylcholine receptor stimulates glucose uptake in L6 skeletal muscle cells by a CaMKK-AMPK-dependent mechanism. Cell Signal. 2010, 22, 1104–1113. [Google Scholar] [CrossRef]
  61. Yan, Y.; Jiang, W.; Liu, L.; Wang, X.; Ding, C.; Tian, Z.; Zhou, R. Dopamine controls systemic inflammation through inhibition of NLRP3 inflammasome. Cell 2015, 160, 62–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Dolma, S.; Selvadurai, H.J.; Lan, X.; Lee, L.; Kushida, M.; Voisin, V.; Whetstone, H.; So, M.; Aviv, T.; Park, N.; et al. Inhibition of dopamine receptor D4 impedes autophagic flux, proliferation, and survival of glioblastoma stem cells. Cancer Cell 2016, 29, 859–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Soll, C.; Jang, J.H.; Riener, M.O.; Moritz, W.; Wild, P.J.; Graf, R.; Clavien, P.A. Serotonin promotes tumor growth in human hepatocellular cancer. Hepatology 2010, 51, 1244–1254. [Google Scholar] [CrossRef]
  64. Fraser, J.; Cabodevilla, A.G.; Simpson, J.; Gammoh, N. Interplay of autophagy, receptor tyrosine kinase signalling and endocytic trafficking. Essays Biochem. 2017, 61, 597–607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Wei, Y.; Zou, Z.; Becker, N.; Anderson, M.; Sumpter, R.; Xiao, G.; Kinch, L.; Koduru, P.; Christudass, C.S.; Veltri, R.W.; et al. EGFR-mediated Beclin 1 phosphorylation in autophagy suppression, tumor progression, and tumor chemoresistance. Cell 2013, 154, 1269–1284. [Google Scholar] [CrossRef] [Green Version]
  66. Yogalingam, G.; Pendergast, A.M. Abl kinases regulate autophagy by promoting the trafficking and function of lysosomal components. J. Biol. Chem. 2008, 283, 35941–35953. [Google Scholar] [CrossRef] [Green Version]
  67. Takeuchi, H.; Kanzawa, T.; Kondo, Y.; Kondo, S. Inhibition of platelet-derived growth factor signalling induces autophagy in malignant glioma cells. Br. J. Cancer 2004, 90, 1069–1075. [Google Scholar] [CrossRef]
  68. Lachmann, A.; Torre, D.; Keenan, A.B.; Jagodnik, K.M.; Lee, H.J.; Wang, L.; Silverstein, M.C.; Ma’ayan, A. Massive mining of publicly available RNA-seq data from human and mouse. Nat. Commun. 2018, 9, 1366. [Google Scholar] [CrossRef] [Green Version]
  69. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  70. Kanehisa, M.; Goto, S.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. Data, information, knowledge and principle: Back to metabolism in KEGG. Nucleic Acids Res. 2014, 42, D199–D205. [Google Scholar] [CrossRef] [Green Version]
  71. Marino, G.; Niso-Santano, M.; Baehrecke, E.H.; Kroemer, G. Self-consumption: The interplay of autophagy and apoptosis. Nat. Rev. Mol. Cell Biol. 2014, 15, 81–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Radogna, F.; Dicato, M.; Diederich, M. Cancer-type-specific crosstalk between autophagy, necroptosis and apoptosis as a pharmacological target. Biochem. Pharmacol. 2015, 94, 1–11. [Google Scholar] [CrossRef] [PubMed]
  73. Liu, J.; Kuang, F.; Kroemer, G.; Klionsky, D.J.; Kang, R.; Tang, D. Autophagy-dependent ferroptosis: Machinery and regulation. Cell Chem. Biol. 2020. [Google Scholar] [CrossRef] [PubMed]
  74. Kast, D.J.; Dominguez, R. The cytoskeleton-autophagy connection. Curr. Biol. 2017, 27, R318–R326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Menzies, F.M.; Fleming, A.; Caricasole, A.; Bento, C.F.; Andrews, S.P.; Ashkenazi, A.; Fullgrabe, J.; Jackson, A.; Jimenez Sanchez, M.; Karabiyik, C.; et al. Autophagy and neurodegeneration: Pathogenic mechanisms and therapeutic opportunities. Neuron 2017, 93, 1015–1034. [Google Scholar] [CrossRef] [Green Version]
  76. Yamamoto, S.; Kuramoto, K.; Wang, N.; Situ, X.; Priyadarshini, M.; Zhang, W.; Cordoba-Chacon, J.; Layden, B.T.; He, C. Autophagy differentially regulates insulin production and insulin sensitivity. Cell Rep. 2018, 23, 3286–3299. [Google Scholar] [CrossRef]
  77. Ma, D.; Li, S.; Molusky, M.M.; Lin, J.D. Circadian autophagy rhythm: A link between clock and metabolism? Trends Endocrinol. Metab. 2012, 23, 319–325. [Google Scholar] [CrossRef] [Green Version]
  78. Hansen, M.; Rubinsztein, D.C.; Walker, D.W. Autophagy as a promoter of longevity: Insights from model organisms. Nat. Rev. Mol. Cell Biol. 2018, 19, 579–593. [Google Scholar] [CrossRef]
  79. Qomaladewi, N.P.; Kim, M.Y.; Cho, J.Y. Rottlerin reduces cAMP/CREB-mediated melanogenesis via regulation of autophagy. Int. J. Mol. Sci. 2019, 20, 2081. [Google Scholar] [CrossRef] [Green Version]
  80. Amaravadi, R.K.; Kimmelman, A.C.; Debnath, J. Targeting autophagy in cancer: Recent advances and future directions. Cancer Discov. 2019, 9, 1167–1181. [Google Scholar] [CrossRef] [Green Version]
  81. Dyczynski, M.; Yu, Y.; Otrocka, M.; Parpal, S.; Braga, T.; Henley, A.B.; Zazzi, H.; Lerner, M.; Wennerberg, K.; Viklund, J.; et al. Targeting autophagy by small molecule inhibitors of vacuolar protein sorting 34 (Vps34) improves the sensitivity of breast cancer cells to Sunitinib. Cancer Lett. 2018, 435, 32–43. [Google Scholar] [CrossRef]
  82. Levine, B.; Mizushima, N.; Virgin, H.W. Autophagy in immunity and inflammation. Nature 2011, 469, 323–335. [Google Scholar] [CrossRef]
  83. Pahari, S.; Negi, S.; Aqdas, M.; Arnett, E.; Schlesinger, L.S.; Agrewala, J.N. Induction of autophagy through CLEC4E in combination with TLR4: An innovative strategy to restrict the survival of Mycobacterium tuberculosis. Autophagy 2019, 1–23. [Google Scholar] [CrossRef]
  84. Homma, K.; Suzuki, K.; Sugawara, H. The Autophagy Database: An all-inclusive information resource on autophagy that provides nourishment for research. Nucleic Acids Res. 2011, 39, D986–D990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Turei, D.; Foldvari-Nagy, L.; Fazekas, D.; Modos, D.; Kubisch, J.; Kadlecsik, T.; Demeter, A.; Lenti, K.; Csermely, P.; Vellai, T.; et al. Autophagy Regulatory Network—A systems-level bioinformatics resource for studying the mechanism and regulation of autophagy. Autophagy 2015, 11, 155–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Deng, W.; Ma, L.; Zhang, Y.; Zhou, J.; Wang, Y.; Liu, Z.; Xue, Y. THANATOS: An integrative data resource of proteins and post-translational modifications in the regulation of autophagy. Autophagy 2018, 14, 296–310. [Google Scholar] [CrossRef] [Green Version]
  87. Nanduri, R.; Kalra, R.; Bhagyaraj, E.; Chacko, A.P.; Ahuja, N.; Tiwari, D.; Kumar, S.; Jain, M.; Parkesh, R.; Gupta, P. AutophagySMDB: A curated database of small molecules that modulate protein targets regulating autophagy. Autophagy 2019, 15, 1280–1295. [Google Scholar] [CrossRef] [PubMed]
  88. Wang, N.N.; Dong, J.; Zhang, L.; Ouyang, D.; Cheng, Y.; Chen, A.F.; Lu, A.P.; Cao, D.S. HAMdb: A database of human autophagy modulators with specific pathway and disease information. J. Cheminform. 2018, 10, 34. [Google Scholar] [CrossRef] [PubMed]
  89. Csermely, P.; Korcsmaros, T.; Kiss, H.J.; London, G.; Nussinov, R. Structure and dynamics of molecular networks: A novel paradigm of drug discovery: A comprehensive review. Pharmacol. Ther. 2013, 138, 333–408. [Google Scholar] [CrossRef] [Green Version]
  90. Ozdemir, E.S.; Halakou, F.; Nussinov, R.; Gursoy, A.; Keskin, O. Methods for Discovering and Targeting Druggable Protein-Protein Interfaces and Their Application to Repurposing. Methods Mol. Biol. 2019, 1903, 1–21. [Google Scholar] [CrossRef]
  91. Gonen, M. Predicting drug-target interactions from chemical and genomic kernels using Bayesian matrix factorization. Bioinformatics 2012, 28, 2304–2310. [Google Scholar] [CrossRef] [PubMed]
  92. Gundersen, G.W.; Jones, M.R.; Rouillard, A.D.; Kou, Y.; Monteiro, C.D.; Feldmann, A.S.; Hu, K.S.; Ma’ayan, A. GEO2Enrichr: Browser extension and server app to extract gene sets from GEO and analyze them for biological functions. Bioinformatics 2015, 31, 3060–3062. [Google Scholar] [CrossRef] [PubMed]
  93. Gyori, B.M.; Bachman, J.A.; Subramanian, K.; Muhlich, J.L.; Galescu, L.; Sorger, P.K. From word models to executable models of signaling networks using automated assembly. Mol. Syst. Biol. 2017, 13, 954. [Google Scholar] [CrossRef] [PubMed]
  94. Pillich, R.T.; Chen, J.; Rynkov, V.; Welker, D.; Pratt, D. NDEx: A Community Resource for Sharing and Publishing of Biological Networks. Methods Mol. Biol. 2017, 1558, 271–301. [Google Scholar] [CrossRef] [PubMed]
  95. Gatliff, J.; East, D.; Crosby, J.; Abeti, R.; Harvey, R.; Craigen, W.; Parker, P.; Campanella, M. TSPO interacts with VDAC1 and triggers a ROS-mediated inhibition of mitochondrial quality control. Autophagy 2014, 10, 2279–2296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Roth, B.L.; Sheffler, D.J.; Kroeze, W.K. Magic shotguns versus magic bullets: Selectively non-selective drugs for mood disorders and schizophrenia. Nat. Rev. Drug Discov. 2004, 3, 353–359. [Google Scholar] [CrossRef] [Green Version]
  97. Ramsay, R.R.; Popovic-Nikolic, M.R.; Nikolic, K.; Uliassi, E.; Bolognesi, M.L. A perspective on multi-target drug discovery and design for complex diseases. Clin. Transl. Med. 2018, 7, 3. [Google Scholar] [CrossRef] [Green Version]
  98. Peters, J.U. Polypharmacology—Foe or friend? J. Med. Chem. 2013, 56, 8955–8971. [Google Scholar] [CrossRef]
  99. Rolf, M.G.; Curwen, J.O.; Veldman-Jones, M.; Eberlein, C.; Wang, J.; Harmer, A.; Hellawell, C.J.; Braddock, M. In vitro pharmacological profiling of R406 identifies molecular targets underlying the clinical effects of fostamatinib. Pharmacol. Res. Perspect. 2015, 3, e00175. [Google Scholar] [CrossRef]
  100. Selent, J.; Lopez, L.; Sanz, F.; Pastor, M. Multi-receptor binding profile of clozapine and olanzapine: A structural study based on the new beta2 adrenergic receptor template. ChemMedChem 2008, 3, 1194–1198. [Google Scholar] [CrossRef] [Green Version]
  101. Sridharan, S.; Jain, K.; Basu, A. Regulation of autophagy by kinases. Cancers (Basel) 2011, 3, 2630–2654. [Google Scholar] [CrossRef] [PubMed]
  102. Nussinov, R.; Zhang, M.; Tsai, C.J.; Jang, H. Calmodulin and IQGAP1 activation of PI3Kalpha and Akt in KRAS, HRAS and NRAS-driven cancers. Biochim. Biophys. Acta Mol. Basis Dis. 2018, 1864, 2304–2314. [Google Scholar] [CrossRef] [PubMed]
  103. Nussinov, R.; Zhang, M.; Tsai, C.J.; Liao, T.J.; Fushman, D.; Jang, H. Autoinhibition in Ras effectors Raf, PI3Kalpha, and RASSF5: A comprehensive review underscoring the challenges in pharmacological intervention. Biophys. Rev. 2018, 10, 1263–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Huang, G.; Zhang, F.; Ye, Q.; Wang, H. The circadian clock regulates autophagy directly through the nuclear hormone receptor Nr1d1/Rev-erbalpha and indirectly via Cebpb/(C/ebpbeta) in zebrafish. Autophagy 2016, 12, 1292–1309. [Google Scholar] [CrossRef] [Green Version]
  105. Torres-Quiroz, F.; Filteau, M.; Landry, C.R. Feedback regulation between autophagy and PKA. Autophagy 2015, 11, 1181–1183. [Google Scholar] [CrossRef] [Green Version]
  106. Chen, Y.; Klionsky, D.J. The regulation of autophagy—Unanswered questions. J. Cell Sci. 2011, 124, 161–170. [Google Scholar] [CrossRef] [Green Version]
  107. Liu, B.; Bhatt, D.; Oltvai, Z.N.; Greenberger, J.S.; Bahar, I. Significance of p53 dynamics in regulating apoptosis in response to ionizing radiation, and polypharmacological strategies. Sci. Rep. 2014, 4, 6245. [Google Scholar] [CrossRef] [Green Version]
  108. Liu, B.; Liu, Q.; Yang, L.; Palaniappan, S.K.; Bahar, I.; Thiagarajan, P.S.; Ding, J.L. Innate immune memory and homeostasis may be conferred through crosstalk between the TLR3 and TLR7 pathways. Sci. Signal. 2016, 9, ra70. [Google Scholar] [CrossRef] [Green Version]
  109. Liu, B.; Gyori, B.; Thiagarajan, P. Statistical model checking based analysis of biological betworks. Autom. Reason. Syst. Biol. Med. 2019. [Google Scholar] [CrossRef] [Green Version]
  110. Liu, B.; Thiagarajan, P.S. Modeling and analysis of biopathways dynamics. J. Bioinform. Comput. Biol. 2012, 10, 1231001. [Google Scholar] [CrossRef] [Green Version]
  111. Ueno, T.; Komatsu, M. Autophagy in the liver: Functions in health and disease. Nat. Rev. Gastroenterol. Hepatol. 2017, 14, 170–184. [Google Scholar] [CrossRef] [PubMed]
  112. Wang, M.M.; Feng, Y.S.; Yang, S.D.; Xing, Y.; Zhang, J.; Dong, F.; Zhang, F. The relationship between autophagy and brain plasticity in neurological diseases. Front. Cell. Neurosci. 2019, 13, 228. [Google Scholar] [CrossRef] [PubMed]
  113. Czaja, M.J.; Ding, W.X.; Donohue, T.M., Jr.; Friedman, S.L.; Kim, J.S.; Komatsu, M.; Lemasters, J.J.; Lemoine, A.; Lin, J.D.; Ou, J.H.; et al. Functions of autophagy in normal and diseased liver. Autophagy 2013, 9, 1131–1158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Schiattarella, G.G.; Hill, J.A. Therapeutic targeting of autophagy in cardiovascular disease. J. Mol. Cell. Cardiol. 2016, 95, 86–93. [Google Scholar] [CrossRef] [Green Version]
  115. Nixon, R.A. The role of autophagy in neurodegenerative disease. Nat. Med. 2013, 19, 983–997. [Google Scholar] [CrossRef]
  116. Garcia, M.L.; Cleveland, D.W. Going new places using an old MAP: Tau, microtubules and human neurodegenerative disease. Curr. Opin. Cell Biol. 2001, 13, 41–48. [Google Scholar] [CrossRef]
  117. Dubey, J.; Ratnakaran, N.; Koushika, S.P. Neurodegeneration and microtubule dynamics: Death by a thousand cuts. Front. Cell. Neurosci. 2015, 9, 343. [Google Scholar] [CrossRef] [Green Version]
  118. Kagan, V.E.; Mao, G.; Qu, F.; Angeli, J.P.; Doll, S.; Croix, C.S.; Dar, H.H.; Liu, B.; Tyurin, V.A.; Ritov, V.B.; et al. Oxidized arachidonic and adrenic PEs navigate cells to ferroptosis. Nat. Chem. Biol. 2017, 13, 81–90. [Google Scholar] [CrossRef]
  119. Wenzel, S.E.; Tyurina, Y.Y.; Zhao, J.; St Croix, C.M.; Dar, H.H.; Mao, G.; Tyurin, V.A.; Anthonymuthu, T.S.; Kapralov, A.A.; Amoscato, A.A.; et al. PEBP1 wardens ferroptosis by enabling lipoxygenase generation of lipid death signals. Cell 2017, 171, 628–641. [Google Scholar] [CrossRef]
  120. Kapralov, A.A.; Yang, Q.; Dar, H.H.; Tyurina, Y.Y.; Anthonymuthu, T.S.; Kim, R.; St Croix, C.M.; Mikulska-Ruminska, K.; Liu, B.; Shrivastava, I.H.; et al. Redox lipid reprogramming commands susceptibility of macrophages and microglia to ferroptotic death. Nat. Chem. Biol. 2020, 16, 278–290. [Google Scholar] [CrossRef]
  121. Steinman, J.; Epperly, M.; Hou, W.; Willis, J.; Wang, H.; Fisher, R.; Liu, B.; Bahar, I.; McCaw, T.; Kagan, V.; et al. Improved total-body irradiation survival by delivery of two radiation mitigators that target distinct cell death pathways. Radiat. Res. 2018, 189, 68–83. [Google Scholar] [CrossRef] [PubMed]
  122. Thermozier, S.; Hou, W.; Zhang, X.; Shields, D.; Fisher, R.; Bayir, H.; Kagam, V.; Yu, J.; Liu, B.; Bahar, I.; et al. Anti-ferroptosis drug enhances total body irradiation mitigation by drugs that block apoptosis and necroptosis. Radiat. Res. 2020, in press. [Google Scholar] [CrossRef]
  123. Lovric, M.; Molero, J.M.; Kern, R. PySpark and RDKit: Moving towards Big Data in Cheminformatics. Mol. Inform. 2019, 38, e1800082. [Google Scholar] [CrossRef] [PubMed]
  124. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  125. Benjamini, Y.; Drai, D.; Elmer, G.; Kafkafi, N.; Golani, I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125, 279–284. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Space of autophagy modulators and their targets. (a) On the left we show three groups of autophagy modulators, and on the right the corresponding targets, which show considerable overlap. The number of targets in each subset are shown by the labels in the Venn diagrams on the right. The numbers on the linkers connecting the drugs to the targets indicate the number of targets associated with each category of autophagy modulators; (b) Distribution of the degrees (number of associated targets) for each category of autophagy modulators.
Figure 1. Space of autophagy modulators and their targets. (a) On the left we show three groups of autophagy modulators, and on the right the corresponding targets, which show considerable overlap. The number of targets in each subset are shown by the labels in the Venn diagrams on the right. The numbers on the linkers connecting the drugs to the targets indicate the number of targets associated with each category of autophagy modulators; (b) Distribution of the degrees (number of associated targets) for each category of autophagy modulators.
Ijms 21 02855 g001
Figure 2. Promiscuity of autophagy modulators. (a) The number of known and predicted targets of activators; (b) The number of known and predicted targets of inhibitors; (c) The number of known and predicted targets of dual-modulators. The drugs are ranked in a descending order of the number of known targets. The contributions of the known and predicted drugs to the total number of targets are shown by the dark and light colors as indicated by the labels. See also Figure S3 for additional data on activators. Insets: Percent distribution of modulators vis-à-vis their promiscuity. The promiscuity is quantified by the degree of each modulator, i.e., the number of targets (abscissas) associated with the modulator.
Figure 2. Promiscuity of autophagy modulators. (a) The number of known and predicted targets of activators; (b) The number of known and predicted targets of inhibitors; (c) The number of known and predicted targets of dual-modulators. The drugs are ranked in a descending order of the number of known targets. The contributions of the known and predicted drugs to the total number of targets are shown by the dark and light colors as indicated by the labels. See also Figure S3 for additional data on activators. Insets: Percent distribution of modulators vis-à-vis their promiscuity. The promiscuity is quantified by the degree of each modulator, i.e., the number of targets (abscissas) associated with the modulator.
Ijms 21 02855 g002
Figure 3. An overview of proteins targeted by autophagy modulators. (a) Frequent targets of autophagy activators (orange), inhibitors (violet), and dual-modulators (blue), ranked by the number of the total known targets in descending order. In light shades are those predicted by our probabilistic matrix factorization (PMF) method. (b) Distribution of targeted proteins by protein families; (c) Tissues enriched in targets of autophagy modulators.
Figure 3. An overview of proteins targeted by autophagy modulators. (a) Frequent targets of autophagy activators (orange), inhibitors (violet), and dual-modulators (blue), ranked by the number of the total known targets in descending order. In light shades are those predicted by our probabilistic matrix factorization (PMF) method. (b) Distribution of targeted proteins by protein families; (c) Tissues enriched in targets of autophagy modulators.
Ijms 21 02855 g003
Figure 4. Gene ontology (GO) annotation and pathway enrichment analysis results based on the targets of autophagy modulators. (a) Enriched GO annotation terms; (b) Enriched KEGG pathways; (c) most common pathways in which autophagy targets take part. The distribution of pathway subclasses targeted by autophagy modulators are shown for each category of autophagy modulators. Note the dominance of signal transduction (red) and endocrine (orange) systems in all categories, the important association of inhibitors with cancers and, to some extent, infectious diseases, and the preferential use of activators and dual modulators in treating neural diseases.
Figure 4. Gene ontology (GO) annotation and pathway enrichment analysis results based on the targets of autophagy modulators. (a) Enriched GO annotation terms; (b) Enriched KEGG pathways; (c) most common pathways in which autophagy targets take part. The distribution of pathway subclasses targeted by autophagy modulators are shown for each category of autophagy modulators. Note the dominance of signal transduction (red) and endocrine (orange) systems in all categories, the important association of inhibitors with cancers and, to some extent, infectious diseases, and the preferential use of activators and dual modulators in treating neural diseases.
Ijms 21 02855 g004
Figure 5. Signal transduction pathways most frequently affected by autophagy modulators and the most promiscuous targets belonging to these pathways. (a) Schematic diagram of the seven autophagy-related (ATG)-signal transduction pathways and highly promiscuous targets participating in those pathways (indicated by the gray arrow pointing from the pathway to the protein). The target proteins are shown in light violet circles, except for a pivotal target (PKA) (shown in yellow). Bold face black arrows from X to Y means that X regulates Y. (b) Target-pathway mapping for the most prominent ATG-signal transduction pathways (abscissas) and ATG-signaling targets (ordinate). The color blocks associated with each protein indicate that it is targeted by autophagy activators (green), inhibitors (red), and dual-modulators (blue). Note that mTOR is also targeted by all three groups of modulators.
Figure 5. Signal transduction pathways most frequently affected by autophagy modulators and the most promiscuous targets belonging to these pathways. (a) Schematic diagram of the seven autophagy-related (ATG)-signal transduction pathways and highly promiscuous targets participating in those pathways (indicated by the gray arrow pointing from the pathway to the protein). The target proteins are shown in light violet circles, except for a pivotal target (PKA) (shown in yellow). Bold face black arrows from X to Y means that X regulates Y. (b) Target-pathway mapping for the most prominent ATG-signal transduction pathways (abscissas) and ATG-signaling targets (ordinate). The color blocks associated with each protein indicate that it is targeted by autophagy activators (green), inhibitors (red), and dual-modulators (blue). Note that mTOR is also targeted by all three groups of modulators.
Ijms 21 02855 g005
Figure 6. A network of pathways connecting intracellular signal transduction events to autophagy modulation. (a) A unified signaling network mediates the effects of autophagy modulators on signaling transduction; (b) The regulatory dependency of signaling transduction modules. Color coding: calcium signaling (red), cAMP signaling (violet), cGMP-PKG signaling (orange), Rap1/Ras/MAPK signaling (blue), PI3K-AKT signaling (green), core autophagy proteins (yellow), and transcription factors (cyan). The core autophagy proteins were grouped (indicated by braces) by their roles in regulating three phases of autophagy: formation of autophagosome, fusion with lysosome, and lysosomal degradation. The key ATG-signal transduction proteins shown in Figure 5a (except MEK) are highlighted by the yellow ellipses. Selected modulators of key targets are listed in gray font. See the code for different types of actions indicated by arrow ends on top.
Figure 6. A network of pathways connecting intracellular signal transduction events to autophagy modulation. (a) A unified signaling network mediates the effects of autophagy modulators on signaling transduction; (b) The regulatory dependency of signaling transduction modules. Color coding: calcium signaling (red), cAMP signaling (violet), cGMP-PKG signaling (orange), Rap1/Ras/MAPK signaling (blue), PI3K-AKT signaling (green), core autophagy proteins (yellow), and transcription factors (cyan). The core autophagy proteins were grouped (indicated by braces) by their roles in regulating three phases of autophagy: formation of autophagosome, fusion with lysosome, and lysosomal degradation. The key ATG-signal transduction proteins shown in Figure 5a (except MEK) are highlighted by the yellow ellipses. Selected modulators of key targets are listed in gray font. See the code for different types of actions indicated by arrow ends on top.
Ijms 21 02855 g006
Table 1. Experimental validation of predicted drug–target interactions.
Table 1. Experimental validation of predicted drug–target interactions.
DrugPredicted Target 1Confidence ScoreBinding Affinity Ki (nM)Reference
olanzapineHTR1F0.6298310[52]
fostamatinibFGFR40.6267350[50]
fluphenazineDRD30.65260.11[52]
fluspirileneDRD30.60850.40[52]
thioridazineDRD30.72031.5[52]
sertindoleDRD30.82652.5[51]
trifluoperazineDRD30.60034.2[51]
prazosinADRA2C0.871310.7[52]
thioridazineHRH10.672216[52]
chlorpromazineCHRM50.841842[52]
sorafenibPDGFRA0.628962[50]
pimozideHTR1A0.604488[52]
verapamilHTR2A0.6419140[52]
fluoxetineHTR2A0.7219148[52]
maprotilineDRD10.8489402[52]
propranololHTR2C0.6593574[52]
nortriptylineHRH20.7008645[52]
acetylcholineCHRM50.8279800[52]
dasatinibFGFR20.61211400[50]
imatinibFLT30.62066300[50]
1 Full names of targets are given in Table S3. See other experimentally validated predictions in Table S3 (rows with binding affinity information in the last column).

Share and Cite

MDPI and ACS Style

Shi, Q.; Pei, F.; Silverman, G.A.; Pak, S.C.; Perlmutter, D.H.; Liu, B.; Bahar, I. Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis. Int. J. Mol. Sci. 2020, 21, 2855. https://doi.org/10.3390/ijms21082855

AMA Style

Shi Q, Pei F, Silverman GA, Pak SC, Perlmutter DH, Liu B, Bahar I. Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis. International Journal of Molecular Sciences. 2020; 21(8):2855. https://doi.org/10.3390/ijms21082855

Chicago/Turabian Style

Shi, Qingya, Fen Pei, Gary A. Silverman, Stephen C. Pak, David H. Perlmutter, Bing Liu, and Ivet Bahar. 2020. "Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis" International Journal of Molecular Sciences 21, no. 8: 2855. https://doi.org/10.3390/ijms21082855

APA Style

Shi, Q., Pei, F., Silverman, G. A., Pak, S. C., Perlmutter, D. H., Liu, B., & Bahar, I. (2020). Mechanisms of Action of Autophagy Modulators Dissected by Quantitative Systems Pharmacology Analysis. International Journal of Molecular Sciences, 21(8), 2855. https://doi.org/10.3390/ijms21082855

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop