Third Extracellular Loop (EC3)-N Terminus Interaction Is Important for Seven-transmembrane Domain Receptor Function: IMPLICATIONS FOR AN ACTIVATION MICROSWITCH REGION*

The canonical heptahelical bundle architecture of seven-transmembrane domain (7TM) receptors is intertwined by three intra- and three extracellular loops, whose local conformations are important in receptor signaling. Many 7TM receptors contain a cysteine residue in the third extracellular loop (EC3) and a complementary cysteine residue on the N terminus. The functional role of such EC3-N terminus conserved cysteine pairs remains unclear. This study explores the role of the EC3-N terminus cysteine pairs on receptor conformation and G protein activation by disrupting them in the chemokine receptor CXCR4, while engineering a novel EC3-N terminus cysteine pair into the complement factor 5a receptor (C5aR), a chemo attractant receptor that lacks it. Mutated CXCR4 and C5aRs were expressed in engineered yeast. Mutation of the cysteine pair with the serine pair (C28S/C274S) in constitutively active mutant CXCR4 abrogated the receptor activation, whereas mutation with the aromatic pair (C28F–C274F) or the salt bridge pair (C28R/C274E), respectively, rescued or retained the receptor activation in response to CXCL12. In this context, the cysteine pair (Cys³⁰ and Cys²⁷²) engineered into the EC3-N terminus (Ser³⁰ and Ser²⁷²) of a novel constitutively active mutant of C5aR restrained the constitutive signaling without affecting the C5a-induced activation. Further mutational studies demonstrated a previously unappreciated role for Ser²⁷² on EC3 of C5aR and its interaction with the N terminus, thus defining a new microswitch region within the C5aR. Similar results were obtained with mutated CXCR4 and C5aRs expressed in COS-7 cells. These studies demonstrate a novel role of the EC3-N terminus cysteine pairs in G protein-coupled receptor activation and signaling.     

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.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Role of the second extracellular loop of human C3a receptor in agonist binding and receptor function

The C3a anaphylatoxin receptor (C3aR) is a G protein-coupled receptor with an unusually large second extracellular loop (e2 loop, approximately 172 amino acids). To determine the function of this unique structure, chimeric and deletion mutants were prepared and analyzed in transfected RBL-2H3 cells. Whereas replacement of the C3aR N-terminal segment with that from the human C5a receptor had minimal effect on C3a binding, substitution of the e2 loop with a smaller e2 loop from the C5a receptor (C5aR) abolished binding of 125I-C3a and C3a-stimulated calcium mobilization. However, as much as 65% of the e2 loop sequence (amino acids 198-308) may be removed without affecting C3a binding or calcium responses. The e2 loop sequences adjacent to the transmembrane domains contain multiple aspartate residues and are found to play an important role in C3a binding based on deletion mutagenesis. Replacement of five aspartate residues in the e2 loop with lysyl residues significantly compromised both the binding and functional capabilities of the C3a receptor mediated by intact C3a or by two C3a analog peptides. These data suggest a two-site C3a-C3aR interaction model similar to that established for C5a/C5aR. The anionic residues near the N and C termini of the C3aR e2 loop constitute a non-effector secondary interaction site with cationic residues in the C-terminal helical region of C3a, whereas the C3a C-terminal sequence LGLAR engages the primary effector site in C3aR.     

On the spatial disposition of the fifth transmembrane helix and the structural integrity of the transmembrane binding site in the opioid and ORL1 G protein-coupled receptor family

Evidence from statistical cluster analyses of a multiple sequence alignment of G protein-coupled receptor seven-helix folds supports the existence of structurally conserved transmembrane (TM) ligand binding sites in the opioid/opioid receptor-like (ORL1) and amine receptor families. Based on the expectation that functionally conserved regions in homologous proteins will display locally higher levels of sequence identity compared with global sequence similarities that pertain to the overall fold, this approach may have wider applications in functional genomics to annotate sequence data. Binding sites in models of the kappa-opioid receptor seven-helix bundle built from the rhodopsin templates of Baldwin et al. (1997) [J. Mol. Biol., 272, 144-164] and Herzyk and Hubbard (1998) [J. Mol. Biol., 281, 742-751] are compared. The Herzyk and Hubbard template is found to be in better accord with experimental studies of amine, opioid and rhodopsin receptors owing to the reduced physical separation of the extracellular parts of TM helices V and VI and differences in the rotational orientation of the N-terminal of helix V that reveal side chain accessibilities in the Baldwin et al. structure to be out of phase with relative alkylation rates of engineered cysteine residues in the TM binding site of the alpha(2A)-adrenergic receptor. TM helix V in the Baldwin et al. template has been remodelled with a different proline kink to satisfy experimental constraints. A recent proposal that rotation of helix V is associated with receptor activation is critically discussed.     

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.     

