Site-bound water and the shortcomings of a less than perfect transition state analogue

Cytosine deamination and the precipitous decline of spontaneous mutation during Earth's history

The hydrolytic deamination of cytosine and 5-methylcytosine residues in DNA appears to contribute significantly to the appearance of spontaneous mutations in microorganisms and in human disease. In the present work, we examined the mechanism of cytosine deamination and the response of the uncatalyzed reaction to changing temperature. The positively charged 1,3-dimethylcytosinium ion was hydrolyzed at a rate similar to the rate of acid-catalyzed hydrolysis of 1-methylcytosine, for which it furnishes a satisfactory kinetic model and a probable mechanism. In agreement with earlier reports, uncatalyzed deamination was found to proceed at very similar rates for cytosine, 1-methylcytosine, cytidine, and cytidine 5'-phosphate, and also for cytosine residues in single-stranded DNA generated from a phagemid, in which we sequenced an insert representing the gene of the HIV-1 protease. Arrhenius plots for the uncatalyzed deamination of cytosine were linear over the temperature range from 90 °C to 200 °C and indicated a heat of activation (ΔH(‡)) of 23.4 ± 0.5 kcal/mol at pH 7. Recent evidence indicates that the surface of the earth has been cool enough to support life for more than 4 billion years and that life has been present for almost as long. If the temperature at Earth's surface is assumed to have followed Newton's law of cooling, declining exponentially from 100 °C to 25 °C during that period, then half of the cytosine-deaminating events per unit biomass would have taken place during the first 0.2 billion years, and <99.4% would have occurred during the first 2 billion years.

Synthesis and conformational analysis of locked carbocyclic analogues of 1,3-diazepinone riboside, a high-affinity cytidine deaminase inhibitor

Cytidine deaminase (CDA) catalyzes the deamination of cytidine via a hydrated transition-state intermediate that results from the nucleophilic attack of zinc-bound water at the active site. Nucleoside analogues where the leaving NH(3) group is replaced by a proton and prevent conversion of the transition state to product are very potent inhibitors of the enzyme. However, stable carbocyclic versions of these analogues are less effective as the role of the ribose in facilitating formation of hydrated species is abolished. The discovery that a 1,3-diazepinone riboside (4) operated as a tight-binding inhibitor of CDA independent of hydration provided the opportunity to study novel inhibitors built as conformationally locked, carbocyclic 1,3-diazepinone nucleosides to determine the enzyme's conformational preference for a specific form of sugar pucker. This work describes the synthesis of two target bicyclo[3.1.0]hexane nucleosides, locked as north (5) and south (6) conformers, as well as a flexible analogue (7) built with a cyclopentane ring. The seven-membered 1,3-diazepinone ring in all the three targets was built from the corresponding benzoyl-protected carbocyclic bis-allyl ureas by ring-closing metathesis. The results demonstrate CDA's binding preference for a south sugar pucker in agreement with the high-resolution crystal structures of other CDA inhibitors bound at the active site.

Analysis of PUGNAc and NAG-thiazoline as transition state analogues for human O-GlcNAcase: mechanistic and structural insights into inhibitor selectivity and transition state poise

O-GlcNAcase catalyzes the cleavage of beta-O-linked 2-acetamido-2-deoxy-beta-d-glucopyranoside (O-GlcNAc) from serine and threonine residues of post-translationally modified proteins. Two potent inhibitors of this enzyme are O-(2-acetamido-2-deoxy-d-glucopyranosylidene)amino-N-phenylcarbamate (PUGNAc) and 1,2-dideoxy-2'-methyl-alpha-d-glucopyranoso[2,1-d]-Delta2'-thiazoline (NAG-thiazoline). Derivatives of these inhibitors differ in their selectivity for human O-GlcNAcase over the functionally related human lysosomal beta-hexosamindases, with PUGNAc derivatives showing modest selectivities and NAG-thiazoline derivatives showing high selectivities. The molecular basis for this difference in selectivities is addressed as is how well these inhibitors mimic the O-GlcNAcase-stabilized transition state (TS). Using a series of substrates, ground state (GS) inhibitors, and transition state mimics having analogous structural variations, we describe linear free energy relationships of log(KM/kcat) versus log(KI) for PUGNAc and NAG-thiazoline. These relationships suggest that PUGNAc is a poor transition state analogue, while NAG-thiazoline is revealed as a transition state mimic. Comparative X-ray crystallographic analyses of enzyme-inhibitor complexes reveal subtle molecular differences accounting for the differences in selectivities between these two inhibitors and illustrate key molecular interactions. Computational modeling of species along the reaction coordinate, as well as PUGNAc and NAG-thiazoline, provide insight into the features of NAG-thiazoline that resemble the transition state and reveal where PUGNAc fails to capture significant binding energy. These studies also point to late transition state poise for the O-GlcNAcase catalyzed reaction with significant nucleophilic participation and little involvement of the leaving group. The potency of NAG-thiazoline, its transition state mimicry, and its lack of traditional transition state-like design features suggest that potent rationally designed glycosidase inhibitors can be developed that exploit variation in transition state poise.

Origin of tight binding of a near-perfect transition-state analogue by cytidine deaminase: implications for enzyme catalysis

Cytidine deaminase (CDA) is a zinc metalloenzyme that catalyzes the hydrolytic deamination of cytidine to uridine. Zebularine (ZEB) binds to CDA, and the binding process leads to a near-perfect transition-state analogue (TSA) inhibitor at the active site with an estimated K(i) value of 1.2 x 10(-)(12) M. The interaction of CDA with the TSA inhibitor has become a paradigm for studying the tight TSA binding by enzymes. The formation of the TSA is catalyzed by CDA by a mechanism that is similar to the formation of the tetrahedral intermediate during the CDA-catalyzed reaction (i.e., through the nucleophilic attack of a Zn-hydroxide group on C(4)). It is believed that the TSA formed at the active site is zebularine 3,4-hydrate. In this paper, it is shown from QM/MM molecular dynamics and free energy simulations that zebularine 3,4-hydrate may in fact be unstable in the enzyme and that a proton transfer from the Zn-hydroxide group to Glu-104 during the nucleophilic attack could be responsible for the very high affinity. The nucleophilic attack by the Zn-hydroxide on C(4) is found to be concerted with two proton transfers. Such concerted process allows the TSA, an alkoxide-like inhibitor, to be stabilized through a mechanism that is similar to the transition-state stabilization in the general acid-base catalysis. It is suggested that the proton transfer from the Zn-hydroxide to Glu-104, which is required to generate the general acid for protonating the leaving ammonia, may play an important role in lowering the activation barrier during the catalysis.

Structural, Thermodynamic, and Kinetic Analyses of Tetrahydrooxazine-derived Inhibitors Bound to β-Glucosidases

The understanding of transition state mimicry in glycoside hydrolysis is increasingly important both in the quest for novel specific therapeutic agents and for the deduction of enzyme function and mechanism. To aid comprehension, inhibitors can be characterized through kinetic, thermodynamic, and structural dissection to build an "inhibition profile." Here we dissect the binding of a tetrahydrooxazine inhibitor and its derivatives, which display Ki values around 500 nm. X-ray structures with both a beta-glucosidase, at 2 A resolution, and an endoglucanase at atomic (approximately 1 A) resolution reveal similar interactions between the tetrahydrooxazine inhibitor and both enzymes. Kinetic analyses reveal the pH dependence of kcat/Km and 1/Ki with both enzyme systems, and isothermal titration calorimetry unveils the enthalpic and entropic contributions to beta-glucosidase inhibition. The pH dependence of enzyme activity mirrored that of 1/Ki in both enzymes, unlike the cases of isofagomine and 1-deoxynojirimycin that have been characterized previously. Calorimetric dissection reveals a large favorable enthalpy that is partially offset by an unfavorable entropy upon binding. In terms of the similar profile for the pH dependence of 1/Ki and the pH dependence of kcat/Km, the significant enthalpy of binding when compared with other glycosidase inhibitors, and the tight binding at the optimal pH of the enzymes tested, tetrahydrooxazine and its derivatives are a significantly better class of glycosidase inhibitor than previously assumed.

Fourier transform ion cyclotron resonance MS reveals the presence of a water molecule in an enzyme transition-state analogue complex

The structures of several powerful inhibitors of hydrolytic enzymes resemble that of the altered substrate in the transition state, except that a hydrogen atom replaces one substituent (typically the leaving group). To test the hypothesis that a water molecule might be present in the gap resulting from this replacement, we examined a transition-state analogue complex formed by Escherichia coli cytidine deaminase by Fourier transform ion cyclotron resonance MS in electrospray mode. Upon nebularization from aqueous solution under conditions (pH 5.6) where the enzyme is active, cytidine deaminase remains dimeric in the vapor phase. In the presence of inhibitor, the enzyme's exact mass can be used to infer the presence at each active site of zinc, 5-fluoro-3,4-dihydrouridine, and a single water molecule.

