Identification of a buried pocket for potent and selective inhibition of Chk1: prediction and verification

Inhibition of the Chk1 kinase by small molecules binding to its active site is a strategy of great therapeutic interest for oncology. We report how computational modelling predicted the binding mode of ligands of special interest to the Chk1 ATP site, for representatives of an indazole series and debromohymenialdisine. These binding modes were subsequently confirmed by X-ray crystallography. The binding mode of a potent indazole derivative involves non-conventional C-H...O and N-H...pi-aromatic interactions with the protein. These interactions are formed in a buried pocket at the periphery of the ATP-binding site, the importance of which has previously been overlooked for ligand design against Chk1. It is demonstrated that filling this pocket can confer ligands with dramatically enhanced affinity for Chk1. Structural arguments in conjunction with assay data explain why targeting this pocket is also advantageous for selective binding to Chk1. Structural overlays of known inhibitors complexed with Chk1 show that only the indazole series utilizes the pocket of interest. Therefore, the analysis presented here should prove helpful in guiding future structure-based ligand design efforts against Chk1.     

Structural basis for Chk1 inhibition by UCN-01

Chk1 is a serine-threonine kinase that plays an important role in the DNA damage response, including G(2)/M cell cycle control. UCN-01 (7-hydroxystaurosporine), currently in clinical trials, has recently been shown to be a potent Chk1 inhibitor that abrogates the G(2)/M checkpoint induced by DNA-damaging agents. To understand the structural basis of Chk1 inhibition by UCN-01, we determined the crystal structure of the Chk1 kinase domain in complex with UCN-01. Chk1 structures with staurosporine and its analog SB-218078 were also determined. All three compounds bind in the ATP-binding pocket of Chk1, producing only slight changes in the protein conformation. Selectivity of UCN-01 toward Chk1 over cyclin-dependent kinases can be explained by the presence of a hydroxyl group in the lactam moiety interacting with the ATP-binding pocket. Hydrophobic interactions and hydrogen-bonding interactions were observed in the structures between UCN-01 and the Chk1 kinase domain. The high structural complementarity of these interactions is consistent with the potency and selectivity of UCN-01.     

X-ray structures of checkpoint kinase 2 in complex with inhibitors that target its gatekeeper-dependent hydrophobic pocket

The serine/threonine checkpoint kinase 2 (Chk2) is an attractive molecular target for the development of small molecule inhibitors to treat cancer. Here, we report the rational design of Chk2 inhibitors that target the gatekeeper-dependent hydrophobic pocket located behind the adenine-binding region of the ATP-binding site. These compounds exhibit IC(50) values in the low nanomolar range and are highly selective for Chk2 over Chk1. X-ray crystallography was used to determine the structures of the inhibitors in complex with the catalytic kinase domain of Chk2 to verify their modes of binding.     

Kinetics, in silico docking, molecular dynamics, and MM-GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP-binding site inhibitors: the role of water molecules examined

With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC?? > 50 μM), indirubin-3'-oxime (IC?? = 144 nM), KT5720 (K(i) = 18.4 nM) and staurosporine (K(i) = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP-binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MMGBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc-staurosporine complex has the greatest number of receptor-ligand hydrogen bonds, while for the indirubin-3'-oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor-enzyme contacts. The MM-GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1-water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schr?dinger induced-fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor-ligand contacts.     

花青素主要成分与HER-2激酶区的分子对接

用分子对接方法预测天然植物化学物质与受体蛋白的相互作用位点并探究作用机制。利用MVD(Molecular Virtual Docker 5.5)软件,以HER-2激酶区为受体模板建立活性位点,与12种花青素成分进行分子对接。结果表明12种化合物均能在同一活性腔中与HER-2激酶区对接(MolDock Score:苷元–105 kJ/mol,单葡糖苷–130 kJ/mol),主要作用力是疏水作用和氢键;该活性腔也是ATP与HER-2激酶区的结合(MolDock Score=–161 kJ/mol)位点,花青素的结合可能会干扰ATP与HER-2之间氢键的形成。提示花青素可能以竞争性结合方式阻碍ATP与HER-2的结合,抑制HER-2磷酸化激活及下游信号通路的激活,从而发挥抑癌活性。     

