Comparative Sequence and Structural Analyses of G-Protein-Coupled Receptor Crystal Structures and Implications for Molecular Models

Background

Up until recently the only available experimental (high resolution) structure of a G-protein-coupled receptor (GPCR) was that of bovine rhodopsin. In the past few years the determination of GPCR structures has accelerated with three new receptors, as well as squid rhodopsin, being successfully crystallized. All share a common molecular architecture of seven transmembrane helices and can therefore serve as templates for building molecular models of homologous GPCRs. However, despite the common general architecture of these structures key differences do exist between them. The choice of which experimental GPCR structure(s) to use for building a comparative model of a particular GPCR is unclear and without detailed structural and sequence analyses, could be arbitrary. The aim of this study is therefore to perform a systematic and detailed analysis of sequence-structure relationships of known GPCR structures.

Methodology

We analyzed in detail conserved and unique sequence motifs and structural features in experimentally-determined GPCR structures. Deeper insight into specific and important structural features of GPCRs as well as valuable information for template selection has been gained. Using key features a workflow has been formulated for identifying the most appropriate template(s) for building homology models of GPCRs of unknown structure. This workflow was applied to a set of 14 human family A GPCRs suggesting for each the most appropriate template(s) for building a comparative molecular model.

Conclusions

The available crystal structures represent only a subset of all possible structural variation in family A GPCRs. Some GPCRs have structural features that are distributed over different crystal structures or which are not present in the templates suggesting that homology models should be built using multiple templates. This study provides a systematic analysis of GPCR crystal structures and a consistent method for identifying suitable templates for GPCR homology modelling that will help to produce more reliable three-dimensional models.     

Structure-based identification of binding sites,native ligands and potential inhibitors for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are the largest family of cell-surface receptors involved in signal transmission. Drugs associated with GPCRs represent more than one fourth of the 100 top-selling drugs and are the targets of more than half of the current therapeutic agents on the market. Our methodology based on the internal coordinate mechanics (ICM) program can accurately identify the ligand-binding pocket in the currently available crystal structures of seven transmembrane (7TM) proteins [bacteriorhodopsin (BR) and bovine rhodopsin (bRho)]. The binding geometry of the ligand can be accurately predicted by ICM flexible docking with and without the loop regions, a useful finding for GPCR docking because the transmembrane regions are easier to model. We also demonstrate that the native ligand can be identified by flexible docking and scoring in 1.5% and 0.2% (for bRho and BR, respectively) of the best scoring compounds from two different types of compound database. The same procedure can be applied to the database of available chemicals to identify specific GPCR binders. Finally, we demonstrate that even if the sidechain positions in the bRho binding pocket are entirely wrong, their correct conformation can be fully restored with high accuracy (0.28 A) through the ICM global optimization with and without the ligand present. These binding site adjustments are critical for flexible docking of new ligands to known structures or for docking to GPCR homology models. The ICM docking method has the potential to be used to "de-orphanize" orphan GPCRs (oGPCRs) and to identify antagonists-agonists for GPCRs if an accurate model (experimentally and computationally validated) of the structure has been constructed or when future crystal structures are determined.     

3D Modeling for TM Receptors: Algorithms and Validations

Summary This paper outlines the basic strategy to build 3D models of transmembrane G-protein coupled receptors (GPCRs) starting from their amino acid sequences in a ‘block-by-block’ manner: (i) locate possible TM helical fragments in the GPCR sequence; (ii) build 3D structures for these helices; (iii) arrange isolated helices across the membrane; (iv) calculate all pairwise helix-helix interactions; (v) assemble helical bundle(s); (vi) restore interhelical loops and N- and C-termini; and (vii) refine the entire 3D structure(s). Computer algorithms and preliminary results for most of the steps are discussed.     

Towards Improved Quality of GPCR Models by Usage of Multiple Templates and Profile-Profile Comparison

