Discovery of novel checkpoint kinase 1 inhibitors by virtual screening based on multiple crystal structures

Incorporating receptor flexibility is considered crucial for improvement of docking-based virtual screening. With an abundance of crystallographic structures freely available, docking with multiple crystal structures is believed to be a practical approach to cope with protein flexibility. Here we describe a successful application of the docking of multiple structures to discover novel and potent Chk1 inhibitors. Forty-six Chk1 structures were first compared in single structure docking by predicting the binding mode and recovering known ligands. Combinations of different protein structures were then compared by recovery of known ligands and an optimal ensemble of Chk1 structures were selected. The chosen structures were used in the virtual screening of over 60?000 diverse compounds for Chk1 inhibitors. Six novel compounds ranked at the top of the hits list were tested experimentally, and two of these compounds inhibited Chk1 activity-the best with an IC(50) value of 9.6 μM. Further study indicated that achieving a better enrichment and identifying more diverse compounds was more likely using multiple structures than using only a single structure even when protein structures were randomly selected. Taking into account conformational energy difference did not help to improve enrichment in the top ranked list.     

Identification of novel inhibitors of mitogen-activated protein kinase phosphatase-1 with structure-based virtual screening

Mitogen-activated protein kinase phosphatase-1 (MKP-1) has proved to be an attractive target for the development of therapeutics for the treatment of cancer. We report the first example for a successful application of the structure-based virtual screening to identify the novel inhibitors of MKP-1. It is shown that the efficiency of virtual screening can be enhanced significantly by the incorporation of a new solvation energy term in the scoring function. The newly found inhibitors have desirable physicochemical properties as a drug candidate and reveal a moderate potency with IC₅₀ values ranging from 20 to 50 μM. Therefore, they deserve a consideration for further development by structure–activity relationship studies to optimize the inhibitory activities. Structural features relevant to the stabilization of the inhibitors in the active site of MKP-1 are discussed in detail.     

Identification of novel inhibitors of tropomyosin-related kinase A through the structure-based virtual screening with homology-modeled protein structure

Tropomyosin-related kinase A (TrkA) is a promising target for the development of cancer and pain therapeutics. Here, we report the first successful example of the use of a structure-based virtual screening to identify novel TrkA inhibitors. The accuracy of the virtual screening was improved by introducing an accurate solvation free energy term into the original AutoDock scoring function. We applied a drug design protocol involving homology modeling, docking analysis of a large chemical library, and enzyme inhibition assays to identify six structurally diverse TrkA inhibitors with K(d) values ranging from 3 to 40 μM. The significant potencies and good physicochemical properties of these drug candidates strongly support their consideration in a development effort that would involve structure-activity relationship (SAR) studies to optimize the inhibitory activities. We also addressed the structural and energetic features associated with binding of the newly identified inhibitors in the ATP-binding site of TrkA. The results indicate that any structural modifications introduced for the purpose of enhancing the activity of TrkA inhibitors should maximize the attractive interactions within the ATP-binding site and simultaneously minimize the desolvation cost for complexation.     

Identification of novel human renin inhibitors through a combined approach of pharmacophore modelling,molecular DFT analysis and in silico screening

Renin is an aspartyl protease of the renin–angiotensin system (RAS) and the first enzyme of the biochemical pathway for the generation of angiotensin II – a potent vasoconstrictor involved in the maintenance of cardiovascular homeostasis and the regulation of blood pressure. High enzymatic specificity of renin and its involvement in the catalysis of the rate-limiting step of the RAS hormone system qualify it as a good target for inhibition of hypertension and other associated diseases. Ligand-based pharmacophore model (Hypo1) was generated from a training set of 24 compounds with renin inhibitory activity. The best hypothesis consisted of one Hydrogen Bond Acceptor (HBA), three Hydrophobic Aliphatic (HY-Al) and one Ring Aromatic (AR) features. This well-validated pharmacophore hypothesis (correlation coefficient 0.95) was further utilized as a 3D query to screen database compounds, which included structures from two natural product repositories. These screened compounds were further analyzed for drug-likeness and ADMET studies. The compounds which satisfied the qualifying criteria were then subjected to molecular docking and Density Functional Theory (DFT) analysis in order to discern their atomic level interactions at the active site of the 3D structure of rennin. The pharmacophore-based modelling that has been used to generate the novel findings of the present study would be an avant-garde approach towards the development of potent inhibitors of renin.     

Identification of four novel phosphorylation sites in estrogen receptor α: impact on receptor-dependent gene expression and phosphorylation by protein kinase CK2

