BioMed Research International

Volume 2014 (2014), Article ID 809816, 13 pages

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

Research Article

Lead Screening for CXCR4 of the Human HIV Infection Receptor Inhibited by Traditional Chinese Medicine

Tzu-Chieh Hung,¹ Wen-Yuan Lee,1,2,3 Kuen-Bao Chen,1,2,4 and Calvin Yu-Chian Chen1,2,5

¹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 City 40447, Taiwan

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

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

Received 22 February 2014; Revised 5 March 2014; Accepted 5 March 2014; Published 5 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.

The acquired immunodeficiency syndrome (AIDS) is a serious worldwide disease caused by the human immunodeficiency virus (HIV) infection. Recent research has pointed out that the G protein-coupled chemokine receptor CXCR4 and the coreceptor C-C chemokine receptor type 5 (CCR5) are important targets for HIV infection. The traditional Chinese medicine (TCM) database has been screened for candidate compounds by simulating molecular docking and molecular dynamics against HIV. Saussureamine C, 5-hydroxy-L-tryptophan, and diiodotyrosine are selected based on the highest docking score. The molecular dynamics is helpful in the analysis and detection of protein-ligand interactions. According to the analysis of docking poses, hydrophobic interactions, hydrogen bond variations, and the comparison of the effect on CXCR4 and CCR5, these results indicate Saussureamine C may have better effect on these two receptors. But for some considerations, diiodotyrosine could make the largest variation and may have some efficacy contrary to expectations.

1. Introduction

Recently, an important topic of the acquired immunodeficiency syndrome (AIDS) had been published in 2013. The G protein-coupled chemokine receptor CXCR4 and the coreceptor C-C chemokine receptor type 5 (CCR5) are important targets for HIV infection [1].

The human immunodeficiency virus (HIV) is a retrovirus which causes the AIDS [2–5]. During the course of this disease, the human immune system becomes compromised, and considerably weaker. The HIV virus is transmitted through a population rapidly by unprotected sexual intercourse [6, 7], contaminated medical equipment [8, 9], vertical infection [10, 11], and bodily fluids.

There were 35.3 million people living with HIV in 2012 and 2.1 million of these patients were adolescents (recorded by WHO). There are still no defined vaccines or drugs available to kill all HIV viruses in patients and then the highly active antiretroviral therapy (HAART) is the standard of care for patients with advanced infection in current treatment [12]. HARRT uses a complex of transcription inhibitors to slow down transcription and then decrease the patient’s total burden of HIV, but this treatment is too expensive and medical costs become a social liability.

Chemokine receptors are critical regulators of cell migration in the context of immune surveillance, inflammation, and development. The one of 19 known human chemokine receptors, the G protein-coupled chemokine receptor CXCR4, is specifically implicated in cancer metastasis and HIV-1 infection [13]. The CCR5 is a receptor for the T-cells that play a central role in cell-mediated immunity against viruses and pathogens. CXCR4 and CCR5 have been defined as coreceptors for the HIV antigen gp120 and then HIV can infect the cell by targeting these receptors [14, 15]. Thus, preventing HIV from targeting the receptor could prevent the virus infection [16, 17]. Based on these observations, the drug IT1t is a CCR5 receptor antagonist, thereby blocking the HIV protein from associating with the receptor.

Computer-aided drug design (CADD) is an in silico simulation technique to screen for novel compounds by their structure and bioactivity from database. The difference from traditional drug design is that CADD has the advantages of both greater speed and lower cost for drug development. The structure-based drug design and ligand-based drug design are two major application areas of CADD. We used CADD to investigate based on structure-based drug design and molecular dynamics [18–21].