Ligand based virtual screening and biological evaluation of inhibitors of chorismate mutase (Rv1885c) from Mycobacterium tuberculosis H37Rv

We have identified new lead candidates that possess inhibitory activity against Mycobacterium tuberculosis H37Rv chorismate mutase by a ligand-based virtual screening optimized for lead evaluation in combination with in vitro enzymatic assay. The initial virtual screening using a ligand-based pharmacophore model identified 95 compounds from an in-house small molecule database of 15,452 compounds. The obtained hits were further evaluated by molecular docking and 15 compounds were short listed based on docking scores and the other scoring functions and subjected to biological assay. Chorismate mutase activity assays identified four compounds as inhibitors of M. tuberculosis chorismate mutase (MtCM) with low K(i) values. The structural models for these ligands in the chorismate mutase binding site will facilitate medicinal chemistry efforts for lead optimization against this protein.     

In silico investigation of potential SRC kinase ligands from traditional Chinese medicine

Src kinase is an attractive target for drug development based on its established relationship with cancer and possible link to hypertension. The suitability of traditional Chinese medicine (TCM) compounds as potential drug ligands for further biological evaluation was investigated using structure-based, ligand-based, and molecular dynamics (MD) analysis. Isopraeroside IV, 9alpha-hydroxyfraxinellone-9-O-beta-D-glucoside (9HFG) and aurantiamide were the top three TCM candidates identified from docking. Hydrogen bonds and hydrophobic interactions were the primary forces governing docking stability. Their stability with Src kinase under a dynamic state was further validated through MD and torsion angle analysis. Complexes formed by TCM candidates have lower total energy estimates than the control Sacaratinib. Four quantitative-structural activity relationship (QSAR) in silico verifications consistently suggested that the TCM candidates have bioactive properties. Docking conformations of 9HFG and aurantiamide in the Src kinase ATP binding site suggest potential inhibitor-like characteristics, including competitive binding at the ATP binding site (Lys295) and stabilization of the catalytic cleft integrity. The TCM candidates have significantly lower ligand internal energies and are estimated to form more stable complexes with Src kinase than Saracatinib. Structure-based and ligand-based analysis support the drug-like potential of 9HFG and aurantiamide and binding mechanisms reveal the tendency of these two candidates to compete for the ATP binding site.     

Development of thioquinazolinones,allosteric Chk1 kinase inhibitors

A high throughput screening campaign was designed to identify allosteric inhibitors of Chk1 kinase by testing compounds at high concentration. Activity was then observed at K_m for ATP and at near-physiological concentrations of ATP. This strategy led to the discovery of a non-ATP competitive thioquinazolinone series which was optimized for potency and stability. An X-ray crystal structure for the complex of our best inhibitor bound to Chk1 was solved, indicating that it binds to an allosteric site ～13 Å from the ATP binding site. Preliminary data is presented for several of these compounds.     

A Targeted Library Screen Reveals a New Inhibitor Scaffold for Protein Kinase D

Protein kinase D (PKD) has emerged as a potential therapeutic target in multiple pathological conditions, including cancer and heart diseases. Potent and selective small molecule inhibitors of PKD are valuable for dissecting PKD-mediated cellular signaling pathways and for therapeutic application. In this study, we evaluated a targeted library of 235 small organic kinase inhibitors for PKD1 inhibitory activity at a single concentration. Twenty-eight PKD inhibitory chemotypes were identified and six exhibited excellent PKD1 selectivity. Five of the six lead structures share a common scaffold, with compound 139 being the most potent and selective for PKD vs PKC and CAMK. Compound 139 was an ATP-competitive PKD1 inhibitor with a low double-digit nanomolar potency and was also cell-active. Kinase profiling analysis identified this class of small molecules as pan-PKD inhibitors, confirmed their selectivity again PKC and CAMK, and demonstrated an overall favorable selectivity profile that could be further enhanced through structural modification. Furthermore, using a PKD homology model based on similar protein kinase structures, docking modes for compound 139 were explored and compared to literature examples of PKD inhibition. Modeling of these compounds at the ATP-binding site of PKD was used to rationalize its high potency and provide the foundation for future further optimization. Accordingly, using biochemical screening of a small number of privileged scaffolds and computational modeling, we have identified a new core structure for highly potent PKD inhibition with promising selectivity against closely related kinases. These lead structures represent an excellent starting point for the further optimization and the design of selective and therapeutically effective small molecule inhibitors of PKD.     