Background

Human S100A12 is a member of the S100 family of EF-hand calcium-modulated proteins that are associated with many diseases including cancer, chronic inflammation and neurological disorders. S100A12 is an important factor in host/parasite defenses and in the inflammatory response. Like several other S100 proteins, it binds zinc and copper in addition to calcium. Mechanisms of zinc regulation have been proposed for a number of S100 proteins e.g. S100B, S100A2, S100A7, S100A8/9. The interaction of S100 proteins with their targets is strongly dependent on cellular microenvironment.

Results

The aim of the study was to explore the factors that influence S100A12 oligomerization and target interaction. A comprehensive series of biochemical and biophysical experiments indicated that changes in the concentration of calcium and zinc led to changes in the oligomeric state of S100A12. Surface plasmon resonance confirmed that the presence of both calcium and zinc is essential for the interaction of S100A12 with one of its extracellular targets, RAGE – the Receptor for Advanced Glycation End products. By using a single-molecule approach we have shown that the presence of zinc in tissue culture medium favors both the oligomerization of exogenous S100A12 protein and its interaction with targets on the cell surface.

Conclusion

We have shown that oligomerization and target recognition by S100A12 is regulated by both zinc and calcium. Our present work highlighted the potential role of calcium-binding S100 proteins in zinc metabolism and, in particular, the role of S100A12 in the cross talk between zinc and calcium in cell signaling.     

Identification and neuroprotective evaluation of a potential c-Jun N-terminal kinase 3 inhibitor through structure-based virtual screening and in-vitro assay

Journal of Computer-Aided Molecular Design - The c-Jun N-terminal kinase 3 (JNK3) signaling cascade is activated during cerebral ischemia leading to neuronal damage. The present study was carried...     

Rapid identification of substrates for novel proteases using a combinatorial peptide library

Fluorogenic substrates for assaying novel proteolytic enzymes could be rapidly identified using an easy, solid-phase combinatorial assay technology. The methodology was validated with leader peptidase of Escherichia coli using a subset of an intramolecularly quenched fluorogenic peptide library. The technique was extended toward the discovery of substrates for a new aspartic protease of pharmaceutical relevance (human napsin A). We demonstrated for the first time known to us that potent fluorogenic substrates can be discovered using extracts of cells expressing recombinant enzyme to screen the peptide library. The straightforward and rapid optimization of protease substrates greatly facilitates the drug discovery process by speeding up the development of high throughput screening assays and thus helps more effective exploitation of the enormous body of information and chemical structures emerging from genomics and combinatorial chemistry technologies.     

Identification of novel S-adenosyl-L-homocysteine hydrolase inhibitors through homology-model-based virtual screening, synthesis, and biological evaluation

The present study describes a successful application of computational approaches to identify novel Leishmania donovani (Ld) AdoHcyase inhibitors utilizing the differences for Ld AdoHcyase NAD(+) binding between human and Ld parasite. The development and validation of the three-dimensional (3D) structures of Ld AdoHcyase using the L. major AdoHcyase as template has been carried out. At the same time, cloning of the Ld AdoHcyase gene from clinical strains, its overexpression and purification have been performed. Further, the model was used in combined docking and molecular dynamics studies to validate the binding site of NAD in Ld. The hierarchical structure based virtual screening followed by the synthesis of five active hits and enzyme inhibition assay has resulted in the identification of novel Ld AdoHcyase inhibitors. The most potent inhibitor, compound 5, may serve as a "lead" for developing more potent Ld AdoHcy hydrolase inhibitors as potential antileishmanial agents.     

Design and synthesis of useful intermediates for novel lipoxygenase substrates through enzymatic resolution

Unnatural lipoxygenase substrates carrying spacing modifiers with a non-ionic hydroxy terminus and with methylene (flanked by the cis,cis diene moiety) pro-(S)-hydrogen can be synthesized from the intermediates 1a–1c. These intermediates are conveniently synthesized via enzymatic resolution with lipase in organic solvents.     

Impaired phosphorylation and mis-localization of Bub1 and BubR1 are responsible for the defective mitotic checkpoint function in Brca2-mutant thymic lymphomas

