Understanding and modulating cyclin-dependent kinase inhibitor specificity: molecular modeling and biochemical evaluation of pyrazolopyrimidinones as CDK2/cyclin A and CDK4/cyclin D1 inhibitors

Summary Cyclin-dependent kinases (CDKs) play a key role in regulating the cell cycle. The cyclins, their activating agents, and endogenous CDK inhibitors are frequently mutated in human cancers, making CDKs interesting targets for cancer chemotherapy. Our aim is the discovery of selective CDK4/cyclin D1 inhibitors. An ATP-competitive pyrazolopyrimidinone CDK inhibitor was identified by HTS and docked into a CDK4 homology model. The resulting binding model was consistent with available SAR and was validated by a subsequent CDK2/inhibitor crystal structure. An iterative cycle of chemistry and modeling led to a 70-fold improvement in potency. Small substituent changes resulted in large CDK4/CDK2 selectivity changes. The modeling revealed that selectivity is largely due to hydrogen-bonded interactions with only two kinase residues. This demonstrates that small differences between enzymes can efficiently be exploited in the design of selective inhibitors.     

吲哚衍生物类VEGFR-2酪氨酸激酶抑制剂的从头设计

以血管内皮生长因子受体-2(VEGFR-2)酪氨酸激酶的晶体结构为基础, 采用从头药物设计方法, 设计了一系列吲哚类化合物, 并用类药性和分子对接进行了筛选, 最后得到10个对接能量较低的化合物分子, 对具有最低结合能的化合物与VEGFR-2酪氨酸激酶的复合物进行了10 ns的分子动力学模拟, 并对其结合模式进行了分析. 这些化合物结构新颖, 可能作为抗肿瘤的先导化合物或候选药物. 本文结果为VEGFR-2酪氨酸激酶抑制剂的进一步改造、 设计及合成提供了理论基础, 并有助于开发高活性和高选择性的抗肿瘤药物.     

p53-MDM2结合抑制剂药效团模型的构建

采用Catalyst软件, 选择5类共24个p53-MDM2结合抑制剂作为训练集, 经计算机建模、构象优化, 由Catalyst系统构建出药效团模型, 并对药效团进行有效性分析, 结合已知的p53-MDM2结合抑制剂的结构信息, 筛选得到含有一个芳环中心、三个疏水中心和一个氢键受体的具有较好预测能力(Correl=0.941, Config=17.530, 吟cost=150.830)的药效团模型.     

The CK2 alpha/CK2 beta interface of human protein kinase CK2 harbors a binding pocket for small molecules

The Ser/Thr kinase CK2 (previously called casein kinase 2) is composed of two catalytic chains (CK2 alpha) attached to a dimer of noncatalytic subunits (CK2 beta). CK2 is involved in suppression of apoptosis, cell survival, and tumorigenesis. To investigate these activities and possibly affect them, selective CK2 inhibitors are required. An often-used CK2 inhibitor is 5,6-dichloro-1-beta-D-ribofuranosylbenzimidazole (DRB). In a complex structure with human CK2 alpha, DRB binds to the canonical ATP cleft, but additionally it occupies an allosteric site that can be alternatively filled by glycerol. Inhibition kinetic studies corroborate the dual binding mode of the inhibitor. Structural comparisons reveal a surprising conformational plasticity of human CK2 alpha around both DRB binding sites. After local rearrangement, the allosteric site serves as a CK2 beta interface. This opens the potential to construct molecules interfering with the CK2 alpha/CK2 beta interaction.     

