Synthesis and application of functionally diverse 2,6,9-trisubstituted purine libraries as CDK inhibitors.

BACKGROUND: Purines constitute a structural class of protein ligands involved in mediating an astonishing array of metabolic processes and signal pathways in all living organisms. Synthesis of purine derivatives targeting specific purine-binding proteins in vivo could lead to versatile lead compounds for use as biological probes or drug candidates. RESULTS: We synthesized several libraries of 2,6, 9-trisubstituted purines using both solution- and solid-phase chemistry, and screened the compounds for inhibition of cyclin-dependent kinase (CDK) activity and human leukemic cell growth. Lead compounds were optimized by iterative synthesis based on structure-activity relationships (SARs), as well as analysis of several CDK-inhibitor cocrystal structures, to afford several interesting compounds including one of the most potent CDK inhibitors known to date. Unexpectedly, some compounds with similar CDK inhibitory activity arrested cellular proliferation at distinctly different phases of the cell cycle and another inhibitor directly induced apoptosis, bypassing cell-cycle arrest. Some of these compounds selectively inhibited growth of cells derived from specific tumors. CONCLUSIONS: 2,6,9-Trisubstituted purines have various and potent biological activities, despite high concentrations of competing endogenous purine ligands in living cells. Purine libraries constitute a versatile source of small molecules that affect distinct biochemical pathways mediating different cellular functions.     

A Combination of Spin Diffusion Methods for the Determination of Protein–Ligand Complex Structural Ensembles

Structure‐based drug design (SBDD) is a powerful and widely used approach to optimize affinity of drug candidates. With the recently introduced INPHARMA method, the binding mode of small molecules to their protein target can be characterized even if no spectroscopic information about the protein is known. Here, we show that the combination of the spin‐diffusion‐based NMR methods INPHARMA, trNOE, and STD results in an accurate scoring function for docking modes and therefore determination of protein–ligand complex structures. Applications are shown on the model system protein kinase A and the drug targets glycogen phosphorylase and soluble epoxide hydrolase (sEH). Multiplexing of several ligands improves the reliability of the scoring function further. The new score allows in the case of sEH detecting two binding modes of the ligand in its binding site, which was corroborated by X‐ray analysis.     

In silico study of porphyrin-anthraquinone hybrids as CDK2 inhibitor

Cyclin-Dependent Kinases (CDKs) are known to play crucial roles in controlling cell cycle progression of eukaryotic cell and inhibition of their activity has long been considered as potential strategy in anti-cancer drug research. In the present work, a series of porphyrin-anthraquinone hybrids bearing meso-substituents, i.e. either pyridine or pyrazole rings were designed and computationally evaluated for their Cyclin Dependent Kinase-2 (CDK2) inhibitory activity using molecular docking, molecular dynamics simulation, and binding free energy calculation. The molecular docking simulation revealed that all six porphyrin hybrids were able to bind to ATP-binding site of CDK2 and interacted with key residues constituted the active cavity of CDK2, while molecular dynamics simulation indicated that all porphyrins bound to CDK2 were stable for 6 ns. The binding free energies predicted by MM-PBSA method showed that most compounds exhibited higher affinity than that of native ligand (4-anilinoquinazoline, DTQ) and the affinity of mono-H₂PyP-AQ was about three times better than that of DTQ, indicating its potential to be advanced as a new CDK2 inhibitor.     

Cyclin-Dependent Kinase 4 and 6 Inhibitors in Cell Cycle Dysregulation for Breast Cancer Treatment

In cell development, the cell cycle is crucial, and the cycle progression’s main controllers are endogenous CDK inhibitors, cyclin-dependent kinases (CDKs), and cyclins. In response to the mitogenic signal, cyclin D is produced and retinoblastoma protein (Rb) is phosphorylated due to activated CDK4/CDK6. This causes various proteins required in the cell cycle progression to be generated. In addition, complexes of CDK1-cyclin A/B, CDK2-cyclin E/A, and CDK4/CDK6-cyclin D are required in each phase of this progression. Cell cycle dysregulation has the ability to lead to cancer. Based on its role in the cell cycle, CDK has become a natural target of anticancer therapy. Therefore, understanding the CDK structures and the complex formed with the drug, helps to foster the development of CDK inhibitors. This development starts from non-selective CDK inhibitors to selective CDK4/CDK6 inhibitors, and these have been applied in clinical cancer treatment. However, these inhibitors currently require further development for various hematologic malignancies and solid tumors, based on the results demonstrated. In drug development, the main strategy is primarily to prevent and asphyxiate drug resistance, thus a determination of specific biomarkers is required to increase the therapy’s effectiveness as well as patient selection suitability in order to avoid therapy failure. This review is expected to serve as a reference for early and advanced-stage researchers in designing new molecules or repurposing existing molecules as CDK4/CDK6 inhibitors to treat breast cancer.     

Improving binding mode and binding affinity predictions of docking by ligand-based search of protein conformations: evaluation in D3R grand challenge 2015

The growing number of protein–ligand complex structures, particularly the structures of proteins co-bound with different ligands, in the Protein Data Bank helps us tackle two major challenges in molecular docking studies: the protein flexibility and the scoring function. Here, we introduced a systematic strategy by using the information embedded in the known protein–ligand complex structures to improve both binding mode and binding affinity predictions. Specifically, a ligand similarity calculation method was employed to search a receptor structure with a bound ligand sharing high similarity with the query ligand for the docking use. The strategy was applied to the two datasets (HSP90 and MAP4K4) in recent D3R Grand Challenge 2015. In addition, for the HSP90 dataset, a system-specific scoring function (ITScore2_hsp90) was generated by recalibrating our statistical potential-based scoring function (ITScore2) using the known protein–ligand complex structures and the statistical mechanics-based iterative method. For the HSP90 dataset, better performances were achieved for both binding mode and binding affinity predictions comparing with the original ITScore2 and with ensemble docking. For the MAP4K4 dataset, although there were only eight known protein–ligand complex structures, our docking strategy achieved a comparable performance with ensemble docking. Our method for receptor conformational selection and iterative method for the development of system-specific statistical potential-based scoring functions can be easily applied to other protein targets that have a number of protein–ligand complex structures available to improve predictions on binding.     

Ligand-biased ensemble receptor docking (LigBEnD): a hybrid ligand/receptor structure-based approach

Ligand docking to flexible protein molecules can be efficiently carried out through ensemble docking to multiple protein conformations, either from experimental X-ray structures or from in silico simulations. The success of ensemble docking often requires the careful selection of complementary protein conformations, through docking and scoring of known co-crystallized ligands. False positives, in which a ligand in a wrong pose achieves a better docking score than that of native pose, arise as additional protein conformations are added. In the current study, we developed a new ligand-biased ensemble receptor docking method and composite scoring function which combine the use of ligand-based atomic property field (APF) method with receptor structure-based docking. This method helps us to correctly dock 30 out of 36 ligands presented by the D3R docking challenge. For the six mis-docked ligands, the cognate receptor structures prove to be too different from the 40 available experimental Pocketome conformations used for docking and could be identified only by receptor sampling beyond experimentally explored conformational subspace.     

Virtual screening for anti-HIV-1 RT and anti-HIV-1 PR inhibitors from the Thai medicinal plants database: a combined docking with neural networks approach

The virtual screening approach for docking small molecules into a known protein structure is a powerful tool for drug design. In this work, a combined docking and neural network approach, using a self-organizing map, has been developed and applied to screen anti-HIV-1 inhibitors for two targets, HIV-1 RT and HIV-1 PR, from active compounds available in the Thai Medicinal Plants Database. Based on nevirapine and calanolide A as reference structures in the HIV-1 RT binding site and XK-263 in the HIV-1 PR binding site, 2,684 compounds in the database were docked into the target enzymes. Self-organizing maps were then generated with respect to three types of pharmacophoric groups. The map of the reference structures were then superimposed on the feature maps of all screened compounds. Only the structures having similar features to the reference compounds were accepted. By using the SOMs, the number of candidates for HIV-1 RT was reduced to six and nine compounds consistent with nevirapine and calanolide A, respectively, as references. For the HIV-1 PR target, there are 135 screened compounds showed good agreement with the XK-263 feature map. These screened compounds will be further tested for their HIV-1 inhibitory affinities. The obtained results indicate that this combined method is clearly helpful to perform the successive screening and to reduce the analyzing step from AutoDock and scoring procedure.     

New approach to molecular docking and its application to virtual screening of chemical databases

This paper describes the validation of a molecular docking method and its application to virtual database screening. The code flexibly docks ligand molecules into rigid receptor structures using a tabu search methodology driven by an empirically derived function for estimating the binding affinity of a protein-ligand complex. The docking method has been tested on 70 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. The lowest energy geometry produced by the docking protocol is within 2.0 A root mean square of the experimental binding mode for 79% of the complexes. The method has been applied to the problem of virtual database screening to identify known ligands for thrombin, factor Xa, and the estrogen receptor. A database of 10,000 randomly chosen "druglike" molecules has been docked into the three receptor structures. In each case known receptor ligands were included in the study. The results showed good separation between the predicted binding affinities of the known ligand set and the database subset.     

Validation of an empirical RNA-ligand scoring function for fast flexible docking using RiboDock®

We report the design and validation of a fast empirical function for scoring RNA-ligand interactions, and describe its implementation within RiboDock, a virtual screening system for automated flexible docking. Building on well-known protein-ligand scoring function foundations, features were added to describe the interactions of common RNA-binding functional groups that were not handled adequately by conventional terms, to disfavour non-complementary polar contacts, and to control non-specific charged interactions. The results of validation experiments against known structures of RNA-ligand complexes compare favourably with previously reported methods. Binding modes were well predicted in most cases and good discrimination was achieved between native and non-native ligands for each binding site, and between native and non-native binding sites for each ligand. Further evidence of the ability of the method to identify true RNA binders is provided by compound selection ('enrichment factor') experiments based around a series of HIV-1 TAR RNA-binding ligands. Significant enrichment in true binders was achieved amongst high scoring docking hits, even when selection was from a library of structurally related, positively charged molecules. Coupled with a semi-automated cavity detection algorithm for identification of putative ligand binding sites, also described here, the method is suitable for the screening of very large databases of molecules against RNA and RNA-protein interfaces, such as those presented by the bacterial ribosome.     

Statistical potential for modeling and ranking of protein-ligand interactions

