Stabilization of the human beta2-adrenergic receptor TM4-TM3-TM5 helix interface by mutagenesis of Glu122(3.41), a critical residue in GPCR structure

G protein-coupled receptor (GPCR) instability represents one of the most profound obstacles in the structural study of GPCRs that bind diffusible ligands. The introduction of targeted mutations at nonconserved residues that lie proximal to helix interfaces has the potential to enhance the fold stability of the receptor helix bundle while maintaining wild-type receptor function. To test this hypothesis, we studied the effect of amino acid substitutions at Glu122^3.41 in the well-studied β₂-adrenergic receptor (β₂AR), which was predicted from sequence conservation to lie at a position equivalent to a tryptophan residue in rhodopsin at the 3,4,5 helix interface among transmembrane (TM) domains 3, 4, and 5. Replacement of Glu122^3.41 with bulky hydrophobic residues, such as tryptophan, tyrosine, and phenylalanine, increases the yield of functionally folded β₂AR by as much as 5-fold. Receptor stability in detergent solution was studied by isothermal denaturation, and it was found that the E122W and E122Y mutations enhanced the β₂AR thermal half-life by 9.3- and 6.7-fold, respectively, at 37 °C. The β₁AR was also stabilized by the introduction of tryptophan at Glu147^3.41, and the effect on protein behavior was similar to the rescue of the unstable wild-type receptor by the antagonist propranolol. Molecular modeling of the E122W and E122Y mutants revealed that the tryptophan ring edge and tyrosine hydroxyl are positioned proximal to the helical break in TM5 introduced by the conserved Pro211^5.50 and may stabilize the helix by interacting favorably with the unpaired carbonyl oxygen of Val206^5.45. Conformational flexibility of TM5 is likely to be a general property of class A GPCRs; therefore, engineering of the TM4-TM3-TM5 interface at the 3.41 position may provide a general strategy for the stabilization of other receptors.     

Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis

The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.     

Molecular dynamics simulations reveal insights into key structural elements of adenosine receptors

The crystal structure of the human A(2A) adenosine receptor, a member of the G protein-coupled receptor (GPCR) family, is used as a starting point for the structural characterization of the conformational equilibrium around the inactive conformation of the human A(2) (A(2A) and A(2B)) adenosine receptors (ARs). A homology model of the closely related A(2B)AR is reported, and the two receptors were simulated in their apo form through all-atom molecular dynamics (MD) simulations. Different conditions were additionally explored in the A(2A)AR, including the protonation state of crucial histidines or the presence of the cocrystallized ligand. Our simulations reveal the role of several conserved residues in the ARs in the conformational equilibrium of the receptors. The "ionic lock" absent in the crystal structure of the inactive A(2A)AR is rapidly formed in the two simulated receptors, and a complex network of interacting residues is presented that further stabilizes this structural element. Notably, the observed rotameric transition of Trp6.48 ("toggle switch"), which is thought to initiate the activation process in GPCRs, is accompanied by a concerted rotation of the conserved residue of the A(2)ARs, His6.52. This new conformation is further stabilized in the two receptors under study by a novel interaction network involving residues in transmembrane (TM) helices TM5 (Asn5.42) and TM3 (Gln3.37), which resemble the conformational changes recently observed in the agonist-bound structure of β-adrenoreceptors. Finally, the interaction between Glu1.39 and His7.43, a pair of conserved residues in the family of ARs, is found to be weaker than previously thought, and the role of this interaction in the structure and dynamics of the receptor is thoroughly examined. All these findings suggest that, despite the commonalities with other GPCRs, the conformational equilibrium of ARs is also modulated by specific residues of the family.     

Relevance of rhodopsin studies for GPCR activation

