|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 39, 28779-28790, September 28, 2007
Mechanism of Activation of a G Protein-coupled Receptor, the Human Cholecystokinin-2 Receptor*From the INSERM, Unit 858, Institut Fédératif de Recherche, 31432 Toulouse, France and the Université de Toulouse 3, 31062 Toulouse, France
Received for publication, January 12, 2007 , and in revised form, June 21, 2007.
G protein-coupled receptors (GPCRs) represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts with the still limited data available concerning their activation mechanism. Here, we investigated the activation mechanism of the cholecystokinin-2 receptor (CCK2R). The three-dimensional structure of inactive CCK2R was homology-modeled on the basis of crystal coordinates of inactive rhodopsin. Starting from the inactive CCK2R modeled structure, active CCK2R (namely cholecystokinin-occupied CCK2R) was modeled by means of steered molecular dynamics in a lipid bilayer and by using available data from other GPCRs, including rhodopsin. By comparing the modeled structures of the inactive and active CCK2R, we identified changes in the relative position of helices and networks of interacting residues, which were expected to stabilize either the active or inactive states of CCK2R. Using targeted molecular dynamics simulations capable of converting CCK2R from the inactive to the active state, we delineated structural changes at the atomic level. The activation mechanism involved significant movements of helices VI and V, a slight movement of helices IV and VII, and changes in the position of critical residues within or near the binding site. The mutation of key amino acids yielded inactive or constitutively active CCK2R mutants, supporting this proposed mechanism. Such progress in the refinement of the CCK2R binding site structure and in knowledge of CCK2R activation mechanisms will enable target-based optimization of nonpeptide ligands.
The superfamily of G protein-coupled receptors (GPCRs)5 represents a large and diverse family of seven transmembrane helix proteins. Because they are involved in all physiological processes, GPCRs are targeted by 50% of marketed drugs. However, these drugs address only a small fraction of the GPCR repertoire. GPCRs represent a major focus in functional genomics programs and drug development research, but their important potential as drug targets contrasts strongly with the very limited data available regarding their activation mechanism. The activation of GPCRs is generally explained by the extended ternary complex model whereby the population of receptors exists as an equilibrium between inactive and active states (1). According to this model, the active conformation of the receptor may be stabilized or induced by agonists. The extended ternary complex model for GPCRs is based on and sustained by the discovery of engineered constitutively active mutants, namely mutants that produce ligand-independent signaling (1). Constitutively active GPCR mutants naturally exist in humans and can cause diseases (2). Interestingly, constitutive GPCR signaling can be inhibited or attenuated by a class of therapeutic agents named inverse agonists (2, 3).
The cholecystokinin (CCK)-2 receptor (CCK2R), previously named the CCKB/gastrin receptor, belongs to family I of GPCRs, which includes rhodopsin (4-6). It binds both cholecystokinin and gastrin, two structurally related neuropeptides with similar high affinities. CCK2R is expressed in the central nervous system and in the gut, where it represents the predominant CCK receptor subtype. CCK2R mediates a wide spectrum of CCK- and gastrin-induced biological effects, including anxiety, pain perception, gastric acid secretion, growth, and differentiation of the gastric mucosa (4, 7). Converging studies have documented that activation of wild-type CCK2R and/or expression of a constitutively active variant may contribute to human diseases (8, 9). These new findings have generated considerable interest in the identification of antagonists of CCK2R. To date, a large panel of chemically distinct CCK2R antagonists have been discovered and used to assess the functions mediated by CCK2R in animals and humans (10). Furthermore, several of these compounds have reached clinical evaluation stages for indications such as anxiety and panic disorders, sleep disorders, drug dependence, pain, gastroesophagus reflux, and gastric secretion disorders (10). However, some reference molecules believed to be pure antagonists turned out to be endowed with some agonist activity in the stomach and pancreas as well as cells expressing CCK1R or CCK2R (11-13). On the other hand, reports on the expression of constitutively active receptors in digestive cancers and on potential therapeutic relevance for neuroendocrine tumors highlighted the need to identify inverse agonists able to reverse basal constitutive activity of CCK2R (14). Target structure-based design of new ligands requires identification of the receptor binding site and determination of the activation mechanism. Toward this aim, our group has progressively delineated the binding site of CCK2R using site-directed mutagenesis combined with molecular modeling approaches (15-17). We identified several aromatic residues belonging to a network presumably critical for CCK2R activation (18). In this last study, differential motion of amino acid side chains within the binding pocket was observed upon in silico binding of different nonpeptide ligands that present distinct pharmacological features. Nevertheless, due to the intrinsic limitations of the modeling procedure used, no global change in CCK2R conformation was observed (18). The aim of the present study was to increase our knowledge of the molecular mechanisms whereby CCK2R is activated. For this purpose, three-dimensional models of inactive and active CCK2R were constructed on the basis of available structural data. These modeled structures were studied using unrestrained molecular dynamics simulations (MDs) in an explicit water-lipid environment. By comparing averaged modeled structures, we identified changes in the relative positions of helices and networks of interacting residues expected to stabilize either active or inactive states of CCK2R. The activation mechanism of the receptor was further investigated using targeted molecular dynamics simulations and site-directed mutagenesis.
Molecular Modeling Modeling of the "Inactive" CCK2R Structure (CCK2R°)—The structure of bovine rhodopsin (Protein Data Bank code 1U19) crystallized in the dark was used as structural template for the homology modeling of the CCK2R receptor (19). Alignment of human CCK2R sequence and bovine rhodopsin sequences was carried out using the ClustalW multiple alignment program (20). Loops connecting transmembrane helices were modeled using the ab initio modeling program MODLOOP (21). The third intracellular loop was not included in the modeled structure. The crude homology model was refined using AMBER8 (available on the World Wide Web) by means of 1000 steps of steepest descent energy minimization followed by 2000 steps of conjugate gradient energy minimization while all C carbons were harmonically restrained to their initial positions. This procedure allowed readjustment of covalent bonds and van der Waals contacts without changing the overall conformation of the protein. The model was then refined by molecular dynamics simulations in an explicit lipid bilayer. It was neutralized by the addition of 17 chloride ions and placed in a rectangular box ( 78 x 80 x 78 Å in size) containing a lipid bilayer (153 molecules of palmitoyloleoylphosphatidylcholine and 8708 molecules of water in addition to the receptor structure), resulting in a final density of 1.0 g/cm. MD simulations at 300 K and constant pressure (1 atm) with anisotropic scaling were then run for 1ns using the SANDER module in AMBER 8. The cut-off distance for the nonbonded interactions was 9 Å, and periodic boundary conditions were applied. Electrostatic interactions were represented using the smooth particle mesh Ewald method with a grid spacing of 1 Å. The coupling constants for the temperature and pressure baths were 1.0 and 0.2 ps, respectively. SHAKE was applied to all bonds involving hydrogen, and an integration step of 2 fs was used throughout. The nonbonded pair list was updated every 10 steps. The simulation protocol involved a series of progressive energy minimizations followed by a 20-ps heating phase and a 250-ps equilibration period, during which secondary structure restraints were applied to the transmembrane helices. These restraints were removed during the following 750 ps of data collection. System coordinates were saved every 2 ps for analysis. Modeling of the "Active" CCK2R Structure (CCK2R*)—The CCK2R° homology model obtained from the rhodopsin crystal was not directly suitable for docking purposes, because the ligand-binding pocket of inactive rhodopsin is remarkably narrow (19). We therefore refined the rhodopsin-based CCK2R° model in order to obtain a model of active CCK2R. Since it has been proposed that the conserved tryptophan belonging to the CWPFF motif in helix VI changes its orientation with respect to the membrane plane during activation, being perpendicular in the inactive form and parallel in the active form, the side chain of Trp-6.48 was manually rotated to adopt the latter conformation (22, 23). We then docked CCK4 into this crude model, taking into account all amino acids and key interactions identified in the course of binding site studies (15, 24, 25) and ran steered molecular dynamics simulations. In fact, steered molecular dynamics were first performed in the presence of CCK4, because contacts were only available for the last four residues of CCK, which correspond to the ligand binding region. Then the N-terminal extension of CCK was added. The simulation was performed in vacuo with a small spherical shell of 450 TIP3P water molecules centered around His-207. The radius of this water shell was 20 Å, ensuring the solvation of the extracellular portion of the helices and the extracellular loops. These water molecules were first energy-minimized, and then waters and protein residues were allowed to relax. In each case, 1000 steps of steepest descent were followed by 2000 steps of conjugate gradient energy minimization. The final coordinate sets were used as input for the subsequent steered MD simulations. SHAKE was used for all bonds, and the integration time step was 2 fs. The cut-off was 12 Å, and a distance-dependent dielectric function was used. The list of nonbonded pairs was updated every 25 steps, and coordinates were saved every 2 ps for further analysis.
The second run of 500-ps steered molecular dynamics simulations consisted of applying distance constraints by means of a harmonic potential as indicated in Table 1. The C
Once the desired distances were attained, CCK4 was substituted by CCK in the binding site. The final model was then refined by 1 ns of unrestrained MDs in an explicit lipid bilayer using the protocol described above for the CCK2R°. Finally, CCK was removed from the complex, and the free active conformation of CCK2R (CCK2R*) was also simulated by 1 ns of MD in a lipid environment. Although CCK2R* is energetically less stable than CCK2R°, the CCK2R* model did not relax back to CCK2R° during MD, because this back conversion would require a microsecond time scale and could not be observed with MD simulations in the nanoseconds time scale.
In Silico Conversion of CCK2R° to CCK2R* Using Targeted Molecular Dynamics Simulations
where E is the biasing potential energy in kcal/mol during the TMD simulation, kr is the force constant, N is the number of atoms, r.m.s.d. is the root mean square deviation from the target structure, and t.r.m.s.d. is the desired root mean square deviation value. The initial value of the desired root mean square deviation was set as the root mean square deviation between the initial and target structures (CCK2R° and CCK2R*, respectively), and during the simulation, the desired root mean square deviation value was gradually decreased to 0. Only the heavy atoms of the protein were considered in the r.m.s. definition. A force constant of 0.50 kcal/mol x Å over 0.5 ns appeared sufficient to find a low energy path leading from the simulated structure to the target structure.
Site-directed Mutagenesis of CCK2R and Transfection of COS-7 Cells
Inositol Phosphate Production Determination
Receptor Binding Assay
Modeled Structure of Active CCK2R (CCK2R*)—Because CCK2R is a member of the family I of rhodopsin-like GPCRs, the three-dimensional model of the inactive state of the CCK2R (CCK2R°) was constructed by homology modeling using coordinates from the inactive rhodopsin crystal isolated in the dark (19). To obtain CCK2R*, the agonist·CCK2R complex was modeled (Fig. 1). The procedure of modeling, which is described in detail under "Experimental Procedures," included the following main steps: rotation of conserved Trp-6.48 side chain (22); docking of CCK4 followed by steered molecular dynamics simulation in vacuum; steered molecular dynamics simulation based on constrained distances indicated in Table 1 and derived from published data on GPCR activation; and refinement in a lipid bilayer (see Table 1 for references).
As illustrated in Fig. 2, after the equilibration period, progression of the r.m.s. deviations of the coordinates of the C
Comparison of the Modeled Structures of CCK2R° and CCK2R*—In order to investigate the activation mechanism of CCK2R, we first compared averaged structures of refined CCK2R° and CCK2R* models. Actually, we assumed that CCK2R naturally exists in an equilibrium between receptors in the inactive (CCK2R°) and active state (CCK2R*) and that specific intramolecular networks stabilize each conformation. In support of this assumption, we previously showed that human wild-type CCK2R presents some degree of constitutive activity (18).
In silico CCK2R° to CCK2R* conversion brought together the upper parts of helices VI and III (Table 2) and caused clockwise rotation of helix VI (observed from extracellular space). Furthermore, in CCK2R*, the cytoplasmic ends of helices VI and VII were more distant from the end of helix III than in CCK2R°, which can be inferred from two distances, one between the Asn-6.55 C
Regarding helix IV, although its movement was small, its extracellular end was brought closer to that of helix VI as a result of its rotation and tilting. Indeed, in CCK2R*, Tyr-4.60 formed a hydrogen bond with Asn-6.55, whereas in CCK2R°, Tyr-4.60 established contacts with Met-4.57 and Leu-5.42 (Fig. 4C). Concerning helix V, comparison of the CCK2R* and CCK2R° structures suggested that this helix underwent a concerted movement to allow helix VI to rotate clockwise during activation. In fact, in CCK2R*, Trp-6.48 and Tyr-6.52 have moved to positions previously occupied by Phe-5.47 and Leu-5.43, respectively, in CCK2R° (Fig. 4D). Finally, helix VII was positioned 2 Å higher in CCK2R* relative to CCK2R°. This motion disrupted a network of interactions present in CCK2R° between Tyr-7.43, His-7.39, and Met-3.32. All of these observations suggested that interhelical contacts described above between helices III, IV, V, VI, and VII may be critical for CCK2R activation. Delineation of the CCK2R° to CCK2R* Conversion by Targeted Molecular Dynamics Simulations—The comparison between modeled structures of CCK2R° and CCK2R* suggested that observed changes in amino acid and helix positions are part of the activation mechanism of CCK2R. However, these data could not answer the questions of how, at a molecular level, the changes involving movements of helices and amino acids actually proceed, in which order they occur, and how they are coordinated. To answer these questions, TMD, an in silico method that can simulate intricate conformational changes in proteins as the pathway between two conformations, was applied to CCK2R (27).
In these in silico experiments, we observed that a force constant of 0.50 kcal mol-1 Å-2 over 500 ps was sufficient to achieve the CCK2R° to CCK2R* conversion, as can be inferred from the final value of the r.m.s. deviation of the backbone atoms. The TMD showed that the rotation of the side chain of Trp-6.48 is a consequence of a completed coordinated tilting and rotation of helix VI and movement of helix V. As a criterion for the tilting of helix VI, two distances, the first one between the C
Experimental Evidence for the Importance of Aromatic/Hydrophobic Residues in the CCK2R Activation—To verify the role of the aforementioned residues in the activation mechanism, they were mutated, and the basal activities of the resulting mutants were determined by measuring inositol phosphate production in transfected COS-7 cells in the absence of agonist stimulation (basal activity) and the presence of CCK. By doing so, we expected that mutation of amino acids that stabilize the receptor in the inactive state should generate constitutively active mutants, whereas mutation of amino acids involved in intramolecular interactions stabilizing the active form should produce mutants devoid of basal activity. As illustrated in Fig. 8, exchange of Val-3.36, Met-4.57, Leu-5.42, Leu-5.43, Phe-5.47, and Trp-6.48 with an Ala yielded constitutively active mutants, indicating that these amino acids are required for stabilization of CCK2R°. Strikingly, mutants V3.36A, L5.42A, L5.43A, and F5.47A displayed robust constitutive activities, since inositol phosphate production in cells transfected with these mutants represented 3-4-fold that of cells transfected with wild-type CCK2R. Mutant W6.48A exhibited a 2.5-fold enhanced activity relative to wild-type CCK2R. Mutant M3.32A and Y7.43A showed a slight or a marked constitutive activity, respectively. For comparison, mutation E3.49A in the conserved (E/D)RY motif, previously found to cause oncogenic activity of the CCK2R, showed a 2-fold increased activity relative to the wild-type CCK2R (26). In contrast, mutation of Tyr-4.60, Phe-6.44, Asn-6.55, and Tyr-6.52 produced CCK2R mutants that exhibited no detectable or significantly lower basal activity than that of wild-type CCK2R, thus ascribing a role in the stability of the active state of the receptor to these amino acids. On the other hand, mutation of Ser-7.45 and His-7.39 did not change basal activity of CCK2R.
The consistency between biological data and structural changes during the CCK2R° to CCK2R* transition was examined. Biological data showing constitutive activity of mutants W6.48A, V3.36A, and M3.32A support the view deduced from analysis of the modeled CCK2R° structure that interactions between Trp-6.48 and Val-3.36/Met-3.32 in CCK2R° stabilize this state. The biological data also suggest that Leu-5.42 and Leu-5.43, which do not interact with Trp-6.48 in CCK2R°, play a role in the maintenance of the inactive state. To explain this result, we suggest that Leu-5.43 and Phe-5.47 maintain CCK2R in its inactive state, because the motion of these two residues is required for, and must precede, rotation of Trp-6.48 and Tyr-6.52. It is also worthy to note that mutation of Ser-7.45, which appeared to be hydrogen-bound with Trp-6.48, did not produce a constitutive CCK2R, indicating that disruption of a hydrogen bond was not sufficient for CCK2R° to CCK2R* conversion. This result fits with the view that an intermediate state was generated by hydrogen bond disruption.
The finding that Y4.60F and N6.55A mutants lack basal activity fits with the view that the hydrogen bonding interaction between these two residues is important for the stability of CCK2R* conformation. The hydrogen bond between Tyr-4.60 and Asn-6.55 in CCK2R* appears possible because the Tyr-4.60 side chain has moved to occupy the position of Leu-5.42 in the CCK2R*, pushing it outside of the binding pocket. As a consequence, interaction between Met-4.57 and Leu-5.42 is disrupted during CCK2R° to CCK2R* conversion. The roles of Met-4.57 and Leu-5.42 in the maintenance of the CCK2R° conformation may be explained in this way. At last, pharmacological experiments were carried out with all mutants to measure their expression level, their binding properties, and their ability to respond to CCK stimulation. A summary of results from these experiments is presented in Table 3.
In this work, using in parallel state-of-the-art molecular modeling and site-directed mutagenesis approaches, the molecular mechanism of CCK2R activation was investigated with respect to the activation of phospholipase C by measuring inositol phosphate production. A number of amino acids that participate in the activation of CCK2R by stabilizing either active or inactive conformations of the receptor were identified, and a global mechanism of activation involving several key amino acids within or near the binding site as well as movements of transmembrane helices IV, V, VI, and VII was deduced. We considered that CCK2R, like other G protein-coupled receptors, naturally exists in an equilibrium between receptors in the inactive conformation and receptors in active conformation, the latter being responsible for constitutive activity of the receptor. Previously, we demonstrated that wild-type human CCK2R displays such a constitutive activity, since it enhances production of inositol phosphates in COS-7 cells to a level that is directly a function of the level of receptor expression at the cell surface (18). As demonstrated with other GPCRs, we presumed that adequate mutation of amino acids critical for the stabilization of the inactive state should generate constitutively active mutants, whereas the mutation of amino acids critical for the stabilization of the active state should yield a receptor with zero or decreased basal activity. According to this model, amino acids that appear to be involved in stabilization of CCK2R° are Met-3.32, Val-3.36, Met-4.57, Leu-5.42, Leu-5.43, Phe-5.47, and Trp-6.48, whereas those that stabilize CCK2R* are Tyr-4.60, Phe-6.44, Tyr-6.52, and Asn-6.55. Importantly, all analyzed mutants retained a good ability to respond to the full agonist CCK. On the other hand, constitutively active mutants displayed an activity that remained moderate relative to the maximal CCK-induced response.
One of the most significant structural changes occurring during CCK2R activation concerns helix VI, which undergoes a clockwise rotation where the extracellular end is tilted toward helix III. Targeted molecular dynamics revealed that helix VI movement is coordinated with rotation of helix V and upper movement of helix VII. As a consequence of these movements, interhelical contacts change from the inactive to active state, the most remarkable being the molecular network around the conserved Trp-6.48. Since the above mechanism of CCK2R activation was deduced from the in silico conversion of CCK2R° to empty CCK2R* and the role of the agonist is to stabilize CCK2R*, the consistency of this activation mechanism with the mode of CCK binding to CCK2R* must be considered. In helix VI, Arg-6.58 and Asn-6.55 participate in CCK binding. In empty CCK2R*, the Arg-6.58 guanidium is associated with water molecules, whereas Asn-6.55 is in tight interaction with the Tyr-4.60 hydroxyl and thus contributes to CCK2R* stabilization. In the presence of CCK, the Arg-6.58 guanidium moves slightly toward the binding cavity to interact with the side chain of CCK-aspartic acid. Interestingly, following CCK binding, the Asn-6.55-Tyr-4.60 bond dissociates, leaving the side chains of these amino acids free to accommodate the CCK amide through hydrogen bonds. As for CCK-Trp, its location between helix III and VII through interactions with His-7.39 and Tyr-7.43 contributes to the stabilization of helix VII in the upper position as found in empty CCK2R*. In the absence of CCK, in CCK2R*, helix VII position is locked by interactions between His-7.39 and Tyr-7.43 with Met-3.32. Furthermore, in CCK2R* an aromatic/hydrophobic pocket is enlarged due to the rotation of amino acids surrounding Trp-6.48. CCK-Phe perfectly accommodates, in this binding pocket formed by Val-3.36, Met4.57 and Tyr4.60, Leu5.46, and Val6.51 and Tyr6.52 and thus contributes to the stabilization of CCK2R*. Several amino acids critical for CCK2R activation, such as Trp-6.48, Leu-5.42, Leu-5.43, and Phe-5.47, are not in contact with CCK, since their side chains are directed outside the binding cavity in CCK2R. All together, these data provide an original picture of how CCK2R can adopt the active conformation on its own, with an energy barrier between active and inactive conformations that is probably small. Furthermore, they show how the peptide agonist, CCK, binds to CCK2R* to stabilize it through specific interactions, some of which concern anchoring points, whereas others serve to maintain helices at critical positions in CCK2R*.
The above global structural changes observed during CCK2R activation are in line with those described or suggested with other GPCRs using biochemical and biophysical techniques, but they appear quantitatively smaller (28-30). Concerning movements of helix VI, pioneering studies using electron paramagnetic resonance analysis of spin-labeled rhodopsin have revealed the proximity of helices III and VI and rigid body movement of helix VI upon light activation (31). This was confirmed by Cys cross-linking and metal-ion binding of these helices, which prevent activation of the rhodopsin-coupled G protein, transducin (31, 32). It is worthy of mention that the strategy of engineered metal-ion binding sites between transmembrane helices applied to a number of peptide receptors confirmed the proximity of the upper part of helices III and VI and their relative motion during receptor activation (33). Interestingly, the maximal efficacy of metals was obtained when the metal ion site was formed by amino acids in positions III:08 and VI:16, which correspond to Met-3.32 and Val-6.51 of CCK2R, respectively. A distance of about 8 Å between the C
Concerning helix V, which undergoes a concerted rotation with helix VII during CCK2R activation, identification of Phe-5.47 as a key residue is in line with the concept that a conserved Phe/Tyr residue in this helix (located at -3 residues of a proline and corresponding to Phe-5.47 in CCK2R) plays a critical role in the activation switch. In ghrelin and With respect to helix VII, its role as an activation or coupling switch has been illustrated. For instance, in the C5a receptor, motion of helix VII has been proposed to explain the activation switch observed with mutants (35). Movement of the cytoplasmic interface of this helix during activation is supported by data showing that the NPXXY motif is accessible to an antibody in metarhodopsin, whereas it is not in rhodopsin maintained in the dark (36). Furthermore, it has been shown that upon formation of metarhodopsin, the salt bridge between Glu-113 in helix III and the protonated Schiff base formed by retinal and Lys-296 in helix VII is disrupted, which requires movement of helix VII relative to helix III (37). In CCK2R, we have noticed an upper movement of helix VII of about 2 Å. Such a movement, if applied to rhodopsin helix VII, would cause disruption of the Glu113-Lys296 salt bridge and would probably establish contact between Lys296 and Glu181 (37). TMD with CCK2R showed that the upper movement of helix VII starts once the Trp-6.48 side chain has begun to rotate and the hydrogen bonding interaction with Ser-7.45 is lost. Disruption of the hydrogen bond may be an important event in the activation process of CCK2R, but the fact that mutation of Ser-7.45 did not generate constitutive activity also indicated that disruption of this hydrogen bond is not sufficient for CCK2R activation. Asn-7.45 in histamine H1 receptor, the cognate residue of Ser-7.45, was reported to play a role in receptor activation (38). Indeed, it was demonstrated that in the H1 receptor, Asn-7.45 constrained Trp-6.48 in the inactive g+ conformation through a hydrogen bonding interaction and that, unlike in CCK2R, disruption of this hydrogen bond by leucine substitution leads to a constitutively active receptor (38). The critical role of helix IV in CCK2R activation is remarkable. This helix does not undergo significant movement during CCK2R°-CCK2R* conversion, but its distance with respect to helix VI seems critical in CCK2R*, as supported by the dramatic decrease of the distance between the two helices in CCK2R* relative to the CCK2R° and the role played by the CCK-amide in the maintenance of this proximity. Helix IV may therefore be considered a part of the rigid body of the receptor and contributes to the correct positioning of helix VI. The importance of helix IV as an anchoring point for peptidic ligand binding has been documented for other peptide receptors. In the CXCR4 chemokine receptor, urotensin II receptor, vasopressin V1a receptor, and melanocortin-4 receptor, the residue equivalent to Tyr-4.60 in CCK2R is critical for agonist binding (39-42). However, with the current work on CCK2R, this is the first example of a key role being ascribed to helix IV in the control of the equilibrium between inactive and active conformations.
Taken together, structural changes observed during CCK2R activation are in line with those described or suggested with other GPCRs using biochemical and biophysical techniques, but they appear quantitatively smaller (28-30). Instead, changes that accompany CCK2R activation are quantitatively more comparable with those deduced from the recent crystallization of a photoactivated, deprotonated intermediate of rhodopsin (43). Furthermore, our data point out that the role played by some conserved amino acids in the activation process is different in CCK2R relative to other members of the GPCR family A. Hence, during activation, CCK2R undergoes conformational changes that are shared by members of GPCR class A but also appears to present divergent features, which may account for its ligand and signaling selectivity. In the future, studies on GPCR mechanism of activation will have to take into account the entire structure of the receptor, including the third intracellular loop, the structure of which is dependent on the relative position of helices V and VI as well as on binding of the It is worthy of note that several amino acids identified as critical in the activation process of CCK2R have been previously mutated. For instance, the exchange of Met-3.32, Val-3.36, and Leu-5.43 was shown to increase basal activity of CCK2R in COS-7 cells (24, 44). Interestingly, in the mentioned studies, partial agonist activity of PD-135,158, a molecule initially believed to be an antagonist, was increased on M4.57A and W6.48A mutants, whereas it disappeared on F5.47A and Y6.52A mutants. Such mutants were previously defined as "a subgroup of activating mutations which facilitates receptor isomerization to the active state and in parallel perturbs ligand receptor interactions" (44). Our data provide a significant advance, since they show how CCK2R can adopt the active conformation on its own and how the peptide agonist, CCK, binds to CCK2R* to stabilize it. Obviously, additional works using biochemical and biophysical approaches are highly desirable to further verify the findings of the present study. Last, this study with CCK2R illustrates how TMD, in contrast to ordinary molecular dynamics simulation, is a powerful molecular modeling technique able to study conformational conversion between inactive and active states of a complex membrane protein such as a GPCR. Previously, TMD has been successfully applied to study conformational transitions of KcsA channel, chymotrypsin, and other proteins (27, 45, 46). However, to our knowledge, analysis of CCK2R activation represents the first successful application of TMD to a GPCR. Based on the findings from this study, TMD appears appropriate to study, at the atomic level, the structural determinants involved in the activation of a GPCR. Interestingly, the intermediary receptor structures that can be identified and extracted from the TMD trajectories can be used as targets for ligand docking studies. Besides providing information of basic scientific interest, it is expected that insight into the activation and inactivation mechanism of GPCRs will be exploited in structure-based drug discovery efforts.
* This work was supported in part by Association pour la Recherche contre le Cancer Grants ARC 3282 and 3756. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 Supported by fellowships from the "Fondation pour la Recherche Medicale."
2 These two authors contributed equally to this work.
3 "Chargée de Recherches" from the Fonds National de la Recherche Scientifique, Belgium. 4 To whom correspondence should be addressed: IFR 31, Institut Louis Bugnard, BP 84225, Unité 858, 31432 Toulouse cedex 4, France. E-mail: fourmyd{at}toulouse.inserm.fr.
5 The abbreviations used are: GPCR, G protein-coupled receptor; CCK, cholecystokinin; CCK2R, cholecystokinin-2 receptor; MD, molecular dynamics; DMEM, Dulbecco's modified Eagle's medium; Ins-P, inositol phosphate; TMD, targeted molecular dynamics; r.m.s., root mean square.
We sincerely thank the Centre Informatique National de l'Enseignement Supérieur (Montpellier, France) for software and hardware support.
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||