Evidence-Based Complementary and Alternative Medicine

Volume 2014 (2014), Article ID 192452, 19 pages

http://dx.doi.org/10.1155/2014/192452

Research Article

In Silico Identification of Potent PPAR- Agonists from Traditional Chinese Medicine: A Bioactivity Prediction, Virtual Screening, and Molecular Dynamics Study

Kuan-Chung Chen¹ and Calvin Yu-Chian Chen2,3

¹School of Pharmacy, China Medical University, Taichung 40402, Taiwan

²School of Medicine, College of Medicine, China Medical University, Taichung 40402, Taiwan

³Department of Biomedical Informatics, Asia University, Taichung 41354, Taiwan

Received 17 December 2013; Accepted 25 January 2014; Published 26 May 2014

Academic Editor: Fuu-Jen Tsai

Copyright © 2014 Kuan-Chung Chen and Calvin Yu-Chian Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The peroxisome proliferator-activated receptors (PPARs) related to regulation of lipid metabolism, inflammation, cell proliferation, differentiation, and glucose homeostasis by controlling the related ligand-dependent transcription of networks of genes. They are used to be served as therapeutic targets against metabolic disorder, such as obesity, dyslipidemia, and diabetes; especially, PPAR- is the most extensively investigated isoform for the treatment of dyslipidemic type 2 diabetes. In this study, we filter compounds of traditional Chinese medicine (TCM) using bioactivities predicted by three distinct prediction models before the virtual screening. For the top candidates, the molecular dynamics (MD) simulations were also utilized to investigate the stability of interactions between ligand and PPAR- protein. The top two TCM candidates, 5-hydroxy-L-tryptophan and abrine, have an indole ring and carboxyl group to form the H-bonds with the key residues of PPAR- protein, such as residues Ser289 and Lys367. The secondary amine group of abrine also stabilized an H-bond with residue Ser289. From the figures of root mean square fluctuations (RMSFs), the key residues were stabilized in protein complexes with 5-Hydroxy-L-tryptophan and abrine as control. Hence, we propose 5-hydroxy-L-tryptophan and abrine as potential lead compounds for further study in drug development process with the PPAR- protein.

1. Introduction

The peroxisome proliferator-activated receptors (PPARs) belonged to the nuclear receptor superfamily of ligand-inducible transcription factors. They are “fatty acid sensors” related to regulation of lipid metabolism, inflammation, cell proliferation, differentiation, and glucose homeostasis by controlling the related ligand-dependent transcription of networks of genes [1–3]. There are three different isoforms of PPARs in mammal, which are PPAR-, PPAR-, and PPAR-/. They have different tissue distributions and responses to different ligands [4–6]. PPARs are used to be served as therapeutic targets against metabolic disorder, such as obesity, dyslipidemia, and diabetes; especially, PPAR- is the most extensively investigated isoform for the treatment of dyslipidemic type 2 diabetes [7–10]. It is a well-known receptor located in fat for antidiabetic insulin sensitizers and has the functions related to adipogenesis, lipogenesis, and glucose homeostasis [11–13]. In rat stroke models, PPAR- has been served as a brain protector against ischemic cerebral infraction [14].

Nowadays, increasing numbers of drug are designed with a target protein against a specific disease [15–18], as increasing numbers of distinct mechanism of diseases have been identified by the researches [19–26]. Recently, the compounds from traditional Chinese herb have been proven to have the therapeutic effects [27–30]. In previous researches, many compounds of traditional Chinese medicine (TCM) have been indicated as potential candidates of lead compounds against cancer [31–34], neuropathic pain [35], stroke [36, 37], and virus infection [38, 39].

In the former study, we aim to detect potential candidates from TCM compounds as agonists targeting PPAR-, PPAR-, and PPAR- [40]. However, a compound which had a higher binding affinity with target protein may not always obtain a higher bioactivity. In this paper, we aimed to focus on the target protein of PPAR- and filter TCM compounds using bioactivities predicted by three distinct prediction models before the virtual screening. The molecular dynamics (MD) simulations were also utilized to investigate the stability of interactions between ligand and PPAR- protein in the docking pose under dynamic conditions. We attempt to identify the potent TCM compounds with higher bioactivities and binding affinity for PPAR- protein and discuss the functional group of these candidates and common binding residues of PPAR- protein in their docking pose.

2. Materials and Methods

2.1. Data Collection

