BioMed Research International

Volume 2014 (2014), Article ID 741703, 14 pages

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

Research Article

In Silico Design of BACE1 Inhibitor for Alzheimer’s Disease by Traditional Chinese Medicine

Hung-Jin Huang,¹ Cheng-Chun Lee,² and Calvin Yu-Chian Chen2,3

¹Department of Chinese Pharmaceutical Sciences and Chinese Medicine Resources, College 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 14 February 2014; Accepted 16 February 2014; Published 8 May 2014

Academic Editor: Chung Y. Hsu

Copyright © 2014 Hung-Jin Huang et al. 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 -site APP cleaving enzyme 1 (BACE1) is an important target for causing Alzheimer’s disease (AD), due to the brain deposition peptide amyloid beta (A) require cleavages of amyloid precursor protein (APP) by BACE1 and -secretase, but treatments of AD still have side effect in recent therapy. This study utilizes the world largest traditional Chinese medicine (TCM) database and database screening to provide potential BACE1 inhibited compound. Molecular dynamics (MD) simulation was carried out to observe the dynamics structure after ligand binding. We found that Triptofordin B1 has less toxicity than pyrimidine analogue, which has more potent binding affinity with BACE1. For trajectory analysis, all conformations are tending to be stable during 5000 ps simulation time. In dynamic protein validation, the residues of binding region are still stable after MD simulation. For snapshot comparison, we found that Triptofordin B1 could reduce the binding cavity; the results reveal that Triptofordin B1 could bind to BACE1 and better than control, which could be used as potential lead drug to design novel BACE1 inhibitor for AD therapy.

1. Introduction

Alzheimer’s disease (AD) is a progressive neurological disease of the central nervous system (CNS) that affects aging patients in the world [1–3]; the causes of AD are not well understood; recent studies indicate that the progression is associated with plaques accumulation and tau protein in the form of neurofibrillary tangles in the cortical region of brain [4–6]. The amyloid hypothesis indicated that amyloid is the initial cause of AD disease contributing to plaques accumulation; one of AD hallmarks is an aggregation of amyloid (A) leading to deposition of -amyloid in the brain [7]. In A reducing approaches, numerous studies demonstrate that amyloid vaccine can remove the amyloid plaques from the brains of the mice and reverse cognitive impairment [8–11], but in human clinical trials, the immunotherapy has side effects during the process of treatment, including autoimmunity [12] and high incidence of meningoencephalitis [13]; clearance of A deposition still has problems for developing AD therapy. Hence, we focus on disrupt formation of A from amyloid precursor protein (APP), cleavage by enzymes for AD prevention. The -secretase is also called BACE1 (-site amyloid precursor protein cleaving enzyme 1), which is an important enzyme in development of AD pathology. BACE1 cleaves transmembrane APP between residues 671 and 672, and carboxy-terminal fragment of APP is cleaved by -secretase, facilitating intramembrane proteolysis by the presenilin 1 (PSEN1) and presenilin 2 (PSEN2) [14, 15]. Subsequently the small 4 kilodalton of amyloid-A1-40 and A1-42 is generated by sequential and -secretase cleavage of APP. Hence, the BACE1 has been recognized as a drug target for curing AD in many studies [16–18].

In this study, in order to design potential lead drugs for BACE1 inhibitor from nature products, computer-aided drug design (CADD) was employed to this research [19, 20], which includes molecular simulation and web server calculation [21, 22]. The nature products from world largest TCM database (TCM Database@Taiwan) were used to investigate more safety drugs [23], because TCM has been wildly used in clinical therapy for thousand years and also used in developing potential illness therapies [24–26]. The approach of CADD combined with TCM database has been wildly used to design new drugs successfully in many cases, including type II diabetes [27], neurotropic pain [28], head and neck cancer [29], hypertension [30], influenza [31–34], inflammation [35], breast cancer [36, 37], HIV virus [38], neuroprotection [39, 40], insomnia [41], erectile dysfunction [42], stroke [43–45], and weight loss [46, 47]. For drug targets, CADD should depend on risk factors study [48–51] and theories [52] to identify direction of each research. Hence, we present small molecular from world largest TCM database to analyze potential nature products by docking study between protein and small TCM compounds. Because BACE1 inhibitors are considered blood brain barrier (BBB) permeability, we also utilize ADMET prediction to evaluate the screening results from docking studies. Besides, we also employed molecular dynamics (MD) simulations to construct dynamic structure of BACE1 with docked ligands, observing the conformation changes over all simulation times.

