Synthesis of 2-,4,-6-, and/or 7-substituted quinoline derivatives as human dihydroorotate dehydrogenase (hDHODH) inhibitors and anticancer agents: 3D QSAR-assisted design

Following our research for human dihydroorotate dehydrogenase (hDHODH) inhibitors as anticancer agents, herein we describe 3D QSAR-based design, synthesis and in vitro screening of 2-,4,-6-, and/or 7-substituted quinoline derivatives as hDHODH inhibitors and anticancer agents. We have designed 2-,4,-6-, and/or 7-substituted quinoline derivatives and predicted their hDHODH inhibitory activity based on 3D QSAR study on 45 substituted quinoline derivatives as hDHODH inhibitors, and also predicted toxicity. Designed compounds were docked into the binding site of hDHODH. Designed compounds which showed good predictive activity, no toxicity, and good docking score were selected for the synthesis, and in vitro screening as hDHODH inhibitors in an enzyme inhibition assay, and anticancer agents in MTT assay against cancer cell lines (HT-29 and MDA-MB-231). Synthesized compounds 7 and 14 demonstrated IC₅₀ value of 1.56?µM and 1.22?µM, against hDHODH, respectively, and these are our lead compounds for the development of new hDHODH inhibitors and anticancer agents.     

Triazoloacridin-6-ones as novel inhibitors of the quinone oxidoreductases NQO1 and NQO2

A range of triazoloacridin-6-ones functionalized at C5 and C8 have been synthesized and evaluated for ability to inhibit NQO1 and NQO2. The compounds were computationally docked into the active site of NQO1 and NQO2, and calculated binding affinities were compared with IC₅₀ values for enzyme inhibition. Excellent correlation coefficients were demonstrated suggesting a predictive QSAR model for this series of structurally similar analogues. From this we have identified some of these triazoloacridin-6-ones to be the most potent NQO2 inhibitors so far reported.     

Inhibitory evaluation of Curculigo latifolia on α-glucosidase,DPP (IV) and in vitro studies in antidiabetic with molecular docking relevance to type 2 diabetes mellitus

The inhibition of α-glucosidase and DPP enzymes capable of effectively reducing blood glucose level in the management of type 2 diabetes. The purpose of the present study is to evaluate the inhibitory potential of α-glucosidase and DPP (IV) activity including with the 2-NBDG uptake assay and insulin secretion activities through in vitro studies. The selected of active compounds obtained from the screening of compounds by LC-MS were docked with the targeted enzyme that involved in the mechanism of T2DM. From the results, root extracts displayed a better promising outcome in α-glucosidase (IC₅₀ 2.72 ± 0.32) as compared with the fruit extracts (IC₅₀ 3.87 ± 0.32). Besides, root extracts also displayed a better activity in the inhibition of DPP (IV), enhance insulin secretion and glucose uptake activity. Molecular docking results revealing that phlorizin binds strongly with α-glucosidase, DPP (IV) and Insulin receptor (IR) enzymes with achieving the lowest binding energy value. The present work suggests several of the compounds have the potential that contribute towards inhibiting α-glucosidase and DPP (IV) and thus effective in lowering post-prandial hyperglycaemia.     

Design,synthesis, docking,ADMET studies,and anticancer evaluation of new 3-methylquinoxaline derivatives as VEGFR-2 inhibitors and apoptosis inducers

Vascular endothelial growth factor receptor-2 (VEGFR-2) plays a critical role in cancer angiogenesis. Inhibition of VEGFR-2 activity proved effective suppression of tumour propagation. Accordingly, two series of new 3-methylquinoxaline derivatives have been designed and synthesised as VEGFR-2 inhibitors. The synthesised derivatives were evaluated in vitro for their cytotoxic activities against MCF-7and HepG2 cell lines. In addition, the VEGFR-2 inhibitory activities of the target compounds were estimated to indicate the potential mechanism of their cytotoxicity. To a great extent, the results of VEGFR-2 inhibition were highly correlated with that of cytotoxicity. Compound 27a was the most potent VEGFR-2 inhibitor with IC₅₀ of 3.2 nM very close to positive control sorafenib (IC₅₀ = 3.12 nM). Such compound exhibited a strong cytotoxic effect against MCF-7 and HepG2, respectively with IC₅₀ of 7.7 and 4.5 µM in comparison to sorafenib (IC₅₀ = 3.51 and 2.17 µM). In addition, compounds 28, 30f, 30i, and 31b exhibited excellent VEGFR-2 inhibition activities (IC₅₀ range from 4.2 to 6.1 nM) with promising cytotoxic activity. Cell cycle progression and apoptosis induction were investigated for the most active member 27a. Also, the effect of 27a on the level of caspase-3, caspase-9, and BAX/Bcl-2 ratio was determined. Molecular docking studies were implemented to interpret the binding mode of the target compounds with the VEGFR-2 pocket. Furthermore, toxicity and ADMET calculations were performed for the synthesised compounds to study their pharmacokinetic profiles     

