Key Features for Designing Phosphodiesterase-5 Inhibitors

Abstract

Phosphodiesterase superfamily is the key regulator of 3′,5′-cyclic guanosine monophosphate (cGMP) decomposition in human body. Phosphodiesterase-5 (PDE-5) inhibitors, sildenafil, vardenafil and tadalafil, are well known oral treatment for males with erectile dysfunction. To investigate the inhibitory effects of traditional Chinese medicine (TCM) compounds to PDE-5, we performed both ligand-based and structure-based studies on this topic. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted to construct three dimensional quantitative structure-activity relationship (3D-QSAR) models of series of known PDE-5 inhibitors. The predictive models had cross-validated, q², and non cross-validated coefficient, r², values of 0.791 and 0.948 for CoMFA and 0.724 and 0.908 for CoMSIA. These two 3D-QSAR models were used to predict activity of TCM compounds. Docking simulations were performed to further analyze the binding mode of training set and TCM compounds. A putative binding model was proposed based on CoMFA and CoMSIA contour maps and docking simulations; formation of pi-stacking, water bridge and specific hydrogen bonding were deemed important interactions between ligands and PDE-5. Of our TCM compounds, engeletin, satisfied our binding model, and hence, emerged as PDE-5 inhibitor candidate.

Using this study as an example, we demonstrated that docking should be conducted for qualitative purposes, such as identifying protein characteristics, rather than for quantitative analyses that rank compound efficacy based on results of scoring functions. Prediction of compound activity should be reserved for QSAR analyses, and scoring functions and docking scores should be used for preliminary screening of TCM database (http://tcm.cmu.edu.tw/index.php).     

Structural and Functional Characterization of the JH2 Pseudokinase Domain of JAK Family Tyrosine Kinase 2 (TYK2)

JAK (Janus family of cytoplasmic tyrosine kinases) family tyrosine kinase 2 (TYK2) participates in signaling through cytokine receptors involved in immune responses and inflammation. JAKs are characterized by dual kinase domain: a tyrosine kinase domain (JH1) that is preceded by a pseudokinase domain (JH2). The majority of disease-associated mutations in JAKs map to JH2, demonstrating its central regulatory function. JH2s were considered catalytically inactive, but JAK2 JH2 was found to have low autoregulatory catalytic activity. Whether the other JAK JH2s share ATP binding and enzymatic activity has been unclear. Here we report the crystal structure of TYK2 JH2 in complex with adenosine 5′-O-(thiotriphosphate) (ATP-γS) and characterize its nucleotide binding by biochemical and biophysical methods. TYK2 JH2 did not show phosphotransfer activity, but it binds ATP and the nucleotide binding stabilizes the protein without inducing major conformational changes. Mutation of the JH2 ATP-binding pocket increased basal TYK2 phosphorylation and downstream signaling. The overall structural characteristics of TYK2 JH2 resemble JAK2 JH2, but distinct stabilizing molecular interactions around helix αAL in the activation loop provide a structural basis for differences in substrate access and catalytic activities among JAK family JH2s. The structural and biochemical data suggest that ATP binding is functionally important for both TYK2 and JAK2 JH2s, whereas the regulatory phosphorylation appears to be a unique property of JAK2. Finally, the co-crystal structure of TYK2 JH2 complexed with a small molecule inhibitor demonstrates that JH2 is accessible to ATP-competitive compounds, which offers novel approaches for targeting cytokine signaling as well as potential therapeutic applications.     

Synthesis and preliminary structure-activity relationship study of 2-aryl-2H-pyrazolo[4,3-c]quinolin-3-ones as potential checkpoint kinase 1 (Chk1) inhibitors

The serine-threonine checkpoint kinase 1 (Chk1) plays a critical role in the cell cycle arrest in response to DNA damage. In the last decade, Chk1 inhibitors have emerged as a novel therapeutic strategy to potentiate the anti-tumour efficacy of cytotoxic chemotherapeutic agents. In the search for new Chk1 inhibitors, a congeneric series of 2-aryl-2?H-pyrazolo[4,3-c]quinolin-3-one (PQ) was evaluated by in-vitro and in-silico approaches for the first time. A total of 30 PQ structures were synthesised in good to excellent yields using conventional or microwave heating, highlighting that 14 of them are new chemical entities. Noteworthy, in this preliminary study two compounds 4e₂ and 4h₂ have shown a modest but significant reduction in the basal activity of the Chk1 kinase. Starting from these preliminary results, we have designed the second generation of analogous in this class and further studies are in progress in our laboratories.     

Synthesis and primary cytotoxicity evaluation of arylmethylenenaphthofuranones derivatives

New series of 2(or 3)-arylmethylenenaphtho[2,1-b]furan-3(or 2)-ones were synthesized, characterized and tested for anticancer properties in vitro. The target compounds were prepared by Knoevenagel coupling between the naphthofuranones 3, 28–30 and formyl derivatives. 2-(4-Oxo-1-benzopyran-3-ylmethylene)naphtho[2,1-b]furan-3-one 36 was the most active compound (IC₅₀ (L1210) = 1.6 μM). These compounds were also evaluated, in an independent manner, as inhibitors of Src protein tyrosine kinase, but only minor activity was observed.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

2-Alkyl(aryl)-quinazolin-4(3H)-thiones, 2-R-(quinazolin-4(3H)-ylthio)carboxylic acids and amides: synthesis,molecular docking,antimicrobial and anticancer properties

In this study, a series of novel 2-alkyl(aryl)-quinazolin-4(3H)-thiones, 2-R-(quinazolin-4(3H)-ylthio)carboxylic acids and amides were synthesized and evaluated for antimicrobial and anticancer activities. Their structure was confirmed by elemental analysis and spectral data (FT-IR, LC-MS, ¹H-NMR). Antimicrobial activity was tested in vitro against Staphylococcus aureus, Enterococcus faecalis, Enterobacter aerogenes, Pseudomonas aeruginosa, Escherichia coli, Klebsiella pneumonia, Candida albicans and NCI in vitro preliminary anticancer activity against nine different cancer types. The most active antibacterial and antifungal compounds were: 2.1, 2.2 and 2.4. The introduction of the carboxylic acid or amide residue into the fourth position of quinazolin-4(3H)-thione resulted in the absence of antimicrobial activity. Substance 3.8 inhibited renal cancer UO-31 line and 2.18 – leukemia CCRF-CEM. The results of in silico molecular docking for DHFR and CK2 kinase had no correlation with in vitro properties, proposing the presence of other biological activity pathways.     

Potential Inhibitors of Angiogenesis. Part I: 3-(Imidazol-4(5)-ylmethylene)indolin-2-ones

The synthesis and pharmacological evaluation of new 3-(imidazol-4(5)-ylmethylene)indolin-2-ones, analogues of SU-5416, are reported. The final compounds 20-51 were obtained by Knoevenagel coupling between the substituted indolin-2-ones 1-15 and either the formylimidazole derivatives 16-18 or 2-formyl-3,5-dimethylpyrrole 19. Methylation at the nitrogen atom of the indolin-2-one and/or imidazole moities was carried out in the presence of the couple NaH/DMF. A Mannich reaction afforded the 1-dimethylaminomethyl derivatives 43 and 48. The antiangiogenic activity of these compounds was evaluated in a three dimensional in vitro rat aortic ring assay. In this test, compound 20 induced a decrease of angiogenesis comparable to that observed with SU-5416; the vascular density indexes at 1 μM were 30±18 and 22±4 % of control, respectively. The compounds were also evaluated, in an independent manner, as inhibitors of the human EGF-receptor tyrosine kinase activity. As expected, only minor activities were observed with four compounds, out of thirty-one, exerting inhibitory effects in the range of 40-55 % at 10 μM concentration.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

De novo design of quinazoline derivatives as CDK2 inhibitors: 3D-QSAR,molecular fragment replacement and Volsurf predictions

Cyclin-dependent kinase 2 (CDK2) has appeared as an important drug target over the years with a multitude of therapeutic potentials. To design compounds with enhanced inhibitory potencies against CDK2, 3D-QSAR and molecular fragment replacement studies were performed on the pyrazolo[4,3-h]quinazoline derivatives, a class of potent CDK2 inhibitors. The contours of 3D-QSAR model revealed important structural features of the inhibitors related to the active site of CDK2. Based on the pyrazolo[4,3-h]quinazoline core, the different substituents at three important points were replaced with diverse molecular fragments. The compounds resulting from fragments assembly with pyrazolo[4,3-h]quinazoline core were then scored with the robust 3D-QSAR model. Furthermore, the absorption, distribution, metabolism and excretion properties of these compounds were predicted by Volsurf to eliminate inappropriate compounds. Thirty-one new potential compounds were finally obtained. These results initiated us to further optimise and design new potential inhibitors.     

Both tissue inhibitors of metalloproteinases-1 (TIMP-1) and TIMP-2 activate Ras but through different pathways

Tissue inhibitors of metalloproteinases-1 (TIMP-1) and TIMP-2 have growth-stimulating activity for a wide range of cell types. Ras, which comprises a family of three members, i.e, Ha-Ras, Ki-Ras, and H-Ras, is known to participate in growth control in all its facets, including cell proliferation, transformation, differentiation, and apoptosis. In this study, we tested the hypothesis that Ras might be involved in the cell growth-promoting activity of TIMPs. Using MG-63 human osteosarcoma cells, we demonstrated that both TIMP-1 and TIMP-2 caused an increase in the Ras-GTP level in a dose-dependent manner. Our previous results indicated that TIMP-1 activity is mediated through the tyrosine kinase (TYK)/mitogen-activated protein kinase (MAPK) pathway. Here, we demonstrated that Ras activation by TIMP-1 was inhibited by a specific TYK inhibitor, herbimycin A, suggesting that the TYK/MAPK signaling pathway was involved in Ras activation by TIMP-1. However, the activation of Ras by TIMP-2 was inhibited by an inhibitor specific for cyclic AMP-dependent protein kinase (PKA), H89, suggesting the involvement of the PKA-mediated pathway. Furthermore, TIMP-2 promoted the formation of a complex between Ras-GTP and phosphoinositide 3-kinase.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

Synthesis,activity evaluation,and docking analysis of barbituric acid aryl hydrazone derivatives as RSK2 inhibitors

The 90 kDa ribosomal S6 kinases (RSKs), especially RSK2, have attracted attention for the development of new anticancer agents. Through structural optimization of the hit compound 1 from our previous study, a series of barbituric acid aryl hydrazone analogues were designed and synthesized as potential RSK2 inhibitors. The most potent one, compound 9, showed a higher activity against RSK2 with an IC₅₀ value of 1.95 μM. To analyze and elucidate their structure-activity relationship, the homology model of RSK2 N-terminal kinase domain was built and molecular docking simulations were performed, which provide helpful clues to design new inhibitors with desired activities.     

Molecular dynamics simulations reveal structural insights into inhibitor binding modes and mechanism of casein kinase II inhibitors

Casein kinase-II, a member of protein kinase family, plays significant role in different cellular processes such as cell growth, differentiation, proliferation, gene expression, and embryogenesis. Being a potent suppressor of apoptosis, it serves as a significant link for its association with various types of malignancies such as colorectal and breast cancer. To overcome its pathological role in various cancerous diseases, CK-II procures great consideration as a therapeutic target. This study aimed at understanding the binding mechanism and structural properties of benzimidazole derivatives by utilizing various computational tools including docking simulation, three-dimensional quantitative structure activity relationships and molecular dynamic simulation. Structure-based 3D-QSAR techniques such as CoMFA and CoMSIA models, were established from the conformations gained by protein–ligand docking approach. The attained models have showed a good extrapolative power for internal as well as external validation. Moreover, MD simulation was carried out to explain the detailed binding mechanism and the comparison of inhibitor’s binding mode with diverse biological activities. A good correlation was observed among docking studies, MD results, and contour map analysis. Interestingly new molecules were designed using detail structural information from MD simulation, showed higher potency of inhibition (pIC₅₀ 7.6–7.7) compare to the most active compound of the series.     

(Thiazol-2-yl)hydrazone derivatives from acetylpyridines as dual inhibitors of MAO and AChE: synthesis,biological evaluation and molecular modeling studies

Abstract

Several (thiazol-2-yl)hydrazone derivatives from 2-, 3- and 4-acetylpyridine were synthesized and tested against human monoamine oxidase (hMAO) A and B enzymes. Most of them had an inhibitory effect in the low micromolar/high nanomolar range, being derivatives of 4-acetylpyridine selective hMAO-B inhibitors also at low nanomolar concentrations. The structure–activity relationship, as confirmed by molecular modeling studies, proved that the pyridine ring linked to the hydrazonic nitrogen and the substituted aryl moiety at C4 of the thiazole conferred the inhibitory effects on hMAO enzymes. Successively, the strongest hMAO-B inhibitors were tested toward acetylcholinesterase (AChE) and the most interesting compound showed activity in the low micromolar range. Our results suggest that this scaffold could be further investigated for its potential multi-targeted role in the discovery of new drugs against the neurodegenerative diseases.     

Designing of benzothiazole derivatives as promising EGFR tyrosine kinase inhibitors: a pharmacoinformatics study

Abstract

Benzothiazole derivatives represent an important class of therapeutic chemical agents and are widely used for interesting biological activities and therapeutic functions including anticancer, antitumor and antimicrobial. In this study, we have performed similarity/substructure-based search of eMolecule database to find out promising benzothiazole derivatives as EGFR tyrosine kinase inhibitors. Several screening criteria that included molecular docking, pharmacokinetics and synthetic accessibility were used on initially derived about 7000 molecules consisting of benzothiazole as major component. Finally, four molecules were found to be promising EGFR tyrosine kinase inhibitors. The best docked pose of each molecule was considered for binding interactions followed by molecular dynamics (MD) and binding energy calculation. Molecular docking clearly showed the final proposed derivatives potential to form a number of binding interactions. MD simulation trajectories undoubtedly indicated that the EGFR protein becomes stable when proposed derivatives bind to the receptor cavity. Strong binding affinity was found for all molecules toward the EGFR which was substantiated by the binding energy calculation using the MM-PBSA approach. Therefore, proposed benzothiazole derivatives may be promising EGFR tyrosine kinase inhibitors for potential application as cancer therapy.

Communicated by Ramaswamy H. Sarma     

3D-QSAR studies of checkpoint kinase 1 inhibitors based on molecular docking and CoMFA

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on a series of substituted 1,4-dihydroindeno[1,2-c]pyrazoles inhibitors, using molecular docking and comparative molecular field analysis (CoMFA). The docking results from GOLD 3.0.1 provide a reliable conformational alignment scheme for the 3D-QSAR model. Based on the docking conformations and alignments, highly predictive CoMFA model was built with cross-validated q ² value of 0.534 and non-cross-validated partial least-squares analysis with the optimum components of six showed a conventional r ² value of 0.911. The predictive ability of this model was validated by the testing set with a conventional r ² value of 0.812. Based on the docking and CoMFA, we have identified some key features of the 1,4-dihydroindeno[1,2-c]pyrazoles derivatives that are responsible for checkpoint kinase 1 inhibitory activity. The analyses may be used to design more potent 1,4-dihydroindeno[1,2-c]pyrazoles derivatives and predict their activity prior to synthesis.