2. Materials and Methods

2.1. Small Molecules and Protein Structure Preparation

The total numbers of TCM compounds from TCM Database@Taiwan were 61,000, and we employed the TCM compounds to search potent ligand as BACE1 inhibitor by docking study. We further used ADMET prediction and Lipinski’s rule of five [53, 54] to estimate drug-likeness of the TCM compounds from docking results; these rules make them a likely oral drug in the human body. For ADMET prediction, we based on BBB penetration, CYP2D6 inhibition, and hepatotoxicity to analyze all docked ligands. The crystal structure of BACE1 was taken from PDB database (PDB code: 4JPE) [55]; the missing atoms and loops were corrected by Prepare Protein module under Accelrys Discovery Studio 2.5.5.9350 (DS 2.5) [56]; residues of BACE1 were protonated in pH 7.4 condition. We also used PONDR-FIT [57] to evaluate unfolded regions on BACE1 sequence for structure validation.

2.2. Docking Study

The volume of BACE1 inhibitor (1M7) in crystal structure of BACE1 was defined as binding site for screening TCM compounds through protein-ligand interaction; different poses of TCM compound were generated by Monte-Carlo techniques; docking study was performed by LigandFit module within DS 2.5. We utilized CHARMm force field [58] to minimize the conformation of each ligand. The energy function is as follows:

Minimization of each docking pose executes 1000 steps of Steepest Descent with an RMS gradient tolerance of 3 and followed by conjugate gradient. The generated conformation of ligands was docked into the defined binding site of BACE1; the ligand poses were calculated by various scoring functions including -PLP1, -PLP2, and -PMF.

2.3. Molecular Dynamics Simulation

The molecular dynamic simulation was carried out by GROMACS 4.5.5 package [59] to simulate the dynamic structure of BACE1 with docked compounds. We utilize charmm27 force field for the simulation system [60]. The distance between the edge of box and protein was set to 1.2 nm. Each protein-ligand system was placed in cubic cell containing water molecular by TIP3P model. Nonbonded interactions include repulsion, dispersion, and Coulomb terms. The repulsion and dispersion terms involve Lennard-Jones interaction [61] and Buckingham potential [62]; the cut-off distance of define van der Waals (VDW) residues was set to 1.4 nm. Long-range electrostatic forces were performed using the PME method [63, 64].

The equation of Lennard-Jones interaction is as follows:

The Buckingham potential is defined as

Topology files and parameters of small compounds in protein-ligand complexes were generated for GROMACS simulation by SwissParam web server [65]. Bonds lengths were constrained by the linear constraint solver (LINCS) algorithm. Na⁺ and Cl⁻ ion were randomly replaced with water molecular to neutralize the simulation systems, and the concentration was set as 0.145 M in solvent system. The energy minimization was used to stabilize the solvent system by Steepest Descent algorithm with 5,000 steps, the follow by equilibration performed under position restraints to equilibrated water molecular in the protein for 1 ns under constant temperature dynamics (NVT type) conditions. In final step, production running for 5000 ps under constant pressure and temperature dynamics (NPT type); all of the temperature simulation system was under 310 K condition. MD conformations are sampled every 20 ps and all frames are analyzed under GROMACS 4.5.5.

3. Results and Discussion

3.1. Docking Results

We utilize PONDR-FIT [57] to understand the amino acids on binding region (GLN60, GLY61, ASP80, ILE158, ILE166, ASP276, GLY278, and THR279) of BACE1 are not disorder structure (Figure 1), and the values of disorder disposition are below 0.5, which indicate that the binding site of BACE1 is order structure, and the ligands binding may not affected by protein structure [66, 67]. For docking analysis, we based on -PLP1, -PLP2, and -PMF to evaluate the docking pose of traditional Chinese medicine (TCM) compounds. From scoring analysis, pyrimidine analogue R-50 (1M7) was regarded as control for comparison, which is synthesis BACE1 inhibitor from Hunt’s study [55]. Top candidates with higher values of scores than 1M7 are shown in Table 1. For ADMET evaluation, all TCM candidates have no CYP2D6 inhibited and hepatotoxicity, suggesting that CYP2D6 may not be affected by these ligands in liver. The 1M7 has hepatotoxicity in ADMET analysis, indicating that our TCM candidates are safer than control. All docked ligands are ranked by -PMF score, due to the prediction of blood-brain barrier (BBB) penetration showing Diterpenoid EF-D with no penetration ability (BBB level = 4); the Triptofordin B1 has -PMF score (194.61) and medium penetration (BBB level = 2), which is better than 1M7 because of the low penetration (BBB level = 3) and low binding score (-PMF = 119.39). Triptofordin B1 is available in Tripterygium wilfordii; the herb extraction has therapeutic effect for SAMP8 mice with AD disease [68]. So we selected Triptofordin B1 for further studies; the chemical scaffolds of TCM candidates and 1M7 are shown in Figure 2. Docking pose of Triptofordin B1 displayed pi-pi interaction with TYR119; close residues include ASP80 and ASP276 (Figure 3(a)). 1M7 binding pose has H-bond with ASP80 and ASP276, but there is no pi interaction presented between residue and ligand. The data reveal that Triptofordin B1 has similar binding position with 1M7 and displayed stronger chemical interaction in BACE1 binding site. In further study, we utilized MD simulation to perform dynamic protein-ligand complexes for variation analysis.

Table 1: Top ten candidates and control.

Name
-PLP1
-PLP2
-PMF
BBB Level
CYP2D6
Hepatotoxicity

Diterpenoid EF-D
79.04
75.67
195.32
4
0
0

Triptofordin B1
68.44
62.28
194.61
2
0
0

Shionoside C
71.67
69.36
193.84
4
0
0

Jangomolide
72.27
67.01
187.36
3
0
0

Vibsanin W
77.98
76.62
184.31
4
0
0

2α,6α-Dihydroxybetulinic acid
59.52
57.72
183.93
2
0
0

Benzoylramanone
63.59
61.05
183.67
2
0
0

Pseurata D
63.28
62.87
180.42
4
0
0

Vibsanin I
78.18
72.43
179.15
4
0
0

1M7
70.70
52.10
119.39
3
0
1

Control.

BBB level (blood brain barrier): high penetration = 1; medium penetration = 2; low penetration = 3; undefined penetration = 4.

CYP2D6: noninhibitor = 0; Inhibitor = 1.

Hepatotoxicity: Non-inhibitor = 0; inhibitor = 1.

Figure 1: Disorder analysis of sequence of BACE1 from result of PONDR-FIT prediction; the value of disorder disposition above 0.5 indicate disorder residues.

Figure 2: Chemical scaffolds of 1M7 (control) and Triptofordin B1.

(a)

(b)

Figure 3: The docking poses of small compounds: (a) Triptofordin B1; (b) 1M7. Small compound and amino acids are colored in green and yellow, respectively.

3.2. Stability Analysis

Structure of BACE1 with docked ligands includes Triptofordin B1 and 1M7 that were carried out by MD simulation, and we use protein structure of BACE1 with no ligand (Apoprotein) for comparison. The analysis result of protein root mean square deviation (RMSD) and radius of gyration (Rg) is shown in Figure 4. 1M7 displayed fluctuation from 500 to 4500 ps and was stable at 0.3 nm of protein RMSD. Triptofordin B1 and Apoprotein show similar trends; the protein RMSD remained stable in the region of 0.3 nm. The radius of gyration (Rg) analysis shows that the compactness of BACE1 with each ligand is less than the Apoprotein structure, because of the docked ligand combined with BACE1. From 3000 to 5000 ps of Rg analysis, the structure tends to be stable around 0.4 nm.