Molecular docking studies on DMDP derivatives as human DHFR inhibitors

Molecular docking is routinely used for understanding drug‐receptor interaction in modern drug design. Here, we describe the docking of 2, 4-diamino-5-methyl-5-deazapteridine (DMDP) derivatives as inhibitors to human dihydrofolate reductase (DHFR). We docked 78 DMDP derivates collected from literature to DHFR and studied their specific interactions with DHFR. A new shape-based method, LigandFit, was used for docking DMDP derivatives into DHFR active sites. The result indicates that the molecular docking approach is reliable and produces a good correlation coefficient (r² = 0.499) for the 73 compounds between docking score and IC₅₀ values (Inhibitory Activity). The chloro substituted naphthyl ring of compound 63 makes significant hydrophobic contact with Leu 22, Phe 31 and Pro 61 of the DHFR active site leading to enhanced inhibition of the enzyme. The docked complexes provide better insights to design more potent DHFR inhibitors prior to their synthesis.     

Exploration of Novel Inhibitors for Class I Histone Deacetylase Isoforms by QSAR Modeling and Molecular Dynamics Simulation Assays

Histone deacetylases (HDAC) are metal-dependent enzymes and considered as important targets for cell functioning. Particularly, higher expression of class I HDACs is common in the onset of multiple malignancies which results in deregulation of many target genes involved in cell growth, differentiation and survival. Although substantial attempts have been made to control the irregular functioning of HDACs by employing various inhibitors with high sensitivity towards transformed cells, limited success has been achieved in epigenetic cancer therapy. Here in this study, we used ligand-based pharmacophore and 2-dimensional quantitative structure activity relationship (QSAR) modeling approaches for targeting class I HDAC isoforms. Pharmacophore models were generated by taking into account the known IC₅₀ values and experimental energy scores with extensive validations. The QSAR model having an external R² value of 0.93 was employed for virtual screening of compound libraries. 10 potential lead compounds (C1-C10) were short-listed having strong binding affinities for HDACs, out of which 2 compounds (C8 and C9) were able to interact with all members of class I HDACs. The potential binding modes of HDAC2 and HDAC8 to C8 were explored through molecular dynamics simulations. Overall, bioactivity and ligand efficiency (binding energy/non-hydrogen atoms) profiles suggested that proposed hits may be more effective inhibitors for cancer therapy.     

Molecular design,synthesis and in vitro biological evaluation of thienopyrimidine–hydroxamic acids as chimeric kinase HDAC inhibitors: a challenging approach to combat cancer

A series of thieno[2,3-d]pyrimidine-based hydroxamic acid hybrids was designed and synthesised as multitarget anti-cancer agents, through incorporating the pharmacophore of EGFR, VEGFR2 into the inhibitory functionality of HDAC6. Three compounds (12c, 15b and 20b) were promising hits, whereas (12c) exhibited potent VEGFR2 inhibition (IC₅₀=185 nM), potent EGFR inhibition (IC₅₀=1.14 µM), and mild HDAC6 inhibition (23% inhibition). Moreover, compound (15c) was the most potent dual inhibitor among all the synthesised compounds, as it exhibited potent EGFR and VEGFR2 inhibition (IC₅₀=19 nM) and (IC₅₀=5.58 µM), respectively. While compounds (20d) and (7c) displayed nanomolar selective kinase inhibition with EGFR IC₅₀= 68 nM and VEGFR2 IC₅₀= 191 nM, respectively. All of the synthesised compounds were screened in vitro for their cytotoxic effect on 60 human NCI tumour cell lines. Additionally, molecular docking studies and ADMET studies were carried out to gain further insight into their binding mode and predict the pharmacokinetic properties of all the synthesised inhibitors.     

Molecular modeling studies,synthesis and biological evaluation of dabigatran analogues as thrombin inhibitors

