Next Article in Journal
Effects of Substituents on the Photophysical/Photobiological Properties of Mono-Substituted Corroles
Next Article in Special Issue
Synthesis of Disubstituted Carboxonium Derivatives of Closo-Decaborate Anion [2,6-B10H8O2CC6H5]−: Theoretical and Experimental Study
Previous Article in Journal
Polylactide-Based Materials: Synthesis and Biomedical Applications
Previous Article in Special Issue
GPU-Enhanced DFTB Metadynamics for Efficiently Predicting Free Energies of Biochemical Systems
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 3
10.3390/molecules28031388
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 Drug Design and Analysis of Dual Amyloid-Beta and Tau Protein-Aggregation Inhibitors for Alzheimer’s Disease Treatment

by
Nisha Job
1,
Venkatesan S. Thimmakondu
2,* and
Krishnan Thirumoorthy
1,*
1
Department of Chemistry, School of Advanced Sciences, Vellore Institute of Technology, Vellore 632014, Tamil Nadu, India
2
Department of Chemistry and Biochemistry, San Diego State University, San Diego, CA 92182, USA
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(3), 1388; https://doi.org/10.3390/molecules28031388
Submission received: 18 November 2022 / Revised: 21 December 2022 / Accepted: 8 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue Advances in the Theoretical and Computational Chemistry)
Browse Figures
Review Reports Versions Notes

Abstract

:
Alzheimer’s disease (AD) is a progressive and irreversible neurodegenerative disorder that gradually leads to the state of dementia. The main features of AD include the deposition of amyloid-beta peptides (Aβ), forming senile plaques, and the development of neurofibrillary tangles due to the accumulation of hyperphosphorylated Tau protein (p-tau) within the brain cells. In this report, seven dual-inhibitor molecules (L1–7) that can prevent the aggregation of both Aβ and p-tau are suggested. The drug-like features and identification of the target proteins are analyzed by the in silico method. L1–7 show positive results in both Blood–Brain Barrier (BBB) crossing and gastrointestinal absorption, rendering to the results of the permeation method. The molecular docking test performed for L1–7 shows binding energies in the range of −4.9 to −6.0 kcal/mol towards Aβ, and −4.6 to −5.6 kcal/mol for p-tau. The drug’s effectiveness under physiological conditions is assessed by the use of solvation models on the investigated systems. Further, the photophysical properties of L1–3 are predicted using TD-DFT studies.

1. Introduction

Alzheimer’s disease (AD) is a progressive and irreversible neurodegenerative disorder that slowly results in the state of dementia [1,2]. Several theories have been proposed to explain the root cause of AD, whilst the exact mechanism is still under debate [3,4]. The amyloid-β (Aβ) cascade hypothesis, cholinergic hypothesis, tau hypothesis, oxidative stress hypothesis, and inflammation hypothesis explain the possible causes of AD in different ways [5,6]. The main features of AD include the deposition of β-amyloid peptides forming senile plaques on the extracellular surface of neurons and the development of neurofibrillary tangles due to the accumulation of hyperphosphorylated Tau protein within the brain cells [7].
Senile plaques, which are made up of Aβ peptides, are found in AD patients. Aβ peptides are the proteolytic product of the amyloid precursor protein (APP). The proteolysis of the amyloid precursor protein produces Aβ-peptides. The splitting of APP in the N-terminus by β-secretase on the cell surface results in the generation of the soluble amyloid precursor protein β-fragment (APPs β). Additionally, the γ-secretase cleaves the C-terminus (C99), which is attached to the cell membrane, resulting in the formation of Aβ-peptides. The aggregation of these peptides causes the formation of insoluble fibrils, which eventually clump together to form the distinctive Aβ plaques as seen in AD patients. However, some other Aβ species mutations favor the longer peptides that terminate in residues 40 (Aβ40) and 42 (Aβ42). The main factor in the formation of senile plaques has been identified as Aβ42 because of its hydrophobic nature and the fact that it aggregates more quickly [8,9].
The p-tau protein is a microtubule-associated protein that is mostly found in the central nervous system (CNS). The functionality of the tau protein includes the stabilization of microtubules in the nerve cells. The alternative splicing of a specific gene called microtubule-associated protein tau (MAPT) produces tau proteins. Normal p-tau regulates gene transcription and cell-cycle activity, and is required for signaling molecules [10]. Under normal circumstances, tau is highly stable; however, different variables can trigger tau aggregation, resulting in neurofibrillary tangles [11]. The microtubule-binding region is known as the projection domain, which plays a vital role in structural stabilization and is generally flanked by tau-phosphosites. The hyperphosphorylation of p-tau causes the microtubule organization to be disrupted [12,13]. In the tau protein, around 80 phosphorylation sites for tyrosine (Tyr), serine (Ser), and threonine (Thr) exist. These sites contain a variety of Tyr, Ser, and Thr phosphatase residues that are involved in phosphorylation. Many kinases are responsible for the regulation of p-tau phosphorylation, including protein kinase A (PKA), protein kinase C (PKC), glycogen synthase kinase-3 (GSK3), and cyclin-dependent kinases (CDKs). In a healthy physiological state, there should be a balance between phosphate and kinase functions leading to normal tau behavior. Most of the hyperphosphorylation regions of p-tau are encoded within three sites (Tyr, Ser, and Thr) that also play a major role in AD progression [14,15]. Hence, the abnormalities of p-tau can be defined as a highly contributing factor to AD.
The current therapeutic strategies employed for the treatment of AD are centered on single-target medications, which are still unable to control the disease’s progression. Studies showed that both Aβ and p-tau influence the cognitive status by synergistically interfering with the neurological pathways. Multi-target drugs offer a single molecular entity to treat several AD-related variables, which has significantly grown in popularity as a potential AD cure [11]. We selected a p-tau inhibitor containing piperidine-substituted 3-nitrobenzaldehyde pharmacophore unit as a reference molecule from the literature, and tried out several combinations to make its derivatives multi-targeted for both Aβ and p-tau to improve their efficiency and to possess Blood–Brain Barrier (BBB) permeability [16]. We employed various in silico approaches, including molecular docking, molecular dynamics, and pharmacokinetic evaluations, for finding new lead molecules for AD treatment. These in silico methods can improve drug efficiency and reduce experimental costs as well. We have also considered only synthetically feasible molecules for conducting the in silico studies.
The first criterion for designing suitable drug molecules for AD treatment is to analyze their inhibitory potency towards Aβ and p-tau aggregation. The second criterion is BBB-crossing and gastrointestinal-absorption ability. Gastrointestinal absorption is vital in terms of drug delivery. For AD treatment, BBB permeability is essential to prevent the aggregation of different proteins inside the human brain. Herein, we designed multiple organic molecules as dual Aβ/p-tau aggregation inhibitors. Based on the available literature, the N-heterocycle substituted 3-nitrobenzaldehyde scaffold was chosen as the core structure. Since these molecules are achievable via different synthetic routes, their experimental analysis is also possible after theoretical confirmation. Each molecule was shortlisted based on target prediction and BBB permeability using SwissADME tools. The seven selected molecules shown in Figure 1 exhibit BBB crossing, higher Aβ, and p-tau binding affinity. We calculated the energies of frontier molecular orbitals (FMOs) in all the molecules to understand their global reactivity descriptors and molecular stability. The seven selected molecules’ 3D structures are depicted in Figure S1 (see Supplementary Materials). The aggregation of the proteins Aβ/p-tau is widely considered the root cause of AD. These seven chosen molecules can all selectively interact with these proteins, helping in their aggregation inhibition. This inhibitory action can only be achieved when the molecules bind in the particular protein position, which avoids aggregation. The primary reason for protein aggregation is the interaction between peptide residues of the protein. To understand these binding mechanisms, we analyze the interaction probability of the designed molecules towards the proteins using molecular docking. The best docking poses are further subjected to molecular dynamics simulation to understand the stability of these interactions. Further, molecules are analyzed using solvation models to mimic physiological conditions. Then, the TD-DFT studies of these molecules are carried out to check their photophysical properties under physiological conditions.