Active site labeling of the gentamicin resistance enzyme AAC(6')-APH(2") by the lipid kinase inhibitor wortmannin

BACKGROUND: Aminoglycoside antibiotic resistance is largely the result of the production of enzymes that covalently modify the drugs including kinases (APHs) with structural and functional similarity to protein and lipid kinases. One of the most important aminoglycoside resistance enzymes is AAC(6')-APH(2"), a bifunctional enzyme with both aminoglycoside acetyltransferase and kinase activities. Knowledge of enzyme active site structure is important in deciphering the molecular mechanism of antibiotic resistance and here we explored active site labeling techniques to study AAC(6')-APH(2") structure and function. RESULTS: AAC(6')-APH(2") was irreversibly inactivated by wortmannin, a potent phosphatidylinositol 3-kinase inhibitor, through the covalent modification of a conserved lysine in the ATP binding pocket. 5'-[p-(Fluorosulfonyl)benzoyl]adenosine, an electrophilic ATP analogue and known inactivator of other APH enzymes such as APH(3')-IIIa, did not inactivate AAC(6')-APH(2"), and reciprocally, wortmannin did not inactivate APH(3')-IIIa. CONCLUSIONS: These distinct active site label sensitivities point to important differences in aminoglycoside kinase active site structures and suggest that design of broad range, ATP binding site-directed inhibitors against APHs will be difficult. Nonetheless, given the sensitivity of APH enzymes to both protein and lipid kinase inhibitors, potent lead inhibitors of this important resistance enzyme are likely to be found among the libraries of compounds directed against other pharmacologically important kinases.     

Molecular insight into the T798M gatekeeper mutation-caused acquired resistance to tyrosine kinase inhibitors in ErbB2-positive breast cancer

Human epidermal growth factor receptor 2 (ErbB2) is an attractive therapeutic target for metastatic breast cancer. The kinase has been clinically observed to harbor a gatekeeper mutation T798M in its active site, which causes acquired resistance to the first-line targeted breast cancer therapy with small-molecule tyrosine kinase inhibitors. Previously, several theories have been proposed to explain the molecular mechanism of gatekeeper mutation-caused drug resistance, such as blocking of inhibitor binding and increasing of ATP affinity. In the current study, the direct binding of three wild type-selective inhibitors (Lapatinib, AEE788 and TAK-285) and two wild type-sparing inhibitors (Staurosporine and Bosutinib) to the wild-type ErbB2 and its T798M mutant are investigated in detail by using rigorous computational analysis and binding affinity assay. Substitution of the polar threonine with a bulky methionine at residue 798 can impair and improve the direct binding affinity of wild type-selective and wild type-sparing inhibitors, respectively. Hindrance effect is responsible for the affinity decrease of wild type-selective inhibitors, while additional nonbonded interactions contribute to the affinity increase of wild type-sparing inhibitors, thus conferring selectivity to the inhibitors for mutant over wild type. The binding affinity of Staurosporine and Bosutinib to ErbB2 kinase domain is improved by 11.9-fold and 2.1-fold upon T798M mutation, respectively. Structural analysis reveals that a nonbonded network of S–π contact interactions (for Staurosporine) or an S-involving halogen bond (for Bosutinib) forms with the sulfide group of mutant Met798 residue.     

7-(2-Anilinopyrimidin-4-yl)-1-benzazepin-2-ones Designed by a “Cut and Glue” Strategy Are Dual Aurora A/VEGF-R Kinase Inhibitors

Although overexpression and hyperactivity of protein kinases are causative for a wide range of human cancers, protein kinase inhibitors currently approved as cancer drugs address only a limited number of these enzymes. To identify new chemotypes addressing alternative protein kinases, the basic structure of a known PLK1/VEGF-R2 inhibitor class was formally dissected and reassembled. The resulting 7-(2-anilinopyrimidin-4-yl)-1-benzazepin-2-ones were synthesized and proved to be dual inhibitors of Aurora A kinase and VEGF receptor kinases. Crystal structures of two representatives of the new chemotype in complex with Aurora A showed the ligand orientation in the ATP binding pocket and provided the basis for rational structural modifications. Congeners with attached sulfamide substituents retained Aurora A inhibitory activity. In vitro screening of two members of the new kinase inhibitor family against the cancer cell line panel of the National Cancer Institute (NCI) showed antiproliferative activity in the single-digit micromolar concentration range in the majority of the cell lines.     

人类胞外信号调节激酶1的同源模建及其抑制剂的对接与结构改造

通过同源模建和分子动力学模拟构建了人类胞外信号调节激酶1(hERK1)的三维结构,并利用profile-3D和procheck方法评估了模型的合理性.对所得的结构使用分子对接程序Affinity和CDOCKER进行了两种抑制剂的对接.结果显示这两种抑制剂与酶的结合方式相似,它们均与残基K36,Q87之间存在氢键作用,二者取代基的不同导致了抑制能力的差别.基于对接结果分析,对已知抑制剂进行结构改造,得到了一个理论上结合能力更强的抑制剂.它在保持与K36和Q87之间氢键的同时,又与残基D93,K96,S135形成了四条氢键,显著提高了与酶的相互作用.对接相互作用能显著下降,MM-PBSA结合自由能降为负值,这些均体现了抑制能力的提高.本工作对于针对该酶的抑制剂设计和相关疾病的新药开发具有理论指导价值.     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

Studying the Binding Modes of Novel 2-Aminopyridine Derivatives as Effective and Selective c-Met Kinase Type 1 Inhibitors Using Molecular Modeling Approaches

The mesenchymal epithelial cell transforming factor c-Met, encoded by c-Met proto-oncogene and known as a high-affinity receptor for Hepatocyte Growth Factor (HGF), is one of the receptor tyrosine kinases (RTKs) members. The HGF/c-Met signaling pathway has close correlation with tumor growth, invasion and metastasis. Thus, c-Met kinase has emerged as a prominent therapeutic target for cancer drug discovery. Recently a series of novel 2-aminopyridine derivatives targeting c-Met kinase with high biological activity were reported. In this study, 3D quantitative structure-activity relationship (QSAR), molecular docking and molecular dynamics simulations (MD) were employed to research the binding modes of these inhibitors.The results show that both the atom-based and docking-based CoMFA (Q² = 0.596, R² = 0.950 in atom-based model and Q² = 0.563, R² = 0.985 in docking-based model) and CoMSIA (Q² = 0.646, R² = 0.931 in atom-based model and Q² = 0.568, R² = 0.983 in docking-based model) models own satisfactory performance with good reliabilities and powerful external predictabilities. Molecular docking study suggests that Tyr1230 and Arg1208 might be the key residues, and electrostatic and hydrogen bond interactions were shown to be vital to the activity, concordance with QSAR analysis. Then MD simulation was performed to further explore the binding mode of the most potent inhibitor. The obtained results provide important references for further rational design of c-Met Kinase type I inhibitors.     

Identification of polyoxometalates as nanomolar noncompetitive inhibitors of protein kinase CK2

Protein kinase CK2 is a multifunctional kinase of medical importance that is dysregulated in many cancers. In this study, polyoxometalates were identified as original CK2 inhibitors. [P2Mo18O62](6-) has the most potent activity. It inhibits the kinase in the nanomolar range by targeting key structural elements located outside the ATP- and peptide substrate-binding sites. Several polyoxometalate derivatives exhibit strong inhibitory efficiency, with IC50 values < or = 10 nM. Furthermore, these inorganic compounds show a striking specificity for CK2 when tested in a panel of 29 kinases. Therefore, polyoxometalates are effective CK2 inhibitors in terms of both efficiency and selectivity and represent nonclassical kinase inhibitors that interact with CK2 in a unique way. This binding mode may provide an exploitable mechanism for developing potent drugs with desirable properties, such as enhanced selectivity relative to ATP-mimetic inhibitors.     

Structure-based 3D QSAR and design of novel acetylcholinesterase inhibitors

The paper describes the construction, validation and application of a structure-based 3D QSAR model of novel acetylcholinesterase (AChE) inhibitors. Initial use was made of four X-ray structures of AChE complexed with small, non-specific inhibitors to create a model of the binding of recently developed aminopyridazine derivatives. Combined automated and manual docking methods were applied to dock the co-crystallized inhibitors into the binding pocket. Validation of the modelling process was achieved by comparing the predicted enzyme-bound conformation with the known conformation in the X-ray structure. The successful prediction of the binding conformation of the known inhibitors gave confidence that we could use our model to evaluate the binding conformation of the aminopyridazine compounds. The alignment of 42 aminopyridazine compounds derived by the docking procedure was taken as the basis for a 3D QSAR analysis applying the GRID/GOLPE method. A model of high quality was obtained using the GRID water probe, as confirmed by the cross-validation method (q² _LOO=0.937, q² _{L50% O}=0.910). The validated model, together with the information obtained from the calculated AChE-inhibitor complexes, were considered for the design of novel compounds. Seven designed inhibitors which were synthesized and tested were shown to be highly active. After performing our modelling study the X-ray structure of AChE complexed with donepezil, an inhibitor structurally related to the developed aminopyirdazines, has been made available. The good agreement found between the predicted binding conformation of the aminopyridazines and the one observed for donepezil in the crystal structure further supports our developed model.     

