| Literature DB >> 34368579 |
Shisanupong Anukanon1, Pornkanok Pongpamorn1, Wareepat Tiyabhorn1, Jaruwan Chatwichien1, Worawat Niwetmarin1, Richard B Sessions2, Somsak Ruchirawat1,3,4, Nopporn Thasana1,3,4.
Abstract
Huperzine A (1, Hup A), a lycodine-type Lycopodium alkaloid isolated from Thai clubmosses Huperzia squarrosa (G. Forst.) Trevis., H. carinata (Desv. ex. Poir.) Trevis., H. phlegmaria (L.), and Phlegmariurus nummulariifolius (Blume) Chambers (Lycopodiaceae), exerts inhibitory activity on acetylcholinesterase, a known target for Alzheimer's disease therapy. This study investigated the structure-activity relationship of C(2)-functionalized and O- or N-methyl-substituted huperzine A derivatives. In silico-guided screening was performed to search for potential active compounds. Molecular docking analysis suggested that substitution at the C(2) position of Hup A with small functional groups could enhance binding affinity with AChE. Consequently, 12 C(2)-functionalized and four O- or N-methyl-substituted compounds were semi-synthesized and evaluated for their eeAChE and eqBChE inhibitory activities. The result showed that 2-methoxyhuperzine A (10) displayed moderate to high eeAChE inhibitory potency (IC50 = 0.16 μM) with the best selectivity over eqBChE (selectivity index = 3633). Notably, this work showed a case of which computational analysis could be utilized as a tool to rationally screen and design promising drug molecules, getting rid of impotent molecules before going more deeply on labor-intensive and time-consuming drug discovery and development processes.Entities:
Year: 2021 PMID: 34368579 PMCID: PMC8340427 DOI: 10.1021/acsomega.1c02875
Source DB: PubMed Journal: ACS Omega ISSN: 2470-1343
Figure 1Structure of (−)-huperzine A (1) and club mosses (a) H. squarrosa, (b) H. carinata, (c) H. phlegmaria, and (d) P. nummulariifolius.
Figure 2Molecular docking of HupA with AChE (4EY5) from BUDE and Sire Waterswap.
Figure 3Design strategy for novel C(2) HupA derivatives targeting AChE.
Molecular Docking Analysis of Compounds 1–17 against hAChE (4EY5) and hBChE (4XII) Using Autodock Software
| compound | estimated free energy of binding (kcal/mol) | estimated inhibition constant, | H-bond interaction residues (distance in Å) | estimated free energy of binding (kcal/mol) | estimated inhibition constant, | H-bond interaction residues (distance in Å) |
|---|---|---|---|---|---|---|
| –12.17 | 1.20 | Tyr337 (2.90), Tyr133 (2.75, 3.15) | –8.08 | 1.20 μM | Glu197 (2.50) | |
| –12.10 | 1.34 | Tyr133 (2.92) | –9.21 | 177.89 | ||
| –10.63 | 16.24 | Tyr133 (2.51), Gly120 (3.15) | –8.87 | 313.41 | Asp70 (2.48), Thr120 (2.36), Tyr332 (2.99) | |
| –12.00 | 1.62 | Tyr133 (2.87) | –8.65 | 453.73 | Glu197 (2.77) | |
| –12.20 | 1.13 | Tyr133 (2.90, 3.15) | 8.33 | 780.19 | Asp70 (2.38), Tyr332 (2.69) | |
| –12.31 | 0.95 | Tyr133 (2.78) | –8.92 | 290.98 | ||
| –10.25 | 30.64 | Tyr133 (2.92), Tyr337 (2.84) | –8.91 | 296.05 | ||
| –12.19 | 1.17 | Trp86 (2.98), Tyr133 (2.87) | –8.65 | 458.89 | Tyr128 (2.43), Glu197 (2.74) | |
| –12.47 | 0.73 | Trp86 (2.91), Tyr133 (2.89) | –8.86 | 321.81 | Tyr128 (2.43), Glu197 (2.78) | |
| –12.75 | 0.45 | Tyr133 (2.84), Glu202 (2.70) | –9.40 | 129.68 | Glu197 (2.51) | |
| –11.30 | 5.23 | Gly120 (2.75), Tyr133 (2.64) | –8.64 | 468.08 | Asp70 (2.36), Tyr (2.68) | |
| –11.04 | 8.13 | –8.89 | 302.29 | Asp70 (2.49), Tyr332 (2.67) | ||
| –8.45 | 643.06 | –8.77 | 370.76 | Glu197 (2.97), His438 (2.77) | ||
| –12.21 | 1.11 | Glu202 (2.53) | –8.99 | 258.05 | Glu197 (2.65) | |
| –8.39 | 713.34 | Ser203 (2.61), His447 (3.22) | –7.29 | 4.50 μM | Ser198 (3.29), His438 (2.68) | |
| –9.11 | 211.03 | Gly120 (2.99) | –7.38 | 3.88 μM | ||
| –7.63 | 2.56 μM | –7.55 | 2.91 μM | Ser79 (3.09) | ||
| donepezil | –12.42 | 0.78 | Phe295 (2.99) | –9.99 | 47.44 | |
Binding free energy and inhibitory constant results obtained from AutoDock 4.2.6 software.[25]
Molecular docking visualization and H-bonding measurement of each compound were obtained from UCSF Chimera Molecular graphics, and analyses were performed with UCSF Chimera.[26]
Donepezil, a known selective AChE inhibitor.[29]
Physiochemical Properties of Compounds 1–17
| compound | MW | HBA | HBD | RB | log BB | ||
|---|---|---|---|---|---|---|---|
| 242.32 | 2 | 2 | 0 | 58.88 | 1.88 | –0.45 | |
| 321.22 | 2 | 2 | 0 | 58.88 | 2.53 | –0.35 | |
| 368.22 | 2 | 2 | 0 | 58.88 | 2.56 | –0.34 | |
| 260.31 | 3 | 2 | 0 | 58.88 | 2.26 | –0.39 | |
| 256.35 | 2 | 2 | 0 | 58.88 | 2.22 | –0.39 | |
| 270.38 | 2 | 2 | 1 | 58.88 | 2.55 | –0.34 | |
| 267.33 | 3 | 2 | 0 | 82.67 | 1.71 | –0.82 | |
| 257.34 | 2 | 3 | 0 | 84.90 | 1.40 | –0.90 | |
| 258.32 | 3 | 3 | 0 | 79.11 | 1.58 | –0.79 | |
| 272.35 | 3 | 2 | 1 | 68.11 | 1.88 | –0.58 | |
| 286.38 | 3 | 2 | 2 | 68.11 | 2.21 | –0.53 | |
| 300.40 | 3 | 2 | 3 | 68.11 | 2.57 | –0.48 | |
| 314.43 | 3 | 2 | 4 | 68.11 | 2.91 | –0.43 | |
| 256.35 | 3 | 1 | 1 | 48.14 | 2.50 | –0.19 | |
| 284.40 | 3 | 0 | 2 | 25.36 | 3.09 | 0.23 | |
| 284.40 | 2 | 0 | 1 | 25.36 | 3.09 | 0.23 | |
| 270.37 | 2 | 1 | 1 | 36.10 | 2.92 | 0.05 | |
| donepezil | 381.51 | 4 | 0 | 6 | 38.77 | 3.85 | 0.15 |
| required parameters | <500 | <10 | <5 | <10 | <90 | 2–5 | >−1.00 |
Calculated using ChemBioDraw Ultra16.0. MW: molecular weight.
Calculated using SwissADME (http://www.swissadme.ch/).[30] HBA: number of hydrogen acceptors; HBD: number of hydrogen donors; RB: number of rotatable bonds; tPSA: total polar surface area; cLog P: log octanol/water partition coefficient.
log BB = −0.0148 × tPSA + 0.152 × cLogP + 0.139.[28]
Donepezil, a known selective AChE inhibitor.[29]
Required parameters necessary to fulfill appropriate physiochemical properties as judged appropriately according to Lipinski’s rules and those important for BBB permeation.[27]
Scheme 1Synthesis of C(2)-Functionalized HupA Derivatives 2–13
Reaction conditions: (a) NBS, DCM, RT, 24 h; (b) NIS, DCM, RT, 48 h; (c) 1-chloromethyl-4-fluoro-1,4-diazoniabicyclo[2.2.2]octane bis(tetrafluoroborate), Pd(dppf)Cl2, THF, reflux, 96 h; (d) CH3MgBr, Pd(dppf)Cl2, THF, reflux, 72 h; (e) Zn(CH3)2, Pd(dppf)Cl2, THF, reflux, 72 h; (f) Zn(CH2CH3)2, Pd(dppf)Cl2, THF, reflux, 72 h; (g) CuCN, K2CO3, DMF, 120 °C, 72 h; (h) NH3 (aq.), CuI, K2CO3, 140 °C, 72 h; (i) H2O, CuI, K2CO3, 130 °C, 72 h; (j) MeOH, CuI, K2CO3, 140 °C, 72 h; (k) EtOH, CuI, K2CO3, 140 °C, 72 h; (l) n-PrOH, CuI, K2CO3, 140 °C, 72 h; (m) n-BuOH, CuI, K2CO3, 140 °C, 72 h.
Scheme 2Synthesis of HupA Derivatives 14–17
Reaction conditions: (a) CH3I, Ag2CO3, DCM, RT, 45 h; (b) CH3I, Ag2CO3, DCM, RT, 72 h.
In Vitro Inhibitory Activity of Compounds 1–17 toward Cholinesterases
| cholinesterase
IC50 (μM)a | |||||||
|---|---|---|---|---|---|---|---|
| compound | R1 | R2 | R3 | R4 | SId | ||
| H | H | H | H | 0.03 ± 0.00 | 68.78 ± 03.15 | 2,273 | |
| Br | H | H | H | 2.56 ± 0.17 | 72.29 ± 04.57 | 28 | |
| I | H | H | H | 14.12 ± 0.84 | 325.22 ± 38.63 | 23 | |
| F | H | H | H | 1.23 ± 0.01 | 63.998 ± 03.63 | 52 | |
| Me | H | H | H | 0.35 ± 0.02 | 486.03 ± 29.71 | 1,388 | |
| Et | H | H | H | 0.57 ± 0.09 | 60.39 ± 02.54 | 105 | |
| CN | H | H | H | 33.88 ± 1.95 | 1,050.88 ± 498.89 | 31 | |
| NH2 | H | H | H | 0.32 ± 0.01 | 524.64 ± 41.56 | 1,634 | |
| OH | H | H | H | 0.53 ± 0.09 | 154.76 ± 04.14 | 293 | |
| OMe | H | H | H | 0.16 ± 0.02 | 599.15 ± 49.76 | 3,633 | |
| OEt | H | H | H | 0.22 ± 0.04 | 253.02 ± 12.24 | 1,154 | |
| O | H | H | H | 10.19 ± 0.50 | 87.27 ± 01.37 | 9 | |
| O | H | H | H | 109.65 ± 9.59 | 16.91 ± 01.23 | 0.2 | |
| H | H | 3.67 ± 0.28 | 903.66 ± 88.20 | 246 | |||
| Me | Me | n.i. | 1,348.75 ± 62.31 | ||||
| Me | Me | Me | n.i. | 258.70 ± 10.81 | |||
| Me | Me | 0.61 ± 0.02 | 223.57 ± 47.78 | 366 | |||
| donepezil | 0.23 ± 0.05 | 18.00 ± 00.79 | 79 | ||||
Mean IC50 (μM) ± SEM, n = 3; n.i. = no inhibition.
AChE (EC 3.1.1.7) from electric eel.
BChE (EC 3.1.1.8) from horse serum.
Selectivity index for AChE is defined as IC50 (BChE)/IC50 (AChE).
Donepezil, a known selective AChE inhibitor used as a positive control.[29]
Figure 4(a, b) Molecular docking and (c, d) binding mode analysis of compounds 10 and 13 toward AChE (4EY5).
Figure 5(a, b) Molecular docking and (c, d) binding mode analysis of compounds 10 and 13 toward BChE (4XII).
Figure 6Correlation of the estimated Ki value according to the docking program and calculated Ki according to in vitro results. Correlation is significant at the p < 0.001 level (two-tailed).
In Vitro Cytotoxicity of HupA Derivatives and Donepezil against IMR90 Cellsa
| %cytotoxicity | %cytotoxicity | ||||
|---|---|---|---|---|---|
| compound | 10 μM | 100 μM | compound | 10 μM | 100 μM |
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | 2.2 ± 2.4 | n.i. | 5.5 ± 9.7 | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | n.i. | n.i. | ||
| n.i. | n.i. | donepezil | n.i. | 29.0 ± 12.5 | |
IMR90 cells were treated with each compound at the indicated concentrations for 24 h. Cell viability was determined by the MTT assay.
The values are shown as mean ± SEM, n = 3. n.i. = no inhibition.
Significantly different from donepezil (100 μM) with p < 0.05.