Biological Network Analysis of PPAR and Related Signaling PathwaysView this Special Issue
Role of PPARs in Progression of Anxiety: Literature Analysis and Signaling Pathways Reconstruction
Peroxisome proliferator-activated receptor (PPAR) group includes three isoforms encoded by PPARG, PPARA, and PPARD genes. High concentrations of PPARs are found in parts of the brain linked to anxiety development, including hippocampus and amygdala. Among three PPAR isoforms, PPARG demonstrates the highest expression in CNS, where it can be found in neurons, astrocytes, and glial cells. Herein, the highest PPARG expression occurs in amygdala. However, little is known considering possible connections between PPARs and anxiety behavior. We reviewed possible connections between PPARs and anxiety. We used the Pathway Studio software (Elsevier). Signal pathways were created according to previously developed algorithms. SNEA was performed in Pathway Studio. Current study revealed 14 PPAR-regulated proteins linked to anxiety. Possible mechanism of PPAR involvement in neuroinflammation protection is proposed. Signal pathway reconstruction and reviewing aimed to reveal possible connection between PPARG and CCK-ergic system was conducted. Said analysis revealed that PPARG-dependent regulation of MME and ACE peptidase expression may affect levels of nonhydrolysed, i.e., active CCK-4. Impairments in PPARG regulation and following MME and ACE peptidase expression impairments in amygdala may be the possible mechanism leading to pathological anxiety development, with brain CCK-4 accumulation being a key link. Literature data analysis and signal pathway reconstruction and reviewing revealed two possible mechanisms of peroxisome proliferator-activated receptors involvement in pathological anxiety: (1) cytokine expression and neuroinflammation mechanism and (2) regulation of peptidases targeted to anxiety-associated neuropeptides, primarily CCK-4, mechanism.
1.1. Anxiety and Anxiety Disorders
Anxiety disorders (including generalized anxiety disorders or panic disorders) are the most widespread mental diseases which are at the same time difficult to treat [1, 2]. Main characteristic of panic disorder is presence of repetitive sudden panic attacks . According to large scale surveys, percentage of population suffering from the anxiety disease throughout lifetime is up to 33.7% [4, 5].
Many researchers report that anxiety disorders cause even more severe decrease in patient’s quality of life and psychosocial functions than other chronic diseases including diabetes, cardiovascular diseases, and lung diseases [6–8]. Both environmental factors and genetic factors are believed to be involved in the development of panic disorders [9, 10].
Last decade studies showed that anxiety and anxiety disorder are associated with amygdala functioning and various types of its dysfunction that lead to a decrease in its activity [11–15]. Hypothalamus is also often associated with anxiety [16, 17]. Thus, it was shown that the anorexigenic neuropeptide CCK4 is able to directly interact with the hypothalamus .
Most modern studies of the molecular basis underlying F40-F48 series (ICD-10) mental disorders focus on polymorphisms in genes encoding the main neurotransmitter system proteins, i.e., catecholamine and GABAergic . However, it is already clear that they are not the primary link in the fine regulation of the formation, severity, and direction of anxious emotions.
1.2. Peroxisome Proliferator-Activated Receptors
Peroxisome proliferator-activated receptor (PPAR) subfamily belongs to nuclear receptor family. Three isoforms encoded by individual genes are known: PPARγ, PPARα, and PPARδ. PPARs are ligand-dependent transcription factors that regulate target gene expression by binding to specific peroxisome proliferator response elements (PPREs) in the enhancer sites of regulated genes. Each receptor binds to its PPRE as a heterodimer with a retinoid X receptor (RXR). Upon binding of the agonist, PPAR conformation changes and stabilizes, after which transcriptional coactivators contribute to activation of target genes .
The PPARs possess the canonical domain structure common to other nuclear receptor family members, including the amino-terminal AF-1 transactivation domain, followed by a DNA-binding domain and a dimerization and ligand binding domain with a ligand-dependent transactivation function AF-2 located at the carboxy-terminal region .
PPARs regulate expression of genes actively involved in lipid and carbohydrate metabolism, vascular homeostasis, tissue repair, cell proliferation and differentiation, and sexual dimorphism .
PPARs are expressed in almost all mammalian tissues and organs. The expression patterns of PPARα, PPARβ/δ, and PPARγ differ, although intersections do occur. A high expression level of all PPARs is observed in tissues with active fatty acid metabolism. PPARβ/δ is constitutively expressed in almost all tissues [22, 23]. PPARγ among all 3 isoforms has the highest expression in the nervous system, where it was found in neurons, astrocytes, and glial cells . Moreover, all three PPA receptors are expressed both in amygdala and in the hypothalamus . PPAR association with various diseases is also shown . However, the association with anxiety is poorly understood.
This study for the first time analyzes the relationship between PPARs and anxiety. Our results will allow to look at molecular-genetic basis of anxiety disorders pathogenesis, underlying mechanisms and related problems from a new angle.
2. Materials and Methods
We used the Pathway Studio® 9 desktop software with ResNet® 14 database and web version of the Pathway Studio software (https://mammalcedfx.pathwaystudio.com) (Elsevier). Additional search of information was performed by using PubMed (http://www.ncbi.nlm.nih.gov/pubmed/), TargetInsights (https://demo.elseviertextmining.com/), and Google Scholar (https://scholar.google.ru/). Signal pathways were created according to previously developed algorithms .
Search algorithm and workflow scheme are presented in Figure 1. Detailed description is given along in Results for better understanding.
3. Results and Discussion
3.1. A Search for Common PPARs and Anxiety Targets
During the first step of our work, we used text-mining and signaling pathway analysis for revealing possible role of PPARs in the development of anxiety.
We used the Pathway Studio® 9 desktop software with ResNet® 14 database and web version of the Pathway Studio software (https://mammalcedfx.pathwaystudio.com) (Elsevier).
The search algorithm for anxiety-protein-PPAR links was as follows: (1)Searching proteins linked to anxiety through concentration change. Linkage type: QuantitativeChange was used in Pathway Studio. This type allows to search for proteins with reported concentration alterations in patients with anxiety. Total of 304 proteins were found(2)304 proteins discovered on the previous step were taken into further analysis by searching linkage between protein and each of the three PPARs independently. We had chosen to search for PPAR connections separately to avoid too complicated pathways. Linkage types: Expression and PromoterBinding were used in Pathway Studio. These two types imply direct PPAR effect on protein concentration through promoter binding or other direct interaction which makes them the most relevant(3)On the next step, all the supporting references were manually revised; at first by reviewing Pathway Studio text-mined “Sentences” section, and after that, if necessary, by studying the whole text for additional details(4)This analysis resulted in three schemes which include all the proteins linked to both anxiety and PPARs. Each scheme was representing one PPAR isoform(5)Three schemes created on the previous step were combined into one for better data representation and reducing number of figures
For describing proteins and their role, search engines highlighted in Materials and Methods were used.
3.2. Common Target Descriptions and Possible Underlying Mechanisms
Only cytokines (IFNG, IL6, and TNF) and vascular endothelial growth factor A (VEGFA) are regulated by all three PPARs. In most cases, connection to PPARα is present. In one case, HSD11B2, only connection to PPARD is present.
Thus, 14 PPAR-regulated proteins with altered concentration in anxiety are known:
3.2.1. BDNF (Brain-Derived Neurotrophic Factor)
Protein’s main function is nerve growth and neuron homeostasis . Both BDNF polymorphisms and concentration alterations are linked to variety of mental and neurological diseases . In our case, reduced BDNF product concentration is observed in both animal models (, see supplementary (available here)) and patients (, see supplementary). The BDNF expression is positively regulated by PPARA (, see supplementary) and PPARG (, see supplementary). If PPARs-independent BDNF expression regulation pathways are unaltered, lower PPARA and PPARG levels are to be expected as well.
3.2.2. CREB1 (CAMP Responsive Element Binding Protein 1)
This is a transcription factor which induces target gene transcription in response to hormonal stimulation via the cAMP pathway. This protein is involved in many cellular processes and is expressed in all tissues (https://www.ncbi.nlm.nih.gov/gene/1385). CREB1 activity is important for brain neuron development and maintenance [33, 34]. Thus, its linkage to mental diseases’ pathogenesis is widely studied [35, 36]. The CREB1 expression is known to be lowered in animal models of anxiety [37, 38]. PPARA and PPARG activate CREB1 transcription by binding to its promoter. Exact in-promoter targets of PPARA are known [39, 40], which points out its important role on CREB1 regulation. PPARG was shown to be able to bind to CREB1 promoter bearing certain nucleotide sequence and block its expression . However, other mechanisms leading to CREB1 activation along with PPARG activation are known, i.e., PKCA (protein kinase C, alpha) activation with further CREB1 activation . Interestingly, CREB1 activation via cAMP-PKCA signaling provokes binding of CREB1 to a cAMP responsive element-like site in PPARG gene promoter region  which allows to speculate of mutual regulation of these factors.
3.2.3. CRP (C-Reactive Protein)
This protein is involved in several host defence-related functions based on its ability to recognize foreign pathogens and damaged host cells and initiate their elimination by interacting with humoral and cellular effector systems in the blood (https://www.ncbi.nlm.nih.gov/gene/1401). CRP linkage to anxiety was shown in several studies. Patients with anxiety were shown to have increased CRP levels [44–46], wherein all authors point out the alterations in patients’ immune status. This can be explained by possible interactions between the central nervous system and the immune system in neuropsychiatric disorders. CRP regulation via PPARs is negative in both PPARA ([47, 48], see supplementary) and PPARG ([49, 50], see supplementary) cases. There is no information on PPARD-CRP interactions. Thus, in case of anxiety, we can observe an increase of CRP levels despite PPARA and PPARG blocking it. The reasons may lie in either other CRP activation signal pathways or tissue effects insofar CRP is expressed in the liver. This is the only protein in this list which lacks expression brain regions of interest.
FOS proteins are members of the transcription factor complex AP-1. FOS proteins were implicated as cell proliferation, differentiation, and transformation regulators. FOS gene is expressed in almost every tissue (https://www.ncbi.nlm.nih.gov/gene/2353). One of the most thorough research of this gene is focused on neuronal plasticity [51, 52]. Animal models of anxiety, fear, and depression show increased FOS in the brain areas linked to the formation of these behavioral reactions in animals [53–55]. Negative effect of PPARA activation on the FOS expression is known [56, 57]. It is also known that PPARG regulates the FOS expression [58–61], but there is no univocal opinion wherever this regulation is positive or negative. It is most likely that several signal pathways involving PPARG and FOS are present. Most of the PPARG-FOS interaction studies are using osteogenesis models, which makes it impossible to speculate about PPARs-FOS-anxiety linkage. This linkage is the least evidence-supported among this list.
The corticosteroid 11-beta-dehydrogenase is a microsomal enzyme complex responsible for the interconversion of glucocorticoid cortisol and its inactive metabolite cortisone; this conversion prevents illicit activation of the mineralocorticoid receptor (https://www.ncbi.nlm.nih.gov/gene/3291). Cortisol is intricately linked to anxiety as increases in its concentration correlates with anxiety and its symptoms [62, 63]. Insofar, as HSD11B2 enzyme is directly involved in processes related to cortisol, its role is pathogenesis of diseases characterized by increased cortisol levels is important . Animal models with HSD11B2 gene knockout have inborn predisposition to increased anxiety . Thus, reduced HSD11B2 activity leads to increased cortisol and anxiety phenotype development. Main organs expressing HSD11B2 are the kidneys, intestines, and salivary glands. Only three of PPAR receptors regulate the HSD11B2 expression wherein this regulation is negative. PPARD is known to bind to HSD11B2 promotor in trophoblasts of the placenta and inhibits its expression [66–68]. We have no knowledge of studies describing PPARD and HSD11B2 interactions in other tissues. However, as far as PPARD is expressed in all tissues, both genes are expressed in the kidney which is the most interesting organ from the point of cortisol accumulation. Thus, in case of HSD11B2, one of peroxisome proliferator-activated receptors (PPARD) may be linked to one of the clinical anxiety manifestations, i.e., increased cortisol.
Interferon gamma is soluble cytokine that is a member of the type II interferon class, secreted by cells of both the innate and adaptive immune systems (https://www.ncbi.nlm.nih.gov/gene/3458). Its expression is limited by immune system cells. Interferon gamma may modulate anxiety and depressive states via its role in brain plasticity . Other study suggests that decreased INFG in patients with anxiety disorder is not a cause but a result of the disease . However, all PPARs regulate its expression negatively via INFG promotor binding and blocking its transcriptional activity ([71–75], see supplementary). Thus, link between PPARs and INFG regulation is unequivocal, but determining if decreased INFG concentration is a disease’s cause, or a consequence is a subject for further studies. It is possible that decrease of INFG concentration is a result of prolonged anxiety.
3.2.7. IL6 (Interleukin 6)
This cytokine plays role in inflammation and maturation of B cells; encoded protein was shown to be an endogenous pyrogen capable of inducing fever in people with autoimmune diseases or infections (https://www.ncbi.nlm.nih.gov/gene/3569). Increase of the IL6 expression in patients with anxiety/depression is well known [76–80]. However, some studies show decrease of IL6 levels in states characterized by pathological anxiety [81, 82]. Thus, IL6 anxiety link is possible but requires further studies. All three PPARs negatively regulate IL6 ([83–85], see supplementary). Further understanding PPARs and IL6 interaction requires evaluation of IL6 role in anxiety.
Leptin is secreted by white adipocytes into the circulation and plays a major role in the regulation of energy homeostasis; it binds to the leptin receptor in the brain, which activates downstream signaling pathways that inhibit feeding and promote energy expenditure (https://www.ncbi.nlm.nih.gov/gene/3952). Leptin concentration is increased in patients with emotional anxiety and reaches its peak in conditions of moderate anxiety . Other studies also show high levels of leptin which correlates with anxiety levels [87–89]. Leptin administration resulted in dose-dependent anxiety decrease [86, 90]. This data points out the positive link between leptin and anxiety. Among three PPARs, two receptors regulate LEP gene expression, i.e., delta and gamma. PPARD binds to the LEP promotor and inhibits its expression [91–93]. Thus, one of the mechanisms of leptin increase in anxiety may be linked to PPARD and PPARG decrease.
Neuropeptide Y is a neuropeptide that is widely expressed in the central nervous system and influences many physiological processes, including cortical excitability, stress response, food intake, circadian rhythms, and cardiovascular function (https://www.ncbi.nlm.nih.gov/gene/4852). Neuropeptide Y deficiency is significantly linked to anxiety development in all animals including fish and human ([94–97], see supplementary). The NPY expression is linked to only one PPAR—gamma; its positive effect on the NPY expression is known in arcuate hypothalamus . We can speculate that NPY decrease may be linked to alterations in PPARG activation/functioning.
Glucocorticoid receptor can function both as a transcription factor that binds to glucocorticoid response elements in the promoters of glucocorticoid responsive genes to activate their transcription and as a regulator of other transcription factors and involved in inflammatory responses, cellular proliferation, and differentiation in target tissues (https://www.ncbi.nlm.nih.gov/gene/2908). Increase of the NR3C1 expression was linked to anxiety increases in rats ; antagonist administration resulted in reduction of anxiety-like behavior in rats . Increased glucocorticoid leads to increased anxiety via activation of receptor NR3C1 [101, 102]. Among all peroxisome proliferator-activated receptors, only PPARA is linked to the NR3C1 expression regulation as it reduces it [103–105]. Wherein mutual gene regulation is known as NR3C1 is capable of activating the PPARA expression . Thus, PPARA is capable of regulating NR3C1-dependent anxiety.
Prolactin is a hormonal growth regulator for many tissues, including cells of the immune system, essential for lactation (https://www.ncbi.nlm.nih.gov/gene/5617). An increase in prolactin levels is associated with anxiety in women during lactation , as well as in paratroopers before jumping . In roman low-avoidance rats with an increased level of anxiety, an increased level of the PRL gene expression in the amygdala was shown . PPARA is the only PPAR able of activating the PRL expression [110–112]. It is also noted that gene is not necessarily activated via PPARA binding to PRL promotor [111, 113]. Thus, one of the mechanisms leading to the increased PRL in anxiety may be explained due to PPARA activity.
Tumor necrosis factor is a multifunctional proinflammatory cytokine that belongs to the tumor necrosis factor superfamily; it is mainly secreted by macrophages and involved in the regulation of a wide spectrum of biological processes including cell proliferation, differentiation, apoptosis, lipid metabolism, and coagulation (https://www.ncbi.nlm.nih.gov/gene/7124). Serum TNF levels are elevated in patients with anxiety symptoms , as well as in patients with generalized anxiety disorder . TNF knockout mice showed a low level of anxiety . Cited works indicate a reliable association of TNF elevation with anxiety. All three PPARs block the expression of TNF, as well as the rest of the cytokines identified in our work ([92, 117–121], see supplementary for additional links). Thus, the explanation of increased TNF in anxiety patients lays either in the malfunction/decreased activity of PPARs or in the activation of another mechanism for regulating the TNF expression.
Translocator protein is a key factor in the flow of cholesterol into mitochondria to permit the initiation of steroid hormone synthesis and interacts with some benzodiazepines (https://www.ncbi.nlm.nih.gov/gene/706). The level of TSPO was significantly reduced in both patients with anxiety and anxiety mice models [122–125] and increased after treatment . Since the benzodiazepines are used in the treatment of anxiety and anxiety disorder , the association of their molecular target TSPO with anxiety is beyond doubt. The TSPO expression is negatively regulated by PPARA ([128, 129], see supplementary). It can be assumed that PPARA activity may adversely affect the development of anxiety due to a decrease in the TSPO expression.
3.2.14. VEGFA (Vascular Endothelial Growth Factor A)
It is a growth factor which induces proliferation and migration of vascular endothelial cells; essential for angiogenesis (https://www.ncbi.nlm.nih.gov/gene/7422). A decrease in VEGFA concentration in patients with anxiety [130, 131] and a high VEGFA level accompanying low anxiety  were reported. The VEGFA expression is regulated by all PPARs: delta and gamma bind to the promoter, activating the translation (, see supplementary), alpha, on the contrary, blocks VEGFA promoter, reducing expression (, see supplementary). Thus, activation and blockade of the VEGFA gene expression are also regulated by transcription factor PPARs. Since the VEGFA expression is reduced in patients with increased anxiety, it can be assumed that in this case, both PPARA activation and a decrease in PPAR delta and gamma activity may occur. Table 2 summarizes above said.
In our analysis, two groups of proteins are the most presented: (i)Inflammatory cytokines—three members (IL6, INFG, and TNF)(ii)Transcription factors—three members (CREB1, FOS, and NR3C1)
Interestingly, all PPARs are linked to the inflammatory cytokine expression, while only PPARA and PPARG are involved linked to transcription factor expression. Transcriptional targets of PPARD are involved in immune response and body homeostasis (HSD11B2, LEP, and VEGFA). It expresses in all brain tissues at the same level , which indicates the probable absence of its participation in the development of most behavioral reactions, including anxiety.
3.3. PPARs, Common Target, and Amygdala
Proteins, identified in the previous step, were reviewed for expression in amygdaloid: (1)Amygdaloid was added to the list, and all the amygdaloid-protein links were extracted(2)All the links except for CellExpression type were excluded(3)All the supporting references were manually revised at first by reviewing Pathway Studio text-mined “Sentences” section, and after that, if necessary, by studying the whole text for additional details
All detected proteins, except for CRP (C-reactive protein), are expressed in amygdala (Figure 3).
Information about the environment is acquired via sensory organs and is transferred to the thalamus nuclei of the limbic system and then to the cortical sections of the sensory analyzer (auditory, visual, tactile cortex): (i)The limbic system responds to the image that the brain has perceived and recognized. In particular, amygdala is responsible for the defensive reaction, fear, and aggression(ii)The amygdala sends a signal to the prefrontal cortex, which evaluates the situation. Its main function in this cascade is to develop a rescue plan in a situation of perceived danger. This path functions inappropriately in case of phobias, which leads to the development of a sense of fear for the stimuli that are not harmful. The amygdala is a key link in the anxiety formation; it is known that the groups of cells of the amygdala are activated when there is fear or aggression. The central nucleus of the amygdala has direct connections with the hypothalamus and brain stem—areas also responsible for fear.
PPARs and retinoid X receptor are expressed in the central nervous system. Delta is widely expressed in all parts of the brain, while alpha and gamma show selective expression: (i)Gamma is not expressed in the structures of the olfactory bulb, in the part of the olfactory cortex, part of the neocortex, some structures of the thalamus, nuclei of the solitary tract, dorsal motor nuclei of the vagus, and Purkinje cells(ii)Alpha is not expressed in the hypothalamus 
All 3 isoforms are actively expressed in the basal ganglia to which amygdala belongs [24, 136, 137]. Thus, PPARs in high concentration are located in the brain areas involved in the formation of anxiety, including those that are widely represented in the hippocampus and amygdala.
3.4. Enrichment Analysis
To further elucidate role of PPARs in anxiety, we conducted Sub-Network Enrichment Analysis (SNEA) using the Pathway Studio® software using list of genes obtained during the first stage of our study. We found considerable number of pathways/gene sets in almost every search subgroup (SNEA: for compounds regulators (branch-drugs), for compounds regulators, for phenotypes and processes (all), with anatomy, with expression regulators) for every gene associated with anxiety disorder.
Table 3 lists top 10 of 100 substances obtained after the SNEA Compounds Regulators Enriched analysis. Compound regulators were chosen as the most representative search subgroup. All these substances (ginsenoside, epinephrine, corticosterone, norepinephrine, Li+, serotonin, n-3 polyunsaturated fatty acid, dehydroepiandrosterone), either have an antistress, antidepressant, neuroprotective, and neurogenic effects or regulate blood pressure, carbohydrate, and lipid metabolism (according to https://pubchem.ncbi.nlm.nih.gov/). Some of them are found in cellular membranes and blood vessel walls. Despite the low Jaccard similarity, which may be explained by relatively low number of searched genes (17), according to enrichment analysis, good ratio of 9/10 of Gene Set Seed related to neurological regulation is observed. This means that 9 substances out of 10 are related to neurological regulation through underlying pathways. This result points out the significant role these genes play in pathogenesis of psychoneurological diseases.
In this set of 10 substances, genes BDNF, FOS, VEGFA, IFNG, IL6, TNF, and PPAR (A, D, G) are present in 100% of underlying pathways. PPARG is present in 80% of them, while PPARA in 70% cases and PPARD in 50%. This tendency is persists in full list of 100 substances obtained after SNEA Compounds Regulators Enriched analysis with BDNF presence in 77%; FOS, 76%; VEGFA, 82%; IFNG, 76%; IL6, 96%; TNF, 96%; and PPAR (A, D, G), 78%. Among PPARs, PPARA is present in 52% and PPARD in 27%; most presented is PPARG in 71%, wherein only 18% out of 100 substances lack any PPARs in their underlying pathway.
This result is an indirect evidence of PPAR role in the anxiety development.
3.5. Neuroinflammation, PPARs, and Anxiety
Recently, the theory of neuroinflammation as the cause of anxiety development was discussed [138–141]. We performed a signal pathway analysis which revealed possible mechanism explaining the role of PPARs in the development of neuroinflammation (Figure 4).
Signal pathway schemes were designed manually using data available in the ResNet database describing links between objects. Annotated schemes available in the Pathway Studio software were used as well.
We suggest that ligand-activated PPARs block the activity of the NF-kappa B family of transcription factors [142, 143]. Coactivators of PPARs-retinoid-X receptor subfamily  and PPARG coactivator 1 alpha (PPARGC1A)  are also involved in this process. NF-kappa B activates the transcription of number of proinflammatory cytokines, including those that change their expression during anxiety, i.e., IL6 , TNF , and IFNG . Thus, PPARs normally block the development of neuroinflammation [149, 150].
3.6. Peroxisome Proliferator-Activated Receptor Gamma (PPARG) and Cholecystokininergic Systems in Anxiety
In recent years, PPARG was shown to be involved in the development of pathological anxiety [151–153]. This surge of research was initiated by Domi et al., whose aim was to study the role of PPARG in the regulation of GABAergic transmission . The work revealed the role of PPARG in the regulation of mood disorders, indicating that weakened signal transmission can contribute to exacerbation of anxiety and the negative effects of stress. Authors suggest that activation of PPARG may be useful in the treatment of psychiatric conditions related to stress and, in particular, anxiety disorders .
Cholecystokinin (CCK) is a neuropeptide which may be found in high concentrations throughout the central nervous system, where it is involved in numerous physiological functions . The role of CCK, especially its smallest functional peptide, CCK-4, in the induction and maintenance of anxiety and major depression is well known [160–162]. The increase in CCK-4 is associated with loss of motivation, anxiety, and panic attacks [159, 163, 164].
The amygdala is firmly connected with the cholecystokininergic system. It was shown that CCK-4 is synthesized and released into amygdala, hippocampal formation, and cerebral cortex [165, 166]. Both CCK-4 receptors are expressed in the amygdala: CCKAR [167, 168] and CCKBR [168, 169].
In this part of the work, we carried out an analysis of signaling pathways in order to identify a possible connection between the cholecystokininergic system and PPARG.
We used the Pathway Studio® 9 desktop software with the ResNet®14 database and web version of the Pathway Studio software (https://mammalcedfx.pathwaystudio.com) (Elsevier). Signal pathways were created according to the algorithms we developed . The results are presented in Figure 5 and are described below.
CCK-4 hydrolysis (Trp-Met-Asp-Phe = WMDF) may be carried out by two enzymes: membrane metalloendopeptidase (MME, also known as NEP, EC 220.127.116.11) and angiotensin I converting enzyme (ACE, EC 18.104.22.168) . MME and ACE cut 2 terminal amino acids (DF) from CCK-4 and CCK-8 which leads to forming of unfunctional peptides (or dipeptides WM and DF in case of CCK-4) [174, 175].
Thus, PPARG-related MME and ACE expression regulation affect levels of active nonhydrolysed CCK-4. Dysregulation of PPARG and the following alteration of peptidase expression in amygdala may be a possible mechanism of pathological anxiety development, with a CCK-4 accumulation as a main cause. This may explain Domi et al. results .
We examined possible associations of PPARs with anxiety. An analysis of the data and signaling pathways available in the literature suggests two mechanisms for the participation of peroxisome proliferator-activated receptors in the formation of anxiety: (1) expression of cytokines and neuroinflammation and (2) regulation of the expression of peptidases, the targets of which are neuropeptides associated with anxiety (CCK-4 in the first place). Further research in these areas will help to better understand the role of PPARs in the development of anxiety.
All supplemental materials (Excel file) are available to download from ResearchGate resource by the link https://www.researchgate.net/publication/340418909_Supplemental_Materials_Role_of_PPARs_in_progression_of_anxiety. Additional data is available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
The reported study was funded by RFBR, project number 19-34-51045.
The supplement file was provided for this study (an excel file), including following information: PPARs and neuroinflammation relation sheet, PPARs and CCK4 relation sheet, PPARs and anxiety common target sheet, and amygdaloid expressed gene sheet. The data on relations between objects on pathways. The downloadable online version can be accessed via ResearchGate: https://www.researchgate.net/publication/340418909_Supplemental_Materials_Role_of_PPARs_in_progression_of_anxiety. (Supplementary Materials)
D. J. Kupfer, “Anxiety and DSM-5,” Dialogues in Clinical Neuroscience, vol. 17, no. 3, pp. 245-246, 2015.View at: Google Scholar
F. Thibaut, “Anxiety disorders: a review of current literature,” Dialogues in Clinical Neuroscience, vol. 19, no. 2, pp. 87-88, 2017.View at: Google Scholar
C. Cackovic, S. Nazir, and R. Marwaha, “Panic disorder (attack),” in StatPearls, StatPearls Publishing, Treasure Island, FL, USA, 2020.View at: Google Scholar
B. Bandelow and S. Michaelis, “Epidemiology of anxiety disorders in the 21st century,” Dialogues in Clinical Neuroscience, vol. 17, no. 3, pp. 327–335, 2015.View at: Google Scholar
S. Chakraborty, S. J. Tripathi, T. R. Raju, and B. S. Shankaranarayana Rao, “Brain stimulation rewarding experience attenuates neonatal clomipramine-induced adulthood anxiety by reversal of pathological changes in the amygdala,” Progress in Neuro-Psychopharmacology & Biological Psychiatry, vol. 103, article 110000, 2020.View at: Publisher Site | Google Scholar
V. Zoete, A. Grosdidier, and O. Michielin, “Peroxisome proliferator-activated receptor structures: ligand specificity, molecular switch and interactions with regulators,” Biochimica et Biophysica Acta (BBA) - Molecular and Cell Biology of Lipids, vol. 1771, no. 8, pp. 915–925, 2007.View at: Publisher Site | Google Scholar
A. P. Nesterova, A. Yuryev, E. A. Klimov et al., Disease Pathways: An Atlas of Human Disease Signaling Pathways, Elsevier, Waltham, 1st edition, 2019.
I. Aksu, B. Baykara, S. Ozbal et al., “Maternal treadmill exercise during pregnancy decreases anxiety and increases prefrontal cortex VEGF and BDNF levels of rat pups in early and late periods of life,” Neuroscience Letters, vol. 516, no. 2, pp. 221–225, 2012.View at: Publisher Site | Google Scholar
C. Sharma, A. Pradeep, R. G. Pestell, and B. Rana, “Peroxisome proliferator-activated receptor gamma activation modulates cyclin D1 transcription via beta-catenin-independent and cAMP-response element-binding protein-dependent pathways in mouse hepatocytes,” The Journal of Biological Chemistry, vol. 279, no. 17, pp. 16927–16938, 2004.View at: Publisher Site | Google Scholar
S. Lim, C. J. Jin, M. Kim et al., “PPARγ gene transfer sustains apoptosis, inhibits vascular smooth muscle cell proliferation, and reduces neointima formation after balloon injury in rats,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 26, no. 4, pp. 808–813, 2006.View at: Publisher Site | Google Scholar
A. Nakanishi and I. Tsukamoto, “n-3 polyunsaturated fatty acids stimulate osteoclastogenesis through PPARγ-mediated enhancement of c-Fos expression, and suppress osteoclastogenesis through PPARγ-dependent inhibition of NFkB activation,” The Journal of Nutritional Biochemistry, vol. 26, no. 11, pp. 1317–1327, 2015.View at: Publisher Site | Google Scholar
M. C. Holmes, C. T. Abrahamsen, K. L. French, J. M. Paterson, J. J. Mullins, and J. R. Seckl, “The mother or the fetus? 11beta-hydroxysteroid dehydrogenase type 2 null mice provide evidence for direct fetal programming of behavior by endogenous glucocorticoids,” The Journal of Neuroscience, vol. 26, no. 14, pp. 3840–3844, 2006.View at: Publisher Site | Google Scholar
A. C. Campos, G. N. Vaz, V. M. Saito, and A. L. Teixeira, “Further evidence for the role of interferon-gamma on anxiety- and depressive-like behaviors: involvement of hippocampal neurogenesis and NGF production,” Neuroscience Letters, vol. 578, pp. 100–105, 2014.View at: Publisher Site | Google Scholar
Y. Okabayashi, S. Nagasaka, G. Kanzaki, N. Tsuboi, T. Yokoo, and A. Shimizu, “Group 1 innate lymphoid cells are involved in the progression of experimental anti-glomerular basement membrane glomerulonephritis and are regulated by peroxisome proliferator-activated receptor α,” Kidney International, vol. 96, no. 4, pp. 942–956, 2019.View at: Publisher Site | Google Scholar
O. Y. Alshogran, A. A. Khalil, A. O. Oweis, S. M. Altawalbeh, and M. A. Y. Alqudah, “Association of brain-derived neurotrophic factor and interleukin-6 serum levels with depressive and anxiety symptoms in hemodialysis patients,” General Hospital Psychiatry, vol. 53, pp. 25–31, 2018.View at: Publisher Site | Google Scholar
E. Avolio, G. Fazzari, M. Mele et al., “Unpredictable chronic mild stress paradigm established effects of pro- and anti-inflammatory cytokine on neurodegeneration-linked depressive states in hamsters with brain endothelial damages,” Molecular Neurobiology, vol. 54, no. 8, pp. 6446–6458, 2017.View at: Publisher Site | Google Scholar
B. Leonard and M. Maes, “Mechanistic explanations how cell-mediated immune activation, inflammation and oxidative and nitrosative stress pathways and their sequels and concomitants play a role in the pathophysiology of unipolar depression,” Neuroscience and Biobehavioral Reviews, vol. 36, no. 2, pp. 764–785, 2012.View at: Publisher Site | Google Scholar
E.-Y. N. Wagner, M.-P. F. Strippoli, V. Ajdacic-Gross et al., “Generalized anxiety disorder is prospectively associated with decreased levels of interleukin-6 and adiponectin among individuals from the community,” Journal of Affective Disorders, vol. 270, pp. 114–117, 2020.View at: Publisher Site | Google Scholar
L. V. d’Uscio, T. He, A. V. R. Santhanam, L.-J. Tai, R. M. Evans, and Z. S. Katusic, “Mechanisms of vascular dysfunction in mice with endothelium-specific deletion of the PPAR-δ gene,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 306, no. 7, pp. H1001–H1010, 2014.View at: Publisher Site | Google Scholar
K.-D. Wagner, M. Benchetrit, L. Bianchini, J.-F. Michiels, and N. Wagner, “Peroxisome proliferator-activated receptor β/δ (PPARβ/δ) is highly expressed in liposarcoma and promotes migration and proliferation,” The Journal of Pathology, vol. 224, no. 4, pp. 575–588, 2011.View at: Publisher Site | Google Scholar
J. T. Garretson, B. J. W. Teubner, K. L. Grove, A. Vazdarjanova, V. Ryu, and T. J. Bartness, “Peroxisome proliferator-activated receptor γ controls ingestive behavior, agouti-related protein, and neuropeptide Y mRNA in the arcuate hypothalamus,” The Journal of Neuroscience, vol. 35, no. 11, pp. 4571–4581, 2015.View at: Publisher Site | Google Scholar
N. Bougarne, R. Paumelle, S. Caron et al., “PPARα blocks glucocorticoid receptor α-mediated transactivation but cooperates with the activated glucocorticoid receptor α for transrepression on NF-κB,” Proceedings National Academy of Sciences United States of America, vol. 106, no. 18, pp. 7397–7402, 2009.View at: Publisher Site | Google Scholar
X. Chen, M. Li, W. Sun et al., “Peroxisome proliferator-activated receptor alpha agonist-induced down-regulation of hepatic glucocorticoid receptor expression in SD rats,” Biochemical and Biophysical Research Communications, vol. 368, no. 4, pp. 865–870, 2008.View at: Publisher Site | Google Scholar
T. Varga, Z. Czimmerer, and L. Nagy, “PPARs are a unique set of fatty acid regulated transcription factors controlling both lipid metabolism and inflammation,” Biochimica et Biophysica Acta (BBA) - Molecular Basis of Disease, vol. 1812, no. 8, pp. 1007–1022, 2011.View at: Publisher Site | Google Scholar
M. Gazouli, Z.-X. Yao, N. Boujrad, J. C. Corton, M. Culty, and V. Papadopoulos, “Effect of peroxisome proliferators on Leydig cell peripheral-type benzodiazepine receptor gene expression, hormone-stimulated cholesterol transport, and steroidogenesis: role of the peroxisome proliferator-activator receptor α,” Endocrinology, vol. 143, no. 7, pp. 2571–2583, 2002.View at: Publisher Site | Google Scholar
O. Petyunina, M. Kopytsya, I. Kuznetsov, and I. Vyshnevska, “Prognostication of clinical outcomes after stemi: the role of vascular endothelial growth factor-A,” Georgian Medical news, vol. 279, pp. 79–86, 2018.View at: Google Scholar
T. Tomita, N. Yasui-Furukori, S. Tsuchimine, A. Kaneda, and S. Kaneko, “Relationships between vascular endothelial growth factor levels and temperament and character inventory traits in healthy Japanese subjects,” Neuropsychobiology, vol. 69, no. 1, pp. 1–5, 2014.View at: Publisher Site | Google Scholar
C. Röhrl, U. Kaindl, I. Koneczny et al., “Peroxisome-proliferator-activated receptors γ and β/δ mediate vascular endothelial growth factor production in colorectal tumor cells,” Journal of Cancer Research and Clinical Oncology, vol. 137, no. 1, pp. 29–39, 2011.View at: Publisher Site | Google Scholar
R. Grau, C. Punzón, M. Fresno, and M. A. Iñiguez, “Peroxisome-proliferator-activated receptor alpha agonists inhibit cyclo-oxygenase 2 and vascular endothelial growth factor transcriptional activation in human colorectal carcinoma cells via inhibition of activator protein-1,” The Biochemical Journal, vol. 395, no. 1, pp. 81–88, 2006.View at: Publisher Site | Google Scholar
H. Shen, E. Oesterling, A. Stromberg, M. Toborek, R. Mac Donald, and B. Hennig, “Zinc deficiency induces vascular pro-inflammatory parameters associated with NF-kappaB and PPAR signaling,” Journal of the American College of Nutrition, vol. 27, no. 5, pp. 577–587, 2008.View at: Publisher Site | Google Scholar
J.-J. Zhou, J.-D. Ma, Y.-Q. Mo et al., “Down-regulating peroxisome proliferator-activated receptor-gamma coactivator-1 beta alleviates the proinflammatory effect of rheumatoid arthritis fibroblast-like synoviocytes through inhibiting extracellular signal-regulated kinase, p38 and nuclear factor-kappaB activation,” Arthritis Research & Therapy, vol. 16, no. 5, p. 472, 2014.View at: Publisher Site | Google Scholar
T. Zimmermann, J. C. Bartsch, A. Beer et al., “Impaired anandamide/palmitoylethanolamide signaling in hippocampal glutamatergic neurons alters synaptic plasticity, learning, and emotional responses,” Neuropsychopharmacology, vol. 44, no. 8, pp. 1377–1388, 2019.View at: Publisher Site | Google Scholar
E. Domi, S. Uhrig, L. Soverchia et al., “Genetic deletion of neuronal PPARγ enhances the emotional response to acute stress and exacerbates anxiety: an effect reversed by rescue of amygdala PPARγ function,” The Journal of Neuroscience, vol. 36, no. 50, pp. 12611–12623, 2016.View at: Publisher Site | Google Scholar
E. Domi, F. F. Caputi, P. Romualdi et al., “Activation of PPARγ attenuates the expression of physical and affective nicotine withdrawal symptoms through mechanisms involving amygdala and hippocampus neurotransmission,” The Journal of Neuroscience, vol. 39, no. 49, pp. 9864–9875, 2019.View at: Publisher Site | Google Scholar
O. Okonkwo and A. Adeyinka, “Biochemistry, cholecystokinin (CCK),” in StatPearls, StatPearls Publishing, Treasure Island, FL, USA, 2020.View at: Google Scholar
M. Kellner, “Experimental panic provocation in healthy man-a translational role in anti-panic drug development?” Dialogues in Clinical Neuroscience, vol. 13, no. 4, pp. 485–493, 2011.View at: Google Scholar
J. Andre, B. Zeau, M. Pohl, F. Cesselin, J.-J. Benoliel, and C. Becker, “Involvement of cholecystokininergic systems in anxiety-induced hyperalgesia in male rats: behavioral and biochemical studies,” The Journal of Neuroscience, vol. 25, no. 35, pp. 7896–7904, 2005.View at: Publisher Site | Google Scholar
F. Noble, S. A. Wank, J. N. Crawley et al., “International Union of Pharmacology. XXI. Structure, distribution, and functions of cholecystokinin receptors,” Pharmacological Reviews, vol. 51, no. 4, pp. 745–781, 1999.View at: Google Scholar
P. Dubreuil, P. Fulcrand, M. Rodriguez, H. Fulcrand, J. Laur, and J. Martinez, “Novel activity of angiotensin-converting enzyme. Hydrolysis of cholecystokinin and gastrin analogues with release of the amidated C-terminal dipeptide,” Biochemical Journal, vol. 262, no. 1, pp. 125–130, 1989.View at: Publisher Site | Google Scholar