Rhodopsin, the dim-light photoreceptor present in the rod cells of the retina, is both a retinal-binding protein and a G protein-coupled receptor (GPCR). Due to this conjunction, it benefits from an arsenal of spectroscopy techniques that can be used for its characterization, while being a model system for the important family of Class A (also referred to as “rhodopsin-like”) GPCRs. For instance, rhodopsin has been a crucial player in the field of GPCR structural biology. Until 2007, it was the only GPCR for which a high-resolution crystal structure was available, so all structure–activity analyses on GPCRs, from structure-based drug discovery to studies of structural changes upon activation, were based on rhodopsin. At present, about a third of currently available GPCR structures are still from rhodopsin. In this review, I show some examples of how these structures can still be used to gain insight into general aspects of GPCR activation. First, the analysis of the third intracellular loop in rhodopsin structures allows us to gain an understanding of the structural and dynamic properties of this region, which is absent (due to protein engineering or poor electron density) in most of the currently available GPCR structures. Second, a detailed analysis of the structure of the transmembrane domains in inactive, intermediate and active rhodopsin structures allows us to detect early conformational changes in the process of ligand-induced GPCR activation. Finally, the analysis of a conserved ligand-activated transmission switch in the transmembrane bundle of GPCRs in the context of the rhodopsin activation cycle, allows us to suggest that the structures of many of the currently available agonist-bound GPCRs may correspond to intermediate active states. While the focus in GPCR structural biology is inevitably moving away from rhodopsin, in other aspects rhodopsin is still at the forefront. For instance, the first studies of the structural basis of disease mutants in GPCRs, or the most detailed analysis of cellular GPCR signal transduction networks using a systems biology approach, have been carried out in rhodopsin. Finally, due again to its unique properties among GPCRs, rhodopsin will likely play an important role in the application of X-ray free electron laser crystallography to time-resolved structural biology in membrane proteins. Rhodopsin, thus, still remains relevant as a model system to study the molecular mechanisms of GPCR activation. This article is part of a Special Issue entitled: Retinal Proteins—You can teach an old dog new tricks.     

Putative Active States of a Prototypic G-Protein-Coupled Receptor from Biased Molecular Dynamics

A major current focus of structural work on G-protein-coupled receptors (GPCRs) pertains to the investigation of their active states. However, for virtually all GPCRs, active agonist-bound intermediate states have been difficult to characterize experimentally owing to their higher conformational flexibility, and thus intrinsic instability, as compared to inactive inverse agonist-bound states. In this work, we explored possible activation pathways of the prototypic GPCR bovine rhodopsin by means of biased molecular dynamics simulations. Specifically, we used an explicit atomistic representation of the receptor and its environment, and sampled the conformational transition from the crystal structure of a photoactivated deprotonated state of rhodopsin to the low pH crystal structure of opsin in the presence of 11-trans-retinal, using adiabatic biased molecular dynamics simulations. We then reconstructed the system free-energy landscape along the predetermined transition trajectories using a path collective variable approach based on metadynamics. Our results suggest that the two experimental endpoints of rhodopsin/opsin are connected by at least two different pathways, and that the conformational transition is populated by at least four metastable states of the receptor, characterized by a different amplitude of the outward movement of transmembrane helix 6.     

Structural basis of activation of bitter taste receptor T2R1 and comparison with Class A G-protein-coupled receptors (GPCRs)

The human bitter taste receptors (T2Rs) are non-Class A members of the G-protein-coupled receptor (GPCR) superfamily, with very limited structural information. Amino acid sequence analysis reveals that most of the important motifs present in the transmembrane helices (TM1-TM7) of the well studied Class A GPCRs are absent in T2Rs, raising fundamental questions regarding the mechanisms of activation and how T2Rs recognize bitter ligands with diverse chemical structures. In this study, the bitter receptor T2R1 was used to systematically investigate the role of 15 transmembrane amino acids in T2Rs, including 13 highly conserved residues, by amino acid replacements guided by molecular modeling. Functional analysis of the mutants by calcium imaging analysis revealed that replacement of Asn-66(2.65) and the highly conserved Asn-24(1.50) resulted in greater than 90% loss of agonist-induced signaling. Our results show that Asn-24(1.50) plays a crucial role in receptor activation by mediating an hydrogen bond network connecting TM1-TM2-TM7, whereas Asn-66(2.65) is essential for binding to the agonist dextromethorphan. The interhelical hydrogen bond between Asn-24(1.50) and Arg-55(2.54) restrains T2R receptor activity because loss of this bond in I27A and R55A mutants results in hyperactive receptor. The conserved amino acids Leu-197(5.50), Ser-200(5.53), and Leu-201(5.54) form a putative LXXSL motif which performs predominantly a structural role by stabilizing the helical conformation of TM5 at the cytoplasmic end. This study provides for the first time mechanistic insights into the roles of the conserved transmembrane residues in T2Rs and allows comparison of the activation mechanisms of T2Rs with the Class A GPCRs.     