Breast cancer susceptibility gene, BRCA2, is a tumor suppressor and individuals who inherit one defected copy of BRCA2 allele experience early onset breast cancer or ovarian cancer accompanied by the loss of the wild type allele. Mouse model for Brca2 mutation shows growth retardation and paradoxical occurrence of thymic lymphomas. Thymic lymphomas from Brca2-mutant mice harbor mutations in p53, Bub1, and BubR1, which function as mitotic checkpoint proteins. Therefore, interplay between Brca2 and mitotic checkpoint has been suggested in the maintenance of genetic fidelity, although it has not been assessed whether the unique mutations in Bub1 and BubR1 found in Brca2-mutant mice are responsible for the abolishment of mitotic checkpoint function. This report demonstrates that Bub1 and BubR1 mutant proteins from Brca2-/- thymic lymphomas have defects in the phosphorylation and kinetochore localization after spindle damage. Thus, the mutations of Bub1 and BubR1 found in Brca2- mutant mice indeed are responsible for the chromosome instability in Brca2-mutated tumors.     

Pharmacophore and docking-based virtual screening approach for the design of new dual inhibitors of Janus kinase 1 and Janus kinase 2

Janus kinase 1 and 2, non-receptor protein tyrosine kinases, are implicated in various cancerous diseases. Involvement of these two enzymes in the pathways that stimulate cell proliferation in cancerous conditions makes them potential therapeutic targets for designing new dual-targeted agents for the treatment of cancer. In the present study, two separate pharmacophore models were developed and the best models for JAK1 (AAADH.25) and JAK2 (ADRR.92) were selected on the basis of their external predictive ability. Both models were subjected to a systematic virtual screening (VS) protocol using a PHASE database of 1.5 million molecules. The hits retrieved in VS were investigated for ADME properties to avoid selection of molecules with a poor pharmacokinetic profile. The molecules considered to be within the range of acceptable limits of ADME properties were further employed for docking simulations with JAK1 and JAK2 proteins to explore the final hits that possess structural features of both pharmacophore models as well as display essential interactions with both of them. Thus, the new molecules obtained in this way should show inhibitory activity against JAK1 and JAK2 and may serve as novel therapeutic agents for the treatment of cancerous disease conditions.     

Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.     

Identification of inhibitors against p90 ribosomal S6 kinase 2 (RSK2) through structure-based virtual screening with the inhibitor-constrained refined homology model

P90 ribosomal S6 kinase 2 (RSK2), which was shown to be overexpressed in human cancers, is a serine/threonine kinase and a potential target for cancer treatment. RSK2 comprises two terminal kinase domains (NTKD and CTKD) that can be inhibited by binding with different types of inhibitors at the ATP binding sites. In the absence of a crystal structure of RSK2, we constructed a model for the 3D structure of the RSK2 NTKD by homology modeling and stepwise constrained refinement with the reported inhibitors using a molecular docking method. Structure-based virtual screening was subsequently performed against a library containing commercially available compounds using the refined model. This resulted in the identification of seven novel RSK2 inhibitors with IC?? values ranging from 2.4 to 14.45 μM.     

Mapping of phosphorylation sites in human MSK1 activated by a novel interaction with MRK‐β

The desire to map reliable phosphorylation signaling network has motivated the development of high‐performance techniques. Targeted biochemical studies and updated methods employing MS techniques are most used in mapping the phosphorylation sites and verifying novel interactions of kinases. Previously, we have established a novel method to efficiently facilitate more comprehensive, accurate phosphorylation site mapping of individual phosphoproteins by using combination of multiple stage MS analysis with target‐decoy database search against the much smaller targeted database. In this study, by applying this method, we have identified the phosphorylation sites in human MSK1 mitogen‐ and stress‐activated protein kinase 1), which has been proved to be a multi‐phosphorylated kinase that plays key roles in various cell functions, activated by a novel interaction with MRK‐β. The results show that this method can find out not only those previously identified active sites in MSK1, but also some novel phosphorylated sites, which correlates with biochemical evidence that, besides p38 and extracellular signal‐regulated kinase, MRK‐β could also activate MSK1 through direct interaction. Hence, we conclude this method is sensitive and reliable as expected and it can be further combined with automated screening and biochemical study in efficiently building up a more comprehensive phosphoprotein network.     

Development of off-line and on-line capillary electrophoresis methods for the screening and characterization of adenosine kinase inhibitors and substrates