2. Results and Discussion

2.1. Prediction of Drug-Target and ADMET Properties

ADMET profile and binding efficiency are the cornerstones of the drug discovery journey. It is paramount to obtain the desired results in these studies for drug candidate molecules before entering the experimental phase [17,18]. In this work, computational tools were used to determine the ADMET features of the candidate molecules. The results obtained from the SwissADME tool are listed in Table 1. The total polar surface area (TPSA) is an important factor in terms of drugs’ bioavailability. The compounds with TPSA > 140 are passively absorbed and have very limited oral bioavailability. The TPSA value of the designed inhibitor molecules in the present work is in the range of 60–80, which is optimal in terms of bioavailability. The log P and log S values are used to measure the lipophilicity and solubility of the molecules that are in the range of −0.4 to 5.6 and −10 to 6, respectively [19]. All seven compounds were within the expected range, which indicates that they have sufficient absorbency and solubility. Skin permeability was evaluated by log Kp values, and the highest skin permeability was shown by L1 among the seven molecules. Lipinski’s rule [18] and Ghose’s rule [20] were used for evaluating the drug-likeness properties of the compounds. The rules showed zero violations against promoting the compounds as potential drug candidates. The toxicity data of the seven inhibitor molecules are listed in Table S1 (see Supplementary Materials). All the molecules are under the class 4 category with median lethal dose (LD50) values between 300 and 2000.
The potential drug candidates were identified by the SwissTargetPrediction. The results showed the inhibitory nature of the compounds toward both Aβ as well as p-tau. Both proteins are responsible for the progression of AD. The gastrointestinal-absorption and BBB-penetration data obtained from the BOILED-Egg method are given in Figure 2. The molecules are represented inside the BOILED-Egg model as either blue, which are substrates of P-glycoprotein, or red, which are not. All seven compounds are located in the yellow region and represented as a red circle in the model, which indicates the BBB-crossing and human gastrointestinal-absorption ability of these molecules. It is also clear that these molecules are not a substrate of P-glycoprotein, which indicates that these compounds cannot be actively effluxed from the CNS.

2.2. Frontier Molecular Orbital and Global Reactivity Descriptors

According to the widely recognized frontier molecular orbital (FMO) theory, chemical reactions are primarily driven by the frontier orbitals, specifically the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) [21]. It is considered a dependable resource for determining the reactivity of small organic molecules. The energy values, hardness, softness, electronegativity, chemical potential, electrophilicity index, and maximum charge transfer index are calculated using HOMO and LUMO [22] energy levels, and are included in Table 2.
The band gap between the HOMO and LUMO is predicted in the range of 6.9–7.6 eV, and all the molecules can be considered biocompatible as reported in the literature, since the band gap is greater than 2.5 eV [23]. L7 has the lowest band-gap value, highest chemical potential, and lowest hardness value. The low chemical softness (σ) and high electrophilicity index (ω) show low cytotoxicity and good biological activity [24]. The lowest σ value of L3 represents the lowest cytotoxicity of all the molecules, and the high ω value (3.3) of L5 shows the highest biological activity [24,25]. The maximum charge transfer index (ΔNmax) is the capacity of an electrophile that can accept a maximum electronic charge from the outer environment and is in the range of 1.1–1.4 for all the molecules [24].

2.3. Molecular Docking Analysis

As per the target prediction results, Aβ and p-tau have the probability to interact with molecules L1–7. Molecular docking studies help to understand the interacting position, affinity, and mechanism of the inhibitor molecules with the target proteins. Two sequence segments of p-tau with six amino-acid residues, 275VQIINK280 at the start of the R2 segment and 306VQIVYK311 at the start of R3, are responsible for p-tau aggregation, as reported in the literature [26]. Since the most potent driver of tau aggregation is the VQIINK segment, we performed docking studies on the 275VQIINK280 segment [27]. We selected the 274KVQIINKKLD283 structure with PDB ID 5V5B from the Research Collaboratory for Structural Bioinformatics (RCSB) for performing molecular docking [28]. In the case of Aβ peptide aggregation, Aβ42 comparatively aggregates faster than other peptides. Because of that, we performed docking studies on the Aβ1–43 segment with PDB ID 2LLM. Based on ten different docking positions, the best poses that have the highest binding score were chosen. The docked compounds and their binding affinity in the ligand-receptor complexes are listed in Table 3. L1–7 show binding energies in the range of −4.6 to −5.6 kcal/mol toward p-Tau, and −4.9 to −6.0 kcal/mol for Aβ.
L1 showed the highest docking score towards p-tau and Aβ of all the investigated molecules. Among the seven molecules, two (L1 and L3) showed a higher binding score with p-tau than when compared to the reference molecule. Figure 3 represents the best binding poses of L1–3 with both proteins. L1–3 molecules showed a higher binding score (−5.3 to −6.0 kcal/mol) with Aβ than other inhibitor molecules, and L2 showed a binding score of −5.1 kcal/mol with p-tau, which very closely trails the reference compound. In further studies, we shortlisted three of these best-performing molecules (L1–3), since the present work aims to find out the lead inhibitor molecule. Figure 4 and Table S2 (available in the Supplementary Materials) give a better insight into the interaction of the three selected molecules with Aβ and p-tau. The images have been generated with Pymol and the Discovery Studio visualizer program [29,30].
As seen in Figure 4, PHE8 is one of the most interactive amino-acid residues of Aβ that can bind through either π–π or π–alkyl interactions with selected inhibitor molecules. In L1, a hydrogen bond is formed between the oxygen atom in the molecule and LYS16 with a bond length of 2.20 Å, as listed in Table S2. L2 shows the second-strongest binding affinity towards Aβ with a score of −5.8, which involves the following interactions: π–σ in LEU5, π–π stacked in PHE8, π–alkyl in PHE8, and alkyl in VAL12.
In the case of p-tau, ASN279, LYS280, LYS281, and LEU282 are the amino acids interacting with the inhibitor molecules. In the three selected inhibitor molecules, the most interactive amino-acid residues are LYS280 and LYS281. Apart from these, ASN279 and LEU282 are also participating in the molecular interactions. L1 is found to be stable at the active site of p-tau by forming π–alkyl (LYS280 and 281), π–σ (LYS 281), and alkyl (LYS280 and LEU282) interactions. The carbonyl oxygen and imine nitrogen of L3 formed a hydrogen bond with ASN279 and LYS281, respectively. Carbon–hydrogen bonding (ASN279), π–alkyl (LYS280 and LYS281), π–σ (LYS281), and alkyl (LEU282) are responsible for the interactions in L3.

2.4. Molecular Dynamics Simulations

Molecular dynamics (MD) simulation is used to evaluate the flexibility or mobility of different parts of the biomolecule and helps to predict how a biomolecular system will react to a perturbation [31]. The best docking poses were subjected to molecular dynamics simulation using the CHARMM36m force field, and the topology of the ligands was generated using the CHARMM-GUI platform [32,33]. The Monte Carlo method is used for ion placing, and molecules are equilibrated using the NVT and NPT isothermal–isobaric ensemble. Subsequently, we executed a 1 ns MD simulation with a 2 fs integrative step. The root-mean-square deviation (RMSD) was used to analyze the conformational stability of biomolecules. The overall stability and structural convergence of the protein–inhibitor complexes along with the apo-protein system were compared. Figure 5 represents the RMSD plot of Aβ’s and p-tau’s protein–inhibitor complexes and their apo forms.
In the case of Aβ, apo and the inhibitor-bound forms deviate only to an acceptable extent, and the deviation range is only in the range of 0.1 to 1.5 nm. Among the selected three molecules (L1–3), the Aβ-L1 complex shows deviation from the apo-Aβ after 0.8 ns. The Aβ-L3 complex shows a slight variation from the apo-protein trend in the range of 0.2 to 0.5 ns. Apo-Aβ and the L3-bound forms show similar RMSD trends throughout the analysis. In the case of apo-p-tau, the RMSD plot of p-tau–L1 and p-tau–L3 showed a deviation in the range of 0.01 to 2.5 nm. The L2-bound p-tau showed a similar RMSD trend compared to others. After 0.7 ns, p-tau–L2 and p-tau–L3 deviate from the RMSD trend of the apo-p-tau. The L3 inhibitor molecule interacting with Aβ is in good agreement with its apo-protein RMSD trend. It also has a similar trend with apo-p-tau over the majority of the time-frame. The given RMSD data are acceptable because of a similar pattern of deviation observed even in the case of the apo-protein. The results show that L2 and L3 show similar RMSD deviation with apo-proteins with tau and Aβ, respectively. In the case of L1, both proteins show variation in the RMSD trends. On comparing L2 and L3 as dual inhibitors, L3 can be a better dual-inhibitor candidate compared to others. L3 shows a similar RMSD trend with apo-Aβ compared to the remaining molecules, while justifying the RMSD trend with p-tau deviates it to a certain extent.

2.5. TD-DFT Analysis and Solvent Effect

The TD-DFT method is used to determine the electronic properties of L1, L2, and L3 in their excited states. The wavelength (λ), the excitation energy, the electronic oscillator strength (f), and the corresponding electronic transition in the gaseous phase are all provided in Table S3 (available in the Supplementary Materials). The optimized geometries of L1, L2, and L3 revealed an intense absorption peak at 291, 299, and 249 nm, respectively. Regarding this absorption peak, L1 showed maximum contribution from the transition HOMO-1→LUMO+1 (40%) at 291 nm with f = 0.3055. In L2, the transitions are mainly attributed to the HOMO-1→LUMO+1 (46%), HOMO-1→LUMO+2 (13%), and HOMO-5→LUMO (13%), HOMO→LUMO (11%), respectively. In the case of L3, HOMO-4→LUMO+1, HOMO-4→LUMO, and HOMO-1→LUMO+2 transitions contribute 22%, 10%, and 8%, respectively. A value of f close to one denotes an allowed transition, whereas an f close to 0.001 denotes a quantum mechanically forbidden transition. The selected molecules show an f value between 0.2 and 0.3, which denotes allowed transitions. The optimized geometries in the water medium of L1, L2, and L3 revealed strong absorption peaks at 315, 333, and 232 nm, respectively (For detailed information, see the supplementary data Table S4). L2 showed the highest absorption peak at 333 nm with maximum contribution from HOMO-1→LUMO (63%). In the case of L1 and L2, the wavelength of the absorption peak is increased in the water medium compared to the gaseous phase. L3 had an intense absorption peak at 232 nm, with f = 0.4720. The f value of all three inhibitor molecules in the water medium is increased compared to the gaseous phase. The obtained electronic oscillator strength values indicate that the molecules can absorb light more efficiently.
The molecular electrostatic potential (MEP) map, which is based on electronic density and chemical reactivity sites, provides details on the total charge distribution of the molecule. The reactive areas of the molecule can therefore be easily determined using the MEP surface. The MEP map given in Figure 6 shows the colors red, orange, yellow, green, and blue, which are related to the distribution of the electronic densities. The difference in MEP in both the gaseous and solvent (water) phases will provide a quantitative perspective of the molecular characteristics under physiological conditions. In L1 and L2, the NO2− group attached to the phenyl ring is distributed over the red region in the gaseous phase and the dark-red region in the water medium. This indicates that, in the water medium, molecules become more nucleophilic and their interactions with the protein may be increased. In L2 and L3, the pyrrolidine ring becomes slightly electron deficient in a water medium. The carbonyl group of L3 in the water medium becomes more nucleophilic than in the gaseous phase. The oxygen of the carbonyl group has an interaction with the ASN279 of tau with a bond length of 2.5 Å, which may be even stronger in the water medium because of the increase in nucleophilicity. The comparison of band-gap energy values of the molecules in the gaseous and solvent phases is listed in Table S5 (available in the supplementary data). The band gap of all three molecules is slightly decreased in the water phase compared to the gaseous phase. Therefore, the solvent phase studies indicate the reactivity of molecules increased under a water medium that simulates physiological conditions. This means the probability of inhibitor–protein interactions is higher in physiological conditions.

3. Materials and Methods

The absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties of the compounds were predicted by using the SwissADME and ProTox-II tools [34,35,36,37]. Target prediction for the inhibitor molecules was recorded with the SwissTargetPrediction. Intestinal absorption and blood–brain barrier (BBB) penetration properties were determined by the BOILED-Egg method [38]. The geometries of the molecules considered in this work have been initially optimized using density functional theory (DFT) at the ωB97XD/6-311++G(2d, 2p) [39,40,41] level of theory. We used (E)-configuration structures in all the computational calculations throughout this work, which is the minimum-energy structure in the potential energy surface. This fact is well supported in the literature [16,42]. The time-dependent density-functional theory (TD-DFT), frontier molecular orbital (FMO), and molecular electrostatic potential (MEP) surface analyses are performed at the same level of theory. The protein-ligand docking process is carried out using AutoDock Vina [43,44]. The binding mechanism between the molecules and their corresponding proteins is examined using molecular dynamics simulation, which is carried out using NAMD/VMD [45,46]. The electronic transitions of the compounds are calculated by using the TD-DFT method [47,48]. All the geometry optimization calculations are performed with the Gaussian suite of programs [49].

