Next Article in Journal
Electrophysiology and Behavioral Assessment of the New Molecule SMe1EC2M3 as a Representative of the Future Class of Triple Reuptake Inhibitors
Next Article in Special Issue
A Second Look at the Crystal Structures of Drosophila melanogaster Acetylcholinesterase in Complex with Tacrine Derivatives Provides Insights Concerning Catalytic Intermediates and the Design of Specific Insecticides
Previous Article in Journal
Comprehensive Analysis of the Copper Exchange Implemented in Ammonia and Protonated Forms of Mordenite Using Microwave and Conventional Methods
Previous Article in Special Issue
Synthesis and In Vitro Screening of Novel Heterocyclic β-d-Gluco- and β-d-Galactoconjugates as Butyrylcholinesterase Inhibitors
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 24
Issue 23
10.3390/molecules24234217
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery of Selective Butyrylcholinesterase (BChE) Inhibitors through a Combination of Computational Studies and Biological Evaluations

by
You Zhou
1,*,
Xin Lu
2,
Hongyu Yang
2,
Yao Chen
3,
Feng Wang
1,
Jifu Li
1,
Zhiran Tang
4,
Xifei Cheng
1,
Yingbin Yang
4,
Li Xu
1 and
Qingyou Xia
1,*
1
State Key Laboratory of Silkworm Genome Biology, College of Biotechnology, Southwest University, Chongqing 400715, China
2
School of Pharmacy, China Pharmaceutical University, Nanjing 210009, China
3
School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing 210023, China
4
School of Life Science, Southwest University, Chongqing 400715, China
*
Authors to whom correspondence should be addressed.
Molecules 2019, 24(23), 4217; https://doi.org/10.3390/molecules24234217
Submission received: 7 November 2019 / Revised: 14 November 2019 / Accepted: 17 November 2019 / Published: 20 November 2019
(This article belongs to the Special Issue Cholinesterase Inhibitors)
Browse Figures
Versions Notes

Abstract

:
As there are increased levels and activity of butyrylcholiesterase (BChE) in the late stage of Alzheimer’s disease (AD), development of selective BChE inhibitors is of vital importance. In this study, a workflow combining computational technologies and biological assays were implemented to identify selective BChE inhibitors with new chemical scaffolds. In particular, a pharmacophore model served as a 3D search query to screen three compound collections containing 3.0 million compounds. Molecular docking and cluster analysis were performed to increase the efficiency and accuracy of virtual screening. Finally, 15 compounds were retained for biological investigation. Results revealed that compounds 8 and 18 could potently and highly selectively inhibit BChE activities (IC50 values < 10 μM on human BChE, selectivity index BChE > 30). These active compounds with novel scaffolds provided us with a good starting point to further design potent and selective BChE inhibitors, which may be beneficial for the treatment of AD.

Graphical Abstract

1. Introduction

Alzheimer’s disease (AD), an age-related and progressive neurodegenerative disorder featuring memory loss and cognitive impairments, is the most common type of dementia among older adults [1]. It affects more than 30 million patients worldwide at present, and the prevalence of AD continues to increase due to population aging [2,3]. Additionally, the current cost of dementia is about a trillion US dollars a year, and this figure will rise to 2 trillion US dollars by 2030 [4]. Thus, AD has become one of the greatest public health issues worldwide and it severely impacts patients and their families [5].
So far, the pathogenesis of AD still remains unclear, while multiple hypotheses, such as the amyloid-β oligomer hypothesis, the cholinergic hypothesis, and the tau hypothesis, have been offered to explain the mechanism of AD development and set the stage for new therapeutic approaches against AD [6,7,8]. Over the past decade, several anti-AD drugs (tacrine, donepezil, and rivastigmine et al.) targeting the cholinergic dysfunction have been clinically employed for the treatment of AD based on the cholinergic hypothesis, according to which the decline of hippocampal and cortical levels of acetylcholine (ACh) contributes substantially to cognitive decline [9,10,11,12]. In addition, several lines of evidence also indicated that both cholinesterases (ChEs), named as acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), play an important role in amyloid-β (Aβ) aggregation during the senile plaque formation and using the inhibitors can decrease senile plaques [13]. Over the past decade, AChE inhibitors have been investigated extensively as anti-AD agents. However, the therapeutic effect of AChE inhibitors is limited, while side effects like nausea and vomiting induced by undesirable inhibition of peripheral ChEs also hinder the long-term application in clinical trials [14,15]. Recently, increasing evidence have indicated that BChE plays a pivotal role in regulating brain ACh level in the late stage of AD. In progressed AD, the AChE level in the brain decreases to 55%–67% of normal values while BChE increases to 165% of normal levels [16,17]. Moreover, AChE knockout mice models indicated that BChE can potentially substitute for AChE, maintaining normal cholinergic pathways in AChE nullizygous animals [18]. In addition, BChE inhibition is not accompanied by the peripheral adverse effects [19]. Therefore, development of highly potent and selective BChE inhibitors that are able to restore Ach levels in the brain, with much reduced peripheral side effects, represents a significant advancement [20].
Up to now, only a limited number of highly selective and potent BChE inhibitors have been reported, as the two isoforms human AChE (hAChE) and human BChE (hBChE) are highly homologous proteins [21] (Figure 1). Stanislav Gobec et al. reported a series of sulfonamide BChE inhibitors in 2016. Among them, compound 1 showed highly selective and nanomolar inhibitory activity against hBChE [22]. Further modification of compound 1 led to the generation of naphthamide compound 2 with slightly improved potency [23]. Additionally, 6-substituted 3,4,5,6-tetrahydroazepino[4,3-b]indol-1(2H)-one (THAI) compound 3 and 2-thiophenyl compound 4 were also highly potent and selective hBChE inhibitors [24,25]. Since BChE-targeting inhibition represents a promising therapeutic approach against AD, the discovery of potent and selective BChE inhibitors is urgent. Over the past decade, multiple X-ray crystallographic structures of complexes between hBChE proteins and hBChE inhibitors with various scaffolds have been determined, which sets the stage for virtual screening [26,27,28]. There are already several reports of using a structure-based virtual screening protocol and successfully obtained potent and highly selective hBChE inhibitors [16,29,30]. After identifying drug hits with novel scaffolds and desired properties by virtual screening, further studies like structure–activity relationship (SAR) studies based on enzyme or cell bioassay and crystallographic study could provide more useful information for structural modification and improvement [31,32,33].
Recently, Sun et al. reported a valid structure-based pharmacophore model which led to the discovery of three compounds displaying IC50 values < 2 μM on hBChE [34]. Considering the effectiveness of the pharmacophore model, further structure-based pharmacophore virtual screening (SBP VS) of three commercial compound collections (ChemBridge, ChEMBL, Vitas-M) was reported in this paper. The employed screening workflow was depicted in Figure 2. With the integration of docking and diversity analysis, 25 compounds were selected and 15 of them were purchased for bioassays. Finally, two compounds with different scaffolds were identified as selective hBChE inhibitors. The new emerging molecules with different scaffolds not only enriched the structural types of selective BChE inhibitors, but also may provide valuable starting structures for medicinal chemists to develop anti-AD drugs.