After TCM compounds from TCM database, Taiwan [41], were filtered by Lipinski et al.’s rule of five [42], a total of 9,029 nonduplicate compounds were prepared by Prepare Ligand module in Discovery Studio 2.5 (DS2.5) to adjust the ionization state to physiological setting for virtual screening. For calculating the pharmacokinetics properties, ADMET Descriptors model in DS2.5 was employed to calculate the aqueous solubility, CYP2D6 binding, hepatotoxicity, and plasma protein binding (PPB) as absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties for each compound.

The X-ray crystallography structure of the human peroxisome proliferator-activated receptor gamma (PPAR-) protein was obtained from RCSB Protein Data Bank with PDB ID: 3K8S [43]. After protein preparation, the chain A of PPAR- protein was used as target protein for virtual screening, and T2384, cocrystallized in PPAR- protein, was used as control.

2.2. Biological Activity Prediction Using Multiple Linear Regression (MLR), Support Vector Machine (SVM), and Bayes Network Toolbox (BNT) Models

For the prediction of biological activity for the TCM compounds, three distinct prediction models were constructed with the pEC₅₀ (log(1/EC₅₀)) value of 20 compounds from Rikimaru et al.’s study [2] as training set. The genetic function approximation module [44] of DS 2.5 was utilized to determine the suitable molecular descriptors for constructing the prediction models, and the fitness of individual model was estimated by square correlation coefficient (). Cross-validation test was used to validate the prediction model. For three distinct prediction models, multiple linear regression and Bayes network toolbox were performed using MATLAB, and support vector machine was performed using LibSVM developed by Chang and Lin [45].

2.3. Docking Simulation

For virtual screening, LigandFit protocol [46] in DS 2.5 was employed to dock each compound into an active site using a shape filter and Monte Carlo ligand conformation generation, and each docked pose was minimized with Chemistry at HARvard Macromolecular Mechanics (CHARMM) force field [47] and evaluated with a set of scoring functions. In addition, LigPlot v.2.2.25 program [48] was employed to identify the interactions between protein and ligand in each docking pose.

2.4. Molecular Dynamics Simulation

Before the molecular dynamics simulation by Gromacs [49], each protein-ligand complex in docking pose has been reprepared. Each ligand was reprepared by SwissParam program [50], and the protein was reprepared with charmm27 force field by Gromacs. The protein-ligand complex was solvated using a water model of TIP3P with a minimum distance of 1.2 Å from the complex and then minimized by steepest descent algorithm [51] with maximum of 5,000 steps. Then a single 10 ps constant temperature (NVT ensemble) equilibration was performed using Berendsen weak thermal coupling method followed by a 40 ns production simulation. For each MD simulation, it adopts the particle mesh Ewald (PME) option with a time step of 2 fs. A series of protocols in Gromacs were employed to analyze the MD trajectories.

3. Results and Discussion

3.1. Biological Activity Predictions

The genetic approximation algorithm determined the six optimum molecular descriptors for constructing prediction models with 20 compounds of training set. The selected descriptors were ES_Sum_sssCH, ES_Count_aaN, BIC, IAC_Mean, CHI_3_P, and JY. These six optimum molecular descriptors can be broadly divided into two groups, which are electronic and special topological descriptors. For electronic topological descriptors, it includes ES_Sum_sssCH, ES_Count_aaN for calculating the sums of the electrotopological state (E-state) values and the counts of each atom type, respectively. For special topological descriptors, BIC and IAC_Mean are bonding information content and mean information of atomic composition, which both belong to Graph-Theoretical InfoContent descriptors [52]. CHI_3_P is a Kier and Hall molecular connectivity index [53]. JY is a Balaban index [54]. According to these selected descriptors, the functional formula of multiple linear regression (MLR) model was constructed as follows:

The support vector machine (SVM) and Bayes network toolbox (BNT) models were also constructed with the identical training set and descriptors. The correlation of predicted and observed activities shown in Figure 1 illustrates the correlation trend and 95% prediction bands for each prediction model. The square correlation coefficients () of training set for MLR, SVM, and BNT models are 0.8442, 0.8536, and 0.7612, respectively. These prediction models are acceptable for predicting activity of PPAR- protein.

(a)

(b)

(c)

Figure 1: Comparative plots of observed versus predicted activity for (a) MLR, (b) SVM, and (c) BNT models. Correlation trend (blue line) and 95% prediction boundaries (enclosed by magenta lines) were shown.

3.2. Docking Simulation

The potent compounds, which have acceptable predicted activities in all three prediction models, have been virtual screening with the target protein. After filtering by the absorption properties, the top TCM candidates ranked by Dock score were listed in Table 1 with their predicted activities and pharmacokinetics properties. Human intestinal absorption model displayed in Figure 2 suggested that the top five TCM candidates may have good absorption.