Applications in structural biology and medicinal chemistry require protein-ligand scoring functions for two distinct tasks: (i) ranking different poses of a small molecule in a protein binding site and (ii) ranking different small molecules by their complementarity to a protein site. Using probability theory, we developed two atomic distance-dependent statistical scoring functions: PoseScore was optimized for recognizing native binding geometries of ligands from other poses and RankScore was optimized for distinguishing ligands from nonbinding molecules. Both scores are based on a set of 8,885 crystallographic structures of protein-ligand complexes but differ in the values of three key parameters. Factors influencing the accuracy of scoring were investigated, including the maximal atomic distance and non-native ligand geometries used for scoring, as well as the use of protein models instead of crystallographic structures for training and testing the scoring function. For the test set of 19 targets, RankScore improved the ligand enrichment (logAUC) and early enrichment (EF(1)) scores computed by DOCK 3.6 for 13 and 14 targets, respectively. In addition, RankScore performed better at rescoring than each of seven other scoring functions tested. Accepting both the crystal structure and decoy geometries with all-atom root-mean-square errors of up to 2 ? from the crystal structure as correct binding poses, PoseScore gave the best score to a correct binding pose among 100 decoys for 88% of all cases in a benchmark set containing 100 protein-ligand complexes. PoseScore accuracy is comparable to that of DrugScore(CSD) and ITScore/SE and superior to 12 other tested scoring functions. Therefore, RankScore can facilitate ligand discovery, by ranking complexes of the target with different small molecules; PoseScore can be used for protein-ligand complex structure prediction, by ranking different conformations of a given protein-ligand pair. The statistical potentials are available through the Integrative Modeling Platform (IMP) software package (http://salilab.org/imp) and the LigScore Web server (http://salilab.org/ligscore/).     

Nonfitting protein–ligand interaction scoring function based on first‐principles theoretical chemistry methods: Development and application on kinase inhibitors

Targeted therapy is currently a hot topic in the fields of cancer research and drug design. An important requirement for this approach is the development of potent and selective inhibitors for the identified target protein. However, current ways to estimate inhibitor efficacy rely on empirical protein–ligand interaction scoring functions which, suffering from their heavy parameterizations, often lead to a low accuracy. In this work, we develop a nonfitting scoring function, which consists of three terms: (1) gas‐phase protein‐ligand binding enthalpy obtained by the eXtended ONIOM hybrid method based on an integration of density functional theory (DFT) methods (XYG3 and ωB97X‐D) and the semiempirical PM6 method, (2) solvation free energy based on DFT‐SMD solvation model, and (3) entropy effect estimated by using DFT frequency analysis. The new scoring function is tested on a cyclin‐dependent kinase 2 (CDK2) inhibitor database including 76 CDK2 protein inhibitors and a p21‐activated kinase 1 (PAK1) inhibitor database including 20 organometallic PAK1 protein inhibitors. From the results, good correlations are found between the calculated scores and the experimental inhibitor efficacies with the square of correlation coefficient R² of 0.76–0.88. This suggests a good predictive power of this scoring function. To the best of our knowledge, this is the first high level theory‐based nonfitting scoring function with such a good level of performance. This scoring function is recommended to be used in the final screening of lead structure derivatives. © 2013 Wiley Periodicals, Inc.     

Ligand-Based Pharmacophore Screening Strategy: a Pragmatic Approach for Targeting HER Proteins

Targeting ErbB family of receptors is an important therapeutic option, because of its essential role in the broad spectrum of human cancers, including non-small cell lung cancer (NSCLC). Therefore, in the present work, considerable effort has been made to develop an inhibitor against HER family proteins, by combining the use of pharmacophore modelling, docking scoring functions, and ADME property analysis. Initially, a five-point pharmacophore model was developed using known HER family inhibitors. The generated model was then used as a query to screen a total of 468,880 compounds of three databases namely ZINC, ASINEX, and DrugBank. Subsequently, docking analysis was carried out to obtain hit molecules that could inhibit the HER receptors. Further, analysis of GLIDE scores and ADME properties resulted in one hit namely BAS01025917 with higher glide scores, increased CNS involvement, and good pharmaceutically relevant properties than reference ligand, afatinib. Furthermore, the inhibitory activity of the lead compounds was validated by performing molecular dynamic simulations. Of note, BAS01025917 was found to possess scaffolds with a broad spectrum of antitumor activity. We believe that this novel hit molecule can be further exploited for the development of a pan-HER inhibitor with low toxicity and greater potential.     

In silico fragment-mapping method: a new tool for fragment-based/structure-based drug discovery

Here, we propose an in silico fragment-mapping method as a potential tool for fragment-based/structure-based drug discovery (FBDD/SBDD). For this method, we created a database named Canonical Subsite–Fragment DataBase (CSFDB) and developed a knowledge-based fragment-mapping program, Fsubsite. CSFDB consists of various pairs of subsite–fragments derived from X-ray crystal structures of known protein–ligand complexes. Using three-dimensional similarity-matching between subsites on one protein and another, Fsubsite compares the surface of a target protein with all subsites in CSFDB. When a local topography similar to the subsite is found on the surface, Fsubsite places a fragment combined with the subsite in CSFDB on the target protein. For validation purposes, we applied the method to the apo-structure of cyclin-dependent kinase 2 (CDK2) and identified four compounds containing three mapped fragments that existed in the list of known inhibitors of CDK2. Next, the utility of our fragment-mapping method for fragment-growing was examined on the complex structure of tRNA-guanine transglycosylase with a small ligand. Fsubsite mapped appropriate fragments on the same position as the binding ligand or in the vicinity of the ligand. Finally, a 3D-pharmacophore model was constructed from the fragments mapped on the apo-structure of heat shock protein 90-α (HSP90α). Then, 3D pharmacophore-based virtual screening was carried out using a commercially available compound database. The resultant hit compounds were very similar to a known ligand of HSP90α. As a result of these findings, this in silico fragment-mapping method seems to be a useful tool for computational FBDD and SBDD.     

Molecular dynamics, density functional, ADMET predictions, virtual screening, and molecular interaction field studies for identification and evaluation of novel potential CDK2 inhibitors in cancer therapy

In this work, we have used molecular dynamics, density functional theory, virtual screening, ADMET predictions, and molecular interaction field studies to design and propose eight novel potential inhibitors of CDK2. The eight molecules proposed showed interesting structural characteristics that are required for inhibiting the CDK2 activity and show potential as drug candidates for the treatment of cancer. The parameters related to the Rule of Five were calculated, and only one of the molecules violated more than one parameter. One of the proposals and one of the drug-like compounds selected by virtual screening indicated to be promising candidates for CDK2-based cancer therapy.     

Fragment-based de novo ligand design by multiobjective evolutionary optimization

GANDI (Genetic Algorithm-based de Novo Design of Inhibitors) is a computational tool for automatic fragment-based design of molecules within a protein binding site of known structure. A genetic algorithm and a tabu search act in concert to join predocked fragments with a user-supplied list of fragments. A novel feature of GANDI is the simultaneous optimization of force field energy and a term enforcing 2D-similarity to known inhibitor(s) or 3D-overlap to known binding mode(s). Scaffold hopping can be promoted by tuning the relative weights of these terms. The performance of GANDI is tested on cyclin-dependent kinase 2 (CDK2) using a library of about 14 000 fragments and the binding mode of a known oxindole inhibitor to bias the design. Top ranking GANDI molecules are involved in one to three hydrogen bonds with the backbone polar groups in the hinge region of CDK2, an interaction pattern observed in potent kinase inhibitors. Notably, a GANDI molecule with very favorable predicted binding affinity shares a 2-N-phenyl-1,3-thiazole-2,4-diamine moiety with a known nanomolar inhibitor of CDK2. Importantly, molecules with a favorable GANDI score are synthetic accessible. In fact, eight of the 1809 molecules designed by GANDI for CDK2 are found in the ZINC database of commercially available compounds which also contains about 600 compounds with identical scaffolds as those in the top ranking GANDI molecules.     

Prediction of multiple binding modes of the CDK2 inhibitors, anilinopyrazoles, using the automated docking programs GOLD, FlexX, and LigandFit: an evaluation of performance

Anilinopyrazoles as CDK2 inhibitors can adopt multiple binding modes depending on the substituents at the 5-position of the pyrazole ring, based on CDK2/cyclin A crystallographic studies. Three commercially available docking programs, FlexX, GOLD, and LigandFit, were tested with 63 anilinopyrazole analogues in an attempt to reproduce the binding modes observed in the crystal structures. Each docking program gave different ligand conformations depending on the scoring or energy functions used. FlexX/drugscore, GOLD/chemscore, and LigandFit/plp were the best combinations of each docking program in reproducing the ligand conformations observed in the crystal structures. The 63 analogues were divided into two groups, type-A and type-B, depending on the substituent at the 5-position of the pyrazole ring. Although an alternate binding mode, observed in a crystal structure of one type-B compound, could not be reproduced with any of the above docking/scoring combinations, GOLD, with a template constraint based on the crystal structure coordinates, was able to reproduce the pose. As for type-A compounds, all docking conditions yielded similar poses to those observed in crystal structures. When predicting activities by scoring programs, the combination of docking with LigandFit/plp and scoring with LIGSCORE1_CFF gave the best correlation coefficient (r=0.60) between experimental pIC50 values and top-ranked rescores of 30 poses of each compound. With regard to type-A compounds, the correlation was 0.69. However, when 11 compounds, whose top-ranked rescored poses did not demonstrate the correct binding modes in reference to the crystal structure, were removed, the correlation rose to 0.75. Consequently, predicting activity on the basis of correct binding modes was found to be reliable.     

Large scale free energy calculations for blind predictions of protein–ligand binding: the D3R Grand Challenge 2015

We describe binding free energy calculations in the D3R Grand Challenge 2015 for blind prediction of the binding affinities of 180 ligands to Hsp90. The present D3R challenge was built around experimental datasets involving Heat shock protein (Hsp) 90, an ATP-dependent molecular chaperone which is an important anticancer drug target. The Hsp90 ATP binding site is known to be a challenging target for accurate calculations of ligand binding affinities because of the ligand-dependent conformational changes in the binding site, the presence of ordered waters and the broad chemical diversity of ligands that can bind at this site. Our primary focus here is to distinguish binders from nonbinders. Large scale absolute binding free energy calculations that cover over 3000 protein–ligand complexes were performed using the BEDAM method starting from docked structures generated by Glide docking. Although the ligand dataset in this study resembles an intermediate to late stage lead optimization project while the BEDAM method is mainly developed for early stage virtual screening of hit molecules, the BEDAM binding free energy scoring has resulted in a moderate enrichment of ligand screening against this challenging drug target. Results show that, using a statistical mechanics based free energy method like BEDAM starting from docked poses offers better enrichment than classical docking scoring functions and rescoring methods like Prime MM-GBSA for the Hsp90 data set in this blind challenge. Importantly, among the three methods tested here, only the mean value of the BEDAM binding free energy scores is able to separate the large group of binders from the small group of nonbinders with a gap of 2.4 kcal/mol. None of the three methods that we have tested provided accurate ranking of the affinities of the 147 active compounds. We discuss the possible sources of errors in the binding free energy calculations. The study suggests that BEDAM can be used strategically to discriminate binders from nonbinders in virtual screening and to more accurately predict the ligand binding modes prior to the more computationally expensive FEP calculations of binding affinity.     

New scoring functions for virtual screening from molecular dynamics simulations with a quantum-refined force-field (QRFF-MD). Application to cyclin-dependent kinase 2

A recently introduced new methodology based on ultrashort (50-100 ps) molecular dynamics simulations with a quantum-refined force-field (QRFF-MD) is here evaluated in its ability both to predict protein-ligand binding affinities and to discriminate active compounds from inactive ones. Physically based scoring functions are derived from this approach, and their performance is compared to that of several standard knowledge-based scoring functions. About 40 inhibitors of cyclin-dependent kinase 2 (CDK2) representing a broad chemical diversity were considered. The QRFF-MD method achieves a correlation coefficient, R(2), of 0.55, which is significantly better than that obtained by a number of traditional approaches in virtual screening but only slightly better than that obtained by consensus scoring (R(2) = 0.50). Compounds from the Available Chemical Directory, along with the known active compounds, were docked into the ATP binding site of CDK2 using the program Glide, and the 650 ligands from the top scored poses were considered for a QRFF-MD analysis. Combined with structural information extracted from the simulations, the QRFF-MD methodology results in similar enrichment of known actives compared to consensus scoring. Moreover, a new scoring function is introduced that combines a QRFF-MD based scoring function with consensus scoring, which results in substantial improvement on the enrichment profile.     

A machine learning-based method to improve docking scoring functions and its application to drug repurposing

Docking scoring functions are notoriously weak predictors of binding affinity. They typically assign a common set of weights to the individual energy terms that contribute to the overall energy score; however, these weights should be gene family dependent. In addition, they incorrectly assume that individual interactions contribute toward the total binding affinity in an additive manner. In reality, noncovalent interactions often depend on one another in a nonlinear manner. In this paper, we show how the use of support vector machines (SVMs), trained by associating sets of individual energy terms retrieved from molecular docking with the known binding affinity of each compound from high-throughput screening experiments, can be used to improve the correlation between known binding affinities and those predicted by the docking program eHiTS. We construct two prediction models: a regression model trained using IC(50) values from BindingDB, and a classification model trained using active and decoy compounds from the Directory of Useful Decoys (DUD). Moreover, to address the issue of overrepresentation of negative data in high-throughput screening data sets, we have designed a multiple-planar SVM training procedure for the classification model. The increased performance that both SVMs give when compared with the original eHiTS scoring function highlights the potential for using nonlinear methods when deriving overall energy scores from their individual components. We apply the above methodology to train a new scoring function for direct inhibitors of Mycobacterium tuberculosis (M.tb) InhA. By combining ligand binding site comparison with the new scoring function, we propose that phosphodiesterase inhibitors can potentially be repurposed to target M.tb InhA. Our methodology may be applied to other gene families for which target structures and activity data are available, as demonstrated in the work presented here.