BioMed Research International

Volume 2014 (2014), Article ID 769867, 15 pages

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

Research Article

In Silico Investigation of Traditional Chinese Medicine Compounds to Inhibit Human Histone Deacetylase 2 for Patients with Alzheimer’s Disease

Tzu-Chieh Hung,¹ Wen-Yuan Lee,1,2,3 Kuen-Bao Chen,1,2,4 Yueh-Chiu Chan,² Cheng-Chun Lee,² and Calvin Yu-Chian Chen1,2,5,6

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

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

³Department of Neurosurgery, China Medical University Hospital, No. 2, Yude Road, North District, Taichung 40447, Taiwan

⁴Department of Anesthesiology, China Medical University Hospital, Taichung 40447, Taiwan

⁵Research Center for Chinese Medicine & Acupuncture, China Medical University, Taichung 40402, Taiwan

⁶Human Genetic Center, Department of Medical Research, China Medical University Hospital, Taichung 40447, Taiwan

Received 20 February 2014; Accepted 5 March 2014; Published 23 June 2014

Academic Editor: Chung Y. Hsu

Copyright © 2014 Tzu-Chieh Hung 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.

Human histone deacetylase 2 (HDAC2) has been identified as being associated with Alzheimer’s disease (AD), a neuropathic degenerative disease. In this study, we screen the world’s largest Traditional Chinese Medicine (TCM) database for natural compounds that may be useful as lead compounds in the search for inhibitors of HDAC2 function. The technique of molecular docking was employed to select the ten top TCM candidates. We used three prediction models, multiple linear regression (MLR), support vector machine (SVM), and the Bayes network toolbox (BNT), to predict the bioactivity of the TCM candidates. Molecular dynamics simulation provides the protein-ligand interactions of compounds. The bioactivity predictions of pIC50 values suggest that the TCM candidatesm, (−)-Bontl ferulate, monomethylcurcumin, and ningposides C, have a greater effect on HDAC2 inhibition. The structure variation caused by the hydrogen bonds and hydrophobic interactions between protein-ligand interactions indicates that these compounds have an inhibitory effect on the protein.

1. Introduction

Alzheimer’s disease (AD) is a neuropathic degenerative disease in which patients will gradually suffer a loss of memory, language, intellect, motor action, and even life. In 2010, it was reported that about 36 million people worldwide suffered from AD [1]. The medical cost of this condition was predicted to be approximately 604 billion USD in 2010 [2]. This huge medical expense becomes a great social burden to an aging society.

Recently, it has been found that Tau protein [3], amyloid-β peptides [4], and human histone deacetylase (HDAC) are major factors in the causation of AD [5]. Human histone deacetylase 2 (HDAC2) is the protein expressed by HDAC2 gene. Some reports have pointed out that HDAC2 is over expressed in AD patients and that this gene negatively regulates memory [6–10]. There are also some references indicating that blocking the HDAC2 gene could be a treatment for AD; furthermore, it has been shown to decrease amyloid-β peptides in mice [5, 11, 12]. HDACs catalyze the acetyl moiety, removing it from the lysine residues of proteins and regulating the level of protein acetylation [13]. The inhibition of HDAC2 has been identified as a mechanism for treating cancer and developing histone deacetylase inhibitors (HDACi) [14]. As shown above, this inhibition mechanism could also be a model for the treatment of AD [6].

Some HDACi studies have indicated a role for chromatin remodeling increasing histone acetylation and enhancing synaptic plasticity and learning behaviors [15–17]. The clinical application of nonselective HDACi in cancer has shown a range of side effects [18, 19]. Suberoylanilide hydroxamic acid (SAHA or vorinostat) is a potent HDACi. SAHA binds to the active site of HDAC where it acts as a chelator for Zinc [13]. SAHA could cross the blood-brain barrier and decrease amyloid β peptides and treat AD and Huntington’s disease (HD) by changes in histone acetylation in the brain [20–22].

Computer-aided drug design (CADD) is an in silico simulation technique for screening novel drug-candidate compounds by structure and prediction of biological activity. The two major application areas of CADD are structure-based drug design and ligand-based drug design. In comparison with traditional drug design, CADD has the advantages of both greater speed and lower cost. We used CADD for molecular simulation based on structure-based drug design, ligand-based drug design, and molecular dynamics [23–28].

Recently, an understanding of personalized medicine and biomedicine has been attracting more and more attention [29]; this department of knowledge could analyze regional diseases [30], clinical diagnosis cases, and disease associated mutations [31]. Traditional Chinese Medicine (TCM) plays an important role in Asia, especially in China, Taiwan, Korea, and Japan. The TCM Database@Taiwan (http://tcm.cmu.edu.tw/) [32] is the largest Traditional Chinese Medicine database in the world. This database contains 2D chemical structures, 3D chemical structures, bioactivity, and molecular information of 61,000 compounds used in Traditional Chinese Medicine. Since 2011, there have been successful discoveries in novel lead compounds from the TCM Database@Taiwan [33–35], including compounds for the putative treatment of AD [36], Parkinson’s Disease [37], insomnia [38], pigmentary disorders [39], and even antivirals [40–44]. Due to the application system of the website [45] and cloud computing platforms [46], the TCM Database@Taiwan is exceptionally helpful for TCM applications and drug design.

In this study, we screen a possible lead compound against HDAC2 from the TCM Database@Taiwan. We use the computational techniques of docking, screening, and ligand-based methods to predict the bioactivity of the selected ligands. Finally, we apply molecular dynamics (MD) simulation to investigate variation from the protein-ligand interactions that may contribute to the evaluation of the effect of HDAC2 inhibition.

2. Materials and Methods

2.1. Data Set

Because the disorder protein plays an important role in drug design, the protein sequence should be submitted to the Database of Protein Disorder (DisProt, http://www.disprot.org/) for disorder prediction [47]. The result of prediction could help define the character of docking site and the efficacy of drug interaction.

A total of 61,000 TCM compounds were downloaded from the TCM database (http://tcm.cmu.edu.tw/). The HDAC2 (PDB ID: 3MAX) crystal structure was obtained from RCSB protein data bank. Based on the research, the SAHA was used as a control [48]. Accelrys Discovery Studio 2.5 (DS 2.5) was used to perform the molecular simulations.

2.2. Molecular Docking

LigandFit is a receptor-rigid docking algorithm program in Discovery Studio 2.5 (DS 2.5). The docking simulation was performed using LigandFit [49] module to dock SAHA and TCM compounds to HDAC2 in the force field of CHARMm [50]. The compounds downloaded from the TCM Database were docked into the SAHA binding site of HDAC2, which was identified from previous research [48].

2.3. Ligand-Based Prediction

We used the SVM, MLR, and BNT approaches as activity prediction models to predict the bioactivity of the TCM compounds. The SVM model, which utilized a nonlinear mapping technique, was built using LibSVM [51]. The MLR, using linear evaluation, was established by MATLAB. The BNT is a directed graphical model based on the study of Bayes probability and network cluster in MATLAB. Fifty compounds [52] were divided randomly into the predictive models as a training set of ten compounds and test set for the other compounds. These compounds were drawn by ChemBioOffice and modified using the prepare ligand module in DS 2.5. The bioactivity (IC50) data of the compounds were converted to pIC50 by the QSAR program and their molecular descriptors were evaluated by calculate molecular properties from 552 different descriptors in Discovery Studio 2.5. The genetic function approximation (GFA) algorithm is based on a genetic algorithm to search all possible QSAR models and use a square correlation coefficient () to estimate the best representative molecular descriptors [53]. After the above processes, the bioactivity of TCM de novo compounds can be predicted through MLR, SVM, and BNT model by MATLAB and LibSVM.

2.4. Molecular Dynamics Simulation

The ligands must be reprepared by using SwissParam (http://swissparam.ch/) [54] before applying MD simulation based on the reference force field [55] of GROMACS 4.5.5 [56]. The HDAC protein combines with ligands as the complex in the buffer (or solution) simulation box. This cubic box provided a minimum distance of 1.2 Å from the complex and was solvated with the TIP3P water model in which sodium and chloride ions were added to neutralize complex charges. The complex was minimized with the steepest descent method for 5,000 steps. The last structure of minimization was transferred for MD simulation. The calculations for electrostatic interactions were based on the particle-mesh Ewald (PME) method [57]. In the PME method, each time step was 2 fs and there were 20,000,000 iterations. Equilibration under the 10 ps constant temperature (NVT ensemble) was based on the Berendsen weak thermal coupling method. The total simulation time for the MD was 40 ns. MD trajectories, RMSD, energy variations of the complex, and the secondary structure database (DSSP) were analyzed using a series of protocols in GROMACS.

3. Result and Discussion

3.1. Molecular Docking

The results of molecular docking are ranked according to the docking score. We selected the pdb data of these docking poses and removed the data that lacked any of the important amino-acids of HDAC2, such as Tyr29, Met35, Phe114, Leu144, Gly154, Phe155, Asp181, His183, Phe210, Asp269, and Leu276, which were identified by [48]. His142 was identified as the active site in HDAC2 based on Uniprot (http://www.uniprot.org/) data. The top twenty ligands were ranked according to the docking score and used for ligand-based predictions of their bioactivity.

3.2. Ligand-Based Prediction

The three prediction models were established from the relations of the test set and training set through MLR, SVM, and BNT (Figure 1). The square correlation coefficients () of these three models were 0.7773, 0.851, and 0.8439, which indicates confidence in these models. The ligands which achieved the top three docking scores from the compounds, having higher bioactivity than the control, were selected and identified as candidate compounds (Table 1).

Table 1

Name
Dock score
SVM
MLR
BNT

Ningposides C
30.26
6.62067
6.88967
5.657169

Monomethylcurcumin
27.295
7.59112
6.106833
6.471242

(−)-Bontl ferulate
21.203
7.51855
7.773727
5.627211

Yakuchinone B
19.121
7.42844
5.443937
5.297163

Spinacetin
15.077
6.82774
5.493229
5.369466

SAHA*
3.51
6.61107
4.949132
4.500999

Control.

(a)

(b)

(c)

Figure 1: Relation of observed activity (pIC50) and predict activity (pIC50). (a) MLR, (b) SVM, and (c) BNT.

The ligand structure (Figure 2) and the docking poses (Figure 3) of these candidate compounds were prepared for MD. From the above description, the docking poses of the candidate compounds indicated that the ligands had interactions with different important amino-acids in the protein. We determined the ligand interactions with the amino-acids by ligplot v.2.2.25 [58] (Figure 4). From this figure, we found that all the ligands can interact with Gly154, His183, and Phe210 of HDAC2. Phe155 and Leu276 are effected by selected ligands (do not have control). This indicates that the docking site could have interactions with important amino-acids and simulate the compound’s target in HDAC2.

(a)

(b)

(c)

(d)

Figure 2: The 2D structure of control and candidate TCM compounds. (a) Control-saha, (b) (−)-Bontl ferulate, (c) monomethylcurcumin, and (d) ningposides C.

(a)

(b)

(c)

(d)

(e)

Figure 3: The docking poses of HDAC2. (a) HDAC2, (b) control-saha (c) (−)-Bontl ferulate, (d) monomethylcurcumin, and (e) ningposides C.

(a)

(b)

(c)

(d)

Figure 4: Ligplot illustrates the hydrophobic interactions. (a) Control-saha, (b) (−)-Bontl ferulate, (c) monomethylcurcumin, and (d) ningposides C.

3.3. Molecular Dynamics Simulation

After 40 ns MD, the total energy and RMSD of the whole complex were recorded, with these data being drawn by OriginPro 8.5 (Figure 5). The amplitude of energy, as illustrated by the total energy of these complexes, tended to the region of −672000~−674000 kcal/mol. The RMSD of the whole complex presented by these protein and ligand interactions will tend to balance ningposides C. From Figure 5, it can be seen that ligand RMSD1 focused on ligand structure variation during the MD process, and ligand RMSD2 described the ligand variation of the whole complex. We analyzed the ligand based on the comparison between ligand RMSD1 and ligand RMSD2.

(a)

(b)

(c)

(d)

Figure 5: Measures of the MD trajectories. (a) Complex RMSD, (b) ligand RMSD focus on complex, (c) ligand RMSD focus on ligand, and (d) the total energy.

The amplitude of monomethylcurcumin in the ligand RMSD1 was the greatest, between 0.1 and 0.2, while the amplitude in the ligand RMSD2 was smooth. From this result, monomethylcurcumin had a strong interaction with HDAC2, but the complex was stable. Ningposides C had the greatest amplitude in ligand RMSD2, and this impossible variation suggested that it was not only the protein and ligand interaction, but the complex might be separated.

The site was closed to the H-bond produced site and will be twisted in torsion analysis, as seen in Figures 6(a)-(1), 6(b)-(8), 6(c)-(12), and 6(c)-(13). However, the ligand still needs interactions with the protein by H-bond production. Although the main structure of the ligand was stable, Figures 6(c)-(14) and 6(d)-(22) described the large variation of ligand structure during the MD process. This phenomenon may be caused by two hexagonal ring structures limiting the flexibility of ligand variation during the interaction.

(a)

(b)

(c)

(d)

Figure 6: Torsion angles during MD. Torsion angle measure is designated by the number which corresponds to the radar chart. The red, green, and blue lines in the radar chart indicate the angle during docking, 0 ns, and the period of MD.

In the clustering based on RMSD variation (Figure 7), monomethylcurcumin had the fewest groups of these compounds, indicating that the variation was lower than in the others. Ningposides C had the greatest RMSD range and the most groups of these compounds. From this result, we suggest that the complex of ningposides C and HDAC2 was unstable and the interaction of this complex seemed to be incessant. This phenomenon was described after the screening of the MD process which recorded ningposides C drifting away from the docking site at 4.12 ns, then retargeting at 14.8 ns. It drifts away from the docking site again at 16.36 ns and targets again at 21.06 ns. This result suggests that ningposides C might not be an appropriate compound due to the possibility of its moving away during MD.

(a)

(b)

(c)

(d)

Figure 7: Clustering of RMSD. (a) Control-saha, (b) (−)-Bontl ferulate, (c) monomethylcurcumin, and (d) ningposides C. In each upper triangle, the color shows the RMSD difference between times on both the -axis and -axis. In the lower triangle, the red triangle shows the same groups based on similar RMSD.

The DSSP describes the type of protein structure (Figure 8). In this figure, SAHA causes the HDAC2 structure to reduce other types (not comprising α helix, β sheet, and turn structure) but increases the percentage of α helix, β sheet, and turn at least 20 ns interaction. (−)-Bontl ferulate, monomethylcurcumin, and ningposides C make different effect from SAHA in that the other type was increased and the α helix was decreased. This result suggests that these candidate compounds may result from different level effects on the amino-acids of HDAC2.

Figure 8: The DSSP of the protein component. DSSP shows the composition variation of the protein structure during MD.

The disorder shows the protein unfolding region and the RMSF describes the protein variation (Figure 9). The disorder analyzed the complete HDAC2 amino-acid sequence and this figure focuses on the region of HDAC2 structure. In this result, the residue with the larger pick is not disorder region. Thus, the disorder region has weak effect on this investigation. The RMSF recorded the amino-acid variation in 40 ns. The amplitude in RMSF among different ligand interactions was similar and this result indicates that the protein structure variation may be similar. This suggestion may be based on the structure that illustrates their matrix (Figure 10). The highest pick in RMSD was Phe210, which is one of the important amino-acids, and this result may indicate that Phe210 will move away from the ligand to inhibit the protein function.

(a)

(b)

(c)

(d)

(e)

Figure 9: The disorder region prediction and RMSF detection. The disorder region could define the efficacy of docking based on the character of docking site and important region. For these RMSF curve is the calculation of the RMSD average for the whole MD focus on each residue.

(a)

(b)

(c)

(d)

Figure 10: HDAC2 structure variation by ligand iteration. (a) Control-saha, (b) (−)-Bontl ferulate, (c) monomethylcurcumin, and (d) ningposides C.

The variation in structure was calculated from the centroid distance from Met35, Gly154, Phe210, and Leu276, which are four important amino-acids in HDAC2. The centroid distance between Met35 and Phe210 and the centroid distance between Gly154 and Leu276 could help us to investigate the HDAC2 variation (Figure 11). The distance between Met35 and Phe210 in SAHA increased from 13.383 to 24.938. The distance between Met35 and Phe210 changed to 32.454 in (−)-Bontl ferulate, 28.667 in monomethylcurcumin, and 23.149 in ningposides C. There was an obvious increase in this distance when the protein and ligand interacted during the MD. On the other hand, the distance between Gly154 and Leu276 in SAHA increased from 10.955 to 15.335, and this distance increased to 12.570, 13.810, and 16.482 during other ligand interactions. These distance variations the candidate ligands made during interactions were similar to SAHA, a result which can be confirmed from Figure 10.

(a)

(b)

Figure 11: The pathway determination and HDAC2 structure variation with SAHA iteration. (a) Pathway and (b) structure variation.

The pathway and structure variation could help to discuss the protein-ligand interaction. The pathway definition is using caver 3.0 to determine the interpath protein path during MD simulation [59] and structure variation comparing the difference between MD 0 ns and MD 40 ns (Figures 11 to 14).

While protein interacts with SAHA, the pathway for ligand was around docking site (Figure 11(a)) and the pole of HDAC2 will be spread (Figure 11(b)). Similar to SAHA, (−)-Bontl ferulate also makes the same variation in both pathway and structure (Figure 12). The pathway for monomethylcurcumin is defined around binding site based on calculation (Figure 13). One pathway for ningposides C is defined inside the protein colored in red, but that is not reasonable in biology. This situation is caused from the structure variation and then the path which is large enough for calculation. Thus, this unreasonable pathway does not take reference for the analysis of the ligand moving.

(a)

(b)

Figure 12: The pathway determination and HDAC2 structure variation with (−)-Bontl ferulate iteration. (a) Pathway and (b) structure variation.

(a)

(b)

Figure 13: The pathway determination and HDAC2 structure variation with monomethylcurcumin iteration. (a) Pathway and (b) structure variation.

(a)

(b)

Figure 14: The pathway determination and HDAC2 structure variation with ningposides C iteration. (a) Pathway and (b) structure variation.

The structure variation of HDAC2 might make ligand move away from docking site if the force of target is much less. If the ningposides C, selected as candidate compounds, away from binding site is unreasonable for drug design, the MD simulation for the investigation of protein-ligand interaction is more important than docking. This situation could present that MD could confirm protein-ligand interaction by computational biology.

4. Conclusion

The analysis of HDAC2 found that the candidate compounds had a HDAC structure variation similar to SAHA. Because these compounds had better bioactivity than SAHA based on SVM, MLR, and BNT, these compounds had a different effect on the amino-acids and caused changes to the component percentage of HDAC2 structure. Finally, from a comparison of HDAC2 structure variation, component, and protein-ligand interactions, we propose that the candidate compounds are appropriate to inhibit HDAC2, including ningposides C due to its unstable interaction.

Conflict of Interests

The authors confirm that there is no conflict of interests to declare.

Authors’ Contribution

Tzu-Chieh Hung, Wen-Yuan Lee, and Kuen-Bao Chen equally contributed to this study.

Acknowledgments

The research was supported by Grants from the National Science Council of Taiwan (NSC102-2325-B039-001 and NSC102-2221-E-468-027-), Asia University (ASIA100-CMU-2, ASIA101-CMU-2, and 102-ASIA-07), and China Medical University Hospital (DMR-103-058 and DMR-102-051). 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. Finally, the authors’ gratitude goes to Tim Williams, Asia University.

References

C. P. Ferri, R. Sousa, E. Albanese, W. S. Ribeiro, and M. Honyashiki, “Alzheimer’s disease international,” in World Alzheimer Report 2009, M. Prince and J. Jackson, Eds., 2009.

A. Wimo and M. Prince, “World Alzheimer report 2010: the global economic impact of dementia,” 2010, http://www.alz.co.uk/research/world-report

J. Giustiniani, M. Sineus, E. Sardin et al., “Decrease of the immunophilin FKBP52 accumulation in human brains of Alzheimer's disease and FTDP-17,” Journal of Alzheimer's Disease, vol. 29, no. 2, pp. 471–483, 2012.

P. E. Cramer, J. R. Cirrito, D. W. Wesson et al., “ApoE-directed therapeutics rapidly clear β-amyloid and reverse deficits in AD mouse models,” Science, vol. 335, no. 6075, pp. 1503–1506, 2012.

J. Gräff, D. Rei, J. S. Guan et al., “An epigenetic blockade of cognitive functions in the neurodegenerating brain,” Nature, vol. 483, no. 7388, pp. 222–226, 2012.

J. S. Guan, S. J. Haggarty, E. Giacometti et al., “HDAC2 negatively regulates memory formation and synaptic plasticity,” Nature, vol. 459, no. 7243, pp. 55–60, 2009.

M. Kilgore, C. A. Miller, D. M. Fass et al., “Inhibitors of class 1 histone deacetylases reverse contextual memory deficits in a mouse model of Alzheimer's disease,” Neuropsychopharmacology, vol. 35, no. 4, pp. 870–880, 2010.

J. Lamb, E. D. Crawford, D. Peck et al., “The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease,” Science, vol. 313, no. 5795, pp. 1929–1935, 2006.

S. Gupta, S. Y. Kim, S. Artis et al., “Histone methylation regulates memory formation,” Journal of Neuroscience, vol. 30, no. 10, pp. 3589–3599, 2010.

L. Welberg, “Learning and memory: HDAC2 is the one,” Nature Reviews Neuroscience, vol. 10, no. 7, p. 471, 2009.

S. H. Chen, H. M. Wu, B. Ossola et al., “Suberoylanilide hydroxamic acid, a histone deacetylase inhibitor, protects dopaminergic neurons from neurotoxin-induced damage,” The British Journal of Pharmacology, vol. 165, no. 2, pp. 494–505, 2012.

Y. M. Sung, T. Lee, H. Yoon et al., “Mercaptoacetamide-based class II HDAC inhibitor lowers Aβ levels and improves learning and memory in a mouse model of Alzheimer's disease,” Experimental Neurology, vol. 239, no. 1, pp. 192–201, 2013.

V. M. Richon, “Cancer biology: mechanism of antitumour action of vorinostat (suberoylanilide hydroxamic acid), a novel histone deacetylase inhibitor,” The British Journal of Cancer, vol. 95, no. 1, pp. S2–S6, 2006.

D. Huang, X. Li, Y. Wei, and Z. Xiu, “A novel series of l-2-benzyloxycarbonylamino-8-(2-pyridyl)-disulfidyloctanoic acid derivatives as histone deacetylase inhibitors: design, synthesis and molecular modeling study,” European Journal of Medicinal Chemistry, vol. 52, pp. 111–122, 2012.

S. Peleg, F. Sananbenesi, A. Zovoilis et al., “Altered histone acetylation is associated with age-dependent memory impairment in mice,” Science, vol. 328, no. 5979, pp. 753–756, 2010.

C. G. Vecsey, J. D. Hawk, K. M. Lattal et al., “Histone deacetylase inhibitors enhance memory and synaptic plasticity via CREB: CBP-dependent transcriptional activation,” The Journal of Neuroscience, vol. 27, no. 23, pp. 6128–6140, 2007.

F. Thiele, S. Young, R. Buchert, and F. Wenzel, “Voxel-based classification of FDG PET in dementia using inter-scanner normalization,” NeuroImage, vol. 77, pp. 62–69, 2013.

N. M. Tsankova, O. Berton, W. Renthal, A. Kumar, R. L. Neve, and E. J. Nestler, “Sustained hippocampal chromatin regulation in a mouse model of depression and antidepressant action,” Nature Neuroscience, vol. 9, no. 4, pp. 519–525, 2006.

W. Renthal, I. Maze, V. Krishnan et al., “Histone deacetylase 5 epigenetically controls behavioral adaptations to chronic emotional stimuli,” Neuron, vol. 56, no. 3, pp. 517–529, 2007.

T. Nuutinen, T. Suuronen, A. Kauppinen, and A. Salminen, “Valproic acid stimulates clusterin expression in human astrocytes: implications for Alzheimer's disease,” Neuroscience Letters, vol. 475, no. 2, pp. 64–68, 2010.

M. Mielcarek, C. L. Benn, S. A. Franklin et al., “SAHA decreases HDAC 2 and 4 levels in vivo and improves molecular phenotypes in the R6/2 mouse model of Huntington's disease,” PLoS ONE, vol. 6, no. 11, Article ID e27746, 2011.

A. Fischer, F. Sananbenesi, X. Wang, M. Dobbin, and L.-H. Tsai, “Recovery of learning and memory is associated with chromatin remodelling,” Nature, vol. 447, no. 7141, pp. 178–182, 2007.

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.

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, p. 844, 2013.

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.

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.

S. C. Basak, “Recent developments and future directions at current computer aided drug design,” Current Computer-Aided Drug Design, vol. 9, no. 1, p. 1, 2013.

C. Y. 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.

W.-L. Liao and F.-J. Tsai, “Personalized medicine: a paradigm shift in healthcare,” BioMedicine, vol. 3, no. 2, pp. 66–72, 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.

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

C. Y. 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.

H. J. Huang, K. J. Lee, H. W. Yu et al., “Structure-based and ligand-based drug design for HER 2 receptor,” Journal of Biomolecular Structure and Dynamics, vol. 28, no. 1, pp. 23–37, 2010.

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

K. C. Chen and C. Y.-C. 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.

M. F. Sun, T. T. Chang, K. C. Chen et al., “Treat Alzheimer's disease by traditional Chinese medicine?” Molecular Simulation, vol. 37, no. 11, pp. 923–931, 2011.

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 & Dynamics, vol. 32, no. 5, pp. 776–791, 2014.

H.-C. Tang and C. Y.-C. Chen, “Investigation of the novel lead of melanocortin 1 receptor for pigmentary disorders,” Evidence-Based Complementary and Alternative Medicine, vol. 2014, Article ID 254678, 13 pages, 2014.

S.-S. Chang, H.-J. Huang, and C. Y. 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.

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.

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. Y. Chen, Y. H. Chang, D. T. Bau et al., “Ligand-based dual target drug design for H1N1: swine flu—a preliminary first study,” Journal of Biomolecular Structure and Dynamics, vol. 27, no. 2, pp. 171–178, 2009.

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.

T. Y. Tsai, K. W. Chang, and C. Y. 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.

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.

J. C. Bressi, A. J. Jennings, R. Skene et al., “Exploration of the HDAC2 foot pocket: synthesis and SAR of substituted N-(2-aminophenyl)benzamides,” Bioorganic and Medicinal Chemistry Letters, vol. 20, no. 10, pp. 3142–3145, 2010.

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 and Modelling, vol. 21, no. 4, pp. 289–307, 2003.

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.

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

J. C. Bressi, R. de Jong, Y. Wu et al., “Benzimidazole and imidazole inhibitors of histone deacetylases: synthesis and biological activity,” Bioorganic and Medicinal Chemistry Letters, vol. 20, no. 10, pp. 3138–3141, 2010.

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.

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.

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.

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.

T. A. Darden and L. G. Pedersen, “Molecular modeling: an experimental tool,” Environmental Health Perspectives, vol. 101, no. 5, pp. 410–412, 1993.

R. A. Laskowski and M. B. Swindells, “LigPlot+: multiple ligand-protein interaction diagrams for drug discovery,” Journal of Chemical Information and Modeling, vol. 51, no. 10, pp. 2778–2786, 2011.

E. Chovancova, A. Pavelka, P. Benes et al., “CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures,” PLoS Computational Biology, vol. 8, no. 10, Article ID e1002708, 2012.