2. Results and Discussion

2.1. Virtual Screening

The construction and validation of the 3D pharmacophore model was reported in the previous paper [34]. Generally, the six-feature model was generated on the basis of 1-BChE complex (PDB ID: 5DYW) [22]. Specifically, the model consists of two hydrogen bond acceptor points, two hydrophobic points, one aromatic ring point, and one positive ionizable point (Figure 3). In this paper, it was used for virtual screening of three commercial compound collections: Vitas-M (1315684), ChEMBL (918887), ChemBridge (777451). We retained 18,498 compounds whose pharmacophore fit value above 3, after the first filtration step. Next, 1502 compounds with favorable physicochemical properties were identified by using drug-like descriptors of Lipinski and Veber rules [35]. In order to reduce the amount of compounds to handle prior to biological experiments, two complementary methods including docking and structural clustering were utilized. On one hand, the CDOCKER module of Discovery Studio 3.0 (DS) was applied to predict the binding modes of all selected virtual hits with BChE active site (PDB ID: 5DYW) to check their interactions with BChE binding site residues [36]. The docked compounds were ranked by the corresponding values of –CDOCKER energy. Careful visual inspection was performed, and the following interactions were considered significant: π–π interaction between the molecule and Trp231, Trp82, and Phe329, H-bonding interaction of the molecule with His438, and occupancy of the molecule in the acyl pocket. A set of 100 compounds were saved. On the other hand, the above-mentioned 1502 compounds were subjected to structural clustering. Ten clusters were produced based on the FCFP_6 fingerprints using the Cluster Ligand module of DS and twenty compounds with diverse scaffolds in each cluster were chosen. In consideration of the binding mode and structural diversity, 25 candidates were selected and 15 commercially available candidates (Figure 4) were purchased from Topscience (www.tsbiochem.com) for biological investigation.

2.2. ChEs Inhibitory Activities of Hit Compounds

Initial screening of the 15 potential inhibitors obtained from in silico studies was performed with human ChEs using a modified Ellman’s assay, and tacrine was used as the reference control (Table 1). The result indicated that compounds 8 and 19 exhibited over 50.0% inhibitory effects on both AChE and BChE at the concentration of 10 μM. Interestingly, compound 18 exhibited selective BChE inhibitory effect (BChE = 58.4% at 10 μM, AChE = 11.1% at 10 μM). Next, the dose-dependent inhibitory activities of compounds 8, 18, and 19 against BChE and AChE were tested at doses ranging from 10−4 to 10−9 M, and their IC50 values were calculated (Figure S1). The result demonstrated that three compounds showed great anti-BChE activities (BChE IC50 < 10 μM). Additionally, compounds 8 and 18 showed much better BChE selective index (SI BChE, AChE IC50/BChE IC50 > 30) than compound 19 (SI BChE = 6). To the best of our knowledge, compounds 8 and 18 were structurally different from the previously reported selective BChE inhibitors, and were used in the follow-up studies.

2.3. Kinetic Studies

As compounds 8 and 18 showed selective BChE inhibitory activity, they were selected to perform enzymatic kinetic studies with BChE in order to gain information about the mode of inhibition and binding. As shown in Figure 5, the patterns clearly indicate both compounds are mixed-type inhibitors: The presence of compounds 8 and 18 reduce the maximum velocity Vm, and increase the Km value. This means that compounds 8 and 18 can bind to the free enzyme, and to the Michaelis complex of the enzyme and substrate. The inhibition constant Ki values of 8 and 18 are shown in Table 2.

2.4. Docking Simulation of Hit Compounds

To verify the binding mode of hit compounds 8 and 18 to BChE, we carried out a docking simulation using CDOCKER module in DS 3.0 and the docking results are shown in Figure 6.
For 8 (Figure 6A), the phenyl group occupies the acyl binding pocket (mainly formed by Trp231, Leu286, and Val288) and interacts with Trp231 via π–π T-shaped interaction. In addition, the positively charged nitrogen of piperidine moiety forms a salt bridge with Asp70 in the peripheral anionic site (PAS) and π–cation interaction with Tyr332, respectively. Basic or permanently charged nitrogens are proposed as important for cation–π interactions with active site residues of ChE and blood-brain barrier (BBB) permeability, and common amidst ChE inhibitors [37,38,39,40]. The 2-methoxybenzyl group fits into the choline binding pocket (mainly formed by Tyr332 and Trp82), which further enhances the binding affinity. Furthermore, the -(CH2)2-N(Me)2 side chain points out of the gorge and stabilizes the U-shaped conformation.
For 18 (Figure 6B), although the structure is obviously different from 8, the binding pattern of 18 with BChE, including the orientation of the binding pose and the key residue for the intermolecular interaction, is similar to that of 8 with BChE. However, the inhibitory activity of 18 towards BChE is lower than 8, this could be explained by the following reasons: (1) The protonation ability of the secondary amine is lower than the tertiary amine; and (2) the 4-substituted benzyl ring of 18 is not completely inserted into acyl binding site, therefore, it is unable to interact with Trp231, which plays an important role in BChE–ligand interaction [41].

2.5. Molecular Dynamics