Hydration free energies and entropies for water in protein interiors

Free energy calculations for the transfer of a water molecule from the pure liquid to an interior cavity site in a protein are presented. Two different protein cavities, in bovine pancreatic trypsin inhibitor (BPTI) and in the I76A mutant of barnase, represent very different environments for the water molecule: one which is polar, forming four water-protein hydrogen bonds, and one which is more hydrophobic, forming only one water-protein hydrogen bond. The calculations give very different free energies for the different cavities, with only the polar BPTI cavity predicted to be hydrated. The corresponding entropies for the transfer to the interior cavities are calculated as well and show that the transfer to the polar cavity is significantly entropically unfavorable while the transfer to the nonpolar cavity is entropically favorable. For both proteins an analysis of the fluctuations in the positions of the protein atoms shows that the addition of a water molecule makes the protein more flexible. This increased flexibility appears to be due to an increased length and weakened strength of protein-protein hydrogen bonds near the cavity.

Iminosugar Glycosidase Inhibitors: Structural and Thermodynamic Dissection of the Binding of Isofagomine and 1-Deoxynojirimycin to β-Glucosidases

The design and synthesis of transition-state mimics reflects the growing need both to understand enzymatic catalysis and to influence strategies for therapeutic intervention. Iminosugars are among the most potent inhibitors of glycosidases. Here, the binding of 1-deoxynojirimycin and (+)-isofagomine to the "family GH-1" beta-glucosidase of Thermotoga maritima is investigated by kinetic analysis, isothermal titration calorimetry, and X-ray crystallography. The binding of both of these iminosugar inhibitors is driven by a large and favorable enthalpy. The greater inhibitory power of isofagomine, relative to 1-deoxynojirimycin, however, resides in its significantly more favorable entropy; indeed the differing thermodynamic signatures of these inhibitors are further highlighted by the markedly different heat capacity values for binding. The pH dependence of catalysis and of inhibition suggests that the inhibitory species are protonated inhibitors bound to enzymes whose acid/base and nucleophile are ionized, while calorimetry indicates that one proton is released from the enzyme upon binding at the pH optimum of catalysis (pH 5.8). Given that these results contradict earlier proposals that the binding of racemic isofagomine to sweet almond beta-glucosidase was entropically driven (Bülow, A. et al. J. Am. Chem. Soc. 2000, 122, 8567-8568), we reinvestigated the binding of 1-deoxynojirimycin and isofagomine to the sweet almond enzyme. Calorimetry confirms that the binding of isofagomine to sweet almond beta-glucosidases is, as observed for the T. maritima enzyme, driven by a large favorable enthalpy. The crystallographic structures of the native T. maritima beta-glucosidase, and its complexes with isofagomine and 1-deoxynojirimycin, all at approximately 2.1 A resolution, reveal that additional ordering of bound solvent may present an entropic penalty to 1-deoxynojirimycin binding that does not penalize isofagomine.

Thermodynamic and extrathermodynamic requirements of enzyme catalysis

An enzyme's affinity for the altered substrate in the transition state (symbolized here as S) matches the value of k(cat)/K(m) divided by the rate constant for the uncatalyzed reaction in water. The validity of this relationship is not affected by the detailed mechanism by which any particular enzyme may act, or on whether changes in enzyme conformation occur on the path to the transition state. It subsumes potential effects of substrate desolvation, H-bonding and other polar attractions, and the juxtaposition of several substrates in a configuration appropriate for reaction. The startling rate enhancements that some enzymes produce have only recently been recognized. Direct measurements of the binding affinities of stable transition-state analog inhibitors confirm the remarkable power of binding discrimination of enzymes. Several parts of the enzyme and substrate, that contribute to S binding, exhibit extremely large connectivity effects, with effective relative concentrations in excess of 10(8) M. Exact structures of enzyme complexes with transition-state analogs also indicate a general tendency of enzyme active sites to close around S in such a way as to maximize binding contacts. The role of solvent water in these binding equilibria, for which Walter Kauzmann provided a primer, is only beginning to be appreciated.

Enzymatic conformational fluctuations along the reaction coordinate of cytidine deaminase