Table 1: Docking results, predicted pEC₅₀, and ADMET properties for top TCM compounds and T2384.

Name
Dock score
Predicted pEC₅₀
CYP2D6a probability
Hepatotoxicity

probability
PPB

levelb

MLR
SVM
BNT

5-Hydroxy-L-tryptophan
148.721
5.89
6.62
6.59
0.069
0.291
0

Abrine
142.592
6.31
6.63
6.44
0.049
0.642
0

Saussureamine C
135.304
5.90
7.27
7.91
0.415
0.450
0

Saussureamine B
124.688
8.85
7.63
8.00
0.356
0.708
0

Saussureamine A
103.030
7.59
7.81
7.64
0.336
0.754
0

* T2384
77.618
7.52
7.06
8.50
0.069
0.953
2

Control.

aInhibition probability of cytochrome P450 2D6 enzyme.

bPlasma protein binding: 0: binding is <90%; 1: binding is >90%; 3: binding is >95%.

Figure 2: Human intestinal absorption model for top TCM compounds and T2384.

For the docking simulation, the binding site of PPAR- protein was defined by the volume and position of control, T2384 (Figure 3(a)). We visually inspected docking poses of top ranked TCM candidates (Figure 3(b)), 5-hydroxy-L-tryptophan, abrine, and saussureamine C interacting with similar PPAR- binding site residues as control (Figures 3(c)-3(d)). Figure 4 displays the structure of T2384 and top three candidates. According to the docking poses shown in Figure 5, T2384 has interactions with residues Phe264 and Phe363, hydrogen bonds (H-bonds) with residues Cys285 and Lys367, and hydrophobic contact with other nine residues.

(a)

(b)

(c)

(d)

Figure 3: Binding sites and common binding residues for PPAR- protein. PPAR- protein with (a) binding site defined by T2384, (b) docking poses of top TCM compounds and T2384 in the binding site. (c) Top view, (d) side view of docking poses with common binding residues.

(a)

(b)

(c)

(d)

Figure 4: Chemical scaffold of control and top three candidates: (a) T2384, (b) 5-hydroxy-L-tryptophan, (c) abrine, and (d) saussureamine C.

Figure 5: Docking pose of PPAR- complexes with T2384, 5-hydroxy-L-tryptophan, abrine, and saussureamine C, respectively. Hydrophobic contacts between PPAR- protein and each compound determined by LIGPLOT program.

Compared with T2384 in PPAR- protein, the top three TCM candidates have been docked with similar docking poses. Due to the molecular size of three TCM compounds, none of them have interaction with Phe264 as T2384. Except saussureamine C, both of 5-hydroxy-L-tryptophan and abrine have interaction with residue Phe363 as control. However, saussureamine C still has hydrophobic contact with residue Phe363. All top three candidates have similar H-bond with residue Lys367 and hydrophobic contacts with some common residues, such as Leu330 and Met364. Except that 5-hydroxy-L-tryptophan has hydrophobic contact instead of H-bond with residue Cys285, the other two candidates have the similar H-bond with residues Cys285 as T2384. In addition, abrine and saussureamine C also have H-bond with Ser289 and Met364, respectively.

3.3. Molecular Dynamics Simulation

The docking poses in the docking simulation illustrate that the top three TCM candidates have similar interactions with the target proteins as T2384. However, the structure of PPAR- protein is fixed during the progress of docking simulation. As this reason, the molecular dynamics (MD) simulations for each protein-ligand complex were performed to investigate the stability of interactions between ligand and target protein in the docking pose under dynamic conditions and investigate the possible variations for each protein-ligand complex after docking.

The root mean square deviations (RMSDs) and radii of gyration for each protein and ligand in the complexes were illustrated in Figure 6. For RMSD, it calculates the deviation of the structure compared with the starting structure over 40 ns of MD simulation. They indicate that all protein-ligand complexes tend to be stable after 30 ns of MD simulation. Radius of gyration, which measures the mass of the atom relative to the center of mass of the complex, is indicative of the compactness of each complex. As shown in Figure 6, there is no significant variation for the compactness of each complex. Figure 7 illustrates the variation of total energy for each protein-ligand complex over the course of 40 ns MD simulation with the average fluctuations in a cycle of 21 frames shown in the center of each graph. Total energy trajectories indicate that these systems were stabilized for PPAR- protein in the complex with T2384 and top three TCM candidates over the course of 40 ns MD simulation. Figure 8 displays the variation of secondary structure of PPAR- protein and secondary structural feature ratio over the course of 40 ns MD simulation for each complex with T2384 and top three TCM candidates. It indicates that docking with three TCM candidates may not cause the significant differences from docking with the control in the secondary structure of PPAR- protein.

