Identification of human tau-tubulin kinase 1 inhibitors: an integrated e-pharmacophore-based virtual screening and molecular dynamics simulation

Computational Design of Novel Tau-Tubulin Kinase 1 Inhibitors for Neurodegenerative Diseases

The tau-tubulin kinase 1 (TTBK1) protein is a casein kinase 1 superfamily member located at chromosome 6p21.1. It is expressed explicitly in the brain, particularly in the cytoplasm of cortical and hippocampal neurons. TTBK1 has been implicated in the phosphorylation and aggregation of tau in Alzheimer's disease (AD). Considering its significance in AD, TTBK1 has emerged as a promising target for AD treatment. In the present study, we identified novel TTBK1 inhibitors using various computational techniques. We performed a virtual screening-based docking study followed by E-pharmacophore modeling, cavity-based pharmacophore, and ligand design techniques and found ZINC000095101333, LD7, LD55, and LD75 to be potential novel TTBK1 lead inhibitors. The docking results were complemented by Molecular Mechanics/Generalized Born Surface Area (MMGBSA) calculations. The molecular dynamics (MD) simulation studies at a 500 ns scale were carried out to monitor the behavior of the protein toward the identified ligands. Pharmacological and ADME/T studies were carried out to check the drug-likeness of the compounds. In summary, we identified a new series of compounds that could effectively bind the TTBK1 receptor. The newly designed compounds are promising candidates for developing therapeutics targeting TTBK1 for AD.

Antisense oligonucleotide-based targeting of Tau-tubulin kinase 1 prevents hippocampal accumulation of phosphorylated tau in PS19 tauopathy mice

Tau tubulin kinase-1 (TTBK1), a neuron-specific tau kinase, is highly expressed in the entorhinal cortex and hippocampal regions, where early tau pathology evolves in Alzheimer's disease (AD). The protein expression level of TTBK1 is elevated in the cortex brain tissues with AD patients compared to the control subjects. We therefore hypothesized that antisense oligonucleotide (ASO) based targeting Ttbk1 could prevent the accumulation of phosphorylated tau, thereby delaying the development of tau pathology in AD. Here we show that in vivo administration of ASO targeting mouse Ttbk1 (ASO-Ttbk1) specifically suppressed the expression of Ttbk1 without affecting Ttbk2 expression in the temporal cortex of PS19 tau transgenic mice. Central administration of ASO-Ttbk1 in PS19 mice significantly reduced the expression level of representative phosphor-tau epitopes relevant to AD at 8 weeks post-dose, including pT231, pT181, and pS396 in the sarkosyl soluble and insoluble fractions isolated from hippocampal tissues as determined by ELISA and pS422 in soluble fractions as determined by western blotting. Immunofluorescence demonstrated that ASO-Ttbk1 significantly reduced pS422 phosphorylated tau intensity in mossy fibers region of the dentate gyrus in PS19 mice. RNA-sequence analysis of the temporal cortex tissue revealed significant enrichment of interferon-gamma and complement pathways and increased expression of antigen presenting molecules (Cd86, Cd74, and H2-Aa) in PS19 mice treated with ASO-Ttbk1, suggesting its potential effect on microglial phenotype although neurotoxic effect was absent. These data suggest that TTBK1 is an attractive therapeutic target to suppress TTBK1 without compromising TTBK2 expression and pathological tau phosphorylation in the early stages of AD.

In silico design and identification of new peptides for mitigating hIAPP aggregation in type 2 diabetes

The aberrant misfolding and self-aggregation of human islet amyloid polypeptide (hIAPP or amylin) into cytotoxic aggregates are implicated in the pathogenesis of type 2 diabetes (T2D). Among various inhibitors, short peptides derived from the amyloidogenic regions of hIAPP have been employed as hIAPP aggregation inhibitors due to their low immunogenicity, biocompatibility, and high chemical diversity. Recently, hIAPP fragment HSSNN18-22 was identified as an amyloidogenic sequence and displayed higher antiproliferative activity to RIN-5F cells. Various hIAPP aggregation inhibitors have been designed by chemical modifications of the highly amyloidogenic sequence (NFGAIL) of hIAPP. In this work, a library of pentapeptides based on fragment HSSNN18-22 was designed and assessed for their efficacy in blocking hIAPP aggregation using an integrated computational screening approach. The binding free energy calculations by molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method identified HSSQN and HSSNQ that bind to hIAPP monomer with a binding affinity of -21.25 ± 4.90 and -19.73 ± 3.10 kcal/mol, respectively, which is notably higher as compared to HSSNN (-11.90 ± 4.12 kcal/mol). The sampling of the non aggregation-prone helical conformation was notably increased from 23.5 ± 3.0 in the hIAPP monomer to 38.1 ± 3.6, and 33.8 ± 3.0% on the incorporation of HSSQN, and HSSNQ, respectively, which indicate reduced aggregation propensity of hIAPP monomer. The pentapeptides, HSSQN and HSSNQ, identified as hIAPP aggregation inhibitors in this work can be further conjugated with various metal chelating peptides to yield more efficacious and clinically relevant multifunctional modulators for targeting various pathological hallmarks of T2D.

CK2 and protein kinases of the CK1 superfamily as targets for neurodegenerative disorders