(a)

(b)

Figure 4: Plots of (a) protein RMSD and (b) radius of gyration from BACE1 during 5000 ps simulation time.

We further analyzed RMSD of each small molecular during MD simulation (Figure 5); ligand RMSD of Triptofordin B1 and 1M7 increases large fluctuation at 2000 ps; the value of ligand RMSD increased from 0.04 to 0.10 nm. Interestingly, 1M7 is decreased from 0.10 nm 0.04 nm after 4500 ps; this finding suggests that the region of 2000 to 4000 ps should be used to analyze the conformation of ligand binding. For total energy analysis, there significant increased values were observed at initial simulation time (Figure 6); the total energy is remained around  kJ/mol for 1M7 and Apoprotein; the Triptofordin B1 was stable at  kJ/mol. These results suggest that all structures of the complexes remained constant after initial simulation time; there is no significant fluctuation among all BACE1 structures.

Figure 5: Plot of ligand RMSD values from BACE1 with docked ligands among 5000 ps simulation times.

(a)  

(b)

(c)

Figure 6: Total energy of BACE1 complexes: (a) Triptofordin B1; (b) 1M7; (c) Apoprotein among 5000 ps simulation times.

3.3. Residues Fluctuation Analysis on the Binding Region

We using root mean squared fluctuation (RMSF) to analyze the fluctuation of residues on protein binding site; the binding region (GLN60, GLY61, ASP80, ILE158, ILE166, ASP276, GLY278, and THR279) shows small flexibility (Figure 7). The largest fluctuation is observed from 425 to 450 residues, because these regions are far away from the binding site, indicating that the flexible amino acids do not affect protein-ligand interaction during MD simulation.

(a)

(b)

(c)

Figure 7: RMSF values of each residue of BACE1 with docked ligand: (a) Triptofordin B1; (b) 1M7 during 5000 ps simulation time.

According to DSSP analysis, the number of helix and beta-sheet remained 150 and 100, respectively (Figure 8); the other secondary type also revealed no distinct changes. Besides, the distance for pair of each residue has no missing plot among all BACE1 structures during 5000 ps, (Figure 9). The results show that structure of BACE1 remained constant during all MD simulations.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 8: DSSP analysis of BACE1 complexes: (a) Triptofordin B1; (b) 1M7; (c) Apoprotein. Number of residues of secondary structure: (d) Triptofordin B1; (e) 1M7 (f) Apoprotein.

(a)

(b)

(c)

Figure 9: Matrix of smallest between each pair of amino acids in BACE1 complexes: (a) Triptofordin B1; (b) 1M7; (c) Apoprotein.

3.4. Movement of Each Ligand Analysis

The mobility of each ligand was analyzed by mean square displacement (MSD) (Figure 10); Triptofordin B1 increased MSD values to 0.3 nm at 2500 ps, and stabilizes until 4000 ps. 1M7 was stable below 0.1 nm and decreased MSD value at 4500 ps. In final simulation after 4500 ps, Triptofordin B1 further increased MSD values to 0.45 nm and tends to be stable to the end time. Here, we further analyze the distance between BACE1 and each ligand among 5000 ps (Figure 11). The distance between 1M7 and BACE1 displayed 1.00 nm before 2000 ps, but Triptofordin B1 increased to 1.50 nm from 2000 to 3500 ps, and the other wild increased distance occur from 4000 to 5000 ps. These results comparing with MSD analysis;the region of 200 to 3500 ps has significant change during dynamics simulation; in the next analysis, we focus on these regions of simulation time for further studies.

Figure 10: Mean square displacement (MSD) of different ligands during 5000 ps simulation times; the value of MSD indicates migration of ligands (Triptofordin B1 and 1M7) from initial site.

Figure 11: Distance between centers of mass of BACE1 and each ligand during 5000 ps simulation times.

