1. Introduction
Alzheimer’s disease (AD) is the most common form of neurodegenerative disease [1–3] with symptoms ranging from intellectual deterioration, cognitive impairment [4, 5] to abnormal behavior, personality changes, depression, and sleep disorders [6, 7]. Neuronal loss, senile plaques in the cerebral cortex, and neurofibrillary tangles in the temporal and frontal lobes are prominent in the brains of AD patients. Deterioration in the limbic system, including pathological changes in the hippocampus and atrophy of pyramidal cells and the amygdale, has also been reported [8, 9]. Formation of
Pioneer research has discovered that the concentration of FKBP52, a FK506-binding protein, is of high density in brain [15] and 40 times higher within the central nervous system than in the immune system [16]. FKBP52 is not only related to immune functions but is also important for CNS protective properties. Interestingly, they can easily bind with highly phosphorylated tau proteins [16], thereby reducing the accumulation of tau proteins [16–18]. FKBP52 was initially discovered as a cochaperone of steroid receptor heterocomplexes [19, 20] and is a member of the FK506-binding protein (FKBP) of immunophilins. FKBP52, which possesses a chaperone function, has a PPIase domain also called FK506-binding domain (FKBD), composed of the first 138 amino acids from the N-terminal, exhibits peptidylprolyl isomerase (PPIase) activity, and plays an important role in regulating tau proteins.
Tacrolimus (FK506) is an immunosuppressant typically used to reduce graft rejection during organ transplants [21]. Research shows that FK506 can also reduce pathological changes in tau proteins [22] and promote neuroregeneration [23, 24]. In addition, FK506 can disrupt steroid receptor heterocomplexes, leading to the release of FKBP52 from the complex and enabling it to function in the nervous system [24–26].
In recent years, most beta-amyloid-centric drug developments have been unsuccessful during Phase III of clinical trials [27, 28]. Currently, the personalized medicine and biomedicine defined an important knowledge for diagnosis and treatment analysis [29, 30], such as rare diseases [31, 32], regional disease [33], signal pathway [34–36], and gene association [37, 38]. The traditional Chinese medicine is an important medicine culture in Asia defined as a personalized medicine, system biology, and biomedicine in medicine practices widely applied as treatments for cancer [39, 40], cardiovascular disease [41, 42], diabetes [43], virus disease [44], inflammation [45], and metabolic disorders [46], indicating its large therapeutic potential for various diseases [47]. In this moment, the cloud computation and system biology could help the screening for TCM application [48, 49]. We have already designed many therapies to different targets [48, 50–57] and have successfully developed method for designing drug in disordered protein [58, 59]. The aim of this research is to screen for potential non-amyloid-beta-centric leads from traditional Chinese medicine targeting FKBP52 and to investigate its mechanisms with the hope of providing important insights for forwarding Alzheimer’s disease.
2. Results
2.1. SVM/MLR/Bayesian Network
The ten representative descriptors determined by GFA were ES_Count_aaCH, ES_Count_sOH, ES_Sum_aaCH, ES_Sum_sssN, Num_Rings6, Molecular_PolarSurfaceArea, CHI_2, Jurs_TPSA, Minimized_Energy, and Shadow_Ylength, suggesting certain relationships between electrotopological properties and polar surface area of ligands with bioactivity (Table
The GFA model with the highest correlation between experimental and predicted pIC50 values
Table 1
Docking scores, bioactivity predictions, and consensus votinga of the TCM candidates and control.
Name | SVM | MLR | Bayesian Network | Dock score | Sum of scores |
---|---|---|---|---|---|
Lathyranoic acid A | 4.81 (0) | 5.12 (0) | 5.45 (0) | 86.20 (1) | 1 |
Daphnetoxin | 5.09 (0) | 8.89 (1) | 7.46 (2) | 84.88 (1) | 4 |
Aurantiamide | 4.89 (0) | 6.42 (0) | 5.43 (0) | 83.89 (1) | 1 |
(6aR,11aR)-9,10-Dimethoxypterocarpan-3-O-beta-D-glucoside | 5.00 (0) | 8.33 (1) | 5.56 (0) | 83.50 (0) | 1 |
Picrasidine M | 5.10 (1) | 7.51 (0) | 5.65 (0) | 69.60 (0) | 1 |
12-O-Acetylphorbol-13-tigliate | 4.91 (0) | 7.40 (0) | 6.15 (0) | 68.88 (0) | 0 |
20-O-(2′E,4′E-Decadienoyl)ingenol | 4.91 (0) | 6.87 (0) | 6.35 (2) | 68.52 (0) | 2 |
Howiinol A II | 5.24 (1) | 6.87 (0) | 5.40 (0) | 68.11 (0) | 1 |
Moellendorffiline | 5.15 (1) | 6.64 (0) | 5.98 (0) | 63.24 (0) | 1 |
Lythrancine II | 5.03 (0) | 8.92 (1) | 6.74 (2) | 62.03 (0) | 3 |
|
|||||
Tacrolimus* | 5.11 | 6.25 | 5.33 | 58.10 | — |
aVoting scores are given in parenthesis.
[figures omitted; refer to PDF]
2.2. Docking Analysis
Summary of docking interactions between FKBP52 and the selected ligands is listed in Table 2. Tacrolimus interacted with FKBP52 via a single H-bond at Tyr113 and three hydrophobic interactions. Higher amounts of total interactions were recorded for TCM candidates. With regard to individual residues, interaction percentages suggest the importance of Arg73, Phe77, and Tyr113 for ligand binding. Interactions with Arg73 were ligand-dependent, whereas Phe77 and Tyr113 were mutual interaction sites for hydrophobic interactions and H-bonds, respectively. Cross comparison reveals that TCM candidates form interactions with all Tacrolimus-interacting residues except Glu85. In addition, each TCM candidate is anchored by interactions with residues not recorded for Tacrolimus. The higher number of interactions observed for TCM candidates might be associated with higher dock scores and suggests higher binding affinity to FKBP52 than to Tacrolimus. Figure 2 was demonstrated by Ligplot, showing the hydrophobic interactions and hydrogen bonds between control compound and TCM candidates. Except for the H-bond between residue Tyr57 and Daphnetoxin, the results were the same for docking analysis.
Table 2
Protein-ligand interactions recorded during docking.
Tyr57 | Asp68 | Arg73 | Phe77 | Glu85 | Val86 | Trp90 | Tyr113 | Lys121 | Phe130 | Totala | |
---|---|---|---|---|---|---|---|---|---|---|---|
Tacrolimus | — | — | Hb | Hb | Hb | — | — | H | — | — | 4 |
Daphnetoxin | — | Hb | Pi | Hb | — | — | — | H | Hb | — | 5 |
20-O-(2′E,4′E-Decadienoyl)ingenol | — | Hb | Hb | Hb | — | Hb | Hb | H | H | — | 7 |
Lythrancine II | H | H | H | Pi/Hb | — | — | — | H | — | Hb | 7 |
|
|||||||||||
Interaction percentageb | 25% | 75% | 100% | 100% | 25% | 25% | 25% | 100% | 50% | 25% |
bPercentage of ligands which interacts with the specified residue.
Hb: hydrophobic interaction.
H: hydrogen bond.
Pi:
[figures omitted; refer to PDF]
2.3. Candidate Selection
Docking scores, predicted bioactivities, and corresponding consensus scores are detailed in Table 1, respectively. In general, candidates with higher dock scores than Tacrolimus were also predicted to have similar if not better bioactivities. Daphnetoxin (TCM origin: Daphne giraldii Nitsche), 20-O-(2′E,4′E-decadienoyl)ingenol (TCM origin: Euphorbia kansui), and Lythrancine II (TCM origin: Lythrum salicaria) were selected as research candidates based on the consensus voting system (Table 1). Structural scaffolds of Tacrolimus and the selected candidates are illustrated in Figure 3.
[figures omitted; refer to PDF]
2.4. Molecular Dynamics Simulation
2.4.1. Stability Profile Analysis
Figure 4 illustrates protein-ligand complex stability in terms of RMSD and total energy. In general, the complex RMSD of Daphnetoxin and 20-O-(2′E,4′E-decadienoyl)ingenol was higher (1.5
[figures omitted; refer to PDF]
2.5. RMS Fluctuation (RMSF) of Protein Residues
Residues contributing to complex structural fluctuations can be assessed by root mean square fluctuations (RMSFs) of each residue (Figure 5). Tacrolimus and Lythrancine II have lower RMSFs and similar trends. 20-O-(2′E,4′E-Decadienoyl)ingenol has the highest RMSF of the four, implying greater fluctuations. Daphnetoxin falls between Tacrolimus, Lythrancine II, and 20-O-(2′E,4′E-decadienoyl)ingenol. Analysis of the RMSF values shows that great differences among the protein complexes were observed at Tyr113, a residue near the binding site.
[figure omitted; refer to PDF]2.6. H-Bond Network during MD Simulation
Next we investigated formation and stability of H-bonds under dynamic conditions. Individual occupancies of detected H-bonds per ligand are detailed in Table
[figures omitted; refer to PDF]
Daphnetoxin forms several discontinuous H-bonds with residues of FKBP52. However, the H-bond occupancies at Tyr113 with O16, O18, and H52 of Daphnetoxin are 31.47%, 2.00%, and 9.89%, respectively. It suggests that the Tyr113 might play an important role in protein-ligand complex stability.
20-O-(2′E,4′E-Decadienoyl)ingenol and the residues of FKBP52 form some unstable H-bonds. Tyr113, O21 and H60 form only transient H-bonds during MD, whereas Asp68 and Lys121 form a stable H-bond with 20-O-(2′E,4′E-decadienoyl)ingenol from 13 to 20 ns.
Lythrancine II and residues of FKBP52 form two stable H-bonds: Ile87 forms a stable H-bond at O36 from 6 to 17.5 ns. Tyr113 and O19 on Lythrancine II form a stable H-bond from 0 to 17.5 ns. The two residues are important for protein-ligand complex stability.
2.7. Torsion MD Simulation
Analysis of torsion helps us understand the stability of ligand and protein binding (Figure 7, Figure S1; see Supplementary Material available online at http://dx.doi.org/10.1155/2014/364819). Torsion (A1) of Tacrolimus changes from significant unstable rotations (0–7.5 ns) to stable rotations (7.5–17 ns), and change to unstable rotations after 17 ns. Comparing to H-bond simulation results, we can assume that Glu110 forms a more stable H-bond. Torsion (A4) is more stable from 2 to 6 ns. During this period, H-bond simulation shows the formation of an H-bond between FKBP52 and Ser118. Torsion (A5) shows some changes, but the small torsion fluctuations at (A2), (A3), (A6), and (A7) show that they have higher stability. Torsion (B8) of Daphnetoxin shows transient stability. Torsion (B9) is stable with small fluctuations. Torsion (B10) has smaller torsion fluctuations from 0 to 10 ns; torsion (B11) also has very small changes. During this period, they separately form H-bonds with Glu85 and Ile87 on FKBP52. Torsion (B12) has small fluctuations from 8 to 18 ns, and torsion (B13) is relatively stable; these locations form H-bonds with Tyr57 and Tyr113. Violent fluctuations at (C14) were observed for 20-O-(2′E,4′E-decadienoyl)ingenol, a phenomenon probably due to the lack of stabilizing bonds with FKBP52. Torsions (C15) and (C16) only form short lived H-bonds which could be attributed to the large torsion changes. Torsion (C17) is relatively stable. With the exception of isolated large fluctuations, torsion (D18) is relatively stable, probably due to the formation of H-bond with Ser118. Torsion (D19) is relatively stable up until 12.5 ns when increased fluctuation was observed. Torsion (D20) is stable, most likely due to the H-bonds formed between Lythrancine II and Tyr57 and Asp68. Torsions (D21) and (D22) had large fluctuations.
[figures omitted; refer to PDF]
2.8. DSSP MD Simulation
Protein functionality is affected by the tertiary structure formed by the secondary structures. The conformation changes of protein secondary structures for each time frame can be computed by the DSSP program. Changes in the secondary protein structure when Tacrolimus, Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, or Lythrancine II is bound to FKBP52 are illustrated in Figure 8. Tacrolimus, Daphnetoxin, and Lythrancine II have an α-helix composition of approximately 6%, but that of 20-O-(2′E,4′E-decadienoyl)ingenol ranges between 6 and 10%. Compositions of β-sheets and turns in the four complexes are approximately 38-39% and 15–20%, but an increase in turn structures was observed in 20-O-(2′E,4′E-decadienoyl)ingenol from 12 ns. From Figure 8 we can observe more significant changes at residue 51. In Tacrolimus, residue 51 is a bend, but in Daphnetoxin, it turns into a turn at 0-1 ns and 12–18 ns. Similarly in 20-O-(2′E,4′E-decadienoyl)ingenol, the residue becomes a turn after 1 ns and after 11 ns. In Lythrancine II, it first appears as a turn from 0 to 2 ns, then turns to bend, and then fluctuates between bends and turns. Comparisons to RMSFs show that RMSF of residue 51 is also larger.
[figures omitted; refer to PDF]
2.9. Solvent Accessible Surface Area (SASA) in MD Simulation
Solvent accessible surface area (SASA) analysis measures the interaction between complexes and solvents. SASAs of the ligand-protein complexes are 140.788 nm/NS2, 141.028 nm/NS2, 141.745 nm/NS2, and 140.004 nm/NS2 for Tacrolimus, Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II, respectively. Average of the ligand SASAs is 9.068 nm/NS2, 5.669 nm/NS2, 6.827 nm/NS2, and 5.978 nm/NS2, respectively. According to the results of proteins and ligands, both the values of 20-O-(2′E,4′E-decadienoyl)ingenol are the greatest, suggesting that it should be accessible for solvents and have more interaction with solvents. In addition, we observed that the SASA values for the four protein complexes (Figure 9(a)) and the ligands (Figure 9(b)) during MD is relatively stable, indicating no significant changes in the protein structure.
[figures omitted; refer to PDF]
2.10. Radius of Gyration (Rg) in MD Simulation
Radius of gyration (Rg) enables one to assess the compactness changes of a ligand-protein complex. Average and maximum values for Tacrolimus, Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II are listed in Table 3. No significant changes were observed for
Table 3
Radius of gyration (Rg) calculated for the control and TCM candidates.
Radius of gyration (Rg; nm) | ||||||||
---|---|---|---|---|---|---|---|---|
Tacrolimus | Daphnetoxin | 20-O-(2′E,4′E-Decadienoyl)ingenol | Lythrancine II | |||||
Average | Maximum | Average | Maximum | Average | Maximum | Average | Maximum | |
Protein | 2.154 | 2.210 | 2.179 | 2.219 | 2.178 | 2.237 | 2.157 | 2.200 |
Mainchain | 2.140 | 2.202 | 2.162 | 2.205 | 2.165 | 2.226 | 2.142 | 2.187 |
Sidechain | 2.167 | 2.216 | 2.160 | 2.202 | 2.190 | 2.247 | 2.170 | 2.212 |
Backbone | 2.137 | 2.200 | 2.160 | 2.202 | 2.163 | 2.223 | 2.140 | 2.185 |
C alpha | 2.141 | 2.204 | 2.165 | 2.207 | 2.167 | 2.227 | 2.144 | 2.189 |
Ligand | 0.541 | 0.575 | 0.378 | 0.384 | 0.570 | 0.558 | 0.394 | 0.413 |
2.11. Cluster MD Simulation
MD simulations from 15 to 20 ns were used for cluster analysis (Figure 10 and Figure
[figures omitted; refer to PDF]
Comparing the conformation at 0 ns and the representative conformation for each complex is shown in Figure 11. Tacrolimus initially form an H-bond with Tyr113, but this H-bond is lost following stabilization. Moreover, the movement of Tacrolimus into the interior of the protein complex (Figure 11(a)) and the higher ligand RMSD (Figure 4) were observed. Therefore, the significant change of Tacrolimus might cause the H-bond to disappear and impact the interaction between Tacrolimus and FKBP52 deeply. The complex RMSD of Daphnetoxin is higher, but the ligand RMSD reach dynamic equilibrium with small fluctuations (Figure 4), it might cause Daphnetoxin to form two H-bonds with Tyr113, and this was maintained throughout MD (Figure 11(b)). The conformation of 20-O-(2′E,4′E-decadienoyl)ingenol and protein shows great differences between 0 ns and 19.58 ns at MD simulation. Especially, the position and the side chain of ligand have obvious changes. Comparing to Figure 4, its complex RMSD and ligand RMSD are higher and show larger fluctuation. As a result, 20-O-(2′E,4′E-decadienoyl)ingenol forms H-bonds with Tyr113 and Lys121 at the beginning of MD, but only the bond with Lys121 was maintained following stabilization (Figure 11(c)). Lythrancine II forms H-bonds with Asp68, Arg73, and Tyr113 at the beginning of MD, but these bonds are lost after stabilization, even though its complex RMSD and ligand RMSD are not higher. It might be attributed to the huge conformation change of the ligand. Similar to Tacrolimus, Lythrancine II was observed to be more embedded in the interior of the protein complex (Figure 11(d)).
[figures omitted; refer to PDF]
3. Discussion
Discrepancies between docking and MD simulation results were observed. During docking, Tacrolimus, Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II were predicted to form H-bonds with Tyr113 of FKBP52. Nonetheless, occupancies of the H-bond with Tyr113 varied greatly. Similarly, RMSF also shows that Tyr113 has great differences among the four complexes. In Table
From other MD simulation analysis we can make the following summarizations. Despite the larger protein complex RMSD of Daphnetoxin compared to Tacrolimus, the H-bond with Tyr113 was presented during both docking and MD simulation, indicating a better affinity towards FKBP52 than Tacrolimus. Summarizing the results from RMSF, SASA, and Rg, the protein complex of Daphnetoxin is stable during MD simulation. 20-O-(2′E,4′E-Decadienoyl)ingenol has higher complex and ligand RMSDs than Tacrolimus and also has more prominent changes in secondary structure changes, RMSF, and Rg. This could be related to the hydrophobic hydrocarbon backbone in its structure. However, 20-O-(2′E,4′E-decadienoyl)ingenol forms more hydrophobic interactions than Tacrolimus and is able to maintain the H-bond with Lys121 following stabilization, suggesting affinity with FKBP52. In addition, 20-O-(2′E,4′E-decadienoyl)ingenol has maximum values of SASA analysis of proteins and ligands, indicating that it should be accessible for solvents and have more interaction with solvents. Though Lythrancine II do not have H-bonds in the stabilized MD simulation conformation, its complex RMSD, ligand RMSD, RMSF, and Rg trends are most similar to Tacrolimus. In addition, movement of Lythrancine II towards the protein interior is also reminiscent of Tacrolimus. Therefore, Lythrancine II shows the most similar ligand-protein conformation with Tacrolimus. From this study, we can also see the importance of MD simulation in computer-aided drug designs.
4. Conclusion
We aimed to utilize structure-based and ligand-based methods to screen for high affinity lead compounds for FKBP52 as alternatives for Tacrolimus to develop better non-amyloid-beta-centric therapies for Alzheimer’s disease. Utilizing docking and highly reliable QSAR models, Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II were selected as potential candidates. Docking, Ligplot, and MD show that Daphnetoxin and Lythrancine II both have better affinity to FKBP52 than Tacrolimus. For 20-O-(2′E,4′E-decadienoyl)ingenol, modifying its long chain hydrocarbon could improve interaction with FKBP52. Based on these results, there is possibility to utilize Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II as backbone structures for modification as a new approach for designing Alzheimer’s disease.
5. Material and Methods
5.1. Data Collection
The FKBP52 crystal structure (PDB ID: 1Q1C) was downloaded from Protein Data Bank [60] and corrected for H-atoms using the Prepare Ligands module in Discovery Studio Client v2.5 (DS2.5; Accelrys Inc., San Diego, CA). Ligand structures and activity data from Gopalakrishnan’s study [61] were used to construct quantitative structure-activity relationship (QSAR) prediction models for FKBP52. TCM ligands used for virtual screening were downloaded from TCM Database@Taiwan [62]. Tacrolimus (FK506), the clinically used immunosuppressant, was used as the control.
5.2. Molecular Docking
Prior to virtual screening, Change Ionization in the Prepare Ligands module (DS 2.5) was applied to adjust ionization states of the downloaded ligands. Downloaded ligands and Tacrolimus were docked to the PPIase domain of FKBP52 under a forcefield of Chemistry at HARvard Molecular Mechanics (CHARMm) using the LigandFit module. Ligands were evaluated based on structural compatibility to the PPIase domain and dock score was selected as the scoring function.
5.3. Support Vector Machine (SVM) and Multiple Linear Regression (MLR) Prediction Models
The 37 ligands [61] were randomly divided into a training set of 30 compounds and a test set of 7 compounds. Molecular properties of each ligand were calculated with Calculate Molecular Properties and the ten representative descriptors most related to bioactivity were determined by Genetic Function Approximation (GFA). Different prediction models were constructed with the representative descriptors and the strength of each model ranked by a square correlation coefficient (
5.4. Bayesian Network
Training/test set groups and representative descriptors determined by GFA were also used to construct the Bayesian Network model. According to distribution characteristics, descriptors and pIC50 were discretized into a maximum of five categories. Linear regression analysis for each pIC50 category in the training dataset was then applied. Banjo [63] was used to discover relationships among the representative descriptors and pIC50 values. The algorithm for predicting pIC50 was written in MATLAB codes integrating Banjo and Bayes Net Toolbox (BNT; https://code.google.com/p/bnt/). Following validation, the algorithm was applied to predict bioactivity of selected TCM candidates.
5.5. Candidate Selection Criteria
All TCM ligands with dock scores lower than Tacrolimus were eliminated. Next, a consensus voting system was used for candidate selection. The three highest scoring candidates in dock score, SVM, or MLR were given a score of “1”; candidates with the three highest Bayesian Network scores were given a score of “2.” Sum of scores was calculated and the three ligands with the highest scores were selected as candidates for further analysis.
5.6. Molecular Dynamics (MD) Simulation
Candidates and Tacrolimus were prepared with SwissParam (http://swissparam.ch/) [64] prior to MD. MD simulation was conducted with GROMACS 4.0.7 under the forcefield of CHARMm27. Distance of protein to box boundaries was 1.2 nm; solvate TIP3P water model was then added to the system. Complex charges were neutralized with sodium and chloride ions using 0.145 M salt model. Simulation was conducted at 310 K under a pressure of 1 bar. Each complex was minimized with 5,000 steps of Steepest Descent, and the final minimized structure was used as the initial structure for MD simulation. Electrostatic interactions were calculated with Particle-Mesh-Ewald (PME) [65]. Equilibration protocol was used to restrain and relax protein-ligand position; first-order kinetics was started from 300 K. Minimized system was used to simulate a five-thousand ps configuration production. MD simulation was conducted for 20 ns with time steps of 2 fs under PME.
Authors’ Contribution
Kai Hsin Liao, Kuen-Bao Chen, and Wen-Yuan Lee contributed equally to this work.
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 (101-ASIA-59, ASIA100-CMU-2, and ASIA101-CMU-2), and China Medical University Hospital (DMR-103-058, DMR-103-001, and DMR-103-096, CMU102-BC-3). 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).
[1] R. L. Adams, P. L. Craig, O. A. Parsons, "Neuropsychology of dementia," Neurologic Clinics, vol. 4 no. 2, pp. 387-404, 1986.
[2] J. L. Cummings, "Subcortical dementia. Neuropsychology, neuropsychiatry, and pathophysiology," British Journal of Psychiatry, vol. 149, pp. 682-697, 1986.
[3] G. N. Savla, B. W. Palmer, "Neuropsychology in Alzheimer's disease and other dementia research," Current Opinion in Psychiatry, vol. 18 no. 6, pp. 621-627, 2005.
[4] B. C. Jost, G. T. Grossberg, "The evolution of psychiatric symptoms in Alzheimer's disease: a natural history study," Journal of the American Geriatrics Society, vol. 44 no. 9, pp. 1078-1081, 1996.
[5] D. J. Selkoe, "Alzheimer's disease is a synaptic failure," Science, vol. 298 no. 5594, pp. 789-791, DOI: 10.1126/science.1074069, 2002.
[6] D. G. Harwood, W. W. Barker, R. L. Ownby, R. Duara, "Relationship of behavioral and psychological symptoms to cognitive impairment and functional status in Alzheimer's disease," International Journal of Geriatric Psychiatry, vol. 15 no. 5, pp. 393-400, 2000.
[7] M. Haupt, A. Kurz, M. Jänner, "A 2-year follow-up of behavioural and psychological symptoms in Alzheimer's disease," Dementia and Geriatric Cognitive Disorders, vol. 11 no. 3, pp. 147-152, 2000.
[8] M. W. Hooper, F. S. Vogel, "The limbic system in Alzheimer's disease. A neuropathologic investigation," The American Journal of Pathology, vol. 85 no. 1, 1976.
[9] R. S. Miletich, T. N. Chase, R. Blesa, "Limbic system dysfunction in Alzheimer's disease," Journal of Neurology Neurosurgery and Psychiatry, vol. 59 no. 4, pp. 450-451, 1995.
[10] J. A. Harris, N. Devidze, L. Verret, K. Ho, B. Halabisky, M. T. Thwin, D. Kim, P. Hamto, I. Lo, G.-Q. Yu, J. J. Palop, E. Masliah, L. Mucke, "Transsynaptic progression of Amyloid- β -induced neuronal dysfunction within the entorhinal-hippocampal network," Neuron, vol. 68 no. 3, pp. 428-441, DOI: 10.1016/j.neuron.2010.10.020, 2010.
[11] H. Zempel, E.-M. Mandelkow, "Linking amyloid- β and tau: amyloid- β induced synaptic dysfunction via local wreckage of the neuronal cytoskeleton," Neurodegenerative Diseases, vol. 10 no. 1–4, pp. 64-72, DOI: 10.1159/000332816, 2012.
[12] V. Nimmrich, U. Ebert, "Is alzheimer's disease a result of presynaptic failure?—Synaptic dysfunctions induced by oligomeric p-amyloid," Reviews in the Neurosciences, vol. 20 no. 1, 2009.
[13] M. A. Busche, X. Chen, H. A. Henning, "Critical role of soluble amyloid-beta for early hippocampal hyperactivity in a mouse model of Alzheimer's disease," Proceedings of the National Academy of Sciences of the United States of America, vol. 109 no. 22, pp. 8740-8745, 2012.
[14] R. B. Maccioni, J. P. Muñoz, L. Barbeito, "The molecular bases of Alzheimer's disease and other neurodegenerative disorders," Archives of Medical Research, vol. 32 no. 5, pp. 367-381, DOI: 10.1016/S0188-4409(01)00316-2, 2001.
[15] J. P. Steiner, T. M. Dawson, M. Fotuhi, C. E. Glatt, A. M. Snowman, N. Cohen, S. H. Snyder, "High brain densities of the immunophilin FKBP colocalized with calcineurin," Nature, vol. 358 no. 6387, pp. 584-587, DOI: 10.1038/358584a0, 1992.
[16] J. Giustiniani, M. Sineus, E. Sardin, O. Dounane, M. Panchal, V. Sazdovitch, C. Duyckaerts, B. Chambraud, E.-E. Baulieu, "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, DOI: 10.3233/JAD-2011-111895, 2012.
[17] B. Chambraud, E. Sardin, J. Giustiniani, O. Dounane, M. Schumacher, M. Goedert, E.-E. Baulieu, "A role for FKBP52 in Tau protein function," Proceedings of the National Academy of Sciences of the United States of America, vol. 107 no. 6, pp. 2658-2663, DOI: 10.1073/pnas.0914957107, 2010.
[18] W. Cao, M. Konsolaki, "FKBP immunophilins and Alzheimer's disease: a chaperoned affair," Journal of Biosciences, vol. 36 no. 3, pp. 493-498, DOI: 10.1007/s12038-011-9080-7, 2011.
[19] C. B. Kang, Y. Hong, S. Dhe-Paganon, H. S. Yoon, "FKBP family proteins: immunophilins with versatile biological functions," NeuroSignals, vol. 16 no. 4, pp. 318-325, DOI: 10.1159/000123041, 2008.
[20] W. B. Pratt, D. O. Toft, "Steroid receptor interactions with heat shock protein and immunophilin chaperones," Endocrine Reviews, vol. 18 no. 3, pp. 306-360, DOI: 10.1210/er.18.3.306, 1997.
[21] S. L. Schreiber, "Chemistry and biology of the immunophilins and their immunosuppressive ligands," Science, vol. 251 no. 4991, pp. 283-287, 1991.
[22] Y. Yoshiyama, M. Higuchi, B. Zhang, S.-M. Huang, N. Iwata, T. Saido, J. Maeda, T. Suhara, J. Q. Trojanowski, V. M.-Y. Lee, "Synapse loss and microglial activation precede tangles in a P301S tauopathy mouse model," Neuron, vol. 53 no. 3, pp. 337-351, DOI: 10.1016/j.neuron.2007.01.010, 2007.
[23] R. B. Birge, S. Wadsworth, R. Akakura, H. Abeysinghe, R. Kanojia, M. Macielag, J. Desbarats, M. Escalante, K. Singh, S. Sundarababu, K. Parris, G. Childs, A. August, J. Siekierka, D. E. Weinstein, "A role for schwann cells in the neuroregenerative effects of a non-immunosuppressive FK506 derivative, JNJ460," Neuroscience, vol. 124 no. 2, pp. 351-366, DOI: 10.1016/j.neuroscience.2003.10.013, 2004.
[24] B. G. Gold, V. Densmore, W. Shou, M. M. Matzuk, H. S. Gordon, "Immunophilin FK506-binding protein 52 (not FK506-binding protein 12) mediates the neurotrophic action of FK506," Journal of Pharmacology and Experimental Therapeutics, vol. 289 no. 3, pp. 1202-1210, 1999.
[25] B. G. Gold, J. E. Villafranca, "Neuroimmunophilin ligands: the development of novel neuroregenerative/ neuroprotective compounds," Current topics in medicinal chemistry, vol. 3 no. 12, pp. 1368-1375, 2003.
[26] S. Barik, "Immunophilins: for the love of proteins," Cellular and Molecular Life Sciences, vol. 63 no. 24, pp. 2889-2900, DOI: 10.1007/s00018-006-6215-3, 2006.
[27] E. Karran, M. Mercken, B. D. Strooper, "The amyloid cascade hypothesis for Alzheimer's disease: an appraisal for the development of therapeutics," Nature Reviews Drug Discovery, vol. 10 no. 9, pp. 698-712, DOI: 10.1038/nrd3505, 2011.
[28] T. E. Golde, L. Petrucelli, J. Lewis, "Targeting A β and tau in Alzheimer's disease, an early interim report," Experimental Neurology, vol. 223 no. 2, pp. 252-266, DOI: 10.1016/j.expneurol.2009.07.035, 2010.
[29] W.-L. Liao, F.-J. Tsai, "Personalized medicine: a paradigm shift in healthcare," BioMedicine, vol. 3 no. 2, pp. 66-72, 2013.
[30] F.-J. Tsai, "Biomedicine brings the future nearer," BioMedicine, vol. 1 no. 1,DOI: 10.1016/j.biomed.2011.10.007, 2011.
[31] F.-J. Tsai, "Rare diseases: a mysterious puzzle," BioMedicine, vol. 3 no. 2, article 65, 2013.
[32] I. C. Chou, W.-D. Lin, C.-H. Wang, "Möbius syndrome in a male with XX/XY mosaicism," BioMedicine, vol. 3 no. 2, pp. 102-104, 2013.
[33] C.-C. Lee, C.-H. Tsai, L. Wan, "Increased incidence of Parkinsonism among Chinese with β -glucosidase mutation in central Taiwan," BioMedicine, vol. 3 no. 2, pp. 92-94, 2013.
[34] D.-Y. Lin, F.-J. Tsai, C.-H. Tsai, C.-Y. Huang, "Mechanisms governing the protective effect of 17 β -estradiol and estrogen receptors against cardiomyocyte injury," BioMedicine, vol. 1 no. 1, pp. 21-28, DOI: 10.1016/j.biomed.2011.10.004, 2011.
[35] C.-H. Wang, W.-D. Lin, D.-T. Bau, "Appearance of acanthosis nigricans may precede obesity: an involvement of the insulin/IGF receptor signaling pathway," BioMedicine, vol. 3 no. 2, pp. 82-87, 2013.
[36] Y.-M. Chang, B. K. Velmurugan, W.-W. Kuo, "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.
[37] I. C. Chou, W.-D. Lin, C.-H. Wang, "Association analysis between Tourette's syndrome and two dopamine genes (DAT1, DBH) in Taiwanese children," BioMedicine, vol. 3 no. 2, pp. 88-91, 2013.
[38] W.-Y. Lin, H.-P. Liu, J.-S. Chang, "Genetic variations within the PSORS1 region affect Kawasaki disease development and coronary artery aneurysm formation," BioMedicine, vol. 3 no. 2, pp. 73-81, 2013.
[39] I. Cohen, M. Tagliaferri, D. Tripathy, "Traditional Chinese medicine in the treatment of breast cancer," Seminars in Oncology, vol. 29 no. 6, pp. 563-574, DOI: 10.1053/sonc.2002.50005, 2002.
[40] S. J. Leuenroth, D. Okuhara, J. D. Shotwell, G. S. Markowitz, Z. Yu, S. Somlo, C. M. Crews, "Triptolide is a traditional Chinese medicine-derived inhibitor of polycystic kidney disease," Proceedings of the National Academy of Sciences of the United States of America, vol. 104 no. 11, pp. 4389-4394, DOI: 10.1073/pnas.0700499104, 2007.
[41] X. Gong, N. J. Sucher, "Stroke therapy in traditional Chinese medicine (TCM): prospects for drug discovery and development," Phytomedicine: International Journal of Phytotherapy and Phytopharmacology, vol. 9 no. 5, pp. 478-484, 2002.
[42] B. La Cour, P. Mølgaard, Z. Yi, "Traditional Chinese medicine in treatment of hyperlipidaemia," Journal of Ethnopharmacology, vol. 46 no. 2, pp. 125-129, DOI: 10.1016/0378-8741(95)01234-5, 1995.
[43] W. L. Li, H. C. Zheng, J. Bukuru, N. de Kimpe, "Natural medicines used in the traditional Chinese medical system for therapy of diabetes mellitus," Journal of Ethnopharmacology, vol. 92 no. 1,DOI: 10.1016/j.jep.2003.12.031, 2004.
[44] T. T. Chang, M. F. Sun, H. Y. Chen, F.-J. Tsai, M. Fisher, J.-G. Lin, C. Y. C. Chen, "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.
[45] K. C. Chen, M. F. Sun, S. C. Yang, S.-S. Chang, H.-Y. Chen, F.-J. Tsai, C. Y. C. Chen, "Investigation into Potent Inflammation Inhibitors from Traditional Chinese Medicine," Chemical Biology and Drug Design, vol. 78 no. 4, pp. 679-688, DOI: 10.1111/j.1747-0285.2011.01202.x, 2011.
[46] J. Yin, H. Zhang, J. Ye, "Traditional Chinese medicine in treatment of metabolic syndrome," Endocrine, Metabolic and Immune Disorders—Drug Targets, vol. 8 no. 2, pp. 99-111, DOI: 10.2174/187153008784534330, 2008.
[47] D. Normile, "The new face of traditional Chinese medicine," Science, vol. 299 no. 5604, pp. 188-190, DOI: 10.1126/science.299.5604.188, 2003.
[48] T. Y. Tsai, K. W. Chang, 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, DOI: 10.1007/s10822-011-9438-9, 2011.
[49] K. W. Chang, T. Y. Tsai, K. C. Chen, S.-C. Yang, H.-J. Huang, T.-T. Chang, M.-F. Sun, H.-Y. Chen, F.-J. Tsai, C. Y. C. Chen, "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.
[50] C. Y. Chen, 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, DOI: 10.1016/j.jmgm.2010.04.002, 2010.
[51] 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, DOI: 10.1016/j.jmgm.2009.08.004, 2009.
[52] K.-C. Chen, 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, DOI: 10.1039/c0sm01548b, 2011.
[53] S. S. Chang, H. J. Huang, 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,DOI: 10.1371/journal.pcbi.1002315, 2011.
[54] S.-C. Yang, S.-S. Chang, H.-Y. Chen, C. Y.-C. Chen, "Identification of potent EGFR inhibitors from TCM Database@Taiwan," PLoS Computational Biology, vol. 7 no. 10,DOI: 10.1371/journal.pcbi.1002189, 2011.
[55] W. I. Tou, C. Y. C. Chen, "In silico investigation of potential Src kinase ligands from traditional chinese medicine," PLoS ONE, vol. 7 no. 3,DOI: 10.1371/journal.pone.0033728, 2012.
[56] W. I. Tou, S. S. Chang, C. C. Lee, C. Y. C. Chen, "Drug design for neuropathic pain regulation from traditional Chinese Medicine," Scientific Reports, vol. 3, 2013.
[57] C. Y. C. Chen, "A novel integrated framework and omproved methodology of computer-aided drug design," Current Topics in Medicinal Chemistry, vol. 13 no. 9, pp. 965-988, 2013.
[58] C. Y. C. Chen, W. L. Tou, "How to design a drug for the disordered proteins," Drug Discovery Today, vol. 18 no. 19-20, pp. 910-915, 2013.
[59] W. I. Tou, C. Y. Chen, "May disordered protein cause serious drug side effect?," Drug Discovery Today, vol. 19 no. 4, pp. 367-372, 2014.
[60] B. Wu, P. Li, Y. Liu, Y. Ding, C. Shu, S. Ye, M. Bartlam, B. Shen, Z. Rao, "3D structure of human FK506-binding protein 52: implication for the assembly of the glucocorticoid receptor/Hsp90/immunophilin heterocomplex," Proceedings of the National Academy of Sciences of the United States of America, vol. 101 no. 22, pp. 8348-8353, DOI: 10.1073/pnas.0305969101, 2004.
[61] R. Gopalakrishnan, C. Kozany, Y. S. Wang, "Exploration of pipecolate sulfonamides as binders of the FK506-binding proteins 51 and 52," Journal of Medicinal Chemistry, vol. 55 no. 9, pp. 4123-4131, 2012.
[62] 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,DOI: 10.1371/journal.pone.0015939, 2011.
[63] J. Yu, V. A. Smith, P. P. Wang, A. J. Hartemink, E. D. Jarvis, "Advances to Bayesian network inference for generating causal networks from observational biological data," Bioinformatics, vol. 20 no. 18, pp. 3594-3603, DOI: 10.1093/bioinformatics/bth448, 2004.
[64] V. Zoete, M. A. Cuendet, A. Grosdidier, O. Michielin, "SwissParam: a fast force field generation tool for small organic molecules," Journal of Computational Chemistry, vol. 32 no. 11, pp. 2359-2368, DOI: 10.1002/jcc.21816, 2011.
[65] T. A. Darden, L. G. Pedersen, "Molecular modeling: an experimental tool," Environmental Health Perspectives, vol. 101 no. 5, pp. 410-412, 1993.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2014 Kai Hsin Liao et al. This is an open access article distributed under the Creative Commons Attribution License (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License. https://creativecommons.org/licenses/by/4.0/
Abstract
Alzheimer’s disease is a neurodegenerative disease that was conventionally thought to be related to the sedimentation of beta-amyloids, but drugs designed according to this hypothesis have generally failed. That FKBP52 can reduce the accumulation of tau proteins, and that Tacrolimus can reduce the pathological changes of tau proteins are new directions away from the long held amyloid-beta-centric concept. Therefore, the screening of traditional Chinese medicine compounds for those with higher affinity towards FKBP52 than Tacrolimus may be a new direction for treating Alzheimer’s disease. This study utilizes ligand-based and structure-based methods as the foundation. By utilizing dock scores and the predicted pIC50 from SVM, MLR, and Bayesian Network, several TCM compounds were selected for further analysis of their protein-ligand interactions. Daphnetoxin has higher affinity and complex structure stability than Tacrolimus; Lythrancine II exhibits the most identical trends in FKBP52 interactions as Tacrolimus, and 20-O-(2′E,4′E-decadienoyl)ingenol may be further modified at its hydrocarbon chain to promote interaction with FKBP52. In addition, we observed the residue Tyr113 of FKBP52 may play a key role in protein-ligand interaction. Our results indicate that Daphnetoxin, 20-O-(2′E,4′E-decadienoyl)ingenol, and Lythrancine II may be starting points for further modification as a new type of non-amyloid-beta-centric drug for Alzheimer’s disease.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer