Complexity in pH-Dependent Ribozyme Kinetics - ACS Publications

Nov 20, 2017 - We previously published simulations of rate−pH profiles for ribozymes in terms of species plots for the general acid and general base...
0 downloads 0 Views 2MB Size
Perspective Cite This: Biochemistry XXXX, XXX, XXX−XXX

pubs.acs.org/biochemistry

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate−pH Profiles Erica A. Frankel†,‡ and Philip C. Bevilacqua*,†,‡,§ †

Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States Center for RNA Molecular Biology, Pennsylvania State University, University Park, Pennsylvania 16802, United States § Department of Biochemistry and Molecular Biology, Pennsylvania State University, University Park, Pennsylvania 16802, United States ‡

S Supporting Information *

ABSTRACT: Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate−pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as “influencers”. Herein, we simulate experimental rate−pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or “dark” pKa shifts, real pKa shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pKas of proton-transferring nucleobases and manifest themselves as “wavy” rate−pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate−pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

G

because like dark matter, their influence can be felt but they have not been observed.31 Second, recent evidence points to the pKa of the general acid or base being influenced by other ribozyme constituents, such as nearby cationic nucleobases or metal ions, acting as influencers.24 Such influencers lead to cooperative thermodynamic models with third-order or higherorder dimensionality, which can give “wavy” rate−pH profiles that can be accounted for by global fitting and statistics.24 We were interested in exploring the mechanistic bases behind ribozyme dark pKas and wavy rate−pH plots. We describe herein how the ribozyme partitions between various states, only a subset of which perform catalysis, i.e., are “functional”. Three kinetics schemes are described that build in complexity, incorporating cooperative interactions and then influencers. All three schemes employ what Cleland originally defined in protein enzymes as “reverse protonation”, namely protonation opposite to expectations for the free catalytic residues.32,33 For instance, A, C, and G are found to be catalytic

eneral acid−base catalysis is a pervasive mechanistic strategy in protein and RNA enzymes.1,2 In recent years, an increasing number of small ribozymes3−8 have been reported to use their nucleobases in catalysis, including the hairpin,9 HDV,10 hammerhead,11 glmS,12 twister,13 VS,14 pistol,15 twister sister,16 and hatchet ribozymes.17,18 Mechanistic roles for nucleobases include proton transfer, electrostatic catalysis, and hydrogen bonding with the transition state.19 The four nucleobases have pKa values far from neutrality and so are not optimal for general acid−base catalysis.20 Although the pKa values of some ribozyme bases have been observed experimentally to shift toward neutrality,21−24 much like in protein enzymes, including but not limited to glycosidases,25 cysteine proteinases,26 and alanine racemase,27,28other ribozyme residues either do not give observable shifts or have small apparent shifts.9,14,29 Such pKa shifts can also be apparent and manifest themselves in coupling RNA folding and activity.29,30 This work was motivated by two considerations of nucleobase pKas and catalysis. First, we wondered whether pKas might be shifted toward neutrality and contribute to rate acceleration but be hidden from traditional experimental approaches. We refer to such pKa shifts as “dark” pKa shifts © XXXX American Chemical Society

Received: August 15, 2017 Revised: November 20, 2017

A

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX

Perspective

Biochemistry in ribozymes in charged states of A+, C+, and G−, respectively, but are expected to be neutral at pH 7 for the free nucleotides. This principle applies to multiple protein enzymes, with more recent examples of a protein arginine deiminase34 and a βlactam synthetase.35 Evidence of reverse protonation in RNA enzymes is manifold: for instance, a cationic general acid of protonated C in the HDV ribozyme is supported by rate−pH profiles,10 inverse thiolate effects,36 Raman crystallography,21 and a crystal structure;37 a cationic general acid of protonated A in the hairpin ribozyme is supported by Raman crystallography;38 while an anionic general base of deprotonated G in the hammerhead and twister ribozymes is supported by rate−pH profiles and a crystal structure.8,11,39 General acid, general base, and influencer species can include not only nucleobases but also the sugar-phosphate backbone, metal ions, and cofactors. For this study, we focus on nucleobase catalysis; however, the same theoretical platform can be readily applied to bimolecular species.24