4. Conclusions

In summary, we report the results of the in silico analysis performed on the potential dual Aβ and p-tau inhibitor molecules for AD treatment. We designed seven synthetically possible inhibitor molecules (L1–7) as dual Aβ and p-tau inhibitors. The binding affinity of these molecules towards both Aβ and p-tau is confirmed by docking analysis, which agrees well with the early pharmacodynamic assessments. All seven molecules exhibited BBB-crossing and gastrointestinal-absorption ability, which is essential for AD treatment. Out of these seven molecules, L1–3 were shortlisted for further investigation due to their higher docking score towards both proteins. LYS280 and LYS281 are the most interacting amino-acid residues in p-tau, and PHE8 is the most interactive amino-acid residue of Aβ. As per the MD simulation data, L3 possesses the highest similarity with the apo-proteins. Hence, L3 could be a lead molecule out of the seven investigated systems. L1 and L2 showed variance in the RMSD analysis from apo-Aβ. Therefore, L1 and L2 are potential molecules only for single-target inhibition of p-tau aggregation. Although the results point to L3 as the lead molecule, in-depth experimental studies are required to confirm the efficiency of all the inhibitor molecules (L1–7). Hence, these molecules can be considered for experimental studies, as they fulfilled the prerequisites of being a drug candidate when assessed theoretically. We put forward these molecules in the experimental analysis, since they can be prepared through different synthetic routes. The results show that the molecules have outstanding pharmacological characteristics with remarkable potential as AD-treating drugs. Further, this work will enlighten theoretical and experimental researchers involved in developing novel pharmacophores. This study can further act as a stepping stone for the auxiliary development of small organic multi-target inhibitor molecules in AD treatment.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28031388/s1, Figure S1. Optimized 3D-tube structures of seven dual inhibitors, Table S1—Toxicity data of all inhibitor molecules, Table S2—Type of interactions between inhibitor molecules and proteins in the complex, Table S3—TD-DFT calculations of the excited state of selected inhibitor molecules in the gaseous phase, Table S4—TD-DFT calculations of the excited state of selected inhibitor molecules in the solvent (water) phase, Table S5—The energy of HOMO and LUMO and their band-gap energy of potential inhibitor molecules in the gaseous and solvent phases, Table S6—Cartesian coordinates of the optimized geometries of inhibitor molecules in the gaseous phase (L1–7) at the ωB97XD/6-311++G(2d,2p) level.

Author Contributions

Conceptualization, N.J., V.S.T. and K.T.; investigation, data curation, visualization, writing and original draft, N.J.; manuscript review and editing, supervision, V.S.T. and K.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by VIT SEED grant SG20210209.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article and supplementary data also available.

Acknowledgments

