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.

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.

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.

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.

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.

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.

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.

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.

Free fatty acid receptor 1 (GPR40) agonists containing spirocyclic periphery inspired by LY2881835

The free fatty acid receptor 1 (FFA1), a G protein-coupled receptor (GPCR) naturally activated by long-chain fatty acids is a novel target for the treatment of metabolic diseases. The basic amine spirocyclic periphery of Eli Lilly's drug candidate LY2881835 for treatment of type 2 diabetes mellitus (which reached phase I clinical trials) inspired a series of novel FFA1 agonists. These were designed to incorporate the 3-[4-(benzyloxy)phenyl]propanoic acid pharmacophore core decorated with a range of spirocyclic motifs. The latter were prepared via the Prins cyclization and subsequent modification of the 4-hydroxytetrahydropyran moiety in the Prins product. Here, we synthesize 19 compounds and test for FFA1 activity. Within this pilot set, a nanomolar potency (EC50=55nM) was reached. Four lead compounds (EC50 range 55-410nM) were characterized for aqueous solubility, metabolic stability, plasma protein binding and Caco-2 permeability. While some instability in the presence of mouse liver microsomes was noted, mouse pharmacokinetic profile of the compound having the best overall ADME properties was evaluated to reveal acceptable bioavailability (F=10.3%) and plasma levels achieved on oral administration.

Atomistic Insights into Rhodopsin Activation from a Dynamic Model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

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.