- Research article
- Open Access
In-silico molecular design of heterocyclic benzimidazole scaffolds as prospective anticancer agents
BMC Chemistry volume 13, Article number: 90 (2019)
Benzimidazole is a valuable pharmacophore in the field of medicinal chemistry and exhibit wide spectrum of biological activity. Molecular docking technique is routinely used in modern drug discovery for understanding the drug-receptor interaction. The selected data set of synthesized benzimidazole compounds was evaluated for its in vitro anticancer activity against cancer cell lines (HCT116 and MCF7) by sulforhodamine B (SRB) assay. Further, molecular docking study of data set was carried out by Schrodinger-Maestro v11.5 using CDK-8 (PDB code: 5FGK) and ER-alpha (PDB code: 3ERT) as possible target for anticancer activity. Molecular docking results demonstrated that compounds 12, 16, N9, W20 and Z24 displayed good docking score with better interaction within crucial amino acids and corelate to their anticancer results. ADME results indicated that compounds 16, N9 and W20 have significant results within the close agreement of the Lipinski’s rule of five and Qikprop rule within the range and these compounds may be taken as lead molecules for the discovery of new anticancer agents.
Benzimidazole is an important structural motif found in extensive number of natural and pharmacologically active compound . The benzimidazole ring itself is an urgent pharmacophore in present day and has been used as privileged scaffolds to synthesize selective drugs of interest in medicinal field including antiulcer , antioxidant , HIV-RT inhibitor , anticancer , antihelmintic , antimicrobial , antihistamine  etc. The selected marketed drugs having benzimidazole moiety (Fig. 1) i.e. veliparib (a), glasdegib (b), liarozole (c), crenolanib (d), abemaciclib (e), pracinostat (f), bendamustine (g) nocodazole (h).
Presently used chemotherapeutic agents restrain the growth of tumour through suppression of DNA replication and transcription. Nevertheless, the attempt of discovering new curative anticancer agents in last decade has led to targets of specific molecular modifications in tumour cells. The new approach now focuses mainly on the development of small biologically active molecules containing significant activity without toxicity related to the usual chemotherapy . The human CDK8 protein which is a part of the RNA polymerase has been one of the proteins responsible for acute lymphoblastic leukaemias.
CDK activity is controlled by association with regulatory subunits (cyclins) and CDK inhibitor proteins, by their phosphorylation state and by ubiquitin-mediated proteolysis. Since the loss of cell cycle control leading to deregulated cell proliferation is one of the hallmarks of cancer, it is anticipated that the inhibition of CDKs will provide an effective approach to control tumor growth and therefore have an impact on cancer therapy. Inhibition of CDKs has been studied by many organizations and has been achieved using a variety of structural templates with varying degrees of selectivity and activities .
Western blotting analysis was used to demonstrate that three hits target CDK8 in HCT 116 colorectal cancer cells. The results showcase the successful application of virtual screening cascades to identify CDK8-targeted scaffolds that can be developed into a drug discovery program. CDK8 is a cyclin-dependent kinase that forms part of the Mediator complex, which itself regulates the transcriptional activity of RNA polymerase II. A number of studies have shown that CDK8 modulates the transcriptional output from distinct transcription factors involved in oncogenic control. These factors include the Wnt/β-catenin pathway, Notch, p53 and TGF-β. CDK8 has been found to be amplified and over expressed in colon cancer. In this context, CDK8 has been reported to act as a colon cancer oncogene. The role of CDK8 in both cellular signaling and colon cancer have relied upon RNAi mediated suppression of CDK8 and on the use of a kinase dead mutant CDK8. In order to more fully investigate the role of CDK8 in colon cancer, we aimed to develop a potent and selective small molecule inhibitor of CDK8 .
CDK-8 is a heterodimeric kinase protein responsible for regulation of cell cycle progression, transcription and other functions. CDKs require cyclin that provides additional sequences for enzymatic potential. All the CDKs (1, 4, 5, 7, 8, 9 and 11) have a two-lobed structure-N-terminal having beta sheets and C-terminal composed of α-helices [12, 13].
Estrogen signaling is essential in the initiation and development of human breast cancer. In the past several decades, extensive efforts have been dedicated to understand the underlying mechanisms of this important signaling pathway in human breast cancer, which have facilitated the development of anti-estrogen therapy, the first targeted therapy for human cancer. Estrogen biology is exceedingly complex and important in the development and function of numerous tissues and physiological phenomena [14, 15]. Computational approaches i.e. molecular docking used for modern drug discovery design for the medicinal drug . Drug molecules might fail during development because of several reasons but as found by the researchers one of the major reasons of failures is related with poor pharmacokinetic: ADME properties . Drug toxicity is the one of the major factors to withdraw drug from the market. Therefore, ADME properties are the crucial determinants for the clinical success of the drug . Now these days, computer based drug design are employed to determine the ADME profile of the compound. ADME modeling has attracted the considerable attention of the pharmaceutical researchers for the drug discovery as they are high-throughput in nature and cost effective .
Recently, a greater emphasis has been given towards the researches on complementary and alternative medicine that deals with cancer management. The present study helps us to understand the interaction between the benzimidazoles and receptors
(CDK-8 and ER-alpha) and also explore their binding mode by in silico molecular design.
Materials and methods
The data set of selected benzimidazole compounds have exhibited better anticancer activity towards human colorectal carcinoma cancer cell line (HCT116) was selected from the earlier study reported by Tahlan et al. [5, 7, 9, 20, 21]. The molecular structures of all reported compounds were drawn in ChemDraw Ultra software 12.0. The data set of selected benzimidazole compounds with their anticancer activity results is shown in Table 1 (S. No. 1 to 18).
Physicochemical and spectral interpretation data of most active compounds
M.pt. °C: > 300; Rf value: 0.52; % yield: 85.87; IR: 3077 (C–H str.), 1439 (C=C str.), 1669 (–CONH str.), 1362 (N=CH str.), 1317 (C–N str.), 730 (C–S str.), 3629 (O–H str.), 625 (C–Br str.); 1H-NMR: 3.84 (s, 2H, CH2), 4.30 (s, 1H, OH), 4.31 (s, 1H, NH of imidazole), 6.77–7.93 (m, 11H, Ar–H), 7.92 (s, 2H, (CONH)2), 7.93 (s, 1H, N=CH); 13C-NMR: 36.76, 111.0, 118.3, 118.8, 119.4, 121.0, 122.8, 125.8, 130.0, 132.1, 135.9, 139.1, 143.2, 150.1, 161.2, 158.1, 167.2; Anal. Calcd. for C23H18N5O3SBr: C, 52.68; H, 3.46; N, 13.36; Found: C, 52.72; H, 3.50; N, 13.40; MS ES + (ToF): m/z 527 [M++1].
M.pt. °C: 122–125; Rf value: 0.80; % yield: 75.42; IR: 3108 (C–H str.), 1606 (C=C str.), 1669 (–CONH str.), 1320 (N=CH str.), 1289 (C–N str.), 665 (C–S str.), 1247 (C–O–C str.), 2835 (C–H str., O–CH3); 1H-NMR: 7.04–8.26 (m, 12H, Ar–H), 4.31 (s, 1H, NH of imidazole), 3.85 (s, 2H, CH2), 7.99 (s, 2H, (–CONH)2), 8.26 (s, 1H, N=CH), 3.85 (s, 3H, CH3); 13C-NMR: 36.6, 56.2, 112.5, 117.7, 120.3, 121.1, 122.2, 123.6, 124.7, 129.5, 129.7, 133.4, 139.5, 150.1, 159.2, 165.9, 167.5; Anal. Calcd. for C24H21N5O3S: C, 62.73; H, 4.61; N, 15.24; Found: C, 62.77; H, 4.65; N, 15.28; MS ES + (ToF): m/z 462 [M++1].
M.pt. °C: 207–210; Rf value: 0.56; % yield: 64.03; IR: 3018 (C–H str.), 1598 (C=C str.), 1668 (–CONH str.), 1360 (N=CH str.), 1281 (C–N str.), 713 (C–S str.), 2915 (C–H str., –CH2–), 2948 (C–H str., CH3); 1H-NMR: 7.16–7.96 (m, 11H, Ar–H), 4.36 (s, 1H, NH of imidazole), 7.97 (s, 2H, (CONH)2), 2.53 (s, 6H, (–CH3)2); 13C-NMR: 39.0, 118.3, 121.6, 125.4, 130.4, 142.8, 149.7, 166.6, 166.8; Anal. Calcd. for C24H22N4O2S: C, 66.96; H, 5.15; N, 13.01; Found: C, 66.99; H, 5.19; N, 13.05; MS ES + (ToF): m/z 431 [M++1].
M.pt. °C: 216–219; Rf value: 0.42; % yield: 83.22; IR: 3096 (C–H str.), 1598 (C=C str.), 1336 (N=CH str.), 1304 (C–N str.), 1664 (–CONH str.), 701 (C–S str.), 2932 (C–H str., CH3), 2826 (C–H str., N–CH3), 2915 (C–H str., –CH2–); 1H-NMR: 7.12–8.33 (m, 13H, Ar–H), 4.31 (s, 1H, NH of imidazole), 7.82 (s, 2H, (CONH)2), 2.51 (q, 2H, CH2); 13C-NMR: 36.1, 119.7, 119.7, 121.4, 123.1, 124.2, 129.0, 131.3, 139.0, 149.7, 166.4, 167.0; Mol. Formula C24H22N4O2S; Elem. Anal. Calcd. C, 66.96; H, 5.15; N, 13.01; Found C, 66.93; H, 5.19; N, 13.04; MS: m/z 431 [M++1].
M.pt. °C: 160–163; Rf value: 0.51; % yield: 78.68; IR: 3087 (C–H str.), 1601 (C=C str.), 1713 (–CO str.), 1354 (C=N str.), 1337 (C–N str.), 1176 (N–N str., hydrazide), 694 (C–S str.), 2843 (C–H str., –CH2–), 684 (C–Br str.), 1623 (conjugated C=C and phenyl subst. C=C); 1H-NMR: 7.13–7.94 (m, 9H, Ar–H), 7.06 (s, 1H, NH of imidazole), 2.00 (s, 2H, CH2), 7.06 (s, 1H, NH), 12.55 (s, 1H, N=CH), 7.19 (s, 1H, Br–C=CH); 13C-NMR: 109.4, 122.1, 123.8, 128.6, 128.9, 130.5, 131.3, 132.2, 132.8, 150.1, 187.8; Mol. Formula C18H15N4OSBr; Elem. Anal. Calcd. C, 52.06; H, 3.64; N, 13.49; Found C, 52.02; H, 3.68; N, 13.45; MS: m/z 416 [M++1].
Preparation of ligand structure
The ligand structures of the data set were prepared by LigPrep module of Schrodinger v11.5. To give the best results, the structures that are docked must be good representations of the actual ligand structures as they would appear in a protein–ligand complex. This means that for Glide docking the structure must meet the following conditions. They must be three-dimensional (3D) form. Glide only modifies the torsional internal coordinates of the ligand during docking, so the rest of the geometric parameters must be optimized beforehand. They must each consist of a single molecule that has no covalent bonds to the receptor, with no accompanying fragments, such as counter ions and solvent molecules. They must have all their hydrogen (filled valences). They must have an appropriate protonation state for physiological pH values (around 7) [22, 23].
Preparation of protein structure
The selected proteins i.e. human cyclin-dependent kinase CDK-8 (PDB code: 5-FGK) and ER-alpha (PDB code: 3ERT) (Figs. 2 and 3) were obtained from the protein data bank (PDB) (http://www.rcsb.org/pdb/home/home.do). The imported typical structure file of protein from the protein data bank is not suitable for immediate use to carry out the molecular docking study. A typical PDB structure file consists of heavy atoms and may include a co-crystallized ligand, water molecules, metal ions and cofactors. The ligand and ligand-receptor complex is suitable for use with other Schrödinger modules. The protein structure was prepared using the protein preparation wizard (preprocessed, optimized and minimized) in the Schrodinger software graphical user interface Maestro v11.5 .
Preparation of grid
Grid generation is done using receptor grid generation module of maestro version 11.5. A grid is generated around the binding site already occupied by the co-crystallized ligand so that co-crystallized ligand can be excluded and new compounds can be attached to the same binding site .
Molecular docking study was applied to investigate the binding mode of compound with selected PDB ID for cancer cell lines. Docking score obtained from GLIDE (maestro v11.5) and binding site was targeted and the grid was created. The active site grid covered the important amino acids interacting with receptor. The 3D structure of the protein was obtained from protein data bank using their specific (PDB code: 5FGK and ER-alpha PDB code: 3ERT). A data set of benzimidazoles was used as ligands and their structures were drawn using the workspace of Maestro and were converted to 3D form for the docking studies. The collected ligands were prepared for docking. Then the prepared ligands were docked into the generated grid in the prepared protein. The best docked pose with lowest glide score value was recorded for each ligand. Extra precision (XP) was performed using Schrödinger-maestro v11.5 (2018-1) [26,27,28].
Theoretical calculation of the ADME properties of data set was calculated by QikProp. Nearly eleven physically significant descriptors and pharmacologically properties of the compounds were analyzed by QikProp. Aqueous solubility of compounds plays a key impact on many ADME associated properties like uptake, distribution, transport and ultimately bioavailability. The benzimidazole derivatives solubility values were within the range. Finally, the Lipinski’s rule of five and Qikprop rule within the range for the benzimidazole derivatives and thus making these derivatives as suitable drug candidates .
Anticancer evaluation (IC50)
The anticancer activity was determined by Sulforhodamine-B (SRB) assay. Briefly, HCT116 and MCF7 cancer cell lines were seeded onto the 96 well plate at 2500 cells/well. The cells were allowed to attach overnight before being exposed to the respective compounds (0.01–100 µg/mL) for 72 h. The highest concentration of each compound tested (100 µg/mL) contained only 0.1% DMSO (non-cytotoxic) SRB assay was then performed whereby the cells were fixed using trichloroacetic acid for 30 min at 4 °C and stained with 0.4% (w/v) SRB mixed with 1% acetic acid for 15 min. After five washes with 1% acetic acid solution, the protein-bound dye was extracted with 10 mM tris base solution. Optical density was read at 570 nm and IC50 of each compound was determined. Anticancer results were presented as mean IC50 of at least triplicates (Tables 3 and 4) .
Results and discussion
Kinase inhibitors are very efficacious for the treatment of cancer especially targeting specific mutations that chiefly drive tumorigenesis. They are categorized according to their capacity to catalyze the transfer of the terminal phosphate of ATP to the substrates that usually contain a serine, threonine or tyrosine residue . Cyclin-dependent kinases (CDKs) are a family of key regulatory proteins that oversee diverse cellular events and their main involvement is in the cell cycle and transcription. Given the fundamental biological roles CDKs perform, it is not surprising that their aberrant activities are a common feature of many diseases, especially cancer. CDKs are a family of serinethreonine protein kinases that govern the initiation, progression and completion of the cell cycle. Activity of the CDKs allows the orderly transition between phases of the cell cycle. Inhibition of cell cycle progression and apoptosis are the most common causes of cell growth inhibition. Cell cycle progression is induced by various cell cycle proteins such as CDKs and cyclins as they are the key regulators of cell cycle . Existing CDK8 X-ray crystal structures have unresolved regions in the vicinity of the ATP-binding site. Thus, homology modeling was used to generate two complete and optimal structural models. The target structure should be determined experimentally by either X-ray crystallography or nuclear magnetic resonance, which can be downloaded from PDB; however, docking has been performed successfully in comparison to homology models or threading. The model should have good quality. It can be tested using validation software such as Molprobity .
Molecular docking results
Molecular docking study was carried out to analyses the binding mode of the compounds against human colorectal carcinoma and breast adenocarcinoma cancer cell lines. Ligand interaction showed the binding mode of compound and standard drugs in the active site of CDK8 (PDB id: 5FGK) have good resolution about 2.36 Å, co-crystallized ligand (5XG) was selected for docking study. Root-mean-square deviation (RMSD) value of docked poses of native co-crystallized ligand is 0.10 Å, R-value free is 0.237 and ER-alpha (PDB id: 3ERT) have good resolution about 1.9 Å, co-crystallized ligand (OHT) was selected for docking study. RMSD value of docked poses of native co-crystallized ligand is 2.0 Å. R-value free is 0.262. Docking study of the data set showed good docking score and interaction with crucial amino acids residues in the binding pocket of the receptor (Table 2). The molecular docking results demonstrated in terms of negative energy value that the lower the binding energy value, best would be the binding affinity with the receptor . Docking results with cdk-8 protein, compounds 12, 16, N9, W20 and Z24 were found to be best molecules and showed better docking score at target site of protein and displayed good to moderate anticancer activity against cancer cell line (HCT116). If we look into the binding mode of compound 12, exhibited good docked score (− 8.907) and formation of hydrogen bond with amino acids i.e. Asp173, Ala100 and Lys52 with oxygen and nitrogen atoms. The binding mode of compound 16 have docking score (− 7.69) and developed of hydrogen bond with amino acids (Lys52 and Ala155) with oxygen atom of OCH3 and nitrogen atom of benzimidazole ring, respectively. The binding mode of compound N9 exhibited good docking score (− 7.425) and prepared hydrogen bond with amino acids (Lys52, Tyr32 and Val27) with oxygen atom and nitrogen atoms of benzimidazole ring, respectively. Compound W20 showed good docking score (− 9.686) and prepared hydrogen bond with amino acids (Ala100 and Ala155) with oxygen and nitrogen atoms of benzimidazole ring, respectively. Compound Z24 displayed the good docking score (− 7.295) and developed hydrogen bond with amino acids (Tyr32 and Val27) with oxygen and nitrogen atoms, respectively. Standard drugs (5-fluorouracil) have good docking score (− 5.79) and formation of hydrogen bond with amino acids residues (Ala100 and Asp98) with oxygen and nitrogen atoms of 5-fluorouracil, respectively. These compounds showed better docked score than the standard drug and the docking results also correlate to their anticancer activity results. The docking results with interacting residues of the docked compounds and standard drug are shown in Table 3; binding surface (3D) and ligand interaction (2D) images are shown in Figs. 4, 5, 6, 7, 8, 9.
Further, the selected data set was docked with ER-alpha receptor of MCF7 (PDB id: 3ERT). Among the docked data set, compounds 12, 16, N9, W20 and Z24 also showed the good docked score and glide energy with crucial amino acids residues in the binding pocket of the receptor. If we look into the binding mode of compound 12 have docking score (− 8.825) and developed hydrogen bond with the crucial amino acids (Val534, Thr347, Leu346 and Cys530). The binding mode of compound 16 have docking score (− 8.986) and developed H-bond with amino acids Thr347, Asp351 and Val534. The binding mode of compound N9 exhibited good docking score (− 6.748) and prepared hydrogen bond with crucial amino acid Asp351. Compound W20 scored docking score (− 7.703) and prepared hydrogen bond with amino acids residues Asp351. Whereas, compound Z24 displayed the good docking score (− 7.275) and prepared hydrogen bond within binding pocket. The standard drug (5-fluorouracil) have docking score (− 3.414) within the binding pocket. The docking results with interacting residues of the compounds and standard drug are shown in Table 4; binding surface (3D) and ligand interaction (2D) images are shown in Figs. 10, 11, 12, 13, 14, 15. Binding mode of five most active compounds is shown in Figs. 16 and 17. Based on the molecular docking analyses these compounds were evaluated for their in vitro anticancer activity against human breast adenocarcinoma cancer cell line (MCF7) by Sulforhodamine-B assay. The anticancer activity result of these compounds is shown in Tables 3 and 4, Figs. 18 and 19.
Molecular docking results suggest that the selected compounds of heterocyclic benzimidazole can act as of great interest in successful chemotherapy. The selected protein data bank i.e. (PDB id: 5FGK and 3ERT) for human colorectal carcinoma and breast adenocarcinoma cancer cell lines may be the good target protein of benzimidazole molecules for their anticancer activity. Based on the docking analysis it is suggested that more structural modifications are required in molecules 12, 16, N9, W20 and Z2 to make them more potent toward cancer cell. The structure activity relationship study based on molecular doking is given in Figs. 20 and 21.
Lipinski’s rule of five is a rule of thumb to evaluate drug likeness or determine if a chemical compound with a certain pharmacological or biological activity has chemical properties and physical properties that would make it a likely orally active drug in humans. The rule describes molecular properties important for a drug’s pharmacokinetics in the human body, including their absorption, distribution, metabolism, and excretion (ADME). The rule is important to keep in mind during drug discovery when a pharmacologically active lead structure is optimized step-wise to increase the activity and selectivity of the compound as well as to ensure drug-like physicochemical properties are maintained as described by Lipinski’s rule which states that (i) no more than 5 hydrogen bond donors, (ii) no more than 10 hydrogen bond acceptors, (iii) a molecular mass less than 500 daltons, (iv) an octanol–water partition coefficient log P not greater than 5 (https://en.wikipedia.org/wiki/Lipinski%27s_rule_of_five). Now these days, computational approaches are employed to determine the ADME of the drug molecules. ADME modeling has attracted the considerable attention of the pharmaceuticals researchers for the drug discovery as they are high-throughput in nature and cost effective . ADME study of the selected compounds was performed using QikProp module of Maestro version 11.5. Around eleven physically relevant and pharmacologically significant parameters of the most active compounds 16, N9 and W20 were determined. The ADME results of the selected heterocyclic benzimidazole compounds 16, N9 and W20 displayed the significant results within the close agreement with the Lipinski’s rule of five and Qikprop rule within the range i.e. Molecular weight of the molecule (mol. MW ≤ 500), Predicted octanol/water partition coefficient (QPlogPo/w = − 2.0 to − 6.5), (QPlogPw = 4.0 to − 45.0), Predicted water/gas partition coefficient (QPlogKp = − 8.0 to − 1.0), Predicted brain/blood partition coefficient (QPlogBB = − 3.0 to − 1.2), donor HB (0.0 to − 6.0), accept HB (2.0 to − 20.0), human oral absorption (1, 2 or 3), percent human oral absorption (0 to 100), Predicted water/gas partition coefficient thus making these compounds as suitable drug candidate. The ADME results are shown in the in Table 5.
In the present work the molecular docking study of the data set of heterocyclic benzimidazole molecules was performed by Maestro version 11.5. In this study we have used CDK8 for human colorectal carcinoma cancer and ER-alpha for breast adenocarcinoma cancer cell lines and identified best anticancer target. Molecular docking results displayed that compounds (12, 16, N9, W20 and Z24) showed the better docked score with moderate to better antiproliferative potency towards cancer cell line within the binding pocket of receptor and comparable to the standard drug. The compounds 16, N9 and W20 also showed good ADME properties within the close agreement of the Lipinski’s rule of five and Qikprop rule within the range and thus making these compounds as suitable drug molecules.
Availability of data and materials
We have presented all our main data in the form of tables and figures.
absorption, distribution, metabolism and excretion
cyclin dependent kinase 8
estrogen receptor alpha
human colorectal carcinoma 116
Michigan Cancer Foundation-7
protein data bank
Abonia R, Cortes E, Insuasty B, Quiroga J, Nogueras M, Cobo J (2011) Synthesis of novel 1,2,5-trisubstituted benzimidazoles as potential antitumor agents. Eur J Med Chem 46:4062–4070
Patil A, Ganguly S, Surana S (2010) Synthesis and antiulcer activity of 2-[5-substituted-1-H-benzo(d)imidazol-2-yl sulfinyl]methyl-3-substituted quinazoline-4-(3H) ones. J Chem Sci 122(3):443–450
Kus C, Ayhan-Kılcıgil G, Ozbey S, Kaynak FB, Kaya M, Coban T, Can-Eke B (2008) Synthesis and antioxidant properties of novel N-methyl-1,3,4-thiadiazol-2-amine and 4-methyl-2H-1,2,4-triazole-3(4H)-thione derivatives of benzimidazole class. Bioorg Med Chem 16:4294–4303
Ziolkowska NE, Michejda CJ, Bujacz GD (2009) Crystal structures of HIV-1 nonnucleoside reverse transcriptase inhibitors: N-benzyl-4-methyl-benzimidazoles. J Mol Strut 930:157–161
Tahlan S, Narasimhan B, Lim SM, Ramasamy K, Mani V, Shah SAA (2018) 2-Mercaptobenzimidazole Schiff bases: design, synthesis, antimicrobial studies and anticancer activity on HCT-116 cell line. Mini-Rev Med Chem. https://doi.org/10.2174/1389557518666181009151008
Hernandez-Covarrubias C, Vilchis-Reyes MA, Yepez-Mulia L, Sanchez-Diaz R, Navarrete-Vazquez G, Hernandez-Campos A, Castillo R, Hernandez-Luis F (2012) Exploring the interplay of physicochemical properties, membrane permeability and giardicidal activity of some benzimidazole derivatives. Eur J Med Chem 52:193–204
Tahlan S, Narasimhan B, Lim SM, Ramasamy K, Mani V, Shah SAA (2018) Design, synthesis, SAR study, antimicrobial and anticancer evaluation of novel 2-mercaptobenzimidazole azomethine derivatives. Mini-Rev Med Chem. https://doi.org/10.2174/1389557518666180903151849
Lavrador-Erb K, Ravula SB, Yu J, Zamani-Kord S, Moree WJ, Petroski RE, Wen J, Malany S, Hoare SRJ, Madan A, Crowe PD, Beaton G (2010) The discovery and structure–activity relationships of 2-(piperidin-3-yl)-1Hbenzimidazoles as selective, CNS penetrating H1-antihistamines for insomnia. Bioorg Med Chem Lett 20:2916–2919
Tahlan S, Ramasamy K, Lim SM, Shah SAA, Mani V, Narasimhan B (2019) 4-(2-(1H-Benzo[d]imidazol-2-ylthio)acetamido)-N-(substituted phenyl) benzamides: design, synthesis and biological evaluation. Chem Cent J 13(12):1–16
Chu X-J, DePinto W, Bartkovitz D, So S-S, Vu BT, Packman K, Lukacs C, Ding Q, Jiang N, Wang K, Goelzer P, Yin X, Smith MA, Higgins BX, Chen Y, Xiang Q, Moliterni J, Kaplan G, Graves B, Lovey A, Fotouhi N (2006) Discovery of [4-amino-2-(1-methanesulfonyl-piperidin-4-ylamino)pyrimidin-5-yl](2,3-difluoro-6-methoxyphenyl)methanone (R547), a potent and selective cyclin-dependent kinase inhibitor with significant in vivo antitumor activity. J Med Chem 49:6549–6560
Koehler MFT, Bergeron P, Blackwood EM, Bowman K, Clark KR, Firestein R, Kiefer JR, Maskos K, McCleland ML, Orren L, Salphati L, Schmidt S, Schneider FV, Wu J, Beresini MH (2016) Development of a potent, specific CDK8 kinase inhibitor which phenocopies CDK8/19 knockout cells. ACS Med Chem Lett 7:223–228
Kalra S, Joshi G, Munshi A, Kumar R (2017) Structural insights of cyclin dependent kinases: implications in design of selective inhibitors. Eur J Med Chem 142:424–458
Kumar S, Lim SM, Ramasamy K, Vasudevan M, Shah SAA, Selvaraj M, Narasimhan B (2017) Synthesis, molecular docking and biological evaluation of bis-pyrimidine Schiff base derivatives. Chem Cent J 11(89):1–16
Wang Z-Y, Yin L (2015) Estrogen receptor alpha-36 (ER-α36): a new player in human breast cancer. Mol Cell Endocrinol 418(3):193–206
Kumar R, Zakharov MN, Khan SH, Miki R, Jang H, Toraldo G, Singh R, Bhasin S, Ravi J (2011) The dynamic structure of the estrogen receptor. J Amino Acids 2011:1–7
Kumar S, Narasimhan B, Lim SM, Ramasamy K, Mani V, Shah SAA (2019) Design, synthesis and biological potential of 5-(2-amino-6-(3/4-bromophenyl)pyrimidin-4-yl) benzene-1,3-diol scaffolds as promising antimicrobial and anticancer agents. Mini-Rev Med Chem. https://doi.org/10.2174/1389557518666181009141924
Wang J, Urban L (2004) The impact of early ADME profiling on drug discovery and development strategy. Drug Discov World Fall 5(4):73–86
Li AP (2001) Screening for human ADME/Tox drug properties in drug discovery. Res Focus 6(7):357–366
Wang Y, Xing J, Xu Y, Zhou N, Peng J, Xiong Z, Liu X, Luo X, Luo C, Chen K, Zheng M, Jiang H (2015) In silico ADME/T modelling for rational drug design. Q Rev Biophys. https://doi.org/10.1017/s0033583515000190:1-28
Tahlan S, Ramasamy K, Lim SM, Shah SAA, Mani V, Narasimhan B (2018) Design, synthesis and therapeutic potential of 3-(2-(1H-benzo[d]imidazol-2-ylthio) acetamido)-N-(substituted phenyl)benzamide analogues. Chem Cent J 12(139):1–12
Tahlan S, Kumar S, Ramasamy K, Lim SM, Shah SAA, Mani V, Pathania R, Narasimhan B (2019) Design, synthesis and biological profile of heterocyclic benzimidazole analogues as prospective antimicrobial and antiproliferative agents. BMC Chem 13(50):1–15
Driessche GVD, Fourches D (2017) Adverse drug reactions triggered by the common HLA-B*57:01 variant: a molecular docking study. J Cheminform 9(13):1–17
Sastry GM, Adzhigirey M, Day T, Annabhimoju R, Sherman W (2013) Protein and ligand preparation: parameters, protocols and influence on virtual screening enrichments. J Comput Aid Mol Des 27(3):221–234
Kumar S, Singh J, Narasimhan B, Shah SAA, Lim SM, Ramasamy K, Mani V (2018) Reverse pharmacophore mapping and molecular docking studies for discovery of GTPase HRas as promising drug target for bis-pyrimidine derivatives. Chem Cent J 12(106):1–11
Sharma V, Sharma PC, Kumar V (2016) In silico molecular docking analysis of natural pyridoacridines as anticancer agents. Adv Chem 2016:1–9
Singh J, Kumar M, Mansuri R, Sahoo GC, Deep A (2016) Inhibitor designing, virtual screening and docking studies for methyltrans-ferase: a potential target against dengue virus. J Pharm Bioallied Sci 8(3):188–194
Friesner RA, Murph RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, Sanschagrin PC, Mainz DT (2006) Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein–ligand complexes. J Med Chem 49:6177–6196
Lenselink EB, Louvel J, Forti AF, van Veldhoven JPD, de Vries H, Mulder-Krieger T, McRobb FM, Negri A, Goose J, Abel R, van Vlijmen HWT, Wang L, Harder E, Sherman W, IJzerman AP, Beuming T (2016) Predicting binding affinities for GPCR ligands using free-energy perturbation. ACS Omega 1:293–304
Skehan P, Storeng R, Scudiero D, Monks A, McMahon J, Vistica D, Warren JT, Bokesch H, Kenney S, Boyd MR (1990) New colorimetric cytotoxicity assay for anticancer-drug screening. J Natl Cancer Inst 82:1107–1112
Bhullar KS, Lagarón NO, McGowan EM, Parmar I, Jha A, Hubbard BP, Rupasinghe HPV (2018) Kinase-targeted cancer therapies: progress, challenges and future directions. Mol Cancer 17(48):1–20
Shaik TB, Malik SM, Seddigid ZS, Routhu SR, Kamala A (2017) Evaluation of anticancer and anti-mitotic properties of quinazoline and quinazolino-benzothiadiazine derivatives. http://dx.doi.org/10.1101/108654
Hernández-Santoyo A, Tenorio-Barajas AY, Altuzar V, Vivanco-Cid H, Mendoza-Barrera C (2013) Protein-protein and protein-ligand docking. Protein engineering-technology and application, Tomohisa Ogawa. IntechOpen. https://doi.org/10.5772/56376
Bassyouni F, El Hefnawi M, El Rashed A, Rehim MA (2017) Molecular modeling and biological activities of new potent antimicrobial, anti-inflammatory and anti-nociceptive of 5-nitro indoline-2-one derivatives. Drug Des 6(2):1–6
Sharma D, Kumar S, Narasimhan B, Ramasamy K, Lim SM, Shah SAA, Mani V (2019) 4-(4-Bromophenyl)-thiazol-2-amine derivatives: synthesis, biological activity and molecular docking study with ADME profile. BMC Chem 13(60):1–16
The authors are thankful to Head, Department of Pharmaceutical Sciences, Maharshi Dayanand University, Rohtak, for providing necessary facilities to carry out this research work.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tahlan, S., Kumar, S., Ramasamy, K. et al. In-silico molecular design of heterocyclic benzimidazole scaffolds as prospective anticancer agents. BMC Chemistry 13, 90 (2019) doi:10.1186/s13065-019-0608-5
- Anticancer activity