Volume 11 Supplement 7
IKKβ inhibitor identification: a multi-filter driven novel scaffold
© Nagarajan et al; licensee BioMed Central Ltd. 2010
Published: 15 October 2010
Nuclear factor kappa B (NF-κB) is a chief nuclear transcription factor that controls the transcription of various genes; and its activation is tightly controlled by Inhibitor kappa B kinase (IKK). The irregular transcription of NF-κB has been linked to auto-immune disorders, cancer and other diseases. The IKK complex is composed of three units, IKKα, IKKβ, and the regulatory domain NEMO, of which IKKβ is well understood in the canonical pathway. Therefore, the inhibition of IKKβ by drugs forms the molecular basis for anti-inflammatory drug research.
The ligand- and structure-based virtual screening (VS) technique has been applied to identify IKKβ inhibitors from the ChemDiv database with 0.7 million compounds. Initially, a 3D-QSAR pharmacophore model has been deployed to greatly reduce the database size. Subsequently, recursive partitioning (RP) and docking filters were used to screen the pharmacophore hits. Finally, 29 compounds were selected for IKKβ enzyme inhibition assay to identify a novel small molecule inhibitor of IKKβ protein.
In the present investigation, we have applied various computational models sequentially to virtually screen the ChemDiv database, and identified a small molecule that has an IC50 value of 20.3μ M. This compound is novel among the known IKKβ inhibitors. Further optimization of the hit compound can reveal a more potent anti-inflammatory agent.
Inhibitor kappa-B kinaseβ (IKKβ) is a serine-threonine protein kinase, which is critically involved in the activation of transcription factor Nuclear Factor kappa B (NF-κB) in response to various inflammatory stimuli . IκB, an inhibitory unit, is responsible for retaining NF-κB in the cytoplasm , for the degradation of IκB by phosphorylation, and for ubiquitination to translocate NF-κB into the nucleolus, leading to transcription initiation . IKKβ plays a crucial role in the way of canonical NF-κB pathway, which phosphorylates IκB protein and thereby translocates NF-κB into the nucleus and initiates pro-inflammatory gene transcription. The canonical NF-κB pathway is well recognized in chronic inflammatory diseases  and inhibition of the IKKβ enzyme by a highly potent inhibitor has remained the primary goal for anti-inflammatory drug discovery.
The IKK complex comprises two catalytic subunits, IKKα and IKKβ, and a regulatory subunit, IKKγ. Although both the catalytic subunits can catalyze the phosphorylation of IκBα, the IKKβ subunit seems to play a dominant role in the canonical pathway. Furthermore, IKKα has a crucial role in mediating p52 activation through the 'non-canonical' pathway . IKKα can form an alternative complex (without IKKβ and IKKγ) and its function is required for the development of the lymphoid organ and the maturation of B cells . Termination of the canonical pathway by inhibiting IKKβ is a potential target in anti-inflammatory drug research.
Recently, the virtual screening (VS) method is playing an increasingly important role in drug discovery. The structure-based method involves docking of small molecules and ranking them based on their score. Every scoring function has its own inherent limitations, and thus, there is a high chance for reporting false positives. In order to minimize the risks of using a structure-based approach, additional filters have been used to enrich the VS scheme. The application of various computational filters in the VS cascade certainly alleviates the difficulties encountered during the initial stages of the drug discovery process. Every model used in the VS scheme has been meticulously validated by test sets that are not included in training the models. In general, the performance of the model is highly dependent on the choice of the ligand that used to train the model.
Results and discussions
3D-QSAR pharmacophore model
Recursive partitioning model
Statistical analysis of the RP model.
Summary of the descriptors that were found to be useful in decision making in the RP model
Sum descriptor for carbon with two single bonds.
Number of hydrogen bond acceptor
Prediction of CYP2D6 inhibition.
Molecular surface area.
The strength and orientation behaviour of a molecule in an electrostatic field.
Number of hydrogen bond donor
Sum descriptor for nitrogen with two aromatic bonds.
The chemical similarity between the VH02 and BMS-34541 provides a basic intuition for the chemical modification of this hit compound. The benzaldehyde moiety of VH02 can be replaced by tiny hydrophobic moieties, whereas, the phenol moiety can be replaced by pyrrole, that can maintain the same distance constraint for nitrogen as that of the BMS compound, to facilitate hydrogen bond formation between the NH group of the ligand and the receptor.
We have developed a filter-driven scaffold model and applied it for the virtual screening of IKKβ inhibitors. Sequential filtering of the database can reduce the false positive rate to a large extent at each stage. The first two models are generated by means of using the known inhibitor information and the third model is a structure-based approach. At the initial level of screening, IKKβ inhibitor-like compounds are retained, and allowed to pass on to the structure-based filter. Docking of several compounds simultaneously to the IKKβ active site revealed the set of compounds that are stable at the ATP binding pocket. In general, identification of lead molecules using a computational modeling approach often relies on approximation and has limited accuracy. Therefore, the VS hits have been validated further by subjecting them to in vitro studies.
The VS approach reported 367 hits; and among these compounds, only 29 have been selected based on encouraging scores, diversity, and commercial availability for the IKKβ inhibition assay. Of the 29 compounds tested, we have identified one hit (VH01) with IC50 20.3 μ M. Despite this inhibition value, this compound is found to be structurally novel among reported IKKβ inhibitors. There are series of similar compounds patented by Zhuravel et al. , which interestingly, also seem to exhibit antitumor activity. Hideshima et al.  have previously explained the use of a small molecule inhibitors of IKKβ and its role in inhibiting the haematological cancer, multiple myeloma. Accordingly, we will focus our attention on the anti-cancer point of view with the identified hit compound. Further optimization of VH01 can lead us to discover more potent compounds that can act as anti-inflammatory as well as anti-cancer agents, and this work currently underway. Although the VH02 compound has not been found to be very potent, its similarity to BMS-345541 has suggested that the screening system could bring out the core features required to be present in the IKKβ inhibitor. Moreover, the VS cascade is not based on serendipity, as it has proven its efficiency in identifying IKKβ inhibitors.
Pharmacophore model generation
The pharmacophore hypothesis modeling was performed using the Catalyst 4.11 (Accelrys, 9685 Scranton Road, San Diego, Calif. 92121) package. A total of 159 compounds collected from the literature [4, 14–17], was made into a library. Subsequently, the library was divided into training and test sets composed of 23 and 136 compounds, respectively. From each scaffold category, a few representative compounds were chosen based on diverse substituents with a wide range of activity data. The 3D-QSAR pharmacophore model known as Hypogen was generated based on 23 IKKβ inhibitors, whose activity data ranges from 3 nM ≤ IC50 ≤ 50000 nM. Detailed information about the pharmacophore can be found elsewhere . The training set compounds were broadly classified into four groups: those with an activity range < 100 nM were classified as highly active (+++); an activity range between > 100 nM to < 1 μM were defined as active (++); compounds with an activity range of > 1 μM to < 10 μM were defined as moderately active (+); and, the compounds having an IC50 value > 10 μM were classified as inactive (-). The same grouping strategy was applied to the test set compounds also. Excluding the training set compounds, the remaining compounds were used as an internal test set to measure the efficiency of the pharmacophore model; no outliers were removed to achieve unrealistic higher correlation values. These compounds also covered a wide range of activity of 4 nM ≤ IC50 ≤ 50000 nM.
For every training set compound, all possible conformers were enumerated and a spreadsheet was prepared with the corresponding activity data and conformers. Additional specifications were made to select desired features, such as hydrogen bond donors, hydrogen bond acceptors, hydrophobes and aromatic rings. The spreadsheet was input to the Catalyst program and in a reasonable time-frame, 10 hypotheses were generated. The best pharmacophore model was selected based on highest correlation, lowest RMSD and the most significant cost values.
Decision tree generation
The RP method of the Cerius2 program was used to generate a decision tree. RP is a classification structure-activity relation (CSAR) method that enables rapid classification of large databases, is non-parametric and captures nonlinear relationships automatically performed based on the Classification and Regression Trees (CART) algorithm . The working principle behind the RP is assembling a set of descriptors, converting them into a data object to reflect the presence or absence of useful features, assembling the data objects into vectors and then into a matrix. Finally, the matrix is divided into two daughter sets, based on the presence/absence of certain useful features. The process is repeated until each member of the matrix has been designated to a terminal node based on the presence/absence of specified features. The RP model is found to be sensitive to the descriptors used, and diversity of the data sets can radically change the property of the decision tree. The method is applicable to structurally unique compounds with activity data to uncover sub-structural rules that govern the biological activity . The RP classification tree is often of great interest to visualize the distribution of potencies at the node and to see how a split at a node divides the potencies at two daughter nodes. This method has been repeatedly used by researchers of bioinformatics and chemoinformatics, either to classify genes or to differentiate active and inactive compounds [20–23]. However, the limitation of the RP method is its inability to extrapolate beyond the range of observed responses. The main objective of incorporating the RP method in the virtual screening process is to rapidly classify unknown compounds based on a small number of readily interpretable descriptors; therefore, for screening compounds.
The recursive partition decision tree model was constructed using a QSAR module of Cerius2 version 4.10.17. The splits were scored using the Gini Impurity scoring function, which minimizes the impurity of the nodes resulting from the split. The tree was set to prune backward through a moderate pruning process, to avoid over splitting. Every node should contain 1% of the samples to qualify for further splits. The knot value was limited to a threshold of 20 per variable and maximum tree depth was set to 10. The best RP tree was generated with these parameters.
Training and test sets of the RP model
Molecular descriptors used for recursive partition model development
MW, Rotlbonds, Hbond acceptor, Hbond Donor
Electrotopological (E-state key)
S_sCH3, S_ssCH2, S_aaCH, S_sssCH, S_tsC, S_dssC, S_aasC, S_sNH2, S_ssNH, S_tN, S_aaN, S_sssN, S_sOH, S_dO, S_ssO, S_sBr
RadOfGyration, Area, Vm, Density, PMI
AlogP98, logP, Molref, Fh2O, Foct
ADMET_Abos, ADMET_Solu, ADMET_Hepatox, ADMET_CYP2D6, ADMET_PPB
A total of 84 compounds were used as an external test set compounds, collected from a different set of published articles, with none of the compounds or similar scaffolds included in the training set. External test set compounds have been reported by two groups [4, 30]. The first set of compounds are derivatives of the imidazothienopyrazine core , with a series of compounds having imidazoquinoxaline  synthesized by same group included in training the model. Another set of compounds reported by Chiristoper et al. , was synthesized based on the benzimidazole core to specifically inhibit IKKμ, but instead inhibited IKKβ. The external test sets were combined to serve as an independent test set to asses the generality of the model. Dependent and independent variables were calculated as explained before.
The templates used to model IKKβ protein
Resolution in Å
Crystal Structure of Human Calmodulin-dependent protein kinase 1 D.
Calmodulin-dependent protein kinase from rat.
Kinase and Ubiquitin-associated domains of MARK3/Par-1
Protein kinase A - fivefold mutant model of Rho-kinase
The FlexX docking program  was employed in the structure-based VS. Prior to docking, hydrogen atoms were added to the protein, and it was minimized using the steepest descent algorithm for about 500 steps. The amino acids Phe26, Val29, Ala42, Lys44, Met65, Val74, Ala76, Glu97, Tyr98, Cys99, Lys106, Val152, and Gly168 and the surrounding residues within the distance range of 6.5 Å were defined as active site. FlexX uses an incremental construction algorithm to place flexible ligands into a fully specified active site, while its empirical scoring function estimates the binding free energy based on physicochemical properties. The FlexX-Pharm  module was used to define the constraints and direct the FlexX docking of several compounds into the specified active site simultaneously. FlexX-Pharm ensures that an interaction is formed between the specified interacting group in the active site and the ligand in a valid docking solution. There are many research groups, who have successfully employed constraints in structure-based VS to increase the enrichment factor [33–36] of active compounds. As we know, most of the ATP competing kinase protein inhibitors make two or three hydrogen bonding interactions with the hinge region [37–39]. Hence, we applied hydrogen bonding as constraints to select compounds that can possibly make hinge interactions. In the docking simulation, two different sets of constraints were applied; namely, 'heavy' and 'light'. The heavy constraint method is very strict in choosing compounds. According to this method, compounds forming three hydrogen bonds with the hinge region (Acceptor (Glu97)-Donor (Cys99)-Acceptor (Cys99)) were alone reported as hits. In the light constraint method, the middle donor interaction is essential and at least one acceptor hydrogen bonding interaction is essential. A maximum of 30 conformers were retained for each compound, passing the constraints criteria. In our previous work , we demonstrated that the f-scoring function was good enough to discriminate IKKβ inhibitors from decoys and so, the same scoring function has been applied in this VS scheme.
In vitro analysis: IKKβ enzyme inhibition assay
IKKβ-TR-FRET reactions for the search of IKKβ inhibitors were carried out based upon the suggestions of the IMAP-TR-FRET system (MDS Analytical Technologies, Sunnyvale, Calif., USA). IKKβ kinase reactions were performed in a reaction buffer (10 mM Tris-HCl, pH 7.2, 10 mM MgCl2, 0.05% NaN3), containing 1 mM DTT and 0.01% Tween-20 (Sigma-Aldrich Co., St. Louis, Mo., USA) to help stabilize the enzyme. The reactions were done at room temperature for 2 h in white standard 384 plates (3572, Corning Life Sciences, Lowell, Mass., USA), using 0.5 μg/ml IKKβ (Millipore Co., Billerica, Mass., USA), 1 μM IKBα-derived substrate (5FAM-GRHDSGLDSMK-NH2; R7574, MDS Analytical Technologies), and 3 μM ATP (Sigma-Aldrich Co.) unless otherwise noted. The total reaction volumes were 20 μl and 10 μM, and compounds were preincubated with the IKKβ enzyme for 10 min before the substrate and ATP were added. For the TR-FRET reaction, 60 μl of the detection mixture (1:600 dilution of IMAP binding reagent and 1:400 dilution of Terbium donor supplied by MDS Analytical Technologies) were added 15 h before reading the plate. The energy transfer signal was measured in a multilabel counter with a TR-FRET option (Victor II, PerkinElmer Oy, Turku, Finland). The counter setting was 340 nm excitation, 100 μs delay, and dual-emission collection for 200 μs at 495 and 520 nm. The energy transfer signal data were used to calculate the percentage inhibition and IC50 values. To monitor the assay system and to compare the hit compounds, Bayer compound (Bayer -5a) was used as a positive control.
List of abbreviations used
3D-Quantitative Structure Activity Relationship
Inhibitor kappa B Kinase
Nuclear Factor kappa B
Matthews correlation coefficient
Hydrogen bond donor
This work is supported by the Korea Institute of Science and Technology (KIST).
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 7, 2010: Ninth International Conference on Bioinformatics (InCoB2010): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S7.
- Mercurio F, Zhu HY, Murray BW, Shevchenko A, Bennett BL, Li JW, Young DB, Barbosa M, Mann M: IKK-1 and IKK-2: Cytokine-activated I kappa B kinases essential for NF-kappa B activation. Science 1997, 278: 860–866. 10.1126/science.278.5339.860View ArticlePubMedGoogle Scholar
- Beg AA, Baldwin AS: The I kappa B proteins: multifunctional regulators of Rel/NF-kappa B transcription factors. Genes Dev 1993, 7: 2064–2070. 10.1101/gad.7.11.2064View ArticlePubMedGoogle Scholar
- Scherer DC, Brockman JA, Chen ZJ, Maniatis T, Ballard DW: Signal-Induced Degradation of I-Kappa-B-Alpha Requires Site-Specific Ubiquitination. Proceedings of the National Academy of Sciences of the United States of America 1995, 92: 11259–11263. 10.1073/pnas.92.24.11259PubMed CentralView ArticlePubMedGoogle Scholar
- Christopher JA, Avitabile BG, Bamborough P, Champigny AC, Cutler GJ, Dyos SL, Grace KG, Kerns JK, Kitson JD, Mellor GW, et al.: The discovery of 2-amino-3,5-diarylbenzamide inhibitors of IKK-alpha and IKK-beta kinases. Bioorganic & Medicinal Chemistry Letters 2007, 17: 3972–3977. 10.1016/j.bmcl.2007.04.088View ArticleGoogle Scholar
- Strnad J, Burke JR: I kappa B kinase inhibitors for treating autoimmune and inflammatory disorders: potential and challenges. Trends in Pharmacological Sciences 2007, 28: 142–148. 10.1016/j.tips.2007.01.005View ArticlePubMedGoogle Scholar
- Lawrence T, Bebien M, Liu GY, Nizet V, Karin M: IKK alpha limits macrophage NF-kappa B activation and contributes to the resolution of inflammation. Nature 2005, 434: 1138–1143. 10.1038/nature03491View ArticlePubMedGoogle Scholar
- Nagarajan S, Ahmed A, Choo H, Cho YS, Oh KS, Lee BH, Shin KJ, Pae AN: 3D QSAR pharmacophore model based on diverse IKKβ inhibitors. Journal of Molecular Modeling 2010. 10.1007/s00894–010–0714–8 10.1007/s00894-010-0714-8Google Scholar
- Lamanna C, Bellini M, Padova A, Westerberg G, Maccari L: Straightforward recursive partitioning model for discarding insoluble compounds in the drug discovery process. Journal of Medicinal Chemistry 2008, 51: 2891–2897. 10.1021/jm701407xView ArticlePubMedGoogle Scholar
- Landis JR, Koch GG: The measurement of observer agreement for categorical data. Biometrics 1977, 33: 159–174. 10.2307/2529310View ArticlePubMedGoogle Scholar
- Murata T, Shimada M, Kadono H, Sakakibara S, Yoshino T, Masuda T, Shimazaki M, Shintani T, Fuchikami K, Bacon KB, et al.: Synthesis and structure-activity relationships of novel IKK-beta inhibitors. Part 2: Improvement of in vitro activity. Bioorganic & Medicinal Chemistry Letters 2004, 14: 4013–4017. 10.1016/j.bmcl.2004.05.040View ArticleGoogle Scholar
- Nagarajan S, Choo H, Cho YS, Oh KS, Lee BH, Shin KJ, Pae AN: IKKβ inhibitors Identification Part II: Ligand and Structure Based Virtual Screening. Bioorganic & Medicinal Chemistry 2010.Google Scholar
- Zhuravel IAI, Kovalenko AV, Borisov SN, Matusovich AV, Il'ya Okun, Tkachenko SE, Chernykh VE: Substituted pyrido[4',3':5,6]pyrano[2,3-d]pyrimidines and combinatorial library of these derivatives with antitumor activity. Book Substituted pyrido[4',3':5,6]pyrano[2,3-d]pyrimidines and combinatorial library of these derivatives with antitumor activity. (Editor ed.^eds.) City 2006, 13.Google Scholar
- Hideshima T, Chauhan D, Richardson P, Mitsiades C, Mitsiades N, Hayashi T, Munshi N, Dang L, Castro A, Palombella V, et al.: NF-kappa B as a therapeutic target in multiple myeloma. Journal of Biological Chemistry 2002, 277: 16639–16647. 10.1074/jbc.M200360200View ArticlePubMedGoogle Scholar
- Baxter A, Brough S, Cooper A, Floettmann E, Foster S, Harding C, Kettle J, McInally T, Martin C, Mobbs M, et al.: Hit-to-lead studies: the discovery of potent, orally active, thiophenecarboxamide IKK-2 inhibitors. Bioorganic & Medicinal Chemistry Letters 2004, 14: 2817–2822. 10.1016/j.bmcl.2004.03.058View ArticleGoogle Scholar
- Beaulieu F, Ouellet C, Ruediger EH, Belema M, Qiu YP, Yang XJ, Banville J, Burke JR, Gregor KR, MacMaster JF, et al.: Synthesis and biological evaluation of 4-amino derivatives of benzimidazoquinoxaline, benzimidazoquinoline, and benzopyrazoloquinazoline as potent IKK inhibitors. Bioorganic & Medicinal Chemistry Letters 2007, 17: 1233–1237. 10.1016/j.bmcl.2006.12.017View ArticleGoogle Scholar
- Morwick T, Berry A, Brickwood J, Cardozo M, Catron K, DeTuri M, Emeigh J, Homon C, Hrapchak M, Jacober S, et al.: Evolution of the thienopyridine class of inhibitors of I kappa B kinase-beta: Part I: Hit-to-lead strategies. Journal of Medicinal Chemistry 2006, 49: 2898–2908. 10.1021/jm0510979View ArticlePubMedGoogle Scholar
- Sugiyama H, Yoshida M, Mori K, Kawamoto T, Sogabe S, Takagi T, Oki H, Tanaka T, Kimura H, Ikeura Y: Synthesis and structure activity relationship studies of benzothieno[3,2-b]furan derivatives as a novel class of IKK beta inhibitors. Chemical & Pharmaceutical Bulletin 2007, 55: 613–624. 10.1248/cpb.55.613View ArticleGoogle Scholar
- De'ath G, Fabricius KE: Classification and regression trees: A powerful yet simple technique for ecological data analysis. Ecology 2000, 81: 3178–3192. 10.1890/0012-9658(2000)081[3178:CARTAP]2.0.CO;2View ArticleGoogle Scholar
- Farmen MWL, Rusinko Christophe III, Young AndrewR, Stanley Sidney: Statistical deconvoluting of mixtures. Book Statistical deconvoluting of mixtures. (Editor ed.^eds.), City 2002., 6434542:Google Scholar
- Huang LT, Gromiha MM, Ho SY: iPTREE-STAB: interpretable decision tree based method for predicting protein stability changes upon mutations. Bioinformatics 2007, 23: 1292–1293. 10.1093/bioinformatics/btm100View ArticlePubMedGoogle Scholar
- Houseman EA, Christensen BC, Yeh RF, Marsit CJ, Karagas MR, Wrensch M, Nelson HH, Wiemels J, Zheng S, Wiencke JK, Kelsey KT: Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. Bmc Bioinformatics 2008., 9: 10.1186/1471-2105-9-365Google Scholar
- Chohan KK, Paine SW, Mistry J, Barton P, Davis AM: A rapid computational filter for cytochrome P450 1A2 inhibition potential of compound libraries. Journal of Medicinal Chemistry 2005, 48: 5154–5161. 10.1021/jm048959aView ArticlePubMedGoogle Scholar
- Burton J, Ijjaali I, Barberan O, Petitet F, Vercauteren DP, Michel A: Recursive partitioning for the prediction of cytochromes P450 2D6 and 1A2 inhibition: Importance of the quality of the dataset. Journal of Medicinal Chemistry 2006, 49: 6231–6240. 10.1021/jm060267uView ArticlePubMedGoogle Scholar
- Crius2: (version 410), Accelrys, Inc, San Diego, CA, USA. 2005.Google Scholar
- Murata T, Shimada M, Sakakibara S, Yoshino T, Masuda T, Shintani T, Sato H, Koriyama Y, Fukushima K, Nunami N, et al.: Synthesis and structure-activity relationships of novel IKK-beta inhibitors. Part 3: Orally active anti-inflammatory agents. Bioorganic & Medicinal Chemistry Letters 2004, 14: 4019–4022. 10.1016/j.bmcl.2004.05.041View ArticleGoogle Scholar
- Bingham AH, Davenport RJ, Gowers L, Knight RL, Lowe C, Owen DA, Parry DM, Pitt WR: A novel series of potent and selective IKK2 inhibitors. Bioorganic & Medicinal Chemistry Letters 2004, 14: 409–412. 10.1016/j.bmcl.2003.10.047View ArticleGoogle Scholar
- Bonafoux D, Bonar S, Christine L, Clare M, Donnelly A, Guzova J, Kishore N, Lennon P, Libby A, Mathialagan S, et al.: Inhibition of IKK-2 by 2-[(aminocarbonyl)amino]-5-acetylenyl-3-thiophenecarboxamides. Bioorganic & Medicinal Chemistry Letters 2005, 15: 2870–2875. 10.1016/j.bmcl.2005.03.090View ArticleGoogle Scholar
- Waelchli R, Bollbuck B, Bruns C, Buhl T, Eder J, Felfel R, Hersperger R, Janser P, Revesz L, Zerwes HG, Schlapbach A: Design and preparation of 2-benzamido-pyrimidines as inhibitors of IKK. Bioorganic & Medicinal Chemistry Letters 2006, 16: 108–112. 10.1016/j.bmcl.2005.09.035View ArticleGoogle Scholar
- Castro AC, Dang LC, Soucy F, Grenier L, Mazdiyasni H, Hottelet M, Parent L, Pien C, Palombella V, Adams J: Novel IKK inhibitors: beta-carbolines. Bioorganic & Medicinal Chemistry Letters 2003, 13: 2419–2422. 10.1016/S0960-894X(03)00408-6View ArticleGoogle Scholar
- Belema M, Bunker A, Nguyen VN, Beaulieu F, Ouellet C, Qu YP, Zhang YH, Martel A, Burke JR, Mclntyre KW, et al.: Synthesis and structure-activity relationship of imidazo(1,2-a)thieno(3,2-e)pyrazines as IKK-beta inhibitors. Bioorganic & Medicinal Chemistry Letters 2007, 17: 4284–4289. 10.1016/j.bmcl.2007.05.031View ArticleGoogle Scholar
- Tripos Inc. SL, MO 63144, USA.
- Hindle SA, Rarey M, Buning C, Lengauer T: Flexible docking under pharmacophore type constraints. Journal of Computer-Aided Molecular Design 2002, 16: 129–149. 10.1023/A:1016399411208View ArticlePubMedGoogle Scholar
- Polgar T, Keseru GM: Virtual screening for beta-secretase (BACE1) inhibitors reveals the importance of protonation states at Asp32 and Asp228. Journal of Medicinal Chemistry 2005, 48: 3749–3755. 10.1021/jm049133bView ArticlePubMedGoogle Scholar
- Polgar T, Magyar C, Simon I, Keseru GM: Impact of ligand protonation on virtual screening against ss-secretase (BACEI). Journal of Chemical Information and Modeling 2007, 47: 2366–2373. 10.1021/ci700223pView ArticlePubMedGoogle Scholar
- Polgar T, Baki A, Szendrei GI, Keseru GM: Comparative virtual and experimental high-throughput screening for glycogen synthase kinase-3 beta inhibitors. Journal of Medicinal Chemistry 2005, 48: 7946–7959. 10.1021/jm050504dView ArticlePubMedGoogle Scholar
- Virag I, Polgar T, Keseru GM: Functional virtual screening of estrogen receptor alpha modulators by FlexX-Pharm. Journal of Molecular Structure-Theochem 2005, 725: 239–242. 10.1016/j.theochem.2005.02.083View ArticleGoogle Scholar
- Aronov AM, McClain B, Moody CS, Murcko MA: Kinase-likeness and kinase-privileged fragments: Toward virtual polypharmacology. Journal of Medicinal Chemistry 2008, 51: 1214–1222. 10.1021/jm701021bView ArticlePubMedGoogle Scholar
- Ghose AK, Herbertz T, Pippin DA, Salvino JM, Mallamo JP: Knowledge based prediction of ligand binding modes and rational inhibitor design for kinase drug discovery. Journal of Medicinal Chemistry 2008, 51: 5149–5171. 10.1021/jm800475yView ArticlePubMedGoogle Scholar
- Liao JJL: Molecular recognition of protein kinase binding pockets for design of potent and selective kinase inhibitors. Journal of Medicinal Chemistry 2007, 50: 409–424. 10.1021/jm0608107View ArticlePubMedGoogle Scholar
- Nagarajan S, Doddareddy MR, Choo H, Cho YS, Oh KS, Lee BH, Pae AN: IKK beta inhibitors identification part I: Homology model assisted structure based virtual screening. Bioorganic & Medicinal Chemistry 2009, 17: 2759–2766. 10.1016/j.bmc.2009.02.041View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.