Systematic identification of the druggable interactions between human protein kinases and naturally occurring compounds in endometriosis

Diverse kinase signaling pathways have been involved in the pathogenesis of endometriosis (EM), which can be modulated either by directly targeting the hub kinases or by indirectly regulating marginal members in the pathways. Here, a systematic kinase–inhibitor interaction profile was created for 8 naturally occurring compounds against 20 human protein kinases. The compounds are all non-sterid that have been reported as pharmacologically active molecular entities potential for EM therapeutics, while the kinases were curated via gene ontology terms enriched from the gene co-citation network with EM. The resulting profile was analyzed at structural, energetic and dynamic levels to identify druggable kinase–compound interactions. The compounds Gossypol, Curcumin and EGCG showed a similar interaction profile across these kinases; they can bind tightly to the top-listed kinases in gene ontology, while the compounds Marrubiin, Apigenin and DIM were predicted to exhibit generally weak affinity for the 20 curated kinases. The JNK kinase, a MAPK family member, was identified as a putative candidate of druggable target for EM therapeutics; the inhibitory activity of eight naturally occurring compounds as well as a sophisticated kinase inhibitor SP600125 against the JNK was tested using enzymatic activity analysis. As might be expected, the Gossypol and EGCG were determined to have high inhibitory activity at namomolar level (IC₅₀ = 55 and 94 nM, respectively), which are comparable with or better than the positive control SP600125 (IC₅₀ = 76 nM), while other tested compounds exhibited weak inhibition (IC₅₀ > 100 nM) or bad potency (IC₅₀ = n.d.) against the kinase.     

Stereoselective synthesis of (3-aminodecahydro-1,4-methanonaphthalen-2-yl)methanols targeted to the C1 domain of protein kinase C

Protein kinase C (PKC) is a widely studied molecular target for the treatment of cancer and other diseases. We have approached the issue of modifying PKC function by targeting the C1 domain in the regulatory region of the enzyme. By using the X-ray crystal structure of the PKCδ C1b domain combined with molecular modeling, we discovered (3-aminodecahydro-1,4-methanonaphthalen-2-yl)methanol as a novel C1 domain ligand. The stereoselective synthesis of this tricyclic γ-amino alcohol was based on two successive Diels-Alder reactions to construct the six continuous stereocenters of the key intermediate.     

Structure-based modelling,scoring, screening,and in vitro kinase assay of anesthetic pkc inhibitors against a natural medicine library

Protein kinase C (PKC) is an intracellular effector of the inositol phosphate-mediated signal transduction pathway. Evidence is emerging that certain general anaesthetics can influence the activity of PKC by interacting with the regulatory domain of the enzyme, and targeting PKC kinase domain is considered as a strategy to modulate the anaesthetic effects. Here, an integrated method was used to perform virtual screening against a large library of natural compounds for the discovery of new and potent PKC modulators. A number of hits were identified and their inhibitory activity against PKC kinase domain was measured by using a standard kinase assay protocol. Three and five compounds were determined to have high and moderate activities with IC₅₀ values at nanomolar and micromolar levels, respectively. These compounds can be considered as promising lead molecular entities to develop efficacious anaesthetic modulators. Structural examination revealed a variety of nonbonded interactions such as hydrogen bonds, cation-π contacts, and hydrophobic forces across the complex interface of PKC with the identified compounds. This study helps to establish an integrative approach to rational kinase inhibitor discovery by efficiently exploiting various existing natural products.     

Virtual screening of protein kinase C inhibitors from natural product library to modulate general anaesthetic effects

Protein kinase C (PKC) plays a key role in neurotransmission in the central nervous system, and targeting PKC domain is considered as a strategy to modulate the anaesthetic effects. In this study, we described a synthetic pipeline to perform high-throughput virtual screening against a large library of 3D structural natural products released recently in order to discover those potential PKC modulators. A total of 100 natural products with top scores were raised, from which 12 promising candidates were tested to determine their inhibitory potencies against PKC. As might be expected, the promiscuous kinase inhibitor staurosporine showed a high PKC inhibitory activity (IC₅₀ = 64 nM), and other two tested compounds, i.e. fisetin and tetrahydropapaverine, were also highly potent with their activities at nanomolar level (IC₅₀ = 370 and 190, respectively).     

Identification of critical chemical features for Aurora kinase-B inhibitors using Hip-Hop, virtual screening and molecular docking