Casein kinases are involved in a variety of signaling pathways, and also in inflammation, cancer, and neurological diseases. Therefore, they are regarded as potential therapeutic targets for drug design. Recent studies have highlighted the importance of the casein kinase 1 superfamily as well as protein kinase CK2 in the development of several neurodegenerative pathologies, such as Alzheimer's disease, Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis. CK1 kinases and their closely related tau tubulin kinases as well as CK2 are found to be overexpressed in the mammalian brain. Numerous substrates have been detected which play crucial roles in neuronal and synaptic network functions and activities. The development of new substances for the treatment of these pathologies is in high demand. The impact of these kinases in the progress of neurodegenerative disorders, their bona fide substrates, and numerous natural and synthetic compounds which are able to inhibit CK1, TTBK, and CK2 are discussed in this review.

The structural, functional, and dynamic effect of Tau tubulin kinase1 upon a mutation: A neuro-degenerative hotspot

Alzheimer's disease (AD) is a progressive disorder that causes brain cells to degenerate and die. AD is one of the common causes of dementia that leads to a decline in thinking, behavioral and social skills that disrupts a person's ability to function independently. Tau-tubulin kinase1 (TTBK1) is a crucial disease regulating AD protein, which is majorly responsible for the phosphorylation and accumulation of tau protein at specific Serine/Threonine residues found in paired helical filaments, suggesting its role in tauopathy. TTBK1 involvement in many diseases and the restricted expression of TTBK1 to the central nervous system (CNS) makes TTBK1 an attractive therapeutic target for tauopathies. The genetic variations in TTBK1 are primarily involved in the TTBK1 pathogenesis. This study highlighted the destabilizing, damaging and deleterious effect of the mutation R142Q on TTBK1 structure through computational predictions and molecular dynamics simulations. The protein deviation, fluctuations, conformational dynamics, solvent accessibility, hydrogen bonding, and the residue-residue mapping confirmed the mutant effect to cause structural aberrations, suggesting overall destabilization due to the protein mutation. The presence of well-defined free energy minima was observed in TTBK1-wild type, as opposed to that in the R142Q mutant, reflecting structural deterioration. The overall findings from the study reveal that the presence of R142Q mutation on TTBK1 is responsible for the structural instability, leading to disruption of its biological functions. The mutation could be used as future diagnostic markers in treating AD.

Pharmacophore- based virtual screening, 3D- QSAR, molecular docking approach for identification of potential dipeptidyl peptidase IV inhibitors

Pharmacophore modeling, molecular docking, and in silico ADME studies have been carried out to determine the binding mode and drug likeliness profile of Pyrrolidine derivatives as Dipeptidyl peptidase IV inhibitors. A five point pharmacophore model (AAADH_1) was generated using 96 compounds with IC₅₀ values ranging from 1.8 to 13000 nM. A statistically significant 3D-QSAR model was generated from the pharmacophore hypothesis. The model had a high correlation coefficient (R² - 0.92), cross validation coefficient (Q² - 0.776) and F value (F - 144) at 6 component Partial least square factor. Pearson r of 0.7124 indicated greater degree of confidence on the model. The accuracy and predictability of the generated model were checked by Enrichment factor, Receiver operating characteristic curves, area under curve, Boltzmann-enhanced discrimination of Receiver operating characteristic and the Robust initial enhancement. To identify novel and potent Dipeptidyl peptidase IV inhibitors, virtual screening was performed using the ligand and database screening. Considering the potent hit molecules on the basis of pharmacophore virtual screening, we have designed new molecules and further subjected to see the interaction with protein. The catalytic domain of Dipeptidyl peptidase IV enzyme in complex with Vildagliptin (PDB Code: 6B1E) was obtained from protein data bank with resolution 1.77 A°. Compound 75 showed the better binding (dock score -7.966) with protein than standard drug vildagliptin (Dock Score -6.554). The hits obtained on virtual screening of the database have provided new chemical starting points for design and development of novel Dipeptidyl peptidase IV inhibitory agents.Communicated by Ramaswamy H. Sarma.

In silico design of peptide inhibitors of tubulin: amyloid-β as a lead compound

Microtubule is one of the most studied targets in cancer research. Stabilizing and destabilizing of the microtubule by targeting its building block tubulin are common mechanisms of microtubule targeting agents. Cancer associates inversely with Alzheimer's disease (AD). So the rate of developing AD is significantly slower in patients with cancer and vice versa. Amyloid-β (Aβ) peptide inhibits tubulin polymerization and induces apoptotic death of cancer cells. We studied the interactions of Aβ with tubulin using protein-protein docking and MD simulation. Aβ bond to the vicinity of the vinblastine binding site and interacted with the H6-H7 loop. Interaction of Aβ with H6-H7 loop blocked nucleotide exchange and may be attributed as a possible reason for blocking of tubulin polymerization. We designed new Aβ-based peptidic inhibitors of tubulin using visual inspection and alanine scanning method. P1 (FRHYHHFFELV) and P9 (HYHHF) bound efficiently to tubulin and also interacted with the H6-H7 loop. Obtained results indicated that proposed peptides could potentially inhibit nucleotide exchange as Aβ.Communicated by Ramaswamy H. Sarma.

Identification of highly active natural thyroid hormone receptor agonists by pharmacophore-based virtual screening

Thyroid hormone receptor (TR) is an important target for the treatment of metabolic diseases. The X-ray crystallographic data for the TR complexed with different ligands were employed to generate feature-based pharmacophore models of the active site of TR receptor. The derived hypothesis was then used to find novel hit compounds through an in silico virtual screening on the Universal Natural Products Database (UNPD). The binding mode and action mechanism of the hit compounds were further investigated by molecular docking and molecular dynamics studies, then compounds possessing similar binding site with the crystal ligand were subjected to binding activity assay. Finally, UNPD19665 and UNPD184785 were proved to be more active than crystal TR ligands, with a binding activity value of 9.82 nM and 12.62 nM, respectively. Communicated by Ramaswamy H. Sarma.