Three-dimensional model for meta-II rhodopsin, an activated G-protein-coupled receptor

Conserved residues in the extracellular loops of short-wavelength cone visual pigments

The role of the extracellular loop region of a short-wavelength sensitive pigment, Xenopus violet cone opsin, is investigated via computational modeling, mutagenesis, and spectroscopy. The computational models predict a complex H-bonding network that stabilizes and connects the EC2-EC3 loop and the N-terminus. Mutations that are predicted to disrupt the H-bonding network are shown to produce visual pigments that do not stably bind chromophore and exhibit properties of a misfolded protein. The potential role of a disulfide bond between two conserved Cys residues, Cys(105) in TM3 and Cys(182) in EC2, is necessary for proper folding and trafficking in VCOP. Lastly, certain residues in the EC2 loop are predicted to stabilize the formation of two antiparallel β-strands joined by a hairpin turn, which interact with the chromophore via H-bonding or van der Waals interactions. Mutations of conserved residues result in a decrease in the level of chromophore binding. These results demonstrate that the extracellular loops are crucial for the formation of this cone visual pigment. Moreover, there are significant differences in the structure and function of this region in VCOP compared to that in rhodopsin.

Simplified modeling approach suggests structural mechanisms for constitutive activation of the C5a receptor

Molecular modeling of conformational changes occurring in the transmembrane region of the complement factor 5a receptor (C5aR) during receptor activation was performed by comparing two constitutively active mutants (CAMs) of C5aR, NQ (I124N/L127Q), and F251A, to those of the wild-type C5aR and NQ-N296A (I124N/L127Q/N296A), which have the wild-type phenotype. Modeling involved comprehensive sampling of various rotations of TM helices aligned to the crystal template of the dark-adapted rhodopsin along their long axes. By assuming that the relative energies of the spontaneously activated states of CAMs should be lower or at least comparable to energies characteristic for the ground states, we selected the plausible models for the conformational states associated with constitutive activation in C5aR. The modeling revealed that the hydrogen bonds between the side chains of D82-N119, S85-N119, and S131-C221 characteristic for the ground state were replaced by the hydrogen bonds D82-N296, N296-Y300, and S131-R134, respectively, in the activated states. Also, conformational transitions that occurred upon activation were hindered by contacts between the side chains of L127 and F251. The results rationalize the available data of mutagenesis in C5aR and offer the first specific molecular mechanism for the loss of constitutive activity in NQ-N296A. Our results also contributed to understanding the general structural mechanisms of activation in G-protein-coupled receptors lacking the "ionic lock", R(3.50) and E/D(6.30). Importantly, these results were obtained by modeling approaches that deliberately simplify many elements in order to explore potential conformations of GPCRs involving large-scale molecular movements.

Complement factor 5a receptor chimeras reveal the importance of lipid-facing residues in transport competence

Residues that mediate helix-helix interactions within the seven transmembranes (TM) of G protein-coupled receptors are important for receptor biogenesis and the receptor switch mechanism. By contrast, the residues directly contacting the lipid bilayer have only recently garnered attention as potential receptor dimerization interfaces. In the present study, we aimed to determine the contributions of these lipid-facing residues to receptor function and oligomerization by systemically generating chimeric complement factor 5a receptors in which the entire lipid-exposed surface of a single TM helix was exchanged with the cognate residues from the angiotensin type 1 receptor. Disulfide-trapping and bioluminescence resonance energy transfer (BRET) studies demonstrated robust homodimerization of both complement factor 5a receptor and angiotensin type 1 receptor, but no evidence for heterodimerization. Despite relatively conservative substitutions, the lipid-facing chimeras (TM1, TM2, TM4, TM5, TM6 or TM7) were retained in the endoplasmic reticulum/cis-Golgi network. With the exception of the TM7 chimera that did not bind ligand, the lipid-facing chimeras bound ligand with low affinity, but similar to wild-type complement factor 5a receptors trapped in the endoplasmic reticulum with brefeldin A. These results suggest that the chimeric receptors were properly folded; moreover, native complement factor 5a receptors are not fully competent to bind ligand when present in the endoplasmic reticulum. BRET oligomerization studies demonstrated energy transfer between the wild-type complement factor 5a receptor and the lipid-facing chimeras, suggesting that the lipid-facing residues within a single TM segment are not essential for oligomerization. These studies highlight the importance of the lipid-facing residues in the complement factor 5a receptor for transport competence.

