(Received for publication, July 31, 1996, and in revised form, December 3, 1996)
From the Laboratory of Pharmacology, Istituto Superiore di Sanità, Viale Regina Elena, 299, 00161 Rome, Italy
We studied the biochemical properties of a genetically engineered neurokinin-1 receptor (NK1R) in which two residues lying on the extracellular edge of the fourth transmembrane domain were replaced by equivalently located elements of the neurokinin-2 receptor (G166C, Y167F NK1R mutant). The mutation produced two effects. The first is enhancement of the apparent binding affinity for heterologous tachykinins (substance K and neurokinin B) and for N- or C-terminal modified analogues of substance P, but not for substance P itself, its full-length analogues, and several peptide and nonpeptide antagonists. Only two antagonists, as exceptions, were found to exhibit a diminished affinity for the mutant receptor. The second effect is a shift in NK1R preference for distinct G protein-mediated signaling pathways. NK1R-mediated phosphoinositide hydrolysis was enhanced both in transiently and permanently transfected cells, while stimulation of cAMP accumulation did not change in transient expression experiments and was reduced in permanently expressing cells.
The effect of the mutation on ligand affinity was not related to any obvious structural commonality, nor to the selectivity for different neurokinin receptors or the agonistic/antagonistic nature of the ligand. However, all ligands responding to the mutation appear to share the ability to induce phosphoinositide signaling more efficiently than cAMP responses when binding to NK1R. We suggest that the mutation shifts the internal equilibria of different functional forms of NK1R. A theoretical analysis according to a multistate allosteric model suggests that the link between binding and biological changes can result from altered stability constants of substates in the conformational space of the receptor.
Tachykinins, as several other families of neuropeptides, display a bipartite distribution of topochemical information on their sequence. The C-terminal half of the molecule is a consensus motif tightly conserved among all known hormones from mammalian and nonmammalian species (1, 2) and, conceivably, constitutes a default conditional element to establish binding affinity for any kind of tachykinin receptor subtype (3, 4). The N-terminal part is instead variable even among the three mammalian tachykinin types, and it is believed to contribute interactions establishing selectivity for specific receptor subtypes (3, 4). The interesting question is whether a similar partition in the molecule of the receptor corresponds to such a sharp division between variant and invariant elements in the molecule of the peptide and, if so, to what relative extent do selective and nonselective interactions contribute to the final binding affinity of each peptide hormone and receptor system.
We focused on two amino acids marking the junction between the fourth putative transmembrane domain (TMD)1 and the second extracellular loop of tachykinin receptors. Replacement of these two residues in the substance P receptor (neurokinin-1 receptor (NK1R)) with those located at equivalent positions in the substance K receptor sequence produced changes in receptor behavior suggesting a specific role of the C-terminal consensus sequence of the peptide in determining sensitivity to the two-residue transmutation (5). It was thus proposed that the two residues introduced by the mutation may induce in the receptor an improved docking configuration for the consensus sequence of tachykinins (5).
Diverse lines of evidence support the notion that residues located between the C-terminal portion of TMD-4 and the N-terminal half of the second extracellular loop play an important role in recognition and binding of tachykinin ligands. First, the two amino acids targeted by the mutation are part of a stretch of five residues that was found to be shared by NK1R and a monoclonal antibody specific for the C-terminal sequence of substance P (5). In a comparative study where three-dimensional models of both macromolecules were examined (6), this sequence appears to be involved in forming part of the surface of a putative pocket displaying strikingly similar charge distributions and may thus be involved in binding the C-terminal sequence of substance P. Second, site-directed mutagenesis of residues in human NK1R closely flanking the two amino acids modified here, such as Gln-165 (7, 8) and Ser-169 (7), results in large changes in the binding affinity for both peptide and nonpeptide ligands. Finally, affinity labeling studies with analogues of substance P bearing photoactivable moieties in the consensus sequence demonstrated alkylation of receptor residues lying in the N-terminal portion of the second extracellular loop and close to the site mutagenized in this study (9-11).
To gain additional understanding of the role of this site of NK1R, we compared the binding affinities of a variety of peptide and nonpeptide ligands. We show evidence suggesting that the main effect of this mutation is an allosteric effect on the conformational equilibrium of the receptor rather than a change in the local docking interactions between the targeted residues and a specific area of the tachykinin sequence. This conformational shift produces two main effects: (a) a change in the binding energy for a chemically heterogeneous group of tachykinin ligands, all of which share the ability to induce polyphosphoinositide (PI) turnover stimulation more efficiently than cAMP accumulation upon activation of NK1R; and (b) an increased ability of the agonist-bound receptor to induce Gq-mediated responses, with no apparent change or even a reduced efficiency in triggering the Gs-mediated pathway.
Monoiodinated 125I-Bolton-Hunter reagent-labeled substance P (2000 Ci/mmol) was purchased from Amersham Corp. All peptides were purchased from Bachem or American Peptide Co. FK224 (12) and FK888 (13) were provided by Dr. Matsuo (Fujisawa Pharmaceuticals Co.). SR 140333 (14) was supplied by Dr. Edmonds-Alt (Senofi Recherche). Men 10930 (15) was provided by Dr. Manzini (A. Menarini Pharmaceuticals). CP 96345 (16), CP 99994 (17), RP 67580 (18), and CGP 49823 (19) were a gift of Dr. Schwartz (University of Copenhagen). All other reagents were from Sigma.
Mutant and Wild-type NK1R GenesThe constructs used in this study were rat NK1R (pCDM8-SPR, kindly donated by Dr. S. Nakanishi) and its G166C,Y167F mutant, prepared by the oligonucleotide-directed M13 method as described previously (5).
Cell CulturePermanently transfected CHO cells were grown in a mixture of Dulbecco's modified Eagle's medium and Ham's F-12 medium (1:1) supplemented with 10% fetal calf serum, 100 units/ml penicillin G, 100 µg/ml streptomycin sulfate, and 100 µg/ml G418 (Life Technologies, Inc.) at 37 °C in a humidified atmosphere of 5% CO2. COS-1 cells and C6 glioma cells were grown under identical conditions, except for the growth medium (only Dulbecco's modified Eagle's medium) and the omission of G418.
Transfection of CellsWild-type and mutant receptors were transfected in CHO cells using Lipofectin (Life Technologies, Inc.), and stably expressing clones were isolated following selection with 400 µg/ml G418 (Geneticin). For transient expression, COS-1 cells were seeded in 24-well plates and transiently transfected with wild-type or mutant receptor cDNA using a DEAE-dextran/chloroquine procedure as described (20). Gradual levels of receptor expression were obtained by transfecting varying amounts of receptor DNA but maintaining the total mass of transfected DNA (2 µg/well) constant with the addition of empty vector.
Phosphoinositide Turnover and cAMP AssayTo compare the
biological activities of mutant NK1R and the wild-type
receptor, we measured phosphoinositide hydrolysis and cAMP
accumulation. Identical procedures were used for transient expression
experiments or permanently expressing cells, except that in the latter
case, the two second messengers were assayed in the same cell extract.
Cells were seeded in 24-well plates and allowed to grow to confluence
prior to the addition of 2 µCi/ml myo-[3H]inositol (80-90 Ci/mmol; DuPont). For
transient tranfections in COS cells, the labeled inositol was added
48 h following transfection. After 24 h of incorporation, the
medium was removed, and cells were incubated in a buffer containing 137 mM NaCl, 2.7 mM KCl, 1.5 mM
KH2PO4, 20 mM NaHepes, 2 mM CaCl2, 1.2 mM MgSO4,
1 mM EGTA, 11.1 mM D-glucose, and
0.1% bovine serum albumin, pH 7.4. This was supplemented with 10 mM LiCl, 100 µM Ro 20-1724, and various concentrations of substance P or other peptides as required by the
experiments. Reactions were conducted for 30 min at 37 °C and
arrested by the removal of the supernatant and the addition of 0.5 ml
of ice-cold 0.1 N HCl to each well. The plates were placed
on ice, and 50 µl of the HCl extract were removed for the determination of cAMP concentration. 1 ml of ice-cold methanol was then
added to each well, taking care not to disturb the monolayer, and the
resulting mixture was carefully collected for the separation of
radiolabeled inositol phosphates. Monolayers were allowed to dry under
a hood and then dissolved into 0.5 ml of a solution prepared by mixing
0.8 M NaOH and 5% SDS (1:2). 200 µl of this extract were
neutralized with 30 µl of glacial acetic acid, mixed with
scintillation mixture (Ultima Flo, Packard Instrument Co.), and counted
in a -counter to determine the levels of the total acid-insoluble
inositol incorporated into the cell membranes. The remaining was used
for protein determinations.
The concentration of radioactive inositol phosphates hydrolyzed was quantified after application of the acid methanolic extract to anion-exchange columns (AG 1-X8, Bio-Rad) and elution into a single fraction as described previously (21). Cyclic AMP levels were determined by radioimmunoassay following acetylation of the sample (22).
Radioreceptor Binding AssaysBinding assays in intact cells were performed as described previously (5) with minor modifications. Cells were seeded in opaque culture plates (Canberra-Packard) and, at confluence, incubated in a reaction buffer with a composition identical to that used for the study of biological activity, supplemented with various concentrations of test peptide and monoiodinated 125I-Bolton-Hunter reagent-labeled substance P (20,000 cpm). Incubations lasted 3 h at 4 °C and were terminated by aspirating the incubation medium and washing the monolayer twice with ice-cold binding buffer. Plates were allowed to drain overnight onto filter paper; then 250 µl of MicroScint (Packard Instrument Co.) were added to each well, and the plates were counted in a Top Count (Canberra-Packard). Binding isotherms of 12 logarithmically spaced concentrations performed in duplicate were routinely used to determine binding parameters.
Data Presentation and AnalysisPhosphoinositide turnover is expressed as the hydrolyzed fraction of total incorporated inositol, calculated as follows: dpm(IP1)/(dpm(Ins) + dpm(IP1)), where dpm(IP1) is the radioactivity determined in the fraction of inositol phosphates eluted from the columns, and dpm(Ins) is that left in the monolayer following the extraction by acidic methanol. The levels of intracellular cAMP are given as pmol/mg of total cell proteins.
Equilibrium binding parameters were computed by analysis of the binding isotherms with the computer program LIGAND (23). Curves for the stimulation of phosphoinositide hydrolysis and cAMP accumulation were fitted using ALLFIT (24) to determine EC50, upper and lower asymptotes, and slopes at midrange. Standard free energy changes were computed as negative natural logarithms of the apparent association binding constants (1/Kd) and are therefore expressed in RT units. All the data are means of several independent experiments as indicated.
Simulations of Binding Data According to an Allosteric Receptor ModelTo interpret and simulate the change in ligand binding
properties induced by the mutation in NK1Rs, we used a
general multistate allosteric model, which is briefly described here.
The receptor (R) exists in n freely
interconverting states. The concentration at equilibrium of each state
(si) is given by a stability constant (j)
describing the first-order transition from any state i 1 to state i. For a sufficient description of this system, microscopic reversibility allows us to consider any state transition as
relative to an arbitrarily chosen reference state
(s0). Therefore, for any state i, the
corresponding stability constant is given by ji = [si]/[s0], and the total
concentration of the receptor existing in n states is thus
as follows (Equation 1).
![]() |
(Eq. 1) |
If a ligand (H) binds to the receptor, it will perturb its state distribution according to the degree of change that the bound molecule imparts to the stability constants of each state. Thus, in the presence of bound ligand, the stability constant of any state i is given by biji = [Hsi]/[Hs0], where the factor b indicates the extent to which the bound ligand alters the stability constant of that state with respect to the unbound receptor.
The overall equilibrium binding affinity of a ligand is the result of the binding affinities of all the states (Equation 2),
![]() |
(Eq. 2) |
To generate experimental binding isotherms representing the concentration of a bound radiolabeled ligand as a function of the total concentration of a second unlabeled ligand, the equilibrium composition of the system must be computed for n number of states and m number of ligands, starting from the equilibrium constants and the total concentrations of receptor and ligands as input parameters. This requires the simultaneous solutions of m + 1 equations describing the concentrations of all free species present (Equations 3 and 4).
![]() |
(Eq. 3) |
![]() |
(Eq. 4) |
Solutions were achieved by an iterative numerical procedure directed to minimize, for each ligand (H) present, the implicit function as follows (Equation 5).
![]() |
(Eq. 5) |
From free species, the concentrations of all bound species can be computed, and the relation between bound and total ligand can be obtained for each ligand. Binding isotherms thus generated were analyzed by LIGAND (23) to compute the apparent Kd and the corresponding equilibrium affinity (1/Kd).
To simulate "allosteric" mutations, we evaluated the effect of changing j parameters either in the absence of any other modification or also in the presence of simultaneous smaller changes in the allosteric b factors of all ligands. To do so, b factors were multiplied by an equally sized matrix of computer-generated random numbers with unitary mean, and a standard deviation equivalent, in energetic units, to one-tenth or one-fifth of the total change that was imposed on the j constants. This introduces scatter in the effect of the mutation that closely resembles that experimentally observed.
Wild-type and mutant NK1Rs were stably transfected into CHO cells, and a pair of wild-type and mutant CHO clones exhibiting comparable levels of receptor expression were selected for the binding studies presented in Table I. To investigate the pharmacological properties of the mutant receptor, we measured equilibrium binding affinities of an ample selection of structurally diverse tachykinin receptor ligands.
|
As previously observed in a transient expression system (5), differences between SP and a number of peptide analogues mark the salient feature of the mutation. In fact, while SP affinity is similar for wild-type and mutant receptors, natural tachykinins selective for other types of neurokinin receptors, such as NKA and NKB, and amino-terminal deleted analogues of both SP and NKA displayed enhanced affinity for the mutant (Table I). Based on these data, we suggested that the main effect of the double-residue replacement could be an improved interaction of the receptor for the consensus pentapeptide C-terminal sequence of tachykinins (FXGLM). The reason why such improvement was evident for heterologous tachykinins and amino-terminal deleted SP analogues, but not for SP itself, might be explained if we assume that the amino-terminal residues, by providing additional interactions at sites others than those targeted by the mutation, could hinder the improvement in affinity for the consensus sequence generated by the mutation itself.
In this study, however, we identified an additional set of ligands responding with increased affinity to the mutation, which suggests that this "message-recognition" hypothesis is not correct. First, we found that septide and its dextro analogue bearing a proline substitution in the C-terminal consensus sequence display greater affinity for the mutant receptor, just like the unmodified hexapeptide SP-(6-11). Second, three C-terminal modified analogues of SP (the methyl ester, the deamidated analogue, and the carboxy-free version of SP-(1-9)) also displayed significantly greater affinity for the mutant receptor than the wild-type receptor. Thus, responsiveness to the mutation requires neither the presence of an intact tachykinin C-terminal consensus core nor the absence of a SP-like amino-terminal sequence.
One clear consequence of the mutation is a change in selectivity of NK1R for natural tachykinin, as the ratios of affinities between SP and either NKA or NKB are much larger in the wild-type receptor than in the G166C,Y167F mutant. Because selectivity may depend partly on differences in sequence and partly on the conformational constraints that those differences impose on the common binding elements of the three peptides (25), we considered whether differences in both selectivity and flexibility among ligands might be the major factor conferring sensitivity to the mutation.
However, the comparison of the binding affinities of a wider range of tachykinin peptide analogues (Table I) indicates that the ability to discriminate between mutant and wild-type receptors is related neither to the conformational flexibility of the peptide ligand nor to its selectivity for different receptor subtypes. For example, septide and D-septide, which differ from SP-(6-11) in a constraining proline substitution at position 9 and enhanced NK1R selectivity, display similar preference for the mutation as their nonconstrained analogues. Also, sarcosine substitution at the same position of the sulfone undecapeptide analogue of SP, although increasing NK1R selectivity, does not bring enhanced affinity for the mutation. Similarly, physalaemin which does not discriminate between mutant and wild-type receptors, does not acquire selectivity for the mutant when presented as a more rigid conformational version in the neurokinin-1-selective antagonist GR 82334, nor do two spantide-like antagonists, either as a full-length undecapeptide or as a C-terminal hexapeptide analogue, display any important difference in affinity between wild-type and mutant receptors. We also tested senktide, which is a more rigid version of the hexapeptide carrying additional substitutions that establish high selectivity for NK3Rs and which does not bind to wild-type NK1R. This peptide did not acquire any apparent affinity for the mutant (no competition at 100 µM for either receptor; data not shown). Taken collectively, these data suggest that conventional structure-activity criteria do not lead to a consistent explanation of why some peptides display enhanced affinity for the mutant receptor, while some others do not.
Table I also illustrates the effect of the mutation on the binding affinities of a representative group of structurally diverse nonpeptide tachykinin antagonists. Most of them showed little difference in Kd between wild-type and mutant receptors, thus behaving just like the agonist SP. But two were exceptions since their binding affinity for the mutant receptor was diminished 5-fold (Table I, last two entries). This suggests that there might be a larger class of tachykinin ligands responding with decreased affinity to the mutation, although we have found only two members thus far. Again, the structural basis for such discrimination remains elusive. There is, in fact, very little apparent structural analogy between CP 96345, a quinuclidine derivative with species selectivity for human and guinea pig NK1R (16), and FK224, a microbial cyclic peptide that was by serendipity found to be an antagonist of NK1R (12).
Free Energy Correlations Suggest a Mutation-induced Conformational ChangeThe finding that no common structural feature satisfactorily predicts the ability of a ligand to detect the change induced by the mutation suggests that the two residues replaced in the receptor may be not primarily involved in direct docking interactions with a particular chemical moiety of the tachykinin sequence, which should be otherwise invariably present in the structure of responding ligands and absent in that of others. Alternatively, they may exert an influence on the entire conformation of the receptor, i.e. acting globally, rather than locally, in changing receptor affinity.
We thus propose a different strategy to analyze the data. The suite of ligands examined in this and previous studies can be viewed as a variation across the conformational space of the "tachykinin pharmacophore" fitting the stereochemical constraints of NK1R. The effect of such variation on binding affinity has been evaluated here for two congeneric forms of the same receptor. Thus, a comparison of binding energy changes due to modification of the ligand with those due to modification of the receptor can provide some insight into the relative extent of the two contributions.
To this end, we first scaled the free energy change resulting from
variation in ligand structure over that resulting from receptor
transmutation. To do so, SP (the highest affinity ligand of the set)
was taken as reference, and we plotted the net difference in binding
energy between each ligand and SP measured for the wild-type receptor
as a function of the net difference due to the mutation. This generated
a scattered distribution of values (Fig. 1A),
among which the points relative to ligands that have increased affinity
for the mutation form a well separated cluster. (Similarly, the two
ligands responding with decreased affinity to mutation appear to be
well separated from the cluster of non-responders, but it is hard to
draw any conclusion in the absence of additional points.) Within the
responder group, there is no significant correlation (r = 0.10, p = 0.75) between the two variables, indicating
that the extent of binding affinity of a ligand for NK1Rs
does not allow a prediction of how that affinity would change in
response to the replacement of the two residues targeted by the
mutation. The lack of overlap between responder and non-responder
groups on the x axis of the plot suggests that, aside from
the ligand-specific effect, a main component of the observed variance
is due to the change in the receptor, i.e. all members of
the responder group appear to differ from all the others primarily by a
shift in binding energy caused by the mutation.
Correlation among binding energies of ligands
for wild-type and mutant NK1Rs. A, the net
differences in free energy changes between each ligand and SP for
binding to wild-type NK1R (ordinate) are plotted
as a function of the net free energy differences between wild-type and
mutant receptors for each ligand (abscissa). Means ± S.D. (bars) were computed after converting into free energy differences (see "Experimental Procedures") the individual affinity values (1/Kd) corresponding to the averaged data and the number of experiments reported in Table I. Different symbols were
used to mark values corresponding to ligands whose affinities are
enhanced (), unchanged (
), or diminished (
) by the mutation. B, the net differences in binding energy relative to SP (see
A) measured for the mutant receptor (Net
Gmut; ordinate) are plotted
versus the corresponding values measured for the wild-type receptor (Net
Gwild type; abscissa).
The symbols identify the same groups of values shown in the scatter
plot of A. The solid lines indicate best fitting
regression lines fitted to the two sets of values, corresponding to
responders (
; slope = 1.02 ± 0.06 and intercept =
2.97 ± 0.51) and non-responders (
; slope = 1.05 ± 0.03 and intercept = 0.078 ± 0.15). The values of the two
ligands with decreased affinity for the mutation (
) were not
included in these calculations (dotted line). To evaluate whether the improvement in the description of the data given by two
separate regression lines compared with a single one traced through all
points is statistically significant, we used analysis of variance and
F statistics according to the extra-sum-of-square principle
(37). The critical F value was computed from the relation ((SS1
SS2)/(df1
df2)) × df2/SS2, where
SS1 and SS2 are the sum of squares of the
residual for the simpler two-parameter model (one regression line) and
for the more complex four-parameter model (two lines), respectively,
while df1 and df2 are the corresponding degrees
of freedom (number of fitted data minus number of total estimated
parameters). The computed value (F = 69.7; degrees of freedom = 2, 20) is highly significant (p < 10
6) and rejects the hypothesis that a single regression
can describe the data equally well. The significance of this test also
indirectly validates the assertion that the two groups of values
form distinct clusters in the graph of A, which can be
viewed, in fact, as a way to display the residuals of the regressions
in B.
To evaluate the magnitude of this shift, we plotted the relation between net binding energies measured for wild-type and mutant receptors (Fig. 1B). The data are best described by two distinct and parallel regression lines corresponding to the groups of responders and non-responders, respectively, with slopes not significantly different from unity. Again, the two ligands displaying a significantly reduced affinity for the mutation seem to belong to a third regression line, although the scarcity of points does not allow proper statistical analysis. The unitary slope of such relations indicates that wild-type and mutant receptors react to changes in ligand structure identically. In this sense, the mutation does not alter the fundamental relationship between structural change in the ligand and change in binding energy for NK1R because that would generate nonparallel lines. However, the line of ligands that respond with increased affinity to the mutation is shifted from that of non-responders by 3 RT units of free energy, indicating that all those ligands have gained, on average, such an extent of binding energy as a result of the change in the receptor. Similarly, by considering a putative line passing through the values of the two antagonists with reduced affinity for the mutant, there may be another subset of ligands displaying an opposite shift with a net loss of 1.5 units of binding energy, although its existence remains speculative.
Although the regression line explains 97% of the variance in the responder group (Fig. 1B), the residuals reflect both experimental error and ligand-specific effects of the mutation, and this analysis does not allow us to distinguish among them. We can, nonetheless, conclude that in the mechanism underlying the change in apparent affinity of the responders, the constant change affecting the whole set of ligands is much more important than smaller contributions due to changes that act differently upon each ligand.
The constancy of such change in binding energy supports the notion that the main effect of the mutation is on the conformational equilibrium of the receptor itself. It also suggests that the class of ligands that is able to "detect" such change must share a common property despite our inability to detect common features in their chemical structure. The next experiments were designed to evaluate this possibility.
Ligands Responding to the Mutation Share Selectivity in Transmembrane SignalingAlthough ligands responding to the mutation are agonists, while antagonists are the majority of non-responders, there is no straightforward relation between mutation-induced change in affinity and intrinsic activity of the ligands. In fact, non-responders also include potent agonists such as SP itself, the [Sar9]sulfone analogue, and physalaemin. Mutagenesis studies of tachykinin receptors (26-28) suggest that peptide and nonpeptide ligands of NK1R might have only partially overlapping binding sites; therefore, replacement of certain residues might only affect the peptide-selective area of the site. But here we find that the peptide SP and several nonpeptide ligands are equally insensitive to the replacement; therefore, the ability to detect the mutation cannot be ascribed to the peptidergic nature of the ligand.
There is, however, an intriguing biological property that septide and
several other ligands of the responder group appear to share. It is
known that the potency of septide in biological assays is much greater
than expected from the binding affinity measured in competition for SP
(reviewed in Ref. 29). In transfected cells, tachykinin receptors
activate at least two different signaling pathways, phosphoinositide
turnover and adenylate cyclase, through the intervention of distinct G
proteins (30). Recent studies show that in cells expressing
NK1Rs, septide, SP-(6-11), NKA, and NKB activate
Gq-mediated PI turnover with a potency similar to that of
SP, despite their low apparent binding affinity. In contrast, the
EC50 of these peptides for stimulation of intracellular accumulation of cAMP is much lower than that of SP (31). This suggests
that ligands responding with increased affinity to the mutation may
have in common the ability to activate more efficiently Gq-
than Gs-dependent signaling. We wondered,
therefore, whether C-terminal analogues of SP, which like septide
respond with increased affinity to the mutation, also exhibit
"septide-like" biological properties. To answer this question, we
compared, in the two transfected CHO clones used for the binding
studies, the relative enhancement of PI turnover and cAMP levels
induced by
[<Glu6]SP-(6-11) 2 and SP
methyl ester. Both ligands displayed little difference in potency or
maximal stimulation in activating PI turnover compared with SP in cells
expressing the wild-type receptor (Fig. 2, top left panel), but clearly differed in their ability to elicit
accumulation of intracellular cAMP (Fig. 2, bottom left
panel): [<Glu6]SP-(6-11) displayed a potency 2 orders of magnitude lower than SP in this response, while the methyl
ester was not only less potent, but also behaved as a partial agonist,
producing 65% of the maximal stimulation observed for SP. In cells
expressing mutant receptors, the relative potency of the three ligands
was not substantially changed (Fig. 2, right panels),
although the cAMP response measured in this case was small and did not
allow reliable comparisons of the concentration-response curves. In
fact, the most striking difference between the two transfected clones
is that in the mutant-expressing cell, the maximal stimulation of PI
turnover is greater, while that of cAMP accumulation is smaller
compared with the wild-type receptor (Fig. 2, right
panels).
Two important conclusions can be drawn from these experiments. First, they indicate that also C-terminal modified analogues of SP share with septide a much better efficiency to activate Gq- rather than Gs-mediated responses. Such a feature may thus be the common denominator of all ligands responding with increased affinity to the mutation. Second, they suggest that the mutation may change the intrinsic ability of the receptor to differentially activate the two G proteins since SP, a ligand that discriminates poorly between the two signaling paths, appears to induce more PI stimulation and less cAMP accumulation in mutant compared with wild-type receptor-expressing cells. This possibility was investigated further.
Differential Ability of Wild-type and Mutant Receptors to Trigger Distinct Signaling PathwaysTo study the relation between
concentration of expressed receptor and extent of activation of the two
different signaling pathways, we first used a transient expression
system. In fact, by calibrating the concentration of transfected
cDNA for the two receptors, their respective ability to activate
second messenger systems can be quantified, despite the lower levels of
expression of the mutant. The expression of wild-type NK1R
over a wide range of densities (Fig. 3A) did
not induce significant changes in PI turnover or intracellular cAMP
concentration in the absence of agonist. Similar data were obtained
with the mutant, indicating that the mutation does not result in a
detectable enhancement of ligand-independent receptor activity.
However, in the presence of saturating concentrations of SP, there was
a clear-cut difference between mutant and wild-type receptors (Fig.
3A). In fact, the unitary increase in PI turnover/pmol of
expressed receptor estimated from a linear regression of the data is
0.025 ± 0.0011 and 0.057 ± 0.0033 for wild-type and mutant
receptors, respectively. Thus, at equal molar units of expressed
receptor, the mutant produces 2-fold more PI turnover stimulation than
the wild-type receptor (Fig. 3A). In contrast, stimulation
of cAMP levels in the same cells did not show a significant difference
between mutant and wild-type receptors (Fig. 3B). Thus, the
mutation can affect the signaling efficiency of NK1R for
only one kind of signal transduction pathway selectively.
To substantiate this finding in a more "physiological" expression
system, we strived to prepare a number of stably transfected cell lines
differing in the number of expressed mutant or wild-type receptors.
Using Chinese hamster ovary cells, we obtained two series of clones
with a reasonably wide range of wild-type or mutant receptor
expression, which were analyzed for both phosphoinositide turnover and
accumulation of cAMP in response to SP. The maximal stimulation of PI
turnover in these cells was related to receptor density (Fig.
4A), but did not show a linear trend as
observed in COS cells. Nonetheless, also in this case, equal levels of PI stimulation occurred at lower expression of the mutant receptor compared with the wild-type receptor. Instead, the mutant was substantially impaired in mediating stimulation of cAMP accumulation (Fig. 4B). Direct assessment of adenylate cyclase activity
in membrane prepared from two selected clones expressing either the wild-type or mutant receptor confirmed that little
Gs-mediated stimulation of enzymatic activity can be seen
by the mutant receptor in this cell system (data not shown).
Additional studies were carried out in transfected C6 cells, a glioma
line that expresses endogenous bombesin and -adrenergic receptors,
coupled to Gq and Gs, respectively, but not
tachykinin receptors. However, only a few mutant-expressing clones were
obtained, which displayed receptor densities far below those measured
in wild-type receptor-expressing cells. Thus, to compare the best mutant-expressing clone with the wild-type receptor-expressing cells,
we divided the extent of maximal PI turnover stimulation in the
presence of SP by the molar concentration of receptors measured on the
cell surface in parallel experiments. Again, the mutant receptor
exhibited an enhanced ability to activate phosphoinositide hydrolysis
(maximal stimulation in fractional PI units/pmol of receptor was
0.025 ± 0.003 in the wild-type receptor and 0.034 ± 0.004 in the mutant receptor). There was no accumulation of intracellular
cAMP in response to SP agonists in transfected glioma cells expressing
either the wild-type or mutant receptor, although all the clones
displayed large and comparable increases in cAMP concentrations in
response to isoproterenol, similar to those measured in nontransfected
C6 cells.
In conclusion, these experiments indicate that
NK1R-mediated stimulation of PI turnover is potentiated by
the mutation in three different cell systems. The effect of the
mutation on cAMP responses is less clear: there is no effect in COS
cells; there is marked impairment in CHO cells; and it cannot be
evaluated in C6 cells, where NK1Rs do not appear to control
cAMP levels. Whether these differences across cell lines reflect a
divergent degree of cross-talk between the two signaling pathways or
differences in stoichiometry between the transfected receptor and the
pool of endogenous G -subunits remains to be investigated. We cannot draw any further conclusions at the moment, except for an obvious consideration of the importance of the surrogate expression host cell
type in the study of recombinant receptors.
The core message of such experiments, however, is that the mutation induces a small but definite shift in receptor preference for the two signaling pathways and, consequently, for the G proteins Gq and Gs, which are involved in these two modes of signal transduction. This implies that the chemical nature of the pair of residues targeted by the mutation can influence either the respective affinities of the receptor for different G proteins or the interconversion of the receptor into divergent "active" forms that have different specificity for distinct G proteins. In either case, the effect is clearly allosteric in nature because the mutated residues, given their location, cannot form direct contact sites for the G protein in the cytosolic region of the receptor.
Allosteric Linkage between Change in Apparent Affinity of the Ligands and Shift in Signaling Property of the ReceptorThere is a phenomenological link between the effects of the mutation on the binding and signaling properties of the receptor: the mutation-induced enhancement of apparent affinity is only "detectable" by agonists that preferentially trigger Gq-mediated responses, and correspondingly, the mutant receptor appears to produce more efficient Gq-mediated signaling even when it is activated by the relatively "transducer-unselective" agonist SP.
The interesting question is whether this link can be interpreted on a more mechanistic and quantitative basis. A ternary complex model describing the interaction between the receptor and two distinct G proteins may be a potential tool to interpret the shift in apparent affinities observed in this study. But it is conceptually difficult to apply such model to the binding data presented here, which were generated in intact cells incubated in isotonic buffer at 0 °C. In fact, at the intracellular concentrations of guanine nucleotides (close to the mM range), this model predicts that the fraction of bound receptor in the G protein-coupled form is negligible at equilibrium, and the measured affinity reflects that of the free receptor state (32). Indeed, the binding isotherms of SP were always consistent with a single class of sites in this study (data not shown).
However, the allosteric concept of ligand-receptor interactions predicts that macroscopic binding affinities are apparent also in the absence of the stabilizing effect of the G protein. Even a monomeric receptor with a single binding site for ligands can be thought of as existing in equilibrium among a large number of different tautomeric states, in each of which the binding affinity of that site will be different (33). The introduction of a ligand in this system shifts the equilibrium toward those conformations that have the highest affinity for that ligand. Therefore, any macroscopic binding constant can always be interpreted as the aggregate of two energetic components: one due to the second-order association process taking place between the two reacting partners and the other resulting from the displacement that the process itself imposes on the intramolecular equilibria within the receptor (34).
Accordingly, even if a mutation does not directly alter "docking" residues within the ligand-binding region of the receptor, it can nonetheless modify "allosterically" the binding affinity of a group of ligands through two distinct mechanisms: (a) by changing the stability of the state (or subset of states) that is commonly induced by that group of ligands or (b) changing the extent to which each ligand alters the stability of those states. The first will produce an equal shift to the affinities of all those ligands; the second affects the affinity of each ligand differently. More realistically, we should think that any mutation may always produce changes by a combination of both effects, but if the change in stability constants of the states is proportionally greater than those affecting the allosteric factors of each ligand, we shall detect a constant shift in the analysis of free energy correlations presented in this study.
To illustrate this by a simplified numerical example, we simulated an allosteric receptor existing in a minimum of four interconverting states (s0 through s3), interacting with four groups of ligands, each consisting of 50 distinct members. Within each group, ligands differ in affinity for the reference state (s0), but share the ability to stabilize preferentially the same receptor state.
To generate the mutation, we applied a net overall free energy change
of 4.5 RT units unequally distributed over the stability constants of the four states so that the constant of
s3 was enhanced by the equivalent of 3 RT free energy units, that of s1 was
diminished by 1.5, while the others remained unaltered (Fig.
5, upper panels). The mutation was simulated
in two ways. First, we only allowed a change in stability constants
(j parameters), and then we also allowed smaller random
changes in the allosteric factors of each ligand (b
parameters). The results of theoretical modeling are plotted in Fig. 5,
exactly as it was done for the experimental data of Fig. 1. Two points
are clear from such simulations.
First, a mutation altering the stability constants of discrete states of the receptor can allosterically change the apparent affinities of all ligands that stabilize those targeted states. If we introduced scatter in the data by allowing the mutation to also change randomly the allosteric factors of ligands, the binding energies of each class of such ligands formed well delimited clusters in plots of wild-type receptor binding energy as a function of net mutation-induced change (Fig. 5, left panels). In contrast, the binding energies of ligands inducing states that were not targeted by the mutation could not be discriminated at all.
Second, the magnitude of change in the stability constants of the targeted states is faithfully reflected in the extent of displacement it imposes on the relation between binding energies for mutant and wild-type receptors. In this simulated example (Fig. 5, right panels), the clusters of binding energies corresponding to ligands selective for s1 and s3 (the stability of which is diminished and enhanced, respectively, by the mutation) generated regression lines with unitary slopes that are shifted 1.5 units to the left and 3 units to the right, respectively, from those of all the other ligands.
According to such a model, the effect of the G166C,Y167F mutation in NK1R can be interpreted as resulting from the enhancement of the intrinsic stability of a state (or collection of substates) that is induced by septide-like ligands and heterologous tachykinins, but not by SP, its full-length analogues, or various antagonists. The same mechanism may also explain the linkage between the change in signaling properties of the receptor and the shift in affinities for ligands sharing increased selectivity for PI signaling if we imagine that the set of functional states induced by septide-like ligands (and favored by the mutation) may preferentially signal via Gq.
In summary, it was shown previously (5) that replacement in NK1R of two NK2R residues located near the extracellular edge of TMD-4 reduces the selectivity of the receptor. We show in this study that the mechanism of increase in affinity for heterologous tachykinins (which underlies such change in selectivity) is not conceivably the result of "local" interactions between the transmutated side chains and specific structural elements of the responding ligands. It consists instead of a broader conformational change that improves the apparent affinity of transducer-selective septide-like tachykinins and also shifts the transducer preference of the receptor itself.
This deduction is supported by three kinds of experimental evidence. First, the analysis of the binding energies of 26 ligands for wild-type and mutant receptors indicates that within the modifications due to the mutation, the constant shift in apparent affinity affecting an entire class of structural heterogeneous tachykinins bears much greater weight than specific changes affecting differently the affinity of each individual ligand. Second, C-terminal modified analogues of SP exhibit the same transducer selectivity as the other members of the group of ligands responding with increased affinity to the mutation, suggesting that this property is shared by all those ligands. Third, the relation between receptor density and biological response evaluated in three different mammalian expression systems indicates that G166C,Y167F NK1R triggers PI signaling more efficiently than the wild-type receptor upon activation by SP.
We interpret these findings by proposing that the G166C,Y167F mutation may enhance the intrinsic stability of a functional substate of NK1R that is preferentially induced upon binding by septide-like agonists, but not by SP-like agonists or nonpeptide ligands. Simulations based on a multistate allosteric receptor model indicate that if a mutation alters individual stability constants of receptor states to a greater extent than the allosteric effect of each ligand, the pattern of change in apparent affinities closely matches that observed in this study. The model therefore supports this interpretation and also provides theoretical background to the analysis of free energy correlations proposed in this paper. Such an approach may be generally useful in the study of site-directed mutagenesis of G protein-coupled receptors and the structure-activity relation of ligand-receptor interactions.
A mutation resulting in selective enhancement of affinity similar to that reported here was described recently for human NK1R (35). Alanine replacement of the highly conserved Tyr-216 located near the cytosolic end of TMD-5 increases the affinity for the heterologous tachykinins NKA and NKB, septide, and C-terminal modified analogues of SP, but not for nonpeptide antagonists and SP itself (35). The increase in binding affinity was attributed to a conformational change in the receptor rather than to a direct interaction between the bound ligands and the mutated side chain (which is very unlikely given the location of the targeted residue in that case). Unlike the G166C,Y167F mutation, the Y216A mutation abolishes PI turnover stimulation, although the possibility that cAMP responsiveness may have been left unchanged was not investigated (35). Considering that TMD-4 and TMD-5 are tethered via a common extracellular loop but are connected to distinct cytosolic segments of the sequence, the finding that different mutations located at opposite edges of the two helices induce similar conformation-mediated changes in ligand binding but divergent effects on receptor activation is provocative. It suggests that the effect on binding of the two mutations may be mediated by a common mechanism, such as, for example, long-range influences producing in both cases a favorable orientation of the second extracellular loop, whereas the differences in receptor activity may reflect the diverse roles that the two intracellular loops play in the process of conversion of the receptor into the active form and in the formation of a docking configuration for the intervening G proteins. Additional studies based on the construction of double-mutant cycles (36) may be useful to evaluate whether there is additivity or cooperativity between perturbations imposed in these two distinct sites of the NK1R transmembrane bundle.
One implication of the results of both studies (this study and Ref. 35) pertains to the mechanism underlying the selectivity of the receptor for different endogenous ligands. As proposed previously (5, 25, 26, 35), the different affinities of tachykinin receptors for their natural ligands may be the result of a selection of optimal matching conformations between common structures in both ligands and receptors, rather than the product of specific local interactions between the variant "address" regions of the peptides and complementary divergent residues in the sequence of the receptors. If selectivity depends on ligand-induced selection of distinct receptor conformers, any mutation that alter the tautomeric equilibria between receptor conformations must also change selectivity. The results presented here confirm this prediction and use the concept of allosteric transition (33) to address this notion on a more theoretical and quantitative basis. In addition, we show here an uncommon example of how receptor mutagenesis can help to identify functional differences among agonists (31) that cannot be easily rationalized by classical pharmacological theories of drug efficacy. Our theoretical analysis indicates that mutations affecting the state distribution of the receptor can dissect discrete clusters within the variation in binding energies of ligands, just like ligands select discrete states within the conformational space of the receptor.
A second implication of this study concerns the mechanism of G protein promiscuity in seven-transmembrane receptors. Here we show that two residues located near the extracellular part of the receptor (the side chains of which are therefore not likely to provide docking sites in the cytosolic binding region for the G protein) are nonetheless important in determining G protein preference since their replacement produces a small but detectable shift toward Gq-mediated PI signaling. This suggests that also G protein selectivity may depend on a conformation-driven selection mechanism as discussed above for the ligands. It could also explain how septide-like agonists, although sharing the same binding site of SP on the receptor, can trigger signals that differ distinctively in transducer selectivity.
We are grateful to Drs. Schwartz, Manzini, Matsuo, and Edmonds-Alt for providing neurokinin-1 ligands.