Next Article in Journal
Structure-Based Virtual Screening towards the Discovery of Novel ULK1 Inhibitors with Anti-HCC Activities
Next Article in Special Issue
Identification of Potential Anti-Neuroinflammatory Inhibitors from Antarctic Fungal Strain Aspergillus sp. SF-7402 via Regulating the NF-κB Signaling Pathway in Microglia
Previous Article in Journal
Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations
Previous Article in Special Issue
Molecular Modeling and Experimental Evaluation of Non-Chiral Components of Bergamot Essential Oil with Inhibitory Activity against Human Monoamine Oxidases
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 9
10.3390/molecules27092626
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

In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease

by
Tijana Bojić
1,*,
Milan Sencanski
2,*,
Vladimir Perovic
2,
Jelena Milicevic
2 and
Sanja Glisic
2
1
Laboratory of Radiobiology and Molecular Genetics-080, Institute of Nuclear Sciences Vinca, National Institute of the Republic of Serbia, University of Belgrade, P.O. Box 522, 11000 Belgrade, Serbia
2
Laboratory for Bioinformatics and Computational Chemistry, Institute of Nuclear Sciences Vinca, National Institute of the Republic of Serbia, University of Belgrade, P.O. Box 522, 11000 Belgrade, Serbia
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(9), 2626; https://doi.org/10.3390/molecules27092626
Submission received: 31 January 2022 / Revised: 31 March 2022 / Accepted: 12 April 2022 / Published: 19 April 2022
(This article belongs to the Special Issue Therapeutic Agents for Neurodegenerative Disorders)
Browse Figures
Versions Notes

Abstract

:
Alzheimer’s disease (AD), a devastating neurodegenerative disease, is the focus of pharmacological research. One of the targets that attract the most attention for the potential therapy of AD is the serotonin 5HT6 receptor, which is the receptor situated exclusively in CNS on glutamatergic and GABAergic neurons. The neurochemical impact of this receptor supports the hypothesis about its role in cognitive, learning, and memory systems, which are of critical importance for AD. Natural products are a promising source of novel bioactive compounds with potential therapeutic potential as a 5HT6 receptor antagonist in the treatment of AD dementia. The ZINC—natural product database was in silico screened in order to find the candidate antagonists of 5-HT6 receptor against AD. A virtual screening protocol that includes both short-and long-range interactions between interacting molecules was employed. First, the EIIP/AQVN filter was applied for in silico screening of the ZINC database followed by 3D QSAR and molecular docking. Ten best candidate compounds were selected from the ZINC Natural Product database as potential 5HT6 Receptor antagonists and were proposed for further evaluation. The best candidate was evaluated by molecular dynamics simulations and free energy calculations.

1. Introduction

