Progress in structure based drug design for G protein-coupled receptors

Structure of the human κ-opioid receptor in complex with JDTic

Opioid receptors (ORs) mediate the actions of endogenous and exogenous opioids for many essential physiological processes including regulation of pain, respiratory drive, mood, and, in the case of κ-opioid receptors (KOR), dysphoria and psychotomimesis. Here we report the crystal structure of the human KOR (hKOR) in complex with the selective antagonist JDTic, arranged in parallel-dimers, at 2.9 angstrom resolution. The structure reveals important features of the ligand binding pocket that contribute to JDTic’s high affinity and subtype-selectivity for hKOR. Modeling of other important KOR-selective ligands, including the morphinan-derived antagonists nor-BNI and GNTI, and the diterpene agonist salvinorin A analog RB-64, reveals both common and distinct features for binding these diverse chemotypes. Analysis of site-directed mutagenesis and ligand structure-activity relationships confirms the interactions observed in the crystal structure, thereby providing a molecular explanation for hKOR subtype-selectivity along with insight essential for the design of hKOR compounds with new pharmacological properties.

Identification of novel adenosine A(2A) receptor antagonists by virtual screening

Virtual screening was performed against experimentally enabled homology models of the adenosine A(2A) receptor, identifying a diverse range of ligand efficient antagonists (hit rate 9%). By use of ligand docking and Biophysical Mapping (BPM), hits 1 and 5 were optimized to potent and selective lead molecules (11-13 from 5, pK(I) = 7.5-8.5, 13- to >100-fold selective versus adenosine A(1); 14-16 from 1, pK(I) = 7.9-9.0, 19- to 59-fold selective).

Fluorescent derivatives of AC-42 to probe bitopic orthosteric/allosteric binding mechanisms on muscarinic M1 receptors

Two fluorescent derivatives of the M1 muscarinic selective agonist AC-42 were synthesized by coupling the lissamine rhodamine B fluorophore (in ortho and para positions) to AC42-NH(2). This precursor, prepared according to an original seven-step procedure, was included in the study together with the LRB fluorophore (alone or linked to an alkyl chain). All these compounds are antagonists, but examination of their ability to inhibit or modulate orthosteric [(3)H]NMS binding revealed that para-LRB-AC42 shared several properties with AC-42. Carefully designed experiments allowed para-LRB-AC42 to be used as a FRET tracer on EGFP-fused M1 receptors. Under equilibrium binding conditions, orthosteric ligands, AC-42, and the allosteric modulator gallamine behaved as competitors of para-LRB-AC42 binding whereas other allosteric compounds such as WIN 51,708 and N-desmethylclozapine were noncompetitive inhibitors. Finally, molecular modeling studies focused on putative orthosteric/allosteric bitopic poses for AC-42 and para-LRB-AC42 in a 3D model of the human M1 receptor.

Emerging role of surface plasmon resonance in fragment-based drug discovery

Surface plasmon resonance (SPR) offers a method of biophysical fragment screening that is fast, efficient, cost effective and accurate. SPR is increasingly being adopted as a secondary assay to validate fragment hits. Recently, technical advances have resulted in the emergence of SPR as a primary screening methodology for fragment-based drug discovery. Moreover, SPR biosensor assays can be developed for a wide range of proteins, including membrane proteins, such as G-protein-coupled receptors. In this review, we discuss the advantages and limitations of SPR fragment screening including experimental consideration of reducing false positive and false negative rates to a minimum. We discuss how ligand efficiency can be used both as a method to eliminate false positives and to understand which fragments in a library may be a source of false negatives.

Discovery of 1,2,4-triazine derivatives as adenosine A(2A) antagonists using structure based drug design

Potent, ligand efficient, selective, and orally efficacious 1,2,4-triazine derivatives have been identified using structure based drug design approaches as antagonists of the adenosine A_2A receptor. The X-ray crystal structures of compounds 4e and 4g bound to the GPCR illustrate that the molecules bind deeply inside the orthosteric binding cavity. In vivo pharmacokinetic and efficacy data for compound 4k are presented, demonstrating the potential of this series of compounds for the treatment of Parkinson’s disease.

Cell-based and in-silico studies on the high intrinsic activity of two boron-containing salbutamol derivatives at the human β₂-adrenoceptor