Modeling molecular mechanisms of binding of the anaphylatoxin C5a to the C5a receptor

This study presents the 3D model of the complex between the anaphylatoxin C5a and its specific receptor, C5aR. This is the first 3D model of a G-protein-coupled receptor (GPCR) complex with a peptide ligand deduced by a molecular modeling procedure analyzing various conformational possibilities of the extracellular loops and the N-terminal segment of the GPCR. The modeling results indicated two very different ways of interacting between C5a and C5aR at the two interaction sites suggested earlier based on the data of site-directed mutagenesis. Specifically, C5a and C5aR can be involved in "mutual-induced fit", where the interface between the molecules is determined by both the receptor and the ligand. The rigid core of the C5a ligand selects the proper conformations of the highly flexible N-terminal segment of C5aR (the first interaction site). At the same time, the binding conformation of the flexible C-terminal fragment of C5a is selected by well-defined interactions with the TM region of the C5aR receptor (the second interaction site). The proposed 3D model of C5a/C5aR complex was built without direct use of structural constraints derived from site-directed mutagenesis reserving those data for validation of the model. The available data of site-directed mutagenesis of C5a and C5aR were successfully rationalized with the help of the model. Also, the modeling results predicted that the full-length C5a and C5a-des74 metabolite would have different binding modes with C5aR. Modeling approaches employed in this study are readily applicable for studies of molecular mechanisms of binding of other polypeptide ligands to their specific GPCRs.     

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

Determinants of Neutralization Resistance in the Envelope Glycoproteins of a Simian-Human Immunodeficiency Virus Passaged In Vivo

In vivo passage of a simian-human immunodeficiency virus (SHIV-89.6) generated a virus, SHIV-89.6P, that exhibited increased resistance to some neutralizing antibodies (G. B. Karlsson et al., J. Exp. Med. 188:1159-1171, 1998). Here we examine the range of human immunodeficiency virus type 1 (HIV-1) neutralizing antibodies to which the passaged virus became resistant and identify envelope glycoprotein determinants of antibody resistance. Compared with the envelope glycoproteins derived from the parental SHIV-89.6, the envelope glycoproteins of the passaged virus were resistant to antibodies directed against the gp120 V3 variable loop and the CD4 binding site. By contrast, both viral envelope glycoproteins were equally sensitive to neutralization by two antibodies, 2G12 and 2F5, that recognize poorly immunogenic structures on gp120 and gp41, respectively. Changes in the V2 and V3 variable loops of gp120 were necessary and sufficient for full resistance to the IgG1b12 antibody, which is directed against the CD4 binding site. Changes in the V3 loop specified complete resistance to a V3 loop-directed antibody, while changes in the V1/V2 loops conferred partial resistance to this antibody. The epitopes of the neutralizing antibodies were not disrupted by the resistance-associated changes. These results indicate that in vivo selection occurs for HIV-1 envelope glycoproteins with variable loop conformations that restrict the access of antibodies to immunogenic neutralization epitopes.     

Identification of ligand effector binding sites in transmembrane regions of the human G protein-coupled C3a receptor

The human C3a anaphylatoxin receptor (C3aR) is a G protein-coupled receptor (GPCR) composed of seven transmembrane alpha-helices connected by hydrophilic loops. Previous studies of chimeric C3aR/C5aR and loop deletions in C3aR demonstrated that the large extracellular loop2 plays an important role in noneffector ligand binding; however, the effector binding site for C3a has not been identified. In this study, selected charged residues in the transmembrane regions of C3aR were replaced by Ala using site-directed mutagenesis, and mutant receptors were stably expressed in the RBL-2H3 cell line. Ligand binding studies demonstrated that R161A (helix IV), R340A (helix V), and D417A (helix VII) showed no binding activity, although full expression of these receptors was established by flow cytometric analysis. C3a induced very weak intracellular calcium flux in cells expressing these three mutant receptors. H81A (helix II) and K96A (helix III) showed decreased ligand binding activity. The calcium flux induced by C3a in H81A and K96A cells was also consistently reduced. These findings suggest that the charged transmembrane residues Arg161, Arg340, and Asp417 in C3aR are essential for ligand effector binding and/or signal coupling, and that residues His81 and Lys96 may contribute less directly to the overall free energy of ligand binding. These transmembrane residues in C3aR identify specific molecular contacts for ligand interactions that account for C3a-induced receptor activation.     

