Next Article in Journal
Design, Synthesis, and Biological Evaluation of 2-Mercaptobenzoxazole Derivatives as Potential Multi-Kinase Inhibitors
Next Article in Special Issue
Piperine–Chlorogenic Acid Hybrid Inhibits the Proliferation of the SK-MEL-147 Melanoma Cells by Modulating Mitotic Kinases
Previous Article in Journal
The Implication of Topoisomerase II Inhibitors in Synthetic Lethality for Cancer Therapy
Previous Article in Special Issue
2-Arylpropionic Acid Pyrazolamides as Cannabinoid CB2 Receptor Inverse Agonists Endowed with Anti-Inflammatory Properties
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pharmaceuticals
Volume 16
Issue 1
10.3390/ph16010095
pharmaceuticals-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

Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies

by
Raquel B. M. de Almeida
1,
Deyse B. Barbosa
1,
Mayra R. do Bomfim
1,
Jéssika A. O. Amparo
2,
Bruno S. Andrade
3,
Silvia L. Costa
2,
Joaquín M. Campos
4,5,
Jorddy N. Cruz
6,
Cleydson B. R. Santos
5,6,*,
Franco H. A. Leite
1 and
Mariana B. Botura
1
1
Department of Health, State University of Feira de Santana, Feira de Santana 44036-900, BA, Brazil
2
Institute of Health Sciences, Federal University of Bahia, Salvador 40170-110, BA, Brazil
3
Department of Biological Sciences, State University of Southwest of Bahia, Jequié 45208-091, BA, Brazil
4
Biosanitary Institute of Granada (ibs.GRANADA), University of Granada, 18071 Granada, Spain
5
Department of Pharmaceutical and Organic Chemistry, Faculty of Pharmacy, Campus of Cartuja, University of Granada, 18071 Granada, Spain
6
Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapá 68902-280, AP, Brazil
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2023, 16(1), 95; https://doi.org/10.3390/ph16010095
Submission received: 5 December 2022 / Revised: 31 December 2022 / Accepted: 3 January 2023 / Published: 9 January 2023
(This article belongs to the Special Issue 10th Brazilian Symposium on Medicinal Chemistry (BrazMedChem_2022))
Browse Figures
Versions Notes

Abstract

:
The enhancement of cholinergic functions via acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) inhibition is considered a valuable therapeutic strategy for the treatment of Alzheimer’s disease. This study aimed to evaluate the in vitro effect of ZINC390718, previously filtered using computational approaches, on both cholinesterases and to characterize, using a molecular dynamics (MD) simulation, the possible binding mode of this compound inside the cholinesterase enzymes. The in vitro cytotoxicity effect was also investigated using a primary astrocyte-enriched glial cell culture. ZINC390718 presented in vitro dual inhibitory activity against AChE at a high micromolar range (IC50 = 543.8 µM) and against BuChE (IC50 = 241.1 µM) in a concentration-dependent manner, with greater activity against BuChE. The MD simulation revealed that ZINC390718 performed important hydrophobic and H-bond interactions with the catalytic residue sites on both targets. The residues that promoted the hydrophobic interactions and H-bonding in the AChE target were Leu67, Trp86, Phe123, Tyr124, Ser293, Phe295, and Tyr341, and on the BuChE target, they were Asp70, Tyr332, Tyr128, Ile442, Trp82, and Glu197. The cytotoxic effect of Z390718, evaluated via cell viability, showed that the molecule has low in vitro toxicity. The in vitro and in silico results indicate that ZINC390718 can be used as chemotype for the optimization and identification of new dual cholinesterase inhibitors.

1. Introduction

Alzheimer’s disease (AD) is a neurodegenerative and irreversible disorder, characterized by a progressive decrease in memory and other cognitive functions that mainly affect the patient’s social life [1]. AD is the most common cause of dementia, especially in the elderly, who account for about 60% to 70% of cases [2]. AD has a complex and multifactorial etiology that involves low levels of acetylcholine, deposits of amyloid plaques, and neurofibrillary tangles of aggregated tau proteins in the brain [3].
Cognitive and behavioral disorders in AD are related to the rupture of cholinergic pathways in the cerebral cortex and in the basal forebrain; they are also associated with decreased levels of the neurotransmitter acetylcholine (Ach) [4]. In this context, studies on AD have shown abnormalities in many neurotransmission systems, and the most prominent and serious damage is that of the cholinergic system, with a selective loss of presynaptic cholinergic neurons that project in the cerebral cortex and hippocampus, leading to the cholinergic hypothesis of the development of AD [5].
Thus, the main approach to treat the disease is to restore the level of Ach through the inhibition of cholinesterase [6]. Cholinesterases are a group of esterase enzymes, and they are responsible for catalyzing the hydrolysis of Ach in choline and acetic acid and are divided into two types: acetylcholinesterase (AChE, EC 3.1.1.7) and butyrylcholinesterase (BuChE, EC 3.1.1.8) [7]. They are homologous glycoproteins and share 50% of sequential identity, but they differ in their distribution in tissues, in their kinetic properties, and in the specificity of their substrates [8].
AChE is found abundantly in the central nervous system (CNS) and is expressed by neurons in skeletal muscles and the erythrocyte membrane, while BuChE is associated with glial cells and is found mainly in blood plasma [8]. Cholinesterase inhibitors, recognized as an official treatment of Alzheimer’s dementia, are also prescribed for dementia without Alzheimer’s; psychiatric disorders; and neurological diseases, such as Huntington’s disease (HD), cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), frontotemporal dementia (FTD), dementia in multiple sclerosis (MS), and progressive supranuclear palsy (PSP) [9].
During the development of AD, AChE levels in the CNS decrease, while BuChE activity increases progressively and, consequently, mainly BuChE controls Ach regulation. In addition, studies employing an AD mouse model have suggested that BuChE may play a role in the etiology of AD due to its contribution to the deposition of amyloid plaques [10]. Thus, an agent that inhibits both AChE and BuChE (a dual inhibitor), such as rivastigmine, provides greater benefits to patients with dementia by decreasing the compensatory effect of Ach hydrolysis by BuChE when AChE levels decrease [11].
Currently, the FDA-approved anticholinesterase drugs are donepezil and galantamine (AChE inhibitors) and rivastigmine (dual inhibitor), which have beneficial effects on the cognitive, functional, and behavioral symptoms of the disease [12]. These drugs are generally well-tolerated, but symptoms such as diarrhea, nausea, vomiting, muscle cramps, decreased blood pressure, insomnia, fatigue, and a loss of appetite have been reported [13].
In this context, much effort has been directed towards the development of new dual cholinesterase inhibitors for the treatment of AD neurological disorders [14]. A combination of in silico and in vitro studies represents a valuable strategy to find new candidate drugs targeting different enzymes related to AD [15].
The toxicity evaluation of bioactive molecules is an important step in the development of drugs. In vitro cytotoxicity assays with cultured cells are widely used for preliminary toxicological studies of chemicals and for drug screening [16,17,18]. For substances with action on the central nervous system, evaluations of cytotoxic effects on neural and glial cells are performed. Astrocyte cultures can be used to evaluate compounds for neurotoxicity or the ability of these compounds to protect astrocytes from disease [19].
Thus, the aim of this study was to identify a new molecule with potential dual inhibition against acetylcholinesterase and butyrylcholinesterase using in silico tools and evaluations of enzymatic inhibition and preliminary in vitro toxicity.

2. Results and Discussion

2.1. In Vitro Anticholinesterase Assays

