Pharmacological characterisation of novel adenosine A3 receptor antagonists.

The adenosine A3 receptor (A3R) belongs to a family of four adenosine receptor (AR) subtypes which all play distinct roles throughout the body. A3R antagonists have been described as potential treatments for numerous diseases including asthma. Given the similarity between (adenosine receptors) orthosteric binding sites, obtaining highly selective antagonists is a challenging but critical task. Here we screen 39 potential A3R, antagonists using agonist-induced inhibition of cAMP. Positive hits were assessed for AR subtype selectivity through cAMP accumulation assays. The antagonist affinity was determined using Schild analysis (pA2 values) and fluorescent ligand binding. Structure-activity relationship investigations revealed that loss of the 3-(dichlorophenyl)-isoxazolyl moiety or the aromatic nitrogen heterocycle with nitrogen at α-position to the carbon of carboximidamide group significantly attenuated K18 antagonistic potency. Mutagenic studies supported by molecular dynamic simulations combined with Molecular Mechanics-Poisson Boltzmann Surface Area calculations identified the residues important for binding in the A3R orthosteric site. We demonstrate that K18, which contains a 3-(dichlorophenyl)-isoxazole group connected through carbonyloxycarboximidamide fragment with a 1,3-thiazole ring, is a specific A3R (< 1 µM) competitive antagonist. Finally, we introduce a model that enables estimates of the equilibrium binding affinity for rapidly disassociating compounds from real-time fluorescent ligand-binding studies. These results demonstrate the pharmacological characterisation of a selective competitive A3R antagonist and the description of its orthosteric binding mode. Our findings may provide new insights for drug discovery.

Discovery of Novel Adenosine Receptor Antagonists through a Combined Structure- and Ligand-Based Approach Followed by Molecular Dynamics Investigation of Ligand Binding Mode.

An intense effort is made by pharmaceutical and academic research laboratories to identify and develop selective antagonists for each adenosine receptor (AR) subtype as potential clinical candidates for "soft" treatment of various diseases. Crystal structures of subtypes A2A and A1ARs offer exciting opportunities for structure-based drug design. In the first part of the present work, Maybridge HitFinder library of 14400 compounds was utilized to apply a combination of structure-based against the crystal structure of A2AAR and ligand-based methodologies. The docking poses were rescored by CHARMM energy minimization and calculation of the desolvation energy using Poisson-Boltzmann equation electrostatics. Out of the eight selected and tested compounds, five were found positive hits (63% success). Although the project was initially focused on targeting A2AAR, the identified antagonists exhibited low micromolar or micromolar affinity against A2A/A3, ARs, or A3AR, respectively. Based on these results, 19 compounds characterized by novel chemotypes were purchased and tested. Sixteen of them were identified as AR antagonists with affinity toward combinations of the AR family isoforms (A2A/A3, A1/A3, A1/A2A/A3, and A3). The second part of this work involves the performance of hundreds of molecular dynamics (MD) simulations of complexes between the ARs and a total of 27 ligands to resolve the binding interactions of the active compounds, which were not achieved by docking calculations alone. This computational work allowed the prediction of stable and unstable complexes which agree with the experimental results of potent and inactive compounds, respectively. Of particular interest is that the 2-amino-thiophene-3-carboxamides, 3-acylamino-5-aryl-thiophene-2-carboxamides, and carbonyloxycarboximidamide derivatives were found to be selective and possess a micromolar to low micromolar affinity for the A3 receptor.

Development of a Predictive Pharmacophore Model and a 3D-QSAR Study for an in silico Screening of New Potent Bcr-Abl Kinase Inhibitors.