Salbutamol is a well-known β(2) adrenoceptor (β(2)AR) partial agonist. We synthesized two boron-containing salbutamol derivatives (BCSDs) with greater potency and efficacy, compared to salbutamol, for inducing β(2)AR-mediated smooth-muscle relaxation in guinea-pig tracheal rings. However, the mechanism involved in this pharmacological effect remains unclear. In order to gain insight, we carried out binding and functional assays for BCSDs in HEK-293T cells transfected with the human β(2)AR (hβ(2)AR). The transfected hβ(2)AR showed similar affinity for BCSDs and salbutamol, but adenosine 3',5'-cyclic phosphate (cAMP) accumulation induced by both BCSDs was similar to that elicited by isoproterenol and greater than that induced by salbutamol. The boron-containing precursors (boric and phenylboronic acids, 100 μM) had no significant effect on salbutamol binding or salbutamol-induced cAMP accumulation. These experimental results are in agreement with theoretical docking simulations on lipid bilayer membrane-embedded hβ(2)AR structures. These receptors showed slightly higher affinity for BCSDs than for salbutamol. An essential change between putative active and inactive conformational states depended on the interaction of the tested ligands with the fifth, sixth and seventh transmembrane domains. Overall, these data suggest that BCSDs induce and stabilize conformational states of the hβ(2)AR that are highly capable of stimulating cAMP production.

Human Serotonin 5-HT2C G Protein-Coupled Receptor Homology Model from the β2 Adrenoceptor Structure: Ligand Docking and Mutagenesis Studies

Activation of the serotonin (5-hydroxytryptamine, 5-HT) 5HT_2C G protein-coupled receptor (GPCR) is proposed as novel pharmacotherapy for obesity and neuropsychiatric disorders. In contrast, activation of the 5-HT_2A and 5-HT_2B GPCRs is associated with untoward hallucinogenic and cardiopulmonary effects, respectively. There is no crystal structure available to guide design of 5-HT_2C receptor-specific ligands. For this reason, a homology model of the 5-HT_2C receptor was built based on the crystal structure of the human β₂ adrenoceptor GPCR to delineate molecular determinants of ligand-receptor interactions for drug design purposes. Computational and experimental studies were carried out to validate the model. Binding of N(CH₃)₂-PAT [(1R, 3S)-(-)-trans-1-phenyl-3-N,N-dimethylamino-1,2,3,4-tetrahydronaphthalene], a novel 5-HT_2C agonist/5-HT_2A/2B inverse agonist, and its secondary [NH(CH₃)-PAT] and primary (NH₂-PAT) amine analogs were studied at the 5-HT_2C wild type (WT) and D3.32A, S3.36A, and Y7.43A 5-HT_2C point-mutated receptors. Reference ligands included the tertiary amines lisuride and mesulergine and the primary amine 5-HT. Modeling results indicated that 5-HT_2C residues D3.32, S3.36, and Y7.43 play a role in ligand binding. Experimental ligand binding results with WT and point-mutated receptors confirmed the impact of D3.32, S3.36, and Y7.43 on ligand affinity.

Microscopic binding of M5 muscarinic acetylcholine receptor with antagonists by homology modeling, molecular docking, and molecular dynamics simulation

By performing homology modeling, molecular docking, and molecular dynamics (MD) simulations, we have developed three-dimensional (3D) structural models of the M5 muscarinic acetylcholine receptor (mAChR) and two complexes for M5 mAChR binding with antagonists SVT-40776 and solifenacin in the environment of lipid bilayer and solvent water. According to the simulated results, each of the antagonists is oriented horizontally in the binding pocket formed by transmembrane helices 2, 3, and 5-7. The cationic headgroup of each of the antagonists interacts with a negatively charged residue, Asp110, through electrostatic and hydrogen-bonding interactions. The simulated results also reveal some significant difference between the binding modes of SVT-40776 and solifenacin. In particular, SVT-40776 is persistently hydrogen bonded with the side chain of residue Tyr458, whereas solifenacin cannot form a similar hydrogen bond with residues around its carbonyl group. Such significant difference in the binding structures is consistent with the fact that SVT-40776 has a much higher binding affinity (K(d) = 0.4 nM) to M5 mAChR than that of solifenacin (K(d) = 31 nM) with the same reeptor. The calculated binding free energy change (-2.3 ± 0.3 kcal/mol) from solifenacin to SVT-40776 is in good agreement with the experimentally derived binding free energy change (-2.58 kcal/mol), suggesting that our modeled M5 mAChR structure and its complexes with the antagonists are reliable. The new structural insights obtained from this computational study are expected to stimulate further biochemical and pharmacological studies on the detailed structures of M5 and other subtypes of mAChRs.

Ligand and decoy sets for docking to G protein-coupled receptors