In this work, 48 thrombin inhibitors based on the structural scaffold of dabigatran were analyzed using a combination of molecular modeling techniques. We generated three-dimensional quantitative structure–activity relationship (3D-QSAR) models based on three alignments for both comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) to highlight the structural requirements for thrombin protein inhibition. In addition to the 3D-QSAR study, Topomer CoMFA model also was established with a higher leave-one-out cross-validation q² and a non-cross-validation r², which suggest that the three models have good predictive ability. The results indicated that the steric, hydrophobic and electrostatic fields play key roles in QSAR model. Furthermore, we employed molecular docking and re-docking simulation explored the binding relationship of the ligand and the receptor protein in detail. Molecular docking simulations identified several key interactions that were also indicated through 3D-QSAR analysis. On the basis of the obtained results, two compounds were designed and predicted by three models, the biological evaluation in vitro (IC₅₀) demonstrated that these molecular models were effective for the development of novel potent thrombin inhibitors.     

Isomeric Mono-, Di-, and Tri-Bromobenzo-1H-Triazoles as Inhibitors of Human Protein Kinase CK2α

To further clarify the role of the individual bromine atoms of 4,5,6,7-tetrabromotriazole (TBBt), a relatively selective inhibitor of protein kinase CK2, we have examined the inhibition (IC₅₀) of human CK2α by the two mono-, the four di-, and the two tri- bromobenzotriazoles relative to that of TBBt. Halogenation of the central vicinal C(5)/C(6) atoms proved to be a key factor in enhancing inhibitory activity, in that 5,6-di-Br₂Bt and 4,5,6-Br₃Bt were almost as effective inhibitors as TBBt, notwithstanding their marked differences in pK_a for dissociation of the triazole proton. The decrease in pK_a on halogenation of the peripheral C(4)/C(7) atoms virtually nullifies the gain due to hydrophobic interactions, and does not lead to a decrease in IC₅₀. Molecular modeling of structures of complexes of the ligands with the enzyme, as well as QSAR analysis, pointed to a balance of hydrophobic and electrostatic interactions as a discriminator of inhibitory activity. The role of halogen bonding remains debatable, as originally noted for the crystal structure of TBBt with CK2α (pdb1j91). Finally we direct attention to the promising applicability of our series of well-defined halogenated benzotriazoles to studies on inhibition of kinases other than CK2.     

Viper venom hyaluronidase and its potential inhibitor analysis: a multipronged computational investigation

Viper venom hyaluronidase (VV-HYA) inhibitors have long been used as therapeutic agents for arresting the local and systemic effects caused during its envenomation. Henceforth, to understand its structural features and also to identify the best potential inhibitor against it the present computational study was undertaken. Structure-based homology modeling of VV-HYA followed by its docking and free energy-based ranking analysis of ligand, the MD simulations of the lead complex was also performed. The sequence analysis and homology modeling of VV-HYA revealed a distorted (β/α)₈ folding as in the case of hydrolases family of proteins. Molecular docking of the resultant 3D structure of VV-HYA with known inhibitors (compounds 1–25) revealed the importance of molecular recognition of hotspot residues (Tyr 75, Arg 288, and Trp 321) other than that of the active site residues. It also revealed that Trp 321 of VV-HYA is highly important for mediating π–π interactions with ligands. In addition, the molecular docking and comparative free energy binding analysis was investigated for the VV-HYA inhibitors (compounds 1–25). Both molecular docking and relative free energy binding analysis clearly confirmed the identification of sodium chromoglycate (compound 1) as the best potential inhibitor against VV-HYA. Molecular dynamics simulations additionally confirmed the stability of their binding interactions. Further, the information obtained from this work is believed to serve as an impetus for future rational designing of new novel VV-HYA inhibitors with improved activity and selectivity.     

Synthesis,in vitro urease inhibition and molecular docking studies of some novel quinazolin-4(3H)-one derivatives containing triazole,thiadiazole and thiosemicarbazide functionalities