and base (Figure 1B and Excel simulation in the Supporting Information). This partition function can then be used to calculate the fraction of the functional forms of the general acid ( f HA+), general base ( f B−), and ribozyme [f R(1)] (see the Supporting Information). The fraction of fully functional ribozyme f R(1) dictates the rate−pH profile of ribozyme rate constant kobs, because kobs = f R(1)k1.2 Notably, f R(1) is simply equal to 1/Q. Because HA+ and B− do not interact in this case, the product of the curves for f HA+ (red) and f B− (blue) also gives the fraction of the functional ribozyme (black), as seen mathematically via f HA+ f B− = Q/Q2 = 1/Q = f R(1) and observed graphically (Figure 1C). In other words, the system is additive in free energy. This noncooperative case has a classic “bell-shaped” rate−pH profile (Figure 1C, black trace), which can be attained experimentally if the pKas of the bases are in the experimentally accessible pH range of approximately 5−10.29 In Figure 1C, we consider two scenarios: the pKas of the general acid and general base are 5 and 9, respectively, like in the ground state hairpin ribozyme,9,38 and the pKas of the general acid and general base are shifted toward neutrality at 6 and 8, respectively. As illustrated in Figure 1C, in each case, the fraction of the functional ribozyme is maximal at neutrality when plotted on the logarithmic axis, albeit with a broader maximum for the unshifted pKas. When the pKa values for the general acid and base are farthest from neutrality, at most only one in every 10000 molecules is functional. In contrast, when the general acid and general base pKa values shift one unit closer to neutrality, the effect is logarithmically additive and the functional population increases 100-fold. Additionally, the pKas of the general acid and base are visible: in the logarithmic plot at both the rounding points and extrapolation lines (see the Supporting Information) and in the linear plots at the inflection points. The significant effect that general acid and base pKa values have on the population of f R(1) underscores the importance of pKa shifting to enzyme activity.



CASE 1. SIMPLE NONCOOPERATIVE MECHANISM The first case is a two-dimensional noncooperative mechanism with two-nucleobase catalysis. This mechanism, which was described by our lab,2 is revisited briefly to frame the influence of cooperative interaction between nucleobases on ribozyme catalysis. In a noncooperative environment where the general acid and base do not interact, the pKa of the general acid remains constant regardless of the protonation state of the general base (top and bottom edges of Figure 1A), and vice versa. This system can be described at any pH using a partition function, Q, that requires only pKa values for the general acid