(a)

(b)

(c)

(d)

Figure 6: RMSDs and radii of gyration for PPAR- protein and ligands over 40 ns MD simulation.

(a)

(b)

(c)

(d)

Figure 7: Total energy of PPAR- complexes with (a) T2384, (b) 5-hydroxy-L-tryptophan, (c) abrine, and (d) saussureamine C over 40 ns MD simulation. The average fluctuations in a cycle of 21 frames were illustrated by yellow (T2384), violet (5-hydroxy-L-tryptophan), wine (abrine), and cyan (saussureamine C) line, respectively.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 8: Secondary structure assignments and secondary structural feature ratio variations of PPAR- complexes over 40 ns MD simulation. Residues 1–65 in -axis correspond to residues 207–271, residues 66-250 in -axis correspond to residues 276–460, and residues 251–262 in -axis correspond to residues 465–476.

The representative structures of each complex after MD simulation were identified by the cluster analysis with a RMSD cutoff of 0.1 nm. In Figure 9, it illustrates the RMSD values and graphical depiction of the clusters over 30–40 ns MD simulation. The representative structures of each complex were identified by middle RMSD structure in the major cluster over 30–40 ns MD simulation, which are 38.88 ns (T2384), 39.86 ns (5-hydroxy-L-tryptophan), 39.80 ns (abrine), and 39.96 ns (saussureamine C), respectively. The snapshots and ligand interaction diagrams for each docking pose of the representative structures are illustrated in Figure 10. For T2384, it maintains the H-bonds with residues Cys285 and Lys367 in a nonstatic condition, which may retain the docking pose of T2384 in the binding pocket of PPAR- protein. In addition, the ligand interaction diagram also indicates that T2384 has interactions with common residues in docking simulation. For 5-hydroxy-L-tryptophan, it keeps the H-bond with residue Lys367 in a nonstatic condition and also has an H-bond with residue Ser289 as the docking pose of abrine in the docking simulation. Similarly, abrine has H-bonds with residues Ser289 and Lys367 as well as has an H-bond and interaction with residue Tyr327. The docking pose of saussureamine C in the docking simulation is not stable in a nonstatic condition. In the representative structures after MD simulation, it has H-bonds with residues His449 and Leu476, as well as a interaction with residue Phe282.

(a)

(b)

(c)

(d)

Figure 9: RMSD values (upper left half) and graphical depiction of the clusters (lower right half) of PPAR- complexes during 30–40 ns MD simulation.

Figure 10: Docking poses of middle RMSD structure in the major cluster during 30–40 ns of MD simulation. Snapshots and ligand interaction diagrams for PPAR- protein complexes with T2384 (38.88 ns), 5-hydroxy-L-tryptophan (39.86 ns), abrine (39.80 ns), and saussureamine C (39.96 ns). For 2D diagrams, residues with magenta cycles are involved in hydrogen-bond, charge, or polar interactions, and residues with green cycles are involved in van der Waals interactions.

The H-bonds occupancies for key residues of PPAR- protein in each complex are shown in Table 2 with cutoff of 0.3 nm. Figure 11 displays the variation of these distances over the course of 40 ns MD simulation. For T2384, the potential H-bonds with key residues of PPAR- protein are formed by its sulfonamide group. 5-Hydroxy-L-tryptophan and abrine form H-bonds with residue Lys367 by the carboxyl group. They form H-bonds with residue Ser289 by the indole group in the beginning of MD simulation, but the H-bond for abrine has shifted from the indole group to secondary amine group after 5 ns of MD simulation. In addition, the carboxyl group of abrine also forms a stable H-bond with residue Tyr327 after 7 ns of MD simulation. For saussureamine C, the docking pose in the docking simulation had changed after MD simulation. The H-bonds formed by the carboxyl group are shifted from residue Lys367 to residue Ser289 after 5 ns of MD simulation. In addition, it forms stable H-bonds with residue His449 by its sulfonamide group and heterocycle group after MD simulation.

Table 2: H-bond occupancy for key residues of PPAR- protein with top three candidates and T2384 overall 40 ns molecular dynamics simulation.

Name
H-bond interaction
Occupancy

T2384
Tyr327:HH/O20
8.70%

Tyr327:OH/H40
5.55%

Lys367:HZ3/O20
39.30%

Lys367:HZ3/O21
0.75%

5-Hydroxy-L-tryptophan
Cys285:O/H21
1.30%

Ser289:HG1/N7
7.05%