3.5. Clustering Analysis for Snapshot Observing

In order to understand the most stable structure during the entire MD simulation for understanding the movement of BACE1, all frames of MD simulation were clustered into different subgroups (Figure 12); the similar MD conformations were grouped into the same cluster. For clustering results, each last group includes last 1000 ps (from 4000 to 5000 ps); hence, we selected the middle flams from each last group for further analysis from all MD complexes (Table 2). Before observing all snapshots from middle frames of last clustering group, we also calculate the distance of H-bonds for each ligand among all simulation times (Figure 13); GLN121 showed decreased distance after 4000 ps for Triptofordin B1 (Figure 13(a)). ASP80 and ASP276 remain revealed low distance with 1M7 (Figure 13(b)), suggesting that GLN121, ASP80, and ASP276 are essential amino acid for ligand binding. In snapshot analysis, we found that Triptofordin B1 could reduce the binding site, because the GLN121 has significant change, and presenting pi interaction with TRP163 (Figure 14), and in ligand channel analysis (Figure 15), we can see that the predicted channel of Triptofordin B1 is shorter than 1M7 and Apoprotein, suggesting that Triptofordin B1 could bind to BACE1 better than 1M7.

Table 2: The middle conformation in each cluster from all MD conformations.

Cluster
Time of middle frame (ps)

Triptofordin B1
1M7
Apoprotein

1
1640
20
1680

2
3540
400
3180

3
3820
1380
4200

4
4300
2340
—

5
—
3760
—

(a)

(b)

(c)

Figure 12: Clustering analyses among 5000 ps simulation times. (a) Triptofordin B1; (b) 1M7 (c) Apoprotein.

(a)

(b)

Figure 13: H-bond distance between residues atoms of BACE1 and ligands during 5000 ps simulation times: (a) Triptofordin B1; (b) 1M7.

(a)

(b)

(c)

Figure 14: The middle structure from each final clustering group: (a) Triptofordin B1; (b) 1M7 (c) Apoprotein.

(a)

(b)

(c)

Figure 15: Predictions of ligand channel of (a) Triptofordin B1, (b) 1M7, and (c) Apoprotein by CAVER 3.0. The similar channels were collected in a group with same color.

4. Conclusion

For ADMET analysis, Triptofordin B1 has more penetration than 1M7 and less toxicity, because 1M7 has hepatotoxicity in ADMET prediction. Three scoring functions, -PLP1, -PLP2, and –PMF, are higher than control. The structure of BACE1 analysis shows that the binding residues have less fluctuation after MD simulations, indicating the each ligand is not affected by protein residues. In migration analysis for Triptofordin B1 and 1M7, the stable region displayed from 3000 to 4000 ps; we utilize clustering analysis to observe this period simulation time. Triptofordin B1 could reduce the binding cavity of BACE1; the results reveal that Triptofordin B1 could bind to BACE1 and better than 1M7, which could be used as potential lead drug to design novel BACE1 inhibitor for AD therapy

Conflict of Interests

The authors declared that there is no conflict of interests.

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 (ASIA100-CMU-2, ASIA101-CMU-2, 102-ASIA-07), and China Medical University Hospital (DMR-103-058, DMR-103-001, 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) and Taiwan Department of Health Cancer Research Center of Excellence (MOHW103-TD-B-111-03).

References

S. H. Mokhtar, M. M. Bakhuraysah, D. S. Cram, and S. Petratos, “The beta-amyloid protein of alzheimer's disease: communication breakdown by modifying the neuronal cytoskeleton,” International Journal of Alzheimer’s Disease, vol. 2013, Article ID 910502, 15 pages, 2013.

M. O. Grimm, J. Mett, C. P. Stahlmann, et al., “Neprilysin and abeta clearance: impact of the APP intracellular domain in NEP regulation and implications in Alzheimer's disease,” Frontiers in Aging Neuroscience, vol. 5, article 98, 2013.

J. R. Petrella, “Neuroimaging and the search for a cure for Alzheimer disease,” Radiology, vol. 269, no. 3, pp. 671–691, 2013.