Benzimidazole inhibitors of the protein kinase CHK2: Clarification of the binding mode by flexible side chain docking and protein–ligand crystallography

Two closely related binding modes have previously been proposed for the ATP-competitive benzimidazole class of checkpoint kinase 2 (CHK2) inhibitors; however, neither binding mode is entirely consistent with the reported SAR. Unconstrained rigid docking of benzimidazole ligands into representative CHK2 protein crystal structures reveals an alternative binding mode involving a water-mediated interaction with the hinge region; docking which incorporates protein side chain flexibility for selected residues in the ATP binding site resulted in a refinement of the water-mediated hinge binding mode that is consistent with observed SAR. The flexible docking results are in good agreement with the crystal structures of four exemplar benzimidazole ligands bound to CHK2 which unambiguously confirmed the binding mode of these inhibitors, including the water-mediated interaction with the hinge region, and which is significantly different from binding modes previously postulated in the literature.     

Discovery of d-amino acid oxidase inhibitors based on virtual screening against the lid-open enzyme conformation

d-Amino acid oxidase (DAAO) inhibitors are typically small polar compounds with often suboptimal pharmacokinetic properties. Features of the native binding site limit the operational freedom of further medicinal chemistry efforts. We therefore initiated a structure based virtual screening campaign based on the X-ray structures of DAAO complexes where larger ligands shifted the loop (lid opening) covering the native binding site. The virtual screening of our in-house collection followed by the in vitro test of the best ranked compounds led to the identification of a new scaffold with micromolar IC₅₀. Subsequent SAR explorations enabled us to identify submicromolar inhibitors. Docking studies supported by in vitro activity measurements suggest that compounds bind to the active site with a salt-bridge characteristic to DAAO inhibitor binding. In addition, displacement of and interaction with the loop covering the active site contributes significantly to the activity of the most potent compounds.     

Synthesis and evaluation of heteroaryl substituted diazaspirocycles as scaffolds to probe the ATP-binding site of protein kinases

With the success of protein kinase inhibitors as drugs to target cancer, there is a continued need for new kinase inhibitor scaffolds. We have investigated the synthesis and kinase inhibition of new heteroaryl-substituted diazaspirocyclic compounds that mimic ATP. Versatile syntheses of substituted diazaspirocycles through ring-closing metathesis were demonstrated. Diazaspirocycles directly linked to heteroaromatic hinge binder groups provided ligand efficient inhibitors of multiple kinases, suitable as starting points for further optimization. The binding modes of representative diazaspirocyclic motifs were confirmed by protein crystallography. Selectivity profiles were influenced by the hinge binder group and the interactions of basic nitrogen atoms in the scaffold with acidic side-chains of residues in the ATP pocket. The introduction of more complex substitution to the diazaspirocycles increased potency and varied the selectivity profiles of these initial hits through engagement of the P-loop and changes to the spirocycle conformation, demonstrating the potential of these core scaffolds for future application to kinase inhibitor discovery.     

Escape from adamantane: Scaffold optimization of novel P2X7 antagonists featuring complex polycycles

The adamantane scaffold, despite being widely used in medicinal chemistry, is not devoid of problems. In recent years we have developed new polycyclic scaffolds as surrogates of the adamantane group with encouraging results in multiple targets. As an adamantane scaffold is a common structural feature in several P2X7 receptor antagonists, herein we report the synthesis and pharmacological evaluation of multiple replacement options of adamantane that maintain a good activity profile. Molecular modeling studies support the binding of the compounds to a site close to the central pore, rather than to the ATP-binding site and shed light on the structural requirements for novel P2X7 antagonists.     

