- 1Institute of Physiologically Active Compounds at Federal Research Center of Problems of Chemical Physics and Medicinal Chemistry, Russian Academy of Sciences, Chernogolovka, Russia
- 2Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia
- 3Institute of Organic Synthesis, Russian Academy of Sciences, Yekaterinburg, Russia
- 4Department of Organic and Biomolecular Chemistry, Ural Federal University, Yekaterinburg, Russia
- 5Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia
- 6Department of Pharmacology of the Institute of Biodesign and Complex System Modeling of Biomedical Science & Technology Park of Sechenov I.M., First Moscow State Medical University, Moscow, Russia
- 7Department of Environmental Health Sciences, University of Michigan, Ann Arbor, MI, United States
- 8Department of Neurology, University of Michigan, Ann Arbor, MI, United States
- 9Center of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, United States
- 10Michigan Institute for Computational Discovery and Engineering, University of Michigan, Ann Arbor, MI, United States
We investigated the inhibitory activities of novel 9-phosphoryl-9,10-dihydroacridines and 9-phosphorylacridines against acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and carboxylesterase (CES). We also studied the abilities of the new compounds to interfere with the self-aggregation of β-amyloid (Aβ42) in the thioflavin test as well as their antioxidant activities in the ABTS and FRAP assays. We used molecular docking, molecular dynamics simulations, and quantum-chemical calculations to explain experimental results. All new compounds weakly inhibited AChE and off-target CES. Dihydroacridines with aryl substituents in the phosphoryl moiety inhibited BChE; the most active were the dibenzyloxy derivative 1d and its diphenethyl bioisostere 1e (IC50 = 2.90 ± 0.23 µM and 3.22 ± 0.25 µM, respectively). Only one acridine, 2d, an analog of dihydroacridine, 1d, was an effective BChE inhibitor (IC50 = 6.90 ± 0.55 μM), consistent with docking results. Dihydroacridines inhibited Aβ42 self-aggregation; 1d and 1e were the most active (58.9% ± 4.7% and 46.9% ± 4.2%, respectively). All dihydroacridines 1 demonstrated high ABTS•+-scavenging and iron-reducing activities comparable to Trolox, but acridines 2 were almost inactive. Observed features were well explained by quantum-chemical calculations. ADMET parameters calculated for all compounds predicted favorable intestinal absorption, good blood–brain barrier permeability, and low cardiac toxicity. Overall, the best results were obtained for two dihydroacridine derivatives 1d and 1e with dibenzyloxy and diphenethyl substituents in the phosphoryl moiety. These compounds displayed high inhibition of BChE activity and Aβ42 self-aggregation, high antioxidant activity, and favorable predicted ADMET profiles. Therefore, we consider 1d and 1e as lead compounds for further in-depth studies as potential anti-AD preparations.
1 Introduction
This paper represents a continuation of our ongoing program of drug discovery for neurodegenerative conditions, with a particular emphasis on Alzheimer’s disease (AD). We have focused on AD for several reasons. For example, among the various types of age-related dementias, AD is the most prevalent. Moreover, AD is devastating for patients and their caregivers. It begins with memory loss, progresses to steady declines in cognitive function, and ultimately results in total disability and death (Calabro et al., 2021). Furthermore, there is no cure; treatment has been limited to partial alleviation of symptoms using the few drugs thus far approved by the FDA (Huang and Mucke, 2012; Cummings et al., 2019).
A major impediment to progress in drug discovery for AD is that the disease does not have a single cause; rather, it is multifactorial (Carreiras et al., 2013; Bachurin et al., 2017). Accordingly, our approach is to combine into single molecules the capabilities to inhibit several of the biological macromolecules or processes involved in the pathogenesis of the disease (Makhaeva et al., 2019c; Bachurin et al., 2023). In the present investigation, we have included the following targets: cholinergic neurotransmission, oxidative stress, and Aβ aggregation, as discussed in turn below.
One of the hallmarks of AD pathology is the demise of cholinergic neurons (Ballard et al., 2005; Hampel et al., 2019). Consequently, during the clinical course of the disease, the activity of acetylcholinesterase (AChE) decreases. In an apparent attempt to compensate for the loss of AChE activity, the activity of butyrylcholinesterase (BChE) increases (Ballard et al., 2005). Nevertheless, in the AD brain, there is an overall deficit of the neurotransmitter, acetylcholine (ACh). Therefore, inhibitors of one or both cholinesterases have been used to increase ACh concentrations in an effort to offset cognitive decline due to impaired cholinergic neurotransmission (Lane et al., 2006; Macdonald et al., 2014). Unfortunately, inhibition of cholinesterase activity on its own is insufficient to halt or retard the underlying neuropathic processes (Martinez and Castro, 2006; Agatonovic-Kustrin et al., 2018).
Among the mechanisms thought to contribute to AD pathogenesis is oxidative stress, which can lead to neuronal death by damaging the lipid and protein constituents of mitochondria, other organelles, and plasma membranes. This situation arises when the concentrations of reactive oxygen species (ROS) in susceptible sites reach critical levels due to their rates of production exceeding their rates of removal (Ortiz et al., 2017; Pohanka, 2018). Although antioxidant processes in the brain decrease with increasing age, these processes are further compromised in AD brain. Consequently, antioxidants are considered to be of adjunctive value in treating AD (Moreira et al., 2005). Indeed, antioxidant properties can be incorporated into drug candidates that contain pharmacophores directed against other targets, such as cholinesterases (Rosini et al., 2011; Ismaili and Romero, 2017; Makhaeva et al., 2020).
Accumulation of Aβ plaques is another process that has been amply linked to AD pathology (Walsh et al., 2005; Castellani et al., 2019; Jeremic et al., 2021). In view of this, under its new accelerated approval program, the FDA has recently approved two monoclonal antibody drugs that could potentially attenuate AD progression by clearing Aβ plaques from the brain (Vaz et al., 2022; van Dyck et al., 2023). However, these approvals have been controversial owing to concerns regarding safety and efficacy (Reardon, 2023). Another approach is to understand the mechanisms of formation of neurotoxic Aβ oligomers and to design compounds to intervene (Jeremic et al., 2021), which was the strategy we employed in our present investigation.
For the current work, we chose acridines as our core structures for building potential multifunctional agents. Our rationale included the fact that in drug discovery for neurodegenerative diseases, acridines are considered to be privileged scaffolds (Bongarzone and Bolognesi, 2011). Moreover, derivatives of acridine have proved to exert therapeutic effectiveness against such diverse illnesses as bacterial infections (Wainwright, 2001), malaria (Girault et al., 2000), trypanosomiasis and leishmaniasis (Gamage et al., 1997), viral infections (Lyakhov et al., 2000), cancer (Denny, 2002), inflammation and diabetes (Sondhi et al., 2010; Mallu et al., 2017) and prion diseases (Korth et al., 2001; Klingenstein et al., 2006; Collinge et al., 2009). They also serve as ideal points of departure for the generation of novel multitarget lead compounds and drug candidates for neurodegenerative and other conditions (Csuk et al., 2009; Antosova et al., 2011; Bongarzone and Bolognesi, 2011). Furthermore, it is well established that numerous compounds containing an acridine moiety display potent inhibitory activity toward AChE and BChE (Kozurkova et al., 2011; Jin, 2014). Tacrine, (9-amino-1,2,3,4-tetrahydroacridine), an especially well-known acridine-based cholinesterase inhibitor, was the first FDA-approved drug for the treatment of the cognitive symptoms of AD (Crismon, 1994). Along with tacrine, the acridine template has been employed for the design and development of new multifunctional molecules as potential therapeutic agents for AD (Makhaeva et al., 2017; Sharma and Piplani, 2017; Chufarova et al., 2018; Lotfi et al., 2020; Tseng et al., 2020; Maciejewska et al., 2022).
In particular, in our previous work, we studied the esterase profile and estimated the antioxidant activity (AOA) of four groups of new 9-substituted acridine derivatives (Makhaeva et al., 2017). We demonstrated that 9-amino-N-methyl-9,10-dihydroacridine derivatives combined effective inhibition of AChE and BChE with a high level of radical-scavenging activity.
The aim of the present investigation was to study the possibility of creating multifunctional agents for AD therapy based on new 9-phosphoryl-9,10-dihydroacridines and 9-phosphorylacridines with various structures of the phosphoryl fragment, obtained using our original strategy of functionalization of acridines by nucleophilic aromatic substitution of hydrogen (SNH) (Chupakhin et al., 1994; Charushin and Chupakhin, 2004; Chupakhin and Charushin, 2017). We were encouraged to pursue 9-phosphoryl derivatives by the observation that connecting acridine with a phosphorus moiety could enhance penetration of acridines through cell membranes (Demkowicz et al., 2016). In addition, recently studied N- and O-phosphorylated tacrines (Przybylowska et al., 2022) demonstrated high anti-AChE and anti-BChE activities, and some of them had reduced hepatotoxicity compared to tacrine.
In summary, we chose to investigate the effects of our new compounds on the components of AD pathogenesis outlined above: cholinergic neurotransmission, oxidative stress, and Aβ aggregation. Accordingly, we determined the following characteristics of the synthesized compounds: the esterase profile, i.e., the inhibitory activities against AChE, BChE, and a structurally related off-target carboxylesterase (CES, EC 3.1.1.1); inhibition of Aβ self-aggregation; and antioxidant activities in the ABTS and FRAP assays. We then analyzed the observed effects using quantum chemical calculations and computational molecular modeling. Finally, we carried out computational predictions of ADMET properties.
2 Materials and methods
2.1 Synthesis of compounds
9-Phosphoryl-9,10-dihydroacridines 1a–f and 9-phosphorylacridines 2a–f were obtained according to Figure 1 as described in our recent publication (Shchepochkin et al., 2021).
2.2 Biological investigations
Experiments were performed in alignment with the standard operating procedures approved by IPAС RAS as described below.
2.2.1 In vitro inhibition of AChE, BChE, and CES activities
Enzymes, substrates, and reference compounds were from Sigma-Aldrich (St. Louis, MO, United States). The activity of enzymes was measured spectrophotometrically, as described in detail in (Makhaeva et al., 2017) using ATCh iodide, BTCh iodide, and 4-NPA as substrates for AChE, BChE, and CES, respectively. Experimental conditions: K,Na-phosphate buffer (100 mM), 25°С, pH 7.5 for AChE and BChE and pH 8.0 for CES assay. Measurements were carried out on a BioRad Benchmark Plus microplate spectrophotometer (Hercules, CA, United States).
Test compounds were dissolved in DMSO; final concentration of solvent in the incubation mixture was 2% (v/v). Initial assessment of inhibitory activity was carried out by determining the degree of enzyme inhibition at a compound concentration of 20 µM. For the most active compounds (inhibition ≥35%), IC50 values were determined.
Mechanism of BChE inhibition was assessed by a detailed analysis of enzyme kinetics with three increasing concentrations of inhibitor and six substrate concentrations as described in detail in Makhaeva et al. (2017).
2.2.2 Inhibition of β-amyloid (1-42) (Aβ42) self-aggregation
Inhibition of Aβ42 self-aggregation by test compounds was studied using the thioflavin T (ThT) fluorescence method (LeVine, 1999; Bartolini et al., 2003; Munoz-Ruiz et al., 2005) with minor modifications as described in detail in (Makhaeva et al., 2023). Lyophilized HFIP-pretreated Aβ42 (AnaSpec Inc., 0.5 mg) was dissolved in DMSO to obtain a stable 500 μM solution. The samples of 50 μM Aβ42 in 215 mM Na-phosphate buffer pH 8.0 were incubated for 24 h at 37°C in the absence or presence of 100 µM test compounds. Myricetin and Propidium iodide in the same concentration were used as references. After incubation, 5 μM ThT in 50 mM glycine-NaOH buffer pH 8.5 was added, and the fluorescence was measured at 440 nm (exc.) and 485 nm (emis.) with a FLUOStar Optima microplate reader (LabTech, Ortenberg, Germany). Blanks consisted of 215 mM Na-phosphate buffer, pH 8.0, 20% (v/v) DMSO or test compounds, respectively.
The inhibition (%) of Aβ42 self-aggregation by the test compounds was calculated using the following equation:
where IFi and IFo are the fluorescence intensities obtained for Aβ42 in the presence or absence of inhibitor, respectively, after subtracting the fluorescence of respective blanks.
2.2.3 Antioxidant activity (AOA)
2.2.3.1 ABTS radical cation scavenging activity assay
Radical scavenging activity of the compounds was evaluated using the ABTS radical cation (2,2ʹ-azinobis-(3-ethylbenzothiazoline-6-sulfonic acid), ABTS•+) decolorization assay (Re et al., 1999), with minor modifications described in detail in Makhaeva et al. (2019a). To prepare the solution of ABTS•+, equal volumes of aqueous solutions of 7 mM ABTS and 2.45 mM potassium persulfate were incubated for 12–16 h at room temperature in the dark. Radical scavenging capacity of the compounds was determined by mixing 10 μL of compound solution in DMSO with 240 μL of ABTS•+ working solution in ethanol (100 μM final concentration). The reduction in ABTS•+ absorbance was measured spectrophotometrically at 734 nm using a xMark UV/VIS microplate spectrophotometer (Bio-Rad, Hercules, CA, United States) for 1 h with an interval of 10 min compared to a standard synthetic antioxidant Trolox (6-hydroxy-2,5,7,8-tetramethylchroman-2-carboxylic acid). Ascorbic acid was used as a positive control.
The antioxidant activity of the compounds was reported as Trolox equivalent antioxidant capacity (TEAC values) as the ratio of the slopes of the concentration−response curves, test compound/Trolox. The IC50 values for the test compounds were also determined.
2.2.3.2 FRAP assay
The ferric reducing antioxidant power (FRAP) assay proposed by Benzie and Strain (1996); Benzie and Strain (1999) modified to be performed in 96-well microplates as described in detail in Makhaeva et al. (2020) was used. 10 μL (0.5 mM) of the tested or reference compounds were mixed with 240 µL of FRAP reagent (2.5 mL of 10 mM TPTZ (2,4,6-tris(2-pyridyl)-s-triazine) solution in 40 mM HCl, 2.5 mL of 20 mM FeCl3 in distilled water and 25 mL of 0.3 M acetate buffer pH 3.6); the absorbance of the mixture was measured spectrophotometrically (λ = 593 nm) with a FLUOStar OPTIMA microplate reader at 600 nm after a 1 h incubation at 37°C against a blank. Trolox was used as a reference compound. The results were expressed as Trolox equivalents (TE)—the ratio of the concentrations of Trolox and the test compound resulting in the same effect on ferric reducing activity.
2.3 Statistical analyses
Plots, linear regressions, and IC50 values were determined using Origin 6.1 for Windows, OriginLab (Northampton, MA, United States). All tests were performed at least in triplicate in three independent experiments. Results are presented as mean ± SEM calculated using GraphPad Prism version 6.05 for Windows (San Diego CA, United States).
2.4 Molecular modeling studies
2.4.1 Molecular docking
The Calculator Plugins of MarvinSketch 21.14.0, ChemAxon (https://www.chemaxon.com, accessed on 27 January 2023) and MolGpKa (Pan et al., 2021) (https://xundrug.cn/molgpka, accessed on 27 January 2023) were used to estimate the pKa values of the ligands. Because the pKa values of the acridine fragments for all considered compounds were estimated to be below 5, all conjugates were utilized for molecular docking in their neutral form. The ligand compounds were optimized using a DFT quantum chemistry method (B3LYP/6-31G*, GAMESS-US (Schmidt et al., 1993) software). For molecular docking, the optimized structures of the ligands were employed, with partial atomic charges obtained from QM results based on the Löwdin scheme (Löwdin, 1970).
The protein targets used for docking included X-ray structures of human AChE co-crystallized with donepezil (PDB: 4EY7, chain A) (Cheung et al., 2012), an optimized X-ray structure of human BChE (PDB: 1P0I) (Nicolet et al., 2003; Masson et al., 2013), and all conformers of the Aβ42 solution NMR structure (PDB ID 1IYT) (Crescenzi et al., 2002).
AutoDock 4.2.6 software (Morris et al., 2009) was utilized to perform molecular docking. The docking grid box was set to cover the entire active site gorge of AChE (22.5 Å × 22.5 Å × 22.5 Å) and BChE (15 Å × 20.25 Å × 18 Å), as well as the entire Aβ42 molecule for all conformers (43.5 Å × 28.5 Å × 54.75 Å), with a grid spacing of 0.375 Å used in all cases. The Lamarckian Genetic Algorithm (LGA) (Morris et al., 1998) was used with 256 runs, 25 × 106 evaluations, 27 × 104 generations, and a population size of 3,000. Figures were created using PyMol (www.pymol.org, accessed on 21 July 2016).
2.4.2 Molecular dynamics simulations
For molecular dynamics (MD) simulations, an initial complex of BChE with the most active compound 1d was taken from molecular docking results. Crystallographic waters found in the initial protein X-ray structure were added, except for ones overlapping with the inhibitor molecule. TIP3P water molecules were added, forming a box with boundaries extending 10 Å beyond the protein using the VMD solvate module. Na+ and Cl− ions were added up to 150 mM ion concentration, and the neutral state of the whole system was achieved using the VMD autoionize module. Force field parameters of the ligand were derived with the CharmGenFF program v. 2.5 (Vanommeslaeghe et al., 2010).
MD simulation was performed with NAMD 3.0b2 software (Phillips et al., 2020) using the CHARMM36 force field (Vanommeslaeghe et al., 2010). During MD simulation, the system was maintained at constant temperature (298 K) and pressure (1 atm, NPT ensemble), using Langevin dynamics and the Nosé-Hoover barostat; the timestep was 1 fs. Periodic boundary conditions and PME electrostatics were applied. Before the productive MD run, a system minimization was conducted during 3,000 steps; next, a 1 ns solvent equilibration run was carried out with the protein and ligand coordinates fixed, and then the structure was fully minimized during 3,000 steps. After that, a 50 ns production run was performed. Analysis of the obtained MD trajectory was done with the VMD 1.9.4a53 package (Humphrey et al., 1996). Molecular modeling figures were prepared with PyMOL.
2.4.3 Quantum-chemical analysis of AOA
Quantum-chemical calculations were performed with the Gaussian 16 (Frisch et al., 2016) and Priroda 19 (Laikov and Ustynyuk, 2005; Laikov, 2020) packages. At the initial stage, the conformational analysis was done in the PBE0 functional (Adamo and Barone, 1999) and TZVP basis set (Schäfer et al., 1994) in the gas phase by the Priroda software. The obtained conformations were optimized using the long-range-corrected functional CAM-B3LYP (Yanai et al., 2004) in the 6-31++G (d,p) basis set (Rassolov et al., 2001) by the Gaussian program. The energies were refined by single-point calculations in the aug-cc-pVTZ basis set (Papajak et al., 2011). Solvent effects were taken into account by applying the Solvation Model based on density (SMD) (Marenich et al., 2009). Six compounds were investigated; namely, three dihydroacridines (1a, 1c, 1e) and three acridines (2a, 2c, 2e).
2.4.4 Prediction of ADMET, physicochemical, and PAINS profiles
Lipophilicity (LogPow) and aqueous solubility (pS) were estimated by the ALogPS 3.0 neural network model implemented in the OCHEM platform (Sushko et al., 2011). Human intestinal absorption (HIA) (Radchenko et al., 2016), blood–brain barrier distribution/permeability (LogBB) (Dyabina et al., 2016; Radchenko et al., 2020), and hERG-mediated cardiac toxicity risk (channel affinity pKi and inhibitory activity pIC50) (Radchenko et al., 2017) were estimated using the integrated online service for the prediction of ADMET properties. This service implements predictive QSAR models based on accurate and representative training sets, fragmental descriptors, and artificial neural networks. The quantitative estimate of drug-likeness (QED) values (Bickerton et al., 2012) were calculated, and the Pan Assay INterference compoundS (PAINS) alerts were checked using RDKit version 2021.09.2 software.
3 Results and discussion
3.1 Chemistry. Synthesis of 9-substituted acridines
Currently, 9-substituted acridines can be obtained by several approaches. The classical method for the preparation of acridine derivatives is based on high-temperature condensation of diphenylamine and carboxylic acid in the presence of zinc chloride (Bernthsen, 1884) (Pathway A, Figure 2). Another method consists of preliminary functionalization of acridine and subsequent substitution of the introduced group by a nucleophile in the presence of transition-metal salts (Liu et al., 2013; Zhilyaev et al., 2022) (Pathway B, Figure 2). The presented strategies do not exclude a number of shortcomings. For example, the desire to introduce auxiliary groups, the use of catalysts and additional ligands, the employment of high temperature processes, and the need to separate reaction products from metal-containing byproducts complicate the preparation of target acridines. In addition, such accompanying negative factors as the high cost of the reaction process and, in some cases, low yields of the final products make the method inefficient and unprofitable. However, an alternative approach to functionalized acridines is nucleophilic aromatic substitution of hydrogen (SNH) (Chupakhin et al., 1994; Charushin and Chupakhin, 2005; Makosza, 2010; Mąkosza, 2011; Mąkosza and Wojciechowski, 2014). This method makes it possible to obtain 9-substituted acridines efficiently and directly, avoiding the disadvantages mentioned above.
The SNH reactions proceed in two stages. The first one is nucleophilic addition to an aromatic substrate to form an intermediate (ϭH-adduct). The ϭH-adducts, e.g., in the case of acridine, are stable compounds, which makes it possible to comprehensively study their physicochemical properties and biological activities (Chupakhin et al., 2017; Volkov et al., 2018). The second stage is characterized by an oxidative process with formation of the target-functionalized acridine (Pathway C, Figure 2).
We used the SNH approach described in detail in a recent publication (Shchepochkin et al., 2021) to obtain the 9-phosphoryl-9,10-dihydroacridines 1a–f and 9-phosphorylacridines 2a–f studied in this work (Figure 1). Briefly, acridine and the corresponding phosphonate were stirred without solvent under an argon atmosphere at 75°C–80°C for 2–7 h. 9-Phosphoryl-9,10-dihydroacridines 1a-f were isolated in 85%–95% yields. The corresponding aromatized compounds 2a-f were obtained in 81%–89% yields by electrochemical oxidation of 1a–f.
Figure 3 illustrates the structures of the previously described compounds, which were the subject of the research in the present work.
FIGURE 3. Structures of the studied acridine derivatives: 9-phosphoryl-9,10-dihydroacridines 1a–f and 9-phosphorylacridines 2a–f.
3.2 Inhibition studies of AChE, BChE and CES. Structure-activity relationships
Table 1 lists the components of the esterase profile of our compounds, i.e., their inhibitory potencies against the enzymes of cholinergic neurotransmission, AChE and BChE, and a structurally related off-target enzyme, CES (Makhaeva et al., 2016; Makhaeva et al., 2019b).
As can be seen from Table 1, all tested compounds weakly inhibited AChE, while a number of compounds demonstrated effective inhibition of BChE.
The most effective BChE inhibitors were found among the series of phosphorylated dihydroacridines. Moreover, the influence of the structure of the phosphorus-containing fragment on the inhibitory activity is clearly seen: dialkylphosphonates weakly inhibited BChE, but the introduction of aryl fragments increased the activity. Thus, diphenylphosphonate 1c was a moderate inhibitor (IC50 = 48.0 ± 3.1 µM), while dibenzyloxy derivative 1d and its diphenethyl bioisostere 1e were effective inhibitors of BChE with IC50 = 2.90 ± 0.23 µM and 3.22 ± 0.25 µM, respectively. The introduction of a chlorine atom into the para-position of the phenethyl fragment of 1e (compound 1f) reduced tenfold the anti-BChE activity (IC50 = 21.7 ± 1.7 µM).
At the same time, aromatized analogs 2 of compounds 1 were largely inactive as BChE inhibitors. An exception was the acridine analog of the most effective dihydroacridine 1d–dibenzyloxy derivative 2d, which inhibited BChE with IC50 = 6.90 ± 0.55 µM, i.e., it was only 2.4 times less effective than the corresponding dihydroacridine.
Weak inhibition of CES, an enzyme that catalyzes the hydrolysis of numerous ester-containing drugs, is a positive property of the studied compounds because this decreases the likelihood of unwanted drug-drug interactions in the event of their use in AD therapy (Tsurkan et al., 2013; Makhaeva et al., 2019b).
3.3 Kinetic studies of BChE inhibition
The most active BChE inhibitor 1d was selected for inhibition kinetics studies. Lineweaver-Burk plots-double reciprocal representations of Michaelis-Menten parameters-were used to assess the type of inhibition. Graphical analysis of the kinetic data on BChE inhibition by the tested compound (Figure 4) demonstrates changes in both Km and Vmax values; this result attests to a mixed type of inhibition. The values obtained for the competitive component (Ki) and the non-competitive component (αKi) of the constants of BChE inhibition by compound 1d were 2.22 ± 0.19 µM and 11.5 ± 0.9 µM, respectively.
FIGURE 4. Double-reciprocal (Lineweaver-Burk) plots of the steady state inhibition of BChE by compound 1d.
3.4 Molecular docking to AChE and BChE
To investigate the structural origins of the observed differences in inhibitory activity, we performed molecular docking to AChE and BChE for the most active dihydroacridine derivatives, namely, the bioisosteric compounds 1d and 1e and their aromatized analogs 2d and 2e.
According to the results of molecular docking, compounds with benzyloxy and phenethyl substituents 1d,e, and 2d,e were able to bind to the AChE active site gorge in a similar manner (Figure 5). However, in this case, the bulky dihydroacridine/acridine and benzyloxy/phenethyl groups were located below amino acid residues Tyr341 and Tyr124, which are the part of a narrow bottleneck (Sussman et al., 1991; Wlodek et al., 1997). This makes trafficking of the ligands to the bottom of the gorge kinetically hindered (Lushchekina and Masson, 2020), thus reducing their inhibitory efficacy.
FIGURE 5. Binding modes of compounds 1d,e (carbon atoms colored green/magenta, respectively) and 2d,e (carbon atoms colored with lighter shades of the same colors) in the active site gorge of human AChE (PDB: 4EY7), according to molecular docking results. Yellow dashes show direct hydrogen bonding interactions. Orange dashes show the hypothetical possibility of water molecule-mediated interactions.
With respect to BChE docking (Figure 6), a significant difference of binding to the enzyme active site at the bottom of the gorge was observed between dihydroacridine/acridine compounds 1d,e and 2d,e with benzyloxy/phenethyl substituents. Due to the wider active site gorge of BChE compared to AChE (Saxena et al., 1997), no kinetic difficulties of such binding were expected.
FIGURE 6. Binding mode of compounds 1d,e (carbon atoms colored green or magenta, respectively) (A), and 2d,e (carbon atoms colored with lighter shades of the same colors) (B) in the active site gorge of human BChE (PDB: 1P0I), according to molecular docking results. Yellow dashes show direct hydrogen bonding interactions. Orange dashes show the hypothetical possibility of water molecule-mediated interactions.
According to the docking results, the dihydroacridine derivatives 1d,e form direct hydrogen bonds with the Tyr332 and D70 side chains (Figure 6A). There is also a theoretical possibility of interaction of the phosphoryl group with the Thr120 side chain oxygen atom through a bridging water molecule (Figure 6A). According to X-ray crystallography (Koellner et al., 2000; Bourne et al., 2010; Berg et al., 2011), and molecular modeling (Nascimento et al., 2017; Makhaeva et al., 2021; van der Westhuizen et al., 2022), such interactions are typical for protein-ligand interactions in the active site gorge of cholinesterases.
Regarding acridines 2d and 2e, no specific interactions were observed, which agrees with their lower experimentally observed inhibitory activity. However, interactions of the acridine fragment of both compounds with the Trp82 main chain oxygen atom are possible (Figure 6B). We observed such a water-mediated interaction for the similar amiridine fragment and the homologous AChE residue Trp86 in molecular dynamics simulations (Makhaeva et al., 2021). As observed with dihydroacridine derivatives, the phosphoryl group of the acridine compounds was located in such a way to suggest that a similar water molecule-bridged interaction occurred with the oxygen atom of the Thr120 side chain (Figure 6B). Regarding the dibenzyloxy group of compound 2d, such interactions are also possible for the oxygen atom (Figure 6B), which is consistent with its higher inhibitory activity compared to compound 2e (Table 1).
3.5 Molecular dynamics simulations of the complex of compound 1d with BChE
We performed MD simulations for the complex of the most active compound 1d with BChE. During the MD trajectory, the ligand remained bound to the enzyme, though it had reoriented, and the resulting binding pose (Figure 7) differed from the docked one (Figure 6A) in specific details. While molecular docking showed that the dihydroacridine group interacted with side chains of the PAS residues Asp70 and Tyr332, these interactions were disrupted shortly after the start of the MD simulation (Figure 8A). Instead, the inhibitor reoriented so that the dihydroacridine group established new interactions with the side chain of Trp82 after 15 ns of simulation (Figure 8B). This residue is located deeper in the gorge than Asp70 and Tyr332 and forms the active site as a part of the cation-binding compartment. The newly established π-cation stacking or T-stacking interactions of the dihydroacridine fragment were maintained during the remainder of the MD trajectory. In addition, T-stacking interactions of one of the ligand’s phenyl rings with the Trp231 indole ring were observed during the entire MD trajectory (Figure 8C). The reorientation of the ligand in the active site gorge of BChE took place during the first 15 ns of the simulation and included temporary displacement of the phosphate group. After this, the phosphate group returned to the position outside the oxyanion hole and close to the Thr120 side chain (Figure 7; Figure 8D). The MD simulation confirmed that their interaction was maintained mainly through frequently exchanging bridging water molecules (Figure 9). However, occasional direct hydrogen bonding is also possible due to thermal motion (Figure 8D).
FIGURE 7. A snapshot from the MD simulation of the complex of compound 1d with BChE, showing the main interactions in the complex after the ligand reorientation discussed in the text. Yellow dashed lines indicate hydrogen bonds, and orange dashes show stacking interactions. The ligand’s atoms referred to in the following figure are labeled. Note the presence of a bridging water molecule whose hydrogen atoms form hydrogen bonds with the Tyr120 Oγ atom and the ligand O2 atom.
FIGURE 8. Interactions along the MD trajectory of the complex of compound 1d with BChE. The evolution of distances between specified atoms or center of masses for the indole rings of Trp residues during the MD simulation is shown. The ligand atoms are indicated according to the labeling in Figure 7. Bold lines show distances averaged over 1 ns of the simulation, and the light trace shadow shows distances recorded every 10 ps. (A) The increase of the shown distances reflects the disruption of interactions between the dihydroacridine group of the ligand and Asp70 and Tyr332 residues and following reorientation; (B) the decrease of the shown distances reflects the establishment of new interactions between the dihydroacridine group of the ligand and the side chain of Trp82 after the ligand reorientation; (C) T-stacking interaction of one of the ligand’s phenyl rings with the Trp231 indole ring was maintained through the whole trajectory and tended to be tighter by the end of the simulation; (D) after the ligand reorientation, the phosphoryl oxygen atom remained close enough to the Thr120 hydroxy group for indirect interactions via bridging water molecules. Only occasionally was the distance between them less than 3Å, which suggests the possibility of transient direct interactions in the form of weak hydrogen bonds.
FIGURE 9. Bridging of the ligand phosphate group and Thr120 side chain by water molecules (LigO2···H1/2-O-H1/2·· Thr120_Oγ) identified in the MD simulation of the complex between compound 1d with BChE. Each line shows the minimum of the two distances between each of the designated water molecule hydrogen atoms and enzyme or ligand oxygen atom. A water molecule was considered bridging when such H-O distances were less than 2 Å simultaneously for the ligand O2 and the enzyme Thr120_Oγ atoms; pairs of these distances for the same water molecule are shown with different shades of the same color. Along the MD trajectory, 133 different water molecules participated in bridging defined in such a way. The plot shows distances for a few water molecules that were most frequently identified as bridging. For clarity, distances for the other type of bridging (LigO2···H1/2-O·· Thr120_Hγ) are not shown.
3.6 Inhibition of β-amyloid (1-42) (Aβ42) self-aggregation
The in vitro inhibitory activity of the test compounds against Aβ42 self-aggregation was determined using the thioflavin T (ThT) fluorimetric assay (LeVine, 1999; Bartolini et al., 2003; Munoz-Ruiz et al., 2005; Bartolini et al., 2011). This widely used procedure is based on specific binding of the fluorescent dye ThT to the β-sheets of assembled amyloid fibrils leading to a significant increase in fluorescence signal (Biancalana and Koide, 2010). Therefore, the decrease in ThT fluorescence in the presence of the studied compound correlates with its activity to inhibit the formation of amyloid aggregates.
As can be seen from Table 1, the ability of the studied compounds to inhibit the self-aggregation of Aβ42 depends both on the structure of their acridine moiety and the substituents at the phosphorus atom; moreover, the structure-activity relationship largely coincides with that observed for BChE inhibition.
Indeed, acridine derivatives 2 inhibit the self-aggregation of Aβ42 very weakly or do not inhibit at all. The ability of dihydroacridines 1 to block Aβ42 self-aggregation depends on the structure of the phosphoryl fragment: dialkoxy derivatives are weak inhibitors, but the introduction of aryl substituents increases activity, which increases upon going from diphenoxy (1с) to dibenzyloxy (1d) substituents. Compound 1d exhibits the maximum inhibitory activity and suppresses the formation of amyloid aggregates by 58.9% ± 4.7%. Its bioisosteric analog, the diphenethyl derivative 1e, also has high activity, inhibiting Aβ42 self-aggregation by 46.9% ± 4.2%. Adding a chlorine atom to the para-position of the phenethyl fragment 1e (compound 1f) slightly reduces the inhibitor activity.
3.7 Molecular docking to Aβ42
Molecular docking did not reveal a direct explanation of the experimentally observed differences in the ability of acridine/dihydroacridine compounds to inhibit Aβ42 self-aggregation. There were slight differences in estimated binding energies (Supplementary Table S1) and number of contacts (Supplementary Table S2) of dihydroacridine derivatives 1d,e compared to acridine derivatives 2d,e with Aβ42 conformers obtained from the NMR structure (PDB: 1IYT) (Crescenzi et al., 2002) [see details in (Makhaeva et al., 2022)]. Although these differences could be considered minor, they are in line with different conformations of the Aβ42 peptide used as molecular docking targets.
Another plausible reason for the experimentally observed substantial differences between active dihydroacridine inhibitors and practically inactive acridine compounds could be related to the different ability of planar aromatic acridines to self-associate through π-π stacking interactions compared to non-planar dihydroacridines. Such self-stacking of acridine derivatives is well known (Evstigneev et al., 2006; Franco Pinto et al., 2022) and could compete with binding to the Aβ42 peptide (Gao et al., 1991; Buurma and Haq, 2008).
3.8 Antioxidant activity (AOA)
Table 2 lists the results of our determinations of the AOA of compounds 1a–f and 2a–f using the ABTS and FRAP assays. In the ABTS test, the radical cation ABTS•+ reacts with an antioxidant via one or both of the following pathways: single-electron transfer (SET) or hydrogen-atom transfer (HAT). In the FRAP test, the ferric 2,4,6-tripyridyl-s-triazine complex [Fe (TPTZ)2]3+ is reduced to the corresponding ferrous complex [Fe (TPTZ)2]2+ by the SET mechanism. In these tests, we used Trolox and ascorbic acid as the reference and positive control compounds, respectively.
All derivatives of 9,10-dihydroacridine 1 demonstrated high radical-scavenging activity in the ABTS test at levels near or slightly above that of Trolox, the standard antioxidant. In contrast, the 9-phosphoryl-acridines 2 were not active in the ABTS test. It should also be noted that compounds 1a–d exhibited a high initial rate of the ABTS•+ binding reaction, comparable to Trolox (the maximum degree of radical binding was achieved within 1 min). This fact is consistent with the SET mechanism for the antiradical activity of these compounds. A slight decrease in the radical-scavenging activity was observed for the phenethyl derivative 1e and its para-Cl-substituted analog 1f, which was also accompanied by a decrease in the reaction rate (the time to reach the maximum degree of radical binding increased from 1 to 30 min).
All dihydroacridine derivatives 1a–f also showed high iron-reducing activity in the FRAP test, at the level or above that of Trolox. Moreover, these compounds were relatively rapid antioxidants, reaching the maximum effect within 4 min. The high activity of dihydroacridines 1a–f in the Fe3+-reducing test is probably directly related to the ease of their electrochemical oxidation to acridines (Shchepochkin et al., 2021). The dimethoxy derivative 1a (1.60 ± 0.01 TE) was the lead compound in this series, and the least active compound was para-Cl-substituted diphenethylphosphoryl dihydroacridine 1f (0.81 ± 0.02 TE). In general, the results of the ABTS and FRAP tests were in good agreement with each other.
Phosphoryl acridines 2a–f had almost no activity in the FRAP test. Only dimethoxy and diethoxy derivatives 2a,b, which are the aromatic analogs of the most active compounds 1a,b in the dihydroacridine series, showed detectable activity.
3.9 Quantum-chemical calculations of AOA
The FRAP assay is a typical SET method (Gulcin, 2020; Spiegel et al., 2020; Florin Danet, 2021), where [Fe (TPTZ)2]3+ is reduced to [Fe (TPTZ)2]2+ by electron transfer from an antioxidant molecule. Unlike [Fe (TPTZ)2]3+, ABTS•+ can be reduced by acquiring both a hydrogen atom and an electron, depending on the antioxidant structure, solvent type, and pH (Ilyasov et al., 2020). In particular, SET is facilitated at acidic pH (Gulcin, 2020). In the ABTS method used in this work, the main solvent is ethanol at pH∼4.5. Under these conditions, the ABTS•+ molecule is deprotonated at both SO3 groups; it is reduced by acquiring an electron from the antioxidant, similarly to [Fe (TPTZ)2]3+ (Ivanova et al., 2020).
In both methods, reduction can occur by either the direct SET mechanism characterized by the ionization potential (IP), or by a hydrogen atom abstraction from the antioxidant molecule, when an electron is transferred to the radical and a proton to a solvent molecule. In the latter case, the actual antioxidant mechanism may be a complex multistage process, but the overall reaction can be characterized by the enthalpy of a hydrogen atom abstraction, i.e., the bond dissociation enthalpy, also known as the bond dissociation energy (BDE).
The quantum mechanical (QM) characteristics of possible antioxidant reactions for dihydroacridines and acridines were calculated taking two items into account: 1) the protonation state of the investigated structures in water and ethanol for the FRAP and ABTS tests, respectively; and 2) the method of synthesizing acridines from dihydroacridines by two-electron transfer electrochemical oxidation (Shchepochkin et al., 2021). The second item suggested a two-stage mechanism for the interaction of the antioxidant with ABTS•+ or FRAP molecules.
The calculated IP and BDE values for acridines were much higher than the corresponding values for dihydroacridines. Moreover, the characteristics of acridines were close to each other, as were the characteristics of dihydroacridines, which agrees with similar experimental results for compounds 2a and 1a, respectively. A detailed description of the approach that was used and all calculated characteristics are given in Supplementary Materials. Here, we briefly summarize the results.
The obtained QM characteristics for the two-stage mechanism of interaction of the antioxidant and ABTS•+ or FRAP molecules allowed us to suggest that each dihydroacridine molecule donates two electrons, yielding the corresponding acridine molecule, and thereby reducing two molecules of ABTS•+ or FRAP. At the same time, the reduction of the second radical molecule was characterized by low values of the BDE, which was consistent with the high AOA of dihydroacridines in both tests. The poor AOA of acridines agrees with high values of their QM characteristics.
The suggested two-stage mechanism of the dihydroacridine antioxidant action is illustrated in Figure 10 for compound 1a in the ABTS test as an example. In path A, the antioxidant (dihydroacridine, AHCHN) sequentially loses an electron and a hydrogen atom to generate the protonated form of the corresponding acridine (AHN+). The latter dissociates to acridine (A); this process is characterized by the proton affinity (PA). In path B, the antioxidant sequentially loses two H-atoms to form the corresponding acridine.
FIGURE 10. Two-stage mechanism of the antioxidant action of dihydroacridines in the ABTS test using compound 1a as an example. The parent dihydroacridine 1a (AHCHN) is circled in green, the final acridine 2a (A) is circled in red, and the protonated form of acridine (AHN+) is circled in blue. The calculated QM characteristics of each stage are shown along the corresponding arrows indicating the reaction. “Solv” is the solvent molecule.
In both cases, the reaction enthalpy of the second stage was low, indicating that the reduction of the second ABTS•+ molecule occurred immediately after the reduction of the first one. For comparison, the IP and the enthalpy of a hydrogen atom abstraction of the phenol group of Trolox, calculated in ethanol, were 98.0 and 73.0 kcal/mol, respectively.
A noticeably higher AOA of compound 1a with methoxy substituents compared to other compounds 1 in the FRAP test and the similar AOA in the entire series of compounds 1 in the ABTS test could be explained by the different access of antioxidant molecules to the ABTS•+ radical versus the [Fe (TPTZ)2]3+ cation. Indeed, [Fe (TPTZ)2]3+, with a three-dimensional propeller structure, is bulkier than the quasi-one-dimensional structure of ABTS•+ (Supplementary Figure S1), which hinders access for bulky antioxidants.
3.10 Prediction of ADMET, physicochemical, and PAINS profiles
Table 3 shows the predicted ADMET properties and physico-chemical characteristics that we computed for compounds 1a–f and 2a–f. As can be seen, with the exception of the phenethyl and chlorophenethyl derivatives 1e–f and 2e–f, we obtained favorable values for intestinal absorption, indicating that the compounds would be expected to be bioavailable when administered orally. Note that it is possible that the unfavorable numbers for compounds 1e–f and 2e–f were likely due to their being outside of the applicability domain of the model. Given the high BBB permeability values, the compounds would be expected to enter the brain to exert effects on CNS targets; this is especially so for the benzyl and phenethyl derivatives. Parameters for cardiac toxicity risk, hERG pKi and pIC50, were within 5.0–7.9 log units, near the center of the possible range of 3–9 log units. Lipophilicities, as assessed by computed LogP values, were in some cases outside of the classical Lipinsky rule of 5 range, but these values were deemed unreliable owing to their being outside of the applicability domain of the model. Otherwise, the molecular weights and predicted aqueous solubility values were near or within acceptable limits for drug-like molecules. Integral quantitative estimates of drug likeness (QED) fell within a range of 0.2–0.9; the most active compounds, 1d,e, scored 0.3 for this parameter. Finally, no alerts were identified by the Pan Assay INterference compoundS (PAINS) filter.
Taken together, the results of the computed ADMET, physico-chemical, and PAINS profiles indicate that the lead compounds were acceptable for being advanced to further optimization and experimental studies of efficacy and safety.
4 Conclusion
Thus, studying the esterase profile of two new groups of acridine derivatives, 9-phosphoryl-9,10-dihydroacridines 1 and 9-phosphoryl-acridines 2, showed that all new compounds weakly inhibited AChE. Dihydroacridines with aryl substituents in the phosphoryl fragment inhibited BChE, with dibenzyloxy derivative 1d and its diphenethyl bioisostere 1e being the most active. Among acridines 2, only one compound 2d, an aromatized analog of 1d, proved to be an effective BChE inhibitor. Using 9-bis(benzyloxy)phosphoryl-9,10-dihydroacridine 1d as an example, a mixed mechanism of BChE inhibition was demonstrated from a kinetics analysis. The results of molecular docking served to explain the experimentally determined parameters of the potency, selectivity, and mechanism of cholinesterases inhibition by the acridine derivatives. Weak inhibition of CES, an off-target esterase, is a positive characteristic of the studied compounds, because this reduces the likelihood of unwanted drug-drug interactions in the case of their use in AD therapy.
The ability of the studied acridine derivatives to inhibit the self-aggregation of Aβ42 depended on both the structure of the acridine moiety and the structure of the substituents at the phosphorus atom. Moreover, the structure-activity relationship largely coincided with that observed for BChE inhibition. Dihydroacridines 1, especially with aryl substituents in the phosphoryl moiety, were able to block Aβ42 self-aggregation. The most active were the dibenzyloxy derivative 1d and its bioisosteric analog, the diphenethyl derivative 1e. On the other hand, acridine derivatives 2 inhibited Aβ42 self-aggregation weakly or not at all. Such a difference may be connected to the different ability to self-associate through π-π stacking interactions of the planar aromatic structure of acridines compared to non-planar dihydroacridines, which could compete with their binding to the Aβ42 peptide.
All 9-phosphoryl-9,10-dihydroacridines 1 demonstrated high ABTS•+-scavenging and iron-reducing activity at the level of Trolox, unlike 9-phosphoryl-acridines 2, which were not active. Quantum-chemical calculations confirmed this difference.
Calculated ADMET parameters of the test compounds predicted favorable intestinal absorption, good blood–brain barrier permeability, and cardiac toxicity potential.
In summary, compounds 1d and 1e with dibenzyloxy- and diphenethyl- substituents were found to display high inhibition of BChE activity and Aβ42 self-aggregation, high antioxidant capability, and favorable predicted ADMET properties. Therefore, these molecules can be considered lead compounds and are recommended for further in-depth studies as promising anti-AD agents.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
GM, SL, OC, and IU contributed to conception and design of the study; AS, MA, IU, and ND synthesized the compounds; NK, ERu, NB, and OS performed experimental investigation of biological activity; SL, TA, ET, ERa, VP, and RR performed computer modeling; NK, AS, TA, SL, ERu, NB, and ERa wrote sections of the manuscript and the first draft of the manuscript; GM, IU, NK, SL, VP, VF, SB, and RR—review and editing of the manuscript; GM, IU, VC—supervision of the project. All authors contributed to the article and approved the submitted version.
Funding
This research was partly supported by grant # 22-13-00298 of the Russian Science Foundation and IPAC RAS State Targets Project # FFSN-2021-0005; quantum-chemical calculations were supported the IBCP RAS State Targets Project # 122041400110-4. The synthesis of the compounds was financially supported by the Russian Foundation for Basic Research (research project # 19-29-08037). Support for RR’s contributions to the computer modeling components of the work was provided in part by a grant from the Alternatives Research and Development Foundation (ARDF) and an Mcubed grant from the University of Michigan.
Acknowledgments
We thank the shared research facilities of the HPC computing resources at Lomonosov Moscow State University (Voevodin et al., 2019) and Joint Supercomputer Center of the Russian Academy of Sciences for supercomputer time, the “Centre for Collective Use of IPAC RAS” for use of equipment cited in the Methods. We are grateful to Prof. A.V. Nemukhin for helpful discussions.
Conflict of interest
RR currently serves as a member of the advisory board of NeuroX1, a startup biotechnology company that is developing a software platform for the discovery and development of drugs for neurodegenerative diseases.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1219980/full#supplementary-material
References
Adamo, C., and Barone, V. (1999). Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 110, 6158–6170. doi:10.1063/1.478522
Agatonovic-Kustrin, S., Kettle, C., and Morton, D. W. (2018). A molecular approach in drug development for Alzheimer's disease. Biomed. Pharmacother. 106, 553–565. doi:10.1016/j.biopha.2018.06.147
Antosova, A., Chelli, B., Bystrenova, E., Siposova, K., Valle, F., Imrich, J., et al. (2011). Structure-activity relationship of acridine derivatives to amyloid aggregation of lysozyme. Biochim. Biophys. Acta 1810, 465–474. doi:10.1016/j.bbagen.2011.01.007
Author Anonymous (2022). 2022 Alzheimer's disease facts and figures. Alzheimers Dement. 18, 700–789. doi:10.1002/alz.12638
Bachurin, S. O., Aksinenko, A. Y., Makhaeva, G. F., and Shevtsova, E. F. (2023). Multipharmacophore strategy in medicinal chemistry for the design of drugs for the treatment of Alzheimer’s and some other neurodegenerative diseases. Russ. Chem. Bull. 72, 130–147. doi:10.1007/s11172-023-3718-0
Bachurin, S. O., Bovina, E. V., and Ustyugov, A. A. (2017). Drugs in clinical trials for alzheimer's disease: The major trends. Med. Res. Rev. 37, 1186–1225. doi:10.1002/med.21434
Ballard, C., Greig, N., Guillozet-Bongaarts, A., Enz, A., and Darvesh, S. (2005). Cholinesterases: Roles in the brain during health and disease. Curr. Alzheimer Res. 2, 307–318. doi:10.2174/1567205054367838
Bartolini, M., Bertucci, C., Cavrini, V., and Andrisano, V. (2003). beta-Amyloid aggregation induced by human acetylcholinesterase: inhibition studies. Biochem. Pharmacol. 65, 407–416. doi:10.1016/s0006-2952(02)01514-9
Bartolini, M., Naldi, M., Fiori, J., Valle, F., Biscarini, F., Nicolau, D. V., et al. (2011). Kinetic characterization of amyloid-beta 1-42 aggregation with a multimethodological approach. Anal. Biochem. 414, 215–225. doi:10.1016/j.ab.2011.03.020
Benzie, I. F., and Strain, J. J. (1999). Ferric reducing/antioxidant power assay: Direct measure of total antioxidant activity of biological fluids and modified version for simultaneous measurement of total antioxidant power and ascorbic acid concentration. Meth. Enzymol. 299, 15–27. doi:10.1016/s0076-6879(99)99005-5
Benzie, I. F., and Strain, J. J. (1996). The ferric reducing ability of plasma (FRAP) as a measure of "antioxidant power": The FRAP assay. Anal. Biochem. 239, 70–76. doi:10.1006/abio.1996.0292
Berg, L., Andersson, C. D., Artursson, E., Hornberg, A., Tunemalm, A. K., Linusson, A., et al. (2011). Targeting acetylcholinesterase: Identification of chemical leads by high throughput screening, structure determination and molecular modeling. PLoS One 6, e26039. doi:10.1371/journal.pone.0026039
Bernthsen, A. (1884). Die acridine. Justus Liebigs Ann. Chem. 224, 1–56. doi:10.1002/jlac.18842240102
Biancalana, M., and Koide, S. (2010). Molecular mechanism of Thioflavin-T binding to amyloid fibrils. Biochim. Biophys. Acta 1804, 1405–1412. doi:10.1016/j.bbapap.2010.04.001
Bickerton, G. R., Paolini, G. V., Besnard, J., Muresan, S., and Hopkins, A. L. (2012). Quantifying the chemical beauty of drugs. Nat. Chem. 4, 90–98. doi:10.1038/nchem.1243
Bongarzone, S., and Bolognesi, M. L. (2011). The concept of privileged structures in rational drug design: Focus on acridine and quinoline scaffolds in neurodegenerative and protozoan diseases. Expert Opin. Drug Discov. 6, 251–268. doi:10.1517/17460441.2011.550914
Bourne, Y., Radic, Z., Taylor, P., and Marchot, P. (2010). Conformational remodeling of femtomolar inhibitor-acetylcholinesterase complexes in the crystalline state. J. Am. Chem. Soc. 132, 18292–18300. doi:10.1021/ja106820e
Buurma, N. J., and Haq, I. (2008). Calorimetric and spectroscopic studies of hoechst 33258: Self-association and binding to non-cognate DNA. J. Mol. Biol. 381, 607–621. doi:10.1016/j.jmb.2008.05.073
Calabro, M., Rinaldi, C., Santoro, G., and Crisafulli, C. (2021). The biological pathways of alzheimer disease: A review. AIMS Neurosci. 8, 86–132. doi:10.3934/Neuroscience.2021005
Carreiras, M. C., Mendes, E., Perry, M. J., Francisco, A. P., and Marco-Contelles, J. (2013). The multifactorial nature of Alzheimer's disease for developing potential therapeutics. Curr. Top. Med. Chem. 13, 1745–1770. doi:10.2174/15680266113139990135
Castellani, R. J., Plascencia-Villa, G., and Perry, G. (2019). The amyloid cascade and alzheimer's disease therapeutics: Theory versus observation. Lab. Investig. 99, 958–970. doi:10.1038/s41374-019-0231-z
Charushin, V. N., and Chupakhin, O. N. (2004). SNH methodology and new approaches to condensed heterocyclic systems. Pure Appl. Chem. 76, 1621–1631. doi:10.1351/pac200476091621
Charushin, V. N., and Chupakhin, O. N. (2005). SNH methodology and new approaches to condensed heterocyclic systems. ChemInform 36. doi:10.1002/chin.200509263
Cheung, J., Rudolph, M. J., Burshteyn, F., Cassidy, M. S., Gary, E. N., Love, J., et al. (2012). Structures of human acetylcholinesterase in complex with pharmacologically important ligands. J. Med. Chem. 55, 10282–10286. doi:10.1021/jm300871x
Chufarova, N., Czarnecka, K., Skibinski, R., Cuchra, M., Majsterek, I., and Szymanski, P. (2018). New tacrine-acridine hybrids as promising multifunctional drugs for potential treatment of Alzheimer's disease. Arch. Pharm. Weinh. 351, e1800050. doi:10.1002/ardp.201800050
Chupakhin, O. N., and Charushin, V. N. (2017). Nucleophilic C–H functionalization of arenes: A new logic of organic synthesis. Pure Appl. Chem. 89, 1195–1208. doi:10.1515/pac-2017-0108
Chupakhin, O. N., Charushin, V. N., and Van Der Plas, H. C. (1994). Nucleophilic aromatic substitution of hydrogen. San Diego: Academic Press, INC. California 92101-4495.
Chupakhin, O. N., Shchepochkin, A. V., and Charushin, V. N. (2017). Atom- and step-economical nucleophilic arylation of azaaromatics via electrochemical oxidative cross C–C coupling reactions. Green Chem. 19, 2931–2935. doi:10.1039/c7gc00789b
Collinge, J., Gorham, M., Hudson, F., Kennedy, A., Keogh, G., Pal, S., et al. (2009). Safety and efficacy of quinacrine in human prion disease (PRION-1 study): A patient-preference trial. Lancet Neurol. 8, 334–344. doi:10.1016/s1474-4422(09)70049-3
Crescenzi, O., Tomaselli, S., Guerrini, R., Salvadori, S., D'ursi, A. M., Temussi, P. A., et al. (2002). Solution structure of the Alzheimer amyloid beta-peptide (1-42) in an apolar microenvironment. Similarity with a virus fusion domain. Eur. J. Biochem. 269, 5642–5648. doi:10.1046/j.1432-1033.2002.03271.x
Crismon, M. L. (1994). Tacrine: First drug approved for alzheimer's disease. Ann. Pharmacother. 28, 744–751. doi:10.1177/106002809402800612
Csuk, R., Barthel, A., Raschke, C., Kluge, R., Strohl, D., Trieschmann, L., et al. (2009). Synthesis of monomeric and dimeric acridine compounds as potential therapeutics in Alzheimer and prion diseases. Arch. Pharm. Weinh. 342, 699–709. doi:10.1002/ardp.200900065
Cummings, J., Lee, G., Ritter, A., Sabbagh, M., and Zhong, K. (2019). Alzheimer's disease drug development pipeline: 2019. Alzheimers Dement. (N Y) 5, 272–293. doi:10.1016/j.trci.2019.05.008
Demkowicz, S., Rachon, J., Daśko, M., and Kozak, W. (2016). Selected organophosphorus compounds with biological activity. Applications in medicine. RSC Adv. 6, 7101–7112. doi:10.1039/c5ra25446a
Denny, W. A. (2002). Acridine derivatives as chemotherapeutic agents. Curr. Med. Chem. 9, 1655–1665. doi:10.2174/0929867023369277
Dyabina, A. S., Radchenko, E. V., Palyulin, V. A., and Zefirov, N. S. (2016). Prediction of blood-brain barrier permeability of organic compounds. Dokl. Biochem. Biophys. 470, 371–374. doi:10.1134/S1607672916050173
Evstigneev, M. P., Evstigneev, V. P., and Davies, D. B. (2006). 1H NMR determination of the self-association of an acridine homodimer and its complexation with ethidium bromide in aqueous solution. J. Biomol. Struct. 784, 162–168. doi:10.1016/j.molstruc.2005.08.033
Florin Danet, A. (2021). “Recent advances in antioxidant capacity assays,” in Antioxidants - benefits, sources, mechanisms of action. Editor V. Waisundara (London: IntechOpen), 1–35.
Franco Pinto, J., Fillion, A., Duchambon, P., Bombard, S., and Granzhan, A. (2022). Acridine-O(6)-benzylguanine hybrids: Synthesis, DNA binding, MGMT inhibition and antiproliferative activity. Eur. J. Med. Chem. 227, 113909. doi:10.1016/j.ejmech.2021.113909
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2016). Gaussian 16 revision C 01 2016. Wallingford, CT, USA: Gaussian Inc.
Gamage, S. A., Figgitt, D. P., Wojcik, S. J., Ralph, R. K., Ransijn, A., Mauel, J., et al. (1997). Structure-activity relationships for the antileishmanial and antitrypanosomal activities of 1'-substituted 9-anilinoacridines. J. Med. Chem. 40, 2634–2642. doi:10.1021/jm970232h
Gao, Q., Williams, L. D., Egli, M., Rabinovich, D., Chen, S. L., Quigley, G. J., et al. (1991). Drug-induced DNA repair: X-Ray structure of a DNA-ditercalinium complex. Proc. Natl. Acad. Sci. U. S. A. 88, 2422–2426. doi:10.1073/pnas.88.6.2422
Girault, S., Grellier, P., Berecibar, A., Maes, L., Mouray, E., Lemière, P., et al. (2000). Antimalarial, antitrypanosomal, and antileishmanial activities and cytotoxicity of bis(9-amino-6-chloro-2-methoxyacridines): Influence of the linker. J. Med. Chem. 43, 2646–2654. doi:10.1021/jm990946n
Gulcin, I. (2020). Antioxidants and antioxidant methods: An updated overview. Arch. Toxicol. 94, 651–715. doi:10.1007/s00204-020-02689-3
Hampel, H., Mesulam, M. M., Cuello, A. C., Khachaturian, A. S., Vergallo, A., Farlow, M. R., et al. (2019). Revisiting the cholinergic hypothesis in alzheimer's disease: Emerging evidence from translational and clinical research. J. Prev. Alzheimers Dis. 6, 2–15. doi:10.14283/jpad.2018.43
Huang, Y., and Mucke, L. (2012). Alzheimer mechanisms and therapeutic strategies. Cell. 148, 1204–1222. doi:10.1016/j.cell.2012.02.040
Humphrey, W., Dalke, A., and Schulten, K. (1996). Vmd: Visual molecular dynamics. J. Mol. Graph. 14, 33–38, 27-28. doi:10.1016/0263-7855(96)00018-5
Ilyasov, I. R., Beloborodov, V. L., Selivanova, I. A., and Terekhov, R. P. (2020). ABTS/PP decolorization assay of antioxidant capacity reaction pathways. Int. J. Mol. Sci. 21, 1131. doi:10.3390/ijms21031131
Ismaili, L., and Romero, A. (2017). “Multitarget-directed antioxidants as therapeutic agents: Putting the focus on the oxidative stress,” in Design of hybrid molecules for drug development. Editor M. Decker (Elsevier Ltd), 5–47.
Ivanova, A., Gerasimova, E., and Gazizullina, E. (2020). Study of antioxidant properties of agents from the perspective of their action mechanisms. Molecules 25, 4251. doi:10.3390/molecules25184251
Jeremic, D., Jiménez-Díaz, L., and Navarro-López, J. D. (2021). Past, present and future of therapeutic strategies against amyloid-β peptides in alzheimer's disease: A systematic review. Ageing Res. Rev. 72, 101496. doi:10.1016/j.arr.2021.101496
Jin, H. (2014). Acetylcholinesterase and butyrylcholinesterase inhibitory properties of functionalized tetrahydroacridines and related analogs. Med. Chem. 4. doi:10.4172/2161-0444.1000213
Klingenstein, R., Lober, S., Kujala, P., Godsave, S., Leliveld, S. R., Gmeiner, P., et al. (2006). Tricyclic antidepressants, quinacrine and a novel, synthetic chimera thereof clear prions by destabilizing detergent-resistant membrane compartments. J. Neurochem. 98, 748–759. doi:10.1111/j.1471-4159.2006.03889.x
Koellner, G., Kryger, G., Millard, C. B., Silman, I., Sussman, J. L., and Steiner, T. (2000). Active-site gorge and buried water molecules in crystal structures of acetylcholinesterase from Torpedo californica. J. Mol. Biol. 296, 713–735. doi:10.1006/jmbi.1999.3468
Korth, C., May, B. C., Cohen, F. E., and Prusiner, S. B. (2001). Acridine and phenothiazine derivatives as pharmacotherapeutics for prion disease. Proc. Natl. Acad. Sci. U. S. A. 98, 9836–9841. doi:10.1073/pnas.161274798
Kozurkova, M., Hamulakova, S., Gazova, Z., Paulikova, H., and Kristian, P. (2011). Neuroactive multifunctional tacrine congeners with cholinesterase, anti-amyloid aggregation and neuroprotective properties. Pharmaceuticals 4, 382–418. doi:10.3390/ph4020382
Laikov, D. N., and Ustynyuk, Y. A. (2005). PRIRODA-04: A quantum-chemical program suite. New possibilities in the study of molecular systems with the application of parallel computing. Russ. Chem. Bull. 54, 820–826. doi:10.1007/s11172-005-0329-x
Lane, R. M., Potkin, S. G., and Enz, A. (2006). Targeting acetylcholinesterase and butyrylcholinesterase in dementia. Int. J. Neuropsychopharmacol. 9, 101–124. doi:10.1017/S1461145705005833
LeVine, H. (1999). Quantification of beta-sheet amyloid fibril structures with thioflavin T. Meth. Enzymol. 309, 274–284. doi:10.1016/s0076-6879(99)09020-5
Liu, Z., Dong, N., Xu, M., Sun, Z., and Tu, T. (2013). Mild Negishi cross-coupling reactions catalyzed by acenaphthoimidazolylidene palladium complexes at low catalyst loadings. J. Org. Chem. 78, 7436–7444. doi:10.1021/jo400803s
Lotfi, S., Rahmani, T., Hatami, M., Pouramiri, B., Kermani, E. T., Rezvannejad, E., et al. (2020). Design, synthesis and biological assessment of acridine derivatives containing 1,3,4-thiadiazole moiety as novel selective acetylcholinesterase inhibitors. Bioorg. Chem. 105, 104457. doi:10.1016/j.bioorg.2020.104457
Löwdin, P.-O. (1970). “On the nonorthogonality problem,” in Advances in quantum chemistry. Editor L. Per-Olov (New York, London: Academic Press), 185–199.
Lushchekina, S. V., and Masson, P. (2020). Slow-binding inhibitors of acetylcholinesterase of medical interest. Neuropharmacology 177, 108236. doi:10.1016/j.neuropharm.2020.108236
Lyakhov, S. A., Suveyzdis, Y. I., Litvinova, L. A., Andronati, S. A., Rybalko, S. L., and Dyadyun, S. T. (2000). Biological active acridine derivatives. Part 4: Synthesis and antiviral activity of some bis-acridinylated diamides. Pharmazie 55, 733–736.
Macdonald, I. R., Rockwood, K., Martin, E., and Darvesh, S. (2014). Cholinesterase inhibition in alzheimer's disease: Is specificity the answer? J. Alzheimers Dis. 42, 379–384. doi:10.3233/JAD-140219
Maciejewska, K., Czarnecka, K., Krecisz, P., Niedzialek, D., Wieczorek, G., Skibinski, R., et al. (2022). Novel cyclopentaquinoline and acridine analogs as multifunctional, potent drug candidates in Alzheimer’s disease. Int. J. Mol. Sci. 23, 5876. doi:10.3390/ijms23115876
Makhaeva, G. F., Elkina, N. A., Shchegolkov, E. V., Boltneva, N. P., Lushchekina, S. V., Serebryakova, O. G., et al. (2019a). Synthesis, molecular docking, and biological evaluation of 3-oxo-2-tolylhydrazinylidene-4,4,4-trifluorobutanoates bearing higher and natural alcohol moieties as new selective carboxylesterase inhibitors. Bioorg. Chem. 91, 103097. doi:10.1016/j.bioorg.2019.103097
Makhaeva, G. F., Kovaleva, N. V., Rudakova, E. V., Boltneva, N. P., Grishchenko, M. V., Lushchekina, S. V., et al. (2023). Conjugates of tacrine and salicylic acid derivatives as new promising multitarget agents for Alzheimer’s disease. Int. J. Mol. Sci. 24, 2285. doi:10.3390/ijms24032285
Makhaeva, G. F., Kovaleva, N. V., Rudakova, E. V., Boltneva, N. P., Lushchekina, S. V., Faingold, I. I., et al. (2020). New multifunctional agents based on conjugates of 4-Amino-2,3-polymethylenequinoline and butylated hydroxytoluene for Alzheimer’s disease treatment. Molecules 25, 5891. doi:10.3390/molecules25245891
Makhaeva, G. F., Lushchekina, S. V., Boltneva, N. P., Serebryakova, O. G., Rudakova, E. V., Ustyugov, A. A., et al. (2017). 9-Substituted acridine derivatives as acetylcholinesterase and butyrylcholinesterase inhibitors possessing antioxidant activity for Alzheimer's disease treatment. Bioorg. Med. Chem. 25, 5981–5994. doi:10.1016/j.bmc.2017.09.028
Makhaeva, G. F., Lushchekina, S. V., Kovaleva, N. V., Astakhova, T. Y., Boltneva, N. P., Rudakova, E. V., et al. (2021). Amiridine-piperazine hybrids as cholinesterase inhibitors and potential multitarget agents for Alzheimer's disease treatment. Bioorg. Chem. 112, 104974. doi:10.1016/j.bioorg.2021.104974
Makhaeva, G. F., Proshin, A. N., Kovaleva, N. V., Rudakova, E. V., Boltneva, N. P., Lushchekina, S. V., et al. (2022). Synthesis and study of the biological activity of thiourea-containing amiridine derivatives as potential multi-target drugs for the treatment of Alzheimer’s disease. Russ. Chem. Bull. 71, 2404–2415. doi:10.1007/s11172-022-3668-y
Makhaeva, G. F., Rudakova, E. V., Kovaleva, N. V., Lushchekina, S. V., Boltneva, N. P., Proshin, A. N., et al. (2019b). Cholinesterase and carboxylesterase inhibitors as pharmacological agents. Russ. Chem. Bull. 68, 967–984. doi:10.1007/s11172-019-2507-2
Makhaeva, G. F., Rudakova, E. V., Serebryakova, O. G., Aksinenko, A. Y., Lushchekina, S. V., Bachurin, S. O., et al. (2016). Esterase profiles of organophosphorus compounds in vitro predict their behavior in vivo. Chem. Biol. Interact. 259, 332–342. doi:10.1016/j.cbi.2016.05.002
Makhaeva, G. F., Shevtsova, E. F., Boltneva, N. P., Lushchekina, S. V., Kovaleva, N. V., Rudakova, E. V., et al. (2019c). Overview of novel multifunctional agents based on conjugates of gamma-carbolines, carbazoles, tetrahydrocarbazoles, phenothiazines, and aminoadamantanes for treatment of Alzheimer's disease. Chem. Biol. Interact. 308, 224–234. doi:10.1016/j.cbi.2019.05.020
Makosza, M. (2010). Nucleophilic substitution of hydrogen in electron-deficient arenes, a general process of great practical value. Chem. Soc. Rev. 39, 2855–2868. doi:10.1039/b822559c
Mąkosza, M. (2011). Nucleophilic substitution of hydrogen in nitroarenes: A new chapter of aromatic chemistry. Synthesis 2011, 2341–2356. doi:10.1055/s-0030-1260668
Mąkosza, M., and Wojciechowski, K. (2014). Synthesis of heterocycles via nucleophilic substitution of hydrogen in nitroarenes. Heterocycles 88, 75. doi:10.3987/rev-13-sr(s)1
Mallu, L., Thirumalai, D., and Asharani, I. V. (2017). One-pot cascade synthesis and in-vitro evaluation of anti-inflammatory and anti-diabetic activities of S-methylphenyl substituted acridine-1,8-diones. Chem. Biol. Drug Des. 90, 520–526. doi:10.1111/cbdd.12973
Marenich, A. V., Cramer, C. J., and Truhlar, D. G. (2009). Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396. doi:10.1021/jp810292n
Martinez, A., and Castro, A. (2006). Novel cholinesterase inhibitors as future effective drugs for the treatment of Alzheimer's disease. Expert Opin. Investig. Drugs 15, 1–12. doi:10.1517/13543784.15.1.1
Masson, P., Lushchekina, S., Schopfer, L. M., and Lockridge, O. (2013). Effects of viscosity and osmotic stress on the reaction of human butyrylcholinesterase with cresyl saligenin phosphate, a toxicant related to aerotoxic syndrome: Kinetic and molecular dynamics studies. Biochem. J. 454, 387–399. doi:10.1042/bj20130389
Moreira, P. I., Siedlak, S. L., Aliev, G., Zhu, X., Cash, A. D., Smith, M. A., et al. (2005). Oxidative stress mechanisms and potential therapeutics in Alzheimer disease. J. Neural Transm. 112, 921–932. doi:10.1007/s00702-004-0242-8
Morris, G. M., Goodsell, D. S., Halliday, R. S., Huey, R., Hart, W. E., Belew, R. K., et al. (1998). Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 19, 1639–1662. doi:10.1002/(sici)1096-987x(19981115)19:14<1639:aid-jcc10>3.0.co;2-b
Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and AutoDock Tools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi:10.1002/jcc.21256
Munoz-Ruiz, P., Rubio, L., Garcia-Palomero, E., Dorronsoro, I., Del Monte-Millan, M., Valenzuela, R., et al. (2005). Design, synthesis, and biological evaluation of dual binding site acetylcholinesterase inhibitors: New disease-modifying agents for alzheimer's disease. J. Med. Chem. 48, 7223–7233. doi:10.1021/jm0503289
Nascimento, E. C. M., Oliva, M., Swiderek, K., Martins, J. B. L., and Andres, J. (2017). Binding analysis of some classical acetylcholinesterase inhibitors: Insights for a rational design using free energy perturbation method calculations with QM/MM MD simulations. J. Chem. Inf. Model. 57, 958–976. doi:10.1021/acs.jcim.7b00037
Nicolet, Y., Lockridge, O., Masson, P., Fontecilla-Camps, J. C., and Nachon, F. (2003). Crystal structure of human butyrylcholinesterase and of its complexes with substrate and products. J. Biol. Chem. 278, 41141–41147. doi:10.1074/jbc.M210241200
Ortiz, G. G., Pacheco Moisés, F. P., Mireles-Ramírez, M., Flores-Alvarado, L. J., González-Usigli, H., Sánchez-González, V. J., et al. (2017). “Oxidative stress: Love and hate history in central nervous system,” in Adv. Protein chem. Struct. Biol. (Academic Press).
Pan, X., Wang, H., Li, C., Zhang, J. Z. H., and Ji, C. (2021). MolGpka: A web server for small molecule pKa prediction using a graph-convolutional neural network. J. Chem. Inf. Model. 61, 3159–3165. doi:10.1021/acs.jcim.1c00075
Papajak, E., Zheng, J., Xu, X., Leverentz, H. R., and Truhlar, D. G. (2011). Perspectives on basis sets beautiful: Seasonal plantings of diffuse basis functions. J. Chem. Theory Comput. 7, 3027–3034. doi:10.1021/ct200106a
Phillips, J. C., Hardy, D. J., Maia, J. D. C., Stone, J. E., Ribeiro, J. V., Bernardi, R. C., et al. (2020). Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 153, 044130. doi:10.1063/5.0014475
Pohanka, M. (2018). Oxidative stress in Alzheimer disease as a target for therapy. Bratisl. Lek. Listy 119, 535–543. doi:10.4149/BLL_2018_097
Przybylowska, M., Dzierzbicka, K., Kowalski, S., Demkowicz, S., Dasko, M., and Inkielewicz-Stepniak, I. (2022). Design, synthesis and biological evaluation of novel N-phosphorylated and O-phosphorylated tacrine derivatives as potential drugs against Alzheimer's disease. J. Enzyme Inhib. Med. Chem. 37, 1012–1022. doi:10.1080/14756366.2022.2045591
QSAR (2022). ADMET prediction service. http://qsar.chem.msu.ru/admet/ (Accessed on June 15, 2022).
Radchenko, E. V., Dyabina, A. S., and Palyulin, V. A. (2020). Towards deep neural network models for the prediction of the blood-brain barrier permeability for diverse organic compounds. Molecules 25, 5901. doi:10.3390/molecules25245901
Radchenko, E. V., Dyabina, A. S., Palyulin, V. A., and Zefirov, N. S. (2016). Prediction of human intestinal absorption of drug compounds. Russ. Chem. Bull. 65, 576–580. doi:10.1007/s11172-016-1340-0
Radchenko, E. V., Rulev, Y. A., Safanyaev, A. Y., Palyulin, V. A., and Zefirov, N. S. (2017). Computer-aided estimation of the hERG-mediated cardiotoxicity risk of potential drug components. Dokl. Biochem. Biophys. 473, 128–131. doi:10.1134/S1607672917020107
Rassolov, V. A., Ratner, M. A., Pople, J. A., Redfern, P. C., and Curtiss, L. A. (2001). 6-31G* basis set for third-row atoms. J. Comput. Chem. 22, 976–984. doi:10.1002/jcc.1058
RDKit (2022). RDKit: Open-source cheminformatics software. [Online]. Available: http://www.rdkit.org (Accessed June 15, 2022).
Re, R., Pellegrini, N., Proteggente, A., Pannala, A., Yang, M., and Rice-Evans, C. (1999). Antioxidant activity applying an improved ABTS radical cation decolorization assay. Free Radic. Biol. Med. 26, 1231–1237. doi:10.1016/s0891-5849(98)00315-3
Reardon, S. (2023). FDA approves Alzheimer's drug lecanemab amid safety concerns. Nature 613, 227–228. doi:10.1038/d41586-023-00030-3
Rosini, M., Simoni, E., Bartolini, M., Tarozzi, A., Matera, R., Milelli, A., et al. (2011). Exploiting the lipoic acid structure in the search for novel multitarget ligands against Alzheimer's disease. Eur. J. Med. Chem. 46, 5435–5442. doi:10.1016/j.ejmech.2011.09.001
Saxena, A., Redman, A. M., Jiang, X., Lockridge, O., and Doctor, B. P. (1997). Differences in active site gorge dimensions of cholinesterases revealed by binding of inhibitors to human butyrylcholinesterase. Biochemistry 36, 14642–14651. doi:10.1021/bi971425+
Schäfer, A., Huber, C., and Ahlrichs, R. (1994). Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr. J. Chem. Phys. 100, 5829–5835. doi:10.1063/1.467146
Schmidt, M. W., Baldridge, K. K., Boatz, J. A., Elbert, S. T., Gordon, M. S., Jensen, J. H., et al. (1993). General atomic and molecular electronic structure system. J. Comput. Chem. 14, 1347–1363. doi:10.1002/jcc.540141112
Sharma, A., and Piplani, P. (2017). Design and synthesis of some acridine-piperazine hybrids for the improvement of cognitive dysfunction. Chem. Biol. Drug Des. 90, 926–935. doi:10.1111/cbdd.13017
Shchepochkin, A. V., Chupakhin, O. N., Demina, N. S., Averkov, M. A., Shimanovskaya, T. Y., Slepukhin, P. A., et al. (2021)., 53. Synthesis, 3791–3798. doi:10.1055/A-1521-3166Electrochemical aromatization of dihydroazines: Effect of chalcogenophosphoryl (CGP) substituents on anodic oxidation of 9-CGP-9,10-dihydroacridine
Sondhi, S. M., Singh, J., Rani, R., Gupta, P. P., Agrawal, S. K., and Saxena, A. K. (2010). Synthesis, anti-inflammatory and anticancer activity evaluation of some novel acridine derivatives. Eur. J. Med. Chem. 45, 555–563. doi:10.1016/j.ejmech.2009.10.042
Spiegel, M., Kapusta, K., Kolodziejczyk, W., Saloni, J., Zbikowska, B., Hill, G. A., et al. (2020). Antioxidant activity of selected phenolic acids-ferric reducing antioxidant power assay and QSAR analysis of the structural features. Molecules 25, 3088. doi:10.3390/molecules25133088
Sushko, I., Novotarskyi, S., Korner, R., Pandey, A. K., Rupp, M., Teetz, W., et al. (2011). Online chemical modeling environment (OCHEM): Web platform for data storage, model development and publishing of chemical information. J. Comput. Aided. Mol. Des. 25, 533–554. doi:10.1007/s10822-011-9440-2
Sussman, J. L., Harel, M., Frolow, F., Oefner, C., Goldman, A., Toker, L., et al. (1991). Atomic structure of acetylcholinesterase from Torpedo californica: A prototypic acetylcholine-binding protein. Science 253, 872–879. doi:10.1126/science.1678899
Tseng, H. J., Lin, M. H., Shiao, Y. J., Yang, Y. C., Chu, J. C., Chen, C. Y., et al. (2020). Synthesis and biological evaluation of acridine-based histone deacetylase inhibitors as multitarget agents against Alzheimer's disease. Eur. J. Med. Chem. 192, 112193. doi:10.1016/j.ejmech.2020.112193
Tsurkan, L. G., Hatfield, M. J., Edwards, C. C., Hyatt, J. L., and Potter, P. M. (2013). Inhibition of human carboxylesterases hCE1 and hiCE by cholinesterase inhibitors. Chem. Biol. Interact. 203, 226–230. doi:10.1016/j.cbi.2012.10.018
Van Der Westhuizen, C. J., Stander, A., Riley, D. L., and Panayides, J. L. (2022). Discovery of novel acetylcholinesterase inhibitors by virtual screening, in vitro screening, and molecular dynamics simulations. J. Chem. Inf. Model. 62, 1550–1572. doi:10.1021/acs.jcim.1c01443
Van Dyck, C. H., Swanson, C. J., Aisen, P., Bateman, R. J., Chen, C., Gee, M., et al. (2023). Lecanemab in early alzheimer's disease. N. Engl. J. Med. 388, 9–21. doi:10.1056/NEJMoa2212948
Vanommeslaeghe, K., Hatcher, E., Acharya, C., Kundu, S., Zhong, S., Shim, J., et al. (2010). CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 31, 671–690. doi:10.1002/jcc.21367
Vaz, M., Silva, V., Monteiro, C., and Silvestre, S. (2022). Role of aducanumab in the treatment of alzheimer's disease: Challenges and opportunities. Clin. Interv. Aging 17, 797–810. doi:10.2147/CIA.S325026
Voevodin, V., Antonov, A., Nikitenko, D., Shvets, P., Sobolev, S., Sidorov, I., et al. (2019). Supercomputer lomonosov-2: Large scale, deep monitoring and fine analytics for the user community. Supercomput. Front. Innov. 6, 4–11. doi:10.14529/jsfi190201
Volkov, P. A., Khrapova, K. O., Telezhkin, A. A., Ivanova, N. I., Albanov, A. I., Gusarova, N. K., et al. (2018). Catalyst-Free phosphorylation of acridine with secondary phosphine chalcogenides: Nucleophilic addition vs SN(H)Ar reaction. Org. Lett. 20, 7388–7391. doi:10.1021/acs.orglett.8b03061
Wainwright, M. (2001). Acridine-a neglected antibacterial chromophore. J. Antimicrob. Chemother. 47, 1–13. doi:10.1093/jac/47.1.1
Walsh, D. M., Klyubin, I., Shankar, G. M., Townsend, M., Fadeeva, J. V., Betts, V., et al. (2005). The role of cell-derived oligomers of Abeta in Alzheimer's disease and avenues for therapeutic intervention. Biochem. Soc. Trans. 33, 1087–1090. doi:10.1042/BST20051087
Wlodek, S. T., Clark, T. W., Scott, L. R., and Mccammon, J. A. (1997). Molecular dynamics of acetylcholinesterase dimer complexed with tacrine. J. Am. Chem. Soc. 119, 9513–9522. doi:10.1021/ja971226d
Yanai, T., Tew, D. P., and Handy, N. C. (2004). A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 393, 51–57. doi:10.1016/j.cplett.2004.06.011
Keywords: 9-phosphoryl-9,10-dihydroacridines, 9-phosphorylacridines, acetylcholinesterase, butyrylcholinesterase, Alzheimer’s disease
Citation: Makhaeva GF, Kovaleva NV, Rudakova EV, Boltneva NP, Lushchekina SV, Astakhova TY, Timokhina EN, Serebryakova OG, Shchepochkin AV, Averkov MA, Utepova IA, Demina NS, Radchenko EV, Palyulin VA, Fisenko VP, Bachurin SO, Chupakhin ON, Charushin VN and Richardson RJ (2023) Derivatives of 9-phosphorylated acridine as butyrylcholinesterase inhibitors with antioxidant activity and the ability to inhibit β-amyloid self-aggregation: potential therapeutic agents for Alzheimer’s disease. Front. Pharmacol. 14:1219980. doi: 10.3389/fphar.2023.1219980
Received: 10 May 2023; Accepted: 18 July 2023;
Published: 09 August 2023.
Edited by:
Syed Nasir Abbas Bukhari, Jouf University, Saudi ArabiaReviewed by:
Jason Schwans, California State University, United StatesEhsan Ullah Mughal, University of Gujrat, Pakistan
Suleyman Akocak, Adiyaman University, Türkiye
Copyright © 2023 Makhaeva, Kovaleva, Rudakova, Boltneva, Lushchekina, Astakhova, Timokhina, Serebryakova, Shchepochkin, Averkov, Utepova, Demina, Radchenko, Palyulin, Fisenko, Bachurin, Chupakhin, Charushin and Richardson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rudy J. Richardson, rjrich@umich.edu