CASE 2. COOPERATIVE MECHANISM: EVIDENCE OF DARK pKas We build off of the preceding noncooperative two-nucleobase catalysis and incorporate an attractive, cooperative interaction between the cationic general acid and anionic general base. In this scenario, the corresponding pKas for those species are perturbed along two edges of the scheme (Figure 2A). A similar scenario has also been described using computational methods for protein enzymes, such as in RNase A.40 However, in that case, the general acid and base do not interact and there is no evidence of dark pKas. The pKa of the general acid is unshifted under conditions where the general base is in its neutral, nonfunctional protonated form (top edge of Figure 2A) but shifts upward toward neutrality by ΔpKcoop units, with a coupling free energy of −1.4ΔpKcoop (kcal/mol) under conditions where the general base is in its anionic functional form (bottom edge of Figure 2A).29 Likewise, the pKa of the general base is unshifted under conditions where the general acid is in its neutral nonfunctional deprotonated form (right edge of Figure 2A) but shifts downward toward neutrality by the same ΔpKcoop units under conditions under which the general acid is protonated (left edge of Figure 2A). The population of the functional form of the ribozyme, f R(1), at any pH can again be described using a partition function, although it is now dependent on both the unshifted and shifted pKa values of the general base and general acid (Figure 2B and

Figure 1. Simulations of pH−species plots for noncooperative general acid−base catalysis. (A) Scheme involving a general acid (GA) and general base (GB). (B) Partition function and fractions of GA, GB, and active ribozyme R(1). (C) Species plots of f HA+ with pKas of 5 and 6, f B− with pKas of 9 and 8, and f R(1), on logarithmic (left) and linear axes (right). B

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX

Perspective

Biochemistry

Figure 2. Simulations of pH−species plots for cooperative general acid−base catalysis. (A) Scheme involving a general acid (GA) and general base (GB). (B) Partition function and fractions of GA, GB, and active ribozyme R(1). (C) Species plots of f HA+ with pKas for the cooperative case of 5 and 7, f B− with pKas for the cooperative case of 9 and 7, and f R(1) on logarithmic (left) and linear axes (right). pKa values for the noncooperative case are 5 and 9 for the general acid and general base, respectively.

∼7 (Figure 2C left, solid red line), where at the general base starts to significantly populate its deprotonated form (Figure 2C left, solid blue line). At this point, the attraction between f HA+ and f B− is enhanced, leading to a slowing of the rate of deprotonation with pH and a rounding point in the general acid species plot near pH 9, the apparent pKa of the general base. Symmetric processes occur with the general base leading to a rounding point in the general base species plot near pH 5, the apparent pKa of the general acid. The net effect is an increase in the population of f R(1) (Figure 2C left, compare dashed and solid black lines) that is similar to that seen in case 1 wherein the pKas of the general acid and base were 6 and 8, respectively. See Figure S1 for further simulations and discussion of the coupling. Common ways to experimentally measure pKa values include spectroscopic methods such as NMR and Raman, which focus on f HA+ or f B− and typically require a ≥5% population for detection23 and are plotted on linear species plots (Figure 2C, right panel). Although the population of f R(1) is now 100-fold greater at any pH than under a noncooperative scenario, the maximum population of f R(1) still occupies only ∼1% of the total population (Figure 2C, left panel). Indeed, this nonadditive effect is completely hidden when looking at a linear species plot (Figure 2C, right panel). Standard spectroscopic methods are generally not sensitive enough to report these ratepromoting pKa shifts. Furthermore, kinetic methods, which report on f R(1) and are typically plotted on logarithmic species plots (Figure 2C, left panel) would also miss these ratepromoting pKa shifts; i.e., the black trace does not capture the

Excel simulation in the Supporting Information). Under these cooperative conditions, the product of the fraction of the general acid (solid red) and base (solid blue) no longer equals the fraction of the functional ribozyme (solid black). This is shown mathematically in Figure 2B and graphically in Figure 2C. In other words, the system is nonadditive in terms of free energy. In Figure 2C, we again plot the fraction of the functional form of the general acid and general base and the fully functional ribozyme along both logarithmic and linear scales. We consider a single scenario in which the unshifted pKas for the general acid and general base are 5 and 9, respectively, with a ΔpKcoop of 2 units (solid lines), and compare that scenario to the noncooperative case from the previous section (dashed lines). The magnitude of ΔpKcoop of 2 units is based on extrapolation of pKa shifts for model compounds (see the Supporting Information). It is visible from the logarithmic plots that even though the apparent pKa values for the general acid and base measured in an experimental rate-pH profile do not change from 5 and 9 (Figure 2C left, rounding points of the solid black line; or Figure 2C right, inflection points of the solid red and blue lines), consideration of the electrostatic attraction between the two species leads to a nonadditive catalytic effect on the fraction of the functional ribozyme: the population of the functional ribozyme increases by 100-fold versus the noncooperative case. It is insightful to consider the origin of this effect. In the cooperative case, the general acid with an apparent pKa of 5 is deprotonated linearly with increasing pH until the pH reaches C

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX

Perspective

Biochemistry

Figure 3. Simulations of pH−species plots for cubic cooperative general acid−base catalysis. (A) Scheme involving a general acid (GA), general base (GB), and influencer (I). (B) Partition function and fractions of GA, GB, and active ribozyme R(1) and R(5). (C) Species plots of f HA+ with a fixed pKa of 7, f B− with pKas of 9 and 7, influencer f HI+ (pink) with a pKas of 4 and 6, and f R(1)+R(5), which represents the two fully functional forms of the ribozyme, on logarithmic (top) and linear (bottom) axes. (D) Same as panel C but with a pKa of GA fixed at 11.

