Journal of Advanced Pharmaceutical Technology & Research

ORIGINAL ARTICLE
Year
: 2020  |  Volume : 11  |  Issue : 4  |  Page : 194--201

In silico docking and comparative ADMET profile of different glycogen synthase kinase 3 beta inhibitors as the potential leads for the development of anti-Alzheimer drug therapy


Manisha Prajapat1, Phulen Sarma1, Nishant Shekhar1, Hardeep Kaur1, Sanjay Singh2, Subodh Kumar1, Harpinder Kaur1, Saniya Mahendiratta1, Amit Raj Sharma3, Sukhmandeep Kaur1, Vidya M Mahalmani1, Bikash Medhi1,  
1 Department of Pharmacology, PGIMER, Chandigarh, India
2 Department of Biotechnology, Thapar Institute of Engineering and Technology, Patiala, Punjab, India
3 Department of Neurology, PGIMER, Chandigarh, India

Correspondence Address:
Prof. Bikash Medhi
Department of Pharmacology, Postgraduate Institute of Medical Education and Research, Chandigarh
India

Abstract

Glycogen synthase kinase 3 beta (GSK3 β) plays a key role in pathologic hyper phosphorylation of tau and plays an important role in the pathogenesis of Alzheimer's disease. In the present study, we have screened a set of potential hits in in silico platform to gain insight regarding binding profile with the target (GSK3 β) from molecular docking, ADME/T, and molecular dynamics (MD) simulations. The three screened compounds 6-BIBEO, 6-BIO, and SB216763 topped the docking score chart when subjected to hard scoring function extraprecision of GLIDE. The active site dynamics study through MD simulations provides insights on residues Asp133, Val135, and Ile62 which are in a state of minimum deviation from their mean special position while they interact with the respective ligands. The same molecules also displayed favorable pharmacokinetic profile, negative Ames test and falls correctly within drug-likeliness rules. These agents can be taken forward further for the development of anti-Alzheimer's drug therapy.



How to cite this article:
Prajapat M, Sarma P, Shekhar N, Kaur H, Singh S, Kumar S, Kaur H, Mahendiratta S, Sharma AR, Kaur S, Mahalmani VM, Medhi B. In silico docking and comparative ADMET profile of different glycogen synthase kinase 3 beta inhibitors as the potential leads for the development of anti-Alzheimer drug therapy.J Adv Pharm Technol Res 2020;11:194-201


How to cite this URL:
Prajapat M, Sarma P, Shekhar N, Kaur H, Singh S, Kumar S, Kaur H, Mahendiratta S, Sharma AR, Kaur S, Mahalmani VM, Medhi B. In silico docking and comparative ADMET profile of different glycogen synthase kinase 3 beta inhibitors as the potential leads for the development of anti-Alzheimer drug therapy. J Adv Pharm Technol Res [serial online] 2020 [cited 2020 Dec 4 ];11:194-201
Available from: https://www.japtr.org/text.asp?2020/11/4/194/297689


Full Text



 Introduction



Alzheimer's disease (AD), the most prevalent form of dementia in the aging world population usually in the age group of 65 years or above.[1] Glycogen synthase kinase 3 beta (GSK3 β) mediated hyperphosphorylation of tau plays a major role in the pathogenesis of AD.[2] GSK3 β is a serine/threonine kinase.[3] Phosphorylation on tyr216 residue generates the active conformation of GSK3 β.[4] GSK3 β favors phosphorylation of prephosphorylated substrate. The primed phosphorylation of residues Ser235 and Ser404 of the tau-protein by other kinases such as CDK5, subsequently aids phosphorylation by GSK-3 on residues Thr231 and Ser400 on tau protein.[5] In this study, we have evaluated 10 different GSK3 β inhibitors (NSC69386, 6BIO, TCG24, Bio-acetoxime, CHIR98014, 6-BIBEO, 6-BIDECO, 6-BIMYEO, LY2090314, SB216763, and SB415286) in in silico platform for the development of potential leads for the treatment of AD.

 Methodology



Retrieval of target structure: The structure of GSK-3 beta, which is target receptor protein of human, was retrieved from protein data bank server (1UV5).[6] The detail of aminoacid sequence is shown in [Figure 1].{Figure 1}

Ligands: We have evaluated 10 different ligands (6-BIO, bio-Acetoxime, TC-G24, CHIR98014, NSC693868, 6-BIBEO, 6-BIDECO, 6-BIMYEO, SB216763, SB415286 and LY2090314). The structure of ligand molecules was retrieved from PubChem[7] [Figure 2]. Since, LY2090314[8] is already in preclinical phase,[9] it was taken as controls.[10],[11]{Figure 2}