Ser289:OG/H26
0.15%

Ser289:OG/H21
14.54%

Tyr327:HH/O13
5.25%

Tyr327:HH/O14
4.55%

Lys367:HZ3/O13
49.93%

Lys367:HZ3/O14
0.90%

His449:HE2/O14
7.60%

Abrine
Cys285:O/H30
1.00%

Ser289:HG1/N7
1.45%

Ser289:OG/H22
3.00%

Ser289:OG/H30
90.05%

Tyr327:HH/O14
83.15%

Lys367:HZ3/O13
41.15%

Lys367:HZ3/O14
37.75%

Saussureamine C
Ser289:HG1/O23
32.40%

Ser289:HG1/O24
30.05%

Tyr327:HH/O14
4.95%

Tyr327:HH/O23
6.20%

Tyr327:HH/O24
14.65%

Tyr327:HH/O28
6.60%

Tyr327:OH/H52
6.65%

Lys367:HZ3/O23
2.65%

Lys367:HZ3/O24
2.70%

His449:HE2/O14
14.65%

His449:HE2/O24
3.50%

His449:HE2/N27
3.90%

His449:HE2/O28
77.60%

Leu476:OT1/H50
13.85%

Leu476:OT2/H50
7.00%

H-bond occupancy cutoff: 0.3 nm.

(a)

(b)

Figure 11: Distances of potential H-bonds between PPAR- protein and each compound during 40 ns MD simulation.

The root mean square fluctuations (RMSFs) shown in Figure 12 illustrate the stability of each residue over 30–40 ns MD simulation. Residues Cys285, Lys367, and His449 are stabilized by all top three TCM candidates and T2384. As abrine forms stable H-bond with residues Ser289 and Tyr327, the RMSFs of Ser289 and Tyr327 are much lower in the complex with abrine than with others. For saussureamine C, as the H-bonds with residue Tyr327 are shifted between the heterocycle group, secondary amine group, and sulfonamide group, it causes the highest value of RMSF for residue Tyr327 in the complex with saussureamine C.

Figure 12: RMSFs for residues 207–476 of PPAR- complexes with each compound over 30–40 ns MD simulation. Common binding residues were illustrated with dash lines.

To consider the variation of each ligand during MD simulation, variation of torsion angles during 40 ns of MD simulation for each ligand in the PPAR- complexes is shown in Figure 13. As T2384 is the cocrystallized compound in the PPAR- protein, the docking pose is stable during 40 ns of MD simulation. For 5-hydroxy-L-tryptophan, the docking pose which is also stable during 40 ns of MD simulation except for the hydroxyl group in the indole ring has a 180-degree shift after MD simulation. For abrine, the variation of torsions 10 and 11 at the initial period of MD simulation may be the reason that the H-bond has shifted from the indole group to secondary amine group, and carboxyl group forms a stable H-bond with residue Tyr327 after MD simulation. Torsions 14 and 16 for saussureamine C indicate that the docking pose of saussureamine C has a fluctuation during 15–30 ns of MD; it can also be seen in the ligand RMSD (Figure 6) and the distance variation with residue Tyr327 (Figure 11). The variation of torsion 19 shows that the sulfonamide group of saussureamine C is flexible over MD simulation.

Figure 13: Variation of ligand torsion angles for each of PPAR- complexes during 40 ns of MD simulation. Red and gray lines represent the ligand torsion angle at docking simulation and first conformation of MD simulation, respectively.

4. Conclusion

This study aims to investigate the potent TCM candidates for PPAR- protein. The biologically activities of candidates were predicted by three distinct prediction models (MLR, SVM and BNT) based on their ligand characteristics. After docking simulation, the docking poses of top TCM compounds ranked by the scoring function were validated by the MD simulation. For the top three TCM candidates, both of 5-hydroxy-L-tryptophan and abrine have an indole ring and carboxyl group to form the H-bonds with the key residues of PPAR- protein. The secondary amine group of abrine also stabilized an H-bond with residue Ser289. The key residues were stabilized in protein complexes with 5-Hydroxy-L-tryptophan and abrine as control. For saussureamine C, the interactions of docking pose in the docking simulation are not stable after MD simulation. Hence, we propose 5-hydroxy-L-tryptophan and abrine as potential lead compounds for further study in drug development process with the PPAR- protein.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The research was supported by Grants from the National Science Council of Taiwan (NSC102-2325-B039-001, NSC102-2221-E-468-027-), Asia University (Asia101-CMU-2, 102-Asia-07), and China Medical University Hospital (DMR-102-105, DMR-103-058, DMR-103-001, and DMR-103-096). This study is also supported in part by Taiwan Department of Health Clinical Trial and Research Center of Excellence (DOH102-TD-B-111-004), Taiwan Department of Health Cancer Research Center of Excellence (MOHW103-TD-B-111-03), and CMU under the Aim for Top University Plan of the Ministry of Education, Taiwan.