We compiled a G protein-coupled receptor (GPCR) ligand library (GLL) for 147 targets, selecting for each ligand 39 decoy molecules, collected in the GPCR Decoy Database (GDD). Decoys were chosen ensuring a ligand-decoy similarity of six physical properties, while enforcing ligand-decoy chemical dissimilarity. The performance in docking of the GDD was evaluated on 19 GPCRs, showing a marked decrease in enrichment compared to bias-uncorrected decoy sets. Both the GLL and GDD are freely available for the scientific community.

In silico screening for agonists and blockers of the β(2) adrenergic receptor: implications of inactive and activated state structures

Ten crystal structures of the β(2) adrenergic receptor have been published, reflecting different signaling states. Here, through controlled-docking experiments, we examined the implications of using inactive or activated structures on the in silico screening for agonists and blockers of the receptor. Specifically, we targeted the crystal structures solved in complex with carazolol (2RH1), the neutral antagonist alprenalol, the irreversible agonist FAUC50 (3PDS), and the full agonist BI-167017 (3P0G). Our results indicate that activated structures favor agonists over blockers, whereas inactive structures favor blockers over agonists. This tendency is more marked for activated than for inactive structures. Additionally, agonists tend to receive more favorable docking scores when docked at activated rather than inactive structures, while blockers do the opposite. Hence, the difference between the docking scores attained with an activated and an inactive structure is an excellent means for the classification of ligands into agonists and blockers as we determined through receiver operating characteristic curves and linear discriminant analysis. With respect to virtual screening, all structures prioritized well agonists and blockers over nonbinders. However, inactive structures worked better for blockers and activated structures worked better for agonists, respectively. Notably, the combination of individual docking experiments through receptor ensemble docking resulted in an excellent performance in the retrieval of both agonists and blockers. Finally, we demonstrated that the induced-fit docking of agonists is a viable way of modifying an inactive crystal structure and bias it toward the in silico recognition of agonists rather than blockers.

Crystal structure-based virtual screening for fragment-like ligands of the human histamine H(1) receptor

The recent crystal structure determinations of druggable class A G protein-coupled receptors (GPCRs) have opened up excellent opportunities in structure-based ligand discovery for this pharmaceutically important protein family. We have developed and validated a customized structure-based virtual fragment screening protocol against the recently determined human histamine H(1) receptor (H(1)R) crystal structure. The method combines molecular docking simulations with a protein-ligand interaction fingerprint (IFP) scoring method. The optimized in silico screening approach was successfully applied to identify a chemically diverse set of novel fragment-like (≤22 heavy atoms) H(1)R ligands with an exceptionally high hit rate of 73%. Of the 26 tested fragments, 19 compounds had affinities ranging from 10 μM to 6 nM. The current study shows the potential of in silico screening against GPCR crystal structures to explore novel, fragment-like GPCR ligand space.

GPCR agonist binding revealed by modeling and crystallography

Despite recent progress in structural coverage of the G-protein-coupled receptor (GPCR) family, high plasticity of these membrane proteins poses additional challenges for crystallographic studies of their complexes with different classes of ligands, especially agonists. The ability to predict computationally the binding of natural and clinically relevant agonists and corresponding changes in the receptor pocket, starting from inactive GPCR structures, is therefore of great interest for understanding GPCR biology and drug action. Comparison of computational models published in 2009 and 2010 with recently determined agonist-bound structures of β-adrenergic and adenosine A(2A) receptors reveals high accuracy of the predicted agonist binding poses (0.8 Å and 1.7 Å respectively) and receptor interactions. In the case of the β(2)AR, energy-based models with limited backbone flexibility have also allowed characterization of side-chain rotations and a finite backbone shift in the pocket region as determinants of full, partial or inverse agonism. Development of accurate models of agonist binding for other GPCRs will be instrumental for functional and pharmacological studies, complementing biochemical and crystallographic techniques.

Diversity and modularity of G protein-coupled receptor structures

G protein-coupled receptors (GPCRs) comprise the most 'prolific' family of cell membrane proteins, which share a general mechanism of signal transduction, but greatly vary in ligand recognition and function. Crystal structures are now available for rhodopsin, adrenergic, and adenosine receptors in both inactive and activated forms, as well as for chemokine, dopamine, and histamine receptors in inactive conformations. Here we review common structural features, outline the scope of structural diversity of GPCRs at different levels of homology, and briefly discuss the impact of the structures on drug discovery. Given the current set of GPCR crystal structures, a distinct modularity is now being observed between the extracellular (ligand-binding) and intracellular (signaling) regions. The rapidly expanding repertoire of GPCR structures provides a solid framework for experimental and molecular modeling studies, and helps to chart a roadmap for comprehensive structural coverage of the whole superfamily and an understanding of GPCR biological and therapeutic mechanisms.