To further understand the binding modes of 8 and 18 with BChE, we explored the interactions between BChE and compounds by carrying out 100 ns molecular dynamics (MD) simulations for the complexes of 8-BChE and 18-BChE. Stable MD simulations trajectories were utilized for data extraction and binding free energy calculation. Time dependencies of root-mean-squared deviation (RMSD) values for the backbone atoms of proteins and compounds during the MD simulation were provided in Figure S2. The total free binding energies were calculated using the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) method. As shown in Table 3, the total binding free energy was −42.90 ± 9.39 and −48.23 ± 6.08 kcal/mol for the complexes of BChE-8 and BChE-18, respectively. The result indicated that the binding of the two compounds with BChE was energetically stable. It is noteworthy that the electrostatic energy (EEL) played a greater role in the binding of the two compounds with BChE than the hydrophobic contacts (van der Waals energy (VDWAALS)).
The contribution of potential hot residues for the binding of 8 and 18 was evaluated with the MM-PBSA method (Figure 7). Meanwhile, energy decomposition of potential hot residues was performed to evaluate which interaction was the dominating factor for the binding free energy (Figure 7). Usually, a residue is supposed to be vital for interaction of proteins with ligands if the interaction energy with ligand is lower than −1 kcal/mol [42]. For 8, none of the determined residues (Asp70, Trp82, Trp231, and Phe329) are the key residue for the binding of 8 with BChE. This result could be explained by the significant conformational change of 8 during MD process, which might be related with the presence of retained water molecules (Figure S2). For 18, Trp82 (−3.30 ± 0.35 kcal/mol), Gly116 (−1.15 ± 0.66 kcal/mol) and Phe329 (−2.61 ± 0.53 kcal/mol) were considered as important residues for binding of 18 with BChE.
The evolution of interatomic distances between functional groups of ligands and the above-mentioned key amino acid residues during MD was also performed to obtain more detailed information. As shown in Figure 8A, the distances between the functional groups of 8 and four monitored amino acid residues in the binding pocket: the piperidine ring-Asp70, the benzene ring-Trp82, the benzene ring-His231, and the benzene ring–Phe329, seemed to be changing greatly during the MD simulation, showing that these interactions are relatively unstable. In contrast, the distances between the benzene ring of 18 and important amino acid residues (Phe329 and Gly116) basically maintained the same during the MD simulation. However, the distances between compound 18 and two amino acid residues, Asp70 and Trp332, change greatly. This might explain why Asp70 and Trp332 played a minor role in the binding of 18 with BChE. All these MD results provide detailed information on the interactions between the hits and BChE.

2.6. Cell Viability Assay

The 3-(4,5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium (MTT) assay was used to evaluate potential cytotoxic effects of hit compounds 8 and 18 on neuroblastoma cell line SH-SY5Y [43]. As indicated in Figure 9, none of compounds was observed to affect cell viability at concentrations of 10 μM and 50 μM. The result indicated that 8 and 18 have preliminary safety on neuroblastoma cell line SH-SY5Y.

2.7. Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) in Silico Prediction

As shown in Table 4, compounds 8 and 18 appeared to have very poor solubility in aqueous media which have already been observed during bioassays. Both compounds were predicted to possess good absorption properties, which means both compounds could successfully enter the blood circulation from the site of administration. Additionally, compounds 8 and 18 were predicted to have high and medium blood–brain barrier (BBB) penetration, respectively, which is vital for AD treatment. The polar surface area (PSA-2D) of both compounds was less than 80. Notably, 8 may bind to CYP2D6, which would affect the efficacy of 8 and result in potential side effects. All compounds may be highly bound to plasma proteins. In this prediction, both compounds were negative in hepatotoxicity. In brief, further biological experiments are required to provide additional data.

3. Materials and Methods

3.1. Virtual Screening

The validated pharmacophore model, consisting of two hydrogen bond acceptor points, two hydrophobic points, one aromatic ring point, and one positive ionizable point, was used to screen three commercial databases using Accelrys Discovery Studio 3.0 (DS, Accelrys, Inc. San Diego, CA, USA). Virtual screening was performed based on the query fit. The drug-likeness of compounds was assessed using the DS 3.0 Filters ligands using Lipinski and Veber Rule protocol. Docking simulation was performed with the CDOCKER module implemented in DS, using the same protein applied for pharmacophore generation. A site sphere (in 10 Å radius) around the native compound was defined as the binding site. Other parameters were set as the defaults. The Cluster Ligand Module was used to cluster compounds. The number of cluster was set as 10 and the predefined set was FCFP_6. Other parameters were maintained as the defaults. Finally, 15 hits were purchased from Topscience (www.tsbiochem.com), with purity >95% (liquid chromatography-mass spectrometry, LC-MS).

3.2. Biological Evaluation

Modified Ellman’s assay was performed to measure the inhibitory effects of purchased compounds on ChEs using a Thermo Fisher Scientific spectrophotometer [44]. AChE (from human erythrocytes, C0663), BChE (from human serum, B4186), 5,5-dithiobis-(2-nitrobenzoic acid) (DTNB, D218200), acetylthiocholine iodide (ATC, A5751), and butyrylthiocholine iodide (BTC, B3253) were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Each compound was dissolved in DMSO and prepared a dilution series of six different concentrations (10−4 to 10−9 M) that the DMSO concentrations lower than 1%. For measurement, a cuvette containing 30 μL of phosphate buffer, 10 μL of AChE (2.5 units/mL) or BChE (2.5 units/mL), and 10 μL of the test compound solution was allowed to stand for 5 min before 10 μL of DTNB was added. After the addition of 20 μL of ATC or BTC, the reaction was initiated and the solution was mixed immediately. Two minutes after substrate addition, the absorption was measured at 412 nm by Thermo Fisher Scientific spectrophotometer (multiskan FC, USA). 10 μL of phosphate-buffered solution replaced the enzyme solution were used to determine the blank value. All the measurements were done in triplicate. The inhibition curve was drawn by plotting the percentage enzyme activity (100% for the reference) versus the logarithm of the test compound concentration. The IC50 values were calculated by GraphPad Prism 6 (GraphPad Software, San Diego, CA, USA), and the data were shown as mean ± standard error of mean (SEM).

3.3. Kinetic Study

Kinetic studies were performed in the same manner as the determination of ChE inhibition as previously described [45]. The substrate (BTC) was used at various concentrations (90, 150, 226, 452, and 904 μM) for each test compound concentration and the enzymatic reaction was extended to four minutes before determining the absorption. The Vmax and Km values of the Michaelis–Menten kinetics were calculated by nonlinear regression from substrate–velocity curves using Graphpad Prism 6 (GraphPad Software, San Diego, CA, USA) (Table S1). Inhibition constants were evaluated from the effect of substrate concentration (S) on the degree of inhibition according to equation:
V0 = Vmax[S]/(Km(1 + [I]/Kic) + [S] (1 + [I]/Kiu))
where S is the substrate BTC, I is the inhibitor, Kic is the enzyme–inhibitor inhibition constant of a complex formed at the catalytic site, Kiu is the Michaelis complex–inhibitor inhibition constant of a complex formed at the peripheral site, Km is the Michaelis complex, and Vm is maximal velocity.

3.4. Binding Mode Prediction

The molecular docking procedure to study the binding mode of compounds 8 and 18 against BChE was performed according to our previously reported protocol [36]. The co-crystal structure of N-((1-benzylpiperidin-3-yl)methyl)-N-(2-methoxyethyl)naphthalene-2-sulfonamide–BChE (PDB ID: 5DYW) was used here.

3.5. Molecular Dynamics