The authors acknowledge VIT for providing the computational resources for this research work.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Haghighijoo, Z.; Zamani, L.; Moosavi, F.; Emami, S. Therapeutic Potential of Quinazoline Derivatives for Alzheimer’s Disease: A Comprehensive Review. Eur. J. Med. Chem. 2022, 227, 113949. [Google Scholar] [CrossRef] [PubMed]
  2. Kumar, A.; Kumar, A. Alzheimer’s Disease Therapy: Present and Future Molecules. Neuromethods 2018, 132, 3–22. [Google Scholar] [CrossRef]
  3. Breijyeh, Z.; Karaman, R. Comprehensive Review on Alzheimer’s Disease: Causes and Treatment. Molecules 2020, 25, 5789. [Google Scholar] [CrossRef]
  4. Abubakar, M.B.; Sanusi, K.O.; Ugusman, A.; Mohamed, W.; Kamal, H.; Ibrahim, N.H.; Khoo, C.S.; Kumar, J. Alzheimer’s Disease: An Update and Insights Into Pathophysiology. Front. Aging Neurosci. 2022, 14, 742408. [Google Scholar] [CrossRef]
  5. Rauk, A. The Chemistry of Alzheimer’s Disease. Chem. Soc. Rev. 2009, 38, 2698–2715. [Google Scholar] [CrossRef] [PubMed]
  6. Bai, R.; Guo, J.; Ye, X.-Y.; Xie, Y.; Xie, T. Oxidative Stress: The Core Pathogenesis and Mechanism of Alzheimer’s Disease. Ageing Res. Rev. 2022, 77, 101619. [Google Scholar] [CrossRef]
  7. dos Santos Picanco, L.C.; Ozela, P.F.; de Fatima de Brito Brito, M.; Pinheiro, A.A.; Padilha, E.C.; Braga, F.S.; de Paula da Silva, C.H.T.; dos Santos, C.B.R.; Rosa, J.M.C.; da Silva Hage-Melim, L.I. Alzheimer’s Disease: A Review from the Pathophysiology to Diagnosis, New Perspectives for Pharmacological Treatment. Curr. Med. Chem. 2016, 25, 3141–3159. [Google Scholar] [CrossRef]
  8. Wang, E.J.; Wu, M.Y.; Lu, J.H. Ferulic Acid in Animal Models of Alzheimer’s Disease: A Systematic Review of Preclinical Studies. Cells 2021, 10, 2653. [Google Scholar] [CrossRef]
  9. Hansson, O.; Zetterberg, H.; Vanmechelen, E.; Vanderstichele, H.; Andreasson, U.; Londos, E.; Wallin, A.; Minthon, L.; Blennow, K. Evaluation of Plasma Aβ40 and Aβ42 as Predictors of Conversion to Alzheimer’s Disease in Patients with Mild Cognitive Impairment. Neurobiol. Aging 2010, 31, 357–367. [Google Scholar] [CrossRef]
  10. Pradeepkiran, J.A.; Reddy, P.H. Phosphorylated Tau Targeted Small-Molecule PROTACs for the Treatment of Alzheimer’s Disease and Tauopathies. Biochim. Biophys. Acta-Mol. Basis Dis. 2021, 1867, 166162. [Google Scholar] [CrossRef]
  11. Malafaia, D.; Albuquerque, H.M.T.; Silva, A.M.S. Amyloid-β and Tau Aggregation Dual-Inhibitors: A Synthetic and Structure-Activity Relationship Focused Review. Eur. J. Med. Chem. 2021, 214, 113209. [Google Scholar] [CrossRef]
  12. Brunden, K.R.; Trojanowski, J.Q.; Lee, V.M.Y. Advances in Tau-Focused Drug Discovery for Alzheimer’s Disease and Related Tauopathies. Nat. Rev. Drug Discov. 2009, 8, 783–793. [Google Scholar] [CrossRef] [PubMed]
  13. Mazanetz, M.P.; Fischer, P.M. Untangling Tau Hyperphosphorylation in Drug Design for Neurodegenerative Diseases. Nat. Rev. Drug Discov. 2007, 6, 464–479. [Google Scholar] [CrossRef] [PubMed]
  14. Alonso, A.D.; Di Clerico, J.; Li, B.; Corbo, C.P.; Alaniz, M.E.; Grundke-Iqbal, I.; Iqbal, K. Phosphorylation of Tau at Thr212, Thr231, and Ser262 Combined Causes Neurodegeneration. J. Biol. Chem. 2010, 285, 30851–30860. [Google Scholar] [CrossRef] [PubMed]
  15. Hosokawa, T.; Saito, T.; Asada, A.; Fukunaga, K.; Hisanaga, S.I. Quantitative Measurement of in Vivo Phosphorylation States of Cdk5 Activator P35 by Phos-Tag SDS-PAGE. Mol. Cell. Proteom. 2010, 9, 1133–1143. [Google Scholar] [CrossRef]
  16. Şahin, S.; Dege, N. A Newly Synthesized Small Molecule: The Evaluation against Alzheimer’s Disease by In Silico Drug Design and Computational Structure Analysis Methods. J. Mol. Struct. 2021, 1236, 130337. [Google Scholar] [CrossRef]
  17. Veber, D.F.; Johnson, S.R.; Cheng, H.-Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular Properties That Influence the Oral Bioavailability of Drug Candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef]
  18. 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. 1997, 23, 3–25. [Google Scholar] [CrossRef]
  19. Xu, Y.; He, Z.; Liu, H.; Chen, Y.; Gao, Y.; Zhang, S.; Wang, M.; Lu, X.; Wang, C.; Zhao, Z.; et al. 3D-QSAR, Molecular Docking, and Molecular Dynamics Simulation Study of Thieno[3,2- b ]Pyrrole-5-Carboxamide Derivatives as LSD1 Inhibitors. RSC Adv. 2020, 10, 6927–6943. [Google Scholar] [CrossRef]
  20. Ho, T.H.; Charlet-B, N.; Poulos, M.G.; Singh, G.; Swanson, M.S.; Cooper, T.A. Muscleblind Proteins Regulate Alternative Splicing. EMBO J. 2004, 23, 3103–3112. [Google Scholar] [CrossRef] [Green Version]
  21. Banerjee, P.; Eckert, A.O.; Schrey, A.K.; Preissner, R. ProTox-II: A Webserver for the Prediction of Toxicity of Chemicals. Nucleic Acids Res. 2018, 46, W257–W263. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, W.; Jiang, Y.; Dostert, K.-H.; O’Brien, C.P.; Riedel, W.; Savara, A.; Schauermann, S.; Tkatchenko, A. Catalysis beyond Frontier Molecular Orbitals: Selectivity in Partial Hydrogenation of Multi-Unsaturated Hydrocarbons on Metal Catalysts. Sci. Adv. 2017, 3, e1700939. [Google Scholar] [CrossRef] [PubMed]
  23. Pangajavalli, S.; Ranjithkumar, R.; Ramaswamy, S. Structural, Hirshfeld, Spectroscopic, Quantum Chemical and Molecular Docking Studies on 6b′, 7′, 8′, 9′-Tetrahydro-2H,6′H-Spiro[Acenaphthylene-1,11′-Chromeno [3,4-a]Pyrrolizine]-2,6′(6a′H,11a′H)-Dione. J. Mol. Struct. 2020, 1209, 127921. [Google Scholar] [CrossRef]
  24. Lotfy, K. Molecular Modeling, Docking and ADMET of Dimethylthiohydantoin Derivatives for Prostate Cancer Treatment. J. Biophys. Chem. 2015, 06, 91–117. [Google Scholar] [CrossRef]
  25. Sahu, V.; Sharma, P.; Kumar, A. Impact of Global and Local Reactivity Descriptors on the Hetero-Diels-Alder Reaction of Enaminothione with Various Electrophiles. J. Chil. Chem. Soc. 2014, 59, 2327–2334. [Google Scholar] [CrossRef]
  26. Namitha, R.; Selvi, G.; Mageswari, D. Theoretical and Experimental Investigations on the Biological Studies on Novel Triazino Quinolines. Int. J. Comput. Intell. Res. 2017, 13, 1021–1025. [Google Scholar]
  27. Fichou, Y.; Al-Hilaly, Y.K.; Devred, F.; Smet-Nocca, C.; Tsvetkov, P.O.; Verelst, J.; Winderickx, J.; Geukens, N.; Vanmechelen, E.; Perrotin, A.; et al. The Elusive Tau Molecular Structures: Can We Translate the Recent Breakthroughs into New Targets for Intervention? Acta Neuropathol. Commun. 2019, 7, 31. [Google Scholar] [CrossRef]
  28. Seidler, P.M.; Boyer, D.R.; Rodriguez, J.A.; Sawaya, M.R.; Cascio, D.; Murray, K.; Gonen, T.; Eisenberg, D.S. Structure-Based Inhibitors of Tau Aggregation. Nat. Chem. 2018, 10, 170–176. [Google Scholar] [CrossRef]
  29. Berman, H.M. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef]
  30. Schrödinger, L.D.W. The PyMOL Molecular Graphics System; Version 2.0; Schrödinger: New York, NY, USA, 2020. [Google Scholar]
  31. Dassault Systèmes BIOVIA. Discovery Studio Modeling Environment; Dassault Systèmes: San Diego, CA, USA, 2017. [Google Scholar]
  32. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef]
  33. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A Web-Based Graphical User Interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef] [PubMed]
  34. Brooks, B.R.; Brooks, C.L.; Mackerell, A.D.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The Biomolecular Simulation Program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef] [PubMed]
  35. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef]
  36. Banerjee, P.; Dehnbostel, F.O.; Preissner, R. Prediction Is a Balancing Act: Importance of Sampling Methods to Balance Sensitivity and Specificity of Predictive Models Based on Imbalanced Chemical Data Sets. Front. Chem. 2018, 6, 362. [Google Scholar] [CrossRef]
  37. Drwal, M.N.; Banerjee, P.; Dunkel, M.; Wettig, M.R.; Preissner, R. ProTox: A Web Server for the In Silico Prediction of Rodent Oral Toxicity. Nucleic Acids Res. 2014, 42, W53–W58. [Google Scholar] [CrossRef]
  38. Daina, A.; Zoete, V. A BOILED-Egg To Predict Gastrointestinal Absorption and Brain Penetration of Small Molecules. ChemMedChem 2016, 11, 1117–1121. [Google Scholar] [CrossRef]
  39. Chai, J.-D.; Head-Gordon, M. Long-Range Corrected Hybrid Density Functionals with Damped Atom–Atom Dispersion Corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615. [Google Scholar] [CrossRef] [PubMed]
  40. Clark, T.; Chandrasekhar, J.; Spitznagel, G.W.; Schleyer, P.V.R. Efficient Diffuse Function-Augmented Basis Sets for Anion Calculations. III. The 3-21+G Basis Set for First-Row Elements, Li-F. J. Comput. Chem. 1983, 4, 294–301. [Google Scholar] [CrossRef]
  41. Krishnan, R.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent Molecular Orbital Methods. XX. A Basis Set for Correlated Wave Functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar] [CrossRef]
  42. Gaenko, A.V.; Devarajan, A.; Gagliardi, L.; Lindh, R.; Orlandi, G. Ab Initio DFT Study of Z–E Isomerization Pathways of N–Benzylideneaniline. Theor. Chem. Acc. 2007, 118, 271–279. [Google Scholar] [CrossRef]
  43. Trott, O.; Olson, A.J. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 2009, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. 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]
  45. Phillips, J.C.; Hardy, D.J.; Maia, J.D.C.; Stone, J.E.; Ribeiro, J.V.; Bernardi, R.C.; Buch, R.; Fiorin, G.; Hénin, J.; Jiang, W.; et al. Scalable Molecular Dynamics on CPU and GPU Architectures with NAMD. J. Chem. Phys. 2020, 153, 044130. [Google Scholar] [CrossRef] [PubMed]
  46. Humphrey, W.; Dalke, A.; Schulten, K. Visual Molecular Dynamics. J. Mol. Graph. 1996, 14.1, 33–38. [Google Scholar] [CrossRef] [PubMed]
  47. Adamo, C.; Jacquemin, D. The Calculations of Excited-State Properties with Time-Dependent Density Functional Theory. Chem. Soc. Rev. 2013, 42, 845–856. [Google Scholar] [CrossRef]
  48. Laurent, A.D.; Adamo, C.; Jacquemin, D. Dye Chemistry with Time-Dependent Density Functional Theory. Phys. Chem. Chem. Phys. 2014, 16, 14334–14356. [Google Scholar] [CrossRef]
  49. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision A.02; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