References

R. M. Evans, G. D. Barish, and Y.-X. Wang, “PPARs and the complex journey to obesity,” Nature Medicine, vol. 10, no. 4, pp. 355–361, 2004.

K. Rikimaru, T. Wakabayashi, H. Abe et al., “A new class of non-thiazolidinedione, non-carboxylic-acid-based highly selective peroxisome proliferator-activated receptor (PPAR) γ agonists: design and synthesis of benzylpyrazole acylsulfonamides,” Bioorganic and Medicinal Chemistry, vol. 20, no. 2, pp. 714–733, 2012.

J. N. Feige, L. Gelman, L. Michalik, B. Desvergne, and W. Wahli, “From molecular action to physiological outputs: peroxisome proliferator-activated receptors are nuclear receptors at the crossroads of key cellular functions,” Progress in Lipid Research, vol. 45, no. 2, pp. 120–159, 2006.

L. Poulsen, M. Siersbaek, and S. Mandrup, “PPARs: fatty acid sensors controlling metabolism,” Seminars in Cell & Developmental Biology, vol. 23, no. 6, pp. 631–639, 2012.

M. Ahmadian, J. M. Suh, N. Hah, et al., “PPARgamma signaling and metabolism: the good, the bad and the future,” Nature Medicine, vol. 19, no. 5, pp. 557–566, 2013.

A. Bugge and S. Mandrup, “Molecular mechanisms and genome-wide aspects of PPAR subtype specific transactivation,” PPAR Research, vol. 2010, Article ID 169506, 12 pages, 2010.

A. R. Miller, “Today's challenges and tomorrow's opportunities: ligands to peroxisome proliferator-activated receptors as therapies for type 2 diabetes and the metabolic syndrome,” Drug Development Research, vol. 67, no. 7, pp. 574–578, 2006.

J. M. Lehmann, L. B. Moore, T. A. Smith-Oliver, W. O. Wilkison, T. M. Willson, and S. A. Kliewer, “An antidiabetic thiazolidinedione is a high affinity ligand for peroxisome proliferator-activated receptor γ (PPARγ),” The Journal of Biological Chemistry, vol. 270, no. 22, pp. 12953–12956, 1995.

N. Kubota, Y. Terauchi, H. Miki et al., “PPARγ mediates high-fat diet-induced adipocyte hypertrophy and insulin resistance,” Molecular Cell, vol. 4, no. 4, pp. 597–609, 1999.

B. M. Forman, J. Chen, and R. M. Evans, “The peroxisome proliferator-activated receptors: ligands and activators,” Annals of the New York Academy of Sciences, vol. 804, pp. 266–275, 1996.

W. W. Cheatham, “Peroxisome proliferator-activated receptor translational research and clinical experience,” The American Journal of Clinical Nutrition, vol. 91, no. 1, pp. 262S–266S, 2010.

A. M. Sharma and B. Staels, “Review: peroxisome proliferator-activated receptor γ and adipose tissue—understanding obesity-related changes in regulation of lipid and glucose metabolism,” The Journal of Clinical Endocrinology and Metabolism, vol. 92, no. 2, pp. 386–395, 2007.

Y.-T. Li, L. Li, J. Chen et al., “7-chloroarctinone-b as a new selective PPARγ antagonist potently blocks adipocyte differentiation,” Acta Pharmacologica Sinica, vol. 30, no. 9, pp. 1351–1358, 2009.

J.-S. Wu, W.-M. Cheung, Y.-S. Tsai et al., “Ligand-activated peroxisome proliferator-activated receptor-γ protects against ischemic cerebral infarction and neuronal apoptosis by 14-3-3ε upregulation,” Circulation, vol. 119, no. 8, pp. 1124–1134, 2009.

C.-L. Jao, S.-L. Huang, and K.-C. Hsu, “Angiotensin I-converting enzyme inhibitory peptides: inhibition mode, bioavailability, and antihypertensive effects,” BioMedicine, vol. 2, no. 4, pp. 130–136, 2012.

M. Lin, S. Tsai, F. Wang, F.-H. Liu, J.-N. Syu, and F.-Y. Tang, “Leptin induces cell invasion and the upregulation of matrilysin in human colon cancer cells,” BioMedicine, vol. 3, no. 4, pp. 174–180, 2013.