2.1.1. Kinetic Parameters

For enzymatic assays, it must be considered that the reactions of the enzymes also depend on factors such as the pH of the substrate and enzyme. The optimal pH of many enzymes is within the physiological range, often between pH 7 and 8 [20]. The normal arterial pH of humans is 7.4, and the blood pH of some other species differs slightly from pH 7.4 [21]. It has been observed that the maximum cholinesterase catalytic range occurs at pH 7.4.
The concentration of the substrate used in kinetic tests directly influences the chance of identifying inhibitors [22]. The Km and Vmax values of the substrate used in the enzymatic reactions with the cholinesterases were determined, as shown in Figure 1.
The concentration of the DTNB followed the value described by Ellman et al. [23] (0.01 M). Even when carrying out tests with varying concentrations, it was observed that this concentration could be maintained without interfering with the reaction of the enzymes. The reading times of the enzymatic reactions with AChE (10 min) and BuChE (20 min) were determined after identifying the decline in the line that represents the end of the total substrate consumption rate.

2.1.2. Dual Anticholinesterase Activity of ZINC390718

In an attempt to assess the ability of ZINC390718 to inhibit the two enzymes, the IC50 values of both targets were calculated (Figure 2). After kinetic parametrization, the enzymatic inhibition activities of the compound ZINC390718 against AChE and BuChE showed a significant effect (p < 0.05) in a concentration-dependent manner for both enzymes (Figure 2). The IC50 concentrations found were 543.8 and 241.1 µM for AChE and BuChE, respectively.
ZINC390718 showed a greater selectivity for BuChE than AChE (selectivity index = 2.25), although for both targets, high IC50 concentration values were observed. In AChE, the estimated volume of the active site is relatively small (302 Å3), being coated with fourteen aromatic amino acid residues, while that of BuChE, with only eight aryl residues, is considerably higher (502 Å3), and this difference in volume can influence the determination of the size of the inhibitor molecules and their accommodation in the active site of each enzyme, thus providing a means of selectivity for the inhibitor [24].
ZINC390718 is a symmetric molecule characterized by the presence of nitrogen groups linked to two rings, as shown in Figure 3. A class of organic compounds with potential anticholinesterase activity is alkaloids, which also contain nitrogen bases as one of their main characteristics [25].
The benzimidazole radical (a benzene fused to the pentagonal ring) present at the ends of ZINC390718 gives it an association with imidazole alkaloids, a group that exhibits many pharmacological activities [26]. The research carried out by Gurjar et al. [27] shows the imidazole analogues as potential cholinesterase inhibitors with a neuroprotective role for Alzheimer’s disease. Hence, ZINC390718 represents a promising bioactive molecule, and future research aiming to improve its molecular arrangement could potentially make it more active against its biological targets.

2.2. In Vitro Neurotoxic Effects

The ZINC390718 molecule was subjected to a primary assessment of its toxicity against glial cells and astrocytes, and the results are shown in Figure 4. Cell viability analyses, carried out by measuring the function of mitochondrial dehydrogenases in converting MTT into crystals of formazan, allowed for the interpretation that, after 24 h of exposure, the study molecule (ZINC390718) had a concentration-dependent effect. The statistically significant interference with cell viability occurred only at the highest concentration tested (100 µM).
The percentage of viable cells was proven to be high for the other concentrations (between 88 and 95%), considering the percentage of the negative control (95%), which suggests a low toxicity. The MTT assay investigates whether cells are metabolically active to provide their dehydrogenase with enough reduced NADH/NADPH cofactors to reduce MTT to blue formazan [28].
In the analysis of the cultures treated with the molecule under study by using phase-contrast microscopy, no changes in the morphology of the treated cells, at different concentrations in relation to the negative control, were identified.
The concentrations evaluated in the toxicity test were lower than the values found for IC50 in the enzymatic assays with the cholinesterases. However, for this cell culture model, the concentration of 100 mM can be considered high. Higher concentrations can interfere with the culture medium and, consequently, cell viability.
Other studies have also used astrocyte cultures to assess the potential for chemical cytotoxicity. In a previous study, an assay was used to test a library with 80 compounds in MTT, and 4 already known cytotoxic compounds (valinomycin; 3,3′, 5,5′-tetrabromobisphenol; deltamethrin; and triphenyl phosphate) were used as controls in populations of neural cells derived from induced pluripotent stem cells (iPSCs) (neuroepithelial stem cells (NSCs), neurons, and astrocytes) at concentrations of 1, 10, and 100 µM [29]. Although 1 µM did not indicate the cytotoxicity of any of the tested compounds already known to be cytotoxic, at a concentration of 10 µM, it was possible to observe a cytotoxicity profile and the loss of cell viability caused, for example, by deltamethrin.
Historically, drug development programs for neurodegenerative diseases generally only target neurons, ignoring the contributions of astrocytes and instead targeting them, offering a new approach to drug development for the treatment of neurological diseases [30]. Preliminary tests on cells of great importance, such as glial cells, for the identification of the adequate functional capacity of the central nervous system, including neurodegenerative diseases, can lead to future evaluations that may bring more conclusive results from the perspective of the cytotoxicity of potentially active molecules in biological targets.

2.3. Molecular Dynamic Simulations