Modeling GPCR active state conformations: the β(2)-adrenergic receptor

The recent publication of several G protein-coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β(2)-adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site-directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site-directed mutagenesis information. Virtual screening experiments show that the models are selective for β-adrenergic agonists over other GPCR ligands, for (R)- over (S)-β-hydroxy agonists and for β(2)-selective agonists over β(1)-selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C-terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C-terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR-G protein model has been used to explain site-directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β(2)-adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs.     

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.     

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism

Background

G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states.

Results

Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 ų reduction in volume.

Conclusions

In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0567-3) contains supplementary material, which is available to authorized users.     

Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors

Activation of G-protein-coupled receptors (GPCRs) is initiated by conformational changes in the transmembrane (TM) helices and the intra- and extracellular loops induced by ligand binding. Understanding the conformational changes in GPCRs leading to activation is imperative in deciphering the role of these receptors in the pathology of diseases. Since the crystal structures of activated GPCRs are not yet available, computational methods and biophysical techniques have been used to predict the structures of GPCR active states. We have recently applied the computational method LITiCon to understand the ligand-induced conformational changes in β₂-adrenergic receptor by ligands of varied efficacies. Here we report a study of the conformational changes associated with the activation of bovine rhodopsin for which the crystal structure of the inactive state is known. Starting from the inactive (dark) state, we have predicted the TM conformational changes that are induced by the isomerization of 11-cis retinal to all-trans retinal leading to the fully activated state, metarhodopsin II. The predicted active state of rhodopsin satisfies all of the 30 known experimental distance constraints. The predicted model also correlates well with the experimentally observed conformational switches in rhodopsin and other class A GPCRs, namely, the breaking of the ionic lock between R135^3.50 at the intracellular end of TM3 (part of the DRY motif) and E247^6.30 on TM6, and the rotamer toggle switch on W265^6.48 on TM6. We observe that the toggling of the W265^6.48 rotamer modulates the bend angle of TM6 around the conserved proline. The rotamer toggling is facilitated by the formation of a water wire connecting S298^7.45, W265^6.48 and H211^5.46. As a result, the intracellular ends of TMs 5 and 6 move outward from the protein core, causing large conformational changes at the cytoplasmic interface. The predicted outward movements of TM5 and TM6 are in agreement with the recently published crystal structure of opsin, which is proposed to be close to the active-state structure. In the predicted active state, several residues in the intracellular loops, such as R69, V139^3.54, T229, Q237, Q239, S240, T243 and V250^6.33, become more water exposed compared to the inactive state. These residues may be involved in mediating the conformational signal from the receptor to the G protein. From mutagenesis studies, some of these residues, such as V139^3.54, T229 and V250^6.33, are already implicated in G-protein activation. The predicted active state also leads to the formation of new stabilizing interhelical hydrogen-bond contacts, such as those between W265^6.48 and H211^5.46 and E122^3.37 and C167^4.56. These hydrogen-bond contacts serve as potential conformational switches offering new opportunities for future experimental investigations. The calculated retinal binding energy surface shows that binding of an agonist makes the receptor dynamic and flexible and accessible to many conformations, while binding of an inverse agonist traps the receptor in the inactive state and makes the other conformations inaccessible.     

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.     

Emerging structural biology of lipid G protein‐coupled receptors