MD simulations were performed using the Particle Mesh Ewald Molecular Dynamics (PMEMD) module in AMBER 16 accelerated by GPU system consisting of the NVIDIA CUDA processor [46]. Five structural water molecules (PDB ID: 5DYW, Chain B: HOH806, HOH765, HOH749, HOH761, and HOH715) were conserved during the preparation of the protein model for molecular docking (Figure S3). The other steps were the same as usual. Briefly, the 10 top ranked binding poses were retained after molecular docking. The structure, from which its ligand had interactions with known crucial amino acids like Trp82, Trp231, and Phe329 et al., was selected for the next MD simulation. The proteins were assigned with the AMBER ff99SB force field, while the ligands were treated with the ANTECHAMBER module and the general AMBER force field [47,48]. All hydrogen atoms of proteins and ligands were added using the reduce module. The systems were solvated in a TIP3P water box in a 9 Å hexahedron. Sodium ions were added to neutralize the systems. Each system was subjected to 1000 steps of steepest-descent energy minimization followed by 1000 steps of conjugate gradient minimization with the purpose of remove possible steric stresses. Both systems were gradually heated from 0 to 300 K using a Langevin thermostat during the initial 60 ps and weak restraints 10 kcal/mol on the protein backbone atoms over 1 ns. Finally, a dynamics simulation of 20 ns NPT ensemble was set at 1 atm and 300 K. After MD simulation, the trajectory was stored every 1 ps for CPPTRAJ analysis [49]. Binding free energies and energy decomposition was determined by using the MM-PBSA method in the AMBER 16 [50]. The distance between chemical functions and residues were monitored within CPPTRAJ module. The chemical group was defined as centroid of multiple atoms. Specifically, the phenyl group was defined using the center of its six carbon atoms. The results were shown in XMGRACE.

3.6. Cell Viability Assay

The human neuroblastoma SH-SY5Y cells (3 × 103, DMEM medium supplemented with 10% fetal bovine serum, volume 0.1 mL) were placed in a 96-well flat-bottomed plate, in a humidified atmosphere of 95% air and 5% CO2 at 37 °C and grown to 80% confluence. Before cell treatment, the complete medium was replaced with reduced-serum medium (BI, 01-052-1ACS) [43]. Then, cells were treated with various concentrations of compounds for 24 h. MTT reagent (D&B, Q108115) (0.5 mg/mL) was added to the wells and the plates were incubated for 4 h at 37 °C. Supernatants were carefully removed, and 150 μL of dimethylsulfoxide (DMSO) were added into each well. The absorbance was measured at 570 nm using the spectrophotometer (Thermo, multiskan FC, USA). Cells were treated in quadruplicate. The values are the mean ± SD of three independent experiments.

3.7. ADMET in Silico Prediction

The ADMET properties (absorption, distribution, metabolism, excretion, and toxicity) and physicochemical properties of compounds 8 and 18 were calculated within ADMET software and Calculate Molecular Properties in Discovery Studio 3.0.

4. Conclusions

Considering the effectiveness of the pharmacophore model built by Sun et al. in the previous study, further pharmacophore-based virtual screening of three commercial compound collections containing 3.0 million compounds was reported in this paper to identify new BChE inhibitors [34]. Molecules passed pharmacophore filters and drug-like filters were further studied by using two complementary methods including docking and cluster to select potential hit compounds for biological evaluation. Satisfyingly, in vitro enzyme inhibition tests confirmed that compounds 8 and 18 could potently and highly selectively inhibit BChE activities (BChE IC50 < 10 μM and SI BChE > 30). In addition, the result of cell viability assay indicated that compounds 8 and 18 have preliminary safety on neuroblastoma cell line SH-SY5Y at concentrations of 10 μM and 50 μM. Furthermore, MD simulation, docking studies, and kinetic studies were performed to analyze the detailed binding modes of representative compounds 8 and 18 with BChE. The results indicated a mixed-type inhibition of the two compounds towards BChE. These active compounds with novel scaffolds provided us with a good starting point to further design potent and selective BChE inhibitors, which may benefit the treatment of AD.

Supplementary Materials

The following are available online at https://www.mdpi.com/1420-3049/24/23/4217/s1, Figure S1. The primary human ChEs IR (%) screening of hit compounds under 10 μM (A and B); The dose-dependent inhibition of hits against human BChE (C) and human AChE (D), respectively. Figure S2: Time dependencies of RMSDs for the heavy atoms (Cα, C, and N) of proteins and ligands. Figure S3. Schematic representation of interaction between ligands (compounds 8 and 18) and BChE before the MD simulation. Figure S4. The binding modes of ligands (compounds 8 and 18) with BChE (PDB ID: 5DYW). Protein models were taken from the last frame of the MD simulation. Table S1. The apparent Vmax and Km values for compounds 8 and 18 in kinetic studies.

Author Contributions

Y.Z., L.X. and Q.X. conceived and designed the experiments; Y.Z., H.Y., F.W., J.L., Z.T. and X.C. performed the experiments; Y.C. and Y.Y. analyzed the data; Y.Z. and X.L. wrote the paper. All the Authors accepted the final version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 81803367; Basic Research and Frontier Exploration Project of Chongqing, grant number cstc2018jcyjAX0715; Fundamental Research Funds for the Central Universities, grant number XDJK2019C103.

Acknowledgments