The molecular dynamics (MD) simulations allowed us to observe the stability and patterns of the established interactions of ZINC390718 with the active sites of the enzymes AChE and BuChE. To ensure reliable data acquisition, one needs to ensure that the system is stable, and one of the main ways to assess structural stability is by using RMSD graphs. This metric is used to analyze the equilibrium and stabilization of the protein structure and its system during the trajectory of MD simulations [31,32,33]. The complexes (ZINC90718 bound to AChE and BuChE) and the APO forms of the proteins were obtained from the RMSD for the main chain (backbone), along with the trajectory (see Figure 5).
In MD simulations, evaluations of the graphs allow for observations of progressive increases in the RMSD values of a protein until the values stabilize. The stability of AChE in the APO form (enzyme without the ligand) was observed from 60 ns (the production phase) of the simulation, where the mean RMSD value was 2.5 ± 0.3 Å; as the enzyme complexed with the ligand (ZINC390718), its production phase was from 15 ns, and the average RMSD value was 2.5 ± 0.3 Å. The stability of BuChE in the APO form was observed from 20 to 95 ns of the simulation, where the mean RMSD value was 2.6 ± 0.3 Å; as the enzyme complexed with the ligand (ZINC390718), its productive phase was 5 ns, and the mean RMSD value was 1.9 ± 0.19 Å.
An RMSD lower than 3.0 Å is considered appropriate when all atoms of the system are used in the evaluation, considering that the equilibrium variation occurs according to the intrinsic characteristics of the protein. The conformational changes that occurred in the AChE and BuChE structures in relation to the initial structures and the identification of the moment when the systems reached equilibrium were also analyzed using the RMSD metric, which mathematically suggests subtle changes between the patterns of the enzymatic systems presented by the APO forms and complexed with the ZINC390718 molecule when compared to each other [34].
The RMSD data are not enough to guarantee the stability of a system, as they do not score the variations in certain regions of the structure. Therefore, an RMSF plot was also generated in an attempt to evaluate the fluctuations of the atoms of each residue during the trajectory. Thus, waste fluctuation graphs were generated for the APO form and compared with those generated for the respective complexes with ZINC390718 during the production phase (Figure 6).
The analysis of the RMSF values demonstrated the atomic fluctuations of the APO form of AChE (RMSF = 0.8 ± 0.4 Å) and its complex (RMSF = 0.9 ± 0.4 Å), as well as the APO form of BuChE (RMSF = 1.0 ± 0.5 Å) and its complex (RMSF = 0.9 ± 0.50 Å). Comparing the regions with the greatest fluctuations between the APO and complex systems (Figure 6), it was observed that the stability pattern remained approximate, with statistically non-significant variations in the values.
The divergences in the findings between AChE and BuChE can be explained by differences in the structures around the main entrance area of the enzyme. In AChE, there are eight aromatic residues: three of them are located at the entrance to the active site (Tyr 72, Trp 286, and Tyr 341) and five are located at nearly the same depth within the slit (Tyr 124, Phe 295, Phe 297, Tyr 337, and Phe 338).
In BuChE, however, six of these residues are altered and substituted for minor residues; that is, Tyr 72 is substituted for Asn 68, Tyr 124 is substituted for Gln 119, Trp 286 is substituted for Ala 277, Phe 295 is substituted for Leu 286, Phe 297 is substituted for Val 288, and Tyr 337 is substituted for Ala 328. These smaller residues significantly extend the entry radius, which makes it possible for large substrates to enter the active site [35].
The data presented herein indicate that the presence of ZINC390718 possibly stabilizes the structures of AChE and BuChE. However, in order to understand the stabilization of the systems, we investigate the intermolecular interactions related to the processes that likely contribute to obtaining a biological response. Thus, hydrogen bonds stand out, since they correspond to the main interactions involved in maintaining the tertiary structure and protein folding, as well as in molecular recognition [36].
Hydrogen interactions have transient characteristics, making it necessary to evaluate their occurrence and permanence during the simulation time. In order to identify whether the stability of the systems, the objective of this study, occurs due to these interactions, Table 1 shows the number and permanence of H bonds formed between the ligand ZINC390718 and the amino acids of the active sites of AChE and BuChE.
The residues that H-bonded with AChE with a permanence greater than 10% were Tyr124 with N1 (16.52%) and N2 (23.36%), Phe295 with O1 (86.61%), Arg296 with O1 (21.37%), and Ser293 with N4 (55.84%). Based on these data, the interactions of the H bond with the amino acids stand out in terms of their relationship with AChE’s active sites. The residue Tyr124 participates in the anionic peripheral active site, and Phe295 is part of the acyl pocket and presents a H interaction with a longer residence time. The residues Tyr124 and Phe295 perform interactions with cannabinoids and molecules that inhibit the enzymatic activity of AChE [37].
AChE also has other important binding sites for ligand-receptor mechanisms: the peripheral anionic site (PAS) (Tyr103, Asp105, Tyr155, Glu316, Trp317, and Tyr372), the acyl pocket (Trp117), the anionic site (Trp267, Phre26, and Phe328), the oxyanionic orifice (Gly152, Gly153, and Ala235), and the aromatic site (Tyr164 and Tyr368) [38].
The interactions of H with BuChE, which remained greater than 10%, were observed with the residues Asp70 (37.9%), Glu197 (12.4%), Tyr332 (38.4%), and Lis339 (33.0%). Asp70 and Tyr332 are residues located at the peripheral site of BuChE; these residues interact with the substrate and change the conformation of the protein to allow for catalysis [39].
The residue Tyr332 has been implicated in binding substrates to BuChE, suggesting that this amino acid residue is part of the P (peripheral) site in this enzyme; anionic Asp70 (Asp74 in AChE) is linked by H interactions, and Tyr332 is the other essential residue of the P site of BuChE [40].
The AChE and BuChE cavities resemble a deep, narrow canyon, about 20 Å long in both enzymes. The catalytic triad consists of serine, histidine, and glutamate residues. The peripheral site that influences the catalytic process is located above the catalytic triad towards the periphery of the active site canyon. This peripheral site contains several aromatic amino acid residues, as well as the aspartate residue (Asp70), an amino acid from BuChE that has been found to have one of the highest percentages of permanence in interactions with H [38].
To select the representative structures of the systems, the most adequate RMSD values were 1.2 Å (AChE) and 1.3 Å (BuChE), as they exhibited the largest number of clusters and significant diversity, and of the total system conformations during the production phase, more than 50% presented in the first three groups for the three systems [41]. The representative structures selected were those that were the most frequent in the groups with the highest number of conformations observed for the cutoff point.
The cutoff point was defined based on the concept that the greater the number of clusters observed in a simulation, the greater the diversity of the conformations adopted by the complex. Moreover, the representative structures were selected via the evaluation of the clusters, using the differences in the maximum RMSD values among the conformations to determine the cutoff point while still considering the native state of the protein as the most stable conformation, which certainly is the most frequent and the most populous group [42].
The representative structures below were selected based on their frequency in the groups with the highest number of conformations observed for the cutoff point (AChE: 26.5 ns; BuChE: 39.15 ns). The set of intermolecular interactions generated between the targets (AChE and BuChE) and the ligand (ZINC390718) responsible for the stabilization of the systems was illustrated using PLIP (Figure 7) after identifying the representative structures for the complexes.
When evaluating the representative structure of the ZINC390718 complex with AChE and BuChE, hydrophobic and hydrogen interactions are observed. The ligand has a benzimidazole radical (a benzene fused to the pentagonal ring) at each of its ends. In the interactions of AChE with the ligand, one of the benzimidazole radicals is involved in the hydrophobic interactions with the residues Leu76, Trp86, and Tyr341, and the N atom establishes a hydrogen interaction with Tyr124.
However, the N opposite to the benzimidazole radical performs a hydrogen interaction with Ser293. Additionally, a hydrophobic interaction is observed between the linker chain and Phe123.
The attractive and non-covalent interactions between the ligand with the residues Tyr 341 and Trp 86 play an important role in the stabilization of the inhibitor in the active site [43]. The peripheral site of AChE is composed of numerous aromatic chains that extend beyond the Tyr337 residue at the catalytic site interface to the entrance of the canyon, contributing to the catalytic efficiency of AChE. This site includes Phe295, Tyr337, and Tyr341 [44].
One end of the region of the linker, the benzimidazole radical, interacts with the Asp70 and Tyr332 residues of BuChE through hydrophobic interactions. The other end, which is also formed by the same benzimidazole radical, interacts hydrophobically with the residues Tyr128, Ile442, Trp82, and Glu197, the latter also interacting with hydrogen as a donor with one of the nitrogen atoms present in this radical.
It has been reported in the literature that the occurrence of hydrogen interactions with the Asp70 residue is important for inhibitory activity against BuChE, because it promotes the blocking of the active site and prevents substrate access to this region; consequently, catalysis does not occur [45].
The interaction with Trp82 (also observed in the structure obtained via molecular docking) is observed in the MD simulations with BuChE inhibitors, and this residue even interacts with the benzene ring of tacrine (cholinesterase inhibitor), playing a crucial role in stabilizing the substrate in the active site [46].
In the studies performed by Zhang et al. [47], the residues, including Asp70, Trp82, and Tyr332, were identified as playing important roles in stabilizing tacrine in the BuChE pocket, where important interactions, such as hydrogen interactions, were observed. Tyr332 and Asp70 residues, which act as gatekeepers, play crucial roles in substrate adjustment.

3. Materials and Methods

3.1. Reagents and Equipment

