Ligand and structure-based methodologies for the prediction of the activity of G protein-coupled receptor ligands

Accurate in silico models for the quantitative prediction of the activity of G protein-coupled receptor (GPCR) ligands would greatly facilitate the process of drug discovery and development. Several methodologies have been developed based on the properties of the ligands, the direct study of the receptor-ligand interactions, or a combination of both approaches. Ligand-based three-dimensional quantitative structure-activity relationships (3D-QSAR) techniques, not requiring knowledge of the receptor structure, have been historically the first to be applied to the prediction of the activity of GPCR ligands. They are generally endowed with robustness and good ranking ability; however they are highly dependent on training sets. Structure-based techniques generally do not provide the level of accuracy necessary to yield meaningful rankings when applied to GPCR homology models. However, they are essentially independent from training sets and have a sufficient level of accuracy to allow an effective discrimination between binders and nonbinders, thus qualifying as viable lead discovery tools. The combination of ligand and structure-based methodologies in the form of receptor-based 3D-QSAR and ligand and structure-based consensus models results in robust and accurate quantitative predictions. The contribution of the structure-based component to these combined approaches is expected to become more substantial and effective in the future, as more sophisticated scoring functions are developed and more detailed structural information on GPCRs is gathered.

Mechanism of Activation of a G Protein-coupled Receptor, the Human Cholecystokinin-2 Receptor

G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.

Modeled structure of a G-protein-coupled receptor:The cholecystokinin-1 receptor

The Cholecystokinin-1 receptor (CCK1R) mediates actions of CCK in areas of the central nervous system and of the gut. It is a potential target to treat a number of diseases. As for all G-protein-coupled receptors, docking of ligands into modeled CCK1R binding site should greatly help to understand intrinsic mechanisms of activation. Here, we describe the procedure we used to progressively build a structural model for the CCK1R, to integrated, and on the basis of site-directed mutagenesis data on its binding site. Reliability of the CCK1R model was confirmed by interaction networks that involved conserved and functionally crucial motifs in G-protein-coupled receptors, such as Glu/Asp-Arg-Tyr and Asn-Pro-Xaa-Xaa-Tyr motifs. In addition, the 3-D structure of CCK1R-bound CCK resembled that determined by NMR in a lipid environment. The derived computational model was also used for revealing binding modes of several nonpeptide ligands and for rationalizing ligand structure-activity relationships known from experiments. Our findings indeed support that our "validated CCK1R model" could be used to study the intrinsic mechanism of CCK1R activation and design new ligands.

A Molecular Mechanism for Sequential Activation of a G Protein-Coupled Receptor

Ligands targeting G protein-coupled receptors (GPCRs) are currently classified as either orthosteric, allosteric, or dualsteric/bitopic. Here, we introduce a new pharmacological concept for GPCR functional modulation: sequential receptor activation. A hallmark feature of this is a stepwise ligand binding mode with transient activation of a first receptor site followed by sustained activation of a second topographically distinct site. We identify 4-CMTB (2-(4-chlorophenyl)-3-methyl-N-(thiazol-2-yl)butanamide), previously classified as a pure allosteric agonist of the free fatty acid receptor 2, as the first sequential activator and corroborate its two-step activation in living cells by tracking integrated responses with innovative label-free biosensors that visualize multiple signaling inputs in real time. We validate this unique pharmacology with traditional cellular readouts, including mutational and pharmacological perturbations along with computational methods, and propose a kinetic model applicable to the analysis of sequential receptor activation. We envision this form of dynamic agonism as a common principle of nature to spatiotemporally encode cellular information.

Rhodopsin and the Others: A Historical Perspective on Structural Studies of G Protein-Coupled Receptors

The role of rhodopsin as a structural prototype for the study of the whole superfamily of G protein-coupled receptors (GPCRs) is reviewed in an historical perspective. Discovered at the end of the nineteenth century, fully sequenced since the early 1980s, and with direct three-dimensional information available since the 1990s, rhodopsin has served as a platform to gather indirect information on the structure of the other superfamily members. Recent breakthroughs have elicited the solution of the structures of additional receptors, namely the beta 1- and beta 2-adrenergic receptors and the A(2A) adenosine receptor, now providing an opportunity to gauge the accuracy of homology modeling and molecular docking techniques and to perfect the computational protocol. Notably, in coordination with the solution of the structure of the A(2A) adenosine receptor, the first "critical assessment of GPCR structural modeling and docking" has been organized, the results of which highlighted that the construction of accurate models, although challenging, is certainly achievable. The docking of the ligands and the scoring of the poses clearly emerged as the most difficult components. A further goal in the field is certainly to derive the structure of receptors in their signaling state, possibly in complex with agonists. These advances, coupled with the introduction of more sophisticated modeling algorithms and the increase in computer power, raise the expectation for a substantial boost of the robustness and accuracy of computer-aided drug discovery techniques in the coming years.