A new series of quinazolinone derivatives containing triazole, thiadiazole, thiosemicarbazide functionalities was synthesized and then screened for their in vitro urease inhibition properties. Most of the compounds showed excellent activity with IC₅₀ values ranging between 1.88 ± 0.17 and 6.42 ± 0.23 µg/mL, compared to that of thiourea (IC₅₀ = 15.06 ± 0.68) and acetohydroxamic acid (IC₅₀ = 21.03 ± 0.94), as reference inhibitors. Among the synthesized molecules, compounds 5c, 5e and 5a showed the best inhibitory effect against urease enzyme with IC₅₀ values of 1.88 ± 0.17 µg/mL, 1.90 ± 0.10 and 1.96 ± 0.07 µg/mL, respectively. Moreover in order to give better understanding of the inhibitory activity of synthesized compounds, molecular docking studies were applied at the target sites of jack bean urease enzyme (JBU). Their binding poses and energy calculations were analyzed using induced fit docking (IFD) and prime-MMGBSA tool. Binding poses of studied compounds were determined using induced fit docking (IFD) algorithms.     

Modification,Biological Evaluation and 3D QSAR Studies of Novel 2-(1,3-Diaryl- 4,5-Dihydro-1H-Pyrazol-5-yl)Phenol Derivatives as Inhibitors of B-Raf Kinase

A series of novel 2-(1,3-diaryl- 4,5-dihydro-1H-pyrazol-5-yl)phenol derivatives (C1–C24) have been synthesized. The B-Raf inhibitory activity and anti-proliferation activity of these compounds have been tested. Compound C6 displayed the most potent biological activity against B-Raf^V600E (IC₅₀ = 0.15 µM) and WM266.4 human melanoma cell line (GI₅₀ = 1.75 µM), being comparable with the positive control (Vemurafenib and Erlotinib) and more potent than our previous best compounds. The docking simulation was performed to analyze the probable binding models and poses while the QSAR model was built to check the previous work as well as to introduce new directions. This work aimed at seeking more potent inhibitors as well as discussing some previous findings. As a result, the introduction of ortho-hydroxyl group on 4,5-dihydro-1H-pyrazole skeleton did reinforce the anti-tumor activity while enlarging the group on N-1 of pyrazoline was also helpful.     

Synthesis and discovery of highly functionalized mono- and bis-spiro-pyrrolidines as potent cholinesterase enzyme inhibitors

Novel mono and bis spiropyrrolidine derivatives were synthesized via an efficient ionic liquid mediated, 1,3-dipolar cycloaddition methodology and evaluated in vitro for their AChE and BChE inhibitory activities in search for potent cholinesterase enzyme inhibitors. Most of the synthesized compounds displayed remarkable AChE inhibitory activities with IC₅₀ values ranging from 1.68 to 21.85 μM, wherein compounds 8d and 8j were found to be most active inhibitors against AChE and BChE with IC₅₀ values of 1.68 and 2.75 μM, respectively. Molecular modeling simulation on Torpedo californica AChE and human BChE receptors, showed good correlation between IC₅₀ values and binding interaction template of the most active inhibitors docked into the active site of their relevant enzymes.     

In silico approaches to identify novel myeloid cell leukemia-1 (Mcl-1) inhibitors for treatment of cancer

Myeloid cell leukemia-1 (Mcl-1) has been a validated and attractive target for cancer therapy. Over-expression of Mcl-1 in many cancers allows cancer cells to evade apoptosis and contributes to the resistance to current chemotherapeutics. Here, we identified new Mcl-1 inhibitors using a multi-step virtual screening approach. First, based on two different ligand-receptor complexes, 20 pharmacophore models were established by simultaneously using ‘Receptor-Ligand Pharmacophore Generation’ method and manual build feature method, and then carefully validated by a test database. Then, pharmacophore-based virtual screening (PB-VS) could be performed by using the 20 pharmacophore models. In addition, docking study was used to predict the possible binding poses of compounds, and the docking parameters were optimized before performing docking-based virtual screening (DB-VS). Moreover, a 3D QSAR model was established by applying the 55 aligned Mcl-1 inhibitors. The 55 inhibitors sharing the same scaffold were docked into the Mcl-1 active site before alignment, then the inhibitors with possible binding conformations were aligned. For the training set, the 3D QSAR model gave a correlation coefficient r² of 0.996; for the test set, the correlation coefficient r² was 0.812. Therefore, the developed 3D QSAR model was a good model, which could be applied for carrying out 3D QSAR-based virtual screening (QSARD-VS). After the above three virtual screening methods orderly filtering, 23 potential inhibitors with novel scaffolds were identified. Furthermore, we have discussed in detail the mapping results of two potent compounds onto pharmacophore models, 3D QSAR model, and the interactions between the compounds and active site residues.