The first crystal structure of a G protein‐coupled receptor (GPCR) was that of the bovine rhodopsin, solved in 2000, and is a light receptor within retina rode cells that enables vision by transducing a conformational signal from the light‐induced isomerization of retinal covalently bound to the receptor. More than 7 years after this initial discovery and following more than 20 years of technological developments in GPCR expression, stabilization, and crystallography, the high‐resolution structure of the adrenaline binding β₂‐adrenergic receptor, a ligand diffusible receptor, was discovered. Since then, high‐resolution structures of more than 53 unique GPCRs have been determined leading to a significant improvement in our understanding of the basic mechanisms of ligand‐binding and ligand‐mediated receptor activation that revolutionized the field of structural molecular pharmacology of GPCRs. Recently, several structures of eight unique lipid‐binding receptors, one of the most difficult GPCR families to study, have been reported. This review presents the outstanding structural and pharmacological features that have emerged from these new lipid receptor structures. The impact of these findings goes beyond mechanistic insights, providing evidence of the fundamental role of GPCRs in the physiological integration of the lipid signaling system, and highlighting the importance of sustained research into the structural biology of GPCRs for the development of new therapeutics targeting lipid receptors.     

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.     

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.     

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment

Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.     

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.     

Active state-like conformational elements in the beta2-AR and a photoactivated intermediate of rhodopsin identified by dynamic properties of GPCRs

G-Protein-coupled receptors (GPCRs) adopt various functionally relevant conformational states in cell signaling processes. Recently determined crystal structures of rhodopsin and the beta 2-adrenergic receptor (beta 2-AR) offer insight into previously uncharacterized active conformations, but the molecular states of these GPCRs are likely to contain both inactive and active-like conformational elements. We have identified conformational rearrangements in the dynamics of the TM7-HX8 segment that relate to the properties of the conserved NPxxY(x)5,6F motif and show that they can be used to identify active state-like conformational elements in the corresponding regions of the new structures of rhodopsin and the beta 2-AR.     

Agonist-bound structures of G protein-coupled receptors

G protein-coupled receptors (GPCRs) play a major role in intercellular communication by binding small diffusible ligands (agonists) at the extracellular surface. Agonist-binding induces a conformational change in the receptor, which results in the binding and activation of heterotrimeric G proteins within the cell. Ten agonist-bound structures of non-rhodopsin GPCRs published last year defined for the first time the molecular details of receptor activated states and how inverse agonists, partial agonists and full agonists bind to produce different effects on the receptor. In addition, the structure of the β(2)-adrenoceptor coupled to a heterotrimeric G protein showed how the opening of a cleft in the cytoplasmic face of the receptor as a consequence of agonist binding results in G protein coupling and activation of the G protein.     

Protein‐based virtual screening of chemical databases. II. Are homology models of g‐protein coupled receptors suitable targets?

The aim of the current study is to investigate whether homology models of G-Protein-Coupled Receptors (GPCRs) that are based on bovine rhodopsin are reliable enough to be used for virtual screening of chemical databases. Starting from the recently described 2.8 A-resolution X-ray structure of bovine rhodopsin, homology models of an "antagonist-bound" form of three human GPCRs (dopamine D3 receptor, muscarinic M1 receptor, vasopressin V1a receptor) were constructed. The homology models were used to screen three-dimensional databases using three different docking programs (Dock, FlexX, Gold) in combination with seven scoring functions (ChemScore, Dock, FlexX, Fresno, Gold, Pmf, Score). Rhodopsin-based homology models turned out to be suitable, indeed, for virtual screening since known antagonists seeded in the test databases could be distinguished from randomly chosen molecules. However, such models are not accurate enough for retrieving known agonists. To generate receptor models better suited for agonist screening, we developed a new knowledge- and pharmacophore-based modeling procedure that might partly simulate the conformational changes occurring in the active site during receptor activation. Receptor coordinates generated by this new procedure are now suitable for agonist screening. We thus propose two alternative strategies for the virtual screening of GPCR ligands, relying on a different set of receptor coordinates (antagonist-bound and agonist-bound states).