Discovery of novel 5α-reductase type II inhibitors by pharmacophore modelling,virtual screening,molecular docking and molecular dynamics simulations

Benign prostatic hyperplasia (BPH) is caused by augmented levels of androgen dihydrotestosterone (DHT) which is involved in the growth of the prostate in humans. 5α-Reductase type II (5αR2) is an intracellular enzyme that catalyses the formation of DHT from testosterone; hence, the inhibition of 5αR2 has emerged as one of the most promising strategies for the treatment of BPH. In this study, a computational approach that integrates ligand-based pharmacophore modelling, virtual screening, molecular docking and molecular dynamics (MD) simulations was adopted to discover novel 5αR2 inhibitors with less side effects. After validating by Fischer's randomisation and Güner–Henry test, the best quantitative pharmacophore model (Hypo1), consisting of two hydrogen-bond acceptors and three hydrophobic features, was subsequently used as a three-dimensional-query in virtual screening to identify potential hits from Maybridge and National Cancer Institute databases. These hits were further filtered by ADMET (absorption, distribution, metabolism, elimination and toxicology) and molecular docking experiments, and their binding stabilities were validated by 10-ns MD simulations. Finally, only one hit was identified as a potential lead based on higher predicted inhibitory activity to 5αR2 compared with the most active inhibitor (finasteride). Our results further suggest that this potential lead could easily be synthesised and has structural novelty, making it a promising candidate for treating BPH.     

E-pharmacophore-based virtual screening to identify GSK-3β inhibitors

Glycogen synthase kinase-3β (GSK-3β) is a serine/threonine kinase which has attracted significant attention during recent years in drug design studies. The deregulation of GSK-3β increased the loss of hippocampal neurons by triggering apoptosis-mediating production of neurofibrillary tangles and alleviates memory deficits in Alzheimer’s disease (AD). Given its role in the formation of neurofibrillary tangles leading to AD, it has been a major therapeutic target for intervention in AD, hence was targeted in the present study. Twenty crystal structures were refined to generate pharmacophore models based on energy involvement in binding co-crystal ligands. Four common e-pharmacophore models were optimized from the 20 pharmacophore models. Shape-based screening of four e-pharmacophore models against nine established small molecule databases using Phase v3.9 had resulted in 1800 compounds having similar pharmacophore features. Rigid receptor docking (RRD) was performed for 1800 compounds and 20 co-crystal ligands with GSK-3β to generate dock complexes. Interactions of the best scoring lead obtained through RRD were further studied with quantum polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area. Comparing the obtained leads to 20 co-crystal ligands resulted in 18 leads among them, lead1 had the lowest docking score, lower binding free energy and better binding orientation toward GSK-3β. The 50?ns MD simulations run confirmed the stable nature of GSK-3β-lead1 docking complex. The results from RRD, QPLD, IFD and MD simulations confirmed that lead1 might be used as a potent antagonist for GSK-3β.     

Optimization of novel reversible Bruton’s tyrosine kinase inhibitors identified using Tethering-fragment-based screens

Since the approval of ibrutinib for the treatment of B-cell malignancies in 2012, numerous clinical trials have been reported using covalent inhibitors to target Bruton’s tyrosine kinase (BTK) for oncology indications. However, a formidable challenge for the pharmaceutical industry has been the identification of reversible, selective, potent molecules for inhibition of BTK. Herein, we report application of Tethering-fragment-based screens to identify low molecular weight fragments which were further optimized to improve on-target potency and ADME properties leading to the discovery of reversible, selective, potent BTK inhibitors suitable for pre-clinical proof-of-concept studies.     

Role of Bruton’s tyrosine kinase in macrophage apoptosis