K.-P. Su, “Inflammation in psychopathology of depression: clinical, biological, and therapeutic implications,” BioMedicine, vol. 2, no. 2, pp. 68–74, 2012.

M. A. Leissring, E. Malito, S. Hedouin et al., “Designed inhibitors of insulin-degrading enzyme regulate the catabolism and activity of insulin,” PLoS ONE, vol. 5, no. 5, Article ID e10504, 2010.

Y. Jiang, X. Li, W. Yang, et al., “PKM2 regulates chromosome segregation and mitosis progression of tumor cells,” Molecular Cell, vol. 53, no. 1, pp. 75–87, 2014.

I. C. Chou, W.-D. Lin, C.-H. Wang, et al., “Association analysis between Tourette's syndrome and two dopamine genes (DAT1, DBH) in Taiwanese children,” BioMedicine, vol. 3, no. 2, pp. 88–91, 2013.

T. Yamamoto, W.-C. Hung, T. Takano, and A. Nishiyama, “Genetic nature and virulence of community-associated methicillin-resistant Staphylococcus aureus,” BioMedicine, vol. 3, no. 1, pp. 2–18, 2013.

C.-H. Wang, W.-D. Lin, D.-T. Bau, I.-C. Chou, C.-H. Tsai, and F.-J. Tsai, “Appearance of acanthosis nigricans may precede obesity: an involvement of the insulin/IGF receptor signaling pathway,” BioMedicine, vol. 3, no. 2, pp. 82–87, 2013.

W. Chen, Z. Wang, C. Jiang, and Y. Ding, “PP2A-mediated anticancer therapy,” Gastroenterology Research and Practice, vol. 2013, Article ID 675429, 10 pages, 2013.

Y.-M. Chang, B. K. Velmurugan, W.-W. Kuo, et al., “Inhibitory effect of alpinate Oxyphyllae fructus extracts on Ang II-induced cardiac pathological remodeling-related pathways in H9c2 cardiomyoblast cells,” BioMedicine, vol. 3, no. 4, pp. 148–152, 2013.

Y. M. Leung, K. L. Wong, S. W. Chen, et al., “Down-regulation of voltage-gated Ca²⁺ channels in Ca²⁺ store-depleted rat insulinoma RINm5F cells,” BioMedicine, vol. 3, no. 3, pp. 130–139, 2013.

S. P. Mahamuni, R. D. Khose, F. Menaa, and S. L. Badole, “Therapeutic approaches to drug targets in hyperlipidemia,” BioMedicine, vol. 2, no. 4, pp. 137–146, 2012.

K. C. Chen, S. S. Chang, F. J. Tsai, and C. Y. Chen, “Han ethnicity-specific type 2 diabetic treatment from traditional Chinese medicine?” Journal of Biomolecular Structure & Dynamics, vol. 31, no. 11, pp. 1219–1235, 2013.

K.-C. Chen, M.-F. Sun, S.-C. Yang et al., “Investigation into potent inflammation inhibitors from traditional Chinese medicine,” Chemical Biology & Drug Design, vol. 78, no. 4, pp. 679–688, 2011.

K. C. Chen, Y. R. Jian, M. F. Sun, T. T. Chang, C. C. Lee, and C. Y. Chen, “Investigation of silent information regulator 1 (Sirt1) agonists from Traditional Chinese Medicine,” Journal of Biomolecular Structure & Dynamics, vol. 31, no. 11, pp. 1207–1218, 2013.

S.-C. Hsu, J.-H. Lin, S.-W. Weng, et al., “Crude extract of Rheum palmatum inhibits migration and invasion of U-2 OS human osteosarcoma cells by suppression of matrix metalloproteinase-2 and -9,” BioMedicine, vol. 3, no. 3, pp. 120–129, 2013.

C.-Y. Chen and C. Y.-C. Chen, “Insights into designing the dual-targeted HER2/HSP90 inhibitors,” Journal of Molecular Graphics & Modelling, vol. 29, no. 1, pp. 21–31, 2010.

S.-C. Yang, S.-S. Chang, H.-Y. Chen, and C. Y.-C. Chen, “Identification of potent EGFR inhibitors from TCM Database@Taiwan,” PLoS Computational Biology, vol. 7, no. 10, Article ID e1002189, 2011.

Y. A. Tsou, K. C. Chen, H. C. Lin, S. S. Chang, and C. Y. C. Chen, “Uroporphyrinogen decarboxylase as a potential target for specific components of traditional Chinese medicine: a virtual screening and molecular dynamics study,” PLoS ONE, vol. 7, no. 11, 2012.