G-protein coupled receptors (GPCRs) are targets of nearly one third of the drugs at the current pharmaceutical market. Despite their importance in many cellular processes the crystal structures are available for less than 20 unique GPCRs of the Rhodopsin-like class. Fortunately, even though involved in different signaling cascades, this large group of membrane proteins has preserved a uniform structure comprising seven transmembrane helices that allows quite reliable comparative modeling. Nevertheless, low sequence similarity between the GPCR family members is still a serious obstacle not only in template selection but also in providing theoretical models of acceptable quality. An additional level of difficulty is the prediction of kinks and bulges in transmembrane helices. Usage of multiple templates and generation of alignments based on sequence profiles may increase the rate of success in difficult cases of comparative modeling in which the sequence similarity between GPCRs is exceptionally low. Here, we present GPCRM, a novel method for fast and accurate generation of GPCR models using averaging of multiple template structures and profile-profile comparison. In particular, GPCRM is the first GPCR structure predictor incorporating two distinct loop modeling techniques: Modeller and Rosetta together with the filtering of models based on the Z-coordinate. We tested our approach on all unique GPCR structures determined to date and report its performance in comparison with other computational methods targeting the Rhodopsin-like class. We also provide a database of precomputed GPCR models of the human receptors from that class.

Availability

GPCRM server and database: http://gpcrm.biomodellab.eu     

Revisiting the homology modeling of G-protein coupled receptors: β1-adrenoceptor as an example

G-protein coupled receptors (GPCRs) are recognized to constitute the largest family of membrane proteins. Due to the disproportion in the quantity of crystal structures and their amino acid sequences, homology modeling contributes a reasonable and feasible approach to GPCR theoretical coordinates. With the brand new crystal structures resolved recently, herein we deliberated how to designate them as templates to carry out homology modeling in four aspects: (1) various sequence alignment methods; (2) protein weight matrix; (3) different sets of multiple templates; (4) active and inactive state of templates. The accuracy of models was evaluated by comparing the similarity of stereo conformation and molecular docking results between models and the experimental structure of Meleagris gallopavo β(1)-adrenergic receptor (Mg_Adrb1) that we desired to develop as an example. Our results proposed that: (1) Cobalt and MAFFT, two algorithms of sequence alignment, were suitable for single- and multiple-template modeling, respectively; (2) Blosum30 is applicable to align sequences in the case of low sequence identity; (3) multiple-template modeling is not always better than single-template one; (4) the state of template is an influential factor in simulating the GPCR structures as well.     

Critical analysis of the successes and failures of homology models of G protein‐coupled receptors

We present a critical assessment of the performance of our homology model refinement method for G protein‐coupled receptors (GPCRs), called LITICon that led to top ranking structures in a recent structure prediction assessment GPCRDOCK2010. GPCRs form the largest class of drug targets for which only a few crystal structures are currently available. Therefore, accurate homology models are essential for drug design in these receptors. We submitted five models each for human chemokine CXCR4 (bound to small molecule IT1t and peptide CVX15) and dopamine D3DR (bound to small molecule eticlopride) before the crystal structures were published. Our models in both CXCR4/IT1t and D3/eticlopride assessments were ranked first and second, respectively, by ligand RMSD to the crystal structures. For both receptors, we developed two types of protein models: homology models based on known GPCR crystal structures, and ab initio models based on the prediction method MembStruk. The homology‐based models compared better to the crystal structures than the ab initio models. However, a robust refinement procedure for obtaining high accuracy structures is needed. We demonstrate that optimization of the helical tilt, rotation, and translation is vital for GPCR homology model refinement. As a proof of concept, our in‐house refinement program LITiCon captured the distinct orientation of TM2 in CXCR4, which differs from that of adrenoreceptors. These findings would be critical for refining GPCR homology models in future. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Structural biology of G protein‐coupled receptor signaling complexes