Macrophages and polymorphonuclear cells (PMNs) rapidly respond to microbial and immune inflammatory stimuli and die during these responses. We have shown earlier that many macrophage and PMN functions are compromised in x-linked immunodeficient (Xid) mice with functional deficiency in Bruton’s tyrosine kinase (Btk). We now report that Btk-deficient macrophages show enhanced susceptibility to apoptotic death on exposure to the microbial and immune inflammatory signals bacterial lipopolysaccharide (LPS) and interferon-gamma (IFNγ) in vitro. In vivo in mixed bone marrow (BM) chimeras Btk deficiency leads primarily to loss of peripheral macrophage numbers without affecting BM development, suggesting a role of inflammation-induced apoptosis in regulating macrophage life span. Surprisingly, Btk deficiency does not affect macrophage apoptosis induced by DNA damage or CD95 engagement. Reactive nitrogen and oxygen species also do not contribute to inflammation-induced apoptosis, but apoptotic process involves loss of mitochondrial potential, shows increased activation of caspase 9 and enhanced loss of Bcl-xL. The lack of pro-survival signaling through the Btk-phosphotidylinositol 3-kinase-Akt pathway, and persistent MEK signaling, lead to enhanced death in Btk-deficient macrophages only downstream of inflammatory triggers. These data underline the complex role of Btk in the regulation of macrophage survival and function.     

Identification of abelson tyrosine kinase inhibitors as potential therapeutics for Alzheimer’s disease using multiple e-pharmacophore modeling and molecular dynamics

Efforts to combat Alzheimer’s disease are focused predominantly on inhibiting the activity of the enzyme(s) that have been identified to be responsible for the production of the amyloid-forming peptide. However, the inherent complexity associated with the network of pathways leading to the disease may involve additional targets for designing effective therapies. Recent experimental findings have identified abelson tyrosine kinase, a non-receptor kinase as a new target for Alzheimer’s. In this work, we employed energy optimized multiple pharmacophore modeling strategy from multiple c-Abl structures bound with ligands in the inactive ATP binding conformation (DFG-out). Virtual screening followed by docking of molecules from ChemBridge resulted in the identification of 10 best scoring molecules. MD simulations of the top three complexes revealed that Compound A, C are the most stable complexes with the most persistent protein–ligand interactions consistent with the calculated binding affinities for the top three compounds. Given the implied role of c-Abl not only in AD but in Parkinson’s disease, the identified compounds may serve as leads for effective neurotherapeutics.     

Purine derivatives as potent Bruton’s tyrosine kinase (BTK) inhibitors for autoimmune diseases

Investigation of various heterocyclic core isosteres of imidazopyrazines 1 & 2 yielded purine derivatives 3 & 8 as potent and selective BTK inhibitors. Subsequent SAR studies of the purine series led to the discovery of 20 as a leading compound. Compound 20 is very selective when screened against a panel of 400 kinases and is a potent inhibitor in cellular assays of human B cell function including B-Cell proliferation and CD86 cell surface expression and exhibited in vivo efficacy in a mouse PCA model. Its X-ray co-crystal structure with BTK shows that the high selectivity is gained from filling a BTK specific lipophilic pocket. However, physical and ADME properties leading to low oral exposure hindered further development.     

Identification of potential PKC inhibitors through pharmacophore designing, 3D-QSAR and molecular dynamics simulations targeting Alzheimer’s disease

Protein kinases are ubiquitously expressed as Serine/Threonine kinases, and play a crucial role in cellular activities. Protein kinases have evolved through stringent regulation mechanisms. Protein kinases are also involved in tauopathy, thus are important targets for developing Anti-Alzheimer’s disease compounds. Structures with an indole scaffold turned out to be potent new leads. With the aim of developing new inhibitors for human protein kinase C, here we report the generation of four point 3D geometric featured pharmacophore model. In order to identify novel and potent PKCθ inhibitors, the pharmacophore model was screened against 80,000,00 compounds from various chemical databases such as., ZINC, SPEC, ASINEX, which resulted in 127 compound hits, and were taken for molecular docking filters (HTVS, XP docking). After in-depth analysis of binding patterns, induced fit docking (flexible) was employed for six compounds along with the cocrystallized inhibitor. Molecular docking study reveals that compound 6F found to be tight binder at the active site of PKCθ as compared to the cocrystal and has occupancy of 90 percentile. MM-GBSA also confirmed the potency of the compound 6F as better than cocrystal. Molecular dynamics results suggest that compound 6F showed good binding stability of active sites residues similar to cocrystal 7G compound. Present study corroborates the pharmacophore-based virtual screening, and finds the compound 6F as a potent Inhibitor of PKC, having therapeutic potential for Alzheimer’s disease. Worldwide, 46.8 million people are believed to be living with Alzheimer’s disease. When elderly population increases rapidly and neurodegenerative burden also increases in parallel, we project the findings from this study will be useful for drug developing efforts targeting Alzheimer’s disease.     

