3D modeling of the activated states of constitutively active mutants of 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.

Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the beta(2) adrenergic receptor (beta(2)AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the beta(2)AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in beta(2)AR, namely, the breaking of the ionic lock between R131(3.50) at the intracellular end of TM3 (part of the DRY motif) and E268(6.30) on TM6, and the rotamer toggle switch on W286(6.48) on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.

Synthesis and spectroscopic characterization of photo-affinity peptide ligands to study rhodopsin-G protein interaction

G protein-coupled receptors (GPCRs) are involved in the control of virtually all aspects of our behavior and physiology. Activated receptors catalyze nucleotide exchange in heterotrimeric G proteins (composed of alpha.GDP, beta and gamma subunits) on the inner surface of the cell membrane. The GPCR rhodopsin and the G protein transducin (G(t)) are key proteins in the early steps of the visual cascade. The main receptor interaction sites on G(t) are the C-terminal tail of the G(t)alpha-subunit and the farnesylated C-terminal tail of the G(t)gamma-subunit. Synthetic peptides derived from these C-termini specifically bind and stabilize the active rhodopsin conformation (R*). Here we report the synthesis of R*-interacting peptides containing photo-reactive groups with a specific isotope pattern, which can facilitate detection of cross-linked products by mass spectrometry. In a preliminary set of experiments, we characterized such peptides derived from the farnesylated G(t)gamma C-terminus (G(t)gamma(60-71)far) in terms of their capability to bind R*. Here, we describe novel peptides with photo-affinity labels that bind R* with affinities similar to that of the native G(t)gamma(60-71)far peptide. Such peptides will enable an improved experimental strategy to probe rhodopsin-G(t) interaction and to map so far unknown interaction sites between both proteins.

Conformational movement of F251 contributes to the molecular mechanism of constitutive activation in the C5a receptor

The activation mechanism of G-protein-coupled receptors triggered upon binding of a ligand represents a very important 'conformational switch' in the biological array of signal transduction. However, the molecular and functional details for this activation switch remain unknown. Random saturation mutagenesis data on the complement factor 5a receptor has provided a large data set of mutants including several constitutively active mutants. In the present study, we employed computational modeling to rationalize the constitutive activity for two constitutively active mutants, NQ (I124N/L127Q) and F251A, and we then made predictions for a series of mutants that either promote or constrain constitutive activity. Biological testing of the site-directed mutants confirmed most of the predictions of the computational modeling. These results support a molecular mechanism of constitutive activity in complement factor 5a receptor mutants that is associated with conformational changes in a network of residues neighboring F251 as the focal point of origin.

Computational modeling and functional analysis of Herpes simplex virus type-1 thymidine kinase and Escherichia coli cytosine deaminase fusion protein

Herpes simplex virus type-1 thymidine kinase (HSV-1TK) and Escherichia coli cytosine deaminase (CD) fusion protein was designed using InsightII software. The structural rationality of the fusion proteins incorporating a series of flexible linker peptide was analyzed, and a suitable linker peptide was chosen for further investigated. The recombinant plasmid containing the coding regions of HSV-1TK and CD cDNA connected by this linker peptide coding sequence was generated and subsequently transfected into the human embryonic kidney 293 cells (HEK293). The Western blotting indicated that the recombinant fusion protein existed as a dimer with a molecular weight of approximately 90 kDa. The toxicity of the prodrug on the recombinant plasmid-transfected human lung cancer cell line NCIH460 was evaluated, which showed that TKglyCD-expressing cells conferred upon cells prodrug sensitivities equivalent to that observed for each enzyme independently. Most noteworthy, cytotoxicity could be enhanced by concurrently treating TKglyCD-expressing cells with prodrugs GCV and 5-FC. The results indicate that we have successfully constructed a HSV-1TKglyCD fusion gene which might have a potential application for cancer gene therapy.

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.