We thank Haopeng Sun from China Pharmaceutical University for providing us with the pharmacophore model used in this study and technical support. We thank Kwon Ho Hong from University of Minnesota for helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Lane, C.A.; Hardy, J.; Schott, J.M. Alzheimer’s disease. Eur. J. Neurol. 2018, 25, 59–70. [Google Scholar] [CrossRef]
  2. Wang, Z.; Hu, J.; Yang, X.; Feng, X.; Li, X.; Huang, L.; Chan, A.S.C. Design, Synthesis, and Evaluation of Orally Bioavailable Quinoline−Indole Derivatives as Innovative Multitarget-Directed Ligands: Promotion of Cell Proliferation in the Adult Murine Hippocampus for the Treatment of Alzheimer’s Disease. J. Med. Chem. 2018, 61, 1871–1894. [Google Scholar] [CrossRef]
  3. Alzheimer’s Disease International. Available online: https://www.alz.co.uk/research/world-report-2018 (accessed on 19 November 2019).
  4. Więckowska, A.; Wichur, T.; Godyń, J.; Bucki, A.; Marcinkowska, M.; Siwek, A.; Więckowski, K.; Zaręba, P.; Knez, D.; Głuch-Lutwin, M.; et al. Novel Multitarget-Directed Ligands Aiming at Symptoms and Causes of Alzheimer’ s Disease. ACS Chem. Neurosci. 2018, 9, 1195–1214. [Google Scholar] [CrossRef]
  5. Wu, Y.; Li, Z.; Huang, Y.; Wu, D.; Luo, H. Novel Phosphodiesterase Inhibitors for Cognitive Improvement in Alzheimer ’ s Disease. J. Med. Chem. 2018, 61, 5467–5483. [Google Scholar] [CrossRef] [PubMed]
  6. Cline, E.N.; Bicca, M.A.; Viola, K.L.; Klein, W.L. The Amyloid-β Oligomer Hypothesis: Beginning of the Third Decade. J. Alzheimer’s Dis. 2018, 64, S567–S610. [Google Scholar] [CrossRef] [PubMed]
  7. Caruso, G.; Distefano, D.A.; Parlascino, P.; Fresta, C.G.; Lazzarino, G.; Lunte, S.M.; Nicoletti, V.G. Receptor-mediated toxicity of human amylin fragment aggregated by short- and long-term incubations with copper ions. Mol. Cell. Biochem. 2017, 425, 85–93. [Google Scholar] [CrossRef] [PubMed]
  8. Julien, C.; Tomberlin, C.; Roberts, C.M.; Akram, A.; Stein, G.H.; Silverman, M.A.; Link, C.D. In vivo induction of membrane damage by β-amyloid peptide oligomers. Acta Neuropathol. Commun. 2018, 6, 131. [Google Scholar] [CrossRef] [PubMed]
  9. Du, X.; Wang, X.; Geng, M. Alzheimer’s disease hypothesis and related therapies. Transl. Neurodegener. 2018, 7, 2. [Google Scholar] [CrossRef]
  10. Sanabria-Castro, A.; Alvarado-Echeverra, I.; Monge-Bonilla, C. Molecular Pathogenesis of Alzheimer’s Disease: An Update. Ann. Neurosci. 2017, 24, 46–54. [Google Scholar] [CrossRef]
  11. Li, Q.; He, S.Y.; Chen, Y.; Feng, F.; Qu, W.; Sun, H.P. Donepezil-based Multi-functional Cholinesterase Inhibitors for Treatment of Alzheimer’s Disease. Eur. J. Med. Chem. 2018, 158, 463–477. [Google Scholar] [CrossRef]
  12. Contestabile, A. The history of the cholinergic hypothesis. Behav. Brain Res. 2011, 221, 334–340. [Google Scholar] [CrossRef] [PubMed]
  13. Hu, J.; Zhang, M.; Yu, L.B.; Ju, Y. Bioorganic & Medicinal Chemistry Letters Synthesis and binding ability of 1, 2, 3-triazole-based triterpenoid receptors for recognition of Hg2+ ion. Bioorg. Med. Chem. Lett. 2010, 20, 4342–4345. [Google Scholar] [PubMed]
  14. Csuk, R.; Albert, S.; Kluge, R.; StrÖhl, D. Resveratrol Derived Butyrylcholinesterase Inhibitors. Arch. Pharm. Chem. Life Sci. 2013, 346, 499–503. [Google Scholar] [CrossRef] [PubMed]
  15. Mishra, P.; Kumar, A.; Panda, G. Anti-cholinesterase hybrids as multi-target-directed ligands against Alzheimer’s disease (1998–2018). Bioorg. Med. Chem. 2019, 27, 895–930. [Google Scholar] [CrossRef] [PubMed]
  16. Dighe, S.N.; Deora, G.S.; Mora, E.D.; Nachon, F.; Chan, S.; Parat, M.; Brazzolotto, X.; Ross, B.P. Discovery and Structure−Activity Relationships of a Highly Selective Butyrylcholinesterase Inhibitor by Structure-Based Virtual Screening. J. Med. Chem. 2016, 59, 7683–7689. [Google Scholar] [CrossRef]
  17. Arendt, T.; Bigl, V.; Walther, F.; Sonntag, M. Decreased Ratio of CSF Acetylcholinesterase to Butyrylcholinesterase Activity in Alzheimer’s Disease. Lancet 1984, 323, 173. [Google Scholar] [CrossRef]
  18. Li, Q.; Yang, H.Y.; Chen, Y.; Sun, H.P. Recent progress in the identification of selective butyrylcholinesterase inhibitors for Alzheimer’s disease. Eur. J. Med. Chem. 2017, 132, 294–309. [Google Scholar] [CrossRef]
  19. Mesulam, M.M.; Guillozet, A.; Shaw, P.; Levey, A.; Duysen, E.G.; Lockridge, O. Acetylcholinesterase knockouts establish central cholinergic pathways and can use butyrylcholinesterase to hydrolyze acetylcholine. Neuroscience 2002, 110, 627–639. [Google Scholar] [CrossRef]
  20. Greig, N.H.; Utsuki, T.; Ingram, D.K.; Wang, Y.; Pepeu, G.; Scali, C.; Yu, Q.; Mamczarz, J.; Holloway, H.W.; Giordano, T.; et al. Selective butyrylcholinesterase inhibition elevates brain acetylcholine, augments learning and lowers Alzheimer β-amyloid peptide in rodent. Proc. Natl. Acad. Sci. USA 2005, 102, 17213–17218. [Google Scholar] [CrossRef]
  21. Fang, L.; Pan, Y.; Muzyka, J.L.; Zhan, C.G. Active site gating and substrate specificity of butyrylcholinesterase and acetylcholinesterase: Insights from molecular dynamics simulations. J. Phys. Chem. B 2011, 115, 8797–8805. [Google Scholar] [CrossRef]
  22. Košak, U.; Brus, B.; Knez, D.; Šink, R.; Žakel, S.; Trontelj, J.; Pišlar, A.; Šlenc, J.; Gobec, M.; Živin, M.; et al. Development of an in-vivo active reversible butyrylcholinesterase inhibitor. Sci. Rep. 2016, 6, 39495. [Google Scholar] [CrossRef] [PubMed]
  23. Košak, U.; Brus, B.; Knez, D.; Žakel, S.; Trontelj, J.; Pišlar, A.; Šink, R.; Jukič, M.; Živin, M.; Podkowa, A.; et al. The Magic of Crystal Structure-Based Inhibitor Optimization: Development of a Butyrylcholinesterase Inhibitor with Picomolar Affinity and in Vivo Activity. J. Med. Chem. 2018, 61, 119–139. [Google Scholar] [CrossRef] [PubMed]
  24. De Candia, M.; Zaetta, G.; Denora, N.; Tricarico, D.; Majellaro, M.; Cellamare, S.; Altomare, C.D. New azepino [4,3-b] indole derivatives as nanomolar selective inhibitors of human butyrylcholinesterase showing protective effects against NMDA-induced neurotoxicity. Eur. J. Med. Chem. 2017, 125, 288–298. [Google Scholar] [CrossRef] [PubMed]
  25. Sawatzky, E.; Wehle, S.; Kling, B.; Wendrich, J.; Bringmann, G.; Sotri, C.A.; Heilmann, J.; Decker, M. Discovery of Highly Selective and Nanomolar Carbamate-Based Butyrylcholinesterase Inhibitors by Rational Investigation into Their Inhibition Mode. J. Med. Chem. 2016, 59, 2067–2082. [Google Scholar] [CrossRef] [PubMed]
  26. Košak, U.; Kenz, D.; Coquelle, N.; Brus, B.; Pišlar, A.; Nachon, F.; Brazzolotto, X.; Kos, J.; Colletier, J.P.; Gobec, S. N-Propargylpiperidines with naphthalene-2-carboxamide or naphthalene-2-sulfonamide moieties: Potential multifunctional anti-Alzheimer’s agents. Bioorg. Med. Chem. 2017, 25, 633–645. [Google Scholar] [CrossRef] [PubMed]
  27. Brus, B.; Košak, U.; Turk, S.; Pišlar, A.; Coquelle, N.; Kos, J.; Stojan, J.; Colletier, J.P.; Gobec, S. Discovery, biological evaluation, and crystal structure of a novel nanomolar selective butyrylcholinesterase inhibitor. J. Med. Chem. 2014, 57, 8167–8179. [Google Scholar] [CrossRef] [PubMed]
  28. Meden, A.; Knez, D.; Jukič, M.; Brazzolotto, X.; Gršič, M.; Pišlar, A.; Zahirović, A.; Kos, J.; Nachon, F.; Svete, J.; et al. Tryptophan-derived butyrylcholinesterase inhibitors as promising leads against Alzheimer’s disease. Chem. Commun. 2019, 55, 3765–3768. [Google Scholar] [CrossRef]
  29. Jiang, C.S.; Ge, Y.X.; Cheng, Z.Q.; Wang, Y.Y.; Tao, H.R.; Zhu, K.K.; Zhang, H. Discovery of New Selective Butyrylcholinesterase (BChE) Inhibitors with Anti-Aβ Aggregation Activity: Structure-Based Virtual Screening, Hit Optimization and Biological Evaluation. Molecules 2019, 24, 2568. [Google Scholar] [CrossRef]
  30. Atatreh, N.; Al Rawashdah, S.; Al Neyadi, S.S.; Abuhamdah, S.M.; Ghattas, M.A. Discovery of new butyrylcholinesterase inhibitors via structure-based virtual screening. J. Enzym. Inhib. Med. Chem. 2019, 34, 1373–1379. [Google Scholar] [CrossRef]
  31. Beesu, M.; Caruso, G.; Salyer, A.C.D.; Khetani, K.K.; Sil, D.; Weerasinghe, M.; Tanji, H.; Ohto, U.; Shimizu, T.; David, S.A. Structure-Based Design of Human TLR8-Specific Agonists with Augmented Potency and Adjuvanticity. J. Med. Chem. 2015, 58, 7833–7849. [Google Scholar] [CrossRef]
  32. Neves, B.J.; Braga, R.C.; Melo-Filho, C.C.; Moreira-Filho, J.T.; Muratov, E.N.; Andrade, C.H. QSAR-based virtual screening: Advances and applications in drug discovery. Front. Pharmacol. 2018, 9, 1275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Beesu, M.; Caruso, G.; Salyer, A.C.D.; Shukla, N.M.; Khetani, K.K.; Smith, L.J.; Fox, L.M.; Tanji, H.; Ohto, U.; Shimizu, T.; et al. Identification of a Human Toll-Like Receptor (TLR) 8-Specific Agonist and a Functional Pan-TLR Inhibitor in 2-Aminoimidazoles. J. Med. Chem. 2016, 59, 3311–3330. [Google Scholar] [CrossRef] [PubMed]
  34. Lu, X.; Yang, H.; Li, Q.; Chen, Y.; Li, Q.; Zhou, Y.; Feng, F.; Liu, W.; Guo, Q.; Sun, H. Expansion of the scaffold diversity for the development of highly selective butyrylcholinesterase (BChE) inhibitors: Discovery of new hits through the pharmacophore model generation, virtual screening and molecular dynamics simulation. Bioorg. Chem. 2019, 85, 117–127. [Google Scholar] [CrossRef] [PubMed]
  35. Lipinski, C.A. Drug-like properties and the causes of poor solubility and poor permeability. J. Pharmacol. Toxicol. Methods 2000, 44, 235–249. [Google Scholar] [CrossRef]
  36. Chen, Y.; Liu, Z.L.; Fu, T.M.; Li, W.; Xu, X.L.; Sun, H.P. Discovery of new acetylcholinesterase inhibitors with small core structures through shape-based virtual screening. Bioorg. Med. Chem. Lett. 2015, 25, 3442–3446. [Google Scholar] [CrossRef]
  37. Rankovic, Z. CNS drug design: Balancing physicochemical properties for optimal brain exposure. J. Med. Chem. 2015, 58, 2584–2608. [Google Scholar] [CrossRef]
  38. Ghose, A.K.; Herbertz, T.; Hudkins, R.L.; Dorsey, B.D.; Mallamo, J.P. Knowledge-based, central nervous system (CNS) lead selection and lead optimization for CNS drug discovery. ACS Chem. Neurosci. 2012, 3, 50–68. [Google Scholar] [CrossRef] [Green Version]
  39. Belinskaya, D.; Juffer, A.; Shestakova, N. The role of electrostatic interactions in the absorption of ligands to the active sites of cholinesterases, as indicated by molecular modeling data. Russ. J. Bioorg. Chem. 2010, 36, 185–190. [Google Scholar] [CrossRef]
  40. Nachon, F.; Carletti, E.; Ronco, C.; Trovaslet, M.; Nicolet, Y.; Jean, L.; Renard, P.-Y. Crystal structures of human cholinesterases in complex with huprine W and tacrine: Elements of specificity for anti- Alzheimer’s drugs targeting acetyl-and butyryl-cholinesterase. Biochem. J. 2013, 453, 393–399. [Google Scholar] [CrossRef] [Green Version]
  41. Rosenberry, T.L.; Id, X.B.; Macdonald, I.R.; Wandhammer, M.; Trovaslet-leroy, M.; Darvesh, S.; Nachon, F. Comparison of the Binding of Reversible Inhibitors to Human Butyrylcholinesterase and Acetylcholinesterase: A Crystallographic, Kinetic and Calorimetric Study. Molecules 2017, 22, 2098. [Google Scholar] [CrossRef] [Green Version]
  42. Zhang, C.; Feng, L.J.; Huang, Y.; Wu, D.; Li, Z.; Zhou, Q.; Wu, Y.; Luo, H.B. Discovery of novel phosphodiesterase-2A inhibitors by structure-based virtual screening structural optimization, and bioassay. J. Chem. Inf. Model. 2017, 57, 355–364. [Google Scholar] [CrossRef] [PubMed]
  43. Sang, Z.; Pan, W.; Wang, K.; Ma, Q.; Yu, L.; Yang, Y.; Bai, P.; Leng, C.; Xu, Q.; Li, X.; et al. European Journal of Medicinal Chemistry Design, synthesis and evaluation of novel ferulic acid- O -alkylamine derivatives as potential multifunctional agents for the treatment of Alzheimer’s disease. Eur. J. Med. Chem. 2017, 130, 379–392. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, Y.; Lin, H.; Yang, H.; Tan, R.; Bian, Y.; Fu, T.; Li, W.; Wu, L.; Pei, Y.; Sun, H. Discovery of new acetylcholinesterase and butyrylcholinesterase inhibitors through structure-based virtual screening. RSC Adv. 2017, 7, 3429–3438. [Google Scholar] [CrossRef] [Green Version]
  45. Lu, X.; He, S.; Li, Q.; Yang, H.; Jiang, X.; Lin, H.; Chen, Y.; Qu, W.; Feng, F.; Bian, Y.; et al. Investigation of multi-target-directed ligands (MTDLs) with (IDO1) inhibition: The design, synthesis of miconazole analogues targeting Alzheimer’s disease. Bioorg. Med. Chem. 2018, 26, 1665–1674. [Google Scholar] [CrossRef]
  46. Case, D.A.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; Homeyer, N.; et al. AMBER 2016; University of California: San Francisco, CA, USA, 2016. [Google Scholar]
  47. Zhang, W.; Hou, T.; Qiao, X.; Xu, X. Parameters for the Generalized Born Model Consistent with RESP Atomic Partial Charge Assignment Protocol. J. Phys. Chem. B 2003, 107, 9071–9078. [Google Scholar] [CrossRef]
  48. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of Multiple Amber Force Fields and Development of Improved Protein Backbone Parameters. Proteins Struct. Funct. Bioinform. 2006, 65, 712–725. [Google Scholar] [CrossRef] [Green Version]
  49. Roe, D.R.; Cheatham, T.E. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef]
  50. Miller, B.R.; Mcgee, T.D.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