Molecular dynamics simulation and molecular docking studies of 1,4-Dihydropyridines as P-glycoprotein’s allosteric inhibitors

P-glycoprotein (P-gp) is a main factor contributing to multidrug resistance. The effect of this transporter protein on limiting the effectiveness of chemotherapy has been shown by various studies. In a previous report, we synthesized some 14-dihydropyridine (DHP) derivatives as inhibitors of human P-gp. In the present study, a computational approach has been exploited to reveal the main interactions between DHPs and P-gp. In order to do this, homology modeling was performed to obtain a model of the protein. Then, molecular dynamics simulation was used to refine the constructed model of P-gp in the presence of the lipids bilayer. Model validation was performed with several tools. Finally, molecular docking followed by MD simulation of ligand–protein complex was employed to elucidate the binding mode and the dynamical changes of protein with/without DHPs bound. The results emphasized that interaction of the residues Gln912, Ser909, Arg905, Ser474, Val472 with DHPs play a crucial role in the inhibitory of these ligands and this was in a relatively good accordance with the results reported in the experimental studies.     

Design,synthesis and biological evaluation of novel 3-substituted pyrazolopyrimidine derivatives as potent Bruton’s tyrosine kinase (BTK) inhibitors

A series of 3-substituted pyrazolopyrimidine derivatives as BTK inhibitors were designed by structure-based drug design and they were synthesized, evaluated by enzyme-based assay and anti-proliferation against Ramos and Raji cells. Most of them displayed good inhibitory activities against both BTK and B-cell lymphoblastic leukemia lines in vitro. Among them, compound 8a exhibited excellent potency (IC₅₀?=?7.95?nM against BTK enzyme, 8.91?μM against Ramos cells and 1.80?μM against Raji cells), with a better hydrophilicity (ClogP?=?3.33). These explorations provided new clues to discover 3-substituted pyrazolopyrimidine derivatives as novel anti-tumor agents.     

Docking,molecular dynamics,binding energy-MM-PBSA studies of naphthofuran derivatives to identify potential dual inhibitors against BACE-1 and GSK-3β

BACE-1 and GSK-3β both are potential therapeutic drug targets for Alzheimer’s disease. Recently, both these targets received attention for designing dual inhibitors. Till now only two scaffolds (triazinone and curcumin) derivatives have been reported as BACE-1 and GSK-3β dual inhibitors. In our previous work, we have reported first in class dual inhibitor for BACE-1 and GSK-3β. In this study, we have explored other naphthofuran derivatives for their potential to inhibit BACE-1 and GSK-3β through docking, molecular dynamics, binding energy (MM-PBSA). These computational methods were performed to estimate the binding affinity of naphthofuran derivatives towards the BACE-1 and GSK-3β. In the docking results, two derivatives (NS7 and NS9) showed better binding affinity as compared to previously reported inhibitors. Hydrogen bond occupancy of NS7 and NS9 generated from MD trajectories showed good interaction with the flap residues Gln73, Thr72 of BACE-1 and Arg141, Thr138 residues of GSK-3β. MM-PBSA and energy decomposition per residue revealed different components of binding energy and relative importance of amino acid involved in binding. The results showed that the binding of inhibitors was majorly governed by the hydrophobic interactions and suggesting that hydrophobic interactions might be the key to design dual inhibitors for BACE1-1 and GSK-3β. Distance between important pair of amino acid residues indicated that BACE-1 and GSK-3β adopt closed conformation and become inactive after ligand binding. The results suggested that naphthofuran derivatives might act as dual inhibitor against BACE-1 and GSK-3β.     