Chronic myelogenous leukemia (CML) is a myeloproliferative disorder, characterized, in most cases, by the presence of the Bcr-Abl fusion oncogene. Bcr-Abl is a constitutively active tyrosine kinase that is responsible for the malignant transformation. Targeting the Bcr-Abl kinase is an attractive treatment strategy for CML. First and second generation Bcr-Abl inhibitors have focused on targeting the ATP-binding domain of the kinase. Mutations in that region are relatively resistant to drug manipulation. Therefore, non-ATP-competitive agents have been recently developed and tested. In the present study, in an attempt to aid the design of new chemotypes with enhanced cytotoxicity against K562 cells, 3D pharmacophore models were generated and 3D-QSAR CoMFA and CoMSIA studies were carried out on the 33 novel Abl kinase inhibitors (E)-α-benzylthio chalcones synthesized by Reddy et al. A five-point pharmacophore with a hydrogen bond acceptor, two hydrophobic groups and two aromatic rings as pharmacophore features, and a statistically significant 3D-QSAR model with excellent predictive power were developed. The pharmacophore model was also used for alignment of the 33 compounds in a CoMFA/CoMSIA analysis. The contour maps of the fields of CoMFA and CoMSIA models were utilized to provide structural insight into how these molecules promote their toxicity. The possibility of using this model for the design of drugs for the treatment of β-thalassemia and sickle cell disease (SCD), since several Bcr-Abl inhibitors are able to promote erythroid differentiation and γ-globin expression in CML cell lines and primary erythroid cells is discussed.

Searching for Novel Janus Kinase-2 Inhibitors Using a Combination of Pharmacophore Modeling, 3D-QSAR Studies and Virtual Screening.

The Janus kinases (JAKs) play a pivotal role in cytokine receptor signaling pathways via activation of downstream signal transducers and activators of transcription (STAT) pathway. Intracellular pathways that include JAKs are critical to immune cell activation and pro-inflammatory cytokine production. Selective inhibitors of JAKs are potentially disease-modifying anti-inflammatory drugs for the treatment of rheumatoid arthritis (RA). Each of the four members of the JAK family plays an individual role in the oncogenesis of the immune system, and therefore, the development of potent and specific inhibitors for each member is needed. Although there is a high sequence homology and structural identity of JAK1 and JAK2, such as a very similar binding mode of inhibitors at the ATPbinding site of enzymes, obvious differences surrounding the JAK1 and JAK2 ATP-binding sites provide a platform for the rational design of JAK2- and JAK1-specific inhibitors. In the present study, a dataset of 33 compounds characterized by a common scaffold of 2-amino-[1,2,4]triazolo[1,5-α]pyridine with well-defined in vitro activity values was computationally explored. Most of the compounds included in the dataset had higher ligand efficiency against JAK2 than JAK1. To improve further the selectivity of these triazolopyridines, Common Pharmacophore Hypotheses (CPHs) were generated and 3D-QSAR studies were carried out on them, in order to comprehend on the molecular features responsible for their selectivity. The proposed computational approach was applied in order to perform an in silico database virtual screening study with the aim to discover novel potent and selective JAK2 inhibitors.

Searching for anthranilic acid-based thumb pocket 2 HCV NS5B polymerase inhibitors through a combination of molecular docking, 3D-QSAR and virtual screening.

A combination of the following computational methods: (i) molecular docking, (ii) 3-D Quantitative Structure Activity Relationship Comparative Molecular Field Analysis (3D-QSAR CoMFA), (iii) similarity search and (iv) virtual screening using PubChem database was applied to identify new anthranilic acid-based inhibitors of hepatitis C virus (HCV) replication. A number of known inhibitors were initially docked into the "Thumb Pocket 2" allosteric site of the crystal structure of the enzyme HCV RNA-dependent RNA polymerase (NS5B GT1b). Then, the CoMFA fields were generated through a receptor-based alignment of docking poses to build a validated and stable 3D-QSAR CoMFA model. The proposed model can be first utilized to get insight into the molecular features that promote bioactivity, and then within a virtual screening procedure, it can be used to estimate the activity of novel potential bioactive compounds prior to their synthesis and biological tests.