For the evaluations of anticholinesterase activities, acetylcholinesterase (AChE) (Electrophorus electricus type VI) and butyrylcholinesterase (BuChE obtained from equine serum), the substrate acetylthiocholine (ACTI), Ellman’s reagent, 5.5′ditiobis (2nitrobenzoic acid) (DTNB), bovine serum albumin (BSA), and the ZINC390718 molecule were obtained commercially from Sigma-Aldrich®. A Multiskan FC microplate reader was used and was obtained from Thermo Scientific® (Waltham, MA, USA).
For the cytotoxicity assay, the following reagents were used: Dulbecco’s modified Eagle’s medium (DMEM, Cultilab, Campinas, SP, Brazil), penicillin, streptomycin, amphotericin B (Sigma-Aldrich®, São Paulo, SP, Brazil), and fetal bovine serum (SFB, Gibco, Grand Island, NY 12657, USA). The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) bromide reduction method was used (Sigma Aldrich, St. Louis, MO, USA).
The ZINC390718 molecule, evaluated in this study, was selected according to in silico data (pharmacophore-based virtual screening and a molecular docking analysis) obtained by our research group [15].

3.2. Determination of the Kinetic Parameters of AChE and BuChE

The standardization of the AChE and BuChE catalytic parameters was evaluated according to Ellman et al. [23] with modifications. The determination of the enzymatic concentration was carried out by keeping the concentration of the used substrate (acetylthiocholine: ACTI) fixed and by varying the enzyme concentrations (0.015, 0.075, and 0.30 U/mL) to evaluate the enzyme reaction rate and substrate consumption. The optical absorbance of the samples was quantified using a Multiskan FC microplate reader (Thermo Scientific®) at a wavelength of 405 nm and a temperature of 25 °C.
After the determination of the enzymatic concentrations, the Michaelis–Menten constant (Km value) for ACTI was determined by varying its concentration (0.0625; 0.1325; 0.2650; 0.5300; 1; 2; 4 mM) for each cholinesterase. The experiments were performed in triplicate, and the Km and Vmax values were obtained using non-linear regression (Michaelis–Menten equation), utilizing the GraphPad Prism program (version 5.0).

3.3. In Vitro Anticholinesterase Activity

The evaluation of anticholinesterase activity was performed as described by Ellman et al. [23] adapted for 96-well microplates, using the kinetic parameters previously determined in Section 2.2.
Pipetting in all wells was initially carried out in the experiments using the following: 140 μL of 0.1 M phosphate buffer at pH 7.4 containing 0.1% bovine serum albumin (BSA); then, in separate triplicates, 20 μL of the sample (molecule ZINC390718 diluted in ethanol), positive control (0.5 mM eserine), and controls (10% ethanol and phosphate buffer); and 20 μL of enzyme in all triplicates (AchE: at a final concentration of 0.15 U/m and BuChE: at a concentration of 0.50 U/mL, both diluted in 0.1 M phosphate buffer). After a 10-minute incubation at room temperature (25 °C), 10 μL of DTNB (10 mM in the well) and 10 μL of ACTI (0.125 mM in the well for AChE, and 0.744 mM in the well for BuChE) were added.
The sample (ZINC390718) was tested at concentrations between 37.5 and 1000 μM and 50 and 500 μM to calculate the IC50 values of AChE and BuChE, respectively. The readings were performed in a cycle from 0 to 11 per minute for AChE and in a cycle from 0 to 21 per minute for BuChE. The optical absorbance of the samples was measured at 405 nm using a Multiskan FC microplate reader (Thermo Scientific®).

3.4. In Vitro Evaluation of Neurotoxic Effects

3.4.1. Primary Astrocyte-Enriched Glial Cell Culture

The glial cell culture was obtained from the cortex of newborn Wistar rats, 0-2 days after their birth. The animals were acquired from the Vivarium of the School of Veterinary Medicine and Animal Sciences of the Federal University of Bahia (UFBA). The procedures were carried out in accordance with the Brazilian guidelines for the production, maintenance, and use of animals for teaching and scientific research activities of the Local Animal Experimentation Ethics Committee (CEUA) of the UFBA Institute of Health Sciences, approved under the registration no. 027/2012, composing the project entitled “Pharmacological bioprospecting of natural products and derivatives in cells of the central nervous system”.
According to the protocol established by Pitanga et al. [48] and Silva et al. [49], glial cells were obtained from the beheading of the animals, followed by the removal of the brain hemispheres aseptically. The meninges and blood vessels were removed from each cortex, and then the material was mechanically dissociated and filtered through a sterile membrane (75 mm). The dissociated cells were grown in bottles with Dulbecco’s modified Eagle’s medium (DMEM), 50% glucose (33 mM), penicillin (100 IU/mL), streptomycin (100 μg/mL), NaHCO3 (2.438 g/L), 6.25 μg/mL amphotericin B, and 10% SFB. The cells were kept in a mixed glial culture for up to 7 days in culture bottles. The microglial cells were isolated from the mixed glial culture by shaking at 1000 RPM for 3 h at 37 °C. After the end of the agitation, the medium was collected, and, thus, the astrocyte-enriched glia remained adhered to the culture bottle. After trypsinization using a trypsin/EDTA solution (0.05% trypsin and 0.02% EDTA in PBS free of Ca2+/Mg2+) and counting in a Neubauer chamber, the cells were seeded at a density of 3.3 × 104 in 96-well plates, previously treated with poly-D-lysine and cultured in a humid chamber with 5% CO2 at 37 °C.

3.4.2. Analysis of Cytotoxicity through the MTT Assay

Cell viability was measured using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium (MTT) bromide reduction method according to Hansen et al. [50]. Cell cultures of normal glial cells were grown in 96-well culture plates, and after 24 h, they were exposed to the molecule in increasing concentrations (0.1, 1, 10, 50, and 100 µM) and to the control (DMSO 0.01%). The experiments were analyzed after 24 h of treatment. After the end of the treatment exposure time, the MTT solution (100 µL/well) diluted in DMEM, at a ratio of 1:4, was added. The plate containing the solution was kept incubated in an oven at 37 °C for 2 h, and then, for the lysis of the cells and complete dissolution of the formazan crystals, a volume of 100 µL/well of the lysis buffer was added, composed of 20% of sodium duodecyl sulfate (SDS) and 50% dimethyl formamide (DMF), pH 4.7.
The cultivation plate was kept at rest on the bench for 18 h, properly protected from light at 37 °C. After this time, the optical absorbance of each sample was quantified using a spectrophotometer (VARIOSKAN FLASH, Thermo Fisher Scientific) at a wavelength of 595 nm. An empty well for each treatment was used as a blank, and its absorbance value was discounted from all others. The experiments were performed in eight replicates, and the results are presented as the percentage of viability (mean and standard deviation) in relation to the control. Cell viability was calculated as a percentage of the negative control, which was considered to be 100%.

3.4.3. Analysis of Cell Morphology Using Phase-Contrast Microscopy

The effects of the molecule on cell morphology were initially observed using phase-contrast microscopy. After 24 h of exposure to the molecule at different concentrations, the cell cultures were analyzed using phase-contrast microscopy with a green filter (Nikon TS-100), and then 3 fields were photographed for each treatment.

3.5. Statistical Analysis of In Vitro Tests

All generated data received adequate statistical treatment using GraphPad Prism software, version 5.0, for Windows (GraphPad Software, San Diego, CA, USA), and they are expressed as means ± standard deviation. The p values adopted as statistically significant in the analyses were those lower than 0.05. The IC50 values for anticholinestease activities were calculated using a non-linear regression analysis.
The selectivity index (IS) of the molecule action against BuChE was determined using the following formula: IC50AChE/IC50BuChE.

3.6. In Silico Evaluation: Molecular Dynamics

