Three Homology Models of PAR2 Derived from Different Templates: Application to Antagonist Discovery

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

BIO-GATS: A Tool for Automated GPCR Template Selection Through a Biophysical Approach for Homology Modeling

G protein-coupled receptors (GPCRs) are the largest family of membrane proteins with more than 800 members. GPCRs are involved in numerous physiological functions within the human body and are the target of more than 30% of the United States Food and Drug Administration (FDA) approved drugs. At present, over 400 experimental GPCR structures are available in the Protein Data Bank (PDB) representing 76 unique receptors. The absence of an experimental structure for the majority of GPCRs demand homology models for structure-based drug discovery workflows. The generation of good homology models requires appropriate templates. The commonly used methods for template selection are based on sequence identity. However, there exists low sequence identity among the GPCRs. Sequences with similar patterns of hydrophobic residues are often structural homologs, even with low sequence identity. Extending this, we propose a biophysical approach for template selection based principally on hydrophobicity correspondence between the target and the template. Our approach takes into consideration other relevant parameters, including resolution, similarity within the orthosteric binding pocket of GPCRs, and structure completeness, for template selection. The proposed method was implemented in the form of a free tool called Bio-GATS, to provide the user with easy selection of the appropriate template for a query GPCR sequence. Bio-GATS was successfully validated with recent published benchmarking datasets. An application to an olfactory receptor to select an appropriate template has also been provided as a case study.

Drug repositioning as an effective therapy for protease-activated receptor 2 inhibition

Proteinase-activated receptor 2 (PAR-2) is a G protein-coupled receptor activated by both trypsin and a specific agonist peptide, SLIGKV-NH2. It has been linked to various pathologies, including pain and inflammation. Several peptide and peptidomimetic agonizts for PAR-2 have been developed exhibiting high potency and efficacy. However, the number of PAR-2 antagonists is smaller. We screened the Food and Drug Administration library of approved compounds to retrieve novel antagonists for repositioning in the PAR-2 structure. The most efficacious compound bicalutamide bound to the PAR-2 binding groove near the extracellular domain as observed in the in silico studies. Further, it showed reduced Ca²⁺ release in trypsin activated cells in a dose-dependent manner. Hence, bicalutamide is a novel and potent PAR-2 antagonist which could be therapeutically useful in blocking multiple pathways diverging from PAR-2 signaling. Further, the novel scaffold of bicalutamide represents a new molecular structure for PAR-2 antagonism and can serve as a basis for further drug development.

Structural Characterization of Agonist Binding to Protease-Activated Receptor 2 through Mutagenesis and Computational Modeling

Protease-activated receptor 2 (PAR2) is a G-protein-coupled receptor that is activated by proteolytic cleavage of its N-terminus. The unmasked N-terminal peptide then binds to the transmembrane bundle, leading to activation of intracellular signaling pathways associated with inflammation and cancer. Recently determined crystal structures have revealed binding sites of PAR2 antagonists, but the binding mode of the peptide agonist remains unknown. In order to generate a model of PAR2 in complex with peptide SLIGKV, corresponding to the trypsin-exposed tethered ligand, the orthosteric binding site was probed by iterative combinations of receptor mutagenesis, agonist ligand modifications, and data-driven structural modeling. Flexible-receptor docking identified a conserved binding mode for agonists related to the endogenous ligand that was consistent with the experimental data and allowed synthesis of a novel peptide (1-benzyl-1H[1,2,3]triazole-4-yl-LIGKV) with functional potency higher than that of SLIGKV. The final model may be used to understand the structural basis of PAR2 activation and in virtual screens to identify novel agonists and competitive antagonists. The combined experimental and computational approach to characterize agonist binding to PAR2 can be extended to study the many other G protein-coupled receptors that recognize peptides or proteins.

PAR2 Modulators Derived from GB88

PAR2 antagonists have potential for treating inflammatory, respiratory, gastrointestinal, neurological, and metabolic disorders, but few antagonists are known. Derivatives of GB88 (3) suggest that all four of its components bind at distinct PAR2 sites with the isoxazole, cyclohexylalanine, and isoleucine determining affinity and selectivity, while the C-terminal substituent determines agonist/antagonist function. Here we report structurally similar PAR2 ligands with opposing functions (agonist vs antagonist) upon binding to PAR2. A biased ligand AY117 (65) was found to antagonize calcium release induced by PAR2 agonists trypsin and hexapeptide 2f-LIGRLO-NH₂ (IC₅₀ 2.2 and 0.7 μM, HT29 cells), but it was a selective PAR2 agonist in inhibiting cAMP stimulation and activating ERK1/2 phosphorylation. It showed anti-inflammatory properties both in vitro and in vivo.