Alzheimer’s disease (AD), a neurodegenerative disease of devastating character, strikes over 40 million people worldwide [1]. It is the progressive and most common form of dementia characterized by cholinergic dysfunction, beta-amyloid plaques, and neurofibrillary tangles clustered in brain structures critical for memory and learning processes. A growing number of world elderly population makes this disease reaching epidemic proportions, causing serious human, social, and economic problems [2]. For that reason, AD is in the focus of pharmacological research. Two FDA approved drug classes for the symptomatic treatment of AD are acetylcholinesterase inhibitors (AChEIs) (donepezil, rivastigmine, and galantamine) and NMDA antagonist (memantine). Both classes of these drugs are administered individually or in combination and result in limited treatment success, partial or no response, high intersubject variability, often poor tolerability, and high costs. The side effects cover the broad palette of neurological and psychiatrical manifestations, gastrointestinal, dermatological and allergic disorders and electrolyte imbalance [3]. All these adverse and side effects raise the necessity for new symptomatic drugs that would have better efficacy-safety profile. One of the receptors that attracts the most attention for the potential individual or adjuvant therapy in AD is 5HT6 receptor [4,5,6,7]. 5HT6 receptor is serotonin postsynaptic G-protein metabotropic receptor coupled with the adenylyl cyclase second messenger system [8], which also activates ERK1/2 via Fyn-dependent pathway [2]. 5HT6 receptor is situated exclusively in CNS, on glutamatergic (hippocampus and throughout cerebral cortex), and GABAergic neurons (striatum, nc. Accumbens and olfactory bulb) [1,9,10]. Both neurotransmitter systems and neuroanatomical structures are constituents of memory and learning brain networks [9]. (Figure 1, adapted from [11]). The precise mechanisms of action of 5HT6 receptor antagonists are still unknown. The procognitive effect of 5HT6 receptor antagonists is probably due to the elevation of acetylcholine and glutamate concentrations in the brain by the deinhibition of GABA neurons [10,12]. Unexpected findings in which both 5HT6 receptor antagonists [13] and agonists [14] enhance the cognition in rodents suggest that 5HT6 receptors activate unknown biochemical pathways or act on distinct neuronal populations, inducing paradoxical procognitive effects [9].
On one side, the fact that (a) 5HT6 receptors regulate the balance between excitatory (glutamate) and inhibitory (GABA) signaling [15] and (b) they modulate synaptic plasticity and long-term potentiation in specific memory and learning regions [16] provides rational arguments for pharmacological targeting of the 5HT6 receptor. The brain tissue-specific distribution of this receptor is also the promising feature that 5HT6 pharmacomodulation will produce less peripheral adverse and side effects [1]. In a review by Khoury et al. [9], an overview was conducted of the current drug therapy for AD with recent clinical trials of the 5-HT6 receptor, with an emphasis on the importance of a better understanding of the mechanism of action of the 5-HT6 receptor and its ligands as partial or inverse agonists against AD [9].
Despite years of study efforts to develop therapies, there is no effective therapy available to cure AD or that considerably inhibits the progression of AD symptoms. Since 2003, although there are more than four hundred trials, no new drugs were approved. Most Phase II clinical trials do not succeed in Phase III often due to serious adverse effects or a lack of efficacy [17]. Eight agents listed in phase III in 2017, including the 5-HT6 inhibitors idalopirdine and intepirdine, failed in clinical trials as of 30 January 2018 [18].
Promising results in phase 2b study using intepirdine as adjunctive therapy to donepezil for AD dementia initiated its further examination in a phase 3 study: MINDSET [19]. However, it failed to achieve the assumed primary outcome as there were no statistically significant differences between intepirdine and placebo on the co-primary endpoints ADAS-Cog and ADCS-ADL. Therefore, the study was terminated as the main trial did not reach its primary outcome [20,21].
Only SUVN-502 achieved significant changes from a baseline in ADAS-Cog, especially in combination with donepezil and memantine [1,9,19,22].
SUVN-502 is a third 5-HT6 receptor antagonist that has been investigated for its safety and effectiveness in the phase IIa trial and in the background of donepezil and memantine [23] in patients with moderate AD. The SUVN-502 treatment, compared with placebo, improved behavioral symptoms as measured by neuropsychiatric inventory scores and reduced statistically significant the agitation/aggression scores. It was suggested that SUVN-502 might have a potential utility in the management of behavioral and psychological symptoms of dementia in AD [24]. In a press release of 16 August 2021, Suven Life sciences announced a Phase 3 Clinical Trial of NCE SUVN-502 (Masupirdine) to treat agitation and aggression in AD dementias. The enrollment of patients was expected from September 2021. This multi-center study is expected to be completed in about 36 months, with the result expected by the end of 2024 [25].
The lack of efficacy of many highly selective and specific drug candidates in treating many central nervous system diseases led to the development of multi-targeted drugs, among which include the multi-targeted drug candidate, AVN-101. AVN-101 is a very potent 5-HT7 receptor antagonist with slightly lower potency toward 5-HT6 [26].
AVN-101 is a potent 5-HT7 receptor antagonist that was well tolerated when taken orally at doses of up to 20 mg daily in the Phase I clinical study. AVN-101 was tested in clinical trials for Alzheimer’s disease [26].
Idalopirdine is a 5HT6 receptor antagonist (Saegis Pharmaceuticals) and was at first developed as a treatment for ‘Cognitive Impairment Associated with Schizophrenia’ (CIAS) [27]. After showing a positive effect on cognition, the drug was entered into clinical trials for use in Alzheimer’s dementia [1].
AVN-211 is a highly selective 5-HT6 receptor antagonist. Avineuro Pharmaceuticals has expressed intentions to start clinical trials of AVN-211 for Alzheimer’s disease in 2015 [28].
SGS-518 is also a potent and selective 5HT6 antagonist that is introduced by Piñeiro-Núñez et al. and Filla et al. [29,30].
The more recent clinical studies of 5HT6R antagonosts are mostly in ongoing or completed phase 1 clinical studies (HEC30654AcOH, ClinicalTrials.gov Identifier: NCT03655873, SB-742457, ClinicalTrials.gov Identifier: NCT00551772, PF-05212377, ClinicalTrials.gov Identifier: NCT01213355) [31,32,33].
Novel therapies prevent or delay the onset, slow the progression, or improve the symptoms associated with AD, and new targets are required for development of new drugs. The importance of finding new dementia treatments was also recognized by the G7 meeting, where it was stated that new efficient medications need to be developed by 2025 [34,35].
Natural products have been the primary source of medicines in all cultures. Success linked with screening progress and growing company collections was evident in providing the leads for many drug classes such as kinases and certain classes of the G-protein-coupled receptors (GPCRs) [36]. Natural products are the source of 25–50% of currently marketed drugs. The percentage of these drugs in anti-cancer and anti-infective diseases is even higher, reaching 66% [37]. Currently, two licensed anti-AD drugs (galantamine and rivastigmine) are based on plant-derived natural products, in addition to the fact that natural products targeting tau aggregation in AD pathogenesis are also reported [38].
Computer-aided or in silico designs have played a key role in the development of small molecule therapeutics for over three decades [39]. In silico techniques used in computational pharmacology allow better understanding and prediction of how drugs and compounds affect biological systems, enabling new methods to investigate the drug space and natural product’s “chemical space” [40,41]. The excellent and thorough review of Kelemen et al. of 5-HT6R as a target in AD with all the data regarding computational modeling tools in the field of 5-HT6R drug design, including ligand-based information, pharmacophore models, quantitative structure-activity relationship methods, crystal structure templates, homology models, along with the main goals and challenges in discovery new 5-HT6 antagonists, is very important and useful guide for the design of new 5-HT6R antagonists [42]. In a recent QSAR study, new potential arylsulfonamide-derived 5-HT6 receptor antagonists were proposed as new anti-AD ligands [43].
Virtual screening such as the in silico analog of high throughput screening of large compound databases is an integral part of the drug discovery process, which significantly reduces the time and cost of discovery of a new drug [44]. The availability of databases is essential for this strategy. One of the largest databases is Super Natural II, the first public, freely available web-based database of natural compounds that contains over 326,000 molecules [45]. The TCM Database@Taiwan is the world’s largest non-commercial database that contains compounds from traditional Chinese medicine with more than 20,000 pure compounds allowing free access to 3D compound structures [46]. The authors that have created and curated this database allowed the incorporation of its structures in ZINC. (https://zinc12.docking.org/pbcs/tcmnp, accessed on 31 March 2022). To date, the ZINC database [47,48] is the largest free 3D molecule database. The ZINC database contains over twenty million commercially available molecules. The database is freely available at http://zinc15.docking.org, accessed on 31 March 2022.
The aim of our study was to in silico investigate the ZINC database of natural products to find the potential candidates 5HT6 receptor antagonist for the treatment of AD dementia.

2. Results and Discussion

First, we have applied EIIP/AQVN-based virtual screening of the ZINC Natural Product database. Since there is currently no approved drug antagonist for the 5HT6 receptor, the selected learning set consisted of 5HT6 receptor antagonists that are now in clinical testing as reported in the literature [42] and clinical candidates from the ChEMBL Target Report Card of Serotonin (5-HT) receptor antagonists [49]. The compounds from the learning set (Table 1, Figure 2) were inside the active domain with AQVN and EIIP values within the intervals of (2.51–2.90) and (0.009–0.09), respectively. The domain occupied by 5HT6 drug candidates was quite wide. In order to increase ligand specificity, we have chosen the two subdomains with drug candidates reported to reach the furthest phase of a clinical trial [42] with EIIP/AQVN values 2.55/0.09 and 2.99/0.009 and the third domain with EIIP/AQVN 2.72/0.055 that occupies serotonin, the natural ligand of 5HT6 receptor. The specified domains with AQVN and EIIP values were selected as a criterion for the selection of compounds representing potential candidate 5HT6 antagonists. By applying the EIIP/AQVN-based virtual screening criterion, 789 candidates were chosen from the ZINC Natural Product database (http://zinc15.docking.org/, accessed on 31 March 2022) [50].
After PCA virtual screening, we selected ten of the most similar compounds to the learning set from each EIIP/AQVN range, which was further subjected to molecular docking. The data on PCA model are available in Supplementary Materials.
Furthermore, after the docking process, we selected ten best candidates (Table 2, Figure 3) based on their docking score, i.e., the best balance between calculated free binding energy and favorable and non-favorable protein-ligand interactions.
The results of the docking process showed that molecules form protein–ligand interactions with the following residues in 5HT6 receptor: Asp 106 (D3.32), as an anchor AA residue, by either salt bridge or hydrogen bonding. Ser 193 (S5.43), Thr 196 (T5.46), Asn 288 (N6.55), and Tyr Y7.42 provide additional hydrogen bonding with the receptor. Residues Trp 102 (W3.28), Val 107 (V3.33), Cys 110(C3.36), Trp 281(W6.48), Phe 284 (F6.51), Phe 285 (F6.52), Asp 303 (D7.36), Trp 307 (W7.40), and L157 from ECL2 form hydrophobic or aromatic interactions with a ligand molecule. Those residues were earlier reported as necessary for antagonistic behaviors [51,52,53,54,55,56,57]. Moreover, other reported interaction residues include Asn 86 (N2.64) [58], Ser185 (S5.36), Thr103 (T3.29) [59], Asp 303 (D7.36), Phe 188 (F5.38), and Ala 157 (A4.56) [60]. As a control, we docked the literature antagonists (Figure 1, Table 1). The 2D binding patterns obtained by GOLD are presented on Figures S1–S6 in Supplementary Materials. Figure 4 and Figure 5 present docked conformations of the compounds ZINC00756618 and ZINC20762773, respectively.
The first and second compounds are the best candidates, and they belong to the AQVN descriptor range 2.55. The third compound is the best among those belonging to the AQVN descriptor range 2.77. As mentioned, the anchor AA residue D3.32 forms either hydrogen bonds or salt bridges with protonated or charged nitrogen atoms, while other residues form appropriate interactions depending on the ligand structure. By comparing the GOLD scores (dimensionless quantity) of our candidates to the known ligands from literature (Table 2), they fall in the almost same range, (32.08–47.13) for the candidates and (32.3–42.36) for the literature compounds. Autodock Vina results confirmed the conformations obtained from GOLD (Figures S7 and S8, Supplementary Material) and binding energies are in good correlation except for a few discrepancies. Molecular dynamics simulation of the 5HT6 best candidate complex showed the stability of the receptor, reaching RMSD values around 4 Å after 12ns (Figure 6). In addition, the stabilities of ligand’s conformation in the binding site (Figure S9) and salt bridge between protonated nitrogen in ligand and D3.32 of the receptor (Figure S10) are presented. For the calculation of binding free energy, we chose the FEP thermodynamic integration approach due to the impossibility to perform metadynamics simulation, as the extracellular loop of the receptor hindered ligand’s pulling out of the binding site. After double annihilation consisting of forward-backward FEP for both bound and unbound ligand states, using SOS energy averaging in VMD, we calculated the binding free energy of the best candidate compound (Figure 7).
ΔΔG = ΔGunbound − ΔGbound = 45.06–29.79 = 15.27 kcal/mol, i.e., −15.27 kcal/mol.
Taking into account the size of the system and error, this is an acceptable correlation with binding free energy obtained from Autodock Vina (−13.13 kcal/mol), which indicates the nanomolar activity of our compound (calculated from relation ΔG = RTlnK, where R = 8.314 J mol−1 K−1, T = 298 K).
The drug-likeness of a compound was assessed according to Lipinski’s Rule of Five [61], which considers molecular weight (<500 Da), the number of hydrogen-bond acceptors (≤10) and donors (≤5), and octanol/water partition coefficient (≤5), and Jorgensen’s rule of three [62], which regards logS (>−5.7), PCaco (>22 nm/s), and primary metabolites (PM) (<7). The violations of these rules are essential for the optimization of biologically active compounds and should not be more than 1.
Table 3 presents predicted the ADMET properties of candidate compounds with the following parameters: molecular weight (MW), number of rotatable bonds (RB), dipole moment (DM), molecular volume (MV), number of hydrogen donors (DHB), number of hydrogen acceptors (AHB), polar surface area (PSA), octanol/water partition coefficient (log P), aqueous solubility (log S), apparent Caco-2 cell permeability (PCaco), number of likely primer metabolic reactions (PM), percentage of human oral absorption (%HOA), violations of rules of three (VRT) and five (VRF), blood–brain barrier permeability parameter (QPlogBB) [62], and Multi-Parameter Optimisation score for drugs targeted at the Central Nervous System (CNS MPO) [63]. ADME parameters are presented along with the violations of Lipinski’s and Jorgensen’s rules. All compounds follow both rules, causing no more than one violation. Thus, according to the ADMET properties, all compounds may possess a good pharmacokinetic profile and may be a good candidate for the treatment of AD dementia. Table 4 presents QPlogBB and CNS MPO values for the literature antagonists for comparison.
Since 1993, the year of the discovery of the 5HT6 receptor and the following period of experiments on 5HT6 receptor antagonists, the 5HT6 receptor deserves special attention in the treatment of cognitive deficits [2,10]. In addition to its procognitive effect, the importance of 5HT6 targeting in AD is emphasized by the fact that this receptor is also involved in affective disorders, anxiety, depression, epilepsy, and obesity, pathologies that often accompany and complicate the clinical status of Alzheimer patients. Even though a substantial number of compounds with antagonist effect on 5HT6 receptors was successful in phase I and II of clinical trials (idalopirdine, intepirdine, latrepirdine, and SAM-531), they did not succeed in Phase III, for various reasons, often owing to serious adverse effects or lack of therapeutic efficacy [17]. As a reason why no 5-HT6R ligand has progressed through all clinical stages and the reason for failure in clinical trials is that candidates are unable to cross the blood–brain barrier [43]. An example of the successful drug on the way to be FDA approved is SUVN-502 [9,12,22]. The triple combination of SUVN-502 with donepezil and memantine synergistically increased the acetylcholine levels in the ventral hippocampus [12]. It is important to stress that SUVN-502 is the closest to EIIP/AQVN values to serotonin, the natural ligand of the 5HT6 receptor. Recently, new benzimidazole-based compounds have been developed as a new class of 5-HT6R antagonists with favorable in vivo properties. Supporting the fact that 5HT6 receptor antagonists represent a step toward the identification of new cognitive enhancers for the treatment of AD, the development of new 5-HT6R antagonists is still required to validate these molecules as a drug class against AD [61].
For in silico approaches or computer-aided drug, the design has provided a large number of opportunities for the identification of novel potential lead compounds against many diseases, including AD, saving the cost and time of performing clinical trials and increasing the success rate of new drugs [64]. These approaches comprise the conformational modeling of small molecules, macromolecules and their docked complexes, calculation descriptors, and the independent molecular properties, followed by model constructions, the design of novel molecules, and the prediction of the potential of bioactive compounds [65].
Natural products are promising potential sources for novel compounds that can meet the criteria of better efficacy and better tolerability in the treatment of AD. In silico approaches could be of significant help as there remains a significant need for new 5-HT6R leads and drug candidates with the potential to improve cognitive deficits in AD.
It was shown earlier for molecular targets in various pathologies that small molecules with similar AQVN and EIIP values interact with the common therapeutic target [66]. This has led to the definition of the EIIP/AQVN criteria for the virtual screening of molecular libraries for compounds with similar treatment properties [67]. Prior studies have verified the EIIP/AQVN approach to be useful in identifying inhibitors against various viral and chronic disease targets [68,69] that were also confirmed experimentally [70,71].
The EIIP/AQVN approach in the previous drug database screening revealed potential candidates for the pharmacotherapy of neurocardiovascular diseases and chronic tinnitus [72,73]. In this study, an in silico EIIP/AQVN-based filter was applied as the first for the selection of 5HT6 receptor antagonists during the virtual screening process. As a result, 789 candidates were chosen from the ZINC—natural product database. This step was followed by PCA analysis and, finally, molecular docking. The screening identified ten best candidates as the most promising 5HT6 receptor antagonist. Our selected candidates show good CNS drug ability properties. There is no known activity for the ten virtual hits in the ZINC database or other databases linked to ZINC.
NPBS (Natural Products & Biological Sources) database is a chemical database connecting natural products and biological sources that is manually curated from the literature [74]. The biological sources cover diverse plant, bacterial, fungal, and marine species. In the NPBS database, we did not find a biological source for any of the ten compounds. Furthermore, we investigated the natural origin of the ten virtual hits in the KNApSAcK Family Database [75], which contains data on compounds linked to information on the plant family and species from which the compound was isolated. However, we did not find any data about plant families or species for the ten proposed hit compounds. The proposed candidates will be evaluated in further experimental testing.

3. Materials and Methods

A library of natural products that contain 1 × 105 compounds was downloaded from the ZINC database in the processed format [48]. The EIIP/AQVN approach was employed for in silico screening of the ZINC—natural product database and then proceeded by PCA analysis and molecular docking for the identification of potential 5HT6 receptor antagonist.

3.1. EIIP/AQVN

The EIIP for organic molecules can be determined by the following simple equation derived from the “general model pseudopotential [67]:
EIIP = 0.25 × Z* × sin (1.04 × π × Z*)/2π,
where Z* is the average quasi valence number (AQVN) determined by the following:
Z* = ∑m (ni × Zi/N),
where Zi is the valence number of the ith atomic component, ni is the number of atoms of the ith component, m is the number of atomic components in the molecule, and N is the total number of atoms. EIIP values calculated according to Equations (1) and (2) are expressed in Rydberg units (Ry).

3.2. Zinc Ligand Database and Preparation

The structures of all candidate ligands were taken from the ZINC database [48] in sdf format, which is already geometrically optimized and protonated at physiological pH. The total number of ligands from all three EIIP/AQVN ranges was 789.
Ligands for the learning sets were downloaded from the CHEMBL database [76], including receptor ID CHEMBL3371 containing 3292 active compounds and containing Ki values in the range from 0.04 nM to 30 μM. The set was randomly sampled, resulting in 100 compound learning set for virtual screening. Molecules were converted from SMILES to sdf format.

3.3. PCA Analysis of Virtual Screening

The obtained learning set from the previous step was imported, and molecules were aligned towards the principal moment of inertia axis and protonated on physiological pH (7.4). The candidate molecules were also imported for the next virtual screening step. MIF descriptors were computed for all molecules. The computation method for descriptor generation was GRID with step 0.5. Applied probes (mapped regions of molecule surface) were DRY (hydrophobic interactions), O (hydrogen bond acceptor), N1 (hydrogen bond donor), and TIP (molecular shape descriptor). The discretization method was AMANDA with a scale factor of 0.55. The encoding method was MACC2, and the weights were the following: DRY: −0.5; O: −2.6; N1: −4.2; and TIP: −0.75. Based on the learning set, the PCA model was built, consisting of five principal components. The explained variance was 64.84%. The virtual screening of compound candidates was carried by the use of the shortest centroid distance. The calculation was carried in Pentacle software version 1.06 for Linux [77,78,79]. From each of the three candidate groups, ten of the most similar molecules were selected and further submitted to molecular docking.

3.4. Receptor Preparation

The homology model of the human 5HT6 receptor was downloaded from the GPCRdb database [80] in its intermediate form. The receptor was protonated at physiological pH 7.4. All receptor amino acid residues (AA) were numbered in Ballesteros-Weinstein nomenclature [81].

3.5. Molecular Docking

Chemscore kinase protocol was used for the molecular docking process [82,83]. The binding site for molecular docking was chosen as a 15 Å radius sphere, with the center in the carbon atom of the D3.32 carboxyl group. The (x,y,z) center of the grid box was (−25.9, 5.82, 18.04). As a scoring function was used, the ChemScore. Search Efficiency was set to 200%. The best ten candidates, based on score and specific receptor–ligand interaction criteria, were selected. The calculations were carried in GOLD software [84]. Molecular docking was also carried in newly developed Autodock Vina 1.2.0 [85] using Autodock4 forcefield. Receptor maps were prepared using Autogrid 4. Exhaustiveness was set to 250. Both receptor and ligands were prepared in ADT Tools 1.5.6 [86]
For intermolecular interaction identification, the hydrogen bonds were recognized if a maximum distance between the donor (D) and a hydrogen atom (H), D · · · H, was 3.4 Å and had an angle D-H-A between 90° and 180°. For a salt bridge, a maximum distance between the charged atoms was 4 Å. For π–alkyl interactions, the maximum distance between the centroid of the aromatic ring and the C-atom of the alkyl group was 4 Å and angle 45°. For π–π interactions, the maximum distance between centroids was 6 Å, and the values of theta and gamma angles were 50° and 35° for stacked and 30° and 55° for T-shaped conformation, respectively. For π–cation interactions, the maximum distance between cation and centroid was 4 Å and angle 40°.

3.6. Molecular Dynamics Simulations and Binding Free Energy Calculations

The 5HT6-ligand complex was inserted into a 70 × 70 Å 2-oleoyl-1-palmitoyl-sn-glyecro-3-phosphocholine (POPC) lipid bilayer, created in VMD 1.9.3 [87]. A 20 Å water layer was added from the both sides of the z-axis. Bad-contact water molecules were removed from the lipid membrane bilayer using the appropriate tcl script. Additionally, the system was neutralized with 0.15 M NaCl, resulting in a 61,203 atoms ensemble. The system was subjected to a 10,000-step energy minimization, 250 NVE ps equilibration, and 40 ns NPT MD production. Pressure and temperature were set to 1 bar and 310 K, respectively, using a Berendsen thermostat, and the applied integration step was 1 fs. In all simulations, periodic boundary conditions with particle-mesh Ewald calculations were implemented. The cut-off was set to 12 Å. A CHARMM36 force field was used [88].
Binding free energy was calculated using Free Energy Perturbation (FEP) approach. Thermodynamic integration was carried for bound and unbound ligand conformation (double annihilation), both forward and backward, using soft-core potential approach. Alchemical Van der Waals shift coefficient was set to 6.0. Each λ step of simulation consisted of 100 ps alchemical equilibration, followed by 400 ps of production. The simulation was set from λ 0 to 1, with step 0.02, in 50 cycles. The simulation took 137,439 s for bound and 27,940 s for unbound form @ Intel i7-10700KF @3.80 GHz and NVIDIA GeForce 1080 Ti.
All simulations were carried in NAMD 3.0 alpha, which is able to perform CUDA accelerated FEP calculations [89]. Results were analyzed in VMD 1.9.3 Analyse FEP simulation extension with an SOS (simple-overlap sampling) estimator.

3.7. ADME/Toxicity Prediction

ADMET parameters of the natural compounds (1–10) were calculated with the help of QikProp software running in normal mode [90,91]. CNS MPO values were calculated using MarvinSketch 22.9 [92].

4. Conclusions

Natural compounds present an important treatment source of novel, potent, and selective candidates for treating the cognitive deficit due to AD. The ZINC Natural Product database was virtually screened for candidate antagonists of 5HT6 receptors. The ten most promising candidate compounds were selected as potential candidates. The results of our study will be confirmed experimentally in further research studies that will show whether these compounds represent suitable preliminary candidates to alleviate the heavy burden of Alzheimer’s disease to patients, families, and the public health system.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27092626/s1, Figure S1. Binding pattern of Intepirdine in the active site of 5HT6 receptor; Figure S2. Binding pattern of Cerlapirdine in the active site of 5HT6 receptor; Figure S3. Binding pattern of Idalopirdine in the active site of 5HT6 receptor; Figure S4. Binding pattern of SG518 in the active site of 5HT6 receptor; Figure S5. Binding pattern of Masupirdine in the active site of 5HT6 receptor; Figure S6. Binding pattern of AVN-211 in the active site of 5HT6 receptor; Figure S7. Binding modes of compound ZINC00756618, obtained by GOLD (green) and Autodock Vina (cyan); Figure S8. Binding modes of compound ZINC2076277, obtained by GOLD (green) and Autodock Vina (cyan); Figure S9. RMSD of ZINC00756618 in the binding site of 5HT6 receptor during 40ns production phase of MD simulation; Figure S10. Salt bridge in 5HT6—ZINC00756618 complex during 40ns production phase simulation, presented as interactions of NH+ (ligand) and oxygens from COO- group from D3.32. Table S1: List of compounds from AQVN group values 2.55 with centroid distance from the learning set; Table S2: List of compounds from AQVN group values 2.72 with centroid distance from the learning set; Table S3: List of compounds from AQVN group values 2.99 with centroid distance the from learning set; Table S4: 5HT6 compounds learning set from ChEMBL database.

Author Contributions

Conceptualization, T.B., M.S., and S.G.; methodology, M.S., S.G., V.P., and J.M.; software, V.P.; validation, S.G., M.S., and J.M.; formal analysis, S.G.; investigation, T.B., S.G., and M.S.; resources, S.G. and M.S.; data curation, S.G.; writing—original draft preparation, T.B., S.G., M.S., and J.M.; writing—review and editing, T.B., S.G., M.S., and J.M.; visualization, M.S.; supervision, S.G.; project administration, T.B., S.G., and M.S.; funding acquisition, S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Andrews, M.; Tousi, B.; Sabbagh, M.N. 5HT6 Antagonists in the Treatment of Alzheimer’s Dementia: Current Progress. Neurol. Ther. 2018, 7, 51–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ferrero, H.; Solas, M.; Francis, P.T.; Ramirez, M.J. Serotonin 5-HT(6) Receptor Antagonists in Alzheimer’s Disease: Therapeutic Rationale and Current Development Status. CNS Drugs 2017, 31, 19–32. [Google Scholar] [CrossRef] [PubMed]
  3. Kanagaratnam, L.; Dramé, M.; Trenque, T.; Oubaya, N.; Nazeyrollas, P.; Novella, J.-L.; Jolly, D.; Mahmoudi, R. Adverse drug reactions in elderly patients with cognitive disorders: A systematic review. Maturitas 2016, 85, 56–63. [Google Scholar] [CrossRef] [PubMed]
  4. Azam, S.; Haque, M.E.; Jakaria, M.; Jo, S.-H.; Kim, I.-S.; Choi, D.-K. G-Protein-Coupled Receptors in CNS: A Potential Therapeutic Target for Intervention in Neurodegenerative Disorders and Associated Cognitive Deficits. Cells 2020, 9, 506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sudoł, S.; Cios, A.; Jastrzębska-Więsek, M.; Honkisz-Orzechowska, E.; Mordyl, B.; Wilczyńska-Zawal, N.; Satała, G.; Kucwaj-Brysz, K.; Partyka, A.; Latacz, G.; et al. The Phenoxyalkyltriazine Antagonists for 5-HT6 Receptor with Promising Procognitive and Pharmacokinetic Properties In Vivo in Search for a Novel Therapeutic Approach to Dementia Diseases. Int. J. Mol. Sci. 2021, 22, 10773. [Google Scholar] [CrossRef] [PubMed]
  6. Li, X.; Gao, L.; Liu, J.; Zhang, H.; Chen, H.; Yang, L.; Wu, M.; Li, C.; Zhu, X.; Ding, Y.; et al. Safety, Tolerability and Pharmacokinetics of the Serotonin 5-HT6 Receptor Antagonist, HEC30654, in Healthy Chinese Subjects. Front. Pharmacol. 2021, 12, 726536. [Google Scholar] [CrossRef] [PubMed]
  7. Di Giovanni, G.; De Deurwaerdère, P. Serotonin research: Crossing scales and boundaries. Neuropharmacology 2020, 181, 108340. [Google Scholar] [CrossRef]
  8. Grimaldi, B.; Bonnin, A.; Fillion, M.P.; Ruat, M.; Traiffort, E.; Fillion, G. Characterization of 5-ht6 receptor and expression of 5-ht6 mRNA in the rat brain during ontogenetic development. Naunyn Schmiedebergs Arch. Pharmacol. 1998, 357, 393–400. [Google Scholar] [CrossRef] [PubMed]
  9. Khoury, R.; Grysman, N.; Gold, J.; Patel, K.; Grossberg, G.T. The role of 5 HT6-receptor antagonists in Alzheimer’s disease: An update. Expert Opin. Investig. Drugs 2018, 27, 523–533. [Google Scholar] [CrossRef]
  10. Ramírez, M.J. 5-HT6 receptors and Alzheimer’s disease. Alzheimer’s Res. Ther. 2013, 5, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Codony, X.; Vela, J.M.; Ramírez, M.J. 5-HT(6) receptor and cognition. Curr. Opin. Pharmacol. 2011, 11, 94–100. [Google Scholar] [CrossRef] [PubMed]
  12. Nirogi, R.; Shinde, A.; Kambhampati, R.S.; Mohammed, A.R.; Saraf, S.K.; Badange, R.K.; Bandyala, T.R.; Bhatta, V.; Bojja, K.; Reballi, V.; et al. Discovery and Development of 1-[(2-Bromophenyl)sulfonyl]-5-methoxy-3-[(4-methyl-1-piperazinyl)methyl]-1H-indole Dimesylate Monohydrate (SUVN-502): A Novel, Potent, Selective and Orally Active Serotonin 6 (5-HT(6)) Receptor Antagonist for Potential Treatment of Alzheimer’s Disease. J. Med. Chem. 2017, 60, 1843–1859. [Google Scholar] [CrossRef] [PubMed]
  13. Bourson, A.; Borroni, E.; Austin, R.H.; Monsma, F.J.J.; Sleight, A.J. Determination of the role of the 5-ht6 receptor in the rat brain: A study using antisense oligonucleotides. J. Pharmacol. Exp. Ther. 1995, 274, 173–180. [Google Scholar] [PubMed]
  14. Nikiforuk, A. The procognitive effects of 5-HT6 receptor ligands in animal models of schizophrenia. Rev. Neurosci. 2014, 25, 367–382. [Google Scholar] [CrossRef] [PubMed]
  15. De Jong, I.E.M.; Mørk, A. Antagonism of the 5-HT(6) receptor—Preclinical rationale for the treatment of Alzheimer’s disease. Neuropharmacology 2017, 125, 50–63. [Google Scholar] [CrossRef] [PubMed]
  16. West, P.J.; Marcy, V.R.; Marino, M.J.; Schaffhauser, H. Activation of the 5-HT(6) receptor attenuates long-term potentiation and facilitates GABAergic neurotransmission in rat hippocampus. Neuroscience 2009, 164, 692–701. [Google Scholar] [CrossRef] [PubMed]
  17. Cummings, J.L.; Morstorf, T.; Zhong, K. Alzheimer’s disease drug-development pipeline: Few candidates, frequent failures. Alzheimer’s Res. Ther. 2014, 6, 37. [Google Scholar] [CrossRef] [Green Version]
  18. Cummings, J.; Lee, G.; Ritter, A.; Zhong, K. Alzheimer’s disease drug development pipeline: 2018. Alzheimer’s Dement. 2018, 4, 195–214. [Google Scholar] [CrossRef] [PubMed]
  19. Lang, F.M.; Mo, Y.; Sabbagh, M.; Solomon, P.; Boada, M.; Jones, R.W.; Frisoni, G.B.; Grimmer, T.; Dubois, B.; Harnett, M.; et al. Intepirdine as adjunctive therapy to donepezil for mild-to-moderate Alzheimer’s disease: A randomized, placebo-controlled, phase 3 clinical trial (MINDSET). Alzheimer’s Dement. Transl. Res. Clin. Interv. 2021, 7, e12136. [Google Scholar] [CrossRef] [PubMed]
  20. Kucwaj-Brysz, K.; Baltrukevich, H.; Czarnota, K.; Handzlik, J. Chemical update on the potential for serotonin 5-HT(6) and 5-HT(7) receptor agents in the treatment of Alzheimer’s disease. Bioorganic Med. Chem. Lett. 2021, 49, 128275. [Google Scholar] [CrossRef] [PubMed]
  21. Huang, L.-K.; Chao, S.-P.; Hu, C.-J. Clinical trials of new drugs for Alzheimer disease. J. Biomed. Sci. 2020, 27, 18. [Google Scholar] [CrossRef] [PubMed]
  22. Nirogi, R.; Abraham, R.; Benade, V.; Medapati, R.B.; Jayarajan, P.; Bhyrapuneni, G.; Muddana, N.; Mekala, V.R.; Subramanian, R.; Shinde, A.; et al. SUVN-502, a novel, potent, pure, and orally active 5-HT6 receptor antagonist: Pharmacological, behavioral, and neurochemical characterization. Behav. Pharmacol. 2019, 30, 16–35. [Google Scholar] [CrossRef] [PubMed]
  23. NIH SUVN-502 with Donepezil and Memantine for the Treatment of Moderate Alzheimer’s Disease-Phase 2a Study 2019, Identifier: NCT02580305. Available online: https://clinicaltrials.gov/ct2/show/NCT02580305 (accessed on 31 March 2022).
  24. Nirogi, R.; Shinde, A.; Jayarajan, P.; Goyal, V.; Bhyrapuneni, G.; Subramanian, R.; Jeta, S.; Ravula, J.; Jasti, V. Potential benefits of Masupirdine (SUVN-502) on Behavioral and Psychological symptoms in patients with moderate Alzheimer’s Disease (5090). Neurology 2020, 94, 5090. [Google Scholar]
  25. Suven Life Sciences Suven Life Sciences Announces Phase 3 Clinical Trial of SUVN-502 (Masupirdine), a 5-HT6 Antagonist for Treatment of Agitation and Aggression in Alzheimer’s Type Dementias. News Release 2021. Available online: http://www.suven.com/pdf/NewRelease16Aug2021.pdf (accessed on 31 March 2022).
  26. Ivachtchenko, A.V.; Lavrovsky, Y.; Okun, I. AVN-101: A Multi-Target Drug Candidate for the Treatment of CNS Disorders. J. Alzheimer’s Dis. 2016, 53, 583–620. [Google Scholar] [CrossRef] [Green Version]
  27. Hesselink, J.M.K. Idalopirdine (LY483518, SGS518, Lu AE 58054) in Alzheimer disease: Never change a winning team and do not build exclusively on surrogates. Lessons Learned from Drug Development Trials. J. Pharmacol. Clin. Res. 2016, 2, 001–004. [Google Scholar] [CrossRef]
  28. Ivachtchenko, A.V.; Lavrovsky, Y.; Ivanenkov, Y.A. AVN-211, Novel and Highly Selective 5-HT6 Receptor Small Molecule Antagonist, for the Treatment of Alzheimer’s Disease. Mol. Pharm. 2016, 13, 945–963. [Google Scholar] [CrossRef]
  29. Piñeiro-Núñez, M.M.; Bauzon, D.D.; Bymaster, F.P.; Chen, Z.; Chernet, M.P.; Clay, M.P.; Crile, R.; Delapp, N.W.; Denny, C.P.; Falcone, J.F.; et al. Discovery and SAR studies of 2,6-difluorobenzenesulfonic acid 1-methyl-3-(1-methylpiperidin-4-yl)-1H-indol-5-yl ester, a novel and potent 5-HT6 antagonist for the treatment of cognitive deficit. In Proceedings of the 229th ACS National Meeting, San Diego, CA, USA, 13–17 March 2005. [Google Scholar]
  30. Filla, S.A.; Flaugh, M.E.; Gillig, J.R.; Heinz, L.J.; Krushinski, J.H.J.; Liu, B.; Pineiro-Nunez, M.M.; Schaus, J.M.; Ward, J.S. Benzenesulfonic Acid Indol-5-yl Esters as Antagonists of the 5-ht6 Receptor. Patent WO2002060871A2, 8 August 2002. [Google Scholar]
  31. Sunshine Lake Pharma Co., Ltd. The Study of a Selective 5-ht6 Receptor Antagonist, HEC30654AcOH, in Healthy Subjects. 2019. Available online: https://clinicaltrials.gov/ct2/show/NCT03655873 (accessed on 31 March 2022).
  32. Pfizer. Scopolamine Challenge Study. 2018. Available online: https://clinicaltrials.gov/ct2/show/NCT01213355 (accessed on 31 March 2022).
  33. GlaxoSmithKline. A Study to Assess the Pharmacokinetics of SB-742457 Formulated as a Capsule and a Tablet in Healthy Elderly Volunteers. 2008. Available online: https://clinicaltrials.gov/ct2/show/NCT00551772 (accessed on 31 March 2022).
  34. Cummings, J.; Lee, G.; Mortsdorf, T.; Ritter, A.; Zhong, K. Alzheimer’s disease drug development pipeline: 2017. Alzheimer’s Dement. 2017, 3, 367–384. [Google Scholar] [CrossRef]
  35. Tewari, D.; Stankiewicz, A.M.; Mocan, A.; Sah, A.N.; Tzvetkov, N.T.; Huminiecki, L.; Horbańczuk, J.O.; Atanasov, A.G. Ethnopharmacological Approaches for Dementia Therapy and Significance of Natural Products and Herbal Drugs. Front. Aging Neurosci. 2018, 10, 3. [Google Scholar] [CrossRef] [Green Version]
  36. Barker, A.; Kettle, J.G.; Nowak, T.; Pease, J.E. Expanding medicinal chemistry space. Drug Discov. Today 2013, 18, 298–304. [Google Scholar] [CrossRef]
  37. Kingston, D.G.I. Modern Natural Products Drug Discovery and Its Relevance to Biodiversity Conservation. J. Nat. Prod. 2011, 74, 496–511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Nisbet, R.M.; Polanco, J.-C.; Ittner, L.M.; Götz, J. Tau aggregation and its interplay with amyloid-β. Acta Neuropathol. 2015, 129, 207–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Sliwoski, G.; Kothiwale, S.; Meiler, J.; Lowe, E.W., Jr. Computational methods in drug discovery. Pharmacol. Rev. 2013, 66, 334–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Hodos, R.A.; Kidd, B.A.; Shameer, K.; Readhead, B.P.; Dudley, J.T. In silico methods for drug repurposing and pharmacology. Wiley Interdiscip. Rev. Syst. Biol. Med. 2016, 8, 186–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Glisic, S.; Sencanski, M.; Perovic, V.; Stevanovic, S.; García-Sosa, A.T. Arginase Flavonoid Anti-Leishmanial in Silico Inhibitors Flagged against Anti-Targets. Molecules 2016, 21, 589. [Google Scholar] [CrossRef] [Green Version]
  42. Kelemen, A.A.; Mordalski, S.; Bojarski, A.J.; Keseru, G.M. Computational modeling of drugs for Alzheimer’s disease: Design of serotonin 5-HT 6 antagonists. In Computational Modeling of Drugs against Alzheimer’s Disease; Springer: New York, NY, USA; Humana Press: Totowa, NJ, USA, 2018; pp. 419–461. [Google Scholar]
  43. Da Silva, A.P.; de Angelo, R.M.; de Paula, H.; Honório, K.M.; da Silva, A.B.F. Drug design of new 5-HT6 antagonists: A QSAR study of arylsulfonamide derivatives. Struct. Chem. 2020, 31, 1585–1597. [Google Scholar] [CrossRef]
  44. McInnes, C. Virtual screening strategies in drug discovery. Curr. Opin. Chem. Biol. 2007, 11, 494–502. [Google Scholar] [CrossRef]
  45. Banerjee, P.; Erehman, J.; Gohlke, B.-O.; Wilhelm, T.; Preissner, R.; Dunkel, M. Super Natural II—A database of natural products. Nucleic Acids Res. 2015, 43, D935–D939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chen, C.Y.-C. TCM Database@Taiwan: The World’s Largest Traditional Chinese Medicine Database for Drug Screening In Silico. PLoS ONE 2011, 6, e15939. [Google Scholar] [CrossRef] [Green Version]
  47. Irwin, J.J.; Sterling, T.; Mysinger, M.M.; Bolstad, E.S.; Coleman, R.G. ZINC: A Free Tool to Discover Chemistry for Biology. J. Chem. Inf. Model. 2012, 52, 1757–1768. [Google Scholar] [CrossRef]
  48. Irwin, J.J.; Shoichet, B.K. ZINC—A Free Database of Commercially Available Compounds for Virtual Screening. J. Chem. Inf. Model. 2005, 45, 177–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. CHembl CHEMBL3371. Available online: https://www.ebi.ac.uk/chembl/target_report_card/CHEMBL3371/ (accessed on 31 March 2022).
  50. Sterling, T.; Irwin, J.J. ZINC 15—Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef] [PubMed]
  51. De la Fuente, T.; Martín-Fontecha, M.; Sallander, J.; Benhamú, B.; Campillo, M.; Medina, R.A.; Pellissier, L.P.; Claeysen, S.; Dumuis, A.; Pardo, L.; et al. Benzimidazole derivatives as new serotonin 5-HT6 receptor antagonists. Molecular mechanisms of receptor inactivation. J. Med. Chem. 2010, 53, 1357–1369. [Google Scholar] [CrossRef] [PubMed]
  52. Mierzejewski, P.; Bienkowski, P. Novel Arylsulfonamide Derivatives with 5-HT6/5-HT7 Receptor Antagonism Targeting Behavioral and Psychological Symptoms of Dementia. J. Med. Chem. 2014, 57, 4543–4557. [Google Scholar]
  53. Pullagurla, M.R.; Westkaemper, R.B.; Glennon, R.A. Possible differences in modes of agonist and antagonist binding at human 5-HT6 receptors. Bioorganic Med. Chem. Lett. 2004, 14, 4569–4573. [Google Scholar] [CrossRef] [PubMed]
  54. Gordillo-Cruz, R.E.; Rentería-Gómez, A.; Islas-Jácome, A.; Cortes-García, C.J.; Díaz-Cervantes, E.; Robles, J.; Gámez-Montaño, R. Synthesis of 3-tetrazolylmethyl-azepino[4,5-b]indol-4-ones in two reaction steps: (Ugi-azide/N-acylation/SN2)/free radical cyclization and docking studies to a 5-Ht6 model. Org. Biomol. Chem. 2013, 11, 6470–6476. [Google Scholar] [CrossRef] [PubMed]
  55. Łażewska, D.; Kurczab, R.; Więcek, M.; Kamińska, K.; Satała, G.; Jastrzębska-Więsek, M.; Partyka, A.; Bojarski, A.J.; Wesołowska, A.; Kieć-Kononowicz, K.; et al. The computer-aided discovery of novel family of the 5-HT6 serotonin receptor ligands among derivatives of 4-benzyl-1,3,5-triazine. Eur. J. Med. Chem. 2017, 135, 117–124. [Google Scholar] [CrossRef] [PubMed]
  56. Staroń, J.; Kurczab, R.; Warszycki, D.; Satała, G.; Krawczyk, M.; Bugno, R.; Lenda, T.; Popik, P.; Hogendorf, A.S.; Hogendorf, A.; et al. Virtual screening-driven discovery of dual 5-HT6/5-HT2A receptor ligands with pro-cognitive properties. Eur. J. Med. Chem. 2019, 185. [Google Scholar] [CrossRef]
  57. Dukat, M.; Mosier, P.D.; Kolanos, R.; Roth, B.L.; Glennon, R.A. Binding of serotonin and N1-benzenesulfonyltryptamine-related analogs at human 5-HT6 serotonin receptors: Receptor modeling studies. J. Med. Chem. 2008, 51, 603–611. [Google Scholar] [CrossRef] [PubMed]
  58. Vera, G.; Lagos, C.F.; Almendras, S.; Hebel, D.; Flores, F.; Valle-Corvalán, G.; Pessoa-Mahana, C.D.; Mella-Raipán, J.; Montecinos, R.; Recabarren-Gajardo, G. Extended N-Arylsulfonylindoles as 5-HT6 Receptor Antagonists: Design, Synthesis & Biological Evaluation. Molecules 2016, 21, 1070. [Google Scholar]
  59. Harris, R.N.; Stabler, R.S.; Repke, D.B.; Kress, J.M.; Walker, K.A.; Martin, R.S.; Brothers, J.M.; Ilnicka, M.; Lee, S.W.; Mirzadegan, T. Highly potent, non-basic 5-HT6 ligands. Site mutagenesis evidence for a second binding mode at 5-HT6 for antagonism. Bioorg. Med. Chem. Lett. 2010, 20, 3436–3440. [Google Scholar] [CrossRef] [PubMed]
  60. Marcinkowska, M.; Bucki, A.; Panek, D.; Siwek, A.; Fajkis, N.; Bednarski, M.; Zygmunt, M.; Godyń, J.; Del Rio Valdivieso, A.; Kotańska, M.; et al. Anti-Alzheimer’s multitarget-directed ligands with serotonin 5-HT6 antagonist, butyrylcholinesterase inhibitory, and antioxidant activity. Arch. Pharm. (Weinh.) 2019, 352, 1900041. [Google Scholar] [CrossRef] [PubMed]
  61. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  62. Kelder, J.; Grootenhuis, P.D.; Bayada, D.M.; Delbressine, L.P.; Ploemen, J.P. Polar molecular surface as a dominating determinant for oral absorption and brain penetration of drugs. Pharm. Res. 1999, 16, 1514–1519. [Google Scholar] [CrossRef] [PubMed]
  63. Wager, T.T.; Hou, X.; Verhoest, P.R.; Villalobos, A. Moving beyond Rules: The Development of a Central Nervous System Multiparameter Optimization (CNS MPO) Approach To Enable Alignment of Druglike Properties. ACS Chem. Neurosci. 2010, 1, 435–449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. González-Vera, J.A.; Medina, R.A.; Martín-Fontecha, M.; Gonzalez, A.; de la Fuente, T.; Vázquez-Villa, H.; García-Cárceles, J.; Botta, J.; McCormick, P.J.; Benhamú, B.; et al. A new serotonin 5-HT6 receptor antagonist with procognitive activity–Importance of a halogen bond interaction to stabilize the binding. Sci. Rep. 2017, 7, 41293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Scotti, L.; Scotti, M.T. In Silico Studies Applied to Natural Products with Potential Activity against Alzheimer’s Disease. In Computational Modeling of Drugs against Alzheimer’s Disease; Roy, K., Ed.; Springer: New York, NY, USA, 2018; pp. 513–531. ISBN 978-1-4939-7404-7. [Google Scholar]
  66. Veljkovic, N.; Glisic, S.; Perovic, V.; Veljkovic, V. The role of long-range intermolecular interactions in discovery of new drugs. Expert Opin. Drug Discov. 2011, 6, 1263–1270. [Google Scholar] [CrossRef]
  67. Veljkovic, N.; Glisic, S.; Prljic, J.; Perovic, V.; Veljkovic, V. Simple and general criterion for “in silico” screening of candidate HIV drugs. Curr. Pharm. Biotechnol. 2013, 14, 561–569. [Google Scholar] [CrossRef] [PubMed]
  68. Tintori, C.; Manetti, F.; Veljkovic, N.; Perovic, V.; Vercammen, J.; Hayes, S.; Massa, S.; Witvrouw, M.; Debyser, Z.; Veljkovic, V.; et al. Novel Virtual Screening Protocol Based on the Combined Use of Molecular Modeling and Electron-Ion Interaction Potential Techniques To Design HIV-1 Integrase Inhibitors. J. Chem. Inf. Model. 2007, 47, 1536–1544. [Google Scholar] [CrossRef] [PubMed]
  69. Veljkovic, V.; Goeijenbier, M.; Glisic, S.; Veljkovic, N.; Perovic, V.R.; Sencanski, M.; Branch, D.R.; Paessler, S. In silico analysis suggests repurposing of ibuprofen for prevention and treatment of EBOLA virus disease. F1000Research 2015, 4, 104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Zhao, Y.; Ren, J.; Harlos, K.; Jones, D.M.; Zeltina, A.; Bowden, T.A.; Padilla-Parra, S.; Fry, E.E.; Stuart, D.I. Toremifene interacts with and destabilizes the Ebola virus glycoprotein. Nature 2016, 535, 169–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Paessler, S.; Huang, C.; Sencanski, M.; Veljkovic, N.; Perovic, V.; Glisic, S.; Veljkovic, V. Ibuprofen as a template molecule for drug design against Ebola virus. Front. Biosci. (Landmark Ed.) 2018, 23, 947–953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Bojić, T.; Perović, V.R.; Senćanski, M.; Glišić, S. Identification of Candidate Allosteric Modulators of the M1 Muscarinic Acetylcholine Receptor Which May Improve Vagus Nerve Stimulation in Chronic Tinnitus. Front. Neurosci. 2017, 11, 636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Bojić, T.; Perović, V.R.; Glišić, S. In silico Therapeutics for Neurogenic Hypertension and Vasovagal Syncope. Front. Neurosci. 2016, 9, 520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Xu, T.; Chen, W.; Zhou, J.; Dai, J.; Li, Y.; Zhao, Y. NPBS database: A chemical data resource with relational data between natural products and biological sources. Database 2020, 2020, baaa102. [Google Scholar] [CrossRef] [PubMed]
  75. Afendi, F.M.; Okada, T.; Yamazaki, M.; Hirai-Morita, A.; Nakamura, Y.; Nakamura, K.; Ikeda, S.; Takahashi, H.; Altaf-Ul-Amin, M.; Darusman, L.K.; et al. KNApSAcK Family Databases: Integrated Metabolite–Plant Species Databases for Multifaceted Plant Research. Plant Cell Physiol. 2012, 53, e1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Bento, A.P.; Gaulton, A.; Hersey, A.; Bellis, L.J.; Chambers, J.; Davies, M.; Krüger, F.A.; Light, Y.; Mak, L.; McGlinchey, S.; et al. The ChEMBL bioactivity database: An update. Nucleic Acids Res. 2014, 42, D1083–D1090. [Google Scholar] [CrossRef] [Green Version]
  77. Durán, Á.; Zamora, I.; Pastor, M. Suitability of GRIND-Based Principal Properties for the Description of Molecular Similarity and Ligand-Based Virtual Screening. J. Chem. Inf. Model. 2009, 49, 2129–2138. [Google Scholar] [CrossRef] [Green Version]
  78. Pastor, M.; Cruciani, G.; McLay, I.; Pickett, S.; Clementi, S. GRid-INdependent Descriptors (GRIND):  A Novel Class of Alignment-Independent Three-Dimensional Molecular Descriptors. J. Med. Chem. 2000, 43, 3233–3243. [Google Scholar] [CrossRef] [PubMed]
  79. Durán, A.; Martínez, G.C.; Pastor, M. Development and validation of AMANDA, a new algorithm for selecting highly relevant regions in Molecular Interaction Fields. J. Chem. Inf. Model. 2008, 48, 1813–1823. [Google Scholar] [CrossRef] [Green Version]
  80. Pándy-Szekeres, G.; Munk, C.; Tsonkov, T.M.; Mordalski, S.; Harpsøe, K.; Hauser, A.S.; Bojarski, A.J.; Gloriam, D.E. GPCRdb in 2018: Adding GPCR structure models and ligands. Nucleic Acids Res. 2018, 46, D440–D446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Ballesteros, J.A.; Weinstein, H. Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors. In Receptor Molecular Biology; Sealfon, S.C., Ed.; Academic Press: Cambridge, MA, USA, 1995; Volume 25, pp. 366–428. ISBN 1043-9471. [Google Scholar]
  82. Baxter, C.A.; Murray, C.W.; Clark, D.E.; Westhead, D.R.; Eldridge, M.D. Flexible docking using Tabu search and an empirical estimate of binding affinity. Proteins 1998, 33, 367–382. [Google Scholar] [CrossRef]
  83. Eldridge, M.D.; Murray, C.W.; Auton, T.R.; Paolini, G.V.; Mee, R.P. Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes. J. Comput. Aided Mol. Des. 1997, 11, 425–445. [Google Scholar] [CrossRef] [PubMed]
  84. Verdonk, M.L.; Cole, J.C.; Hartshorn, M.J.; Murray, C.W.; Taylor, R.D. Improved protein-ligand docking using GOLD. Proteins 2003, 52, 609–623. [Google Scholar] [CrossRef] [PubMed]
  85. Eberhardt, J.; Santos-Martins, D.; Tillack, A.F.; Forli, S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. J. Chem. Inf. Model. 2021, 61, 3891–3898. [Google Scholar] [CrossRef] [PubMed]
  86. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  87. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  88. Huang, J.; MacKerell, A.D.J. CHARMM36 all-atom additive protein force field: Validation based on comparison to NMR data. J. Comput. Chem. 2013, 34, 2135–2145. [Google Scholar] [CrossRef] [Green Version]
  89. Chen, H.; Maia, J.D.C.; Radak, B.K.; Hardy, D.J.; Cai, W.; Chipot, C.; Tajkhorshid, E. Boosting Free-Energy Perturbation Calculations with GPU-Accelerated NAMD. J. Chem. Inf. Model. 2020, 60, 5301–5307. [Google Scholar] [CrossRef]
  90. Schrödinger LLC. QikProp. 2021. Available online: https://www.schrodinger.com/products/qikprop (accessed on 31 March 2022).
  91. Jorgensen, W.L.; Duffy, E.M. Prediction of drug solubility from structure. Adv. Drug Deliv. Rev. 2002, 54, 355–366. [Google Scholar] [CrossRef]
  92. ChemAxon Ltd. MarvinSketch. 2022. Available online: https://chemaxon.com/products/marvin (accessed on 31 March 2022).
Molecules 27 02626 g001 550
Figure 1. Neurochemical and biochemical mechanisms mediating 5-HT6 receptor functions. In addition to the activation of cAMP signalling pathways, 5-HT6 receptors activate the extracellular signal regulated kinase1/2 (ERK1/2) via Fyn-dependent pathway. The proposed neurochemical circuitry for 5-HT6 receptors to influence cognition involves modulation of cholinergic and/or glutamatergic activity through GABAergic interneurons. AC: adenylate cyclase; ACH: acetylcholine; Glut: glutamate. (+): stimulation; (—): inhibition. (adapted and modified from [11]).
Figure 1. Neurochemical and biochemical mechanisms mediating 5-HT6 receptor functions. In addition to the activation of cAMP signalling pathways, 5-HT6 receptors activate the extracellular signal regulated kinase1/2 (ERK1/2) via Fyn-dependent pathway. The proposed neurochemical circuitry for 5-HT6 receptors to influence cognition involves modulation of cholinergic and/or glutamatergic activity through GABAergic interneurons. AC: adenylate cyclase; ACH: acetylcholine; Glut: glutamate. (+): stimulation; (—): inhibition. (adapted and modified from [11]).
Molecules 27 02626 g001
Molecules 27 02626 g002 550
Figure 2. Structural formula of 5HT6 receptor ligand from literature.
Figure 2. Structural formula of 5HT6 receptor ligand from literature.
Molecules 27 02626 g002
Molecules 27 02626 g003 550
Figure 3. Structural formula of selected 5HT6 receptor candidates from ZINC database.
Figure 3. Structural formula of selected 5HT6 receptor candidates from ZINC database.
Molecules 27 02626 g003
Molecules 27 02626 g004 550
Figure 4. Homology model of 5HT6 receptor with docked best candidate compound ZINC00756618 and interacting AA residues. Green lines: hydrogen bonds; purple: aromatic interactions; grey: hydrophobic interactions; orange: anion—PI interaction.
Figure 4. Homology model of 5HT6 receptor with docked best candidate compound ZINC00756618 and interacting AA residues. Green lines: hydrogen bonds; purple: aromatic interactions; grey: hydrophobic interactions; orange: anion—PI interaction.
Molecules 27 02626 g004
Molecules 27 02626 g005 550
Figure 5. Homology model of 5HT6 receptor with docked second best candidate compound ZINC20762773 and interacting AA residues. Green lines: hydrogen bonds; purple: aromatic interactions; grey: hydrophobic interactions; orange: salt bridge or anion—PI interaction.
Figure 5. Homology model of 5HT6 receptor with docked second best candidate compound ZINC20762773 and interacting AA residues. Green lines: hydrogen bonds; purple: aromatic interactions; grey: hydrophobic interactions; orange: salt bridge or anion—PI interaction.
Molecules 27 02626 g005
Molecules 27 02626 g006 550
Figure 6. RMSD of 5HT6 backbone atoms during the production phase of complex with ZINC00756618.
Figure 6. RMSD of 5HT6 backbone atoms during the production phase of complex with ZINC00756618.
Molecules 27 02626 g006
Molecules 27 02626 g007 550
Figure 7. FEP results for bound (red line) and unbound (black line) conformation of the 5HT6–ZINC00756618 complex.
Figure 7. FEP results for bound (red line) and unbound (black line) conformation of the 5HT6–ZINC00756618 complex.
Molecules 27 02626 g007
Table
Table 1. List of literature compounds for 5HT6 receptor, with their docking scores.
Table 1. List of literature compounds for 5HT6 receptor, with their docking scores.
NoNameGold ScoreAutodock Vina Binding Energy (kcal/mol)AQVNEIIP
1Intepirdine42.36−5.942.9090.009
2Cerlapirdine (SAM-531)36.63−8.132.880.000
3Idalopirdine 34.59−9.1932.550.089
4SGS51833.92−9.3662.780.035
5Masupirdine (SUVN-502)35.55−10.522.790.032
6AVN-21132.30−5.3353.080.073
Table
Table 2. List of the best candidate compounds for 5HT6 receptor, based on their docking scores.
Table 2. List of the best candidate compounds for 5HT6 receptor, based on their docking scores.
NoNameGold ScoreAutodock Vina Binding Energy (kcal/mol)AQVNEIIP
1ZINC0075661843.17−13.132.550.090
2ZINC2076277339.05−11.572.550.090
3ZINC2050190535.01−11.472.550.090
4ZINC0876497134.65−9.112.720.055
5ZINC7232453534.48−10.72.720.055
6ZINC0052622334.38−9.1352.990.009
7ZINC0876499333.48−9.3712.720.055
8ZINC1289277032.94−8.9952.720.055
9ZINC1827704932.65−9.0722.720.055
10ZINC1932360632.08−10.042.720.055
Table
Table 3. Calculated absorption, distribution, metabolism, elimination, and toxicity (ADMET) parameters of the compounds.
Table 3. Calculated absorption, distribution, metabolism, elimination, and toxicity (ADMET) parameters of the compounds.
CompoundMWRBDMMVDHBAHBPSAlogPlogSPCacoPM%HOAVRFVRTQPlogBBCNS MPO
ZINC00756618390.563.21304.914.4531.55.24−5.6814168100110.333.22
ZINC20762773390.568.11287.32.253.7586.83.98−4.9134849500−0.924.81
ZINC20501905390.556.71317.706.546.93.48−3.34125485000.293.22
ZINC08764971352.485.91133.52465.43.74−3.921475510000−0.554.73
ZINC72324535357.536.31182.21345.34.82−5.971440410001−0.284.53
ZINC00526223342.432.61102.32469.63.99−5.73839310001−0.754.41
ZINC08764993386.888.71202.52463.34.46−5.391724510000−0.393.94
ZINC12892770352.485.11179.72468.24.03−4.731159510000−0.694.58
ZINC18277049357.406.81087.122.562.24.61−5.691541810001−0.164.08
ZINC19323606357.404.71082.31643.72.55−2.48175282000.844.14
MW: Molecular weight; RB: Number of rotatable bonds; DM: computed dipole moment; MV: total solvent-accessible volume; DHB: estimated number of hydrogen-bond donors; AHB: estimated number of hydrogen-bond acceptors; PSA: van der Waals surface area of polar nitrogen and oxygen atoms and carbonyl carbon atoms; logP: predicted octanol/water partition coefficient; log S: predicted aqueous solubility; PCaco: predicted apparent Caco-2 cell permeability; PM: number of likely metabolic reactions; % HOA: predicted human oral absorption percentage; VRF: number of violations of Lipinski rule of five (the rules are as follows: MW< 500, log P < 5, DHB ≤ 5, AHB ≤ 10, positive PSA value); VRT: number of violations of Jorgensen rule of three (the rules are as follows: log S > −5.7, PCaco > 22 nm/s, PM < 7));QP log BB, brain/blood (good values to range 95% of drugs are from −3 to 1.2); CNS MPO. The result is a decimal number between 0 and 6, with “ideal” molecules scoring 6.
Table
Table 4. List of literature compounds for 5HT6 receptor, with their calculated QPlogBB and CNS MPO values.
Table 4. List of literature compounds for 5HT6 receptor, with their calculated QPlogBB and CNS MPO values.
NoNameQPlogBBCNS MPO
1Intepirdine−0.255.36
2Cerlapirdine (SAM-531)−0.744.89
3Idalopirdine 0.533.07
4SGS5180.133.41
5Masupirdine (SUVN-502)0.464.52
6AVN-211−0.295.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bojić, T.; Sencanski, M.; Perovic, V.; Milicevic, J.; Glisic, S. In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease. Molecules 2022, 27, 2626. https://doi.org/10.3390/molecules27092626

AMA Style

Bojić T, Sencanski M, Perovic V, Milicevic J, Glisic S. In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease. Molecules. 2022; 27(9):2626. https://doi.org/10.3390/molecules27092626

Chicago/Turabian Style

Bojić, Tijana, Milan Sencanski, Vladimir Perovic, Jelena Milicevic, and Sanja Glisic. 2022. "In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease" Molecules 27, no. 9: 2626. https://doi.org/10.3390/molecules27092626

Article Metrics

Back to TopTop