profiles, and thus, fractions of functional ribozyme plots are wavy, too. The coupling free energy between HI+ and B− is mathematically and graphically equivalent to that between HA+ and B− in case 2, with a coupling of ΔpKcoop of 2 units. Because the general acid does not couple in the case presented here, coupling between HI+ and B− is the same on the left and right sides of the cube, which is termed direct coupling.43,44 The effect of the influencer makes physical sense in terms of the ribozyme structure. The active sites of ribozymes have nucleobases, which may be charged, in the proximity of the catalytic species.24,45,46 It is thus likely that the ionization states of nucleobases near the active site may affect the ionization state of the general acid and general base. An example of this was recently reported in the hammerhead ribozyme where C17+ interacts with the G12− base24 and has also been seen in oligoalanine models (Table S1).47 Inspection of Figure 3 reveals that there is a near-plateau region for f R(1)+R(5), which represent the two states of the ribozyme with a fully functional general acid and general base, between pH ∼4 and 6 (top graph of Figure 3C). In the case of

wavy data present in the general acid and general base species plots. This leads to the notion that the rate-promoting pKa shifts in Figure 2 are dark pKa shifts and are indecipherable relative to noncooperative models (Figure S2). With this said, recent NMR relaxation−dispersion approaches are allowing users to see rare or “excited” states of RNA, especially where there is a significant chemical shift change between states, as might be expected for altered base protonation states.41,42 Thus, experimental elucidation of these rare pKas in active sites of ribozymes may be on the horizon.



CASE 3. NUCLEOBASE INFLUENCER: WAVY RATE−pH PROFILES To promote wavy behavior in f R(1), we need an attractive cooperative interaction between an influencer and either a general acid or a general base. In the case provided below, the anionic general base, B−, and cationic influencer, HI+, interact favorably, while for the sake of simplicity, the general acid and general base do not couple. In this way, the population of the general base but not the general acid is wavy in its rate−pH D

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX

Perspective

Biochemistry the general acid having a pKa of 7, bell-shaped behavior is observed (Figure 3C top, black line), while in the case of the general acid having a pKa of 11, the high-pH leg of the bell is missing in the experimentally accessible pH range of approximately 5−10 (Figure 3D top, black line). The width of the plateau can be increased by increasing the coupling (ΔpKa) between HI+ and B− (Excel simulations in the Supporting Information). General acid−base catalysis is a ubiquitous catalytic strategy in small ribozymes and protein enzymes alike, utilizing side chains as well as metal ions and cofactors. Depending on the species employed as the general acid and base, the observed rate−pH profiles can take on different shapes. Herein, we explored rate−pH profiles for three cases, motivated by simple pairwise interaction between the anionic general base and cationic general acid, or between the anionic general base and cationic influencer: a noncooperative thermodynamic square, a cooperative thermodynamic square, and a cooperative thermodynamic cube. Relative to the noncooperative case, cooperativity in the thermodynamic square (case 2) leads to rate enhancement without experimentally observable changes in rate−pH profiles for kinetics or signal−pH profiles for spectroscopy. The extent to which such “dark pKas”, which could also include pKas shifted away from neutrality, exist and promote ribozyme or protein enzyme reactivity awaits further exploration. Relative to the noncooperative case, cooperativity in the thermodynamic cube (case 3) leads to “wavy” rate−pH profiles in which the pKa of the functional form of the ribozyme, f R(1) + f R(5), is itself a function of pH. Studies of an enolase protein enzyme also revealed a pH-dependent pKa, with a line shape reminiscent of the black line in Figure 3C top panel, which those authors described as “wavelike”.33 The magnitude of the wave in their study parallels that seen in our study if we set the coupling to one pKa unit. Both this “wavelike” behavior and the “reverse protonation” mentioned above show that similar physical principles apply to enzymes made out of amino acids and nucleotides. A difference is that reverse protonation seems to be the exception in protein enzymes but the rule in RNA enzymes given that the pKas of the free nucleobases are far from neutrality. In case 3, we considered the scenario in which coupling between the influencer and general base is direct, but coupling could also be indirect and depend on the protonation state of the general acid, as well as additional influencers. Such additional complexity could lead to multiple plateaus and the appearance of fractional slopes in rate−pH profiles. Perturbations in the pKas of titratable groups that lead to complexities in the rate profiles such as these have also been supported on a theoretical level for protein enzymes.48−50 Global fitting is an approach that we recently introduced to fit rate−pH profiles that are parametric in Mg2+ concentration that have wavy features.24 Interactions among the general acid, general base, and influencer have been taken herein as electrostatic but could be of other forms that depend on charge state such as cation−π interactions.24,51,52 Calculations and new experimental methods may lend insight into the phenomena described herein.