G protein‐coupled receptors (GPCRs) constitute the largest family of cell surface receptors that mediate numerous cell signaling pathways, and are targets of more than one‐third of clinical drugs. Thanks to the advancement of novel structural biology technologies, high‐resolution structures of GPCRs in complex with their signaling transducers, including G‐protein and arrestin, have been determined. These 3D complex structures have significantly improved our understanding of the molecular mechanism of GPCR signaling and provided a structural basis for signaling‐biased drug discovery targeting GPCRs. Here we summarize structural studies of GPCR signaling complexes with G protein and arrestin using rhodopsin as a model system, and highlight the key features of GPCR conformational states in biased signaling including the sequence motifs of receptor TM6 that determine selective coupling of G proteins, and the phosphorylation codes of GPCRs for arrestin recruitment. We envision the future of GPCR structural biology not only to solve more high‐resolution complex structures but also to show stepwise GPCR signaling complex assembly and disassembly and dynamic process of GPCR signal transduction.     

Signal transmission via G protein-coupled receptors in the light of rhodopsin structure determination

G protein-coupled receptors (GPCRs) transducing diverse external signals to cells via activation of heterotrimeric GTP-binding (G) proteins, estimated to mediate actions of 60% of drugs, had been resistant to structure determination until summer 2000. The first atomic-resolution experimental structure of a GPCR, that of dark (inactive) rhodopsin, thus provides a trustworthy 3D prototype for antagonist-bound forms of this huge family of proteins. In this work, our former theoretical GPCR models are evaluated against the new experimental template. Subsequently, a working hypothesis regarding the signal transduction mechanism by GPCRs is presented.     

Conserved activation pathways in G-protein-coupled receptors

GPCRs (G-protein-coupled receptors) are seven-transmembrane helix proteins that transduce exogenous and endogenous signals to modulate the activity of downstream effectors inside the cell. Despite the relevance of these proteins in human physiology and pharmaceutical research, we only recently started to understand the structural basis of their activation mechanism. In the period 2008-2011, nine active-like structures of GPCRs were solved. Among them, we have determined the structure of light-activated rhodopsin with all the features of the active metarhodopsin-II, which represents so far the most native-like model of an active GPCR. This structure, together with the structures of other inactive, intermediate and active states of rhodopsin constitutes a unique structural framework on which to understand the conserved aspects of the activation mechanism of GPCRs. This mechanism can be summarized as follows: retinal isomerization triggers a series of local structural changes in the binding site that are amplified into three intramolecular activation pathways through TM (transmembrane helix) 5/TM3, TM6 and TM7/TM2. Sequence analysis strongly suggests that these pathways are conserved in other GPCRs. Differential activation of these pathways by ligands could be translated into the stabilization of different active states of the receptor with specific signalling properties.     

Fractal dimension as a measure of surface roughness of G protein-coupled receptors: implications for structure and function

Protein surface roughness is a structural property associated with ligand-protein and protein-protein binding interfaces. In this work we apply for the first time the concept of surface roughness, expressed as the fractal dimension, to address structure and function of G protein-coupled receptors (GPCRs) which are an important group of drug targets. We calculate the exposure ratio and the fractal dimension for helix-forming residues of the β(2) adrenergic receptor (β(2)AR), a model system in GPCR studies, in different conformational states: in complex with agonist, antagonist and partial inverse agonists. We show that both exposure ratio and roughness exhibit periodicity which results from the helical structure of GPCRs. The pattern of roughness and exposure ratio of a protein patch depends on its environment: the residues most exposed to membrane are in general most rough whereas parts of receptors mediating interhelical contacts in a monomer or protein complex are much smoother. We also find that intracellular ends (TM3, TM5, TM6 and TM7) which are relevant for G protein binding and thus receptor signaling, are exposed but smooth. Mapping the values of residual fractal dimension onto receptor 3D structures makes it possible to conclude that the binding sites of orthosteric ligands as well as of cholesterol are characterized with significantly higher roughness than the average for the whole protein. In summary, our study suggests that identification of specific patterns of roughness could be a novel approach to spot possible binding sites which could serve as original drug targets for GPCRs modulation.     

Structure modeling of all identified G protein-coupled receptors in the human genome