Y. A. Tsou, K. C. Chen, S. S. Chang, Y. R. Wen, and C. Y. Chen, “A possible strategy against head and neck cancer: in silico investigation of three-in-one inhibitors,” Journal of Biomolecular Structure & Dynamics, vol. 31, no. 12, pp. 1358–1369, 2013.

W. I. Tou, S. S. Chang, C. C. Lee, and C. Y. C. Chen, “Drug design for neuropathic pain regulation from traditional Chinese medicine,” Scientific Reports, vol. 3, article 844, 2013.

K.-C. Chen, K.-W. Chang, H.-Y. Chen, and C. Y.-C. Chen, “Traditional Chinese medicine, a solution for reducing dual stroke risk factors at once?” Molecular BioSystems, vol. 7, no. 9, pp. 2711–2719, 2011.

T.-T. Chang, K.-C. Chen, K.-W. Chang et al., “In silico pharmacology suggests ginger extracts may reduce stroke risks,” Molecular BioSystems, vol. 7, no. 9, pp. 2702–2710, 2011.

H. J. Huang, Y. R. Jian, and C. Y. Chen, “Traditional Chinese medicine application in HIV: an in silico study,” Journal of Biomolecular Structure & Dynamics, vol. 32, no. 1, pp. 1–12, 2014.

S.-S. Chang, H.-J. Huang, and C. Y.-C. Chen, “Two birds with one stone? Possible dual-targeting H1N1 inhibitors from traditional Chinese medicine,” PLoS Computational Biology, vol. 7, no. 12, Article ID e1002315, 2011.

K. C. Chen, S. S. Chang, H. J. Huang, T. L. Lin, Y. J. Wu, and C. Y. Chen, “Three-in-one agonists for PPAR-alpha, PPAR-gamma, and PPAR-delta from traditional Chinese medicine,” Journal of Biomolecular Structure & Dynamics, vol. 30, no. 6, pp. 662–683, 2012.

C. Y.-C. Chen, “TCM Database@Taiwan: the world's largest traditional Chinese medicine database for drug screening in silico,” PLoS ONE, vol. 6, no. 1, Article ID e15939, 2011.

C. A. Lipinski, F. Lombardo, B. W. Dominy, and P. J. Feeney, “Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings,” Advanced Drug Delivery Reviews, vol. 46, no. 1–3, pp. 3–26, 2001.

Y. Li, Z. Wang, N. Furukawa et al., “T2384, a novel antidiabetic agent with unique peroxisome proliferator-activated receptor γ binding properties,” The Journal of Biological Chemistry, vol. 283, no. 14, pp. 9168–9176, 2008.

D. Rogers and A. J. Hopfinger, “Application of genetic function approximation to quantitative structure-activity relationships and quantitative structure-property relationships,” Journal of Chemical Information and Computer Sciences, vol. 34, no. 4, pp. 854–866, 1994.

C. Chang and C.-J. Lin, “LIBSVM: a library for support vector machines,” ACM Transactions on Intelligent Systems and Technology, vol. 2, no. 3, pp. 1–27, 2011.

C. M. Venkatachalam, X. Jiang, T. Oldfield, and M. Waldman, “LigandFit: a novel method for the shape-directed rapid docking of ligands to protein active sites,” Journal of Molecular Graphics & Modelling, vol. 21, no. 4, pp. 289–307, 2003.

B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus, “CHARMM: a program for macromolecular energy minimization and dynamics calculations,” Journal of Computational Chemistry, vol. 4, pp. 187–217, 1983.

A. C. Wallace, R. A. Laskowski, and J. M. Thornton, “LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions,” Protein Engineering, vol. 8, no. 2, pp. 127–134, 1995.

B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, “GRGMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation,” Journal of Chemical Theory and Computation, vol. 4, no. 3, pp. 435–447, 2008.

V. Zoete, M. A. Cuendet, A. Grosdidier, and O. Michielin, “SwissParam: a fast force field generation tool for small organic molecules,” Journal of Computational Chemistry, vol. 32, no. 11, pp. 2359–2368, 2011.

R. Fletcher, Optimization, Academic Press, New York, NY, USA, 1969.

D. Bonchev, Information Theoretic Indices for Characterization of Chemical Structures, Research Studies Press, 1983.

L. B. Kier and L. H. Hall, Molecular Connectivity in Chemistry and Drug Research, Academic Press, New York, NY, USA, 1976.

A. T. Balaban, “Highly discriminating distance-based topological index,” Chemical Physics Letters, vol. 89, no. 5, pp. 399–404, 1982.