Fast and convenient CE assays were developed for the screening of adenosine kinase (AK) inhibitors and substrates. In the first method, the enzymatic reaction was performed in a test tube and the samples were subsequently injected into the capillary by pressure and detected by their UV absorbance at 260 nm. An MEKC method using borate buffer (pH 9.5) containing 100 mM SDS (method A) was suitable for separating alternative substrates (nucleosides). For the CE determination of AMP formed as a product of the AK reaction, a phosphate buffer (pH 7.5 or 8.5) was used and a constant current (95 microA) was applied (method B). The methods employing a fused-silica capillary and normal polarity mode provided good resolution of substrates and products of the enzymatic reaction and a short analysis time of less than 10 min. To further optimize and miniaturize the AK assays, the enzymatic reaction was performed directly in the capillary, prior to separation and quantitation of the product employing electrophoretically mediated microanalysis (EMMA, method C). After hydrodynamic injection of a plug of reaction buffer (20 mM Tris-HCl, 0.2 mM MgCl2, pH 7.4), followed by a plug containing the enzyme, and subsequent injection of a plug of reaction buffer containing 1 mM ATP, 100 microM adenosine, and 20 microM UMP as an internal standard (I.S.), as well as various concentrations of an inhibitor, the reaction was initiated by the application of 5 kV separation voltage (negative polarity) for 0.20 min to let the plugs interpenetrate. The voltage was turned off for 5 min (zero-potential amplification) and again turned on at a constant current of -60 microA to elute the products within 7 min. The method employing a polyacrylamide-coated capillary of 20 cm effective length and reverse polarity mode provided good resolution of substrates and products. Dose-response curves and calculated K(i) values for standard antagonists obtained by CE were in excellent agreement with data obtained by the standard radioactive assay.     

A CE‐based assay for human protein kinase CK2 activity measurement and inhibitor screening

A new assay for protein kinase CK2 activity determination based on the quantification of a phosphorylated substrate was developed. The common CK2 substrate peptide RRRDDDSDDD, conjugated with the fluorophore 5‐[(2‐aminoethyl)amino]naphthalene‐1‐sulfonic acid at the C‐terminus served as the analyte. By means of CZE using 2 mol/L acetic acid as electrolyte and UV detection at 214 nm, the non‐phosphorylated and the phosphorylated peptide variants could be resolved within 6 min from a complex assay mixture. By this means, activity of human CK2 could be monitored by a kinetic, as well as an endpoint, method. Inhibition of human recombinant CK2 holoenzyme by 6‐methyl‐1,3,8‐trihydroxyanthraquinone and 4,5,6,7‐tetrabromobenzotriazole resulted in IC₅₀ values of 1.33 and 0.27 μM, respectively, which were similar to those obtained with the standard radiometric assay. These results suggest that the CE/UV strategy described here is a straightforward assay for CK2 inhibitor testing.     

Identification and characterization of the novel nuclease activity of human phospholipid scramblase 1

Background

Human phospholipid scramblase 1 (hPLSCR1) was initially identified as a Ca²⁺ dependent phospholipid translocator involved in disrupting membrane asymmetry. Recent reports revealed that hPLSCR1 acts as a multifunctional signaling molecule rather than functioning as scramblase. hPLSCR1 is overexpressed in a variety of tumor cells and is known to interact with a number of protein molecules implying diverse functions.

Results

In this study, the nuclease activity of recombinant hPLSCR1 and its biochemical properties have been determined. Point mutations were generated to identify the critical region responsible for the nuclease activity. Recombinant hPLSCR1 exhibits Mg²⁺ dependent nuclease activity with an optimum pH and temperature of 8.5 and 37 °C respectively. Experiments with amino acid modifying reagents revealed that histidine, cysteine and arginine residues were crucial for its function. hPLSCR1 has five histidine residues and point mutations of histidine residues to alanine in hPLSCR1 resulted in 60 % loss in nuclease activity. Thus histidine residues could play a critical role in the nuclease activity of hPLSCR1.

Conclusions

This is the first report on the novel nuclease activity of the multi-functional hPLSCR1. hPLSCR1 shows a metal dependent nuclease activity which could play a role in key cellular processes that needs to be further investigated.

     

Synthesis of 5-substituted-1H-indol-2-yl-1H-quinolin-2-ones: a novel class of KDR kinase inhibitors

[reaction: see text] A number of approaches for the synthesis of the 1H-indol-2-yl-1H-quinolin-2-one ring system found in the potent and selective KDR kinase inhibitor 1 are described. The preparation and reaction of trimethylsilylnitrobenzene 26 with 2-methoxy-3-quinolinecarboxaldehyde 28 afforded alcohol 30, which was the key intermediate for the preparation of the target compounds. Conversion of alcohol 30 to either nitroketone 36 or nitrostyrene 45 set the stage for reductive cyclization and the formation of indole 25. The quinolin-2-one functionality was unmasked in the last step to provide compound 1 in 56-60% overall yield from readily available starting materials.