Molecular docking: All in silico evaluations were carried out using Schrödinger Maestro suite 2019.

Pharmacokinetics properties of ligands

Admet SAR,[12] AMDET labs[13] and Swiss ADME server[14] were used for the evaluation of pharmacokinetic properties.

Drug-likeliness

Ligands were evaluated for drug likeness using Lipinski, Ghose, Veber, Egan, Muegge criteria using Swiss ADME software.

Molecular dynamics simulations

The overall three best performers and the control were further evaluated in molecular dynamics (MD) studies using a three step process of “system building,” minimization, and MD simulation using the Desmond module of Schrodinger Inc (simulation time-50 ns, ensemble class nonproliferation treaty, temperature-300K and pressure = 1 bar).

 Results and Discussion



Chemical structure

Ligand serial number, name, and their chemical structure are illustrated in [Figure 2].

Docking profile of the ligands

According to docking score 7, compound shown good binding profile in docking (6BIBEO > 6BIO > CHIR98014 > SB415286>NSC693868 > LY2090314 > SB216763 > TC-G24). Docking score data are showed in [Table 1] and [Figure 3], [Figure 4], [Figure 5], [Figure 6].{Table 1}{Figure 3}{Figure 4}{Figure 5}{Figure 6}

ADMET profile

Comparative ADMET profiles of the different agents are showed in [Table 2].{Table 2}

Physicochemical properties

Log P value (distribution coefficient P)

In our study, all the compounds had Log P value below 5. The compound NSC693868 had poor lipid bi-layer permeability (Log P = 0.558) compared to other ligands (6BIO, 6BEBIO, BIOACETOXIME, 6-BIDECO, CHIR98014, TC-G24 and SB216763).

LogD7.4

logD7.4 value was in the low in case of NSC693868 and SB415286, highlighting their hydrophilicity. Apart from these two ligands, LogD7.4 value was between 2 and 3.2. None of the compounds had LogD7.4 value higher than 3.5 [Table 2].

Absorption

All the ligands were found to be positive for human intestinal absorption. However, only 3 of the ligands (NSC693868, TC-G24, and SB216763) were permeable through Caco-2. None of the ligands were substrate of P-gp; however, most of the drugs were P-gp inhibitor except NSC693868, CHIR98014, and SB415286. None of the ligands inhibited renal organic cation transporter [Table 2].

Distribution

These compounds were distributed in three main sub-cellular regions that are plasma membrane, lysosome, and mitochondria. 6-BIO and BIOACETOXIME showed distribution in plasma membrane, NSC693868 in lysosome and other remaining drugs in mitochondria. The plasma protein binding were <80% in case of NSC6938, 80%–90% in case of 6-BIO, 6-BIMYEO, CHIR98014 and LY2090314 and was more than 90% in case of 6-BIBEO, Bio-acetoxime, 6-BIDECO, TCG-24, SB216763, and SB415286. All of the ligands were permeable through blood–brain barrier except 6 BIDECO and SB415286 [Table 2].

Excretion

The half-life of all the ligands was in between 0.8 and 1.9 h. The highest T1/2 life is 1.932 h was in case of SB216763, and highest clearance rate was 1.87 ml/min/kg which was found with NSC693868 [Table 2].

Toxicity

All the ligands that were evaluated were weak hERG channel inhibitors, except CHIR98014, which showed strong inhibition. Five ligands showed toxicity in Ames test (BIOACETOXIME, NSC693868, CHIR98014, TC-G24, and SB415286), whereas the rest of the compounds did not show toxicity. None of the ligands were potential carcinogens. All the ligands showed high fish toxicity and TP toxicity and low honey bee toxicity. Skin sensitization was seen in case of CHIR98014, TCG24. Predicted rat LD50 of the compounds ranged from 2.4 to 2.6 mol/kg [Table 2].

Metabolism profile

All of the ligands were CYP3A4 substrate except compound NSC693868 and CHIR98014, which were non substrates. None of the ligands were substrate of CYP2C9 and CYP2D6. Only six ligands (6-BIO, 6-BIBEO, BIOACETOXIME, 6-BIDECO, SB216763, and LY2090314) were inhibitors of CYP3A4, whereas rest were noninhibitors [Table 3].{Table 3}

Drug-likeliness of the ligands

All the selected ligands followed these rules expect ligand CHIR98014 (violated all 5 rules) and LY2090314 (violated Lipinski rule and Ghose rule) [Table 4].{Table 4}

Selection of ligand for further molecular dynamics simulations

Upon comparing and integrating the knowledge acquired from docking and ADMET score, 3 ligands (6-BIO, 6-BIBEO, and SB216763) were selected out to be good hit and were taken forward for MD simulations. LY2090314 was taken as control.

Molecular dynamic simulation

Root mean square deviation (RMSD): The respective RMSD observed were 1.8Š for 6-BIO, 1.9Š for 6-BIBEO, 1.75Š for SB216763, and 2.0 for LY2090314. These RMSD values were well-within the acceptable range of 0–3Š and also RMSD progression equilibrates when it approaches the end of trajectory hence implying to a stable protein-ligand complex formation which can be inferred for positive interaction of 1UV5 with all three ligands (6BIO, 6BIBEO, SB216763, and LY2090314) [Figure 7]b, [Figure 7]e, [Figure 7]h and [Figure 7]k. Details of RMSF values are showed in [Figure 7]a, [Figure 7]d, [Figure 7]g and [Figure 7]J.{Figure 7}

Protein-Ligand Interactions: Residues Asp133 and Val135 were found to be predominantly important residues exhibiting the high percentage of H-bonding with the three selected candidate (6-BIO, 6-BIBEO, and SB216763) compounds. Residue Asp133, Val135, and Ile62 are the major residues involved in the core-binding cavity showing predominant interaction with the ligand. Moreover, compounds 6-BIBEO (Ile62, Val70, Ala83, and Leu188), 6-BIO (Val70, Ala83, and Leu188) SB216763 (Ile62, Val70, Ala83, Val110, Leu132, and Leu188), and LY2090314 (Ile62, Phe67, Val70, Lys85, and Leu188) exhibit some degree of hydrophobic interactions. Moreover, water bridges were observed with core active site residues in 6-BIBEO (Arg141 and Ile62), in 6-BIO (Gln185 and Tyr140) and in SB216763 (Ile62, Tyr134, Pro136, Tyr138, Val135, Thr138, Arg141, Tyr140, and Gln185) and LY2090314 (Lys85, Asp133, Val135, Pro136, Thr138, Arg141, Lys183, Gln185, and Asp200) [Figures 7c, [Figure 7]f, [Figure 7]i, [Figure 7]l and [Figure 8].{Figure 8}

 Discussion



GSK3 β has a key role in hyperphosphorylation of tau and plays an important role in regulation of intra-neuronal hyperphosphorylated tau level.[15] In this study, we have targeted GSK3 β for in silico identification of possible hits for the development of anti-Alzheimer's therapy. Among the 10 selected ligands, 6-BIBEO, 6-BIO, CHIR98014, SB415286, NSC693868, LY2090314, SB216763, and TC-G24 showed good binding profile, evaluated in terms of docking score. logD7.4 value was in the lower side in case of NSC693868 and SB415286, highlighting their hydrophilicity. AD being a disease of the central nervous system, BBB permeability is a major factor. In our study, all the ligands were able to pass the BBB except two ligands (6-BIDECO and SB415286).

CHIR98014 was found to be a strong inhibitor to the hERG, whereas hERG inhibitory profiles of other drugs were comparable to control. Ames test was positive for BIOACETOXIME, NSC693868, CHIR98014, TCG24, and SB415286.

Taking in account, the pharmacokinetics result and target binding profile, 3 of the compounds (6-BIBEO > 6-BIO > SB216763) were found to be the most suitable agents for further development process.

After performing MD simulation of these selected compound with the common target GSK3 β, we can interpret from the data providing insight on RMSD of GSK3 β-ligand complex was 1.8Š for 6-BIO, 1.9Š for 6-BIBEO, 1.75Š for SB216763, and 2.0 for LY2090314 with respect to its C-alpha position. The trajectory frames are well under the scale of 3Š and stabilized as it propagated further. The root mean square fluctuation (RMSF) and ligand-contact ratio of GSK3-Bresidues showing with the respective ligands emphasize that SB216763 shows greater percentage interaction with the core residues of target site, with higher number of water bridge formation, while 6-BIO and 6-BIBEO were dynamically similar in behavior and also formed similar interaction profile with the target site, of which H-bond was a major part. The positive control LY2090314 showed uneven interaction pattern with the GSK3 β residues but showed great affinity by the means of water bridge formation [Figure 7]l. Ligand RMSF details shown in [Figure 8], on the basis of which 6-BIBEO emerges out to be most stable at all trajectory frames.

 Conclusion



Among the ten ligands evaluated, 6-BIO, 6-BIBEO, and SB216763 needs further evaluation as probable anti-Alzheimer's drugs considering the in silico ADME parameters, toxicity, blood brain barrier permeability, docking scores, and MD simulation.

Acknowledgment

The authors acknowledge Prof. Prasad V. Bharatam and Mr. Gurudutt Dubey of NIPER, Mohali and Mr. Vinod Devaraji and Dr. Prajwal Nadekar from Schrodinger, Bangalore.

fFinancial support and sponsorship

Nil.

Conflicts of interest

There are no conflicts of interest.

References

1Weller J, Budson A. Current understanding of Alzheimer's disease diagnosis and treatment. F1000Res 2018 Jul 31;7:1161.
2Chong FP, Ng KY, Koh RY, Chye SM. Tau Proteins and Tauopathies in Alzheimer's Disease. Cell Mol Neurobiol 2018;38:965-80.
3Jantrapirom S, Nimlamool W, Chattipakorn N, Chattipakorn S, Temviriyanukul P, Inthachat W, et al. Liraglutide Suppresses Tau Hyperphosphorylation, Amyloid Beta Accumulation through Regulating Neuronal Insulin Signaling and BACE-1 Activity. International Journal of Molecular Sciences. 2020;21:1725.
4Das TK, Jana P, Chakrabarti SK, Abdul Hamid MRW. Curcumin Downregulates GSK3 and Cdk5 in Scopolamine-Induced Alzheimer's Disease Rats Abrogating Aβ40/42 and Tau Hyperphosphorylation. J Alzheimers Dis Rep 3:257-67.
5Hanger DP, Noble W. Functional Implications of Glycogen Synthase Kinase-3-Mediated Tau Phosphorylation. International Journal of Alzheimer's Disease. 2011;2011:1-11.
6Meijer L, Skaltsounis A-L, Magiatis P, Polychronopoulous P, Knockaert M, Leost M, et al. Gsk-3-Selective Inhibitors Derived from Tyrian Purple Indurubins. ChemBiol. 2003;10:1255.
7PubChem. PubChem [Internet]. Available from: https://pubchem.ncbi.nlm.nih.gov/. [Last cited on 2020 Sep 30].
8PubChem. 3-[6-Fluoro-10-(piperidine-1-carbonyl)-1,10-diazatricyclo[6.4.1.04,13]trideca-2,4,6,8(13)-tetraen-3-yl]-4-imidazo[1,2-a]pyridin-3-ylpyrrole-2,5-dione [Internet]. Available from: https://pubchem.ncbi.nlm.nih.gov/compound/10029385. [Last cited on 2020 Sep 30].
9Vasdev N, Garcia A, Stableford WT, Young AB, Meyer JH, Houle S, et al. Synthesis and ex vivo evaluation of carbon-11 labelled N-(4-methoxybenzyl)-N′-(5-nitro-1,3-thiazol-2-yl)urea ([11C]AR-A014418): A radiolabelled glycogen synthase kinase-3β specific inhibitor for PET studies. Bioorganic & Medicinal Chemistry Letters. 2005;15:5270-3.
10Pandey MK, DeGrado TR. Glycogen Synthase Kinase-3 (GSK-3)-Targeted Therapy and Imaging. Theranostics 2016;6:571-93.
11Zamek-Gliszczynski MJ, Abraham TL, Alberts JJ, Kulanthaivel P, Jackson KA, Chow KH, et al. Pharmacokinetics, Metabolism, and Excretion of the Glycogen Synthase Kinase-3 Inhibitor LY2090314 in Rats, Dogs, and Humans: A Case Study in Rapid Clearance by Extensive Metabolism with Low Circulating Metabolite Exposure. Drug Metab Dispos 2013;41:714-26.
12admetSAR @ LMMD [Internet]. Available from: http://lmmd.ecust.edu.cn/admetsar1. [Last cited on 2019 May 30].
13Home-ADMElab: ADMET Prediction,ADMET Predictor,QSAR,ADMET Database [Internet]. Available from: http://admet.scbdd.com/ [Last cited on 2020 Sep 30].
14SwissADME [Internet]. Available from: http://www.swissadme.ch/. [Last cited on 2020 Sep 30].
15Llorens-Marítin M, Jurado J, Hernández F, Ã vila J. GSK-3Î2, a pivotal kinase in Alzheimer disease. Front Mol Neurosci [Internet]. 2014;7. Available from: http://journal.frontiersin.org/article/10.3389/fnmol.2014.00046/abstract. [Last cited on 2020 Sep 30].