MedPharmRes
University of Medicine and Pharmacy at Ho Chi Minh City
Original article

Structure-based virtual screening of plant-derived natural compounds as potential PPARα agonists for the treatment of dyslipidemia

Phuong Thuy Viet Nguyena,*https://orcid.org/0000-0002-0233-8692, Truong Le Quang Voahttps://orcid.org/0000-0003-4350-1342, Thao Anh Nguyenahttps://orcid.org/0000-0002-7265-9096, Phuong Ngoc Khanh Hoahttps://orcid.org/0000-0002-1747-9199, Bao Hoang Gia Nguyenahttps://orcid.org/0000-0001-9872-9390
aFaculty of Pharmacy, University of Medicine and Pharmacy at Ho Chi Minh City, Ho Chi Minh city 700000, Viet Nam
*Address correspondence: Phuong Thuy Viet Nguyen at the Faculty of Pharmacy, University of Medicine and Pharmacy at Ho Chi Minh City, Ho Chi Minh city 700000, Viet Nam; Tel: +084-919520708; E-mail: ntvphuong@ump.edu.vn

© Copyright 2021 MedPharmRes. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Mar 31, 2021; Revised: Jul 12, 2021; Accepted: Jul 27, 2021

Published Online: Dec 31, 2021

Abstract:

Background:

Nowadays, metabolic disorders such as dyslipidemia have become serious health problems in the modern world. PPARs are regulators of numberous metabolic pathways, hence there has been a huge increase in the development and use of the PPARs agonists, especially PPARα agonists as main therapeutic of dyslipidemia.

Objectives:

The study aimed to explore potential plant-derived natural compounds as PPARα agonist agent for drug discovery of dyslipidemia.

Methods:

Structure-based virtual screening through molecular docking was conducted for 142 bioactive compounds from 29 medicinal plants on the main binding site of PPARα (PDB ID: 5HYK). Binding affinities and binding interactions between the ligands and PPARα were investigated.

Results:

Screening results showed that 34 compounds had strong binding affinities into the PPARα (binding affinities of less than -8.0 kcal.mol-1), including 20 flavonoid, 4 terpenoid and 10 alkaloid compounds. Flavonoid was found as the best group which fitted well in the binding site of the PPARα. Top compounds were identified, including formononetin from Thermopsis alterniflora (-10.2 kcal.mol-1), diosmetin from Musa spp. (-10.1 kcal.mol-1), luteolin from Elsholtzia ciliate (-9.9 kcal.mol-1); steviol from Stevia rebaudiana (-9.4 kcal.mol-1); and tuberocrooline from Stemona tuberosa (-10.5 kcal.mol-1), respectively. These compounds showed the potential agonistic activities due to forming the hydrogen bonds as well as hydrophobic interactions with four key residues of the receptor such as Ser280, Tyr314, His440 and Tyr464.

Conclusions:

These potential natural compounds may provide useful information in the drug design and discovery for anti-dyslipidemia agents.

Keywords: structure-based virtual screening; molecular docking; PPARα; natural compounds; dyslipidemia

1. INTRODUCTION

Nowadays, with the change of living standards and lifestyle, metabolic disorders such as dyslipidemia (or often hyperlipidemia) and obesity have become serious health problems in the modern world [1]. These diseases can lead to high incidence of morbidity and mortality in both men and women as well as cause an economic burden to the society [1]. Dyslipidemia is characterized by a high levels of total or low-density lipoprotein (LDL) cholesterol, elevated triglycerides and/ or a low levels of high-density lipoprotein (HDL) cholesterol [2]. Dyslipidemia is considered as the biggest contributing factor to the development of athrosclerosis and cardiovascular diseases which is the first cause of death in both developed and developing countries [2]. According to the World Health Organisation, there are about 50% of patients of ischemic heart disease associated with dyslipidemia and more than 4 million deaths every year [3].

Peroxisome Proliferator-Activated Receptors (PPARs) are the nuclear hormone receptors superfamily (class II) [4] including three subtypes, PPARα, PPARβ or δ, and PPARγ, respectively [5]. The PPARs are ligand-activated transcription factors which play a crucial role in the regulation of metabolic processes, energy homeostatis by inducing or repressing target genes [6]. Each subtype displays distinct functions in different cell types [5]. The PPARα was known for controlling and regulating lipid metabolism and inflammation [7]. This receptor is expressed in tissues which have high rate of fatty acid catabolism and in tissues involved in lipid oxidation such as cardiac muscle, liver, kidney, skeletal and adrenal glands [8]. The PPARδ appears in skeletal muscle and adipose tissues which is best known for skin homeostatis [8] and for regulatation of cholesterol, adipogenesis metabolism, and colon cancer [9]. The PPARγ is present in adipose tissue, vascular smooth muscles and immune cells repsonible for energy storage [5] by inducing lipogenesis and fat storage in the tissue as well as improving insulin sensitivity in skeletal muscle [10].

All three PPARs are regulators of numberous metabolic pathways, hence there is a huge increase in the development and use of the PPARs agonists, especially PPARα agonists as main therapeutics of dyslipidemia during the last decade [8]. Number of current marketed drugs for dyslipidemia treatment targets the PPARs. For example, fibrates drugs are activators of PPARα to reduce triacylglycerols used for the treatment of dyslidemia and thiazolidinediones as activatiors of PPARγ for treating hyperglycemia in the type-2 diabetes. Beside, there are other molecular targets for anti-dyslipidemia such as inhibitions of cholesteryl ester transfer protein, cholesterol absorption by binding to cholesterol transporter NPC1L1 (Niemann-pick C1-like1) protein, cholesterol O-acyltransferase enzymes involved in re-esterification of absorbed cholesterol, cholesterol-metabolizing cytochrome P450 or activtion of AMP-activated protein kinase and omega-3 fatty acids, etc [8]. However, there is still rarely drug compounds for those targets [8]. Thus, the PPARs, in particular PPARα has been the main target for anti-dyslipidemia.

