Abstract
A restricted version of the ternary complex model for receptor-G protein complex formation has recently been proposed. Known as the two-state model, this model proposes that in the context of agonist and G protein interactions, only two thermodynamic states exist for the receptor: active (R*) and inactive (R). One form of this model suggests that only the R* state of the receptor is capable of interacting with and subsequently activating G proteins. We directly tested the kinetic aspects of a strict two-state receptor model in a cell line containing the native β2-adrenergic receptor that is capable of inducing Gs expression. We examined adenylyl cyclase activity in the presence of limiting GTP levels and concluded that there exists a different rate of heterotrimer dissociation (i.e., HR*G yields HR* + G*) for different β2-agonists. This finding is inconsistent with a strict two-state model in which R* is a characteristic of the receptor that is independent of the identity of the agonist. It implies that agonist activation of adenylyl cyclase is more complicated than a simple two-state model.
The concept of agonist-induced receptor conformation has been the subject of much debate. Although the order of the interaction between agonists, receptors, and G proteins is relatively undisputed, the stoichiometry and nature of the proteins during activation remain very much in question. One model in particular that seems to be gathering widespread acceptance for G protein-coupled receptors is the two-state model for receptor activation (1-3). The two-state model is a restricted version of the ternary complex model (4) and has been extremely useful for explaining the experimental observations involving β2-adrenergic receptor activation of adenylyl cyclase. Moreover, its stark simplicity makes it easy to set up arguments that attempt to predict the behavior of other features of the system.
In physicochemical terms, the two-state model cannot be exactly true. Molecules exist in a very large number of conformational states, and trying to describe the thermodynamic properties of a receptor with just two conformational states must be an approximation. That given, the usefulness of the model depends on just how reliable its predictions are. In terms of receptor-G protein interactions, a strict two-state model is as shown here: Equation 1where G is the G protein, R is receptor, and H is agonist. It is supposed that the receptor must exist in one of two conformations (R and R*), with both H and G binding with greater affinity to the R* form. Given this condition, there automatically is cooperativity in binding of the G and H to the R (which explains the GTP shift). Moreover, the different efficiencies of different agonists to promote adenylyl cyclase activation is easily explained by supposing that stronger agonists favor binding to the R* form to a greater extent than do weaker agonists and hence more readily form the HR*G ternary complex. In addition, inverse agonists tend to stabilize the inactive R state (5, 6). It must be emphasized that in a strict two-state model, the affinity of G for R is identical to that for HR and its affinity for R* is identical to that for HR* because R and R* represent unique defined conformational states (or collection of states). Given the above set-up, it is possible to develop the thermodynamics of the system in terms of the relative concentrations of the two conformations: when free and when bound to G protein, to agonist or to both.
The classic development of the thermodynamics of the H, R, and G interactions in terms of the ternary complex model (4) shown below is not a one-state but rather a multistate approach. The use of a different dissociation constants for the binding of G to HR than for the binding to R is an acknowledgment that the receptor in HR is in a different conformational state than in R Equation 2In this equation, therefore, unlike that for the two-state model, R, HR, and so on represent the totality of possible conformations and not a single conformation or a defined group of conformations in a constant proportion. Thus, although it is superficially simpler than Eq. 1, Eq. 2 is thermodynamically complete and, unlike Eq. 1, it is guaranteed to be correct if the ternary complex model is correct as drawn. Eq. 2 may be compared to a true one-state model shown below in Eq. 3. In this equation, because the receptor can exist in only one conformation, there can be no cooperativity of binding between G protein and agonist. Equation 3The scheme proposed by Samama et al. (1) (shown below in Eq. 4), although described in that publication as a two-state scheme, is in fact a multistate model. The fact that different affinities are used for the binding of H to the R* and the R*G state means that a different conformation is assumed for R* when bound and when not bound to G protein. In other words, there are at least three states. The use of this scheme enables one to be completely accurate with the thermodynamics, but it prevents one from making the simplistic twostate arguments, which are the attraction of the model. Equation 4For the two-state model to be useful in successfully making predictions not possible by the classic thermodynamic approach, additional assumptions beyond the schemes discussed above are necessary. If the assumptions are limited to proposing that agonists bind more strongly to the conformation that also more strongly binds G protein, then no argument can be drawn beyond those already possible for the classic thermodynamic scheme. Similarly, negative agonists that reduce the amount of adenylyl cyclase activity below the basal level found in the absence of ligand can be equally well explained (and predicted) by the classic thermodynamic approach. Where the strict two-state model potentially offers advantages is in quantification of predictions and as a potential for a quantitative relationship between thermodynamic properties (e.g., GTP shifts) and kinetic features (e.g., efficiencies of agonists in adenylyl cyclase activation).
In the current study, we investigated the specific hypothesis that the fraction of ternary complex in the activating conformation (i.e., in the HR*G form) determines the rate at which GDP/GTP exchange occurs. In the strict two-state model, the HR*G conformation is the same for all agonists, which differ only in the ratio of HR*G to HRGtotal [where HRGtotal = HRG (inactive R) + HR*G (active R)]. The rationale for the test of this hypothesis is to determine the rate of breakdown of the ternary complex for four different β2 agonists. It will be shown that in all the cases the thermodynamics requires that >85% of the ternary complex must be in the HR*G rather than the HRG form. Proportionality of the rate of breakdown for the ternary complex with the fraction of ternary complex that is HR*G would therefore require that all the rates be similar. This did not occur. We will therefore conclude that the two-state model in its strictest form cannot be applied in all cases.
Experimental Procedures
Materials.
Molecular biology reagents, Dulbecco’s modified Eagle’s medium, and geneticin were obtained from GIBCO BRL (Gaithersburg, MD). Tris Base, GTP, and guanosine-5′-O-(3-thio)triphosphate were from Boehringer-Mannheim Biochemicals (Indianapolis, IN). [α-32P]ATP, Na-125I, and [2,8-3H]cAMP were from DuPont-New England Nuclear (Boston, MA). Dexamethasone and the remaining reagents from Sigma Chemical (St. Louis, MO).
Cell culture of an inducible Gsα cell line.
The establishment and characterization of a stably transfected murine S49 cyc− T cell lymphoma capable of inducing Gs protein after dexamethasone treatment have been previously described (7). Briefly, the S49 cyc− cell line, lacking Gsα mRNA and protein (8), was electroporated with the 7.7-kb pMMTV · Gsα neo vector (a generous gift from J. Gonzales; described in detail in Ref. 9). The vector contains the cDNA encoding rat Gsαlong linked downstream of the dexamethasone-inducible mouse mammary tumor virus long terminal repeat promoter. The vector also contains the selection marker for neomycin resistance, neomycin phosphotransferase, constitutively driven by the human β-globin promoter.
Transfected drug-resistant cell lines (S49*cyc−) were maintained in stock tissue culture flasks (Corning Glassworks, Corning, NY) at 37° in HEPES-buffered Dulbecco’s modified Eagle’s medium supplemented with penicillin, streptomycin, 10% heat-treated horse serum, and 200 μg/ml geneticin to maintain selective pressure. Cells were expanded for Gsα gene induction experiments into eight individual preconditioned 2-liter roller bottles (Corning) from a single stock source. The final cell density was 1 × 106 cells/ml when the volume was made up to 2 liters with fresh Dulbecco’s modified Eagle’s medium plus 10% horse serum media. Gsα protein induction was initiated when the S49*cyc− cells were incubated with 5 μmdexamethasone (final, added from a 10 mg/ml stock prepared in 95% ethanol) for times of 1–24 hr. Transfected control cells underwent no treatment.
Membrane preparations.
Cell membranes were prepared and isolated as follows. Cells were washed twice with an excess of buffer A (137 mm NaCl, 5.36 mm KCl, 1.1 mmKH2PO4, and 1.08 mmNa2HPO4, pH 7.2) by centrifugation at 600 × g. The cells were then resuspended in ice-cold cell lysis buffer B (20 mm HEPES, 150 mm NaCl, 5 mm NaH2PO4, 1 mm EDTA, and 1 mm benzamidine, pH 7.4; buffer B also contained 10 μg/ml trypsin inhibitor and 10 μg/ml leupeptin to protect Gsα from possible proteolysis) and placed in a Parr bomb (at 500 p.s.i.) for 25 min. The disrupted cells were centrifuged for 5 min at 600 × g to pellet nuclear debris. The supernatant was layered onto a 23% and 43% sucrose step gradient in HE buffer (20 mm HEPES and 1 mm EDTA, pH 8.0) and centrifuged at 25,000 rpm in a Beckman Instruments (Columbia, MD) SW 28 rotor for 45 min at 4°. The membrane fraction was collected as a band at the sucrose interface. The membranes were immediately frozen in liquid N2 and stored at −80°. Membrane concentrations were determined using the BioRad (Hercules, CA) assay (10).
β2 Agonist competition binding: GTP shift analysis.
GTP shift binding analyses were used to analyze the total amount of ternary complex (HRGtotal) present in the S49*cyc− membrane preparations containing different Gs levels as [GTP] → 0. Binding analyses were carried out in 500-μl reactions with a single 80 pm concentration of the radiolabeled β2 antagonist 125I-CYP [prepared according to the protocol of Barovsky et al. (11) with modifications by Hoyer et al. (12)] in the presence of increasing β2 agonist concentrations and the following final concentrations of the reagents: 1 mm EDTA, pH 7.4, 20 mm HEPES, pH 7.4, 10 μm phentolamine, 0.3 mm MgCl2, and 20–50 μg of cell membranes (diluted with HE, pH 8.0). The binding reaction was conducted in the presence and absence of 10 μmguanosine-5′-O-(3-thio)triphosphate for 55 min at 30°. In some experiments, binding was conducted under adenylyl cyclase assay conditions (see below) to examine the effect of varied GTP concentrations on the area of the GTP shift. In all experiments, the nonspecific 125I-CYP binding was determined in the presence of 10 μm alprenolol. Reactions were terminated with the addition of 2.5 ml of ice-cold stop buffer (50 mm Tris-Cl, pH 7.4, 10 mm MgCl2) followed by the immediate filtration of the solution [with the use of a Millipore (Bedford, MA) 1225 vacuum filtration apparatus] onto Whatman (Clifton, NJ) GF/C filters. The reaction tubes were rinsed with an additional 2.5 ml of ice-cold stop buffer. The GF/C filters were washed four or five times with 2.5 ml of stop buffer, dried, and removed to scintillation vials, in which 125I-CYP activity was counted on a Beckman Gamma 4000 System Counter for 1 min. Activity, measured as counts per minute, was converted into units of fmol/mg with the aid of a spread-sheet program using Lotus 1-2-3. The specific binding of 125I-CYP for each concentration was determined as mean values of triplicate measurements for total binding less the mean values of triplicate measurements for nonspecific binding. Data were plotted as log [agonist] versus the normalized fraction of agonist bound and were subsequently analyzed using our recently developed Scatchard method (15) to generate agonist affinity constants (to R and RG, respectively) and [Gs]total and HRGtotallevels. The use of two-component nonlinear regression analyses (GraphPAD) to analyze the −GTP curve for HRGtotal levels resulted in values similar to those obtained according to the Scatchard method (±2.5%). The area for the GTP shifts were determined by the trapezoid method. Comparisons of data ± standard deviations were performed using single-factor analysis of variance (p < 0.05).
Adenylyl cyclase assays.
Adenylyl cyclase assays, conducted with saturating agonist levels and decreasing GTP concentrations, were used to compare the rate of the ternary complex breakdown for different agonists. The premise for using these assays is described in greater detail below and in Discussion. Briefly, it is known experimentally that at normal (high) GTP concentration, the formation of the heterotrimer complex (HRGtotal) is rate limiting, and once formed, it interacts with GTP extremely rapidly (13). For a strict two-state model, the rate at which the HR*G complex dissociates should be independent of the agonist used to stabilize the R* conformation; therefore, the rate of HRGtotal breakdown should be proportional to the fraction of HR*G formed by each agonist. The rate of HR*G dissociation was made rate limiting by reducing GTP and allowing the HRGtotal complex to accumulate.
Adenylyl cyclase activity was examined in S49*cyc−membranes containing different Gs levels using the conditions for adenylyl cyclase assays as described by Clark et al. (14), with some modifications. Briefly, adenylyl cyclase activity in plasma membranes (∼0.2 μg/μl, final) was brought to steady state for 3 min at 30° in a 1800-μl volume in a 15 × 100 mm borosilicate glass tube with the following reagents (given as final concentrations): 40 mm HEPES, pH 7.7, 1 mm EDTA, 0.3 mm MgCl2, 8 mm creatine phosphate, 16 units/ml creatine phosphokinase, 0.05 mm ATP, and 0.1 mm3-isobutyl-1-methylxanthine. The preincubation period was necessary to bring the HRGtotal levels for each agonist to equilibrium. [For any agonist to be effective, the t 1/2 of HRGtotal formation cannot be less than thet 1/2 for adenylyl cyclase inactivation (t 1/2 = 15 sec for k −1in S49 cells (16)]. Each assay tube consisted of a different GTP concentration (0, 10, 30, 100, and 300 nm, and 1, 3, and 10 μm), so that for a single agonist, a series of eight independent time courses was required to produce a single GTP dose-response curve. The final concentrations of agonist required to saturate the β2-adrenergic receptors (i.e., occupy >90% of receptor sites with agonist) were 10 μm for epinephrine (Kd = 1100 nm), 2.2 μm for isoproterenol (Kd = 188 nm), 2.8 μm for fenoterol (Kd = 400 nm), and 13 μm for dobutamine (Kd = 1400 nm).
Time course assays, spanning a 3-min period, were initiated with the addition of 40 × 106 cpm [α-32P]ATP in a 200-μl aliquot containing the concentration of reagents described above. Time course assays were selected over a typical 10-min adenylyl cyclase assay to examine the linearity of the response at low GTP concentrations. The reactions were quenched every 18 sec by removing 100-μl aliquots into 500 μl of ice-cold stop buffer. The isolation and determination of [32P]cAMP activity were performed as described by Salomon et al. (15). GTP-dependent adenylyl cyclase dose-response curves were produced using linear regression analysis (GraphPAD, San Diego, CA) to generate the slopes (i.e., cyclase activity) from the linear cAMP accumulation versus time plots. These data were subsequently transformed to Eadie-Hofstee plots to examine the rate of HRGtotal breakdown (see below). All comparisons of data were performed using single-factor analysis of variance (p < 0.05, Excel).
Calculating the limiting HRGtotal and minimum HR*G levels.
The limiting levels of HRGtotal, as [GTP] → 0, were calculated from the Gs:R ratios obtained from GTP shift analyses for each β2 agonist. The Gs:R ratios, representing the maximum level of receptor-G protein interaction in the presence of saturating agonist, were obtained for each S49*cyc− membrane preparation containing variable Gs levels using a newly developed Scatchard method (16, and data not shown); however, analysis of HR versus HRGtotal plots for the −GTP curves with a rectangular hyperbola will yield similar results (±2.5%). With the Gs:R level for an agonist, one can determine the limiting HRGtotal levels by using the following relationship (16):
Equation 5where K
RG represents the dissociation constant between HR and G. The values used forK
RG (unitless when [R]total is set to unity and all concentrations are expressed as fractions or multiples of [R]total) were 0.007 ± 0.002 for epinephrine, 0.011 ± 0.006 for isoproterenol, 0.007 ± 0.003 for fenoterol, and 0.029 ± 0.016 for dobutamine, and [Gs]total:R ranged from 0.246 to 0.692 when epinephrine was used in GTP shift studies (16). Under saturating agonist conditions, [HR] ≅ [R]total, and the limiting values for HRGtotal are determined by solving the following quadratic equation:
Determination of the rate of HRGtotal breakdown.
In a strict two-state model, the rate of HRGtotal breakdown should be proportional to the fraction of HR*G formed by each agonist (i.e., once HR*G is formed the rate of breakdown should be agonist independent). The experimental rate of HRGtotal breakdown was determined by fitting the Eadie-Hofstee plots (with reversed axes so that the x-axis represented υ and the y-axis represented v/s) for cyclase activity to a quartic solution that kinetically describes the GTP-dependent activation of adenylyl cyclase (Appendix ). A formal derivation reveals that the y-axis in the present case is
Results
Examination of adenylyl cyclase response with varying GTP concentrations.
The kinetics of adenylyl cyclase activity, under reduced GTP concentrations, were examined in S49*cyc−membranes containing variable Gs protein levels as described in Experimental Procedures. At low GTP levels, the system allows heterotrimeric (HRGtotal) complexes to accumulate, which results in a change in the rate-limiting step of adenylyl cyclase activation from the rate of heterotrimer formation to the rate of active heterotrimer dissociation. The build-up of heterotrimer complexes at low GTP concentrations was confirmed by conducting GTP shift binding analyses with subsaturating GTP levels (4) using adenylyl cyclase assay conditions and demonstrating that the area of the GTP shift was proportional to the level of GTP in the assay (data not shown).
Fig. 1 compares the GTP dependence for adenylyl cyclase activity in the presence of saturating concentrations of four different β2 agonists. In the absence of GTP, there was a small but detectable adenylyl cyclase activity, which increased from 1.4 ± 0.4 pmol/min/mg for control membranes with minimal Gsexpression to 6.73 ± 2.2 pmol/min/mg at high Gslevels. This increase in spontaneous activity presumably reflects the increased Gs levels and the presence of endogenous guanine nucleotides associated with the membrane preparation. Under saturating agonist conditions, the presence of GTP resulted in increased response until the response became saturated at ∼3–10 μm GTP. Linear regression analysis of the curves indicated that adenylyl cyclase activity was always linear with time (r > 0.9), even under low GTP conditions, during the 3-min course of the assay.
The magnitude of the adenylyl cyclase response was dependent on the increasing Gs levels in the membrane and the type of agonist used to stimulate activity. The rank order of efficacy (epinephrine ≥ isoproterenol > fenoterol > dobutamine) can be seen clearly in a comparison of the normalized adenylyl cyclase activity with respect to increasing GTP concentrations (Fig. 2). Fig. 2 was generated by computing the slopes of the individual linear responses shown in Fig. 1 using linear regression analysis for each GTP concentration. Sigmoid nonlinear regression curve analysis of the GTP response curves also revealed some agonist dependency of the EC50 value for GTP, but in most cases the EC50 value for GTP was ∼100–150 nm. The data shown in Fig. 2 were transformed to Eadie-Hofstee plots to analyze the rate of HRGtotalbreakdown (see below).
Determination of the limiting HRGtotal, HR*G, and rate of HRG breakdown for each agonist.
The test of a true two-state model for receptor activation lies in the assumption that the experimental rate of HRGtotal dissociation should be proportional to the fraction of activated HR*G complexes. This is predicted to be so if the R* conformation is independent of the type of agonist. Fig. 3 compares the relative levels of HRGtotal, the minimum HR*G, and rate of ternary complex breakdown for different β2 agonists.
The limiting level of heterotrimer complex (HRGtotal) formed in the absence of GTP was calculated for each agonist using GTP shift binding assays (Ref. 16 and data not shown) as described in Experimental Procedures. With a value of 1.0 for epinephrine, the normalized HRGtotal values (± standard error) for isoproterenol, fenoterol, and dobutamine are 0.92 ± 0.04, 0.92 ± 0.04, and 0.82 ± 0.09, respectively. The HRGtotal calculated levels for isoproterenol, fenoterol, and dobutamine were consistently less than that for epinephrine regardless of the Gs protein level in the plasma membrane (data not shown). An analysis of variance yielded no statistical difference for the effects of Gs protein levels on the measured differences in HRGtotal levels between the agonists and epinephrine. These data demonstrate that the HRGtotal levels formed in the absence of GTP seem to be relatively similar from the strong agonist epinephrine to the weak agonist dobutamine.
The values for the minimum HR*G level for each of the agonists were calculated as described in Appendix and Experimental Procedures. Fig. 3 demonstrates that the minimum fraction of active HR*G heterotrimer complexes makes up a significant proportion (>85%) of the HRGtotal pool. The HR*G levels for each agonist were normalized relative to the HR*G levels for epinephrine. The relative minimum HR*G values (± standard error) for isoproterenol, fenoterol, and dobutamine are 0.920 ± 0.001, 0.915 ± 0.002, and 0.764 ± 0.006, respectively. The 15% difference between HR*G and HRGtotal for dobutamine is thought to be due to the difficulty in measuring K d2, the high agonist affinity component; it is highly dependent on the HRG 171 HR*G distribution (in the HRGtotal pool) in that even a small contribution of HRG (complex with inactive R) leads to a substantial affect on the measured K d2. The use of a two-component nonlinear curve regression analysis to obtain a value for the high affinity binding component (data not shown) revealed that the difference between K d1, the low affinity component, and K d2 averaged 28 ± 21-fold for dobutamine (16 experiments), 215 ± 70-fold for fenoterol (16 experiments), 232 ± 82-fold for isoproterenol (eight experiments), and 400 ± 144-fold for epinephrine (24 experiments). The K d2 values determined in this fashion also seemed to decrease with increasing Gs levels (r = −0.94), which may reflect a complicated receptor-G protein stoichiometry or the difficulty in determining the high agonist affinity using a two-component curve-fitting analysis.
A comparison between the minimum HR*G levels and the relative rate of ternary complex breakdown is also illustrated in Fig. 3. The rate of ternary complex breakdown was determined from they-intercept value of Eadie-Hofstee plots representing GTP-dependent adenylyl cyclase activity for each agonist (see above). The Eadie-Hofstee plots were used to fit a quartic solution describing receptor-G protein/adenylyl cyclase interactions in the presence of decreasing GTP (Appendix ). Using the conditions for fitting as described in Experimental Procedures, the x-axis represented the number of active G* s-C complexes, a number that is limited by C in this system (16). Therefore in practice, when Gs ≤ C, the plots were linear; when Gs ≥ C, the plot tails downward as it approaches the x-axis.
In our simulations of fitting the quartic solution to adenylyl cyclase response data (data not shown), only a change ink 2 levels could result in differenty-intercepts for agonists of different efficacy (i.e., different k 1 values). Using the quartic solution to analyze experimental data revealed that the rates of breakdown (normalized average ± standard error) for isoproterenol, fenoterol, and dobutamine were 0.906 ± 0.005, 0.738 ± 0.015, and 0.505 ± 0.011, respectively. No statistical significance could be demonstrated for the effects of Gsprotein levels on the measured differences in the rates of heterotrimer dissociation between the agonists and epinephrine. In addition, no significant difference could be demonstrated between the minimum HR*G levels and the rate of heterotrimer breakdown for either epinephrine or isoproterenol. However, a significant difference could be demonstrated for fenoterol (p = 4 × 10−3, eight experiments) and dobutamine (p = 1.5 × 10−4; eight experiments). These data demonstrate that the rate of breakdown is not similar for different β2agonists (even when corrected for disparities in HR*G levels between the agonists), implying that agonist activation of receptors is more complicated than a simple two-state model.
Discussion
Two-state receptor activation model and the rate of HRGtotal breakdown.
The extended ternary complex model has provided an extremely useful thermodynamic description of agonist/receptor-G protein interactions (1-3). This model proposes that the receptor exists in one of two conformations: an inactive state, R, which displays low agonist affinity and presumably does not couple to G (2), and an active state, R*, which displays high agonist affinity and couples to the G protein. The inactive state, R, can isomerize into the active state R*, either spontaneously or with the presence of an agonist. The agonist then stabilizes the R* form, and in general it is accepted that stronger agonists have a greater ability to stabilize the active state.
The two-state model of receptor conformation offers advantages over the more general classic ternary complex model only if it can make verifiable predictions that the latter cannot. There are no qualitative predictions that can be made by the two-state model that are not also obvious from the classic model. The advantage, if there is one, is to use the two-state system to make quantitative predictions based on the amounts of the active form (R*, HR*, or HR*G) present in the system. Essentially, this devolves to relating the quantities (or the relative quantities) of these forms to effector activation or to the steady state of effector activation under conditions in which the nature of the agonist or antagonist is changed.
The scheme for the two-state ternary complex model as illustrated by Samama et al. (1), in fact, depicts a multistate model in which
1) We give measurements that show the actual accumulation of ternary complex at [GTP] → 0 for all the agonists used. In combination with 2) below, we can therefore calculate the maximum differences between the agonists for the steady state concentrations of [HR*G].
2) If the two-state model is an adequate description of the thermodynamics of ternary complex formation, then all of the agonists used are primarily in the HR*G state rather than in the HRG (R inactive) state in the ternary complex (HRGtotal).
3) We measure the adenylyl cyclase activity for all agonists as [GTP] → 0. As [GTP] → 0, the rate of breakdown of the ternary complex is the rate-determining step in adenylyl cyclase activation. Measurement of adenylyl cyclase activity under these conditions therefore gives estimates for the relative rates of ternary complex conversion to active G* s.
4) The strict two-state model requires proportionality between the amount of HR*G in the system at steady state and the rate of conversion to active G* s.
It is known experimentally that at low values of [GTP], the ternary complex accumulates. This is the basis for the theory of the GTP shift for agonist binding (4). In any process involving more than a single step, the intermediate species before the rate-determining step accumulates; therefore, it can reasonably be supposed that when GTP levels are very low, the rate at which adenylyl cyclase is activated becomes directly dependent on the rate at which the ternary complex breaks down. In the current report, however, we performed a complete kinetic analysis over the full range of GTP concentrations from very low to saturating concentrations of GTP. A complete theoretical analysis for the relationship between GTP concentration and adenylyl cyclase activity is given for the shuttle model of cyclase activation (Appendix ). However, although the overall equation is somewhat different when developed for a precoupled model (18), the conclusion that the rate of breakdown of the ternary complex can be inferred from the adenylyl cyclase activity is the same in both cases. It should be pointed out that the exact form of the relationship between adenylyl cyclase activity and GTP concentration is not crucial to the test of the two-state model performed in this report. The solution for that relationship was derived using a Gs-to-C shuttle model as the primary mechanism for adenylyl cyclase activation because it is interesting in its own right and because it allows the mathematical demonstration that as [GTP] → 0, the adenylyl cyclase activity is directly proportional to the rate of ternary complex breakdown (Appendix ).
Experimental conditions were chosen such that similar levels of HRGtotal were used to stimulate adenylyl cyclase activity for four different β2-adrenergic agonists: the strong agonists epinephrine and isoproterenol, the moderate agonist fenoterol, and the weak agonist dobutamine. If a strict two-state model were manifest for receptor conformation, then similar levels of HR*G are predicted for each agonist when examining adenylyl cyclase activity under limiting GTP concentrations. Appendix provides the formal thermodynamic justification of the statement above by showing that it is possible to use GTP shift binding data to calculate the minimum percentage of HRGtotal that must be in the HR*G state. Fig.3 reveals the similarity in the HR*G levels between the agonists and shows that the weak agonist dobutamine is capable of producing >80% of the HR*G complexes as epinephrine under saturating agonist conditions. Because the percentage is very high even for the weak agonist dobutamine, the subsequent arguments are facilitated.
If the properties of the ternary complex depended proportionately only on the fraction in the active HR*G form, then there should be little difference between the ternary complex that included the weak agonist dobutamine, the stronger agonist fenoterol, and the very strong agonists epinephrine and isoproterenol. It would then be expected that the rate of HRGtotal breakdown would be similar for all of the agonists. Significant differences in the rates of breakdown would demonstrate that a simple two-state model is inappropriate for this process. Fig. 3 demonstrates that the relative rates of heterotrimer breakdown for fenoterol and dobutamine are significantly different from the minimum HR*G levels determined for each of the two agonists. This suggests that the rate of HRGtotal breakdown is dependent on the type of agonist, a finding that is inconsistent with a single conformation for the activated state of the receptor.
It is unlikely that the agonist-dependent differences in the rate of heterotrimer breakdown are due to 1) the ability of a particular agonist to form the ternary complex or 2) the agonist-dependent variation in the rates of agonist association and dissociation to the β2-adrenergic receptor. Our experiments were designed to minimize the effects described above by bringing adenylyl cyclase activity to steady state during a 3-min incubation period before the addition of labeled ATP and subsequent data collection. Therefore, problems resulting from the delay in reaching full activity on the addition of stimulating agonist do not arise. Similarly, the preincubation ensures that agonist binding has also attained a steady state before measurements are taken. Thus, the rates at which agonist binds and dissociates from the receptor should also not be a significant factor.
A possible explanation for these data that would preserve a strict two-state model requires the selective activation of the inhibitory heterotrimeric G protein, Gi. If Gi were preferentially activated by dobutamine and fenoterol, then the reduced adenylyl cyclase activity at low GTP, under stimulation by those agonists, would be the result of Gi-mediated inhibition rather than a slower rate of breakdown of the ternary complex. In the current circumstances, this is unlikely because Gi-coupled α-adrenergic receptors have not been demonstrated in these cells. In addition, the low level of β-adrenergic receptors makes it unlikely that a low-efficiency interaction between the β-adrenergic receptor and Gi could have a significant effect.
Our data therefore are not consistent with a strict two-state model as described by Samama et al. (1). Our data do not, however, invalidate the classic ternary complex model (4), but they do suggest that receptor activation of Gs is more complicated than agonist stabilization of an active R* state. Indeed, other investigators have suggested that each agonist stabilizes a unique set of receptor conformations (known as the multistate model) with a range of G protein-activating abilities (17, 18). However, caution should be taken in extrapolating data from solubilized receptor systems and mutant receptors to biological systems. It should also be noted that although other models for receptor activation are possible, it may be that a more realistic model of agonist/receptor/G protein interactions might only complicate matters while only marginally improving the resolution of the two-state model. In any case, these findings have important implications not only in the study of the mechanisms involved in hormone-activated receptors but also for future studies involving the search pharmaceutical agents that may selectively activate different forms of a receptor.
Appendix
This section describes the derivation of a complete kinetic analysis of steady state adenylyl cyclase activity with respect to varying GTP concentrations in terms of the Cassel-Selinger cycle (19). The scheme shown below shows the cycle of Gs activation and inactivation after receptor activation. In this version, a shuttle model for G* s to adenylyl cyclase catalytic unit is assumed, although an analysis of the process in terms of a precoupled model for Gs and adenylyl cyclase would give almost identical values for the estimates of the rate of breakdown from experimental data. The description of active adenylyl cyclase is given as the solution of a quartic equation; however, the quartic itself can be factored into two quadratic equations, which simplifies matters somewhat. Equation 9In the scheme above, HR is the agonist-occupied receptor, Gs is the inactivated stimulatory G protein, G* s is activated G protein, C is inactive adenylyl cyclase, and G* sC is the active adenylyl cyclase complex. The kinetic constants governing the cycle are represented byk 1, the rate of heterotrimer (HRG) formation;k 10, the rate of dissociation of an inactive heterotrimer; k 2, the rate of active heterotrimer (HRG) dissociation; k 3, the rate of adenylyl cyclase activation by activated G* s, andk −1, the rate of G* s and G* sC inactivation via the intrinsic GTPase mechanism [because adenylyl cyclase does not seem to act as a GAP for G* s (20) k −1, the rates for G* s and G* sC are identical].
In the presence of GTP, we want to determine the steady state rate of adenylyl cyclase, G*
sC activity. From the conservation of mass, the levels of R, Gs, and C are defined as
Relating the rate of heterotrimer breakdown to they-intercept of the Eadie-Hofstee plot.
From eq. 27, we find that
The thermodynamics of receptor/G protein/agonist interactions in terms of the two-state model.
In this section, we use the two-state model to calculate the minimum fraction of receptors that must be in the R* state in the ternary complex to achieve the measured GTP shifts. If a receptor exists in two distinct conformational states that have different affinities for a ligand (Ka for state A and Kb for state B), then the measured dissociation constant (Kd ) for the ligand is given by Kd =Kafa + Kbfb , wherefa is the fraction of receptor in state A in the receptor/ligand complex, and fb is the fraction of receptor in state B in the ligand/receptor complex. Because this is a two-state model, fa plus fb is equal to unity. The above equations is derived from the first principles in Appendix .
If we adhere strictly to the two-state model, then when a ligand binds to a given receptor with increased affinity as a result of a change in conditions, according to the above equation, which is completely general, all changes in Kd must be interpreted in terms of changes in fa and fb, becauseKa and Kb refer to the binding to the two defined states.
In the absence of GTP, the β-adrenergic receptor forms a ternary complex with Gs and agonist (4). In this ternary complex, the affinity of the agonist for the receptor is increased relative to its binding in the presence of GTP when ternary complex does not accumulate. If B is the high affinity form of the receptor, thenKb < Ka , and the greater measured affinity for ligand in the absence of GTP requires there be a greater contribution of the B form to the ternary complex. In fact, the measured greater affinity places limits on the fraction of each conformation that is thermodynamically allowable.
If K
d1 is the dissociation constant of ligand from the receptor alone and K
d2 is the dissociation constant from the ternary complex, then we may writeK
d1 = K
afa1 +K
bfb1 and K
d2 =K
afa2 + K
bfb2. Noting that fa1 = 1 − fb1 andfa2 = 1 − fb2, these two equations may be combined and rearranged to give
It can be seen by inspection that any significant difference inK d1 and K d2 must lead to large proportion of receptors (>0.9) in the fb2 state. In the specific case of this report, even the rather weak agonist dobutamine requires that ≥80% (based on the differences in affinity obtained using a nonlinear two-component curve fit with the high and low values of the curve constant in addition to the low affinity portion of the curve) of the ternary complex be in the HR*G form (i.e., the active form) in terms of the two-state model.
Appendix
In the following system, only a single dissociation constant is observed provided the conversion time between the different state of the receptor R ↔ R* is rapid compared with the time allowed for equilibration of ligand binding.
Equation 44In this scheme, the dissociation for agonist bindingKD
is given as
Footnotes
- Received December 26, 1996.
- Accepted April 4, 1997.
-
Send reprint requests to: Roger Barber, Ph.D., The University of Texas-Houston Medical School, Department of Integrative Biology, Pharmacology, and Physiology, 6431 Fannin, P.O. Box 20708, Houston, TX 77225-0334. E-mail:rbarber{at}farmr1.med.uth.tmc.edu
-
This work was supported by National Institutes of Health Grant RR07710.
- The American Society for Pharmacology and Experimental Therapeutics