Methods, derivations of equations, and supporting figures (PDF) Simulation spreadsheet (XLSX)

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Phone: (814) 863-3812. ORCID

Philip C. Bevilacqua: 0000-0001-8074-3434 Funding

NASA Exobiology, Grant 80NSSC17K0034. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank Jamie Bingaman, Kyle Messina, Drew Veenis, and Scott Showalter for helpful comments.



ABBREVIATIONS GA, general acid; GB, general base; f HA+, fraction of functional general acid; f B−, fraction of functional general base; f R(1), fraction of functional ribozyme.



REFERENCES

(1) Jencks, W. P. (1975) Binding energy, specificity, and enzymic catalysis: the circe effect. Adv. Enzymol. Relat. Areas Mol. Biol. 43, 219− 410. (2) Bevilacqua, P. C. (2003) Mechanistic considerations for general acid-base catalysis by RNA: revisiting the mechanism of the hairpin ribozyme. Biochemistry 42, 2259−2265. (3) Hertel, K. J., Herschlag, D., and Uhlenbeck, O. C. (1994) A kinetic and thermodynamic framework for the hammerhead ribozyme reaction. Biochemistry 33, 3374−3385. (4) Rupert, P. B., and Ferre-D’Amare, A. R. (2001) Crystal structure of a hairpin ribozyme-inhibitor complex with implications for catalysis. Nature 410, 780−786. (5) Nakano, S.-i., Proctor, D. J., and Bevilacqua, P. C. (2001) Mechanistic characterization of the HDV genomic ribozyme: assessing the catalytic and structural contributions of divalent metal ions within a multichannel reaction mechanism. Biochemistry 40, 12022−12038. (6) Wilson, T. J., McLeod, A. C., and Lilley, D. M. J. (2007) A guanine nucleobase important for catalysis by the VS ribozyme. EMBO J. 26, 2489−2500. (7) Klein, D. J., Been, M. D., and Ferre-D’Amare, A. R. (2007) Essential role of an active-site guanine in glmS ribozyme catalysis. J. Am. Chem. Soc. 129, 14858−14859. (8) Liu, Y., Wilson, T. J., McPhee, S. A., and Lilley, D. M. J. (2014) Crystal structure and mechanistic investigation of the twister ribozyme. Nat. Chem. Biol. 10, 739−744. (9) Kath-Schorr, S., Wilson, T. J., Li, N.-S., Lu, J., Piccirilli, J. A., and Lilley, D. M. J. (2012) General acid−base catalysis mediated by nucleobases in the hairpin ribozyme. J. Am. Chem. Soc. 134, 16717− 16724. (10) Nakano, S., Chadalavada, D. M., and Bevilacqua, P. C. (2000) General acid-base catalysis in the mechanism of a hepatitis delta virus ribozyme. Science 287, 1493−1497. (11) Han, J., and Burke, J. M. (2005) Model for general acid-base catalysis by the hammerhead ribozyme: pH-activity relationships of G8 and G12 variants at the putative active site. Biochemistry 44, 7864− 7870. (12) Klein, D. J., Been, M. D., and Ferré-D’Amaré, A. R. (2007) Essential role of an active-site guanine in glmS ribozyme catalysis. J. Am. Chem. Soc. 129, 14858−14859.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.biochem.7b00784. E

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX

Perspective