The three-dimensional structures of human PPARα in complex with several agonists were solved (PDB codes: 5HYK, 4BCR, 3SP6, 3KDT, 3G8I, 3FEI, 3ET1, 2REW, 2GTK, 1I7G). The structure of human PPARα is very similar to both of the PPARγ and PPARδ [5]. The structure includes four functional domains (A/B, C, D, E/F) (5), in which there are the ligand-independent activation function 1 (AF1) in the A/B domain; the conserved central DNA binding domain in the C domain; docking site for cofactors in the D domain and the ligand binding domain in the E region [5]. The PPARs have multiple binding sites depending on various binding ligands [11]. There are two known binding sites, the main binding pocket for full agonistic and the second one for partial agonistic activities on the PPARα [12].

The ligand binding domain of PPARα is a 12-helix which forms a large hydrophobic pocket [5] (T-shaped cavity) with the volume of about 1300 Å3 [13]. The central pocket spans the region located in between the C-terminal helix 12 (forming one side of a second activation function, AF2 helix) and the 3-stranded antiparallel β sheet [13]. After ligation with the agonist such as fibrate drugs, PPARα undergoes conformation changes [6, 14]. The conserved hydrogen bonds involving the AF2 helix were formed between the protein and the agonists which was found as the most structural differences between the apo protein and the agonist-bound protein [13]. The full agonists stabilise an active conformation of the AF2 helix and enhances heterodimerization process with the retinoid X receptors, promoting recruitment of nuclear receptor coactivators and gene transcription [12, 13]. In addition, a loss of this stabilizating interactions suggested the partial agonist [14]. Without ligands, the PPARα binds to promoters of the target genes leading to transcription repress [6].

However, fibrate drugs as PPARα agonists demonstrated to cause some side effects such as hepatomegaly or liver dysfunction, gastrointestinal disturbance, increase of creatinine levels, etc. [15]. As the results, there is necessary to explore novel compounds for the anti-hyperlipidemic activity with less toxicities. Some medicinal plants with bioactive compounds have been reported for anti-lipid effects as rich source for new effective and safe medicines [16-20] despite of required sufficient evidence for their activities. For example, anthraquinones from Rheum officinale; ginsenoside, ginseng, and polysaccharides from Radix ginseng; triterpenes from Rhizoma alismatis,... reduced triglycerides, LDL-C [19]. Furthermore, some compounds were reported as PPARα agonists such as picrasidine C (an alkaloid from the root of Picrasma quassioides); bixin (a carotenoid from the pericarp of the seeds of Bixa orellana); naringenin (flavanones from the dried, immature fruit of Citrus aurantium); secoiridoid excelside B and some metabolites from Fraxinus excelsior L.etc [20]. Therefore, this study aimed to discover the plant-derived natural compounds against PPARα as PPARα agonists to assist in drug design for anti-dyslipidemia. By using molecular docking for screening, 142 investigated compounds from 29 medicinal plants were docked into the PPARα. The results of ligand binding affinities and their interactions with this receptor were combined to select the potential natural compounds as potential lipid lowering agents.

2. MATERIALS AND METHOD

Molecular docking was applied in in silico screening to select natural compound with high binding affinity into the PPARα. Initially, the diverse bioactive compounds belonging to four main groups, flavonoid, terpenoid, alkaloid and saponin, respectively with their lower-lipid effects were selected for this study. To assist for drug discovery in the next step, the medicinal plants were then searched for containing the compounds. In total, there were 142 natural compounds from 29 medicinal plants including alkaloid (51 compounds), flavonoid (36 compounds), terpenoid (33 compounds) and saponin (22 compounds) were chosen [20-27]. Molecular docking was conducted to investigate binding interactions and to reveal potential natural compounds for treatment of lipid disorders using AutoDock Vina software version 1.1.2 (an open-source program) [28].