J. D. Weinstein, E. R. Gonzalez, R. D. Egleton, and D. A. Hunt, “A paradigm shift for evaluating pharmacotherapy for Alzheimer's disease: the 10-patient screening protocol,” The Consultant Pharmacist, vol. 28, no. 7, pp. 443–454, 2013.

L. Ghezzi, E. Scarpini, and D. Galimberti, “Disease-modifying drugs in Alzheimer's disease,” Drug Design, Development and Therapy, vol. 7, pp. 1471–1479, 2013.

C. Duyckaerts, B. Delatour, and M.-C. Potier, “Classification and basic pathology of Alzheimer disease,” Acta Neuropathologica, vol. 118, no. 1, pp. 5–36, 2009.

H. W. Querfurth, “Calcium ionophore increases amyloid β peptide production by cultured cells,” Biochemistry, vol. 33, no. 15, pp. 4550–4561, 1994.

D. Morgan, D. M. Diamond, P. E. Gottschall et al., “A β peptide vaccination prevents memory loss in an animal model of Alzheimer's disease,” Nature, vol. 408, no. 6815, pp. 982–985, 2000.

J. McLaurin, R. Cecal, M. E. Kierstead et al., “Therapeutically effective antibodies against amyloid-β peptide target amyloid-β residues 4-10 and inhibit cytotoxicity and fibrillogenesis,” Nature Medicine, vol. 8, no. 11, pp. 1263–1269, 2002.

B. Philippe, J.-P. Brion, A.-F. Macq, and J.-N. Octave, “A new monoclonal antibody against the anionic domain of the amyloid precursor protein of Alzheimer's disease,” NeuroReport, vol. 5, no. 3, pp. 289–292, 1993.

D. M. Wilcock, A. Rojiani, A. Rosenthal et al., “Passive immunotherapy against Aβ in aged APP-transgenic mice reverses cognitive deficits and depletes parenchymal amyloid deposits in spite of increased vascular amyloid and microhemorrhage,” Journal of Neuroinflammation, vol. 1, no. 1, article 24, 2004.

T. Wisniewski and U. Konietzko, “Amyloid-β immunisation for Alzheimer's disease,” The Lancet Neurology, vol. 7, no. 9, pp. 805–811, 2008.

S. Rasool, R. Albay, H. Martinez-Coria, et al., “Vaccination with a non-human random sequence amyloid oligomer mimic results in improved cognitive function and reduced plaque deposition and micro hemorrhage in Tg2576 mice,” Molecular Neurodegeneration, vol. 7, article 37, 2012.

D. J. Selkoe and M. S. Wolfe, “Presenilin: running with Scissors in the membrane,” Cell, vol. 131, no. 2, pp. 215–221, 2007.

M.-C. Yin, “Anti-glycative potential of triterpenes: a mini-review,” BioMedicine, vol. 2, no. 1, pp. 2–9, 2012.

S. Butini, E. Gabellieri, M. Brindisi, et al., “A stereoselective approach to peptidomimetic BACE1 inhibitors,” European Journal of Medicinal Chemistry, vol. 70, pp. 233–247, 2013.

Y. Zou, L. Xu, W. Chen, et al., “Discovery of pyrazole as C-terminus of selective BACE1 inhibitors,” European Journal of Medicinal Chemistry, vol. 68, pp. 270–283, 2013.

S. Eketjall, J. Janson, F. Jeppsson, et al., “AZ-4217: a high potency BACE inhibitor displaying acute central efficacy in different in vivo models and reduced amyloid deposition in Tg2576 mice,” The Journal of Neuroscience, vol. 33, no. 24, pp. 10075–10084, 2013.

H.-J. Huang, H. W. Yu, C.-Y. Chen et al., “Current developments of computer-aided drug design,” Journal of the Taiwan Institute of Chemical Engineers, vol. 41, no. 6, pp. 623–635, 2010.

C. Y. Chen, “A novel integrated framework and improved methodology of computer-aided drug design,” Current Topics in Medicinal Chemistry, vol. 13, no. 9, pp. 965–988, 2013.