Theoretical study of the human bradykinin-bradykinin B2 receptor complex

The interaction of bradykinin (BK) with the bradykinin B2 receptor (B2R) was analyzed by using molecular modeling (MM) and molecular dynamics (MD) simulations. A homology model for B2R has been generated and the recently determined receptor-bound solid-state NMR spectroscopic structure of BK (Lopez et al., Angew. Chem. 2008, 120, 1692-1695; Angew. Chem. Int. Ed. 2008, 47, 1668-1671) has been modeled into the binding pocket of the receptor to probe the putative ligand-receptor interface. The experimental hormone structure fitted well into the binding pocket of the receptor model and remained stable during the MD simulation. We propose a parallel orientation of the side chains for Arg1 and Arg9 in BK that is bound to B2R. The MD simulation study also allows the conformational changes that lead to the activated form of B2R to be analyzed. The hydrogen bond between N140 (3.35) and W283 (6.48) is the key interaction that keeps the receptor in its inactive form. This hydrogen bond is broken during the MD simulation due to rotation of transmembrane helix 3 (TM3) and is replaced by a new hydrogen bond between W283 (6.48) and N324 (7.45). We propose that this interaction is specific for the activated form of the bradykinin B2 receptor. Additionally, we compared and discussed our putative model in the context of the structural model of the partially activated rhodopsin (Rh*) and with the known biochemical and structural data.

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.

How do receptors activate G proteins?

Heterotrimeric G proteins couple the activation of heptahelical receptors at the cell surface to the intracellular signaling cascades that mediate the physiological responses to extracellular stimuli. G proteins are molecular switches that are activated by receptor-catalyzed GTP for GDP exchange on the G protein alpha subunit, which is the rate-limiting step in the activation of all downstream signaling. Despite the important biological role of the receptor-G protein interaction, relatively little is known about the structure of the complex and how it leads to nucleotide exchange. This chapter will describe what is known about receptor and G protein structure and outline a strategy for assembling the current data into improved models for the receptor-G protein complex that will hopefully answer the question as to how receptors flip the G protein switch.

In silico study of the human rhodopsin and meta rhodopsin II/S-arrestin complexes: impact of single point mutations related to retina degenerative diseases

We propose two models of the human S-arrestin/rhodopsin complex in the inactive dark adapted rhodopsin and meta rhodopsin II form, obtained by homology modeling and knowledge based docking. First, a homology model for the human S-arrestin was built and validated by molecular dynamics, showing an average root mean square deviation difference from the pattern behavior of 0.76 A. Then, combining the human S-arrestin model and the modeled structure of the two human rhodopsin forms, we propose two models of interaction for the human S-arrestin/rhodopsin complex. The models involve two S-arrestin regions related to the N domain (residues 68-78; 170-182) and a third constituent of the C domain (248-253), with the rhodopsin C terminus (330-348). Of the 22 single point mutations related to retinitis pigmentosa and congenital night blindness located in the cytoplasmatic portion of rhodopsin or in S-arrestin, our models locate 16 in the interaction region and relate two others to possible dimer formation. Our calculations also predict that the light activated complex is more stable than the dark adapted rhodopsin and, therefore, of higher affinity to S-arrestin.

Defining the interface between the C-terminal fragment of alpha-transducin and photoactivated rhodopsin

A novel combination of experimental data and extensive computational modeling was used to explore probable protein-protein interactions between photoactivated rhodopsin (R*) and experimentally determined R*-bound structures of the C-terminal fragment of alpha-transducin (Gt(alpha)(340-350)) and its analogs. Rather than using one set of loop structures derived from the dark-adapted rhodopsin state, R* was modeled in this study using various energetically feasible sets of intracellular loop (IC loop) conformations proposed previously in another study. The R*-bound conformation of Gt(alpha)(340-350) and several analogs were modeled using experimental transferred nuclear Overhauser effect data derived upon binding R*. Gt(alpha)(340-350) and its analogs were docked to various conformations of the intracellular loops, followed by optimization of side-chain spatial positions in both R* and Gt(alpha)(340-350) to obtain low-energy complexes. Finally, the structures of each complex were subjected to energy minimization using the OPLS/GBSA force field. The resulting residue-residue contacts at the interface between R* and Gt(alpha)(340-350) were validated by comparison with available experimental data, primarily from mutational studies. Computational modeling performed for Gt(alpha)(340-350) and its analogs when bound to R* revealed a consensus of general residue-residue interactions, necessary for efficient complex formation between R* and its Gt(alpha) recognition motif.

