Constrained Bonding Environment in the Michaelis Complex of

Aug 7, 2012 - Constrained Bonding Environment in the Michaelis Complex of Trypanosoma cruzi Uridine ... Inverse 1′-3H and 5′-3H BIE's indicate a c...
0 downloads 0 Views 523KB Size
Rapid Report pubs.acs.org/biochemistry

Constrained Bonding Environment in the Michaelis Complex of Trypanosoma cruzi Uridine Phosphorylase Rafael G. Silva, D. Randal Kipp, and Vern L. Schramm* Department of Biochemistry, Albert Einstein College of Medicine of Yeshiva University, 1300 Morris Park Avenue, Bronx, New York 10461, United States S Supporting Information *

and ribose 1-phosphate,12 and KIE studies using phosphate or arsenate as the nucleophile indicated that the reaction proceeds via an expanded SN2-like transition state.13,14 If the C1′−N1 bond of uridine becomes longer in the Michaelis complex toward the transition state, the anomeric carbon develops partial sp2 character with an increased degree of freedom in the C1′−H1′ bond out-of-plane bending mode. Experimentally, this gives rise to a normal 1′-3H BIE. Scheme S1 of the Supporting Information depicts the atom numbering of uridine. Tritium BIE’s were measured for the 1′-H and 5′-H2 positions of uridine, using [5′-14C]uridine as the light substrate in competitive radiolabel experiments. Briefly, 20 μM TcUP was incubated with 50 μM uridine ([5′-14C]uridine and either [1′-3H]- or [5′-3H2]uridine) and 50 mM (NH4)2SO4 (pH 7.5) in 100 mM HEPES (pH 7.5). The incubation mixture (100 μL) was loaded onto an ultrafiltration apparatus, and an argon pressure of 30 psi was applied until approximately half the volume had passed through a semipermeable membrane (10 kDa molecular mass exclusion). At that time, 30 μL was removed from the top and bottom of the ultrafiltration wells. Scintillation counting and spectral deconvolution were performed as previously reported.13 BIE’s were calculated with eq 1, where 14CT and 14CB are the total 14C counts in the top and bottom wells, respectively, and 3HT and 3HB are the total 3H counts in the top and bottom wells, respectively.15 BIE’s are expressed as averages and their standard errors for at least 10 replicates in two independent experiments. The effect of sulfate ion on the 1′-3H BIE was assessed by excluding (NH4)2SO4 from the incubation mixture. Enzyme was omitted in control experiments.

