The primary goal of rational drug discovery is the identification of selective ligands which act on single or multiple drug targets to achieve the desired clinical outcome through the exploration of total chemical space. To identify such desired compounds, computational approaches are necessary in predicting their drug-like properties. G Protein-Coupled Receptors (GPCRs) represent one of the largest and most important integral membrane protein families. These receptors serve as increasingly attractive drug targets due to their relevance in the treatment of various diseases, such as inflammatory disorders, metabolic imbalances, cardiac disorders, cancer, monogenic disorders, etc. In the last decade, multitudes of three-dimensional (3D) structures were solved for diverse GPCRs, thus referring to this period as the “golden age for GPCR structural biology.” Moreover, accumulation of data about the chemical properties of GPCR ligands has garnered much interest toward the exploration of GPCR chemical space. Due to the steady increase in the structural, ligand, and functional data of GPCRs, several cheminformatics approaches have been implemented in its drug discovery pipeline. In this review, we mainly focus on the cheminformatics-based paradigms in GPCR drug discovery. We provide a comprehensive view on the ligand– and structure-based cheminformatics approaches which are best illustrated via GPCR case studies. Furthermore, an appropriate combination of ligand-based knowledge with structure-based ones, i.e., integrated approach, which is emerging as a promising strategy for cheminformatics-based GPCR drug design is also discussed.

Janus kinase 1 (JAK1) plays a key role in most cytokine-mediated inflammatory and autoimmune responses through JAK/STAT signaling; thus, JAK1 inhibition is a promising therapeutic strategy for several diseases. Analysis of the binding modes of current JAK inhibitors to JAK isoforms allowed the design of N-alkyl-substituted 1-H-pyrrolo­[2,3-b] pyridine carboxamide as a JAK1-selective scaffold, and the synthesis of various methyl amide derivatives provided 4-((cis-1-(4-chlorobenzyl)-2-methylpiperidin-4-yl)­amino)-N-methyl-1H-pyrrolo­[2,3-b]­pyridine-5-carboxamide (31g) as a potent JAK1-selective inhibitor. In particular, the (S,S)-enantiomer of 31g (38a) exhibited excellent potency for JAK1 and selectivity over JAK2, JAK3, and TYK2. On investigating the effect of 31g on hepatic fibrosis, it was found that it reduces the proliferation and fibrogenic gene expression of TGF-β-induced hepatic stellate cells (HSCs). Specifically, 31g significantly inhibited TGF-β-induced migration of HSCs at 0.25 μM in wound-healing assays.

Glutaminyl cyclase (QC) has been implicated in the formation of toxic amyloid plaques by generating the N-terminal pyroglutamate of β-amyloid peptides (pGlu-Aβ) and thus may participate in the pathogenesis of Alzheimer's disease (AD). We designed a library of glutamyl cyclase (QC) inhibitors based on the proposed binding mode of the preferred substrate, Aβ. An in vitro structure-activity relationship study identified several excellent QC inhibitors demonstrating 5- to 40-fold increases in potency compared to a known QC inhibitor. When tested in mouse models of AD, compound 212 significantly reduced the brain concentrations of pyroform Aβ and total Aβ and restored cognitive functions. This potent Aβ-lowering effect was achieved by incorporating an additional binding region into our previously established pharmacophoric model, resulting in strong interactions with the carboxylate group of Glu327 in the QC binding site. Our study offers useful insights in designing novel QC inhibitors as a potential treatment option for AD.

Potent and selective A3 adenosine receptor (AR) agonists were identified by replacement of 4′-oxo- or 4′-thionucleosides with bioisosteric selenium. Unlike previous agonists, 4′-seleno analogues preferred a glycosidic syn conformation and a South sugar puckering, as shown in the X-ray crystal structure of 5′-N-methylcarbamoyl derivative 3p. Among compounds tested, N6-3-iodobenzyl analogue 3d was found to be the most potent A3AR full agonist (Ki = 0.57 nM), which was ≥ 800- and 1900-fold selective for A1 and A2AARs, respectively. In the N6-cycloalkyl series, 2-Cl analogues generally exhibited better hA3AR affinity than 2-H analogues, while 2-H > 2-Cl in the N6-3-halobenzyl series. The N7-isomers 3t and 3u were much weaker in binding than corresponding N9-isomers, but compound 3t lacked A3AR activation, appearing to be a weak antagonist. 2-Cl-N6-3-iodobenzyl analogue 3p inhibited chemoattractant-induced migration of microglia/monocytes without inducing cell death at ≤ 50 μM. This suggests the potential for the development of 4′-selenonucleoside A3AR agonists as novel anti-stroke agents.

A series of 1-substituted 3-(t-butyl/trifluoromethyl)pyrazole C-region analogues of 2-(3-fluoro-4-methylsulfonamidophenyl)propanamides were investigated for hTRPV1 antagonism. The structure activity relationship indicated that the 3-chlorophenyl group at the 1-position of pyrazole was the optimized hydrophobic group for antagonistic potency and the activity was stereospecific to the S-configuration, providing exceptionally potent antagonists 13S and 16S with K=0.1nM. Particularly significant, 13S exhibited antagonism selective for capsaicin and NADA and not for low pH or elevated temperature. Both compounds also proved to be very potent antagonists for rTRPV1, blocking in vivo the hypothermic action of capsaicin, consistent with their in vitro mechanism. The docking study of compounds 13S and 16S in our hTRPV1 homology model indicated that the binding modes differed somewhat, with that of 13S more closely resembling that of GRT12360.