3.6.1. Ligand Topology

The best pose of ZINC390718 obtained via the molecular docking assays was submitted to the ATB 3.0 server (https://atb.uq.edu.au/index.py, accessed on 5 July 2021) for the generation of its topology. The parameters of atomic charge, bond length, torsional angles, and dihedrals were obtained using the GROMOS96 54A7 force field [51].

3.6.2. Molecular Dynamics Simulations

MD simulations were performed using the GROMACS 5.1.2 package [52], in which the GROMOS96 54A7 force field parameters, a temperature of 298.15 K, and a pressure of 1.013 bar were adopted.
The 3D structures were obtained from the PDB (PDB ID for AChE: 4M0E, and for BuChE: 4BDS), from which the crystallographic ligand, non-structural water molecules, and crystallization artifacts were removed. The non-modeled regions of the structures were built with the help of the SWISS-MODEL server [53].
The protonation states of the acidic and basic residues of the targets were adjusted in the pdb2gmx module implemented in GROMACS 5.1.2 according to pH 7.4 for AChE and BuChE [54]. For this, the pKa values of the residues were evaluated on the H ++ server (http://biophysics.cs.vt.edu/index.php, accessed on 5 July 2021) [55]. To solvate the systems, a dodecahedral box with the water model SPC-E [56] was used, with a minimum distance of 1.4 nm from the edge of the box. For neutralization, in systems involving AChE (APO and complex), seven Na2+ ions were added, while in systems with BuChE, four Cl− ions were added.
All systems (APO and complex forms) were minimized in two steps: initially by using the Steepest Descent (SD) algorithm with 10,000 cycles and, later, by using the Conjugated Gradient (GC) algorithm with 1000 cycles. After the minimization steps, the balancing step (t = 1.0 ns) and, finally, the production dynamics (APO systems: t = 100.0 ns; complex systems: t = 50.0 ns) were carried out. The MD data were obtained from the isothermal–isobaric ensemble (NPT—number of particles, constant temperature, and pressure).
The electrostatic and hydrophobic interactions were described using the Particle Mesh Ewald (PME) method [57], with a cutting radius equal to 0.9 nm. The stability of the system was evaluated based on the analysis of the variation in the Root Mean Square Deviation (RMSD) and the Root Mean Square Fluctuation (RMSF), calculated through the modules rms and rmsf available in GROMACS 5.1.2.50.
Additionally, the number and permanence of hydrogen bonds in the complexes were calculated using the g_hbond module, implemented in the GROMACS 5.1.2 package and the HbMap2Grace program (http://lmdm.biof.ufrj.br/software/hbmap2grace/index.html, accessed on 5 July 2021). The hydrogen bonds that showed a distance between the donor and acceptor that was lower than or equal to 3.5 Å and an angle between the donor/acceptor and the H atom that was lower than or equal to 60° were identified, as well as the bonds whose permanence time was lower than 10% of the total simulation time [58].

3.6.3. Clustering Analysis and Selection of Average Structure

A representative structure of the ligand-macromolecule complex during the MD production phase was obtained by grouping similar conformations with the aid of the g_cluster module available in the GROMACS 5.1.2 package, using the joining method GROMOS [59], with cutoff points defined from the RMSD values of 1.0; 1.1; 1.2; 1.3; 1.4; 1.5 Å. The average structure of the most populous group was then selected for the analysis of the main interactions established between the prioritized ligand and the targets.

4. Conclusions

The ZINC390718 molecule had dual in vitro inhibitory activity against AChE and BuChE, acting in a concentration-dependent manner and being more active against BuChE. These in vitro effects are in agreement with the in silico findings. The MD simulation revealed that ZINC390718 performed important hydrophobic and H-interactions with residues that are part of AChE and BuChE catalytic sites. The residues Tyr341, Phe338, Tyr337, Val365, Val 294, Trp 86, Phe 295, and Gly342 participated in interactions with AChE, and the residues Asp70, Tyr332, Tyr128, Ile442, Trp82, and Glu197 participated in interactions with BuChE. These results support the ability of ZINC390718 to interact at the active sites of AChE and BuChE.

Author Contributions

The manuscript consists of contributions of all the aforementioned authors. All authors contributed to the conceptualization of the project and methodology. Conceptualization, R.B.M.d.A., F.H.A.L. and M.B.B.; Data curation, D.B.B., J.A.O.A., B.S.A., S.L.C. and J.M.C.; Formal analysis, R.B.M.d.A., D.B.B., M.R.d.B., J.A.O.A., S.L.C., J.M.C., J.N.C., C.B.R.S. and M.B.B.; Funding acquisition, J.M.C. and C.B.R.S.; Investigation, R.B.M.d.A., M.R.d.B., J.A.O.A. and S.L.C.; Methodology, D.B.B., M.R.d.B., J.A.O.A., B.S.A., F.H.A.L. and M.B.B.; Software, B.S.A. and F.H.A.L.; Supervision, F.H.A.L. and M.B.B.; Validation, M.R.d.B. and B.S.A.; Visualization, D.B.B., S.L.C., J.N.C. and C.B.R.S.; Writing—original draft, R.B.M.d.A.; Writing—review and editing, J.M.C., J.N.C., C.B.R.S., F.H.A.L. and M.B.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundação de Amparo à Pesquisa do Estado da Bahia (FAPESB, BOL0346/2017) and Programa de Pós-Graduação em Biotecnologia of Universidade Estadual de Feira de Santana (UEFS).

Institutional Review Board Statement

The animal study protocol was approved by the Local Animal Experimentation Ethics Committee (CEUA) of the UFBA Institute of Health Sciences under the registration no. 027/2012.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors would like to thank the Programa de Pós-Graduação em Biotecnologia of Universidade Estadual de Feira de Santana (UEFS) and the Fundação de Amparo à Pesquisa do Estado da Bahia (FAPESB). We also thank Abilio Borghi for assistance with the language review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ambure, P.; Bhat, J.; Puzyn, T.; Roy, K. Identifying natural compounds as multi-target-directed ligands against Alzheimer’s disease: An in silico approach. J. Biomol. Struct. Dyn. 2019, 37, 1282–1306. [Google Scholar] [CrossRef] [PubMed]
  2. WHO. Vector-Borne Diseases; National Academies Press: Washington, DC, USA, 2008. [CrossRef]
  3. Long, J.M.; Holtzman, D.M. Alzheimer Disease: An Update on Pathobiology and Treatment Strategies. Cell 2019, 179, 312–339. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, X.Q.; Mobley, W.C. Exploring the pathogenesis of Alzheimer disease in basal forebrain cholinergic neurons: Converging insights from alternative hypotheses. Front. Neurosci. 2019, 13, 446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ferreira-Vieira, T.H.; Guimaraes, I.M.; Silva, F.R.; Ribeiro, F.M. Alzheimer’s disease: Targeting the Cholinergic System. Curr. Neuropharmacol. 2016, 14, 101–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Da Silva, V.B.; De Andrade, P.; Kawano, D.F.; Morais, P.A.B.; De Almeida, J.R.; Carvalho, I.; Taft, C.A.; Da Silva, C.H.T.D.P. In silico design and search for acetylcholinesterase inhibitors in Alzheimer’s disease with a suitable pharmacokinetic profile and low toxicity. Future Med. Chem. 2011, 3, 947–960. [Google Scholar] [CrossRef]
  7. Ciriello, R.; Lo Magro, S.; Guerrieri, A. Assay of serum cholinesterase activity by an amperometric biosensor based on a co-crosslinked choline oxidase/overoxidized polypyrrole bilayer. Analyst 2018, 143, 920–929. [Google Scholar] [CrossRef] [Green Version]
  8. Hampel, H.; Mesulam, M.M.; Cuello, A.C.; Farlow, M.R.; Giacobini, E.; Grossberg, G.T.; Khachaturian, A.S.; Vergallo, A.; Cavedo, E.; Snyder, P.J.; et al. The cholinergic system in the pathophysiology and treatment of Alzheimer’s disease. Brain 2018, 141, 1917–1933. [Google Scholar] [CrossRef]
  9. Li, Y.; Hai, S.; Zhou, Y.; Dong, B.R. Cholinesterase inhibitors for rarer dementias associated with neurological conditions. Cochrane Database Syst. Rev. 2015, 2015, CD009444. [Google Scholar] [CrossRef] [Green Version]
  10. Greig, N.H.; Lahiri, D.K.; Sambamurti, K. Butyrylcholinesterase: An important new target in Alzheimer’s disease therapy. Int. Psychogeriatr. 2002, 14, 77–91. [Google Scholar] [CrossRef]
  11. Eldufani, J.; Blaise, G. The role of acetylcholinesterase inhibitors such as neostigmine and rivastigmine on chronic pain and cognitive function in aging: A review of recent clinical applications. Alzheimer’s Dement. Transl. Res. Clin. Interv. 2019, 5, 175–183. [Google Scholar] [CrossRef]
  12. Larik, F.A.; Shah, M.S.; Saeed, A.; Shah, H.S.; Channar, P.A.; Bolte, M.; Iqbal, J. New cholinesterase inhibitors for Alzheimer’s disease: Structure activity relationship, kinetics and molecular docking studies of 1-butanoyl-3-arylthiourea derivatives. Int. J. Biol. Macromol. 2018, 116, 144–150. [Google Scholar] [CrossRef]
  13. Li, S.Y.; Jiang, N.; Xie, S.S.; Wang, K.D.G.; Wang, X.B.; Kong, L.Y. Design, synthesis and evaluation of novel tacrine-rhein hybrids as multifunctional agents for the treatment of Alzheimer’s disease. Org. Biomol. Chem. 2014, 12, 801–814. [Google Scholar] [CrossRef] [PubMed]
  14. Sharma, K. Cholinesterase inhibitors as Alzheimer’s therapeutics. Mol. Med. Rep. 2019, 20, 1479–1487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Mascarenhas, A.M.S.; de Almeida, R.B.M.; de Araujo Neto, M.F.; Mendes, G.O.; da Cruz, J.N.; Santos, C.B.R.; Botura, M.B.; Leite, F.H.A. Pharmacophore-based virtual screening and molecular docking to identify promising dual inhibitors of human acetylcholinesterase and butyrylcholinesterase. J. Biomol. Struct. Dyn. 2020, 39, 6021–6030. [Google Scholar] [CrossRef] [PubMed]
  16. Almeida, V.M.; Dias, Ê.R.; Souza, B.C.; Cruz, J.N.; Santos, C.B.R.; Leite, F.H.A.; Queiroz, R.F.; Branco, A. Methoxylated flavonols from Vellozia dasypus Seub ethyl acetate active myeloperoxidase extract: In vitro and in silico assays. J. Biomol. Struct. Dyn. 2022, 40, 7574–7583. [Google Scholar] [CrossRef] [PubMed]
  17. Rego, C.M.A.; Francisco, A.F.; Boeno, C.N.; Paloschi, M.V.; Lopes, J.A.; Silva, M.D.S.; Santana, H.M.; Serrath, S.N.; Rodrigues, J.E.; Lemos, C.T.L.; et al. Inflammasome NLRP3 activation induced by Convulxin, a C-type lectin-like isolated from Crotalus durissus terrificus snake venom. Sci. Rep. 2022, 12, 4706. [Google Scholar] [CrossRef]
  18. Santos, C.B.R.; Santos, K.L.B.; Cruz, J.N.; Leite, F.H.A.; Borges, R.S.; Taft, C.A.; Campos, J.M.; Silva, C.H.T.P. Molecular modeling approaches of selective adenosine receptor type 2A agonists as potential anti-inflammatory drugs. J. Biomol. Struct. Dyn. 2021, 39, 3115–3127. [Google Scholar] [CrossRef]
  19. Ganesan, A.; Coote, M.L.; Barakat, K. Molecular dynamics-driven drug discovery: Leaping forward with confidence. Drug Discov. Today 2017, 22, 249–269. [Google Scholar] [CrossRef] [PubMed]
  20. Bisswanger, H. Enzyme assays. Perspect. Sci. 2014, 1, 41–55. [Google Scholar] [CrossRef]
  21. Hubbell, J.A.E.; Muir, W.W. Monitoring Anesthesia. In Equine Anesthesia: Monitoring and Emergency Therapy, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2008; pp. 149–170. [Google Scholar] [CrossRef]
  22. Robinson, P.K. Enzymes: Principles and biotechnological applications. Essays Biochem. 2015, 59, 1–41. [Google Scholar] [CrossRef]
  23. Ellman, G.L.; Courtney, K.D.; Andres, V.; Featherstone, R.M. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961, 7, 88–95. [Google Scholar] [CrossRef]
  24. Darvesh, S.; McDonald, R.S.; Darvesh, K.V.; Mataija, D.; Conrad, S.; Gomez, G.; Walsh, R.; Martin, E. Selective reversible inhibition of human butyrylcholinesterase by aryl amide derivatives of phenothiazine. Bioorganic Med. Chem. 2007, 15, 6367–6378. [Google Scholar] [CrossRef]
  25. Eguchi, R.; Ono, N.; Hirai Morita, A.; Katsuragi, T.; Nakamura, S.; Huang, M.; Altaf-Ul-Amin, M.; Kanaya, S. Classification of alkaloids according to the starting substances of their biosynthetic pathways using graph convolutional neural networks. BMC Bioinform. 2019, 20, 380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Jin, Z. Muscarine, imidazole, oxazole and thiazole alkaloids. Nat. Prod. Rep. 2009, 26, 382–445. [Google Scholar] [CrossRef] [PubMed]
  27. Gurjar, A.S.; Darekar, M.N.; Yeong, K.Y.; Ooi, L. In silico studies, synthesis and pharmacological evaluation to explore multi-targeted approach for imidazole analogues as potential cholinesterase inhibitors with neuroprotective role for Alzheimer’s disease. Bioorganic Med. Chem. 2018, 26, 1511–1522. [Google Scholar] [CrossRef] [PubMed]
  28. Berridge, M.V.; Tan, A.S. Characterization of the Cellular Reduction of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT): Subcellular Localization, Substrate Dependence, and Involvement of Mitochondrial Electron Transport in MTT Reduction. Arch. Biochem. Biophys. 1993, 303, 474–482. [Google Scholar] [CrossRef] [PubMed]
  29. Pei, Y.; Peng, J.; Behl, M.; Sipes, N.S.; Shockley, K.R.; Rao, M.S.; Tice, R.R.; Zeng, X. Comparative neurotoxicity screening in human iPSC-derived neural stem cells, neurons and astrocytes. Brain Res. 2016, 1638, 57–73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Gorshkov, K.; Aguisanda, F.; Thorne, N.; Zheng, W. Astrocytes as targets for drug discovery. Drug Discov. Today 2018, 23, 673–680. [Google Scholar] [CrossRef]
  31. Tsay, J.J.; Su, S.C.; Yu, C.S. A multi-objective approach for protein structure prediction based on an energy model and backbone angle preferences. Int. J. Mol. Sci. 2015, 16, 15136–15149. [Google Scholar] [CrossRef] [Green Version]
  32. dos Santos, K.L.B.; Cruz, J.N.; Silva, L.B.; Ramos, R.S.; Neto, M.F.A.; Lobato, C.C.; Ota, S.S.B.; Leite, F.H.A.; Borges, R.S.; da Silva, C.H.T.P.; et al. Identification of novel chemical entities for adenosine receptor type 2a using molecular modeling approaches. Molecules 2020, 25, 1245. [Google Scholar] [CrossRef] [Green Version]
  33. de Melo Lima, A.; Siqueira, A.S.; Möller, M.L.S.; de Souza, R.C.; Cruz, J.N.; de Lima, A.R.J.; da Silva, R.C.; Aguiar, D.C.F.; da Silva Gonçalves Vianez, J.L., Jr.; Gonçalves, E.C. In silico improvement of the cyanobacterial lectin microvirin and mannose interaction. J. Biomol. Struct. Dyn. 2020, 40, 1064–1073. [Google Scholar] [CrossRef] [PubMed]
  34. Sargsyan, K.; Grauffel, C.; Lim, C. How Molecular Size Impacts RMSD Applications in Molecular Dynamics Simulations. J. Chem. Theory Comput. 2017, 13, 1518–1524. [Google Scholar] [CrossRef] [PubMed]
  35. 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] [PubMed] [Green Version]
  36. Hubbard, R.E.; Kamran Haider, M. Hydrogen Bonds in Proteins: Role and Strength. eLS, 2010. [Google Scholar] [CrossRef]
  37. Furqan, T.; Batool, S.; Habib, R.; Shah, M.; Kalasz, H.; Darvas, F.; Kuca, K.; Nepovimova, E.; Batool, S.; Nurulain, S.M. Cannabis constituents and acetylcholinesterase interaction: Molecular docking, in vitro studies and association with CNR1 RS806368 and ACHE RS17228602. Biomolecules 2020, 10, 758. [Google Scholar] [CrossRef]
  38. Kiametis, A.S.; Silva, M.A.; Romeiro, L.A.S.; Martins, J.B.L.; Gargano, R. Potential acetylcholinesterase inhibitors: Molecular docking, molecular dynamics, and in silico prediction. J. Mol. Model. 2017, 23, 67. [Google Scholar] [CrossRef]
  39. Bajda, M.; Łątka, K.; Hebda, M.; Jończyk, J.; Malawska, B. Novel carbamate derivatives as selective butyrylcholinesterase inhibitors. Bioorganic Chem. 2018, 78, 29–38. [Google Scholar] [CrossRef]
  40. Rosenberry, T.L.; Brazzolotto, X.; 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]
  41. Pita, S.S.R.; Pascutti, P.G. Pharmacophore analysis of Trypanosoma cruzi trypanothione reductase (TR) complexed with peptide mimetic inhibitors. Rev. Virtual Quim. 2012, 4, 788–804. [Google Scholar] [CrossRef]
  42. Torda, A.E.; van Gunsteren, W.F. Algorithms for clustering molecular dynamics configurations. J. Comput. Chem. 1994, 15, 1331–1340. [Google Scholar] [CrossRef]
  43. Jang, C.; Yadav, D.K.; Subedi, L.; Venkatesan, R.; Venkanna, A.; Afzal, S.; Lee, E.; Yoo, J.; Ji, E.; Kim, S.Y.; et al. Identification of novel acetylcholinesterase inhibitors designed by pharmacophore-based virtual screening, molecular docking and bioassay. Sci. Rep. 2018, 8, 14921. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Szegletes, T.; Mallender, W.D.; Thomas, P.J.; Rosenberry, T.L. Substrate binding to the peripheral site of acetylcholinesterase initiates enzymatic catalysis. Substrate inhibition arises as a secondary effect. Biochemistry 1999, 38, 122–133. [Google Scholar] [CrossRef]
  45. Senol, F.S.; Ślusarczyk, S.; Matkowski, A.; Pérez-Garrido, A.; Girón-Rodríguez, F.; Cerón-Carrasco, J.P.; den-Haan, H.; Peña-García, J.; Pérez-Sánchez, H.; Domaradzki, K.; et al. Selective in vitro and in silico butyrylcholinesterase inhibitory activity of diterpenes and rosmarinic acid isolated from Perovskia atriplicifolia Benth. and Salvia glutinosa L. Phytochemistry 2017, 133, 33–44. [Google Scholar] [CrossRef] [PubMed]
  46. 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. Bioorganic Chem. 2019, 85, 117–127. [Google Scholar] [CrossRef]
  47. Zhang, Z.; Fan, F.; Luo, W.; Zhao, Y.; Wang, C. Molecular Dynamics Revealing a Detour-Forward Release Mechanism of Tacrine: Implication for the Specific Binding Characteristics in Butyrylcholinesterase. Front. Chem. 2020, 8, 730. [Google Scholar] [CrossRef] [PubMed]
  48. Pitanga, B.P.S.; Silva, V.D.A.; Souza, C.S.; Junqueira, H.A.; Fragomeni, B.O.N.; Nascimento, R.P.; Silva, A.R.; de Fátima, D.; Costa, M.; El-Bachá, R.S.; et al. Assessment of neurotoxicity of monocrotaline, an alkaloid extracted from Crotalaria retusa in astrocyte/neuron co-culture system. Neurotoxicology 2011, 32, 776–784. [Google Scholar] [CrossRef]
  49. Silva, V.D.A.; Cuevas, C.; Muñoz, P.; Villa, M.; Ahumada-Castro, U.; Huenchuguala, S.; Dos Santos, C.C.; De Araujo, F.M.; Ferreira, R.S.; Da Silva, V.B.; et al. Autophagy protects against neural cell death induced by piperidine alkaloids present in Prosopis juliflora (Mesquite). An. Acad. Bras. Cienc. 2017, 89, 247–261. [Google Scholar] [CrossRef] [Green Version]
  50. Hansen, M.B.; Nielsen, S.E.; Berg, K. Re-examination and further development of a precise and rapid dye method for measuring cell growth/cell kill. J. Immunol. Methods 1989, 119, 203–210. [Google Scholar] [CrossRef]
  51. Schmid, N.; Eichenberger, A.P.; Choutko, A.; Riniker, S.; Winger, M.; Mark, A.E.; Van Gunsteren, W.F. Definition and testing of the GROMOS force-field versions 54A7 and 54B7. Eur. Biophys. J. 2011, 40, 843–856. [Google Scholar] [CrossRef]
  52. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindah, E. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  53. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; De Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Quinn, B. Preparation and Maintenance of Live Tissues and Primary Cultures for Toxicity Studies. In Biochemical Ecotoxicology: Principles and Methods; Elsevier: Amsterdam, The Netherlands, 2014; pp. 33–47. [Google Scholar] [CrossRef]
  55. Mermelstein, D.J.; McCammon, J.A.; Walker, R.C. pH-dependent conformational dynamics of beta-secretase 1: A molecular dynamics study. J. Mol. Recognit. 2019, 32, e2765. [Google Scholar] [CrossRef] [PubMed]
  56. Berendsen, H.J.C.; Grigera, J.R.; Straatsma, T.P. The missing term in effective pair potentials. J. Phys. Chem. 1987, 91, 6269–6271. [Google Scholar] [CrossRef]
  57. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  58. Araujo, J.S.C.; de Souza, B.C.; Costa Junior, D.B.; de Mattos Oliveira, L.; Santana, I.B.; Duarte, A.A.; Lacerda, P.S.; dos Santos Junior, M.C.; Leite, F.H.A. Identification of new promising Plasmodium falciparum superoxide dismutase allosteric inhibitors through hierarchical pharmacophore-based virtual screening and molecular dynamics. J. Mol. Model. 2018, 24, 220. [Google Scholar] [CrossRef]
  59. Daura, X.; Van Gunsteren, W.F.; Mark, A.E. Folding-unfolding thermodynamics of a β-heptapeptide from equilibrium simulations. Proteins Struct. Funct. Genet. 1999, 34, 269–280. [Google Scholar] [CrossRef]