In silico binding free energy predictability by using the linear interaction energy (LIE) method: bromobenzimidazole CK2 inhibitors as a case study

Protein kinase CK2 is essential for cell viability, and its control regards a broad series of cellular events such as gene expression, RNA, and protein synthesis. Evidence of its involvement in tumor development and viral replication indicates CK2 as a potential target of antineoplastic and antiviral drugs. In this study the Linear Interaction Energy (LIE) Method with the Surface Generalized Born (SGB) continuum solvation model was used to study several bromobenzimidazole CK2 inhibitors. This methodology, developed by Aqvist, finds a plausible compromise between accuracy and computational speed in evaluating binding free energy (DeltaGbind) values. In this study, two different free binding energy models, named "CK2scoreA" and "CK2scoreB", were developed using 22 inhibitors as the training set in a stepwise approach useful to appropriately select both the tautomeric form and the starting binding position of each inhibitor. Both models are statistically acceptable. Indeed, the better one is characterized by a correlation coefficient (r2) of 0.81, and the predictive accuracy was 0.65 kcal/mol. The corresponding validation, using an external test set of 16 analogs, showed a correlation coefficient (q2) of 0.68 and a prediction root-mean-square error of 0.78 kcal/mol. In this case, the LIE approach has been proved to be an efficient methodology to rationalize the difference of activity, the key interactions, and the different possible binding modes of this specific class of potent CK2 inhibitors.     

Small molecule inhibitors of the MDM2-p53 interaction discovered by ensemble-based receptor models

Five nonpeptide, small-molecule inhibitors of the human MDM2-p53 interaction are presented, and each inhibitor represents a new scaffold. The most potent compound exhibited a Ki of 110 +/- 30 nM. These compounds were identified using our multiple protein structure (MPS) method which incorporates protein flexibility into a receptor-based pharmacophore model that identifies appropriate hotspots of binding. Docking the inhibitors with an induced-fit docking protocol suggested that the inhibitors mimicked the three critical binding residues of p53 (Phe19, Trp23, and Leu26). Docking also predicted a new orientation of the scaffolds that more fully fills the binding cleft, enabling the inhibitors to take advantage of additional hydrogen-bonding possibilities not explored by other small molecule inhibitors. One inhibitor in particular was proposed to probe the hydrophobic core of the protein by taking advantage of the flexibility of the binding cleft floor. These results show that the MPS technique is a promising advance for structure-based drug discovery and that the method can truly explore broad chemical space efficiently in the quest to discover potent, small-molecule inhibitors of protein-protein interactions. Our MPS technique is one of very few ensemble-based techniques to be proven through experimental verification of the discovery of new inhibitors.     

Identification of novel nt-MGAM inhibitors for potential treatment of type 2 diabetes: Virtual screening,atom based 3D-QSAR model,docking analysis and ADME study

In this study, a virtual screening procedure was applied to identify new potential nt-MGAM inhibitors as a possible medication for type 2 diabetes. To this aim, a series of salacinol analogues were first investigated by docking analysis for their binding to the X-ray structure of the biological target nt-MGAM. Key interactions for ligand binding into the receptor active site were identified which shared common features to those found for other known inhibitors, which strengthen the results of this study. 3D QSAR model was then built and showed to be statistically significant and with a good predictive power for the training (R² = 0.99, SD = 0.17, F = 555.3 and N = 27) and test set (Q² = 0.81, Pearson(r) = 0.92, RMSE = 0.52, N = 08). The model was then used to virtually screen the ZINC database with the aim of identifying novel chemical scaffolds as potential nt-MGAM inhibitors. Further, in silico predicted ADME properties were investigated for the most promising molecules. The outcome of this investigation sheds light on the molecular characteristics of the binding of salacinol analogues to nt-MGAM enzyme and identifies new possible inhibitors which have the potential to be developed into drugs, thus significantly contributing to the design and optimization of therapeutic strategies against type 2 diabetes.     

Docking and 3D-QSAR Analysis on a Series of Pyridone-based EZH2 Inhibitors