T.-Y. Tsai, K.-W. Chang, and C. Y.-C. Chen, “IScreen: world's first cloud-computing web server for virtual screening and de novo drug design based on TCM database@Taiwan,” Journal of Computer-Aided Molecular Design, vol. 25, no. 6, pp. 525–531, 2011.

K.-W. Chang, T.-Y. Tsai, K.-C. Chen et al., “iSMART: an integrated cloud computing web server for traditional Chinese medicine for online virtual screening, de novo evolution and drug design,” Journal of Biomolecular Structure and Dynamics, vol. 29, no. 1, pp. 243–250, 2011.

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.

P.-C. Lin, P.-Y. Liu, S.-Z. Lin, and H.-J. Harn, “Angelica sinensis: a Chinese herb for brain cancer therapy,” BioMedicine, vol. 2, no. 1, pp. 30–35, 2012.

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.

S. C. Hsu and J. G. Chung, “Anticancer potential of emodin,” BioMedicine, vol. 2, no. 3, pp. 108–116, 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 and Dynamics, vol. 31, no. 11, pp. 1219–1235, 2013.

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

Y. A. Tsou, K. C. Chen, H. C. Lin, S. S. Chang, and C. Y. 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, Article ID e50087, 2012.

W. I. Tou and C. Y.-C. Chen, “In silico investigation of potential Src kinase ligands from traditional chinese medicine,” PLoS ONE, vol. 7, no. 3, Article ID e33728, 2012.

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.

C.-H. Lin, T.-T. Chang, M.-F. Sun et al., “Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database,” Journal of Biomolecular Structure and Dynamics, vol. 28, no. 4, pp. 471–482, 2011.

T.-T. Chang, M.-F. Sun, H.-Y. Chen et al., “Screening from the world's largest TCM database against H1N1 virus,” Journal of Biomolecular Structure and Dynamics, vol. 28, no. 5, pp. 773–786, 2011.

S.-S. Chang, H.-J. Huang, and C. Y.-C. Chen, “High performance screening, structural and molecular dynamics analysis to identify H1 inhibitors from TCM Database@Taiwan,” Molecular BioSystems, vol. 7, no. 12, pp. 3366–3374, 2011.

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

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.

S.-C. Yang, S.-S. Chang, and C. Y.-C. Chen, “Identifying HER2 inhibitors from natural products database,” PLoS ONE, vol. 6, no. 12, Article ID e28793, 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.

K. C. Chen, Y. R. Jian, M. F. Sun, et al., “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.

H.-J. Huang, C.-C. Lee, and C. Y.-C. Chen, “Pharmacological chaperone design for reducing risk factor of Parkinson's Disease from Traditional Chinese Medicine,” Evidence-Based Complementary and Alternative Medicine, vol. 2014, Article ID 830490, 12 pages, 2014.

H. Y. Chen, S. S. Chang, Y. C. Chan, and C. Y. Chen, “Discovery of novel insomnia leads from screening traditional Chinese medicine database,” Journal of Biomolecular Structure and Dynamics, vol. 32, no. 5, pp. 776–791, 2014.

C. Y.-C. Chen, “Computational screening and design of Traditional Chinese Medicine (TCM) to block phosphodiesterase-5,” Journal of Molecular Graphics and Modelling, vol. 28, no. 3, pp. 261–269, 2009.

K.-C. Chen and C. Yu-Chian Chen, “Stroke prevention by traditional Chinese medicine? A genetic algorithm, support vector machine and molecular dynamics approach,” Soft Matter, vol. 7, no. 8, pp. 4001–4008, 2011.

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.

K. Y. Chen, S. S. Chang, and C. Y. Chen, “In silico identification of potent pancreatic triacylglycerol lipase inhibitors from traditional Chinese medicine,” PLoS ONE, vol. 7, no. 9, Article ID e43932, 2012.

K. C. Chen, S. S. Chang, H. J. Huang, et al., “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.