Catalytic transitions in the human MDR1 P-glycoprotein drug binding sites

Multidrug resistance proteins that belong to the ATP-binding cassette family like the human P-glycoprotein (ABCB1 or Pgp) are responsible for many failed cancer and antiviral chemotherapies because these membrane transporters remove the chemotherapeutics from the targeted cells. Understanding the details of the catalytic mechanism of Pgp is therefore critical to the development of inhibitors that might overcome these resistances. In this work, targeted molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp. Crystal structures of homologues in four different conformations were used as intermediate targets in the dynamics simulations. Transitions from conformations that were wide open to the cytoplasm to transition state conformations that were wide open to the extracellular space were studied. Twenty-six nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses. Coupled movement of nucleotide binding domains (NBDs) and transmembrane domains (TMDs) that form the drug binding cavities were observed. Pronounced twisting of the NBDs as they approached each other as well as the quantification of a dramatic opening of the TMDs to the extracellular space as the ATP hydrolysis transition state was reached were observed. Docking interactions of 21 known transport ligands or inhibitors were analyzed with each of the 26 transitional structures. Many of the docking results obtained here were validated by previously published biochemical determinations. As the ATP hydrolysis transition state was approached, drug docking in the extracellular half of the transmembrane domains seemed to be destabilized as transport ligand exit gates opened to the extracellular space.     

ATP and AMP Mutually Influence Their Interaction with the ATP-binding Cassette (ABC) Adenylate Kinase Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) at Separate Binding Sites

Cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel in the ATP-binding cassette (ABC) transporter protein family. In the presence of ATP and physiologically relevant concentrations of AMP, CFTR exhibits adenylate kinase activity (ATP + AMP ⇆ 2 ADP). Previous studies suggested that the interaction of nucleotide triphosphate with CFTR at ATP-binding site 2 is required for this activity. Two other ABC proteins, Rad50 and a structural maintenance of chromosome protein, also have adenylate kinase activity. All three ABC adenylate kinases bind and hydrolyze ATP in the absence of other nucleotides. However, little is known about how an ABC adenylate kinase interacts with ATP and AMP when both are present. Based on data from non-ABC adenylate kinases, we hypothesized that ATP and AMP mutually influence their interaction with CFTR at separate binding sites. We further hypothesized that only one of the two CFTR ATP-binding sites is involved in the adenylate kinase reaction. We found that 8-azidoadenosine 5′-triphosphate (8-N₃-ATP) and 8-azidoadenosine 5′-monophosphate (8-N₃-AMP) photolabeled separate sites in CFTR. Labeling of the AMP-binding site with 8-N₃-AMP required the presence of ATP. Conversely, AMP enhanced photolabeling with 8-N₃-ATP at ATP-binding site 2. The adenylate kinase active center probe P¹,P⁵-di(adenosine-5′) pentaphosphate interacted simultaneously with an AMP-binding site and ATP-binding site 2. These results show that ATP and AMP interact with separate binding sites but mutually influence their interaction with the ABC adenylate kinase CFTR. They further indicate that the active center of the adenylate kinase comprises ATP-binding site 2.     

Identification of novel, selective and potent Chk2 inhibitors

A series of isothiazole carboxamidine compounds were synthesized and discovered as novel and selective inhibitors for Chk2. They are not active against the related Chk1 kinase. The structure-activity relationship studies were performed on the scaffold, and enzymatic kinetic analysis showed they are simple ATP competitive inhibitors with K(i) values as low as 11 nM for Chk2. Computer modeling studies were employed to comprehend the mechanism of action and SAR of these compounds.     

The structure of gelsolin bound to ATP

Calcium activation of the actin-modifying properties of gelsolin is sensitive to ATP. Here, we show that soaking calcium-free gelsolin crystals in ATP-containing media results in ATP occupying a site that spans the two pseudosymmetrical halves of the protein. ATP binding involves numerous polar and hydrophobic contacts and is identical for the two copies of gelsolin related by non-crystallographic symmetry within the crystal. The gamma-phosphate of ATP participates in several charge-charge interactions consistent with the preference of gelsolin for ATP, as a binding partner, over ADP. In addition, disruption of the ATP-binding site through Ca2+ activation of gelsolin reveals why ATP binds more tightly to the inactive molecule, and suggests how the binding of ATP may modulate the sensitivity of gelsolin to calcium ions. Similarities between the ATP and PIP2 interactions with the C-terminal half of gelsolin are evident from their overlapping binding sites and in that both molecules bind more tightly in the absence of calcium ions. We propose a model for how PIP2 may bind to calcium-free gelsolin based on the ATP-binding site.     