Enhancer of Zeste homolog 2(EZH2) is closely correlated with malignant tumor and regarded as a promising target to treat B-cell lymphoma. In our research, the molecular docking and three-dimensional quantitative structure-activity relationships(3D-QSAR) studies were performed on a series of pyridone-based EZH2 compounds. Molecular docking allowed us to study the critical interactions at the binding site of EZH2 protein with inhibitors and identify the practical conformations of ligands in binding pocket. Moreover, the docking-based alignment was applied to derive the reliable 3D-QSAR models. Comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA) provided available ability of visualization. All the derived 3D-QSAR models were considered to be statistically significant with respect to the internal and external validation parameters. For the CoMFA model, q~2 = 0.649, r~2 = 0.961 and r~2 pred = 0.877. For the CoMSIA model, q~2 = 0.733, r~2 = 0.980 and r~2 pred = 0.848. With the above arguments, we extracted the correlation between the biological activity and structure. Based on the binding interaction and 3D contour maps, several new potential inhibitors with higher biological activity predicted were designed, which still awaited experimental validation. These theoretical conclusions could be helpful for further research and exploring potential EZH2 inhibitors.     

Exploration of a binding mode of benzothiazol-2-yl acetonitrile pyrimidine core based derivatives as potent c-Jun N-terminal kinase-3 inhibitors and 3D-QSAR analyses

C-Jun N-terminal kinase (JNK) is a therapeutic target for inhibitors which may provide clinical benefit in the pathogenesis of rheumatoid arthritis (RA) as well as in various apoptosis-related disorders. The benzothiazol-2-yl acetonitrile derivatives, recently reported by Pascale et al. (J. Med. Chem. 2005, 48, 4596-4607), are the first generation JNK inhibitors of this class. To understand inhibitory mechanisms and elucidate pharmacophoric properties of these derivatives molecular docking and 3D-QSAR studies were performed on a set of 44 compounds. Ligand Fit module of Cerius2 (4.9) was employed to locate the binding orientations of all the compounds within the JNK-3 ATP binding site. A good correlation (r2=0.810) between the calculated binding free energies (-PMF score) and the experimental inhibitory activities suggests that the identified binding conformations of these potential inhibitors are reliable. Based on the binding conformations, robust and highly predictive 3D-QSAR models were developed with conventional r2 0.886 and 0.802, full cross-validation r2 0.980 and 0.788, and predictive r2 0.965 and 0.968 for MFA and MSA, respectively. The interaction mode was demonstrated taking into consideration inhibitor conformation, hydrogen bonding, and electrostatic interaction. The 3D-QSAR model built in this study will provide clear guidelines for a novel inhibitor design based on the benzothiazole derivatives against JNK-3 for the treatment of inflammatory disorders.     

Investigation of indirubin derivatives: a combination of 3D-QSAR,molecular docking,and ADMET towards the design of new DRAK2 inhibitors

Kinase-related apoptosis-inducing kinase 2 (DRAK2) is a serine/threonine kinase and belongs to the death-associated protein kinase DPAK family, which is responsible for induction of apoptosis in many cell types. Thus, DRAK2 is regarded as a promising target for treatment of autoimmune diseases. To investigate the binding between DRAK2 and indirubin inhibitors and design potent inhibitors, a three-dimensional quantitative structure-activity relationship (3D-QSAR) and molecular docking were performed. Comparative Molecular Similarity Indices Analysis (CoMSIA) was developed using 33 molecules having pIC₅₀ ranging from 8.523 to 5.000 (IC₅₀ in nM). The best CoMSIA model gave a significant coefficient of determination (R²?=?0.93), as well as a (leave-one-out cross-validation coefficient Q² of 0.81. The predictive ability of this model was evaluated by external validation using a test set of eight compounds and yielded a predicted coefficient of determination R²_test of 0.94. The contour maps could provide structural features to improve inhibitory activity. Good consistency between contour maps and molecular docking strongly suggests that the molecular modeling is reliable. Based on these satisfactory results, we designed several new DRAK2 inhibitors and their inhibitory activities were predicted using different models, which are developed on different training and test sets. Additionally, these newly designed inhibitors showed promising results in the preliminary in silico ADMET evaluations compared to the best inhibitor from the studied dataset. This study could be useful in lead identification and optimization for early drug discovery of DRAK2 inhibitors.