Analysis of the crystal structures for cytidine deaminase complexed with substrate analog 3-deazacytidine, transition-state analog zebularine 3,4-hydrate, and product uridine establishes significant changes in the magnitude of atomic-scale fluctuations along the (approximate) reaction coordinate of this enzyme. Differences in fluctuations between the substrate analog complex, transition-state analog complex, and product complex are monitored via changes in corresponding crystallographic temperature factors. Previously, we reported that active-site conformational disorder is substantially reduced in the transition-state complex relative to the two ground-state complexes. Here, this result is statistically corroborated by crystallographic data for fluorinated zebularine 3,4-hydrate, a second transition-state analog, and by multiple regression analysis. Multiple regression explains 70% of the total temperature factor variation through a predictive model for the average B-value of an amino acid as a function of the catalytic state of the enzyme (substrate, transition state, product) and five other physical and structural descriptors. Furthermore, correlations of atomic fluctuation magnitudes throughout the body of each complex are quantified through an auto-correlation function. The transition-state analog complex shows the greatest correlations between temperature factor magnitudes for spatially separated atoms, underscoring the strong ability of this reaction-coordinate species to "organize" enzymatic fluctuations. The catalytic significance for decreased atomic-scale motions in the transition state is discussed. A thermodynamic argument indicates that the significant decreases in local enzymatic conformational entropy at the transition state result in enhanced energetic stabilization there.

Catalysis by entropic effects: the action of cytidine deaminase on 5,6-dihydrocytidine

In neutral solution, 5,6-dihydrocytidine undergoes spontaneous deamination (k25 approximately 3.2 x 10(-5) s(-1)) much more rapidly than does cytidine (k25 approximately 3.0 x 10(-10) s(-1)), with a more favorable enthalpy of activation (DeltaDeltaH# = -8.7 kcal/mol) compensated by a less favorable entropy of activation (TDeltaDeltaS# = -1.8 kcal/mol at 25 degrees C). E. coli cytidine deaminase enhances the rate of deamination of 5,6-dihydrocytidine (kcat/k(non) = 4.4 x 10(5)) by enhancing the entropy of activation (DeltaDeltaH# = 0 kcal/mol; TDeltaDeltaS# = +7.6 kcal/mol, at 25 degrees C). Binding of the competitive inhibitor 3,4,5,6-tetrahydrouridine (THU), a stable analogue of 5,6-dihydrocytidine in the transition state for its deamination, is accompanied by a release of enthalpy (DeltaH = -7.1 kcal/mol, TDeltaDeltaS = +2.2 kcal/mol) that approaches the estimated enthalpy of binding of the actual substrate in the transition state for deamination of 5,6-dihydrocytidine (DeltaH = -8.1 kcal/mol, TDeltaDeltaS = +6.0 kcal/mol). Thus, the shortcomings of THU in capturing all of the binding affinity expected of an ideal transition-state analogue reflect a less favorable entropy of association. That difference may arise from the analogue's inability to displace a water molecule from the "leaving group site" at which ammonia is generated in the normal reaction. The effect on binding of removing the 4-OH group from the transition-state analogue THU, to form 3,4,5,6-tetrahydrozebularine (THZ) (DeltaDeltaH = -2.1 kcal/mol, TDeltaDeltaS = -4.4 kcal/mol), is mainly entropic, consistent with the inability of THZ to displace water from the "attacking group site". These results are consistent with earlier indications [Snider, M. J., and Wolfenden, R. (2001) Biochemistry 40, 11364] that site-bound water plays a prominent role in substrate activation and inhibitor binding by cytidine deaminase.

The depth of chemical time and the power of enzymes as catalysts

The fastest known reactions include reactions catalyzed by enzymes, but the rate enhancements that enzymes produce had not been fully appreciated until recently. In the absence of enzymes, these same reactions are among the slowest that have ever been measured, some with half-times approaching the age of the Earth. This difference provides a measure of the proficiencies of enzymes as catalysts and their relative susceptibilities to inhibition by transition-state analogue inhibitors. Thermodynamic comparisons between spontaneous and enzyme-catalyzed reactions, coupled with structural information, suggest that in addition to electrostatic and H-bonding interactions, the liberation of water molecules from an enzyme's active site into bulk solvent sometimes plays a prominent role in determining the relative binding affinities of the altered substrate in the ground state and transition state. These comparisons also indicate a high level of synergism in the action of binding determinants of both the substrate and the enzyme, that are not directly involved in the chemical transformation of the substrate but contribute to the rate of its transformation at an enzyme's active site.