Method To Assess Packing Quality of Transmembrane α-Helices in Proteins. 1. Parametrization Using Structural Data

Integral membrane proteins (MPs) are pharmaceutical targets of exceptional importance. Modern methods of three-dimensional protein structure determination often fail to supply the fast growing field of structure-based drug design with the requested MPs' structures. That is why computational modeling techniques gain a special importance for these objects. Among the principal difficulties limiting application of these methods is the low quality of the MPs' models built in silico. In this series of two papers we present a computational approach to the assessment of the packing "quality" of transmembrane (TM) alpha-helical domains in proteins. The method is based on the concept of protein environment classes, whereby each amino acid residue is described in terms of its environment polarity and accessibility to the membrane. In the first paper we analyze a nonredundant set of 26 TM alpha-helical domains and compute the residues' propensities to five predefined classes of membrane-protein environments. Here we evaluate the proposed approach only by various test sets, cross-validation protocols and ability of the method to delimit the crystal structure of visual rhodopsin, and a number of its erroneous theoretical models. More advanced validation of the method is given in the second article of this series. We assume that the developed "membrane score" method will be helpful in optimizing computer models of TM domains of MPs, especially G-protein coupled receptors.

Predisposition of the dark state of rhodopsin to functional changes in structure

As the only member of the family of G-protein-coupled receptors for which atomic coordinates are available, rhodopsin is widely studied for insight into the molecular mechanism of G-protein-coupled receptor activation. The currently available structures refer to the inactive, dark state, of rhodopsin, rather than the light-activated metarhodopsin II (Meta II) state. A model for the Meta II state is proposed here by analyzing elastic network normal modes in conjunction with experimental data. Key mechanical features and interactions broken/formed in the proposed model are found to be consistent with the experimental data. The model is further tested by using a set of Meta II fluorescence decay rates measured to empirically characterize the deactivation of rhodopsin mutants. The model is found to correctly predict 93% of the experimentally observed effects in 119 rhodopsin mutants for which the decay rates and misfolding data have been measured, including a systematic analysis of Cys-->Ser replacements reported here. Based on the detailed comparison between model and experiments, a cooperative activation mechanism is deduced that couples retinal isomerization to concerted changes in conformation, facilitated by the intrinsic dynamics of rhodopsin. A global hinge site is identified near the retinal-binding pocket that ensures the efficient propagation of signals from the central transmembrane region to both cytoplasmic and extracellular ends. The predicted activation mechanism opens the transmembrane helices at the critical G-protein binding cytoplasmic domain. This model provides a detailed, mechanistic description of the activation process, extending experimental observations and yielding new insights for further tests.

Docking studies of agonists and antagonists suggest an activation pathway of the A3 adenosine receptor

Structural determinants of ligand efficacy in the human A(3) adenosine receptor (AR) were studied using pharmacophore and docking analyses of various categories of A(3) selective ligands: inverse agonist, neutral antagonist (nonnucleoside and nucleoside), and agonist (partial and full). The homology modeling of GPCRs was adapted to provide two templates: the rhodopsin-based resting state for antagonist binding and a putative Meta I state, conformationally altered at a key residue (W6.48), for agonist binding. The preferential binding domains and/or local conformational changes associated with docking of three high affinity A(3)AR ligands were compared: inverse agonist PSB-11 1 ((R)-8-ethyl-4-methyl-2-phenyl-imidazo[2,1-i]purin-5-one); neutral antagonist MRE-3008F20 7 (5-[[(4-methoxyphenyl)amino]carbonyl]amino-8-methyl-2-(2-furyl)pyra-zolo[4,3-e]1,2,4-triazolo[1,5-c]pyrimidine), and full agonist Cl-IB-MECA 21 (2-chloro-N(6)-(3-iodobenzyl)-5'-N-methylcarboxamidoadenosine) to define a distinct recognition mode for each. Ribose-containing agonists were more hydrophilic than nonnucleoside antagonists, and H-bonding ability at the ribose 3'- and 5'-positions was required for agonism. From the receptor perspective, common requirements for activation included the destabilization of H-bond networks at W6.48 and H7.43, the specific interactions of the ribose moiety in its putative hydrophilic pocket at T3.36, S7.42, and H7.43, the stabilization of the complex by inward movement of F5.43, and the characteristic rotation of W6.48. By analogy, outward rotation of the W6.48 side-chain upon activation of an internally-crosslinking mutant M(3) muscarinic receptor was indicated by constrained molecular dynamics (MD). Our results are consistent with an anti-clockwise rotation (from the extracellular view) of transmembrane domains 3, 5, 6, and 7, as proposed for other Family A GPCRs. Thus, the putative conformational changes associated with A(3)AR activation indicate a shared mechanism of GPCR activation similar to rhodopsin.

A comprehensive structure-function map of the intracellular surface of the human C5a receptor. I. Identification of critical residues

G protein-coupled receptors are one of the largest protein families in nature; however, the mechanisms by which they activate G proteins are still poorly understood. To identify residues on the intracellular face of the human C5a receptor that are involved in G protein activation, we performed a genetic analysis of each of the three intracellular loops and the carboxyl-terminal tail of the receptor. Amino acid substitutions were randomly incorporated into each loop, and functional receptors were identified in yeast. The third intracellular loop contains the largest number of preserved residues (positions resistant to amino acid substitutions), followed by the second loop, the first loop, and lastly the carboxyl terminus. Surprisingly, complete removal of the carboxyl-terminal tail did not impair C5a receptor signaling. When mapped onto a three-dimensional structural model of the inactive state of the C5a receptor, the preserved residues reside on one half of the intracellular surface of the receptor, creating a potential activation face. Together these data provide one of the most comprehensive functional maps of the intracellular surface of any G protein-coupled receptor to date.

Interactions between Conserved Residues in Transmembrane Helices 2 and 7 during Angiotensin AT1Receptor Activation

Site-directed mutagenesis studies and independent molecular modeling studies were combined to investigate the network of inter-residue interactions within the transmembrane region of the angiotensin AT(1a) receptor. Site-directed mutagenesis was focused on residues Tyr292, Asn294, Asn295, and Asn298 in transmembrane helix 7, and the conserved Asp74 in helix 2 and other polar residues. Functional interactions between pairs of residues were evaluated by determining the effects of single and double-reciprocal mutations on agonist-induced AT(1a) receptor activation. Replacement of Tyr292 by aspartate in helix 7 abolished radioligand binding to both Y292D and D74Y/Y292D mutant receptors. Reciprocal mutations of Asp74/Asn294, Ser115/Asn294, Ser252/Asn294, and Asn298/Sen115 caused additive impairment of function, suggesting that these pairs of residues make independent contributions to AT(1a) receptor activation. In contrast, mutations of the Asp74/Tyr298 pair revealed that the D74N/N298D reciprocal mutation substantially increased the impaired inositol phosphate responses of the D74N and N298D receptors. Extensive molecular modeling yielded 3D models of the TM region of the AT(1) receptor and the mutants as well as of their complexes with angiotensin II, which were used to rationalize the possible reasons of impairing of function of some mutants. These data indicate that Asp74 and Asn298 are not optimally positioned for direct strong interaction in the resting conformation of the AT(1a) receptor. Balance of interactions between residues in helix 2 (as D74) and helix 7 (as N294, N295 and N298) in the AT(1) receptors, however, has a crucial role both in determining their functional activity and levels of their expression.

G-protein coupled receptor structure

Because of their central role in regulation of cellular function, structure/function relationships for G-protein coupled receptors (GPCR) are of vital importance, yet only recently have sufficient data been obtained to begin mapping those relationships. GPCRs regulate a wide range of cellular processes, including the senses of taste, smell, and vision, and control a myriad of intracellular signaling systems in response to external stimuli. Many diseases are linked to GPCRs. A critical need exists for structural information to inform studies on mechanism of receptor action and regulation. X-ray crystal structures of only one GPCR, in an inactive state, have been obtained to date. However considerable structural information for a variety of GPCRs has been obtained using non-crystallographic approaches. This review begins with a review of the very earliest GPCR structural information, mostly derived from rhodopsin. Because of the difficulty in crystallizing GPCRs for X-ray crystallography, the extensive published work utilizing alternative approaches to GPCR structure is reviewed, including determination of three-dimensional structure from sparse constraints. The available X-ray crystallographic analyses on bovine rhodopsin are reviewed as the only available high-resolution structures for any GPCR. Structural information available on ligand binding to several receptors is included. The limited information on excited states of receptors is also reviewed. It is concluded that while considerable basic structural information has been obtained, more data are needed to describe the molecular mechanism of activation of a GPCR.

Comparison of the retinitis pigmentosa mutations in rhodopsin with a functional map of the C5a receptor

We compare the known retinitis pigmentosa (RP) mutations in rhodopsin with mutational data obtained for the complement factor 5a receptor (C5aR), a member of the rhodopsin-like family of G protein-coupled receptors (GPCRs). We have performed genetic analyses that define residues that are required for C5aR folding and function. The cognate residues in rhodopsin are not preferentially mutated in RP, suggesting that the predominant molecular defect in RP involves more than simple misfolding or inactivation. Energy calculations are performed to elucidate the structural effects of the RP mutations. Many of these mutations specifically disrupt the environment of the retinal prosthetic group of rhodopsin, and these do not correspond to essential residues in C5aR. This may be because a retinal group is present in rhodopsin but not in C5aR. Another subset of RP mutations is more generally important for receptor structure, and these mutations correlate with essential residues of C5aR.

Modeling activated states of GPCRs: the rhodopsin template

Activation of G Protein-Coupled Receptors (GPCRs) is an allosteric mechanism triggered by ligand binding and resulting in conformational changes transduced by the transmembrane domain. Models of the activated forms of GPCRs have become increasingly necessary for the development of a clear understanding of signal propagation into the cell. Experimental evidence points to a multiplicity of conformations related to the activation of the receptor, rendered important physiologically by the suggestion that different conformations may be responsible for coupling to different signaling pathways. In contrast to the inactive state of rhodopsin (RHO) for which several high quality X-ray structures are available, the structure-related information for the active states of rhodopsin and all other GPCRs is indirect. We have collected and stored such information in a repository we maintain for activation-specific structural data available for rhodopsin-like GPCRs, http://www.physiology.med.cornell.edu/GPCRactivation/gpcrindex.html . Using these data as structural constraints, we have applied Simulated Annealing Molecular Dynamics to construct a number of different active state models of RHO starting from the known inactive structure. The common features of the models indicate that TM3 and TM5 play an important role in activation, in addition to the well-established rearrangement of TM6. Some of the structural changes observed in these models occur in regions that were not involved in the constraints, and have not been previously tested experimentally; they emerge as interesting candidates for further experimental exploration of the conformational space of activated GPCRs. We show that none of the normal modes calculated from the inactive structure has a dominant contribution along the path of conformational rearrangement from inactive to the active forms of RHO in the models. This result may differentiate rhodopsin from other GPCRs, and the reasons for this difference are discussed in the context of the structural properties and the physiological function of the protein.

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.

Equilibrium between Metarhodopsin-I and Metarhodopsin-II Is Dependent on the Conformation of the Third Cytoplasmic Loop

Rhodopsin is a G-protein-coupled receptor (GPCR) that is the light detector in the rod cells of the eye. Rhodopsin is the best understood member of the large GPCR superfamily and is the only GPCR for which atomic resolution structures have been determined. However, these structures are for the inactive, dark-adapted form. Characterization of the conformational changes in rhodopsin caused by light-induced activation is of wide importance, because the metarhodopsin-II photoproduct is analogous to the agonist-occupied conformation of other GPCRs, and metarhodopsin-I may be similar to antagonist-occupied GPCR conformations. In this work we characterize the interaction of antibody K42-41L with the metarhodopsin photoproducts. K42-41L is shown to inhibit formation of metarhodopsin-II while it stabilizes the metarhodopsin-I state. Thus, K42-41L recognizes an epitope accessible in dark-adapted rhodopsin and metarhodopsin-I that is lost upon formation of metarhodopsin-II. Previous work has shown that the peptide TGALQERSK is able to mimic the K42-41L epitope, and we have now determined the structure of the K42-41L-peptide complex. The structure demonstrates a central role for elements of the rhodopsin C3 loop, particularly Gln238 and Glu239, in the interaction with K42-41L. Geometric constraints taken from the antibody-bound peptide were used to model the epitope on the rhodopsin surface. The resulting model suggests that K42-41L locks the C3 loop into an extended conformation that is intermediate between two compact conformations seen in crystal structures of dark-adapted rhodopsin. Together, the structural and functional data strongly suggest that the equilibrium between metarhodopsin-I and metarhodopsin-II is dependent upon the conformation of the C3 loop. The biological implications of this model and its possible relations to dimeric and multimeric complexes of rhodopsin are discussed.

Adenosine receptor modelling. A1/A2a selectivity

Three-dimensional models of the A(1) and A(2a) adenosine receptors (AR) were constructed by means of a homology procedure, using bovine rhodopsin as a template. In order to validate the two models, a docking analysis of selective agonists was carried out. The study shows that A(1)/A(2a) selectivity is mainly influenced by the different ability of the two receptors to give lipophilic interactions, instead of giving different H bonds. The binding site cavity of the A(1)AR is smaller than that of the A(2a)AR, and for this reason, less bulky ligands like CPA are able to give close interactions with the A(1)AR, unlike larger ligands such as CGS-21680. The different dimensions of the binding site cavity could be due to the presence of three residues of proline, which cause a different rearrangement of the TM, thus modifying the side chain disposition inside the inter-helix channel.

Bridging the gaps in 3D structure of the inositol 1,4,5-trisphosphate receptor-binding core

Calcium concentration is strictly regulated in all cells. The inositol 1,4,5-trisphosphate receptor (IP(3)R), which forms a homotetrameric Ca2+ release channel in the endoplasmic reticulum, is one of the key molecules responsible for this regulation. The opening of this channel requires binding of two intracellular messengers, which are inositol 1,4,5-trisphosphate (IP(3)) and Ca2+. To promote the Ca2+-channel gating and release from the endoplasmic reticulum, IP(3) binds to the amino-terminal region of IP(3)R. Recently, the crystal structure of IP(3)R-binding core in complex with its ligand was presented [I. Bosanac, J.R. Alattia, T.K. Mai, J. Chan, S. Talarico, F.K. Tong, K.I. Tong, F. Yoshikawa, T. Furuichi, M. Iwai, T. Michikawa, K. Mikoshiba, M. Ikura, Structure of the inositol 1,4,5-trisphosphate receptor binding core in complex with its ligand, Nature 420 (2002) 696-700; I. Bosanac, H. Yamazaki, T. Matsu-ura, T. Michikawa, K. Mikoshiba, M. Ikura, Crystal structure of the ligand-binding suppressor domain of type 1 inositol 1,4,5-trisphosphate receptor, Mol. Cell 17 (2005) 193-203]. The space positions of residues 289-301 (segment A), 320-350 (segment B), 373-386 (segment C), and 529-545 (segment D) were not determined by the X-ray crystallography. To bridge these gaps, the computer modeling of physiologically meaningful low-energy 3D structures of the segments A-D of the inositol 1,4,5-trisphosphate receptor has been carried out by using a hierarchical conformational search algorithm combining two approaches: knowledge-based homology modeling and ab initio conformational search strategy. The structure analysis suggests a Ca2+-binding site of high affinity formed by residues 296-335, several low-energy regular secondary structure units within the segment B, and a number of hinge regions within the segments A-D, important for the receptor functioning.

Computational studies of the primary phototransduction event in visual rhodopsin

This Account addresses recent advances in the elucidation of the detailed molecular rearrangements due to the primary photochemical event in rhodopsin, a prototypical G-protein-coupled receptor (GPCR) responsible for the signal transmission cascade in the vertebrate vision process. The reviewed studies provide fundamental insight on long-standing problems regarding the assembly and function of the individual residues and bound water molecules that form the rhodopsin active site, a center that catalyzes the 11-cis/all-trans isomerization of the retinyl chromophore in the primary step of the phototransduction mechanism. Emphasis is placed on the authors' recent computational studies, based on state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, addressing the structural refinement of the retinyl chromophore binding site in high-resolution X-ray structures of bovine visual rhodopsin, the energy storage mechanism, and the molecular origin of spectroscopic changes due to the primary photochemical event.

Molecular mechanisms of constitutive activity: mutations at position 111 of the angiotensin AT1 receptor

A possible molecular mechanism for the constitutive activity of mutants of the angiotensin type 1 receptor (AT1) at position 111 was suggested by molecular modeling. This involves a cascade of conformational changes in spatial positions of side chains along transmembrane helix (TM3) from L112 to Y113 to F117, which in turn, results in conformational changes in TM4 (residues I152 and M155) leading to the movement of TM4 as a whole. The mechanism is consistent with the available data of site-directed mutagenesis, as well as with correct predictions of constitutive activity of mutants L112F and L112C. It was also predicted that the double mutant N111G/L112A might possess basal constitutive activity comparable with that of the N111G mutant, whereas the double mutants N111G/Y113A, N111G/F117A, and N111G/I152A would have lower levels of basal activity. Experimental studies of the above double mutants showed significant constitutive activity of N111G/L112A and N111G/F117A. The basal activity of N111G/I152A was higher than expected, and that of N111G/Y113A was not determined due to poor expression of the mutant. The proposed mechanism of constitutive activity of the AT(1) receptor reveals a novel nonsimplistic view on the general problem of constitutive activity, and clearly demonstrates the inherent complexity of the process of G protein-coupled receptor (GPCR) activation.

Modeling flexible loops in the dark-adapted and activated states of rhodopsin, a prototypical G-protein-coupled receptor

Conformational possibilities of flexible loops in rhodopsin, a prototypical G-protein-coupled receptor, were studied by modeling both in the dark-adapted (R) and activated (R*) states. Loop structures were built onto templates representing the R and R* states of the TM region of rhodopsin developed previously (G. V. Nikiforovich and G. R. Marshall. 2003. Biochemistry. 42:9110). Geometrical sampling and energy calculations were performed for each individual loop, as well as for the interacting intracellular loops IC1, IC2, and IC3 and the extracellular loops EC1, EC2, and EC3 mounted on the R and R* templates. Calculations revealed that the intra- and extracellular loops of rhodopsin possess low-energy structures corresponding to large conformational movements both in the R and R* states. Results of these calculations are in good agreement with the x-ray data available for the dark-adapted rhodopsin as well as with the available experimental biophysical data on the disulfide-linked mutants of rhodopsin. The calculated results are used to exemplify how the combined application of the results of independent calculations with emerging experimental data can be used to select plausible three-dimensional structures of the loops in rhodopsin.

Molecular Characterization of a Purified 5-HT4 Receptor

Serotonin 5-HT(4(a)) receptor, a G-protein-coupled receptor (GPCR), was produced as a functional isolated protein using Escherichia coli as an expression system. The isolated receptor was characterized at the molecular level by circular dichroism (CD) and steady-state fluorescence. A specific change in the near-UV CD band associated with the GPCR disulfide bond connecting the third transmembrane domain to the second extracellular loop (e2) was observed upon agonist binding to the purified receptor. This is a direct experimental evidence for a change in the conformation of the e2 loop upon receptor activation. Different variations were obtained depending whether the ligand was an agonist (partial or full) or an inverse agonist. In contrast, antagonist binding did not induce any variation. These observations provide a first direct evidence for the fact that free (or antagonist-occupied), active (partial- or full agonist-occupied) and silent (inverse agonist-occupied) states of the receptor involve different arrangements of the e2 loop. Finally, ligand-induced changes in the fluorescence emission profile of the purified receptor confirmed that the partial agonist stabilized a single, well-defined, conformational state and not a mixture of different states. This result is of particular interest in a pharmacological perspective since it directly demonstrates that the efficacy of a drug is likely due to the stabilization of a ligand-specific state rather than selection of a mixture of different conformational states of the receptor.

Architecture of the Human Urotensin II Receptor: Comparison of the Binding Domains of Peptide and Non-Peptide Urotensin II Agonists

The human urotensin II receptor (h-UTR) is a member of the family of rhodopsin-like G-protein-coupled receptors (GPCRs) involved in the modulation of the functionality of many tissues and organs. Recently the urotensin-II (UII) neuropeptide, which is a potent vasoconstrictor in mammals and it is postulated to play a central role in cardiovascular homeostasis, has been identified as an agonist of the UII receptor. To elucidate the receptor's molecular recognition, a h-UTR model was constructed by homology modeling using the 2.6 A crystal structure of bovine rhodopsin as a template and subsequently refined by molecular dynamics simulations. The molecular recognition of h-UTR was probed by automated docking of P5U, a potent UII peptide agonist, as well as of the non-peptide compounds 1-4. We believe that this new model of the h-UTR provides the means for understanding the ligand's potency and for facilitating the design of novel and more potent UII ligands.

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.

Random mutagenesis of the M3 muscarinic acetylcholine receptor expressed in yeast: identification of second-site mutations that restore function to a coupling-deficient mutant M3 receptor

The M(3) muscarinic receptor is a prototypical member of the class A family of G protein-coupled receptors (GPCRs). To gain insight into the structural mechanisms governing agonist-mediated M(3) receptor activation, we recently developed a genetically modified yeast strain (Saccharomyces cerevisiae) which allows the efficient screening of large libraries of mutant M(3) receptors to identify mutant receptors with altered/novel functional properties. Class A GPCRs contain a highly conserved Asp residue located in transmembrane domain II (TM II; corresponding to Asp-113 in the rat M(3) muscarinic receptor) which is of fundamental importance for receptor activation. As observed previously with other GPCRs analyzed in mammalian expression systems, the D113N point mutation abolished agonist-induced receptor/G protein coupling in yeast. We then subjected the D113N mutant M(3) receptor to PCR-based random mutagenesis followed by a yeast genetic screen to recover point mutations that can restore G protein coupling to the D113N mutant receptor. A large scale screening effort led to the identification of three such second-site suppressor mutations, R165W, R165M, and Y250D. When expressed in the wild-type receptor background, these three point mutations did not lead to an increase in basal activity and reduced the efficiency of receptor/G protein coupling. Similar results were obtained when the various mutant receptors were expressed and analyzed in transfected mammalian cells (COS-7 cells). Interestingly, like Asp-113, Arg-165 and Tyr-250, which are located at the cytoplasmic ends of TM III and TM V, respectively, are also highly conserved among class A GPCRs. Our data suggest a conformational link between the highly conserved Asp-113, Arg-165, and Tyr-250 residues which is critical for receptor activation.

Novel amino acid derived natural products from the ascidian Atriolum robustum: identification and pharmacological characterization of a unique adenosine derivative

Investigation of the methanolic extract of the Australian ascidian Atriolum robustum led to the isolation and characterization of five new amino acid derived structures (1-5). The structures were elucidated employing spectroscopic techniques (NMR, MS, UV, and IR). The absolute stereochemistry of 1 and 2 was established by chemical degradation, derivatization, and chiral GC-MS analysis. Structures 4 and 5 are complex nucleosides containing rare methylthioadenosine and methylsulfinyladenosine moieties, respectively. In radioligand binding studies the 5'-deoxy-5'-methylthioadenosine-2',3'-diester 4 exhibited affinity for A(1) and A(3) adenosine receptors with K(i) values below 10 microM. Its affinity was somewhat lower for A(2A) (K(i) = 17 microM) and much lower for A(2B) adenosine receptors. Analytical experiments using capillary electrophoresis showed that compound 4 was stable under the conditions of radioligand binding studies. Incubation with carboxylesterase resulted in slow hydrolysis of the adenosine derivative to 5'-deoxy-5'-methylthioadenosine (MTA), which was about 10-fold more potent at adenosine receptors than compound 4. Thus, the 2',3'-diester derivative 4 may act as a lipophilic prodrug of MTA in addition to its own adenosine receptor activity. GTP shift experiments indicated that the adenosine derivative was a partial agonist at A(1) adenosine receptors of rat brain cortical membranes. Compound 4 inhibited cAMP accumulation in Chinese hamster ovary (CHO) cell membranes recombinantly expressing the human A(3) adenosine receptor, thus indicating that the adenosine derivative also acted as a partial agonist at A(3)ARs. Homology models of the A(1) and the A(3) adenosine receptors in their putative active and inactive conformations were built and used for docking of the sterically demanding compound 4. It was found that this ligand fit well into the binding pockets of both receptor subtypes because of its highly flexible structure, although in somewhat different binding modes.