The repertoire of G protein-coupled receptors in the human parasite Schistosoma mansoni and the model organism Schmidtea mediterranea.

Background

G protein-coupled receptors (GPCRs) constitute one of the largest groupings of eukaryotic proteins, and represent a particularly lucrative set of pharmaceutical targets. They play an important role in eukaryotic signal transduction and physiology, mediating cellular responses to a diverse range of extracellular stimuli. The phylum Platyhelminthes is of considerable medical and biological importance, housing major pathogens as well as established model organisms. The recent availability of genomic data for the human blood fluke Schistosoma mansoni and the model planarian Schmidtea mediterranea paves the way for the first comprehensive effort to identify and analyze GPCRs in this important phylum.

Results

Application of a novel transmembrane-oriented approach to receptor mining led to the discovery of 117 S. mansoni GPCRs, representing all of the major families; 105 Rhodopsin, 2 Glutamate, 3 Adhesion, 2 Secretin and 5 Frizzled. Similarly, 418 Rhodopsin, 9 Glutamate, 21 Adhesion, 1 Secretin and 11 Frizzled S. mediterranea receptors were identified. Among these, we report the identification of novel receptor groupings, including a large and highly-diverged Platyhelminth-specific Rhodopsin subfamily, a planarian-specific Adhesion-like family, and atypical Glutamate-like receptors. Phylogenetic analysis was carried out following extensive gene curation. Support vector machines (SVMs) were trained and used for ligand-based classification of full-length Rhodopsin GPCRs, complementing phylogenetic and homology-based classification.

Conclusions

Genome-wide investigation of GPCRs in two platyhelminth genomes reveals an extensive and complex receptor signaling repertoire with many unique features. This work provides important sequence and functional leads for understanding basic flatworm receptor biology, and sheds light on a lucrative set of anthelmintic drug targets.

Unraveling the Structure and Function of G Protein-Coupled Receptors Through NMR Spectroscopy

G protein-coupled receptors (GPCRs) are a large superfamily of signaling proteins expressed on the plasma membrane. They are involved in a wide range of physiological processes and, therefore, are exploited as drug targets in a multitude of therapeutic areas. In this extent, knowledge of structural and functional properties of GPCRs may greatly facilitate rational design of modulator compounds. Solution and solid-state nuclear magnetic resonance (NMR) spectroscopy represents a powerful method to gather atomistic insights into protein structure and dynamics. In spite of the difficulties inherent the solution of the structure of membrane proteins through NMR, these methods have been successfully applied, sometimes in combination with molecular modeling, to the determination of the structure of GPCR fragments, the mapping of receptor-ligand interactions, and the study of the conformational changes associated with the activation of the receptors. In this review, we provide a summary of the NMR contributions to the study of the structure and function of GPCRs, also in light of the published crystal structures.

Modeled structure of the whole regulator G-protein signaling-2

There is an increasing interest towards the mechanism by which regulators of G-protein signaling regulate signals of G-protein-coupled receptors. RGS2 is a regulator of Gq protein signaling (RGS), the N-terminal region of which is known to contain determinants for G protein-coupled receptor recognition, but its structure is still unknown. To understand the molecular basis for this recognition, the three-dimensional model of RGS2, including N-terminal region and RGS box, was modeled. For this, RGS4 box structure and data from circular dichroism study of RGS2 N-terminal region were used. Then, membrane-targeting activity of the RGS2 amphipathic helix contained in the N-terminal region was investigated. Furthermore, in cellulo study provided first evidence that an internal sequence within the N-terminal region of RGS2 is involved in RGS2 regulation of cholecystokinin receptor-2 signal. RGS2 modeled structure can now serve to study molecular recognition of RGS2 by signaling molecules. © 2006 Elsevier Inc. All rights reserved.

Evidence for a Direct and Functional Interaction between the Regulators of G Protein Signaling-2 and Phosphorylated C Terminus of Cholecystokinin-2 Receptor