Sample Availability: Samples of compounds 5–19 are available from the authors.
Molecules 24 04217 g001 550
Figure 1. Recently reported selective human butyrylcholinesterase (hBChE) inhibitors.
Figure 1. Recently reported selective human butyrylcholinesterase (hBChE) inhibitors.
Molecules 24 04217 g001
Molecules 24 04217 g002 550
Figure 2. An overview of the structure-based pharmacophore virtual screening (SBP VS) protocol applied to identify selective butyrylcholinesterase (BChE) inhibitors.
Figure 2. An overview of the structure-based pharmacophore virtual screening (SBP VS) protocol applied to identify selective butyrylcholinesterase (BChE) inhibitors.
Molecules 24 04217 g002
Molecules 24 04217 g003 550
Figure 3. The pharmacophore model bound with 1-BChE complex. The Compound is represented in yellow stick mode. In the pharmacophore model, the hydrogen bond acceptor, hydrophobic portion, aromatic ring, cation, and excluded volumes are colored green, cyan, orange, red, and grey, respectively.
Figure 3. The pharmacophore model bound with 1-BChE complex. The Compound is represented in yellow stick mode. In the pharmacophore model, the hydrogen bond acceptor, hydrophobic portion, aromatic ring, cation, and excluded volumes are colored green, cyan, orange, red, and grey, respectively.
Molecules 24 04217 g003
Molecules 24 04217 g004 550
Figure 4. Chemical structure of the compounds purchased for in vitro tests.
Figure 4. Chemical structure of the compounds purchased for in vitro tests.
Molecules 24 04217 g004
Molecules 24 04217 g005 550
Figure 5. Representative plot of BChE activity and the effect of substrate concentration (90–904 μM) in the absence of inhibitor and in the presence of 8 and 18 (0.5–2 μM). (A) Substrate-velocity curves of BChE inhibition by compound 8; (B) Substrate-velocity curves of BChE inhibition by compound 18.
Figure 5. Representative plot of BChE activity and the effect of substrate concentration (90–904 μM) in the absence of inhibitor and in the presence of 8 and 18 (0.5–2 μM). (A) Substrate-velocity curves of BChE inhibition by compound 8; (B) Substrate-velocity curves of BChE inhibition by compound 18.
Molecules 24 04217 g005
Molecules 24 04217 g006 550
Figure 6. Binding mode predictions for compound 8 (A) and 18 (B) with BChE domain (PDB ID: 5DYW). Compounds were shown in green stick mode; key residues were shown in yellow stick mode; the solvent accessible surface (SAS) of proteins contoured by white.
Figure 6. Binding mode predictions for compound 8 (A) and 18 (B) with BChE domain (PDB ID: 5DYW). Compounds were shown in green stick mode; key residues were shown in yellow stick mode; the solvent accessible surface (SAS) of proteins contoured by white.
Molecules 24 04217 g006
Molecules 24 04217 g007 550
Figure 7. Residue contributions and energy decomposition of potential hot residues of 8 (A) and 18 (B); for all energies, unit is kcal/mol.
Figure 7. Residue contributions and energy decomposition of potential hot residues of 8 (A) and 18 (B); for all energies, unit is kcal/mol.
Molecules 24 04217 g007
Molecules 24 04217 g008 550
Figure 8. The evolution of interatomic distances between chemical functions of the ligands and key amino acid residues in the binding pocket during molecular dynamics (MD). (A) The evolution of interatomic distances between chemical functions of 8 and key amino acid residues during MD. Asp70 represented the distance between N atom of the piperidine ring and Asp70; Trp82 represented the distance between the methoxy-substituted benzene ring and Trp82; Trp231 represented the distance between the unsubstituted benzene ring and Trp231; (B) The evolution of interatomic distances between chemical functions of 18 and key amino acid residues during MD. Asp70 represented the distance between N atom of the secondary amine group and Asp70; Gly116 represented the distance between the methoxy-substituted benzene ring and Gly116; Phe329 represented the distance between the methoxy-substituted benzene ring and Phe329; Tyr332 represented the distance between N atom of the secondary amine group and Tyr332.
Figure 8. The evolution of interatomic distances between chemical functions of the ligands and key amino acid residues in the binding pocket during molecular dynamics (MD). (A) The evolution of interatomic distances between chemical functions of 8 and key amino acid residues during MD. Asp70 represented the distance between N atom of the piperidine ring and Asp70; Trp82 represented the distance between the methoxy-substituted benzene ring and Trp82; Trp231 represented the distance between the unsubstituted benzene ring and Trp231; (B) The evolution of interatomic distances between chemical functions of 18 and key amino acid residues during MD. Asp70 represented the distance between N atom of the secondary amine group and Asp70; Gly116 represented the distance between the methoxy-substituted benzene ring and Gly116; Phe329 represented the distance between the methoxy-substituted benzene ring and Phe329; Tyr332 represented the distance between N atom of the secondary amine group and Tyr332.
Molecules 24 04217 g008
Molecules 24 04217 g009 550
Figure 9. The cytotoxicity of hit compounds on SH-SY5Y cells.
Figure 9. The cytotoxicity of hit compounds on SH-SY5Y cells.
Molecules 24 04217 g009
Table
Table 1. The inhibitory activities against cholinesterases (ChEs) of the hits from virtual screening.
Table 1. The inhibitory activities against cholinesterases (ChEs) of the hits from virtual screening.
Compound BChEAChE
IR a (%)IC50 b (μM) IR c (%)IC50 (μM)
57.2 ± 0.6nd. d−0.31 ± 0.5nd.
68.5 ± 0.3nd.−1.5 ± 0.5nd.
716.3 ± 1.1nd.0.6 ± 0.6nd.
868.6 ± 0.71.1 ± 0.658.5 ± 1.243.2 ± 17.6
915.5 ± 1.6nd.16.0 ± 1.5nd.
109.9 ± 1.0nd.7.8 ± 0.7nd.
1114.8 ± 1.3nd.−0.7 ± 0.7nd.
12−1.8 ± 1.1nd.1.1 ± 1.0nd.
1320.1 ± 1.2nd.11.3 ± 1.3nd.
143.4 ± 0.4nd.10.9 ± 0.8nd.
15−0.6 ± 0.5nd.0.6 ± 1.0nd.
1626.4 ± 1.1nd.38.7 ± 1.7nd.
1711.8 ± 1.2nd.2.9 ± 0.5nd.
1858.4 ± 0.96.3 ± 2.011.1 ± 1.5nd.
19
Tacrine		61.2 ± 1.8
100 		2.4 ± 1.0
0.003 ± 0.004		53.2 ± 0.6
95.2 ± 0.3		13.8 ± 6.0
0.01 ± 0.003
All data are shown as mean ± SEM of three experiments. SEM = standard error of mean. a Inhibition ratio (IR) against AChE at 10 μM. b IC50 values represent the concentration of inhibitor required to decrease enzyme activity by 50%. c Inhibition ratio (IR) against BChE at 10 μM. d nd = not determined.
Table
Table 2. The inhibition constants for the inhibition of BChE by compounds 8 and 18.
Table 2. The inhibition constants for the inhibition of BChE by compounds 8 and 18.
CompoundKic aKiu b
80.88 ± 0.07 µM3.61 ± 0.24 µM
180.93 ± 0.13 µM2.31 ± 0.32 µM
All data are shown as mean ± SEM of three experiments. a Kic is the inhibition constant for the competitive part of inhibition. b Kiu is the inhibition constant for the uncompetitive part of inhibition.
Table
Table 3. Predicted binding free energies (kcal/mol) for bindings of 8 or 18 with BChE by the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) method.
Table 3. Predicted binding free energies (kcal/mol) for bindings of 8 or 18 with BChE by the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) method.
Energy Terms (kcal/mol)818
VDWAALS a−43.3 ± 5.8 −52.2 ± 3.5
EEL b−51.3 ± 16.8−142 ± 24.0
EGB c57.6 ± 15.0153 ± 21.3
ESURF d−5.9 ± 1.0−7.1 ± 0.7
DELTA G gas e−94.6 ± 19.1−195 ± 24.6
DELTA G solv f51.7 ± 14.2146 ± 20.8
DELTA TOTAL g−42.9 ± 9.4−48.2 ± 6.1
a van der Waals energy. b Electrostatic energy. c Polar solvation energy. d Non-polar solvation energy. e Total gas phase free energy. f Total solvation free energy. g Total binding free energy.
Table
Table 4. Predicted pharmacokinetic properties of hits.
Table 4. Predicted pharmacokinetic properties of hits.
Comp.AlogP98 aPSA-2D bSolubility
Level c		Absorption
Level d		BBB
Level e		PPB fCYP2D6 gHepatotoxic h
83.5353.589301truetruefalse
183.51463.209302truefalsefalse
a AlogP98: Lipophilicity descriptor. b PSA-2D: Polar surface area. c Solubility Level: (0, good; 1, moderate; 2, poor; 3, very poor). d Absorption Level: (0, good; 1, moderate; 2, poor; 3, very poor). e blood–brain barrier (BBB) Level: (0, very high blood–brain barrier penetration; 1, high; 2, medium; 3, low). f PPB Prediction: PPB refers to plasma protein binding. The classification describing whether a compound is highly bound (>= 90% bound) to plasma proteins using the cutoff Bayesian score of −2.209 (obtained by minimizing the total number of false positives and false negatives). g CYP2D6 Prediction: The classification describing whether a compound is a cytochrome P450 2D6 (CYP2D6) inhibitor using the cutoff Bayesian score of 0.161 (obtained by minimizing the total number of false positives and false negatives). h Hepatotoxic Prediction: The classification describing whether a compound is hepatotoxic using the cutoff Bayesian score of −4.154 (obtained by minimizing the total number of false positives and false negatives).