I. C. Chou, W. D. Lin, C. H. Wang, et al., “Möbius syndrome in a male with XX/XY mosaicism,” BioMedicine, vol. 3, no. 2, pp. 102–104, 2013.

W.-Y. Lin, H.-P. Liu, J.-S. Chang, et al., “Genetic variations within the PSORS1 region affect Kawasaki disease development and coronary artery aneurysm formation,” BioMedicine, vol. 3, no. 2, pp. 73–81, 2013.

C. C. Lee, C. H. Tsai, L. Wan, et al., “Increased incidence of Parkinsonism among Chinese with α-glucosidase mutation in central Taiwan,” BioMedicine, vol. 3, no. 2, pp. 92–94, 2013.

F. J. Tsai, “Rare diseases: a mysterious puzzle,” BioMedicine, vol. 3, no. 2, article 65, 2013.

C. Y.-C. Chen, “Weighted equation and rules—a novel concept for evaluating protein-ligand interaction,” Journal of Biomolecular Structure and Dynamics, vol. 27, no. 3, pp. 271–282, 2009.

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.

C. A. Lipinski, “Lead- and drug-like compounds: the rule-of-five revolution,” Drug Discovery Today: Technologies, vol. 1, no. 4, pp. 337–341, 2004.

K. W. Hunt, A. W. Cook, R. J. Watts, et al., “Spirocyclic beta-site amyloid precursor protein cleaving enzyme 1 (BACE1) inhibitors: from hit to lowering of cerebrospinal fluid (CSF) amyloid beta in a higher species,” Journal of Medicinal Chemistry, vol. 56, no. 8, pp. 3379–3403, 2013.

Accelrys, Discovery Studio Client v2.5, Accelrys, San Diego, Calif, USA, 2009.

B. Xue, R. L. Dunbrack, R. W. Williams, A. K. Dunker, and V. N. Uversky, “PONDR-FIT: a meta-predictor of intrinsically disordered amino acids,” Biochimica et Biophysica Acta: Proteins and Proteomics, vol. 1804, no. 4, pp. 996–1010, 2010.

B. R. Brooks, C. L. Brooks III, A. D. Mackerell Jr. et al., “CHARMM: the biomolecular simulation program,” Journal of Computational Chemistry, vol. 30, no. 10, pp. 1545–1614, 2009.

S. Pronk, S. Pall, R. Schulz, et al., “GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit,” Bioinformatics, vol. 29, no. 7, pp. 845–854, 2013.

U. D. Priyakumar and A. D. MacKerell, “Comparison of the CHARMM27, AMBER4. 1 and BMS nucleic acid force fields via free energy calculations of base flipping,” Abstracts of Papers of the American Chemical Society, vol. 230, pp. U1391–U1392, 2005.

J. E. Jones, “On the determination of molecular fields. II. From the equation of state of a gas,” Proceedings of the Royal Society of London Series A, vol. 106, no. 738, pp. 463–477, 1924.

R. A. Buckingham, “The classical equation of state of gaseous helium, neon and argon,” Proceedings of the Royal Society of London A: Mathematical and Physical Sciences, vol. 168, no. 933, pp. 264–283, 1938.

T. Darden, D. York, and L. Pedersen, “Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems,” The Journal of Chemical Physics, vol. 98, no. 12, pp. 10089–10092, 1993.

U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen, “A smooth particle mesh Ewald method,” The Journal of Chemical Physics, vol. 103, no. 19, pp. 8577–8593, 1995.

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.

W. I. Tou and C. Y. Chen, “May disordered protein cause serious drug side effect?” Drug Discovery Today, vol. 19, no. 4, pp. 367–372, 2014.

C. Y. Chen and W. I. Tou, “How to design a drug for the disordered proteins?” Drug Discovery Today, vol. 18, no. 19-20, pp. 910–915, 2013.

J. Nie, M. Zhou, C. Lü et al., “Effects of triptolide on the synaptophysin expression of hippocampal neurons in the AD cellular model,” International Immunopharmacology, vol. 13, no. 2, pp. 175–180, 2012.