Understanding the thermodynamic stability of an RNA hairpin and its mutant.

The thermodynamic stability of RNA hairpin loops has been a subject of considerable interest in the recent past (Wimberly et al., 1991). There have been experimental reports indicating that the hairpins with a C(UUCG)G loop sequence are thermodynamically very stable (Wimberly et al., 1991). We used the solution structure of GGAC(UUCG)GUCC (Cheong et al., 1990; Varani et al., 1991) as the starting conformation in our attempt to understand its thermodynamic stability. We carried out molecular dynamics/free energy simulations to understand the basis for the destabilization of the C(UUCG)G loop by mutating cytosine (C7)-->uracil. Because of the limited length of simulation and the presence of kinetic barriers (solvent intervention) to the uracil-->cytosine mutation, all of our computed free energy differences are based on multiple forward simulations. Based on these calculations we find that the cytosine-->uracil mutation in the loop destabilizes it by approximately 1.5kcal/mol relative to that of the reference state, an A-form RNA but with cytosine (C7) looped out. This is the same sign and magnitude as that observed in the thermodynamic studies carried out by Varani et al.(1991). We have carried out free energy component analysis to understand the effect of mutating the cytosine residue to uracil on the thermodynamic stability of the C(UUCG)G hairpin loops. Our calculations show that the most significant contribution to the stability is from the phosphate group linking U5 and U6, which favors the cytosine residue over uracil by about 6.0 kcal/mol. The residues U5, U6, and G8 in the loop region also contribute significantly to the stability. The contributions from the salt and solvent compensate each other, indicating the dynamic nature of interactions of the environment with the nucleic acid system and the coupling between these two components.     

Light-driven activation of beta 2-adrenergic receptor signaling by a chimeric rhodopsin containing the beta 2-adrenergic receptor cytoplasmic loops

Structure-function studies of rhodopsin indicate that both intradiscal and transmembrane (TM) domains are required for retinal binding and subsequent light-induced structural changes in the cytoplasmic domain. Further, a hypothesis involving a common mechanism for activation of G-protein-coupled receptor (GPCR) has been proposed. To test this hypothesis, chimeric receptors were required in which the cytoplasmic domains of rhodopsin were replaced with those of the beta(2)-adrenergic receptor (beta(2)-AR). Their preparation required identification of the boundaries between the TM domain of rhodopsin and the cytoplasmic domain of the beta(2)-AR necessary for formation of the rhodopsin chromophore and its activation by light and subsequent optimal activation of beta(2)-AR signaling. Chimeric receptors were constructed in which the cytoplasmic loops of rhodopsin were replaced one at a time and in combination. In these replacements, size of the third cytoplasmic (EF) loop critically determined the extent of chromophore formation, its stability, and subsequent signal transduction specificity. All the EF loop replacements showed significant decreases in transducin activation, while only minor effects were observed by replacements of the CD and AB loops. Light-dependent activation of beta(2)-AR leading to Galphas signaling was observed only for the EF2 chimera, and its activation was further enhanced by replacements of the other loops. The results demonstrate coupling between light-induced conformational changes occurring in the transmembrane domain of rhodopsin and the cytoplasmic domain of the beta(2)-AR.     

Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism

The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.     

The second cytoplasmic loop of metabotropic glutamate receptor functions at the third loop position of rhodopsin.