Share and Cite

MDPI and ACS Style

Zhou, Y.; Lu, X.; Yang, H.; Chen, Y.; Wang, F.; Li, J.; Tang, Z.; Cheng, X.; Yang, Y.; Xu, L.; et al. Discovery of Selective Butyrylcholinesterase (BChE) Inhibitors through a Combination of Computational Studies and Biological Evaluations. Molecules 2019, 24, 4217. https://doi.org/10.3390/molecules24234217

AMA Style

Zhou Y, Lu X, Yang H, Chen Y, Wang F, Li J, Tang Z, Cheng X, Yang Y, Xu L, et al. Discovery of Selective Butyrylcholinesterase (BChE) Inhibitors through a Combination of Computational Studies and Biological Evaluations. Molecules. 2019; 24(23):4217. https://doi.org/10.3390/molecules24234217

Chicago/Turabian Style

Zhou, You, Xin Lu, Hongyu Yang, Yao Chen, Feng Wang, Jifu Li, Zhiran Tang, Xifei Cheng, Yingbin Yang, Li Xu, and et al. 2019. "Discovery of Selective Butyrylcholinesterase (BChE) Inhibitors through a Combination of Computational Studies and Biological Evaluations" Molecules 24, no. 23: 4217. https://doi.org/10.3390/molecules24234217

APA Style

Zhou, Y., Lu, X., Yang, H., Chen, Y., Wang, F., Li, J., Tang, Z., Cheng, X., Yang, Y., Xu, L., & Xia, Q. (2019). Discovery of Selective Butyrylcholinesterase (BChE) Inhibitors through a Combination of Computational Studies and Biological Evaluations. Molecules, 24(23), 4217. https://doi.org/10.3390/molecules24234217

Article Metrics

Back to TopTop