Molecules 28 01388 g001 550
Figure 1. Molecular structures of seven dual inhibitors (L1–7).
Figure 1. Molecular structures of seven dual inhibitors (L1–7).
Molecules 28 01388 g001
Molecules 28 01388 g002 550
Figure 2. BOILED-Egg model of L1–7 to predict gastrointestinal absorption and BBB penetration. (HIA—human gastrointestinal absorption; PGP+—substrate of P-glycoprotein; PGP−—non-substrate of P-glycoprotein).
Figure 2. BOILED-Egg model of L1–7 to predict gastrointestinal absorption and BBB penetration. (HIA—human gastrointestinal absorption; PGP+—substrate of P-glycoprotein; PGP−—non-substrate of P-glycoprotein).
Molecules 28 01388 g002
Molecules 28 01388 g003 550
Figure 3. Docking poses of L1–3 in the binding site of Aβ and p-tau.
Figure 3. Docking poses of L1–3 in the binding site of Aβ and p-tau.
Molecules 28 01388 g003
Molecules 28 01388 g004 550
Figure 4. Two-dimensional representation of interactions between L1–3 and proteins.
Figure 4. Two-dimensional representation of interactions between L1–3 and proteins.
Molecules 28 01388 g004
Molecules 28 01388 g005 550
Figure 5. Root-mean-square deviation (RMSD) plots of Aβ and p-tau in both apo (in yellow) and with inhibitor-bound forms (L1, L2, and L3 are in red, green, and blue respectively) in the molecular simulation of 1 ns.
Figure 5. Root-mean-square deviation (RMSD) plots of Aβ and p-tau in both apo (in yellow) and with inhibitor-bound forms (L1, L2, and L3 are in red, green, and blue respectively) in the molecular simulation of 1 ns.
Molecules 28 01388 g005
Molecules 28 01388 g006 550
Figure 6. Molecular electrostatic potential map of L1–3.
Figure 6. Molecular electrostatic potential map of L1–3.
Molecules 28 01388 g006
Table
Table 1. In silico drug-likeness properties of L1–7 obtained using the SwissADME tool.
Table 1. In silico drug-likeness properties of L1–7 obtained using the SwissADME tool.
MoleculesMWTPSALog PESOL
Log S		Log Kp (cm s−1)ViolationsSynthetic AccessibilityTarget
LipinskiGhoseMAPTAmyloid Beta
L1359.4261.424.15−5.56−4.75003.05yesyes
L2345.3961.423.88−5.27−4.92002.99yesyes
L3323.3578.492.91−4.11−5.76002.89yesyes
L4337.3778.493.11−4.41−5.59002.96yesyes
L5306.32762.51−4.02−5.92002.83yesyes
L6309.3661.423.15−4.42−5.33002.89yesyes
L7305.3363.113.11−4.42−5.46002.76yesyes
MW—molecular weight; TPSA—topological polar surface area; log P—lipophilicity of arithmetic mean of the values predicted by the five proposed methods; ESOL log S—water solubility by estimated solubility method; log Kp—permeability coefficient; MAPT—microtubule-associated protein tau.
Table
Table 2. The computed global reactivity descriptors for L1–7.
Table 2. The computed global reactivity descriptors for L1–7.
Parameters/DescriptorsL1L2L3L4L5L6L7
EHOMO (eV)−7.73−7.62−8.43−8.12−8.30−7.77−8.06
ELUMO (eV)−0.74−0.49−0.86−0.62−1.31−0.43−1.13
Energy band gap (eV)
(ΔE = ELUMO − EHOMO)		6.997.137.577.506.997.346.93
Ionization potential
(I = −EHOMO) (eV)		7.737.628.438.128.307.778.06
Electron affinity (eV)
(A = −ELUMO)		0.730.480.850.621.300.421.12
Chemical hardness
(η = (I − A)/2) 		3.493.563.783.753.493.673.46
Chemical softness
(σ = 1/2 η)		0.140.140.130.130.140.130.14
Electronegativity
(χ = (I + A)/2)		4.234.054.644.374.804.104.59
Chemical potential
(μ = −(I + A)/2)		−4.23−4.05−4.64−4.37−4.80−4.10−4.59
Electrophilicity index
(ω = μ2/2η)		2.562.302.842.553.302.293.04
Maximum charge transfer index
(ΔNmax = −μ/η)		1.211.131.221.161.371.111.32
Table
Table 3. The molecular docking scores of L1–7 with Aβ and p-tau.
Table 3. The molecular docking scores of L1–7 with Aβ and p-tau.
MoleculesDocking Score (kcal/mol)
Amyloid Betap-tau
L1−6.0−5.6
L2−5.8−5.1
L3−5.3−5.4
L4−5.2−4.6
L5−5.0 −5.0
L6−5.0 −5.0
L7−4.9−5.2
Reference molecule [16]-−5.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Job, N.; Thimmakondu, V.S.; Thirumoorthy, K. In Silico Drug Design and Analysis of Dual Amyloid-Beta and Tau Protein-Aggregation Inhibitors for Alzheimer’s Disease Treatment. Molecules 2023, 28, 1388. https://doi.org/10.3390/molecules28031388

AMA Style

Job N, Thimmakondu VS, Thirumoorthy K. In Silico Drug Design and Analysis of Dual Amyloid-Beta and Tau Protein-Aggregation Inhibitors for Alzheimer’s Disease Treatment. Molecules. 2023; 28(3):1388. https://doi.org/10.3390/molecules28031388

Chicago/Turabian Style

Job, Nisha, Venkatesan S. Thimmakondu, and Krishnan Thirumoorthy. 2023. "In Silico Drug Design and Analysis of Dual Amyloid-Beta and Tau Protein-Aggregation Inhibitors for Alzheimer’s Disease Treatment" Molecules 28, no. 3: 1388. https://doi.org/10.3390/molecules28031388

Article Metrics

Back to TopTop