Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Computational Predictions of G Protein-Coupled Receptor Structures and Binding Sites

G protein-coupled receptors (GPCRs) are the largest family of proteins within the human genome. They consist of seven transmembrane (TM) helices, with a N-terminal region of varying length and structure on the extracellular side, and a C-terminus on the intracellular side. GPCRs are involved in transmitting extracellular signals to cells, and as such are crucial drug targets. Designing pharmaceuticals to target GPCRs is greatly aided by full-atom structural information of the proteins. In particular, the TM region of GPCRs is where small molecule ligands (much more bioavailable than peptide ligands) typically bind to the receptors. In recent years nearly thirty distinct GPCR TM regions have been crystallized. However, there are more than 1,000 GPCRs, leaving the vast majority of GPCRs with limited structural information. Additionally, GPCRs are known to exist in a myriad of conformational states in the body, rendering the static x-ray crystal structures an incomplete reflection of GPCR structures. In order to obtain an ensemble of GPCR structures, we have developed the GEnSeMBLE procedure to rapidly sample a large number of variations of GPCR helix rotations and tilts. The lowest energy GEnSeMBLE structures are then docked to small molecule ligands and optimized. The GPCR family consists of five subfamilies with little to no sequence homology between them: class A, B1, B2, C, and Frizzled/Taste2. Almost all of the GPCR crystal structures have been of class A GPCRs, and much is known about their conserved interactions and binding sites. In this work we particularly focus on class B1 GPCRs, and aim to understand that family’s interactions and binding sites both to small molecules and their native peptide ligands. Specifically, we predict the full atom structure and peptide binding site of the glucagon-like peptide receptor and the TM region and small molecule binding sites for eight other class B1 GPCRs: CALRL, CRFR1, GIPR, GLR, PACR, PTH1R, VIPR1, and VIPR2. Our class B1 work reveals multiple conserved interactions across the B1 subfamily as well as a consistent small molecule binding site centrally located in the TM bundle. Both the interactions and the binding sites are distinct from those seen in the more well-characterized class A GPCRs, and as such our work provides a strong starting point for drug design targeting class B1 proteins. We also predict the full structure of CXCR4 bound to a small molecule, a class A GPCR that was not closely related to any of the class A GPCRs at the time of the work.

Estudio de una serie musteriense de La Quina (Charente, Francia): el ejemplo de la trinchera C (excavaciones DR. L. y G. Henri-Martin)

[ES] El análisis de una serie industrial correspondiente a los cinco niveles arqueológicos de la trinchera C de La Quina (Charente, Francia), proveniente de las excavaciones del Dr. L. Henri-Martin y de G. Henri-Martin, pone de manifiesto, junto a la caracterización del conjunto estratigráfico por la presencia de raederas, el progresivo incremento de los denticulados en el devenir temporal. Este hecho asociado a otros particulares fenómenos secundarios (en el alargamiento de las formas, en la tipología de las raederas, en la elaboración del retoque, ... ) traducen una complicación diacrónica creciente en la dinámica de este complejo musteroide, presagiando, quizás, su relativa afinidad con el proceso leptolitizante de las series industriales evolucionadas del Paleolítico medio.

Subcellular localization and characterization of G protein-coupled receptor homolog from lymphocystis disease virus isolated in China

G protein-coupled receptors (GPCRs) constitute a large superfamily involved in various types of signal transduction pathways, and play an important role in coordinating the activation and migration of leukocytes to sites of infection and inflammation. Viral GPCRs, on the other hand, can help the virus to escape from host immune surveillance and contribute to viral pathogenesis. Lymphocystis disease virus isolated in China (LCDV-C) contains a putative homolog of cellular GPCRs, LCDV-C GPCR. In this paper, LCDV-C GPCR was cloned, and the subcellular localization and characterization of GPCR protein were investigated in fish cells. LCDV-C GPCR encoded a 325-amino acid peptide, containing a typical seven-transmembrane domain characteristic of the chemokine receptors and a conserved DRY motif that is usually essential for receptor activation. Transient transfection of GPCR-EGFP in fathead minnow (FHM) cells and epithelioma papulosum cyprini (EPC) cells indicated that LCDV-C GPCR was expressed abundantly in both the cytoplasm and nucleoplasm. Transient overexpression of GPCR in these two cells cannot induce obvious apoptosis. FHM cells stably expressing GPCR showed enhanced cell proliferation and significant anchorage-independent growth. The effects of GPCR protein on external apoptotic stimuli were examined. Few apoptotic bodies were observed in cells expressing GPCR treated with actinomycin D (ActD). Quantitative analysis of apoptotic cells indicated that a considerable decrease in the apoptotic fraction of cells expressing GPCR, compared with. the control cells, was detected after exposure to ActD and cycloheximide. These data suggest that LCDV-C GPCR may inhibit apoptosis as part of its potential mechanism in mediating cellular transformation.

Trafficking of G protein-coupled receptors.