G protein-coupled receptors identified so far are classified into at least three major families based on their amino acid sequences. For the family of receptors homologous to rhodopsin (family 1), the G protein activation mechanism has been investigated in detail, but much less for the receptors of other families. To functionally compare the G protein activation mechanism between rhodopsin and metabotropic glutamate receptor (mGluR), which belong to distinct families, we prepared a set of bovine rhodopsin mutants whose second or third cytoplasmic loop was replaced with either the second or third loop of Gi/Go- or Gq-coupled mGluR (mGluR6 or mGluR1). Among these mutants, the mutants in which the second or third loop was replaced with the corresponding loop of mGluR exhibited no G protein activation ability. In contrast, the mutant whose third loop was replaced with the second loop of Gi/Go-coupled mGluR6 efficiently activated Gi but not Gt: this activation profile is almost identical with those of the mutant rhodopsins whose third loop was replaced with those of the Gi/Go-coupled receptors in family 1 [Yamashita et al. (2000) J. Biol. Chem. 275, 34272-34279]. The mutant whose third loop was replaced with the second loop of Gq-coupled mGluR1 partially retained the Gi coupling ability of rhodopsin, which is in contrast to the fact that all the rhodopsin mutants having the third loops of Gq-coupled receptors in family 1 exhibit no detectable Gi activation. These results strongly suggest that the molecular architectures of rhodopsin and mGluR are different, although the G protein activation mechanism involving the cytoplasmic loops is common.     

Structural studies on site-directed mutants of domain 3 of Xenopus laevis oocyte 5 S ribosomal RNA

Base substitutions have been introduced into the highly conserved sequences of loops D and E within domain 3 of Xenopus laevis oocyte 5 S rRNA. The effects of these mutations on the solution structure of this 5 S rRNA have been studied by means of probing with nucleases, and with chemical reagents under native and semi-denaturing conditions. The data obtained with these mutants support the graphic model of Xenopus oocyte 5 S rRNA proposed by Westhof et al. In particular, our results rule out the existence of long-range base-pairing interactions between loop C and either loop D or loop E. The data also confirm that loops D and E in the wild-type 5 S RNA adopt unusual secondary structures and illustrate the importance of nucleotide sequence in the formation of intrinsic local loop conformations via non-canonical base-pairs and specific base-phosphate contacts. Consistent with this conclusion is our observation that the domain 3 fragment of Xenopus oocyte 5 S rRNA adopts the same conformation as the corresponding region in the full-length 5 S rRNA.     

Vasopressin V2 Receptor/Bioligand Interactions

We predict some essential interactions between the V2 vasopressin renal receptor (V2R) and its agonists [Arg⁸]vasopressin (AVP) and [D-Arg⁸]vasopressin (DAVP), and the non-peptide antagonist OPC-31260. V2R controls antidiuresis and belongs to the superfamily of heptahelical transmembrane (7TM) G-protein-coupled receptors (GPCRs). The receptor was built, the ligands were docked and the structures relaxed using advanced molecular modeling techniques. Docked agonists and antagonists appear to prefer similar V2R compartments. A number of receptor amino acid residues are indicated, mainly in the TM3–TM7 helices, as potentially important in ligand binding. Many of these residues are invariant for either the GPCR superfamily or the subfamily of related (vasopressin V2R, V1aR and V1bR and oxytocin OR) receptors. Moreover, some of the equivalent residues in V1aR have already been found critical for ligand affinity [Mouillac et al., J. Biol. Chem., 270 (1995) 25771].     

The conformational and dynamic basis for ligand binding reactivity in hemoglobin Ypsilanti (beta 99 asp-->Tyr): origin of the quaternary enhancement effect

Hemoglobin Ypsilanti (HbY) is a stable tetrameric hemoglobin that binds oxygen with little or no cooperativity and with high affinity [Doyle, M. L., et al. (1992) Proteins: Struct., Funct., Genet. 14, 351-362]. It displays an especially large quaternary enhancement effect. An X-ray crystallographic study [Smith, F. R., et al. (1991) Proteins: Struct., Funct., Genet. 10, 81-91] of the carboxy derivative of this hemoglobin (COHbY) revealed a new quaternary structure that partially resembles the recently described R2 structure [Silva, M. M., et al. (1992) J. Biol. Chem. 267, 17248-17256]. Very little is known about either the solution phase conformations of the liganded and deoxy forms of HbY or the molecular basis for the large quaternary enhancement effect (Doyle et al., 1992). In this study, near-IR absorption, Soret-enhanced Raman, and UV (229 nm) resonance Raman spectroscopies are used to probe the liganded and deoxy derivatives of HbY in solution. Nanosecond time-resolved near-IR absorption measurements are used to expose the relaxation properties of the photoproduct of COHbY. Time-resolved (Soret band) absorption is used to generate the geminate and solvent phase ligand rebinding curves for photodissociated COHbY. The spectroscopic results indicate that COHbY has an R-like conformation with respect to both the proximal heme pocket and the hinge region of the alpha 1 beta 2 interface. The deoxy derivative of HbY has spectroscopic features that are very similar to those observed for species assigned to the deoxy R or half-liganded R conformations of human adult hemoglobin (HbA). The 10 ns to 100 micros relaxation properties of the photoproduct of COHbY are distinctly different from those of HbA in that for HbY, little if any tertiary or quaternary relaxation is observed. The near-absence of relaxation in the HbY photoproduct explains the differences in the geminate and solvent phase CO recombination between HbA and HbY. The impact of the conformational and relaxation properties of HbY on the geminate rebinding process forms the basis of a model that accounts for the large quaternary enhancement effect reported for HbY (Doyle et al., 1992). In addition, the spectroscopic data and the X-ray crystallographic results explain the slow relaxation for HbY and the near-absence of cooperative ligand binding for this protein based on the behavior of the penultimate tyrosines.     