ABSTRACT: The transition state for the Trypanosoma cruzi uridine phosphorylase (TcUP) reaction has an expanded SN2 character. We used binding isotope effects (BIE's) to probe uridine distortion in the complex with TcUP and sulfate to mimic the Michaelis complex. Inverse 1′-3H and 5′-3H BIE's indicate a constrained bonding environment of these groups in the complex. Quantum chemical modeling identified a uridine conformer whose calculated BIE's match the experimental values. This conformer differs in sugar pucker and uracil orientation from the unbound conformer and the transition-state structure. These results support ground-state stabilization in the Michaelis complex.

G

round-state destabilization has been proposed as a strategy employed by enzymes, including triosephosphate isomerase, enoyl-CoA hydratase, and alkaline phosphatase, to accomplish catalysis.1−4 The process consists of enzymeinduced electronic strain of substrates in the Michaelis complex toward their transition-state configurations. However, some theoretical and experimental evidence has questioned this hypothesis in certain enzymes, for instance, orotidine 5′decarboxylase and lysozyme.5−7 Kinetic isotope effects (KIE’s) measured by the internal competition method in enzymatic reactions report on V/K, probing events from free substrate in solution up to, and including, the first irreversible step, often assumed to be the release of the first product under initial rate conditions. KIE’s on V/K can be corrected for kinetic complexities, such as commitment factors, to yield intrinsic KIE’s, which reflect differences in bond vibration frequencies between the free substrate and the transition state for the chemical step, providing no information about the bonding environment in the Michaelis complex.2,8 Conversely, equilibrium binding isotope effects (BIE’s) arise from distinct bonding environments in the free substrate and in the Michaelis complex, prior to the chemical step.9 When both types of isotope effects are available for an enzyme reaction, the degree of distortion in the enzyme-bound substrate can be compared with that at the transition state, allowing an evaluation of ground-state destabilization.10,11 Here, the hypothesis of ground-state destabilization in the reaction catalyzed by Trypanosoma cruzi uridine phosphorylase (TcUP) (EC 2.4.2.3) is investigated by BIE analysis. TcUP catalyzes the reversible phosphorolysis of uridine to yield uracil © 2012 American Chemical Society

BIE = (14C T/14 C B − 1)/( 3HT/3 HB − 1)

(1)

BIE values for 1′-3H of 0.947 ± 0.007 (−5.3%) and for 5′-3H2 of 0.979 ± 0.006 (−2.1%) were obtained for the formation of the TcUP−uridine−sulfate ternary complex. A 1′-3H BIE of 0.953 ± 0.009 (−4.7%) was measured in the absence of (NH4)2SO4, indicating that the nucleophile mimic exerts a negligible influence on the inverse isotope effect. The sulfate ion is known to occupy the phosphate binding sites in ribosyltransferases and is employed as a nonreactive mimic of phosphate in ternary complex crystal structures and BIE experiments.10,11,16 The lack of chemical activity of this complex was confirmed by incubating 50 mM (NH4)2SO4 Received: July 9, 2012 Revised: August 6, 2012 Published: August 7, 2012 6715

dx.doi.org/10.1021/bi300914q | Biochemistry 2012, 51, 6715−6717

Biochemistry

Rapid Report

(pH 7.5), 10 μM TcUP, and 100 μM uridine in 100 mM HEPES (pH 7.5) for 90 min, followed by reverse-phase highperformance liquid chromatography separation as previously described.13 No uracil could be detected under this condition (Figure 1), which resembles that of the ultrafiltration experiments.

Figure 1. Reverse-phase chromatograms of incubation mixtures containing TcUP, uridine, and sulfate (), uracil and sulfate (−−−), and uridine and sulfate (···), in HEPES (pH 7.5).

Inverse BIE’s occur when the bond to the tritium isotope is more constrained in the enzyme-bound state. Therefore, atoms H1′ and H(2)5′ experience a more constrained bonding environment when bound to TcUP. These values contrast with the intrinsic KIE’s for 1′-3H of 1.132 (13.2%) and for 5′-3H2 of 1.041 (4.1%) for uridine arsenolysis by TcUP, which established increased degrees of bond vibrational freedom for these atoms at the transition state.13 Hence, ternary complex interactions contribute isotope effects in the opposite direction to those arising from transition-state formation. This finding is contrary to those found with Plasmodium falciparum and human orotate phosphoribosyltransferases (OPRT), where half of the normal intrinsic KIE values for 1′-3H are contributed by BIE’s,10 and with human thymidine phosphorylase (TP), where the 5′-3H2 BIE for thymidine is equivalent to and thus accounts for the entire intrinsic KIE.17 In summary, the experimental BIE data for TcUP support ground-state stabilization in both binary and ternary complexes. Density-functional theory (DFT) calculations were performed to search for a conformational geometry of uridine in the Michaelis complex that could account for the BIE’s. Several conformations of uridine were optimized in vacuo at the B3LYP level of theory with the 6-31G** basis set, as implemented in Gaussian 09,18 by varying the O4′−C1′− N1−C6 dihedral angle with either 2′-endo-3′-exo or 2′-exo-3′endo ribosyl pucker. All optimizations resulted in ground-state structures containing no imaginary frequencies. Isotope effects for each structure were calculated from scaled vibrational frequencies using ISOEFF9819 and compared with the experimental BIE values. The 2′-endo-3′-exo and 2′-exo-3′-endo uridine conformations yielded distinct BIE’s throughout the range of O4′-C1′− N1−C6 dihedrals (Figure 2). The lowest-energy conformation of free uridine, used as the reference upon which all BIE’s were calculated, adopts a 2′-endo-3′-exo pucker and an O4′−C1′− N1−C6 torsion angle of −117°13 (Figure 3, left). A 2′-exo-3′endo conformer with an O4′−C1′−N1−C6 dihedral angle of −75° (Figure 3, middle) generated calculated BIE’s for 1′-3H of

Figure 2. Dependence of calculated 1′-3H (blue) and 5′-3H2 (pink) BIE’s on O4′−C1′−N1−C6 dihedral angles of 2′-endo-3′-exo (top) and 2′-exo-3′-endo (bottom) conformers of uridine.

Figure 3. Stick models of DFT-optimized conformations of uridine as the free substrate (left), in the Michaelis complex (middle), and at the transition state for arsenolysis (right), and their respective O4′−C1′− N1−C6 dihedral angles.

0.949 (−5.1%) and for 5′-3H2 of 0.972 (−2.8%), in excellent agreement with the experimental values of −5.3 and −2.1%, respectively. Inverse BIE’s were also calculated from other conformers, but none of them matched the experimental ones as closely (Table S1 of the Supporting Information). The 2′-exo-3′-endo ribosyl pucker differs from that of free uridine and also from that of the KIE-derived transition-state structure for the TcUP-catalyzed arsenolysis of uridine. At the transition state, a flatter 2′-endo-3′-exo sugar is predicted13 (Figure 3, right). The potential energy of the conformer proposed to be present in the Michaelis complex is approximately 3 kcal mol−1 higher than that of free uridine (Figure S1 of the Supporting Information), an amount of energy stabilized, for example, by a modest H-bond.20 This could be explained by the loss of an intramolecular H-bond between O(H)5′ and O2. This H-bond was predicted by DFT optimizations of ground-state uridine, as previously reported,13 and is likely disrupted as the O4′−C1′−N1−C6 torsion angle increases from −117° to −75° upon enzyme binding. At the transition state, this angle is further increased to −23°, favoring H-bond interactions with the side chain of an active site 6716

dx.doi.org/10.1021/bi300914q | Biochemistry 2012, 51, 6715−6717

Biochemistry



glutamine residue conserved in all known uridine phosphorylases.13,21 At present, structural data for predicting interactions with the −75° dihedral conformation of uridine are lacking. Distinct atomic charges are a prominent feature of transitionstate structures of N-ribosyltransferase reactions.22 Groundstate destabilization in the Michaelis complex is therefore expected to alter the charge distribution of bound uridine toward that of the transition state. Accordingly, natural bond orbital (NBO) charges were calculated for the proposed TcUPbound uridine conformer and compared with those reported for free substrate and transition-state structures (Table 1).13

a

position

FU

BU

TS

0.25a 0.22 −0.47a −0.59a −0.78 0.22 0.22 −0.66

0.25 0.25 −0.48 −0.61 −0.76 0.23 0.21 −0.64

0.36a 0.29 −0.62a −0.51a −0.78 0.22 0.19 −0.71

*E-mail: [email protected]. Phone: (718) 4302813. Funding

This work was supported by National Institutes of Health Grants GM41916 and AI049512 to V.L.S. Notes

The authors declare no competing financial interest.



REFERENCES

(1) Jencks, W. P. (1987) Catalysis in Chemistry and Enzymology, Chapter 5, Dover Publications, New York. (2) Anderson, V. E. (2005) Arch. Biochem. Biophys. 433, 27−33. (3) Belasco, J. G., and Knowles, J. R. (1980) Biochemistry 19, 472− 477. (4) Andrews, L. D., Deng, H., and Herschlag, D. (2011) J. Am. Chem. Soc. 133, 11621−11631. (5) Warshel, A., Strajbl, M., Villa, J., and Florian, J. (2000) Biochemistry 39, 14728−14738. (6) Miller, B. G., Butterfoss, G. L., Short, S. A., and Wolfenden, R. (2001) Biochemistry 40, 6227−6232. (7) Warshel, A., and Levitt, M. (1976) J. Mol. Biol. 103, 227−249. (8) Ruszczycky, M. W., and Anderson, V. E. (2006) J. Theor. Biol. 243, 328−342. (9) Schramm, V. L. (2007) Curr. Opin. Chem. Biol. 11, 529−536. (10) Zhang, Y., and Schramm, V. L. (2011) Biochemistry 50, 4813− 4818. (11) Murkin, A. S., Birck, M. R., Rinaldo-Matthis, A., Wuxian, S., Taylor, E. A., Almo, S. C., and Schramm, V. L. (2007) Biochemistry 46, 5038−5049. (12) Paege, L. M., and Schlenk, F. (1952) Arch. Biochem. Biophys. 40, 42−49. (13) Silva, R. G., Vetticatt, M. J., Merino, E. F., Cassera, M. B., and Schramm, V. L. (2011) J. Am. Chem. Soc. 133, 9923−9931. (14) Silva, R. G., and Schramm, V. L. (2011) Biochemistry 50, 9158− 9166. (15) Murkin, A. S., Tyler, P. C., and Schramm, V. L. (2008) J. Am. Chem. Soc. 130, 2166−2167. (16) Canduri, F., dos Santos, D. M., Silva, R. G., Mendes, M. A., Basso, L. A., Palma, M. S., de Azevedo, W. F., Jr., and Santos, D. S. (2004) Biochem. Biophys. Res. Commun. 313, 907−914. (17) Birck, M. R., and Schramm, V. L. (2004) J. Am. Chem. Soc. 126, 6882−6883. (18) Frisch, M. J., et al. (2009) Gaussian 09, revision A.02, Gaussian, Inc., Wallingford, CT. (19) Anisimov, V., and Paneth, P. (1999) J. Math. Chem. 26, 75−86. (20) Lin, J., and Frey, P. A. (2000) J. Am. Chem. Soc. 122, 11258− 11259. (21) Larson, E. T., Mudeppa, D. G., Gillespie, J. R., Mueller, N., Napuli, A. J., Arif, J. A., Ross, J., Arakaki, T. L., Lauricella, A., Detitta, G., Luft, J., Zucker, F., Verlinde, C. L., Fan, E., Van Voorhis, W. C., Buckner, F. S., Rathod, P. K., Hol, W. G., and Merritt, E. A. (2010) J. Mol. Biol. 396, 1244−1259. (22) Schramm, V. L. (2011) Annu. Rev. Biochem. 80, 703−732.

Previously published values.13

The only continuous trends in charge distribution among the three structures concern the subtle accumulation of negative charge at the N1 position and the small decrease in electron density at the H1′ position of bound uridine, both of which become more pronounced at the transition state. Position C1′ conserves its charge as uridine binds TcUP and loses electron density only at the transition state as the N-ribosidic bond electrons migrate to the uracil moiety. Positions O4′, O5′, O2, pro-R-H5′, and pro-S-H5′ show altered NBO charges in the Michaelis complex in comparison with those of the free substrate. However, those charges change in the opposite direction as the system reaches the transition state. The inverse experimental 1′-3H and 5′-3H2 BIE’s presented in this work contrast with the normal KIE’s measured for TcUP-catalyzed arsenolysis of uridine. The uridine conformer whose calculated 1′- 3 H and 5′- 3 H 2 BIE’s match the experimental values representing the conformation of uridine in the Michaelis complex does not incorporate characteristics of the transition-state structure. Neither the ribosyl pucker nor the NBO charge distribution pattern is altered to more closely resemble the transition state. These observations are consistent with ground-state stabilization in the Michaelis complex of TcUP. Ground-state stabilization here refers to local effects on the bound substrate, such as bonding environment, ribosyl pucker, and NBO charges, in contrast with those in the enzymatic transition state. Compared with previous reports of BIE’s for human and P. falciparum OPRT10 and human TP,17 this work exemplifies the diverse strategies employed by different enzymes to catalyze similar chemical transformations.



AUTHOR INFORMATION

Corresponding Author

Table 1. Natural Bond Orbital Charges of Free Uridine (FU), Bound Uridine (BU), and the Transition State (TS) C1′ H1′ N1 O4′ O5′ pro-R-H5′ pro-S-H5′ O2

Rapid Report

ASSOCIATED CONTENT

S Supporting Information *

Detailed description of materials and methods, Scheme S1, Table S1, Figure S1, and atomic coordinates for all optimized ground-state geometries of uridine. This material is available free of charge via the Internet at http://pubs.acs.org. 6717

dx.doi.org/10.1021/bi300914q | Biochemistry 2012, 51, 6715−6717