Pharmacophore-based virtual screening of catechol-o-methyltransferase (COMT) inhibitors to combat Alzheimer’s disease

Alzheimer’s disease (AD) is one of the most significant neurodegenerative disorders and its symptoms mostly appear in aged people. Catechol-o-methyltransferase (COMT) is one of the known target enzymes responsible for AD. With the use of 23 known inhibitors of COMT, a query has been generated and validated by screening against the database of 1500 decoys to obtain the GH score and enrichment value. The crucial features of the known inhibitors were evaluated by the online ZINC Pharmer to identify new leads from a ZINC database. Five hundred hits were retrieved from ZINC Pharmer and by ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering by using FAF-Drug-3 and 36 molecules were considered for molecular docking. From the COMT inhibitors, opicapone, fenoldopam, and quercetin were selected, while ZINC63625100_413 ZINC39411941_412, ZINC63234426_254, ZINC63637968_451, and ZINC64019452_303 were chosen for the molecular dynamics simulation analysis having high binding affinity and structural recognition. This study identified the potential COMT inhibitors through pharmacophore-based inhibitor screening leading to a more complete understanding of molecular-level interactions.     

3D-QSAR,molecular docking and molecular dynamics studies of a series of RORγt inhibitors

The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R²_cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R²_pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity.     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Conversion of carbazole carboxamide based reversible inhibitors of Bruton’s tyrosine kinase (BTK) into potent,selective irreversible inhibitors in the carbazole,tetrahydrocarbazole, and a new 2,3-dimethylindole series

Incorporation of a suitably-placed electrophilic group transformed a series of reversible BTK inhibitors based on carbazole-1-carboxamide and tetrahydrocarbazole-1-carboxamide into potent, irreversible inhibitors. Removal of one ring from the core of these compounds provided a potent irreversible series of 2,3-dimethylindole-7-carboxamides having excellent potency and improved selectivity, with the additional advantages of reduced lipophilicity and molecular weight.     

Benzofuran–appended 4-aminoquinazoline hybrids as epidermal growth factor receptor tyrosine kinase inhibitors: synthesis,biological evaluation and molecular docking studies

A series of 2-arylbenzo[b]furan–appended 4-aminoquinazoline hybrids were prepared and evaluated for cytotoxicity in vitro against the human lung cancer (A549), colorectal adenocarcinoma (Caco-2), hepatocellular carcinoma (C3A) and cervical cancer (HeLa) cell lines. Compounds 10d and 10j exhibited significant cytotoxicity against the C3A and Caco-2 cell lines and induced apoptosis in these cell lines. Likewise, compounds 10d and 10e exhibited significant inhibitory activity towards epidermal growth factor receptor-tyrosine kinase phosphorylation (IC₅₀ values of 29.3?nM and 31.1?nM, respectively) against Gefitinib (IC₅₀?=?33.1?nM). Molecular docking of compounds 10 into EGFR-TK active site suggests that they bind to the region of EGFR like Gefitinib does.

     

Mitigation of reactive metabolite formation for a series of 3-amino-2-pyridone inhibitors of Bruton’s tyrosine kinase (BTK)

Reactive metabolites have been putatively linked to many adverse drug reactions including idiosyncratic toxicities for a number of drugs with black box warnings or withdrawn from the market. Therefore, it is desirable to minimize the risk of reactive metabolite formation for lead molecules in optimization, in particular for non-life threatening chronic disease, to maximize benefit to risk ratio. This article describes our effort in addressing reactive metabolite issues for a series of 3-amino-2-pyridone inhibitors of BTK, e.g. compound 1 has a value of 459 pmol/mg protein in the microsomal covalent binding assay. Parallel approaches were taken to successfully resolve the issues: establishment of a predictive screening assay with correlation association of covalent binding assay, identification of the origin of reactive metabolite formation using MS/MS analysis of HLM as well as isolation and characterization of GSH adducts. This ultimately led to the discovery of compound 7 (RN941) with significantly reduced covalent binding of 26 pmol/mg protein.