Use of an in situ disulfide cross-linking strategy to study the dynamic properties of the cytoplasmic end of transmembrane domain VI of the M3 muscarinic acetylcholine receptor

The ligand-induced activation of G protein-coupled receptors (GPCRs) is predicted to involve pronounced conformational changes on the intracellular surface or the receptor proteins. A reorientation of the cytoplasmic end of transmembrane domain VI (TM VI) is thought to play a key role in GPCR activation and productive receptor/G protein coupling. Disulfide cross-linking studies with solubilized, Cys-substituted mutant versions of bovine rhodopsin and the M3 muscarinic acetylcholine receptor suggested that the cytoplasmic end of TM VI is conformationally highly flexible, even in the absence of activating ligands (Farrens, D. L., et al. (1996) Science 274, 768-770; Zeng, F. Y., et al. (1999) J. Biol. Chem. 274, 16629-16640). To test the hypothesis that the promiscuous disulfide cross-linking pattern observed in these studies was caused by the use of solubilized receptor proteins endowed with increased conformational flexibility, we employed a recently developed in situ disulfide cross-linking strategy that allows the detection of disulfide bonds in Cys-substituted mutant M3 muscarinic receptors present in their native membrane environment. Specifically, we used membranes prepared from transfected COS-7 cells to analyze a series of double Cys mutant M3 receptors containing one Cys residue within the sequence K484(6.29) to S493(6.38) at the cytoplasmic end of TM VI and a second Cys residue at the cytoplasmic end of TM III (I169C(3.54)). This analysis revealed a disulfide cross-linking pattern that was strikingly more restricted than that observed previously with solubilized receptor proteins, both in the absence and in the presence of the muscarinic agonist, carbachol. Carbachol stimulated the formation of disulfide bonds in only two of the 10 analyzed mutant muscarinic receptors, I169C(3.54)/K484C(6.29) and I169C(3.54)/A488C(6.33), consistent with an agonist-induced rotation of the cytoplasmic end of TM VI. These findings underline the usefulness of analyzing the structural and dynamic properties of GPCRs in their native lipid environment.     

Predictive modelling of topology and loop variations in dimeric DNA quadruplex structures

We have used a combination of simulated annealing (SA), molecular dynamics (MD) and locally enhanced sampling (LES) methods in order to predict the favourable topologies and loop conformations of dimeric DNA quadruplexes with T2 or T3 loops. This follows on from our previous MD simulation studies on the influence of loop lengths on the topology of intramolecular quadruplex structures [P. Hazel et al. (2004) J. Am. Chem. Soc., 126, 16405–16415], which provided results consistent with biophysical data. The recent crystal structures of d(G4T3G4)2 and d(G4BrUT2G4) (P. Hazel et al. (2006) J. Am. Chem. Soc., in press) and the NMR-determined topology of d(TG4T2G4T)2 [A.T. Phan et al. (2004) J. Mol. Biol., 338, 93–102] have been used in the present study for comparison with simulation results. These together with MM-PBSA free-energy calculations indicate that lateral T3 loops are favoured over diagonal loops, in accordance with the experimental structures; however, distinct loop conformations have been predicted to be favoured compared to those found experimentally. Several lateral and diagonal loop conformations have been found to be similar in energy. The simulations suggest an explanation for the distinct patterns of observed dimer topology for sequences with T3 and T2 loops, which depend on the loop lengths, rather than only on G-quartet stability.