Biochemistry (13) Eiler, D., Wang, J., and Steitz, T. A. (2014) Structural basis for the fast self-cleavage reaction catalyzed by the twister ribozyme. Proc. Natl. Acad. Sci. U. S. A. 111, 13028−13033. (14) Smith, M. D., Mehdizadeh, R., Olive, J. E., and Collins, R. A. (2008) The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK (a). RNA 14, 1942−1949. (15) Harris, K. A., Lünse, C. E., Li, S., Brewer, K. I., and Breaker, R. R. (2015) Biochemical analysis of pistol self-cleaving ribozymes. RNA 21, 1852−1858. (16) Liu, Y., Wilson, T. J., and Lilley, D. M. J. (2017) The structure of a nucleolytic ribozyme that employs a catalytic metal ion. Nat. Chem. Biol. 13, 508−513. (17) Li, S., Lünse, C. E., Harris, K. A., and Breaker, R. R. (2015) Biochemical analysis of hatchet self-cleaving ribozymes. RNA 21, 1845−1851. (18) Lee, K.-Y., and Lee, B.-J. (2017) Structural and biochemical properties of novel self-cleaving ribozymes. Molecules 22, E678. (19) Lilley, D. M. J. (2011) Mechanisms of RNA catalysis. Philos. Trans. R. Soc., B 366, 2910−2917. (20) Izatt, R. M., Christensen, J. J., and Rytting, J. H. (1971) Sites and thermodynamic quantities associated with proton and metal ion interaction with ribonucleic acid, deoxyribonucleic acid, and their constituent bases, nucleosides, and and nucleotides. Chem. Rev. 71, 439−481. (21) Gong, B., Chen, J. H., Chase, E., Chadalavada, D. M., Yajima, R., Golden, B. L., Bevilacqua, P. C., and Carey, P. R. (2007) Direct measurement of a pK(a) near neutrality for the catalytic cytosine in the genomic HDV ribozyme using Raman crystallography. J. Am. Chem. Soc. 129, 13335−13342. (22) Gong, B., Klein, D. J., Ferré-D’Amaré, A. R., and Carey, P. R. (2011) The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate. J. Am. Chem. Soc. 133, 14188− 14191. (23) Thaplyal, P., and Bevilacqua, P. C. (2014) Experimental approaches for measuring pKa’s in RNA and DNA. Methods Enzymol. 549, 189−219. (24) Frankel, E. A., Strulson, C. A., Keating, C. D., and Bevilacqua, P. C. (2017) Cooperative interactions in the hammerhead ribozyme drive pKa shifting of G12 and its stacked base C17. Biochemistry 56, 2537− 2548. (25) McIntosh, L. P., Hand, G., Johnson, P. E., Joshi, M. D., Korner, M., Plesniak, L. A., Ziser, L., Wakarchuk, W. W., and Withers, S. G. (1996) The pKa of the general acid/base carboxyl group of a glycosidase cycles during catalysis: a 13C-NMR study of bacillus circulans xylanase. Biochemistry 35, 9958−9966. (26) Pinitglang, S., Watts, A. B., Patel, M., Reid, J. D., Noble, M. A., Gul, S., Bokth, A., Naeem, A., Patel, H., Thomas, E. W., Sreedharan, S. K., Verma, C., and Brocklehurst, K. (1997) A classical enzyme active center motif lacks catalytic competence until modulated electrostatically. Biochemistry 36, 9968−9982. (27) Harris, T. K., and Turner, G. J. (2002) Structural basis of perturbed pKa values of catalytic groups in enzyme active sites. IUBMB Life 53, 85−98. (28) Sun, S., and Toney, M. D. (1999) Evidence for a two-base mechanism involving tyrosine-265 from arginine-219 mutants of alanine racemase. Biochemistry 38, 4058−4065. (29) Moody, E. M., Lecomte, J. T., and Bevilacqua, P. C. (2005) Linkage between proton binding and folding in RNA: a thermodynamic framework and its experimental application for investigating pKa shifting. RNA 11, 157−172. (30) Knitt, D. S. (1996) pH dependencies of the Tetrahymena ribozyme reveal an unconventional origin of an apparent pKa. Biochemistry 35, 1560−1570. (31) Tyson, N. D. (2017) Astrophysics for people in a hurry, 1st ed., W. W. Norton & Co., London. (32) Cleland, W. W. (1977) Determining the chemical mechanisms of enzyme-catalyzed reactions by kinetic studies. Adv. Enzymol. Relat. Areas Mol. Biol. 45, 273−387.