Signaling of G protein-coupled receptors (GPCRs) is regulated by different mechanisms. One of these involves regulators of G protein signaling (RGS), which are diverse and multifunctional proteins that bind to active G alpha subunits of G proteins and act as GTPase-activating proteins. Little is known about the molecular mechanisms that govern the selective use of RGS proteins in living cells. We first demonstrated that CCK2R-mediated inositol phosphate production, known to be G(q-)dependent, is more sensitive to RGS2 than to RGS4 and is insensitive to RGS8. Both basal and agonist-stimulated activities of the CCK2R are regulated by RGS2. By combining biochemical functional, and in silico structural approaches, we demonstrate that a direct and functional interaction occurs between RGS2 and agonist-stimulated cholecystokinin receptor-2 (CCK2R) and identified the precise residues involved: phosphorylated Ser434 and Thr439 located in the C-terminal tail of CCK2R and Lys62, Lys63, and Gln67, located in the N-terminal domain of RGS2. These findings confirm previous reports that RGS proteins can interact with GPCRs to modulate their signaling and provide a molecular basis for RGS2 recognition by the CCK2R.

Discovery of novel Agonists and antagonists of the free fatty acid receptor 1 (FFAR1) using virtual screening

The G-protein-coupled receptor free fatty acid receptor 1 (FFAR1), previously named GPR40, is a possible novel target for the treatment of type 2 diabetes. In an attempt to identify new ligands for this receptor, we performed virtual screening (VS) based on two-dimensional (2D) similarity, three-dimensional (3D) pharmacophore searches, and docking studies by using the structure of known agonists and our model of the ligand binding site, which was validated by mutagenesis. VS of a database of 2.6 million compounds followed by extraction of structural neighbors of functionally confirmed hits resulted in identification of 15 compounds active at FFAR1 either as full agonists, partial agonists, or pure antagonists. Site-directed mutagenesis and docking studies revealed different patterns of ligand-receptor interactions and provided important information on the role of specific amino acids in binding and activation of FFAR1.

P2Y(1) antagonists: Combining receptor-based modeling and QSAR for a quantitative prediction of the biological activity based on consensus scoring

P2Y(1) is an ADP-activated G protein-coupled receptor (GPCR). Its antagonists impede platelet aggregation in vivo and are potential antithrombotic agents. Combining ligand and structure-based modeling we generated a consensus model (LIST-CM) correlating antagonist structures with their potencies. We docked 45 antagonists into our rhodopsin-based human P2Y(1) homology model and calculated docking scores and free binding energies with the Linear Interaction Energy (LIE) method in continuum-solvent. The resulting alignment was also used to build QSAR based on CoMFA, CoMSIA, and molecular descriptors. To benefit from the strength of each technique and compensate for their limitations, we generated our LIST-CM with a PLS regression based on the predictions of each methodology. A test set featuring untested substituents was synthesized and assayed in inhibition of 2-MeSADP-stimulated PLC activity and in radioligand binding. LIST-CM outperformed internal and external predictivity of any individual model to predict accurately the potency of 75% of the test set.

Identification of Galanin and Its Receptor GalR1 as Novel Determinants of Resistance to Chemotherapy and Potential Biomarkers in Colorectal Cancer

Purpose: A major factor limiting the effective clinical management of colorectal cancer (CRC) is resistance to chemotherapy. Therefore, the identification of novel, therapeutically targetable mediators of resistance is vital.Experimental design: We used a CRC disease-focused microarray platform to transcriptionally profile chemotherapy-responsive and nonresponsive pretreatment metastatic CRC liver biopsies and in vitro samples, both sensitive and resistant to clinically relevant chemotherapeutic drugs (5-FU and oxaliplatin). Pathway and gene set enrichment analyses identified candidate genes within key pathways mediating drug resistance. Functional RNAi screening identified regulators of drug resistance.

Results: Mitogen-activated protein kinase signaling, focal adhesion, cell cycle, insulin signaling, and apoptosis were identified as key pathways involved in mediating drug resistance. The G-protein-coupled receptor galanin receptor 1 (GalR1) was identified as a novel regulator of drug resistance. Notably, silencing either GalR1 or its ligand galanin induced apoptosis in drug-sensitive and resistant cell lines and synergistically enhanced the effects of chemotherapy. Mechanistically, GalR1/galanin silencing resulted in downregulation of the endogenous caspase-8 inhibitor FLIP(L), resulting in induction of caspase-8-dependent apoptosis. Galanin mRNA was found to be overexpressed in colorectal tumors, and importantly, high galanin expression correlated with poor disease-free survival of patients with early-stage CRC.

Conclusion: This study shows the power of systems biology approaches to identify key pathways and genes that are functionally involved in mediating chemotherapy resistance. Moreover, we have identified a novel role for the GalR1/galanin receptor-ligand axis in chemoresistance, providing evidence to support its further evaluation as a potential therapeutic target and biomarker in CRC. Clin Cancer Res; 18(19); 5412–26. © 2012 AACR.