Pharmaceuticals 16 00095 g001 550
Figure 1. Determination of the kinetic parameters of AChE and BuChE under steady-state conditions. (A) Km and Vmax values of AChE against the substrate (ACTI). (B) Km and Vmax values of BuChE against the substrate (ACTI). Values were calculated using non-linear regression.
Figure 1. Determination of the kinetic parameters of AChE and BuChE under steady-state conditions. (A) Km and Vmax values of AChE against the substrate (ACTI). (B) Km and Vmax values of BuChE against the substrate (ACTI). Values were calculated using non-linear regression.
Pharmaceuticals 16 00095 g001
Pharmaceuticals 16 00095 g002 550
Figure 2. Concentration–response curve of the in vitro anticholinesterase effect of ZINC390718: (A) acetylcholinesterase (AChE) and (B) butyrylcholinesterase (BuChE).
Figure 2. Concentration–response curve of the in vitro anticholinesterase effect of ZINC390718: (A) acetylcholinesterase (AChE) and (B) butyrylcholinesterase (BuChE).
Pharmaceuticals 16 00095 g002
Pharmaceuticals 16 00095 g003 550
Figure 3. The 2D chemical structure of ZINC390718, C19H18N4O, with nitrogen groups (blue) and carbonyl (red) highlighted.
Figure 3. The 2D chemical structure of ZINC390718, C19H18N4O, with nitrogen groups (blue) and carbonyl (red) highlighted.
Pharmaceuticals 16 00095 g003
Pharmaceuticals 16 00095 g004 550
Figure 4. Percentage of cell viability in glial cell cultures of rats after 24 h of exposure to the Z90718 molecule (MTT assay). (*) Significant differences from negative control (p < 0.05).
Figure 4. Percentage of cell viability in glial cell cultures of rats after 24 h of exposure to the Z90718 molecule (MTT assay). (*) Significant differences from negative control (p < 0.05).
Pharmaceuticals 16 00095 g004
Pharmaceuticals 16 00095 g005 550
Figure 5. RMSD (backbone) of the APO structures and AChE and BuChE complexed with ZINC390718 during molecular dynamics of 100 ns (APO forms) and 50 ns (complexed forms). (a) AChE APO form; (b) AChE complexed with ZINC390718; (c) BuChE APO form; (d) BuChE complexed with ZINC390718.
Figure 5. RMSD (backbone) of the APO structures and AChE and BuChE complexed with ZINC390718 during molecular dynamics of 100 ns (APO forms) and 50 ns (complexed forms). (a) AChE APO form; (b) AChE complexed with ZINC390718; (c) BuChE APO form; (d) BuChE complexed with ZINC390718.
Pharmaceuticals 16 00095 g005
Pharmaceuticals 16 00095 g006 550
Figure 6. RMSF (Å) of backbone Cα atoms of the APO structures and complexes of AChE (A) and BuChE (B) with ZINC390718 during the production phase (AChE apo: 15 to 100 ns; AChE complex: 15 to 50 ns) (BuChE apo: 20 to 95 ns; BuChE complex: 5 to 50 ns). Highlights in blue correspond to regions of the active site.
Figure 6. RMSF (Å) of backbone Cα atoms of the APO structures and complexes of AChE (A) and BuChE (B) with ZINC390718 during the production phase (AChE apo: 15 to 100 ns; AChE complex: 15 to 50 ns) (BuChE apo: 20 to 95 ns; BuChE complex: 5 to 50 ns). Highlights in blue correspond to regions of the active site.
Pharmaceuticals 16 00095 g006
Pharmaceuticals 16 00095 g007 550
Figure 7. ZINC04541733 interactions at AChE and BuChE binding sites obtained from the representative structure of the MD simulation. White: ligand carbons; blue: nitrogen; red: oxygen; white sphere: aromatic center; solid blue line: hydrogen bond; gray dashed line: hydrophobic interaction.
Figure 7. ZINC04541733 interactions at AChE and BuChE binding sites obtained from the representative structure of the MD simulation. White: ligand carbons; blue: nitrogen; red: oxygen; white sphere: aromatic center; solid blue line: hydrogen bond; gray dashed line: hydrophobic interaction.
Pharmaceuticals 16 00095 g007
Table
Table 1. Permanence time of hydrogen interactions (Hbond) in the active sites of AChE and BuChE during the production phase and identification of the pairs involved.
Table 1. Permanence time of hydrogen interactions (Hbond) in the active sites of AChE and BuChE during the production phase and identification of the pairs involved.
AChE
ResidueProteinZINC390718Permanence (%)Total
Tyr124OHN116.5239.88
OHN223.36
Gln291NE2N30.5710.83
ON410.26
Ser293NN30.2855.68
ON455.84
Phe295NO186.6186.61
Arg296NO121.3721.37
Tyr341OHN14.2710.53
OHN20.85
ON45.41
BuChE
Asp70OD1N437.9275.17
OD2N437.25
Glu197OE1N212.4215.08
OE2N22.66
Tyr332OHN42.8844.12
OHO10.22
OHN12.66
OHN238.36
Lis339NZN333.0433.04
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de Almeida, R.B.M.; Barbosa, D.B.; do Bomfim, M.R.; Amparo, J.A.O.; Andrade, B.S.; Costa, S.L.; Campos, J.M.; Cruz, J.N.; Santos, C.B.R.; Leite, F.H.A.; et al. Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies. Pharmaceuticals 2023, 16, 95. https://doi.org/10.3390/ph16010095