Recently, more attention has been given to personalized medicine and biomedicine [22, 23]. By this knowledge, people could discover the association from the mutation [24, 25], pathway [26, 27], the cause for special disease [28–30], and even the case from clinical diagnosis [31] with disease. Traditional Chinese medicine (TCM) is an identified personalized medicine and this clinical diagnosis has an important role in Asia, especially in China, Taiwan, Korea, and Japan. In 2011, the TCM Database@Taiwan (http://tcm.cmu.edu.tw/) [32] which is the largest traditional Chinese medicine database in the world was established. In this TCM database, both of 2D and 3D chemical structures, bioactivity, and molecular information for over 61,000 compounds of traditional Chinese medicinal herbs could be generated. Until today, there has been successful novel drug discovery from the TCM Database@Taiwan, such as cancer treatment [33–36], stroke prevention [37], EGFR inhibition [38], inflammation inhibition [39], pain relief [19], and antivirals [40–43]. Since the application system of the website [44] and the cloud computing platform [45], the TCM Database@Taiwan could be rigorous and valuable for TCM application and drug design.

In this research, we screen a possible compound against HIV from the TCM Database@Taiwan based on molecular docking. Finally, we use molecular dynamics (MD) simulation to investigate the protein-ligand interactions that may contribute to evaluate the effect of human HIV receptor inhibition.

2. Materials and Methods

2.1. Data Set

A total of 61,000 TCM compounds were downloaded from the TCM database (http://tcm.cmu.edu.tw/). The CXCR4 (PDB ID: 3ODU) crystal structure was generated from RCSB Protein Data Bank (PDB). Based on the literature, IT1t was defined as a control [13]. Then, the Accelrys Discovery Studio 2.5 (DS 2.5) was used for the molecular docking selection.

2.2. Disorder Protein Detection

Because the disorder protein plays an important role in drug design, we take the crystal structure to predict the disorder region by the Database of Protein Disorder (DisProt: http://www.disprot.org/). Based on the prediction, we can decide the character of the docking site and assess the efficacy of the drug [46, 47].

After a comparison of the disorder regions and the defined docking sites, we could evaluate drug efficacy from the protein-ligand interaction.

2.3. Molecular Docking

The docking platform dock IT1t and TCM compounds to CXCR4 in the force field of CHARMm [48] by using LigandFit [49], which is a receptor-rigid docking algorithm in Discovery Studio 2.5 (DS 2.5). The docking site of CXCR4 was defined by the research around Trp94, Asp97, Trp102, Val112, His113, Tyr116, Arg183, Cys186, Asp187, Arg188, and Glu288 [1, 13]. The complexes of candidate compounds with CXCR4 were selected for hydrophobic interactions by Ligplus [50, 51].

2.4. Molecular Dynamics Simulation

These selected ligands must take preparation by using SwissParam (http://swissparam.ch/) [52] before MD simulation based on the reference force field [53] of GROMACS 4.5.5 [54]. The CXCR4 combines with ligands as the complex set into the full buffer (or solution) simulation box. This cubic box was set with a minimum distance of 1.2 Å from the complex and in this simulation box was solvated with the TIP3P water model to regulate the sodium and chloride ion to neutralize complex charges. The minimization applies the steepest descent method for 5,000 steps in the beginning. Then the last structure was transferred to MD simulation. The equilibration was based on the Berendsen weak thermal coupling method under the 100 ps constant temperature (PER ensemble). The Particle-Mesh Ewald (PME) [55] method was used to calculate the electrostatic interactions in 2 fs per time and the numbers of steps were 5,000,000 times and then accomplished 10 ns simulation time of MD. Gromacs 4.5.5 also has protocol to analyze MD trajectories, RMSD, and energy variations.

3. Results and Discussion

3.1. The Detection of Disorder Protein

The disorder protein is defined as unfolding protein. For this character, while the drug is docking to the disorder region, the complex will stabilize with difficulty. There are some references [46, 47] that indicate that the disorder region is not any defined domain; therefore, therefore the drug targeting the disorder region may have lower side effects than a drug interacting with the widespread domain. Thus, the disorder region can be defined as a hard work for drug design. The disorder regions of CXCR4 are defined as having a disposition of over than 0.5 (Figure 1). This result indicates that the important amino acids are not disorder regions; thus, the complex selected based on docking could have an influence on CXCR4.

Figure 1: The disorder and binding site detection. The green curve in the figure is the disorder disposition of each amino acid and the red lines are the residues of the important amino acids.

3.2. Molecular Docking

After molecular docking and ranking by docking score, the top three TCM compounds are defined as candidate compounds which are Saussureamine C, 5-hydroxy-L-tryptophan, and diiodotyrosine derived from the TCM herbs Saussurea lappa Clarke, Mucuna pruriens (L) D., and Ox Thyroid of Bos taurus domesticus Gmelin (or Bubalus bubalis L.), respectively (Table 1). The top ranking compound, Saussureamine C is know for anti-ulcer treatment [56] and the herb Saussurea lappa Clarke can inhibit breast cancer migration [57], treat cardiovascular disease [58, 59], be used against hepatotoxic activity [60], and inhibit cytotoxic T lymphocytes [61]. The second ranked herb, Mucuna pruriens, can prevent Parkinson’s disease based on antioxidation [62, 63]. The third ranked compound, diiodotyrosine from the herb Ox Thyroid of Bos taurus domesticus Gmelin, has been related to the thyroid [64–66] and pH-sensitive pore-forming [67]. As reported in the literature, most of these compounds can have an effect on immunity, especially on cancer. For the above reference, we suggest that these compounds can have an effect on T-cell receptors, such as CXCR4. For our previous research about the coreceptor CCR5, the top and second compounds are the same. We make the sequence align between CXCR4 and CCR5 (Figure 2). This result indicates that the most important amino acids both of CXCR4 and CCR5 are similar and presents that the docking site and binding domain designed by important amino acids are similar. For this result, the drug has an influence on CXCR4 which may also have an effect on CCR5.

Table 1: Scoring functions of the top three compounds and the inhibitors of CXCR4.

Compounds
Herbs
-PLP1
-PLP2
Dock score

Saussureamine C
Saussurea lappa Clarke
58.25
49.34
207.651

5-Hydroxy-L-tryptophan
Mucuna pruriens
21.5
20.69
205.102

Diiodotyrosine
Ox thyroid of Bos  taurus  domesticus  Gmelin (or Bubalus  bubalis  L.)
34.37
30.49
203.933

IT1t*

73.42
72.93
71.873

Control.

Figure 2: The sequence analysis between two HIV influence receptors. The deep blue of sequence means these two regions are the same, low blue means similar, and white means different.

The structures of the candidate compounds and control were screened from TCM database (Figure 3). Then, the docking poses, the docking site, and the amino acid neighbors by ligands are presented (Figure 4). From this result, we observe that Asp97 and Asp187 are defined as the amino acids that can interact with all the selected ligands; thus, these amino acids may play important roles in target function of CXCR4.

(a)

(b)

(c)

(d)

Figure 3: The structure of control and candidate TCM compounds. (a) IT1t, (b) Saussureamine C, (c) 5-hydroxy-L-tryptophan, and (d) diiodotyrosine.

(a)

(b)

(c)

(d)

(e)

Figure 4: The docking poses of ligands. (a) The crystal structure of CXCR4 and the docking site, (b) IT1t, (c) Saussureamine C, (d) 5-hydroxy-L-tryptophan, and (e) diiodotyrosine.

The hydrophobic interaction can be analyzed by Ligplus (Figure 5). This result shows that the signed deep red amino acids are at a high frequency while ligands target in docking site. Most of these amino acids have been defined as important amino acids in the literature; thus, this hydrophobic interaction analysis is credible to present the selected compounds effect on CXCR4.

(a)

(b)

(c)

(d)

Figure 5: Ligplot illustrates the protein-ligand interactions. (a) IT1t, (b) Saussureamine C, (c) 5-hydroxy-L-tryptophan, and (d) diiodotyrosine. The deep red color of the hydrophobic interactions indicates a high frequency in all ligand interactions.

3.3. Molecular Dynamics Simulation

The RMSD and total energy of a complex during MD simulation were calculated (Figure 6). The total energy is in the range −2308~−2294*10³ kJ/mol and tends to −2300*10³ kJ/mol. Although the compounds Saussureamine C and 5-hydroxy-L-tryptophan have a high variation in ligand RMSD, the amplitude in both ligand RMSD and protein RMSD is more gentle after MD 8 ns. Thus, we suggest that these two ligands may make the complex balance quickly.

Figure 6: Measures of the MD trajectories. (a) Complex RMSD, (b) ligand RMSD, and (c) the total energy.

The clustering based on RMSD variation could be calculated (Figure 7). In this result, the complex with compound Saussureamine C or 5-hydroxy-L-tryptophan and the protein structure position will be similar to the same group after MD 5 ns. This result confirms our suggestion and presents that these two ligands will be fit to CXCR4.

(a)

(b)

(c)

(d)

Figure 7: The clustering based on RMSD variation while protein interaction. (a) IT1t, (b) Saussureamine C, (c) 5-Hydroxy-L-tryptophan, (d) Diiodotyrosine.

The calculations of RMSD in each residue during the whole MD, the root mean square fluctuations (RMSF), show the variation in CXCR4 (Figure 8). In this result, the pick sites of residue in four complexes are similar, and then the pick in residues 50 to 300 could help to define the relation in interaction. The similar pick sites present these compounds effect on the same amino acids and these amino acids may play important function in CXCR4. The largest amplitude of the complex with diiodotyrosine may indicate that this compound will have stronger effect on protein making the complex unstable.

Figure 8: The RMSF of each residue of protein.

The H-bond occupancy and structure variation were calculated for the protein-ligand interactions (Figures 9–12). In the complex with control, the H-bond occupancy is less than 10% (Figure 9(a)), but there was still variation in both the position and the composition (Figures 9(b)-8). These variations, we suggest, are due to the complex having more hydrophobic interactions than others might have. This is an important function, while protein-ligand interactions then inhibit the influence of the G protein.

(a)

(b)

Figure 9: The variation of control and CXCR4 complex in MD simulation. (a) H-bond variation, (b) structure variation. The (1)–(9) red color indicates the difference through MD.

(a)

(b)

Figure 10: The variation of Saussureamine C and CXCR4 complex in MD simulation. (a) H-bond variation, (b) structure variation. The (1)–(7) red color indicates the difference through MD.

(a)

(b)

Figure 11: The variation of 5-hydroxy-L-tryptophan and CXCR4 complex in MD simulation. (a) H-bond variation, (b) structure variation. The (1)–(5) red color indicates the difference through MD.

(a)

(b)

Figure 12: The variation of diiodotyrosine and CXCR4 complex in MD simulation. (a) H-bond variation, (b) structure variation. The (1)–(10) red color indicates the difference through MD.

In Figure 10(a), the different atoms of Glu288 interact with ligand after 1000 ps. Besides composition variation, the structure variations of CXCR4 with Saussureamine C are similar to the control (Figure 10(b)). Although the hydrophobic interaction are fewer than the control, the binding site Glu288, for chemokine (defined by Uniprot), will be targeted by Saussureamine C to replace the force of the hydrophobic interaction.

The high H-bond occupancy for CXCR4 with 5-hydroxy-L-tryptophan is not only Glu288 (Figure 11(a)). The Asp97 and His113 had been defined as chemokine binding site by the UniProt; the ligand interacts with these functional sites may cause the structure to loose the helix, which will effect the function of the G protein (Figure 11(b)).

The larger variation in both H-bond and structure for the Diiodotyrosine complex might indicate that this complex does not tend to balance (Figure 12). This situation might indicate that diiodotyrosine is a compound with stronger force and long term interaction for CXCR4. For this possible suggestion, diiodotyrosine might not be bad for the inhibition of CXCR4 but using diiodotyrosine it may have some consideration or supporting measures.

The difference from CCR5 is that the structure variation in CXCR4 is more intense, while the protein-ligand interaction and the ligand are the same. Accordingly, we suggest the ligand inhibits these two receptors which might be intense interaction to change CXCR4 structure and strong target to prevent other influence on CCR5.

The pathway definition is according to the calculation of caver 3.0 [68] to determine the interpath in protein while interacting (Figure 13). In this result, we could find most pathways around the docking site, and only diiodotyrosine is different. This phenomenon may be caused by larger structural variations making some pathways hidden to the protein as the virus enters the human cell. Although we could not identify the existence of this hypothesized pathway, we note that Diiodotyrosine is different from the other selected compounds.

(a)

(b)

(c)

(d)

Figure 13: The pathway calculation of CXCR4 complex in MD simulation. (a) IT1t, (b) Saussureamine C, (c) 5-hydroxy-L-tryptophan, and (d) diiodotyrosine.

4. Conclusion

Based on the above discussion, we have seen TCM compounds Saussureamine C, 5-hydroxy-L-tryptophan, and diiodotyrosine can have an effect on CXCR4 against HIV infection. The control interacts with CXCR4 by more hydrophobic interactions but with other compounds on the basis of H-bond, a stronger force in the interaction. The structural variations in CXCR4, with the same compounds, being more intense than in CCR5 might be due to the different modes of protein inhibition. Although Saussureamine C is the best selection from drug design, diiodotyrosine might have some efficacy contrary to expectations from the largest variation.

Conflict of Interests

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

Authors’ Contribution

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

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, 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) and 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. We also thank Pei-Chen Chang for help in the analysis of these data. Finally, the authors’ gratitude goes to Dr. Tim Williams, Asia University.

References

Q. Tan, Y. Zhu, J. Li, et al., “Structure of the CCR5 chemokine receptor-HIV entry inhibitor maraviroc complex,” Science, vol. 341, no. 6152, pp. 1387–1390, 2013.

F. B. Sinoussi, J. C. Chermann, and F. Rey, “Isolation of a T-lymphotropic retrovirus from a patient at risk for acquired immune deficiency syndrome (AIDS),” Science, vol. 220, no. 4599, pp. 868–871, 1983.

M. W. Cloyd, J. J.-Y. Chen, and L. Wang, “How does HIV cause AIDS? The homing theory,” Molecular Medicine Today, vol. 6, no. 3, pp. 108–111, 2000.

R. A. Weiss, “How does HIV cause AIDS?” Science, vol. 260, no. 5112, pp. 1273–1279, 1993.

D. C. Douek, M. Roederer, and R. A. Koup, “Emerging concepts in the immunopathogenesis of AIDS,” Annual Review of Medicine, vol. 60, pp. 471–484, 2009.

S. Weller and K. Davis, “Condom effectiveness in reducing heterosexual HIV transmission,” The Cochrane Database of Systematic Reviews, no. 1, 2002.

B. Leynaert, A. M. Downs, and I. De Vincenzi, “Heterosexual transmission of human immunodeficiency virus: variability of infectivity throughout the course of infection,” American Journal of Epidemiology, vol. 148, no. 1, pp. 88–96, 1998.

D. M. Bell, “Occupational risk of human immunodeficiency virus infection in healthcare workers: an overview,” The American Journal of Medicine, vol. 102, no. 5 B, pp. 9–15, 1997.

J. L. Cleveland, L. Barker, B. F. Gooch, E. M. Beltrami, and D. Cardo, “Use of HIV postexposure prophylaxis by dental health care personnel: an overview and updated recommendations,” Journal of the American Dental Association, vol. 133, no. 12, pp. 1619–1626, 2002.

J. C. Forbes, A. M. Alimenti, J. Singer et al., “A national review of vertical HIV transmission,” AIDS, vol. 26, no. 6, pp. 757–763, 2012.

K. M. Little, P. H. Kilmarx, A. W. Taylor, C. E. Rose, E. D. Rivadeneira, and S. R. Nesheim, “A review of evidence for transmission of HIV from children to breastfeeding women and implications for prevention,” Pediatric Infectious Disease Journal, vol. 31, no. 9, pp. 938–942, 2012.

D. Finzi, M. Hermankova, T. Pierson et al., “Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy,” Science, vol. 278, no. 5341, pp. 1295–1300, 1997.

B. Wu, E. Y. T. Chien, C. D. Mol et al., “Structures of the CXCR4 chemokine GPCR with small-molecule and cyclic peptide antagonists,” Science, vol. 330, no. 6007, pp. 1066–1071, 2010.

E. A. Berger, P. M. Murphy, and J. M. Farber, “Chemokine receptors as HIV-1 coreceptors: roles in viral entry, tropism, and disease,” Annual Review of Immunology, vol. 17, pp. 657–700, 1999.

E. Coakley, C. J. Petropoulos, and J. M. Whitcomb, “Assessing chemokine co-receptor usage in HIV,” Current Opinion in Infectious Diseases, vol. 18, no. 1, pp. 9–15, 2005.

V. M. Muniz-Medina, S. Jones, J. M. Maglich et al., “The relative activity of “function sparing” HIV-1 entry inhibitors on viral entry and CCR5 internalization: is allosteric functional selectivity a valuable therapeutic property?” Molecular Pharmacology, vol. 75, no. 3, pp. 490–501, 2009.

W.-T. Choi, R. Nedellec, M. Coetzer et al., “CCR5 mutations distinguish N-terminal modifications of RANTES (CCL5) with agonist versus antagonist activity,” Journal of Virology, vol. 86, no. 18, pp. 10218–10220, 2012.

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.-C. Chen, “Drug design for neuropathic pain regulation from traditional Chinese medicine,” Scientific Reports, vol. 3, p. 844, 2013.

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

W.-L. Liao and F.-J. Tsai, “Personalized medicine: a paradigm shift in healthcare,” BioMedicine, vol. 3, no. 2, pp. 66–72, 2013.

F.-J. Tsai, “Biomedicine brings the future nearer,” BioMedicine, vol. 1, no. 1, p. 1, 2011.

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.

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.-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–81, 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.

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

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.

Y.-T. Chang, W.-D. Lin, Z.-N. Chin et al., “Nonketotic hyperglycinemia: a case report and brief review,” BioMedicine, vol. 2, no. 2, pp. 80–82, 2012.

C.-H. Wang, W.-D. Lin, and F.-J. Tsai, “Craniofacial dysmorphism, what is your diagnosis?” BioMedicine, vol. 2, no. 2, pp. 49–50, 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.

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.

W. I. Tou and C. Y.-C. Chen, “Traditional Chinese medicine as dual guardians against hypertension and cancer?” Journal of Biomolecular Structure and Dynamics, vol. 30, no. 3, pp. 299–317, 2012.

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.

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

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.

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

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.

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.

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

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

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

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.

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.

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.

M. Yoshikawa, S. Hatakeyama, Y. Inoue, and J. Yamahara, “Saussureamines A, B, C, D, and E, new anti-ulcer principles from Chinese Saussureae radix,” Chemical and Pharmaceutical Bulletin, vol. 41, no. 1, pp. 214–216, 1993.

Y. K. Choi, S.-G. Cho, and S.-M. Woo, “Saussurea lappa clarke-Derived costunolide prevents TNF alpha induced breast cancer cell migration and invasion by inhibiting NF-κB activity,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 936257, 10 pages, 2013.

O. P. Upadhyay, J. K. Ojha, H. S. Bajpai, and A. K. Hathwal, “Study of kustha (saussurea lappa, clarke) in ischaemic heart disease,” Ancient Science of Life, vol. 13, no. 1-2, pp. 11–18, 1993.

T. S. M. Saleem, N. Lokanath, A. Prasanthi, M. Madhavi, G. Mallika, and M. N. Vishnu, “Aqueous extract of Saussurea lappa root ameliorate oxidative myocardial injury induced by isoproterenol in rats,” Journal of Advanced Pharmaceutical Technology and Research, vol. 4, no. 2, pp. 94–100, 2013.

S. Yaeesh, Q. Jamal, A. J. Shah, and A. H. Gilani, “Antihepatotoxic activity of Saussurea lappa extract on D-galactosamine and lipopolysaccharide-induced hepatitis in mice,” Phytotherapy Research, vol. 24, no. 2, pp. S229–S232, 2010.

M. Taniguchi, T. Kattaoka, H. Suzuki et al., “Costunolide and dehydrocostus lactone as inhibitors of killing function of cytotoxic T lymphocytes,” Bioscience, Biotechnology and Biochemistry, vol. 59, no. 11, pp. 2064–2067, 1995.

S. K. Yadav, J. Prakash, S. Chouhan, and S. P. Singh, “Mucuna pruriens seed extract reduces oxidative stress in nigrostriatal tissue and improves neurobehavioral activity in paraquat-induced Parkinsonian mouse model,” Neurochemistry International, vol. 62, no. 8, pp. 1039–1047, 2013.

S. Ovallath and P. Deepa, “The history of parkinsonism: descriptions in ancient Indian medical literature,” Movement Disorders, vol. 28, no. 5, pp. 566–568, 2013.

A. Phatarphekar, J. M. Buss, and S. E. Rokita, “Iodotyrosine deiodinase: a unique flavoprotein present in organisms of diverse phyla,” Molecular BioSystems, vol. 10, no. 1, pp. 86–92, 2014.

K. B. Paul, J. M. Hedge, C. Macherla et al., “Cross-species analysis of thyroperoxidase inhibition by xenobiotics demonstrates conservation of response between pig and rat,” Toxicology, vol. 312, no. 1, pp. 97–107, 2013.

G. Toth and B. Noszal, “Thyroid hormones and their precursors I. Biochemical properties,” Acta Pharmaceutica Hungarica, vol. 83, no. 2, pp. 35–45, 2013.

D. H. Haas and R. M. Murphy, “Design of a pH-sensitive pore-forming peptide with improved performance,” Journal of Peptide Research, vol. 63, no. 1, pp. 9–16, 2004.

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.