(33) Sims, P. A., Larsen, T. M., Poyner, R. R., Cleland, W. W., and Reed, G. H. (2003) Reverse protonation is the key to general acid− base catalysis in enolase. Biochemistry 42, 8298−8306. (34) Knuckley, B., Bhatia, M., and Thompson, P. R. (2007) Protein arginine deiminase 4: evidence for a reverse protonation mechanism. Biochemistry 46, 6578−6587. (35) Raber, M. L., Arnett, S. O., and Townsend, C. A. (2009) A conserved tyrosyl-glutamyl catalytic dyad in evolutionarily linked enzymes: carbapenam synthetase and beta-lactam synthetase. Biochemistry 48, 4959−4971. (36) Das, S. R., and Piccirilli, J. A. (2005) General acid catalysis by the hepatitis delta virus ribozyme. Nat. Chem. Biol. 1, 45−52. (37) Chen, J.-H., Yajima, R., Chadalavada, D. M., Chase, E., Bevilacqua, P. C., and Golden, B. L. (2010) A 1.9 Å crystal structure of the HDV ribozyme precleavage suggests both Lewis acid and general acid mechanisms contribute to phosphodiester cleavage. Biochemistry 49, 6508−6518. (38) Liberman, J. A., Guo, M., Jenkins, J. L., Krucinska, J., Chen, Y., Carey, P. R., and Wedekind, J. E. (2012) A transition-state interaction shifts nucleobase ionization toward neutrality to facilitate small ribozyme catalysis. J. Am. Chem. Soc. 134, 16933−16936. (39) Martick, M., and Scott, W. G. (2006) Tertiary contacts distant from the active site prime a ribozyme for catalysis. Cell 126, 309−320. (40) Dissanayake, T., Swails, J. M., Harris, M. E., Roitberg, A. E., and York, D. M. (2015) Interpretation of pH−activity profiles for acid− base catalysis from molecular simulations. Biochemistry 54, 1307−1313. (41) Kimsey, I. J., Petzold, K., Sathyamoorthy, B., Stein, Z. W., and Al-Hashimi, H. M. (2015) Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes. Nature 519, 315−320. (42) Szymanski, E. S., Kimsey, I. J., and Al-Hashimi, H. M. (2017) Direct NMR evidence that transient tautomeric and anionic states in dG.dT Form Watson-Crick-like base pairs. J. Am. Chem. Soc. 139, 4326−4329. (43) Horovitz, A., and Fersht, A. R. (1990) Strategy for analysing the co-operativity of intramolecular interactions in peptides and proteins. J. Mol. Biol. 214, 613−617. (44) Di Cera, E. (1998) Site-specific thermodynamics: Understanding cooperativity in molecular recognition. Chem. Rev. 98, 1563− 1592. (45) Wilcox, J. L., and Bevilacqua, P. C. (2013) pKa shifting in double-stranded RNA is highly dependent upon nearest neighbors and bulge positioning. Biochemistry 52, 7470−7476. (46) Gottstein-Schmidtke, S. R., Duchardt-Ferner, E., Groher, F., Weigand, J. E., Gottstein, D., Suess, B., and Wohnert, J. (2014) Building a stable RNA U-turn with a protonated cytidine. RNA 20, 1163−1172. (47) Cantor, C. R., and Schimmel, P. R. (1980) Biophysical chemistry, W. H. Freeman, San Francisco. (48) Onufriev, A., Case, D. A., and Ullmann, G. M. (2001) A novel view of pH titration in biomolecules. Biochemistry 40, 3413−3419. (49) Ullmann, G. M. (2003) Relations between protonation constants and titration curves in polyprotic acids: A critical view. J. Phys. Chem. B 107, 1263−1271. (50) Klingen, A. R., Bombarda, E., and Ullmann, G. M. (2006) Theoretical investigation of the behavior of titratable groups in proteins. Photochem. Photobiol. Sci. 5, 588−596. (51) Reiss, C. W., Xiong, Y., and Strobel, S. A. (2017) Structural basis for ligand binding to the guanidine-I riboswitch. Structure 25, 195− 202. (52) Huang, L., Wang, J., and Lilley, D. M. J. (2017) The structure of the guanidine-II riboswitch,. Cell Chem. Biol. 24, 695−702.

F

DOI: 10.1021/acs.biochem.7b00784 Biochemistry XXXX, XXX, XXX−XXX