This study was performed to find the selective chemical features for Aurora kinase-B inhibitors using the potent methods like Hip-Hop, virtual screening, homology modeling, molecular dynamics and docking. The best hypothesis, Hypo1 was validated toward a wide range of test set containing the selective inhibitors of Aurora kinase-B. Homology modeling and molecular dynamics studies were carried out to perform the molecular docking studies. The best hypothesis Hypo1 was used as a 3D query to screen the chemical databases. The screened molecules from the databases were sorted based on ADME and drug like properties. The selective hit compounds were docked and the hydrogen bond interactions with the critical amino acids present in Aurora kinase-B were compared with the chemical features present in the Hypo1. Finally, we suggest that the chemical features present in the Hypo1 are vital for a molecule to inhibit the Aurora kinase-B activity.     

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling,virtual screening based on multiple conformations,and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5’-phosphate and 3’-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3′-hydroxyl group and 5′-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5′-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.     

Homology modeling and molecular dynamics study of GSK3/SHAGGY-like kinase

Although the GSK3/SHAGGY-like kinase is a highly conserved serine/threonine kinase implicated in many signaling pathways in eukaryotes, the lack of knowledge of its three-dimensional (3D) structure has hindered efforts to understand the binding specificities of substrate and catalytic mechanism. To understand the structure-activity relationships, the protein 3D structure was built by using homology modeling based on the known X-ray diffraction structure of Glycogen synthase kinase-3beta (Gsk3beta) and the model structure was further refined using unrestrained molecular dynamics simulations. The research indicates that the general 3D organization of the GSK3/SHAGGY-like kinase is a typical kinase family and comprises an N-terminal domain of beta-sheet and a larger C-terminal domain mainly constituted by alpha-helix. In order to understand the molecular interactions between the natural substrate-ATP and GSK3/SHAGGY-like kinase, a 3D model of the complex ATP-GSK3/SHAGGY-like kinase is developed by molecular docking program, which is helpful to guide the experimental realization and the new mutant designs as well. One important finding is that the identification of the key binding-site residue of Lys69 which plays an important role in the catalysis of GSK3/SHAGGY-like kinase and this is in consistent with experimental observation.     

Computer-aided active-site-directed modeling of the Herpes Simplex Virus 1 and human thymidine kinase

Summary Thymidine kinase (TK), which is induced by Herpes Simplex Virus 1 (HSV1), plays a key role in the antiviral activity of guanine derivatives such as aciclovir (ACV). In contrast, ACV shows only low affinity to the corresponding host cell enzyme. In order to define the differences in substrate binding of the two enzymes on molecular level, models for the three-dimensional (3-D) structures of the active sites of HSV1-TK and human TK were developed. The reconstruction of the active sites started from primary and secondary structure analysis of various kinases. The results were validated to homologous enzymes with known 3-D structures. The models predict that both enzymes consist of a central core -sheet structure, connected by loops and -helices very similar to the overall structure of other nucleotide binding enzymes. The phosphate binding is made up of a highly conserved glycine-rich loop at the N-terminus of the proteins and a conserved region at the C-terminus. The thymidine recognition site was found about 100 amino acids downstream from the phosphate binding loop. The differing substrate specificity of human and HSV1-TK can be explained by amino-acid substitutions in the homologous regions.To achieve a better understanding of the structure of the active site and how the thymidine kinase proteins interact with their substrates, the corresponding complexes of thymidine and dihydroxypropoxyguanine (DHPG) with HSV1 and human TK were built. For the docking of the guanine derivative, the X-ray structure of Elongation Factor Tu (EF-Tu), co-crystallized with guanosine diphosphate, was taken as reference. Fitting of thymidine into the active sites was done with respect to similar interactions found in thymidylate kinase. To complement the analysis of the 3-D structures of the two kinases and the substrate enzyme interactions, site-directed mutagenesis of the thymidine recognition site of HSV1-TK has been undertaken, changing Asp¹⁶² in the thymidine recognition site into Asn. First investigations reveal that the enzymatic activity of the mutant protein is destroyed.     

人类蛋白激酶组的多靶点分子对接系统

蛋白激酶在信号转导、基因转录和蛋白翻译等生物过程起关键性作用,因而与大量人类疾病密切相关。所以,蛋白激酶的抑制剂筛选是抗肿瘤药物开发的热点,正在向基于全激酶组的高通量多靶点筛选模式发展。为了降低大规模实验筛选的成本,提高成功率,本文构建人类蛋白激酶组的多靶点分子对接系统,对抑制剂-激酶组的相互作用进行预测。我们首先利用同源模建方法,对人类激酶组约500个激酶变异体的催化域进行结构建模;接着以催化域结构模型为受体,用已知激酶抑制剂进行分子对接,对抑制剂与各激酶变异体的结合亲和力进行了定量计算。结果显示,本文所建立的多靶点分子对接系统可以准确预测抑制剂与激酶变异体的相互作用,结合自由能的计算值与实验值有很强的相关性。所以,该分子对接系统可用于多靶点激酶抑制剂的计算筛选,为激酶抑制剂开发与抗肿瘤药物设计提供理论依据。     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.