Structure- and Ligand-Based Virtual Screening Identifies New Scaffolds for Inhibitors of the Oncoprotein MDM2

A major challenge in the field of ligand discovery is to identify chemically useful fragments that can be developed into inhibitors of specific protein-protein interactions. Low molecular weight fragments (with molecular weight less than 250 Da) are likely to bind weakly to a protein’s surface. Here we use a new virtual screening procedure which uses a combination of similarity searching and docking to identify chemically tractable scaffolds that bind to the p53-interaction site of MDM2. The binding has been verified using capillary electrophoresis which has proven to be an excellent screening method for such small, weakly binding ligands.     

Targeting cysteine rich C1 domain of Scaffold protein Kinase Suppressor of Ras (KSR) with anthocyanidins and flavonoids – a binding affinity characterization study

Kinase Suppressor of Ras (KSR) is a molecular scaffold that interacts with the core kinase components of the ERK cascade, Raf, MEK, ERK to provide spatial and temporal regulation of Ras-dependent ERK cascade signaling. Interruption of this mechanism can have a high influence in inhibiting the downstream signaling of the mutated tyrosine kinase receptor kinase upon ligand binding. Still none of the studies targeted to prevent the binding of Raf, MEK binding on kinase suppressor of RAS. In that perspective the cysteine rich C1 domain of scaffold proteins kinase suppressor of Ras-1 was targeted rather than its ATP binding site with small ligand molecules like flavones and anthocyanidins and analyzed through insilico docking studies. The binding energy evaluation shows the importance of hydroxyl groups at various positions on the flavone and anthocyanidin nucleus. Over all binding interaction shows these ligands occupied the potential sites of cysteine rich C1 domain of scaffold protein KSR.     

Probing molecular docking in a charged model binding site

A model binding site was used to investigate charge-charge interactions in molecular docking. This simple site, a small (180A(3)) engineered cavity in cyctochrome c peroxidase (CCP), is negatively charged and completely buried from solvent, allowing us to explore the balance between electrostatic energy and ligand desolvation energy in a system where many of the common approximations in docking do not apply. A database with about 5300 molecules was docked into this cavity. Retrospective testing with known ligands and decoys showed that overall the balance between electrostatic interaction and desolvation energy was captured. More interesting were prospective docking scre"ens that looked for novel ligands, especially those that might reveal problems with the docking and energy methods. Based on screens of the 5300 compound database, both high-scoring and low-scoring molecules were acquired and tested for binding. Out of 16 new, high-scoring compounds tested, 15 were observed to bind. All of these were small heterocyclic cations. Binding constants were measured for a few of these, they ranged between 20microM and 60microM. Crystal structures were determined for ten of these ligands in complex with the protein. The observed ligand geometry corresponded closely to that predicted by docking. Several low-scoring alkyl amino cations were also tested and found to bind. The low docking score of these molecules owed to the relatively high charge density of the charged amino group and the corresponding high desolvation penalty. When the complex structures of those ligands were determined, a bound water molecule was observed interacting with the amino group and a backbone carbonyl group of the cavity. This water molecule mitigates the desolvation penalty and improves the interaction energy relative to that of the "naked" site used in the docking screen. Finally, six low-scoring neutral molecules were also tested, with a view to looking for false negative predictions. Whereas most of these did not bind, two did (phenol and 3-fluorocatechol). Crystal structures for these two ligands in complex with the cavity site suggest reasons for their binding. That these neutral molecules do, in fact bind, contradicts previous results in this site and, along with the alkyl amines, provides instructive false negatives that help identify weaknesses in our scoring functions. Several improvements of these are considered.