G protein-coupled receptors (GPCRs) play an integral role in the signal transduction of an enormous array of biological phenomena, thereby serving to modulate at a molecular level almost all components of human biology. This role is nowhere more evident than in cardiovascular biology, where GPCRs regulate such core measures of cardiovascular function as heart rate, contractility, and vascular tone. GPCR/ligand interaction initiates signal transduction cascades, and requires the presence of the receptor at the plasma membrane. Plasma membrane localization is in turn a function of the delivery of a receptor to and removal from the cell surface, a concept defined most broadly as receptor trafficking. This review illuminates our current view of GPCR trafficking, particularly within the cardiovascular system, as well as highlights the recent and provocative finding that components of the GPCR trafficking machinery can facilitate GPCR signaling independent of G protein activation.

Purification, crystallization and preliminary X-ray diffraction studies of a complex between G protein-coupled receptor kinase 2 and Gbeta1gamma2.

G protein-coupled receptor kinase 2 (GRK2) phosphorylates activated G protein-coupled receptors (GPCRs), which ultimately leads to their desensitization and/or downregulation. The enzyme is recruited to the plasma membrane via the interaction of its carboxyl-terminal pleckstrin-homology (PH) domain with the beta and gamma subunits of heterotrimeric G proteins (Gbetagamma). An improved purification scheme for GRK2 has been developed, conditions under which GRK2 forms a complex with Gbeta(1)gamma(2) have been determined and the complex has been crystallized in CHAPS detergent micelles. Crystals of the GRK2-Gbetagamma complex belong to space group C2 and have unit-cell parameters a = 187.0, b = 72.1, c = 122.0 A, beta = 115.2 degrees. A complete data set has been collected to 3.2 A resolution with Cu Kalpha radiation.

The role of β-arrestins in the termination and transduction of G-protein-coupled receptor signals

β-arrestins are versatile adapter proteins that form complexes with most G-protein-coupled receptors (GPCRs) following agonist binding and phosphorylation of receptors by G-protein-coupled receptor kinases (GRKs). They play a central role in the interrelated processes of homologous desensitization and GPCR sequestration, which lead to the termination of G protein activation. β-arrestin binding to GPCRs both uncouples receptors from heterotrimeric G proteins and targets them to clathrincoated pits for endocytosis. Recent data suggest that β-arrestins also function as GPCR signal transducers. They can form complexes with several signaling proteins, including Src family tyrosine kinases and components of the ERK1/2 and JNK3 MAP kinase cascades. By recruiting these kinases to agonist-occupied GPCRs, β-arrestins confer distinct signaling activities upon the receptor. β-arrestin-Src complexes have been proposed to modulate GPCR endocytosis, to trigger ERK1/2 activation and to mediate neutrophil degranulation. By acting as scaffolds for the ERK1/2 and JNK3 cascades, β-arrestins both facilitate GPCR-stimulated MAP kinase activation and target active MAP kinases to specific locations within the cell. Thus, their binding to GPCRs might initiate a second wave of signaling and represent a novel mechanism of GPCR signal transduction.

Receptor and G betagamma isoform-specific interactions with G protein-coupled receptor kinases.

The G protein-coupled receptor (GPCR) kinases (GRKs) phosphorylate and desensitize agonist-occupied GPCRs. GRK2-mediated receptor phosphorylation is preceded by the agonist-dependent membrane association of this enzyme. Previous in vitro studies with purified proteins have suggested that this translocation may be mediated by the recruitment of GRK2 to the plasma membrane by its interaction with the free betagamma subunits of heterotrimeric G proteins (G betagamma). Here we demonstrate that this mechanism operates in intact cells and that specificity is imparted by the selective interaction of discrete pools of G betagamma with receptors and GRKs. Treatment of Cos-7 cells transiently overexpressing GRK2 with a beta-receptor agonist promotes a 3-fold increase in plasma membrane-associated GRK2. This translocation of GRK2 is inhibited by the carboxyl terminus of GRK2, a known G betagamma sequestrant. Furthermore, in cells overexpressing both GRK2 and G beta1 gamma2, activation of lysophosphatidic acid receptors leads to the rapid and transient formation of a GRK/G betagamma complex. That G betagamma specificity exists at the level of the GPCR and the GRK is indicated by the observation that a GRK2/G betagamma complex is formed after agonist occupancy of the lysophosphatidic acid and beta-adrenergic but not thrombin receptors. In contrast to GRK2, GRK3 forms a G betagamma complex after stimulation of all three GPCRs. This G betagamma binding specificity of the GRKs is also reflected at the level of the purified proteins. Thus the GRK2 carboxyl terminus binds G beta1 and G beta2 but not G beta3, while the GRK3 fusion protein binds all three G beta isoforms. This study provides a direct demonstration of a role for G betagamma in mediating the agonist-stimulated translocation of GRK2 and GRK3 in an intact cellular system and demonstrates isoform specificity in the interaction of these components.

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.

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.

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.

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.

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.