AMA Style

de Almeida RBM, Barbosa DB, do Bomfim MR, Amparo JAO, Andrade BS, Costa SL, Campos JM, Cruz JN, Santos CBR, Leite FHA, et al. Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies. Pharmaceuticals. 2023; 16(1):95. https://doi.org/10.3390/ph16010095

Chicago/Turabian Style

de Almeida, Raquel B. M., Deyse B. Barbosa, Mayra R. do Bomfim, Jéssika A. O. Amparo, Bruno S. Andrade, Silvia L. Costa, Joaquín M. Campos, Jorddy N. Cruz, Cleydson B. R. Santos, Franco H. A. Leite, and et al. 2023. "Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies" Pharmaceuticals 16, no. 1: 95. https://doi.org/10.3390/ph16010095

APA Style

de Almeida, R. B. M., Barbosa, D. B., do Bomfim, M. R., Amparo, J. A. O., Andrade, B. S., Costa, S. L., Campos, J. M., Cruz, J. N., Santos, C. B. R., Leite, F. H. A., & Botura, M. B. (2023). Identification of a Novel Dual Inhibitor of Acetylcholinesterase and Butyrylcholinesterase: In Vitro and In Silico Studies. Pharmaceuticals, 16(1), 95. https://doi.org/10.3390/ph16010095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop