Structure of the PAP ternary complex
PAP(D154A) was crystallized with MgATP and a 5-mer poly(A) oligonucleotide. The mutation of Asp154, a catalytically essential active site residue, renders the enzyme inactive and was introduced to trap the closed, substrate-bound complex. This mutation resulted in only one of the catalytically-required metals being bound at the active site, but the mutation does not appear to significantly perturb the expected positions of the bound substrates. To our knowledge, this is the first crystal structure of a polymerase wherein one of the three catalytically-required aspartic acids has been mutated. There are, however, many instances where structures have been determined with only metal B at the active site, and in many of these the missing metal appears to cause a relatively small change in the nucleotide conformation. Our structure suggests that the strategy of mutating one of the catalytic aspartic acids may be generally useful, particularly in cases where, as in PAP, using unnatural substrates to trap the enzyme-substrate complex proves problematic. The structure was refined to 1.8 Å resolution (), and clear electron density was observed for both MgATP and RNA (Supplemental Figure 1
). In the ternary complex, the structure of individual domains are similar to those seen earlier, but these domains have undergone significant movement about both of the previously defined hinge regions; one between the N- and middle domains (residues around 40 and around 190) and the other between the middle and C-terminal domains (residues around 353) (Balbo et al., 2007
). As shown in , when compared to the most open structure of yeast PAP (molecule A of PDB code 1FA0) the N-terminal domain is rotated 23.4° towards the middle domain. The trajectory of this rotation is similar to that seen earlier, but this structure is 14.7° more closed than the most closed state reported previously (PDB code 2HHP). The C-terminal domain is in an intermediately closed state, near the average position of this domain observed in the various crystal structures. The translational components associated with movements about both hinges are very small.
Figure 1 Structure of PAP in complex with ATP and RNA. (A) The structure of the closed, ternary complex of PAP is shown in purple, and ATP and RNA are shown in white and yellow, respectively. The middle domains of the earlier PAP structures have all been superimposed (more ...)
In earlier PAP crystal structures, a large cleft (~30 Å long and ~35 Å deep) was observed between the N- and C-terminal domains (Balbo et al., 2007
; Bard et al., 2000
; Martin and Keller, 2004
; Martin et al., 2000
). In these structures, the active site is located at the bottom of this cleft, near the interface of the N- and middle domains. In the closed state, the N- and C-terminal domains interact, closing off the top of the cleft, but leaving openings at both ends near the bottom (). The closed structure reveals extensive contacts involving both substrates and residues within the central cavity. In fact, the substrates themselves largely mediate contact across the domains in the closed state. The enzyme-substrate interface buries ~2470 Å2
on the RNA and 1090 Å2
on the protein); interactions between the N- and C- terminal domains bury an additional ~360 Å2
of accessible surface area. The latter interface has good shape complementarity and includes a salt bridge between Arg125 in the N-terminal domain and Glu373 in the C-terminal domain, located near the top of the cleft.
The PAP-bound RNA substrate adopts an extended conformation, and there are no contacts between the RNA bases. In this respect, the PAP-RNA complex is quite different from the nucleic acid-containing complexes of two other template-independent polymerases, CCA adding enzyme and terminal deoxynucleotidyltransferase. All three enzymes belong to the polymerase beta family, but contain no homology outside of the catalytic palm domain. In both other structures, the active site cleft into which the nucleic acid substrate binds is wider to accommodate the single-stranded oligonucleotide chains, which are much more solvent exposed and generally form base-stacking interactions (Delarue et al., 2002
; Tomita et al., 2004
). This is in stark contrast to the extensive protein-RNA interactions observed in PAP ternary complex. The trajectory of the RNA leaving the active site also differs when PAP is compared to these other polymerases. PAP is also distinct from other ssRNA-binding proteins. Though the C-terminal domain of PAP contains an RNA recognition motif (RRM), the way in which this domain interacts with the RNA is completely dissimilar to that seen in other RRM proteins. RNA is generally bound to a particular face of RRM domains (Conte et al., 2000
), but in PAP, this RNA-binding surface is largely occluded by C-terminal domain residues 419 through 439. Furthermore, while ssRNA does interact with the C-terminal domain, it does so through residues which are on the opposite side of the RRM’s central beta sheet relative to the RNA-binding parts of the other structures.
Other than those of bovine PAP, the structures most similar to yeast PAP are those of the trypanosomal terminal uridylyltransferases, Tb
RET2 and Tb
TUT4 (Deng et al., 2005
; Stagno et al., 2007
). These proteins add a uradine residue to the 3′ end of a single-stranded RNA primer in a template-independent reaction, and they are homologous to PAP throughout the PAP’s N-terminal and middle domains. TUTases lack PAP’s C-terminal domain, but, like PAP, they are capable of significant domain closure. There is no crystal structure of a TUTase with the 3′ end position occupied. As discussed below, superimposition of PAP with TUTases reveals a very similar incoming nucleotide orientation, and as there is space in the TUTase structures for a nucleotide similar to that at position -1 in PAP, it would appear that the PAP structure is a good model for the TUTase ternary complex (Supplemental Figure 2
There are 24 protein-RNA interactions that are mediated by an extensive network of highly ordered water molecules, most of which are buried (as defined in Methods) within the interior of the complex ( and ). Overall, 58 water molecules are completely buried, and 31 of these are within the central, substrate-binding cavity. The degree to which the individual water molecules are buried is indicated in , and all of the solvent molecules are shown in Supplementary Figure 3
. These ordered waters are mostly hydrogen bonded to the polar atoms of the adenine bases as well as to the phosphodiester backbone. Consequently, the poly(A) substrate remains largely solvated in the enzyme ternary complex. Recognition of the 2′ and 3′ hydroxyl groups of the ribose of ATP is also mediated by bridging waters, which seems to be a strategy employed in many structurally unrelated ribonucleotide binding proteins (Babor et al., 2002
). As discussed below, some of the waters, particularly those that interact with the base of ATP, could help ensure substrate specificity.
Figure 2 Determinants of ATP and RNA binding. Interacting residues from the N-terminal, middle, and C-terminal domains are colored yellow, green, and cyan, respectively. Residues within cloud-shaped bubbles on either side of the bases indicate hydrophobic/van (more ...)
Figure 3 Stereo views of the substrate binding sites of PAP. (A–C) Detailed substrate interactions formed in closed, ternary complex. ATP (yellow carbons) and the 3′ end (blue carbons) are shown along with PAP with the surrounding amino acids (green (more ...)
Substrate binding determinants
In the sections below, the individual bases will be referred to by their position on the chain, with -1 being the 3′ end of the RNA and -5 denoting the 5′ end. The nucleotide at position -5 leaves the protein and interacts with F511 and K478 in a neighboring molecule. We do not believe this interaction is physiologically important. Poly(A) residue -4 is located at the surface of the enzyme, partially exposed to the solvent. The adenine bases at positions -1, -2, and -3 are completely buried within the closed structure, and PAP interacts with the 2′-OH groups of each of these nucleotides. These observations are consistent with earlier studies indicating that DNA primers are poor substrates for the enzyme (Edmonds, 1990
) and that efficient polyadenylation requires a primer of at least three residues (Zhelkovsky et al., 1998
PAP does not exhibit great specificity with respect to the sequence of its RNA substrate (Butler and Platt, 1988
; Butler et al., 1990
). This is consistent with the requirement for initiating polyadenylation at non-ploy(A) sequences that can occur in freshly cleaved pre-mRNA substrates (Zhao et al., 1999
). There are, however, three base-specific interactions with the RNA observed in the crystal structure. Two of these are between the protein and the RNA base at position -4, where the side chains of Lys 392 and Glu487 make hydrogen bonds with the N1 and N6 atoms of the base, respectively. The -4 nucleotide lies on the surface of the C-terminal domain and exhibits relatively high B-factors, suggesting that it may not be as well anchored to the protein as the bases at positions -1, -2 and -3. In addition to the two interactions mentioned above, the hydrophobic face of the -4 base sits atop C-terminal domain residues Leu388 and Leu 491, and the 2′ hydroxyl makes a hydrogen bond with the side chain of His314. As the base-specific interactions we observe at position -4 are at the protein surface, there is likely little energetic penalty associated with their loss when other nucleotides occupy this position because hydrogen bonds to the base can easily be replaced by those to water.
The nucleotide at position -3 interacts with all three domains of the polymerase. The side chain of Gln294 makes a hydrogen bond with the N6 position of the adenine ring, potentially providing base specificity at this position; however, this residue is not conserved. The -3 base is sandwiched between highly conserved residue Phe140, from the N-terminal domain, and His314, from the middle domain. The side chains of both of these residues were exposed to solvent in earlier PAP structures and appeared highly flexible as indicated by their B-factors and in some cases the absence of side chain electron density. Conserved residue Arg387 from the C-terminal domain also interacts with the -3 phosphate moiety, as does His314. The latter residue is replaced with tyrosine in higher eukaryotes, thus presumably retaining both its base stacking and hydrogen bonding functions. Finally, the 2′ hydroxyl makes a hydrogen bond with the side chain of Asn315. As is typical for the remaining nucleotides in this structure, this nucleoside makes six water-mediated interactions. Many of these interactions, as well as the protein atoms to which these water molecules are coordinated, are detailed in and .
The base of the -2 nucleotide is sandwiched between the N- and C-terminal domains. One side of its hydrophobic surface interacts with Val136; the other contacts Gly380 and Ser384 from helix N. The 2′ hydroxyl of this nucleotide makes a hydrogen bond with Lys145, and an additional hydrogen bond is made between the phosphate moiety and the main chain amide of Phe140. All of these residues are very highly conserved, and again there are well ordered, buried water molecules which make additional interactions with the nucleotide.
The base of the 3′-terminal residue (-1 position) of the poly(A) substrate is sandwiched between the side chain of Val141 and the base of the incoming ATP. In addition, hydrogen bonds are formed between the 2′ hydroxyl and the side chain hydroxyl of Tyr87 and between the backbone phosphate oxygen and the δN of Asn226. In the two-metal ion dependent nucleotidyltransferases, the 3′-terminal residue is also bound through coordination of the 3′ hydroxyl group to a second, so-called catalytic Mg2+
ion (Batra et al., 2006
; Sawaya et al., 1994
; Sawaya et al., 1997
). The mutation of Asp154 abrogates binding of this metal, and this interaction is not observed in the present structure. Nonetheless, the positions of the MgATP and the poly(A) terminus observed in the structure generally conform to this mechanism. The position of the -1 nucleotide coincides well with the nucleotide (3′-dAMP) bound at the poly(A) subsite in the earlier structure of yeast PAP (Bard et al., 2000
) () and with the nucleotides observed in structures of DNA pol β (Sawaya et al., 1997
) (). Furthermore, despite the absence of the second Mg2+
, the 3′ hydroxyl is located 3.2 Å from the α-phosphorus atom of ATP, as if poised for nucleophilic attack. This all supports the overall correctness of the nucleotide conformation in this structure.
Figure 4 Comparison of the active site with other structures. (A) The alpha carbons of the N-terminal domain of the RNA-bound PAP molecule were superimposed with those of the two crystallographically-independent molecules from PDB code 1FA0. In the earlier structure (more ...)
The triphosphate moiety of ATP exhibits tridentate coordination to metal B, and this Mg2+
exhibits octahedral geometry involving coordination to the side chain carboxylates of Asp100 and Asp102 and an ordered water, in addition to the triphosphate (). This structural arrangement is an essential characteristic of this class of enzymes (Sawaya et al., 1997
). Additionally, the γ-phosphate of ATP makes direct contact with the side chains of Ser89, Lys215, and Tyr224. The latter two interactions are formed only in the closed enzyme state. A backbone amide from Ser89 forms an additional hydrogen bond to the non-bridging oxygen of the β-phosphate. The non-bridging α-phosphate oxygen is exposed to solvent molecules and does not interact with active site residues, a result consistent with previous studies (Balbo et al., 2007
). The coordination environment about this oxygen is important because a partial negative charge develops here in the transition state (Burgers and Eckstein, 1979
; Steitz and Steitz, 1993
). Exposure of this atom to solvent and the absence of an appropriately positioned acidic residue is evidence that water-derived protons function in charge stabilization during catalysis.
Structural basis for nucleotide base recognition
A central question regarding the mechanism of PAP concerns the structural determinants for ATP specificity. The present structure confirms that domain closure results in the formation of new interactions, both direct and water-mediated, between the enzyme and the nucleotide substrate. It also suggests that there are multiple factors determining nucleotide base recognition. Overall, the nucleotide binding site exhibits good surface complementarity to the MgATP, but there are no direct hydrogen bonds observed between the enzyme and the adenine base. The surface of the binding site of the adenine-ribose moiety in PAP is defined primarily by two features: (1) non-polar interactions above and below the plane of the adenine base and to the sugar, and (2) water mediated interactions to polar groups in the plane of the base and to the sugar hydroxyl groups. These results seem to conform to general principles utilized in many, structurally unrelated adenylate-binding proteins (Moodie et al., 1996
). There is also suggestive evidence that N226 may participate in adenine recognition via an interaction with the N6 group.
As shown in , a base stacking interaction occurs between the ATP and the 3′-terminus of the poly(A) substrate. The side chain of Val234 makes van der Waals contact with the ribose-base moiety of ATP, on the opposite face as the base-stacking interaction. Nucleotide bases form stacking interactions both in solution and in nucleic acids (Broom et al., 1967
; Ts’o et al., 1969
; Ts’o et al., 1962
), with the order of stacking preference: purine-purine > purine-pyrimidine > pyrimidine-pyrimidine. The theoretical free energy change realized upon purine-purine base stacking is 2–6 kcal/mol, and A-A stacking is more favorable than A–C stacking by 0.5 – 2 kcal/mol (Friedman and Honig, 1995
; Norberg and Nilsson, 1995
). As discussed below, these energetic values are consistent with the lower catalytic efficiency exhibited by PAP for nucleotidyltransfer reactions utilizing the pyrimidine substrate, CTP (Balbo et al., 2005
). In vivo
, the terminal nucleotidyl residue formed by cleavage at the poly(A) site is often, but not always, adenylate (Zhao et al., 1999
). Previously, we studied the role of base stacking in substrate binding and the effect of alternative 3′-terminal bases by measuring the kinetics of the reaction: A17
-C + MgATP
-CA + MgPPi
(Balbo et al., 2005
). The catalytic efficiency for this reaction was 10-fold lower than the reaction with A18
and MgATP, supporting the assertion that base stacking is significant.
Base stacking might provide a substantial amount of substrate binding free energy, but it does not explain specificity between ATP and GTP. PAP exhibits 800-fold lower catalytic efficiency when GTP is utilized as a substrate instead of ATP (this work), strongly suggesting that other factors are involved in base recognition. Discrimination between adenine and guanine is largely based on differential electrostatic properties and shape (Moodie et al., 1996
; Nobeli et al., 2001
). In the present structure, three water molecules, which are completely buried within the interior of the complex, coordinate the N3, N6, and N7 positions of the incoming nucleotides base. Substitution of the N6 of adenine with the O6 of guanine replaces a H-bond donor with a H-bond acceptor, and results in protonation of the N1 atom, altering its H-bonding properties as well. These perturbations would be expected to have significant consequences on base recognition even if molecular recognition at any of these positions is mediated by a buried water. Also, the adenine near the C2 position is packed in a small pocket formed by the side chains of Thr304, Met310, and Ala312. Therefore, steric exclusion might contribute to negative selection against the 2-amino group of GTP.
Finally, the proximity of Asn226 to the adenine moiety of ATP has led to previous speculation of its involvement in base recognition (Balbo et al., 2007
; Bard et al., 2000
). In the present structure, the δO of Asn226 is only 3.57 Å from the N6 position of the ATP. A small change in the nucleotide or side chain conformation, perhaps by restoration of the Asp154 side chain and the second Mg2+
, could align these groups and promote an interaction. Indeed, no other active site residue is close enough to the nucleotide to make a base-specific contact without a large structural rearrangement of the active site. To investigate the role of Ans226 in nucleotide specificity, a kinetic analysis of GTP utilization by PAP (w.t.) and the N226A mutant was performed (). As is the case for CTP (Balbo et al., 2005
), substrate selectivity against GTP is manifested in Vmax
, rather than Km
. reports the apparent ΔΔGcat
contributed by each mutated residue as measured for a given reaction. The small discrepancy between the ΔΔGcat
for the reactions utilizing ATP and GTP (−1.29 − (−0.74) = −0.54 kcal/mol) is consistent with a small degree of antagonism between N226 and the guanine base, as would be expected if this residue is involved in the negative selection of GTP. However, due to the limitations in interpreting kinetic studies of this nature, such speculation should be made cautiously. Presently, stronger evidence indicating an ATP-Asn226 interaction is unavailable.
Steady State Kinetic Parameters for PAP and Mutants
Figure 5 Kinetics and mutangenesis of PAP. (A) Total apparent free energy contributions to catalysis (Δ Gcat) from mutated residues. Values are derived from nucleotidyltransfer (utilizing ATP, CTP or GTP) or pyrophosphorolysis (PPi) reactions. (B) Relative (more ...)
Kinetics of PAP mutants
Kinetic studies were undertaken to investigate the relationship between substrate binding, domain movement, and catalysis. We describe below a model wherein the enzyme converts between open and closed conformations during the catalytic cycle residues in the central cavity function to stabilize the closed conformation of the enzyme-substrate ternary complex. Full active site assembly occurs only upon formation of the closed state, thus proper substrate recognition, which induces formation of the closed state, will promote catalysis. This results in nucleotide specificity appearing in the Vmax
term, rather than Km
. Previously, yeast and bovine PAP have both been the subject of mutagenesis/activity studies (Martin et al., 1999
; Martin et al., 2004
; Zhelkovsky et al., 2004
), however these studies did not consider such a model. Rather than attempt to reinterpret the existing data, which is complicated by the use of different reaction conditions, different methods of assaying activity, and intrinsic differences between the yeast and mammalian enzymes, we have chosen to perform a highly detailed analysis of a subset of mutations.
Here, the mutants N189A, K215A, Y224F, and N226A were characterized; these residues are all located in the active site cleft and are or have been implicated in substrate binding and recognition. Asn189 was previously suggested to be important in ATP specificity (Martin and Keller, 2004
). It is located at the hinge between the N-terminal and middle-domains, and while it does not interact with either substrate, it does make a tertiary contact to the carbonyl oxygen of Tyr307. This bridging interaction between the N- and middle-domains occurs only in the closed state. Lys215 and Tyr224 bind the γ-phosphate of ATP in the closed state and were observed to interact with the product, MgPPi
, in the open state (Bard et al., 2000
). Asn226 is in the middle domain. It does not interact with substrates in the open state, but in the closed state, it contacts a backbone phosphate of the poly(A) and possibly also the ATP base (see discussion above).
The steady state kinetic mechanism of PAP is rapid equilibrium random, and the overall rate is limited by the chemical step (Balbo et al., 2005
). The steady state kinetic parameters (Eqs. 1
) were evaluated for these mutant enzymes in both the forward (polyadenylation) and reverse (pyrophosphorolysis) directions and for the cytidylyltransfer reaction utilizing the alternative nucleotide substrate, CTP. The results are reported in alongside the published results (Balbo et al., 2005
) for wild type PAP. For each mutant, difference free energy changes relative to wild-type PAP were calculated for each of the reactions performed () including nucleotidyltransfer utilizing ATP, CTP or GTP and pyrophosphorolysis. These are defined as ΔΔGcat
)]. These values measure the apparent, total
free energy contribution of the residue to both substrate binding and catalysis. Agreement among the values for a given mutation reflects consistency between the different sets of experiments. Additionally, for each experiment involving the nucleotidyltransfer reactions with ATP or CTP, the values, ΔΔGES
(see Methods), were calculated and plotted against the total difference in free energy change (ΔΔGcat
) relative to a reference reaction (the cytidylyltransfer reaction catalyzed by K215A) (). This treatment of the data illustrates the relationship between the total free energy change for catalysis, ΔΔGcat
, and the individual energy contributions to this term from substrate binding (ΔΔGES
) and the chemical (rate determining) step (ΔΔG‡
Certain generalizations can be made upon inspection of the data in and . With regard to adenylyltransfer activity, all mutations primarily exert their effect on the “Km
” terms, Ka
. For the rapid equilibrium random mechanism employed by PAP, these terms reflect the equilibrium dissociation constants of the poly(A) (Ka
) and ATP (Kb
) substrates from the enzyme-substrate ternary
complex (see (Balbo et al., 2005
) for a thorough discussion). In contrast, the inhibition constants, Kia
, reflect dissociation of poly(A) and ATP, respectively, from the binary
complex. Importantly, the observation that a given mutation (see N226A, Y224F) has an essentially equal effect on both Ka
, but essentially no effect on the Ki
terms, indicates that these mutations do something other than simply disrupt substrate “binding.” We argue below that these mutations affect the equilibrium between the open- and closed-domain forms of the enzyme. Even N189A, a residue that does not directly contact substrates, but bridges the N- and middle domains in the closed state, exhibits a Km
effect. Furthermore, these mutations have a comparatively small effect on V1
, suggesting that the catalytic machinery was uncompromised. This result is consistent with the location of these residues, which are at the active site, but away from the chemically reactive parts of the substrate(s).
In contrast, for the pyrophosphorolysis reaction, () the mutations primarily affect the maximal velocity (V2). For N189A and N226A, the Ki and Km terms were very similar to those of wild type. Interestingly, the Y224F (and K215A) mutant(s) showed lower affinity for MgPPi in the binary complex (Kip), but the total free energy change (calculated from the product K′aKip) upon enzyme-product ternary complex formation by Y224F is similar to that of wild type (ΔΔG = 0.4 kcal/mol) due to a compensatory increase in affinity for An in the ternary complex. Attempts to fit the pyrophosphorolysis data for K215A were unsuccessful; they did suggest, however, that the Km values were much greater than 10 mM (not shown).
The cytidylyltransfer data exhibit more scatter than those for the reactions involving ATP. Even so, the total free energy change upon mutation calculated from these data is mostly consistent with those from the other experiments (). The largest discrepancy (0.5–0.6 kcal/mol) between these and the other data sets involves K215A and Y224F and suggests that the compound effect of mutating either K215 or Y224 and utilization of the CTP is not additive. This could possibly reflect the binding of CTP in an alternative (i.e., non-ATP-like) binding mode. Aside from this complication, the important generalization from this data set, particularly the experiments with wild type, N189A and N226A, is that the ΔΔGcat
is primarily accounted for in the maximal velocity rather than substrate binding (), a qualitatively different result from the adenylyltransfer data set. As also illustrates, when CTP is utilized by w.t. PAP, the apparent “binding” energy is similar to that when ATP is utilized, and the main difference is in the reaction velocity. This was previously noted (Balbo et al., 2005
) and indicates differences in the mechanisms for utilization of ATP and CTP. As mentioned above, the N226A mutant also exhibits a velocity effect for the guanylyltransfer reaction, mirroring the result for the CTP reaction. A general mechanism accounting for this phenomenon is discussed in detail below.
Mechanism of poly(A) polymerase
The kinetic results can be rationalized by the mechanism described in Scheme 1, where A
= poly(A), B
= MgATP, P
, and E
= PAP; the designation, E
′, refers to the enzyme when the poly(A) substrate is bound in the product-binding mode. This mechanism amends our previous model (Balbo et al., 2005
) by introducing conformational change steps (domain reorientation) on either side of the chemical step. The equilibrium constants describing the open-to-closed conformational change for the enzyme-substrate and enzyme-product ternary complexes are Kc1
, respectively, where each Kc
is the ratio of closed enzyme to open enzyme. The central step in which the chemical transformation of substrates occurs involves closed complexes and is denoted by the dashed box. In the sections below, we describe the predicted kinetic behavior of the system in Scheme 1 and explain how the mutagenesis results are consistent with such a model.
Because the steady state kinetic mechanism of PAP is rapid equilibrium random (Balbo et al., 2005
), relatively simple relationships exist between the kinetic parameters and the thermodynamic equilibrium constants in Scheme 1. It follows from both the previous kinetic study and the present work that domain movement is also rapid. The kinetic effect of internal conformational change steps such as those proposed here was previously described in the Tetrahymena
ribozyme system, and these references (Narlikar and Herschlag, 1998
; Narlikar et al., 1997
) should be consulted for a more detailed description of the kinetic features of this type of system. Briefly, the conformational change steps (Kc1
in Scheme 1) will affect the “Km
” terms for each substrate or product. In each case, the Km
term equals the equilibrium constant (or function of constants) describing dissociation of the substrate or product from the appropriate enzyme ternary complex in the open state
divided by the term, (1 + Kc
). This termcauses a reduction in Km
is significantly greater than 1 (closed conformation is favored), and has a negligible effect when Kc
is significantly less than 1 (open conformation is favored). Definitions for all kinetic parameters are given in .
Steady State Kinetic Parameters of PAP
The adenylyltransfer reaction catalyzed by wild type PAP can be described by the free energy diagram in . For reactions having a single, rate limiting step, the Michaelis complex corresponds to the lowest-energy ground state enzyme substrate complex that exists prior to the highest-energy transition state barrier (i.e. the rate limiting step). In this model, the closed enzyme conformation is of lower energy than the open conformation, and the Michaelis complex corresponds to EAB(closed)
. This idea is supported by previous fluorescence quenching results suggesting that MgATP and poly(A), specifically, induce domain closure (Balbo et al., 2007
). The free energy required to stabilize the closed, active state is derived from the numerous and extensive interactions formed in the closed state (tertiary interactions) as a result of binding both ATP and poly(A). Residues that contribute tertiary interactions will lower the free energy of both the closed ground state and the transition state by exactly the same amount, since these interactions occur in both states. This phenomenon has been called “uniform binding” (Albery and Knowles, 1976
). As shown in , the effect of mutating residues that donate tertiary interactions is to destabilize the closed state, and, hence, lower Kc1
and increase Km
. In this situation, the lowest energy ground state remains EAB(closed)
and the height of the energy barrier from EAB(closed)
to the transition state is unchanged, thus Vmax
is unaffected by the mutation. The mutagenesis results for the adenylyltransfer reaction () are consistent with this mechanism. Furthermore, as is diagramed in , the mutation does not affect the free energy of the enzyme product ternary complex (E
), thus the free energy barrier from E
to the transition state increases and the mutants exhibit a velocity effect in the reverse direction. Thus, the pyrophosphorolysis results are also consistent with this mechanism. One interpretation of these results is that the Michaelis complex for the reverse reaction corresponds to the open enzyme conformation, E
. The situation diagramed in represents a different mechanism in which the lowest-energy ground state (the Michaelis complex) is EAB(open)
. This occurs when there is insufficient free energy available from tertiary interactions in the closed state (i.e., Kc1
< 1). The data for the reactions utilizing either CTP or GTP are consistent with the mechanism of . Substitution of the nucleotide base from adenine to cytidine results in a loss of 3.4 kcal/mol (ΔΔGcat
) of available free energy (Balbo et al., 2005
); and from adenine to guanine, 4.0 kcal/mol (this work). As a result, the open conformation is favored significantly over the closed conformation. Because the Michaelis complex for this reaction corresponds to the open conformation, the mutants, particularly N226A and K215A, which affect the energy of the closed state, exhibit a velocity effect, as described by the model in .
The main catalytic features of PAP as determined from kinetic experiments, including the site directed mutagenesis (this work) and previous work with alternative substrates (Balbo et al., 2005
), are summarized as free energy diagrams in . In , the results for three “mutations” in the forward reaction are compared to the adenylyltransfer reaction catalyzed by wild type PAP: Y224F, K215A, and cytidylyltransfer catalyzed by wild type PAP (effectively, a “mutation” of the adenine base of ATP to cytidine). shows only the diagram for wild type adenylyltransfer reaction with the relevant free energy terms indicated. The terms ΔGcat
, and ΔG‡
are related to the steady state kinetic parameters V/Kia
, and V, respectively (see Methods). Also indicated are terms relating to the intrinsic binding energy, ΔGIB
, and strain energy, ΔGstrain
, where strain represents the energetic penalties incurred upon formation of the ground state. The relationship between intrinsic binding energy and strain has been widely recognized (Fersht, 1974
; Herschlag, 1988
; Jencks, 1975
), and can be represented by: ΔGES
. For poly(A) polymerase, ΔGstrain
represents the free energy utilized in the ground state destabilization mechanism. We have demonstrated that specificity for ATP (versus CTP or GTP) is manifested primarily in the Vmax
-fold effect), whereas only a small (3-fold) effect on Km
is observed. We proposed an induced fit mechanism for nucleotide specificity that was described as a composite of two effects: uniform binding and ground state destabilization (Balbo et al., 2005
). The model of Scheme 1 offers an elaboration of this mechanism in structural terms, specifically that ground state destabilization occurs in the closed conformation. In contrast, the cytidylyl- and guanylyltransfer reactions do not exhibit ground state destabilization (Balbo et al., 2005
); this is explained in part by perturbation of Kc1
so that the open complex is favored in the ground state, thus the energetic penalty from ground state destabilization is not incurred in the formation of the Michaelis complex when alternative NTPs are utilized. The polyadenylation results for N189A, N226A, and Y224F, suggest these are classic uniform binding mutations (Albery and Knowles, 1976
). The ΔΔGcat
for these mutants are ≤ 2 kcal/mol, thus leaving sufficient free energy available to induce the full ground state destabilization effect. In contrast, the K215A mutant, though generally similar to Y224F (), exhibits a measurable effect on ΔG‡
(). We speculate this is due to weaker induction of strain in the Michaelis complex, and could suggest that Lys215 plays a catalytic role, possible by donation of a positive charge in the ground state.
The mechanism for poly(A) polymerase must account for a number of experimental observations including the previously demonstrated rapid equilibrium random mechanism, domain motion, nucleotide specificity, and the kinetic effects of PAP mutation. According to the induced fit mechanism of PAP, the binding of both MgATP and poly(A) induce a conformational change to form the closed enzyme state. Analogous domain movements have been demonstrated in other polymerases (Doublie et al., 1999
), and indeed, rigid body motion such as domain reorientation or flap closure is a common feature in enzymes. In many cases, domain movement is associated with the full assembly of the active site and conversion of the enzyme to a catalytically active state (Hammes, 2002
). Numerous tertiary interactions between residues in the central cleft and substrates, as revealed in the present structure, provide the free energy necessary to stabilize the enzyme in the closed state. Furthermore, a substantial portion of the intrinsic binding energy derived from ATP base recognition is utilized to promote the reaction velocity via a ground state destabilization mechanism. This is the basis for manifestation of substrate specificity in the velocity term (Vmax), rather than the Km
. Upon substrate binding and formation of the E
ternary complex, domain closure is induced and adenylyltransfer is then catalyzed in the closed enzyme state. Because incorrect nucleotide substrates do not induce the conformational change, their utilization is inefficient, and substrate dissociation (of the incorrect nucleotide) is likely to occur. The details of how catalysis is promoted in the closed, fully assembled active site, including the chemical nature of the proposed ground state destabilization mechanism, remains an important question in the enzymology of PAP and all polymerase beta family enzymes.
PAP functions in vivo
as part of a multi-protein complex which includes a number of RNA-binding proteins (Colgan and Manley, 1997
; Zhao et al., 1999
). These proteins tether the elongating poly(A) substrate to PAP, the effect of which is to increase the local concentration of RNA at the site of polyadenylation. PAP is distributive (Lingner et al., 1991
; Wahle, 1991
), but appears processive in the presence of other polyadenylation factors (Preker et al., 1997
). The kinetic mechanism provides an explanation for how the intrinsic properties of the polymerase allow for efficient catalysis in the context of RNA tethering in vivo
. Efficient polyadenylation is governed, in large part, by the thermodynamics of substrate/product binding and domain movement. Poly(A) can bind in either of two, mutually exclusive binding modes (as either a substrate or product), and it is essential that the product poly(A) is released efficiently after each round of catalysis. The high values for K′a
indicate that product dissociation readily occurs even under the cellular condition of locally high poly(A) concentration. One new insight provided by this study (particularly the kinetics of pyrophosphorolysis by the mutants) is that the open conformation of the enzyme-product ternary complex (E
) is more stable than the closed conformation. This suggests that upon conversion of the enzyme-bound substrates to products, domain opening is induced, facilitating dissociation of products from the enzyme after each round of catalysis. Furthermore, we speculate that the large number of buried waters in the closed complex serves to lower the kinetic barrier to domain opening and closing, which must occur rapidly during successive rounds of adenylyltransfer, by minimizing the extent of substrate desolvation upon Michaelis complex formation. These factors together prevent unproductive reaction in the reverse direction and ensure efficient catalysis and forward flux under cellular conditions.