G protein-coupled receptors (GPCRs), encoded by about 5% of human genes, comprise the largest family of integral membrane proteins and act as cell surface receptors responsible for the transduction of endogenous signal into a cellular response. Although tertiary structural information is crucial for function annotation and drug design, there are few experimentally determined GPCR structures. To address this issue, we employ the recently developed threading assembly refinement (TASSER) method to generate structure predictions for all 907 putative GPCRs in the human genome. Unlike traditional homology modeling approaches, TASSER modeling does not require solved homologous template structures; moreover, it often refines the structures closer to native. These features are essential for the comprehensive modeling of all human GPCRs when close homologous templates are absent. Based on a benchmarked confidence score, approximately 820 predicted models should have the correct folds. The majority of GPCR models share the characteristic seven-transmembrane helix topology, but 45 ORFs are predicted to have different structures. This is due to GPCR fragments that are predominantly from extracellular or intracellular domains as well as database annotation errors. Our preliminary validation includes the automated modeling of bovine rhodopsin, the only solved GPCR in the Protein Data Bank. With homologous templates excluded, the final model built by TASSER has a global C(alpha) root-mean-squared deviation from native of 4.6 angstroms, with a root-mean-squared deviation in the transmembrane helix region of 2.1 angstroms. Models of several representative GPCRs are compared with mutagenesis and affinity labeling data, and consistent agreement is demonstrated. Structure clustering of the predicted models shows that GPCRs with similar structures tend to belong to a similar functional class even when their sequences are diverse. These results demonstrate the usefulness and robustness of the in silico models for GPCR functional analysis. All predicted GPCR models are freely available for noncommercial users on our Web site (http://www.bioinformatics.buffalo.edu/GPCR).     

3D Modeling for TM Receptors: Algorithms and Validations

This paper outlines the basic strategy to build 3D models of transmembrane G-protein coupled receptors (GPCRs) starting from their amino acid sequences in a block-by-block manner: (i) locate possible TM helical fragments in the GPCR sequence; (ii) build 3D structures for these helices; (iii) arrange isolated helices across the membrane; (iv) calculate all pairwise helix-helix interactions; (v) assemble helical bundle(s); (vi) restore interhelical loops and N- and C-termini; and (vii) refine the entire 3D structure(s). Computer algorithms and preliminary results for most of the steps are discussed.     

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.     

3D modeling of the activated states of constitutively active mutants of rhodopsin

The activated (R*) states in constitutively active mutants (CAMs) of G-protein-coupled receptors (GPCRs) are presumably characterized by lower energies than the resting (R) states. If specific configurations of TM helices differing by rotations along the long transmembrane axes possess energies lower than that in the R state for pronounced CAMs, but not for non-CAMs, these particular configurations of TM helices are candidate 3D models for the R* state. The hypothesis was studied in the case of rhodopsin, the only GPCR for which experimentally determined 3D models of the R and R* states are currently available. Indeed, relative energies of the R* state were significantly lower than that of the R state for the rhodopsin mutants G90D/M257Y and E113Q/M257Y (strong CAMs), but not for G90D, E113Q, and M257Y (not CAMs). Next, the developed build-up procedure successfully identified few similar configurations of the TM helical bundle of G90D/M257Y and E113Q/M257Y as possible candidates for the 3D model of the R* state of rhodopsin, all of them being in good agreement with the model suggested by experiment. Since constitutively active mutants are known for many of GPCRs belonging to the large rhodopsin-like family, this approach provides a way for predicting possible 3D structures corresponding to the activated states of the TM regions of many GPCRs for which CAMs have been identified.     

Constitutive activation of angiotensin II type 1 receptor alters the orientation of transmembrane Helix-2

A key step in transmembrane (TM) signal transduction by G-protein-coupled receptors (GPCRs) is the ligand-induced conformational change of the receptor, which triggers the activation of a guanine nucleotide-binding protein. GPCRs contain a seven-TM helical structure essential for signal transduction in response to a large variety of sensory and hormonal signals. Primary structure comparison of GPCRs has shown that the second TM helix contains a highly conserved Asp residue, which is critical for agonist activation in these receptors. How conformational changes in TM2 relate to signal transduction by a GPCR is not known, because activation-induced conformational changes in TM2 helix have not been measured. Here we use modification of reporter cysteines to measure water accessibility at specific residues in TM2 of the type 1 receptor for the octapeptide hormone angiotensin II. Activation-dependent changes in the accessibility of Cys76 on TM2 were measured in constitutively activated mutants. These changes were directly correlated with measurement of function, establishing the link between physical changes in TM2 and function. Accessibility changes were measured at several consecutive residues on TM2, which suggest that TM2 undergoes a transmembrane movement in response to activation. This is the first report of in situ measurement of TM2 movement in a GPCR.     

NMR structure of the second intracellular loop of the alpha 2A adrenergic receptor: evidence for a novel cytoplasmic helix

A major, unresolved question in signal transduction by G protein coupled receptors (GPCRs) is to understand how, at atomic resolution, a GPCR activates a G protein. A step toward answering this question was made with the determination of the high-resolution structure of rhodopsin; we now know the intramolecular interactions that characterize the resting conformation of a GPCR. To what degree does this structure represent a structural paradigm for other GPCRs, especially at the cytoplasmic surface where GPCR-G protein interaction occurs and where the sequence homology is low among GPCRs? To address this question, we performed NMR studies on approximately 35-residue-long peptides including the critical second intracellular loop (i2) of the alpha 2A adrenergic receptor (AR) and of rhodopsin. To stabilize the secondary structure of the peptide termini, 4-12 residues from the adjacent transmembrane helices were included and structures determined in dodecylphosphocholine micelles. We also characterized the effects on an alpha 2A AR peptide of a D130I mutation in the conserved DRY motif. Our results show that in contrast to the L-shaped loop in the i2 of rhodopsin, the i2 of the alpha 2A AR is predominantly helical, supporting the hypothesis that there is structural diversity within GPCR intracellular loops. The D130I mutation subtly modulates the helical structure. The spacing of nonpolar residues in i2 with helical periodicity is a predictor of helical versus loop structure. These data should lead to more accurate models of the intracellular surface of GPCRs and of receptor-mediated G protein activation.     

A transmembrane helix-bundle from G-protein coupled receptor CB2: biosynthesis, purification, and NMR characterization

The cannabinoid receptor subtype 2 (CB2) is a member of the G-protein coupled receptor (GPCR) superfamily. As the relationship between structure and function for this receptor remains poorly understood, the present study was undertaken to characterize the structure of a segment including the first and second transmembrane helix (TM1 and TM2) domains of CB2. To accomplish this, a transmembrane double-helix bundle from this region was expressed, purified, and characterized by NMR. Milligrams of this hydrophobic fragment of the receptor were biosynthesized using a fusion protein overexpression strategy and purified by affinity chromatography combined with reverse phase HPLC. Chemical and enzymatic cleavage methods were implemented to remove the fusion tag. The resultant recombinant protein samples were analyzed and confirmed by HPLC, mass spectrometry, and circular dichroism (CD). The CD analyses of HPLC-purified protein in solution and in DPC micelle preparations suggested predominant alpha-helical structures under both conditions. The 13C/15N double-labeled protein CB2(27-101) was further verified and analyzed by NMR spectroscopy. Sequential assignment was accomplished for more than 80% of residues. The 15N HSQC NMR results show a clear chemical shift dispersion of the amide nitrogen-proton correlation indicative of a pure double-labeled polypeptide molecule. The results suggest that this method is capable of generating transmembrane helical bundles from GPCRs in quantity and purity sufficient for NMR and other biophysical studies. Therefore, the biosynthesis of GPCR transmembrane helix bundles represents a satisfactory alternative strategy to obtain and assemble NMR structures from recombinant "building blocks."     

Structure and pH sensitivity of the transmembrane segment 3 of rhodopsin

Activation of G protein-coupled receptors (GPCRs) originates in ligand-induced protein conformational changes that are transmitted to the cytosolic receptor surface. In the photoreceptor rhodopsin, and possibly other rhodopsin-like GPCRs, protonation of a carboxylic acid in the conserved E(D)RY motif at the cytosolic end of transmembrane helix 3 (TM3) is coupled to receptor activation. Here, we have investigated the structure of synthetic peptides derived from rhodopsin TM3. Polarized FTIR spectroscopy reveals a helical structure of a 31-mer TM3 peptide reconstituted into PC vesicles with a large tilt of 40-50 degrees of the helical axis relative to the membrane normal. Helical structure is also observed for the TM3 peptide in detergent micelles and depends on pH, especially in the C-terminal sequence. In addition, the fluorescence emission of the single tyrosine of the D(E)RY motif in the TM3 peptide exhibits a pronounced pH sensitivity that is abolished when Glu is replaced by Gln, demonstrating that protonation of the conserved Glu side chain affects the structure in the environment of the D(E)RY motif of TM3. The pH regulation of the C-terminal TM3 structure may be an intrinsic feature of the E(D)RY motif in other class I receptors, allowing the coupling of protonation and conformation of membrane-exposed residues in full-length GPCRs.     

Selective reconstitution of human D4 dopamine receptor variants with Gi alpha subtypes

G protein-coupled receptors (GPCRs) are seven-transmembrane (TM) helical proteins that bind extracellular molecules and transduce signals by coupling to heterotrimeric G proteins in the cytoplasm. The human D4 dopamine receptor is a particularly interesting GPCR because the polypeptide loop linking TM helices 5 and 6 (loop i3) may contain from 2 to 10 similar direct hexadecapeptide repeats. The precise role of loop i3 in D4 receptor function is not known. To clarify the role of loop i3 in G protein coupling, we constructed synthetic genes for the three main D4 receptor variants. D4-2, D4-4, and D4-7 receptors contain 2, 4, and 7 imperfect hexadecapeptide repeats in loop i3, respectively. We expressed and characterized the synthetic genes and found no significant effect of the D4 receptor polymorphisms on antagonist or agonist binding. We developed a cell-based assay where activated D4 receptors coupled to a Pertussis toxin-sensitive pathway to increase intracellular calcium concentration. Studies using receptor mutants showed that the regions of loop i3 near TM helices 5 and 6 were required for G protein coupling. The hexadecapeptide repeats were not required for G protein-mediated calcium flux. Cell membranes containing expressed D4 receptors and receptor mutants were reconstituted with purified recombinant G protein alpha subunits. The results show that each D4 receptor variant is capable of coupling to several G(i)alpha subtypes. Furthermore, there is no evidence of any quantitative difference in G protein coupling related to the number of hexadecapeptide repeats in loop i3. Thus, loop i3 is required for D4 receptors to activate G proteins. However, the polymorphic region of the loop does not appear to affect the specificity or efficiency of G(i)alpha coupling.     

TESS: a geometric hashing algorithm for deriving 3D coordinate templates for searching structural databases. Application to enzyme active sites.

It is well established that sequence templates such as those in the PROSITE and PRINTS databases are powerful tools for predicting the biological function and tertiary structure for newly derived protein sequences. The number of X-ray and NMR protein structures is increasing rapidly and it is apparent that a 3D equivalent of the sequence templates is needed. Here, we describe an algorithm called TESS that automatically derives 3D templates from structures deposited in the Brookhaven Protein Data Bank. While a new sequence can be searched for sequence patterns, a new structure can be scanned against these 3D templates to identify functional sites. As examples, 3D templates are derived for enzymes with an O-His-O "catalytic triad" and for the ribonucleases and lysozymes. When these 3D templates are applied to a large data set of nonidentical proteins, several interesting hits are located. This suggests that the development of a 3D template database may help to identify the function of new protein structures, if unknown, as well as to design proteins with specific functions.