Preparation of protein: The 3D crystal structure of the complex PPARα was retrieved from the Protein Data Bank (PDB ID: 5HYK – resolution: 1.83 Å, https://www.rcsb.org) [6]. The structure of PPARα has a co-crystallized ligand-2-methyl-2-[4-(naphthalen-1-yl)phenoxy]propanoic acid (or AL26-29) as the full agonist [6]. The residues making up the main binding pocket are Cys276, Thr279, Ser280, Tyr314, Ile317, Met330, His440 and Tyr464 [6] which covered the co-crystallised-ligand. Using AutoDock Tools 1.5.7rc1, this receptor was prepared for docking; with removing all water molecules and heteroatoms and adding polar hydrogen atoms and Kollman charges into the protein structure. Redocking was carried out by extracting the co-crystallized ligand from the experimental structure and docking the ligand into the binding pocket of PPARα encompassing the native ligand.

Preparation of ligand: Ligands (142 compounds) were prepared in 2D by Chem3D Ultra program (CDX file). Open Babel was employed to convert all 2D structures to 3D structures whose energies were subsequently minimized with the YASARA Energy Minimization server (https://www.yasara.org).

Docking parameters: The docking parameters included the coordinate parameters of the center x, y, z of 7.57 Å; 32.34 Å; 23.882 Å. The grid box was centered on the native ligand with the dimension of 24 x 24 x 24 (Å)3; spacing distance = 1 Å; with default exhaustiveness = 8.

Evaluation of docking results: The docking results were evaluated by the ligand binding affinities (kcal.mol-1), binding pose and the possible interactions between the key residues for biological activities in the target such as Tyr314 for maintaining protein active conformation of protein and for selectivity ligands. In addition, the fibrate drugs as the PPARα agonists including fenofibrate, gemfibrozil, clofibrate, bezafibrate, ciprofibrate were used for docking into the active site of the receptor as reference drugs. The lower binding affinity, the better ligand. BIOVIA Discovery Studio Visualizer 2020 (a free version, downloaded from the website https://discover.3ds.com/) was used for visualisation, assisting in analysis and creating images of 3D models of the protein and the binding mode of protein-ligand.

3. RESULTS

To evaluate the docking protocol, redocking of the co-crystallized ligand into the main binding site of PPARα was conducted. The results showed that this native ligand (the naphthalenic derivatives) bound well into the PPARα with good binding affinity (-11.3 kcal.mol-1) and the root mean square deviation (RMSD) between the docked structure and the native one, using only movable heavy atoms (i.e., only ligand atoms, not hydrogen), was 1.46 Å (less than 2.00 Å). The binding mode and interactions of this ligand with the PPARα were mostly similar to the experimental structure. For example, the polar head of the ligand could form the hydrogen bonds with the side chains of key residues of the active site of the PPARα such as Tyr464 on the H12 helix responsible for agonistic activities towards PPARα and Tyr314, and with the OH of Ser280 [6]. The ligand also interacted with the PPARα through hydrophobic interactions with Phe273 and Phe351. These results demonstrated the reliability of the docking program, so the docking protocol could be then used for screening process.

Molecular docking

All of 142 natural compounds of the medicinal plants were located into the main binding pocket of PPARα (PDB id: 5HYK) in docking which showed their potential activities against PPARα. These compounds were classified into the subgroups: alkaloid (51 compounds), terpenoid (33 compounds), flavonoid (36 compounds) and saponin (22 compounds) for binding analysis. The results were also compared with those results obtained with the agonisit ligands such as fibrate drugs as reference compounds of PPARα. It showed that flavonoid was the best group of binding well into the PPARα compared to the other groups. Notably, total of 34/142 natural compounds had high affinity for binding into the PPARα (binding affinities of less than -8.0 kcal.mol-1), including 20 flavonoid compounds, 4 terpenoid compounds and 10 alkaloid compounds (Table 1). Of which, the top compounds for each group included: flavonoid group: formononetin from Thermopsis alterniflora (-10.2 kcal.mol-1), diosmetin from Musa spp. (-10.1 kcal.mol-1), luteolin from Elsholtzia ciliate (-9.9 kcal.mol-1); one terpenoid compound: steviol from Stevia rebaudiana (-9.4 kcal.mol-1); and one alkaloid compound: tuberocrooline from Stemona tuberosa (-10.5 kcal.mol-1) (Figure 1 and 2). These compounds strongly accommodated and formed good interactions with the residues of the PPARα binding site as the reference drugs did. Structure and binding affinity relationships of the ligands were also taken for analysis. However, there were four compounds, namely saponarin, vincenin-2, stemonone and tuberostemoenone did not show negative binding affinities as the ligands did not accommodate fully in the binding pocket.

Table 1. Screening results of 142 natural compounds based on molecular docking from 29 medicinal plants into the PPARα (PDB id: 5HYK) with the binding affinities (kcal.mol-1).
No Medicinal plant Natural compound Binding affinity (kcal.mol-1)
1 Momordica charantia Cucurbitaceae Momordicine 28 -1.8
Gypsogenin -3.7
Goyaglicoside A -6.9
Goyaglicoside C -7.2
Goyaglicoside E -8.2
Goyaglicoside F -7.2
Vicine -6.9
Goyaglicoside H -6.0
2 Stevia rebaudiana Asteraceae Steviol -9.4
3 Ganoderma lucidum Ganodermataceae Ganoderic acid A -1.2
Ganoderic acid B -1.1
Ganoderic acid D -1.5
Ganoderic acid E -3.8
Ganoderic acid F 4.0
Ganoderic acid G -0.3
Ganoderic acid H -5.3
Ganoderic acid I -3.9
Ganoderic acid J -2.7
Ganoderic acid K 5.9
Ganoderic acid L -0.8
Ganoderic acid M -0.2
Ganoderic acid N -5.0
Acid ganodermic S -0.5
Acid ganodermic P2 -0.2
4 Cynara scolymus Asteraceae Cynaropicrin -8.9
Dehydrocynaropicrin -8.3
Grossheimin -7.0
5 Sesamum indicum Pedaliaceae 3-epibartogenic acid -2.3
6 Celastrus hindsii Celastraceae Celasdin A -2.1
Celasdin B -0.9
Celasdin C -6.0
Maytenfolone A -0.6
7 Cynachium stauntonii Apocynaceae Hancolupenone -1.4
Hancolupenol -2.1
8 Folium nelumbinis Nuciferae Nuciferin -5.3
Pronuciferine -4.0
Roemerine -7.4
2-hydroxy-1-methoxyaporphine -6.9
Dehydronuciferine -2.6
N-nornuciferine -3.5
Astragalin -3.8
9 Erythrina orientalis Fabaceae Soyasapogenol B -3.2
10 Codonopsis pilosula Campanulaceae Radicamine A -8.2
Codonopyrrolidiums A -5.9
Codonopyrrolidiums B -7.1
Codonopsinols A -7.9
Codonopsinols B -7.6
Codonopsinols C -7.4
Codonopiloside A -7.5
Tryptophan -7.7
Perlolyrine -9.1
Nicotinic acid -5.5
11 Stemona tuberosa Stemonaceae Isostenine (neostenine) -5.4
Tuberostemonine H -2.3
Tuberostemonine N -2.6
Tuberostemonine K -6.1
Neotuberostemonol -4.2
Epi-Bisdehydroneotuberostemonine J (aka epibisdehydrotuberostemonine J) -8.5
9a-Bisdehydrotuberostemonine -1.0
9a-Bisdehydrotuberostemonine A -7.0
Tridehydrotuberostemonine -2.0
Bisdehydroneostemoninine -9.1
Bisdehydrostemoninine A -8.5
Bisdehydrostemoninine B -8.5
Bisdehydrostemoninine -8.1
Isobisdehydrostemoninine -7.3
Oxystemoninine -3.7
Stemoenonine -5.0
9a-O-Methylstemoenonine -4.3
Oxystemoenonine -4.0
Tuberostemospironine -7.5
10-Hydroxycroomine -7.5
6-Hydroxycroomine (aka 6ahydroxycroomine) -7.6
Dehydrocroomine -5.8
Tuberospironine -7.6
Neotuberostemoninol -6.7
Sessilifoliamide F -2.5
Tuberostemoline -0.7
Tuberocrooline -10.5
1,9a-seco-Stemoenonine -6.2
Tuberostemoenone 0.9
Croomine -8.9
Stemoninoamide -9.8
Stemotinine -6.1
Tuberostemonone 2.9
Neotuberostemonine (aka tuberostemonine LG) -1.0
12 Platycodon grandiflorum Campanulaceae Platycodin A -4.5
Platycodin C -4.8
Platycodin D -3.9
Deapioplatycodin D -4.1
13 Kochia scoparia Chenopodiaceae Momoridin Ic -0.9
14 Aesculus turbinate Aesculaceae Escin Ia -2.5
Escin IIa -2.2
Escin Ib -3.4
Escin IIb -3.9
15 Fructus Momordicae grosvenorii Cucurbitaceae Mogroside IV -1.2
Mogroside V -0.7
16 Schefflera heptaphylla Araliaceae Silphioside F -2.1
Copteroside B -3.7
Gypsogenin 3-O-D-glucuronide -0.5
17 Schefflera octophylla Araliaceae Sessiloside -6.2
Chiisanoside -7.1
18 Gynostemma pentaphyllum Cucurbitaceae Rutin -3.4
Damulin A -7.4
Damulin B -6.7
19 Lysimachia foenum-graecum Primulaceae Foenumoside B -1.3
20 Dioscorea nipponica Makino Dioscoreaceae Trillin -2.9
21 Panax ginseng Araliaceae Ginsenoside Rb1 -5.5
22 Scutellaria baicalensis Lamiaceae Formononetin -10.2
Daidzein -9.7
Chrysin -9.6
Isoscitellarein -9.4
Daidzin -9.1
Apigenin-7-glucoside -8.9
Luteolin-7-rutinoside -7.2
Puerarin -5.9
23 Elsholtzia ciliate Lamiaceae Luteolin -9.9
Apigenin -9.4
Kumatakenin -8.1
Linarin -6.5
24 Musa spp. Musaceae Epicatechin -8.9
Gallocatechin -8.9
Dimer procyanidin -7.7
Semilicoisoflavone -7.7
25 Glycyrrhiza uralensis Fabaceae Chalcone -8.9
Isoliquiritigenin -8.6
Ononin -8.6
Licoflavonol -8.4
Isoangustone -8.2
Licochalcone -7.7
Glyasperin C -7.6
7-O-methylluteone -7.1
Licorisoflavan A -6.9
Dehydroglyasperin -6.3
26 Cynanchum Apocynaceae Quercetin -9.2
27 Thermopsoside alternilora Fabaceae Thermopsoside -8.8
28 Cucumis sativus Cucurbitaceae Saponarin 4-O-glucoside -8.7
Isovitexin -7.2
Isoorientin -6.8
Saponarin 1.5
Vicenin-2 2.5
29 Senna alata Fabaceae Diosmetin -10.1
Download Excel Table
mpr-5-4-34-g1
Figure 1. Binding modes of 3 top flavonoid compounds, namely formononetin from Thermopsis alterniflora (-10.2 kcal.mol-1), diosmetin from Musa spp. (-10.1 kcal.mol-1), luteolin from Elsholtzia ciliate (-9.9 kcal.mol-1); respectively into the PPARα (PDB id: 5HYK) with green line represented for hydrogen bond, pink line for hydrophobic contacts
Download Original Figure
mpr-5-4-34-g2
Figure 2. Binding modes of 2 top compounds, one terpenoid compound: steviol from Stevia rebaudiana (-9.4 kcal.mol-1); and one alkaloid compound: tuberocrooline from Stemona tuberosa (-10.5 kcal.mol-1), respectively into the PPARα (PDB id: 5HYK) with green line represented for hydrogen bond, pink line for hydrophobic contacts
Download Original Figure
Binding modes and interactions analysis

Flavonoid compounds such as formononetin (-10.2 kcal.mol-1) and diosmetin (-10.1 kcal.mol-1) had good binding affinities to the binding pocket of PPARα which showed the potential agonistic activities. This was explained by the presence of -OH groups in the structure of these flavonoids created the hydrogen bonds with four key residues of the receptor such as Ser280, Tyr314, His440 and Tyr464 [6] as well as the aromatic rings made hydrophobic interactions with the target residues which led to the tightly attachment of the compounds into the binding site of PPARα. These results were compatible with the experimental results that PPAR agonists interact to form similar interactions with polar residues, especially with Tyr464 localized in AF2 which was important for stabilizing the active conformation of protein [6]. In addition, the three other flavonoids such as luteolin (-9.9 kcal.mol-1), daidzein (-9.7 kcal.mol-1), and chrysin (-9.6 kcal.mol-1) were also good ones towards the binding affinities on the PPARα. The ligands with the simple backbone could attach the hydrophobic pocket better than the isoflavon and flavon compounds with the presence of glucosides or alkyl groups, for example: luteolin (-9.9 kcal.mol-1) > luteolin-7-rutinoside (-7.2 kcal.mol-1); daidzein (-9.7 kcal.mol-1) > daidzin (-9.1 kcal.mol-1); and apigenin (-9.4 kcal.mol-1) > apigenin-7-glucoside (-8.9 kcal.mol-1). The results could be related to the experimental resuls that the aglycone penetrates easily due to the high lipophilicity and low molecular weight [29]. The other isoflavon and flavon compounds with the presence of glucoside or alkyl groups such as daidzin, thermopsoside, ononin, isoangustone, apigenin-7-glucoside, saponarin 4-O-glucoside, licoflavonol, kumatakenin did not show strong binding affinities (Table 1). Structures of the eight compounds with the good binding affinities of less than - 8.0 kcal.mol-1 and their binding modes into the PPARα (PDB id: 5HYK) were presented in Figure 3 and 4.

mpr-5-4-34-g3
Figure 3. Structures of 8 flavonoid compounds belonging to isoflavon and flavon with the good binding affinities of less than -8.0 kcal.mol-1 into the PPARα (PDB id: 5HYK)
Download Original Figure
mpr-5-4-34-g4
Figure 4. Binding modes of 8 flavonoid compounds belonging to isoflavon and falvon with the good binding affinities of less than -8.0 kcal.mol-1 into the PPARα (PDB id: 5HYK)
Download Original Figure

Alkaloid compounds bound to the PPARα with binding affinities ranging from -1.0 to -10.5 kcal.mol-1. Ten alkaloid compounds including tuberocrooline (-10.5 kcal.mol-1 had good affinities on the PPARα due to forming the hydrophobic interactions between the aromactic rings containing nitrogen and some important amino acids of binding site such as His440 (for example, perlolyrine: -9.1 kcal.mol-1). The compounds also generated the hydrogen bonds through their side chains having some functional groups (OH or COOH groups) with residues Ser280, Tyr314 and His440 of PPARα (for example, radicamine A: -8.2 kcal.mol-1). These interactions were shared by all PPARα agonist like fibrates and the other ligands reported in the PDB [6, 12], with data showed in Table 2. The alkaloid possess the simple aromactic ring like radicamine A and perlolysine could attach further in the binding pocket than the ligands having bulky aromatic groups. However, replacing pyrrole ring into pyrrolidine ring made the compound not going deeper into the cavity or extending the length of side chain of pyrrole which led to reduce the binding affinities, such as the case of compound 9a-bisdehydrotuberostemonie with the pyrrole ring (-7.0 kcal.mol-1) > neotuberostemonol with the pyrrolidine ring (-4.2 kcal.mol-1), and stemoninoamide (-9.8 kcal.mol-1) > bisdehydrostemonine A and B (-8.5 kcal.mol-1). When pyrrolidine opens, for example tuberocrooline (-10.5 kcal.mol-1), the compound made the C=O group of tetrahydrofuran interacted with Tyr314 and also formed more hydrogen bond with the residue Ala455 of target, so this compound showed strongly affinity than, croomine (-8.9 kcal.mol-1). All of Figure 5, 6 and 7 illustrated 10 top alkaloid compounds, namely radicamine A, perlolyrine, epibisdehydroneotuberostemonine J, bisdehydroneostemoninine, stemoninoamide, bisdehydrostemoninine A, bisdehydrostemoninine B, bisdehydrostemoninine, tuberocrooline, croomine, respectively with their good binding affinities of less than -8.0 kcal.mol-1 and their binding modes into the PPARα.

Table 2. Binding interactions of 15 compounds including 10 top alkaloid compounds, namely radicamine A, perlolyrine, epibisdehydroneotuberostemonine J, bisdehydroneostemoninine, stemoninoamide, bisdehydrostemoninine A, bisdehydrostemoninine B, bisdehydrostemoninine, tuberocrooline, croomine, and 5 reference compounds (fenofibrate, gemfibrozil, clofibrate, bezafibrate and ciprofibrate), respectively and the PPARα (PDB id: 5HYK)
No Compound name Binding affinity (kcal.mol-1) Hydrogen bonds Hydrophobic interactions
1 Radicamine A -8.2 Ser280, His440, Tyr464 Ile354, Ile447, Ala454
2 Perlolyrine -9.1 Phe273, Cys276, Ile354, His440, Val444, Ile447, Ala454
3 Epi-Bisdehydroneotuberostemonine J -8.5 Ser280, Tyr314 Ile354, Val444, Ile447, Lys448, Leu456
4 Bisdehydroneostemoninine -9.1 Ala455
5 Bisdehydrostemoninine A -8.5 Ser280, His440 Ile354, Val444, Ile447, Leu456
6 Bisdehydrostemoninine B -8.5 Cys276 Ile354, Val444, Ile447, Leu456
7 Bisdehydrostemoninine -8.1 His440, Leu456 Phe273, Ile354, Val444, Ile447, Leu456
8 Tuberocrooline -10.5 Ser280, Tyr314, Ala455
9 Croomine -8.9 Ser280, His440 Ile354, Val444, Ile447
10 Stemoninoamide -9.8 Ile354, His440, Val444, Ile447, Leu456
11 Fenofibrate -10.2 Ser280, His440 Phe273, Ile354, Val444, Ile447, Ala454, Leu456 Val270, Phe351,
12 Gemfibrozil -8.2 His440 Val444, Ile447, Lys448, Ala454, Ala455, Leu456
13 Clofibrate -7.3 Phe273, Cys276, Ile354, Met355
14 Bezafibrate -9.1 Ser280, Tyr314, His440, Val444 Ile354, Val444, Ala454 Val270, Phe273,
15 Ciprofibrate -8.6 Ser280, Tyr314, His440 Phe351, Val444, Ile447, Ala454
Download Excel Table
mpr-5-4-34-g5
Figure 5. The molecular structures of 10 top alkaloid compounds, namely radicamine A, perlolyrine, epibisdehydroneotuberostemonine J, bisdehydroneostemoninine, stemoninoamide, bisdehydrostemoninine A, bisdehydrostemoninine B, bisdehydrostemoninine, tuberocrooline, croomine, respectively with their binding affinities into the PPARα (PDB id: 5HYK)
Download Original Figure
mpr-5-4-34-g6
Figure 6. Binding modes of 10 top alkaloid compounds, namely radicamine A, perlolyrine, epibisdehydroneotuberostemonine J, bisdehydroneostemoninine, stemoninoamide, bisdehydrostemoninine A, bisdehydrostemoninine B, bisdehydrostemoninine, tuberocrooline, croomine, respectively with the good binding affinities of less than -8 kcal.mol-1 into the PPARα (PDB id: 5HYK)
Download Original Figure
mpr-5-4-34-g7
Figure 7. Binding interactions of 10 top alkaloid compounds, namely radicamine A, perlolyrine, epibisdehydroneotuberostemonine J, bisdehydroneostemoninine, stemoninoamide, bisdehydrostemoninine A, bisdehydrostemoninine B, bisdehydrostemoninine, tuberocrooline, croomine, respectively and the residues at the binding site of the PPARα (PDB id: 5HYK) with green line represented for hydrogen bond, pink line for hydrophobic contacts
Download Original Figure

For terpenoids and alkaloids, the compounds belonging to these groups could get into the active site of the PPARα, however, their binding affinities were not as good as the other groups. Of which, three terpenoid compounds, namely steviol, cynarpicrine and dehydrocynaropicrine and goyaglycoside E interacted well with the residues Tyr314, Tyr464, Ile447 at the active site of the PPARα by hydrogen bonding with their affinities of -9.4; -8.9, -8.3 and -8.2 kcal.mol-1, respectively (presented in Figure 8). It could be seen that the presence of OH group in the structures such as steviol, cyanopicrin, dehydrocyanopicrin increased the binding affinities. The results agreed with the experimental data that -OH group was essential for the activities in terpenoids (20). Only one saponin compound, diosgenin displayed good binding affinity of -7.5 kcal.mol-1. The saponin ginsenosides, such as ginsenoside Rb1 did not show high binding affinity (-5.5 kcal.mol-1). The experimental results also reported that the compound responsible for the bioactivities of ginseng which demonstrated the inhibition of PPARα in vitro and in vivo with a consequent increased serum concentrations of total cholesterol, triglycerides, and HDL cholesterol [30].

mpr-5-4-34-g8
Figure 8. Binding modes of 4 top terpenoid compounds, namely steviol, cynaropicrin, dehydrocynaropicrin, goyaglicoside E, respectively with the good binding affinities of less than -8.0 kcal.mol-1 into the PPARα (PDB id: 5HYK)
Download Original Figure

Moreover, the docking results of fibrate drugs as the PPARα agonists (fenofibrate, gemfibrozil, clofibrate, bezafibrate, ciprofibrate) into the PPARα as reference drugs showed that these drugs bound closely to the binding site of the receptor with high binding affinities: fenofibrate (-10.2 kcal.mol-1), gemfibrozil (-8.2 kcal.mol-1), clofibrate (-7.3 kcal.mol-1), bezafibrate (-9.1 kcal.mol-1), ciprofibrate (-8.6 kcal.mol-1) (Table 2). Comparison of docking results of the 142 natural compounds and the reference drugs, it showed that most of the investigated compounds bound to the main binding pocket and strongly mimicking the observed interactions like the reference ones did. The reference drugs interacted well with the receptor by hydrogen bonds and hydrophobic interactions with some important residues such as Tyr464 important for recruiting co-activators (6) and Phe273, Tyr314, Phe351 and Tyr464. In particular, the standard polar interactions with Ser280, Tyr314 and Tyr464 presented in the all agonists complexed to the receptor (6) which were found in all of the top investigated compounds.

Overall, top compounds selected from each group including formononetin, diosmetin, luteolin, daidzein, chrysin (flavonoid compounds); steviol (terpenoid compound) and tuberrocrooline (alkaloid compound) possess high binding affinities. They occupied the same position and shared ability to form interactions with the residues of binding pocket of the PPARα. However, taking into account of forming hydrogen bonds network with the key residues Tyr464, His440, Tyr314 and Ser280 for agonistic activities, the following compounds, namely formononetin from Thermopsis alterniflora (-10.2 kcal.mol-1), diosmetin from Musa spp. (-10.1 kcal.mol-1), luteolin from Elsholtzia ciliate (-9.9 kcal.mol-1); steviol from Stevia rebaudiana (-9.4 kcal.mol-1); and tuberocrooline from Stemona tuberosa (-10.5 kcal.mol-1) were selected. They could be used as the potential PPARα agonists for further biological testings for their anti-dyslipidemia effects.

4. DISCUSSION

The natural compounds from medicinal plants have been of interest in drug discovery for the treatment of many diseases, including dyslipidemia in recent years. Therefore, in this study, structure-based virtual screening through molecular docking was conducted for 142 investigated natural compounds into the PPARα (PDB id: 5HYK) to discover potential compounds for anti-dyslipidemia. Through the criteria including ligand binding affinities, binding modes and binding interactions, flavonoid compounds performed the most potential compounds. These compounds formed favorable bindings with the PPARα by hydrogen bonds with key residues of the target, Ser280, Tyr314 and Tyr464, crucial for agonistic activities. Of which, two flavonoids, formononetin from Thermopsis alterniflora (-10.2 kcal.mol-1), diosmetin from Musa spp. (-10.1 kcal.mol-1) and luteolin from Elsholtzia ciliate (-9.9 kcal.mol-1) were the top compounds of high binding affinities which suggesting they might be potential compounds as PPARα agonists for further testing toward anti-dyslipidemia effects.

This results were consistent with the reported results from the experiments that formononetin was the most potent activators of PPARα with the value of EC50 of less than 1.0 mol/L [20, 31]. Diosmetin demonstrated inhibition of fat accumulation in liver and epididymal tissues, however diosmetin had little effect on the parameters related to dyslipidemia [32]. Diosmetin (3’,5,7-trihydroxy-4’-methoxy-flavone), a flavone aglycone have been displayed to have several pharmacological effects such as anti-inflammatory and hepatoprotective role in an in vivo model [33] or antihyperglycemic activities [34]. Another study reported that diosmetin treatment inhibited fat accumulation in liver and epididymal tissues and improved glucose intolerance by lowering glucose levels during a glucose tolerance test [32]. Furthermore, luteolin has been shown to have anti-inflammatory activity as well as to protect against tumor necrosis factor-alpha and especially to prevent lipid accumulation and decrease total cholesterol, triglyceride and LDL-cholesterol levels [35]. Therefore, the findings in this study suggested the mechanism of these compounds toward to PPARα.

Moreover, two other compounds, one terpenoid compound: steviol from Stevia rebaudiana (-9.4 kcal.mol-1); and one alkaloid compound: tuberocrooline from Stemona tuberosa (-10.5 kcal.mol-1) were also identified as the hit natural compounds for PPARα agonists. Steviol, the terpene aglycone, the final metabolite of all steviol glycosides was reported with some bioactivites such as antidiabetic effects on streptozotocin-induced diabetic mice [36], antihyperglycemic, antihypertensive, antioxidant and anti-inflammatory effects, etc [37]. Steviol causes a decrease in glucose accumulation in intestinal ring tissue, liver and kidney and also enhances insulin secretion [36]. Tuberocrooline has not been demonstrated for anti-dyslipidemia, but this compound was known for the antitussive activity [38]. Furthermore, the formononetin has been proved for anticancer activity [39], diosmetin for the biological activity on alloxan diabetic rats [40], steviol for antioxidant capacity [41] and tuberocrooline for the acetylcholinesterase inhibitory activity [42].

Therefore, the good compounds were selected, namely formononetin from Thermopsis alterniflora, diosmetin from Musa spp., luteolin from Elsholtzia ciliate; steviol from Stevia rebaudiana; and tuberocrooline from Stemona tuberosa, respectively for further steps in drug discovery and design of PPARα agonists for anti-hyperlipidemia activities. Combination of the experimental results of these compounds for dyslipidemia effects and the results obtained in this study, the current research also suggested the mechanism of these compounds towards to the PPARα for agonistic activities. However, there were still some limits in this study. It is required more number of investigated compounds and partial agonistic on the PPARα or dual agonist/ pan-agonistic activities combining on PPARγ could be investigated in the screening process of better anti-dyslipidemia and also for PPARα agonists compounds.

Conclusion

Metabolic disorders such as dyslipidemia have been in the list of serious health problems in the modern world. This disease is considered as the major contributing factor to the mortality of men and women in both developed and developing countries. With the aims of identifying the potential plant-derived natural compounds as PPARα agonists for assissiting in drug discovery of anti-hyperlipidemia, 142 compounds were docked into the structure of PPARα (PDB ID: 5HYK). The screening results discovered 34 compounds including 20 flavonoid, 4 terpenoid and 10 alkaloid compounds strongly accommodated and formed good interactions with the residues of the PPARγ binding site as the reference drugs did. Flavonoid was the best compounds attaching into the binding site of the PPARα. Top compounds for anti-dyslipidemia were identified, namely formononetin from Thermopsis alterniflora, diosmetin from Musa spp., luteolin from Elsholtzia ciliate; steviol from Stevia rebaudiana; and tuberocrooline from Stemona tuberosa. These compounds formed favorable interactions with the binding site of PPARα through hydrogen bonds and hydrophobic interactions with key residues of the target, Ser280, Tyr314 and Tyr464, leading to agonistic activities. Combination of these results and the experimental results of these compounds for dyslipidemia effects, the current study also suggested the mechanism of the compounds as the PPARα agonists. Therefore, this study provided useful information for further drug discovery and design of potential lipid-lowering agents as the replacement for the current PPARα agonists like fibrate drugs.

LIST OF ABBREVIATIONS

LDL:

low-density lipoprotein

HDL:

high-density lipoprotein

PPARs:

Peroxisome Proliferator-Activated Receptors

NPC1L1:

Niemann-pick C1-like1

AF1:

activation function 1

AF2:

activation function 2

RMSD:

root mean square deviation.

FUNDING

The authors would like to thank University of Medicine and Pharmacy at Ho Chi Minh City for the financial support.

CONFLICT OF INTEREST

The authors declare that there is no conflict of interest.

ACKNOWLEDGEMENTS

The authors would like to thank Department of Pharmaceutical Information Technology, Faculty of Pharmacy, University of Medicine and Pharmacy at Ho Chi Minh City for their support and funding of conducting this research.

ETHICAL CONSIDERATIONS

Not available.

AUTHORS’ CONTRIBUTIONS

The authors have participated in completing the work including TLQV for docking of terpenoids; TAN for docking of saponins; PNKH for docking of alkaloids; BHGN for docking of fibrates and falvonoid compounds; and PTVN for conducting the project, supervising the process and drafting, editing and proofreading manuscript.

REFERENCES

1.

Pitsavos C, Panagiotakos D, Weinem M, Stefanadis C. Diet, exercise and the metabolic syndrome. Rev Diabet Stud. 2006 Fall; 3(3):118-126

2.

Nelson RH. Hyperlipidemia as a risk factor for cardiovascular disease. Prim Care. 2013 Mar; 40(1):195-211

3.

Moini J, Ahangari R, Miller C, Samsam M. Global Health Complications of Obesity. 1st ed USA: Elsevier Academic Press. 2020; Chapter 18, Perspective on economics and obesityp p. 411-423

4.

Bain DL, Heneghan AF, Connaghan-Jones KD, Miura MT. Nuclear receptor structure: implications for function. Annual review of physiology. 2007 Mar; 69:201-20

5.

Christofides A, Konstantinidou E, Jani C, Boussiotis VA. The role of peroxisome proliferator-activated receptors (PPAR) in immune responses. Metabolism. 2021 Aug; 2021114(154338):1-13

6.

Capelli D, Cerchia C, Montanari R, Loiodice F, Tortorella P, Laghezza A, et al. Structural basis for PPAR partial or full activation revealed by a novel ligand binding mode. Sci Rep. 2016 Oct; 20166(34792):1-28

7.

Harmon GS, Lam MT, Glass CK. PPARs and lipid ligands in inflammation and metabolism. Chem rev. 2011 Oct; 12111(10):6321-40

8.

Mahamuni SP, Khose RD, Menaa F, Badole SL. Therapeutic approaches to drug targets in hyperlipidemia. BioMedicine. 2012 Dec; 2(4):137-146

9.

Coll T, Rodrïguez-Calvo R, Barroso E, Serrano L, Eyre E, Palomer X, et al. Peroxisome proliferator-activated receptor (PPAR) beta/delta: a new potential therapeutic target for the treatment of metabolic syndrome. Curr Mol Pharmacol. 2009 Jan; 2(1):46-55

10.

Zieleniak A, Wójcik M, Woźniak LA. Structure and physiological functions of the human peroxisome proliferator-activated receptor gamma. Arch Immunol Ther Exp. 2008 Sep-Oct; 56(5):331-45

11.

Kuwabara N, Oyama T, Tomioka D, Ohashi M, Yanagisawa J, Shimizu T, et al. Peroxisome proliferator-activated receptors (PPARs) have multiple binding points that accommodate ligands in various conformations: phenylpropanoic acid-type PPAR ligands bind to PPAR in different conformations, depending on the subtype. J Med Chem. 2012 Jan; 2655(2):893-902

12.

Bernardes A, Souza PC, Muniz JR, Ricci CG, Ayers SD, Parekh NM, et al. Molecular mechanism of peroxisome proliferator-activated receptor α activation by WY14643: a new mode of ligand recognition and receptor stabilization. J Mol Biol. 2013 Aug; 23425(16):2878-93

13.

Cronet P, Petersen JFW, Folmer R, Blomberg N, Sjöblom K, Karlsson U, et al. Structure of the PPARα and -γ Ligand Binding Domain in Complex with AZ 242; Ligand Selectivity and Agonist Activation in the PPAR Family. Structure. 2001 Aug; 9(8):699-706

14.

Nettles KW, Greene GL. Ligand control of coregulator recruitment to nuclear receptors. Annu Rev Physiol. 2005 Mar; 67:309-33

15.

Zodda D, Giammona R, Schifilliti S. Treatment Strategy for Dyslipidemia in Cardiovascular Disease Prevention: Focus on Old and New Drugs. Pharmacy (Basel). 2018 Jan; 216(1) 10:1-35

16.

Rouhi-Boroujeni H, Rouhi-Boroujeni H, Heidarian E, Mohammadizadeh F, Rafieian-Kopaei M. Herbs with anti-lipid effects and their interactions with statins as a chemical anti-hyperlipidemia group drugs: A systematic review. ARYA Atheroscler. 2015; 11(4):244-251.

17.

Dhaliya SA, Surya AS, Tomy DV, Carla B, Kumar AR, Christudas S. A review of hyperlipidemia and medicinal plants. IntJAPSBMS. 2013 Oct-Dec; 2(4):219-237.

18.

Ji X, Shi S, Liu B, Shan M, Tang D, Zhang W, et al. Bioactive compounds from herbal medicines to manage dyslipidemia. Biomed Pharmacother. 2019 Oct; 118(109338):1-12

19.

Guo M, Liu Y, Gao Z-Y, Shi D-z. Chinese Herbal Medicine on Dyslipidemia: Progress and Perspective. Evid-Based Complementary Altern Med. 2014 Feb; 163036:1-12

20.

Rigano D, Sirignano C, Taglialatela-Scafati O. The potential of natural products for targeting PPARα. Acta Pharm Sin B. 2017 Jul; 7(4):427-438

21.

Ma C, Wang J, Chu H, Zhang X, Wang Z, Wang H, et al. Purification and characterization of aporphine alkaloids from leaves of Nelumbo nucifera Gaertn and their effects on glucose consumption in 3T3-L1 adipocytes. Int J Mol Sci. 2014; 15(3):3481-3494

22.

Xiao J, Tian B, Xie B, Yang E, Shi J, Sun Z. Supercritical fluid extraction and identification of isoquinoline alkaloids from leaves of Nelumbo nucifera Gaertn. Euro Food Res Technol. 2010 May; 231(3):407-414

23.

Mongalo NI, Mashele SS, Makhafola TJ. Ziziphus mucronata Willd. (Rhamnaceae): it’s botany, toxicity, phytochemistry and pharmacological activities. Heliyon. 2020 Apr; 6(4)e03708:1-58

24.

He J-y, Ma N, Zhu S, Komatsu K, Li Z-Y, Fu W-M. The genus Codonopsis (Campanulaceae): A review of phytochemistry, bioactivity and quality control. J Nat Med. 2015 Jan; 69(1):1-21

25.

Pilli R, Rosso G, Oliveira MdC. The Chemistry of Stemona Alkaloids: An Update. Nat Prod Rep. 2010 Nov; 27:1908-37

26.

Hussain MS, Jahan N, Rashid MM, Hossain MS, Chen U, Rahman N. Antihyperlipidemic screening and plasma uric acid reducing potential of Momordica charantia seeds on Swiss albino mice model. Heliyon. 2019 May; 5(5)e01739:1-18

27.

Hu X-Y, Qin F-Y, Lu X-F, Zhang L-S, Cheng Y-X. Three New Polyynes from Codonopsis pilosula and Their Activities on Lipid Metabolism. Molecules. 2018 Apr; 23(4) 887:1-15

28.

Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS, et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem. 2009 Dec; 30(16):2785-91

29.

Murota K, Nakamura Y, Uehara M. Flavonoid metabolism: the interaction of metabolites and gut microbiota. Biosci Biotechnol Biochem. 2018 Apr; 82(4):600-610

30.

Yoon M, Lee H, Jeong S, Kim JJ, Nicol CJ, Nam KW, et al. Peroxisome proliferator-activated receptor alpha is involved in the regulation of lipid metabolism by ginseng. Br J Pharmacol. 2003 Apr; 138(7):1295-302

31.

Shen P, Liu MH, Ng TY, Chan YH, Yong EL. Differential effects of isoflavones, from Astragalus membranaceus and Pueraria thomsonii, on the activation of PPARalpha, PPARgamma, and adipocyte differentiation in vitro. The J Nutr. 2006 Apr; 136(4):899-905

32.

Chung S, Kim HJ. Comparative study of the effects of diosmin and diosmetin on fat accumulation, dyslipidemia, and glucose intolerance in mice fed a high-fat high-sucrose diet. Foof Sci Nutr. 2020 Nov; 8(11):5976-5984

33.

Yang Y, Gong X-B, Huang L-G, Wang Z-X, Wan R-Z, Zhang P, et al. Diosmetin exerts anti-oxidative, anti-inflammatory and anti-apoptotic effects to protect against endotoxin-induced acute hepatic failure in mice. Oncotarget. 2017 May; 8(19):30723-30733

34.

Shakthi Deve A, Sathish Kumar T, Kumaresan K, Rapheal VS. Extraction process optimization of polyphenols from Indian Citrus sinensis - as novel antiglycative agents in the management of diabetes mellitus. J Diabetes Metab Disord. 2014 Jan; 713(1) 11:1-21

35.

Li J, Dong JZ, Ren YL, Zhu JJ, Cao JN, Zhang J, et al. Luteolin decreases atherosclerosis in LDL receptor-deficient mice via a mechanism including decreasing AMPK-SIRT1 signaling in macrophages. Exp Ther Med. 2018 Sep; 16(3):2593-2599

36.

Myint KZ, Chen J-M, Zhou Z-Y, Xia Y-M, Lin J, Zhang J. Structural dependence of antidiabetic effect of steviol glycosides and their metabolites on streptozotocin-induced diabetic mice. J Sci Food Agric. 2020 Aug; 100(10):3841-3849

37.

Chatsudthipong V, Muanprasat C. Stevioside and related compounds: therapeutic benefits beyond sweetness. Pharmacol Ther. 2009 Jan; 121(1):41-54

38.

Lin L-G, Pak-Ho Leung H, Zhu J-Y, Tang C-P, Ke C-Q, Rudd JA, et al. Croomine- and tuberostemonine-type alkaloids from roots of Stemona tuberosa and their antitussive activity. Tetrahedron. 2008 Oct; 64(44):10155-10161

39.

Tay K-C, Tan LT-H, Chan CK, Hong SL, Chan K-G, Yap WH, et al. Formononetin: A Review of Its Anticancer Potentials and Mechanisms. Front Pharmacol. 2019 Jul; 10 820:1-52

40.

Michael HN, Salib JY, Eskander EF. Bioactivity of diosmetin glycosides isolated from the epicarp of date fruits, Phoenix dactylifera, on the biochemical profile of alloxan diabetic male rats. Phytother Res. 2013 May; 27(5):699-704

41.

Karimi M, Hashemi J, Ahmadi A, Abbasi A, Esfahani M. Study on the bioactivity of steviol and isosteviol in stevia (Stevia rebaudiana Bertoni). Acta Physiol Plant. 2014 Dec; 36(12):3243-3248

42.

Ramli RA, Lie W, Pyne SG. Alkaloids from the Roots of Stichoneuron caudatum and Their Acetylcholinesterase Inhibitory Activities. J Nat Prod. 2014 Apr; 77(4):894-901