Transition State Structure of the Solvolytic Hydrolysis of NAD+

Biomimetic α-selective ribosylation enables two-step modular synthesis of biologically important ADP-ribosylated peptides

The α-type ADP-ribosylated peptides represent a class of important molecular tools in the field of protein ADP-ribosylation, however, they are difficult to access because of their inherent complicated structures and the lack of effective synthetic tools. In this paper, we present a biomimetic α-selective ribosylation reaction to synthesize a key intermediate, α-ADP-ribosyl azide, directly from native β-nicotinamide adenine dinucleotide in a clean ionic liquid system. This reaction in tandem with click chemistry then offers a two-step modular synthesis of α-ADP-ribosylated peptides. These syntheses can be performed open air in eppendorf tubes, without the need for specialized instruments or training. Importantly, we demonstrate that the synthesized α-ADP-ribosylated peptides show high binding affinity and desirable stability for enriching protein partners, and reactivity in post-stage poly ADP-ribosylations. Owing to their simple chemistry and multidimensional bio-applications, the presented methods may provide a powerful platform to produce general molecular tools for the study of protein ADP-ribosylation.

ADP-ribosylated peptides are important molecular tools, however, the lack of effective synthetic methods hinders the access to their complex structures. Here, the authors present a biomimetic α-selective ribosylation reaction coupled with click chemistry ultimately providing a two-step modular synthesis of α-ADP-ribosylated peptides.

Influence of Dielectric Environment upon Isotope Effects on Glycoside Heterolysis: Computational Evaluation and Atomic Hessian Analysis

Isotope effects depend upon the polarity of the bulk medium in which a chemical process occurs. Implicit solvent calculations with molecule-shaped cavities show that the equilibrium isotope effect (EIE) for heterolysis of the glycosidic bonds in 5'-methylthioadenosine and in 2-(p-nitrophenoxy)tetrahydropyran, both in water, are very sensitive in the range 2 ≤ ε ≤ 10 to the relative permittivity of the continuum surrounding the oxacarbenium ion. However, different implementations of nominally the same PCM method can lead to opposite trends being predicted for the same molecule. Computational modeling of the influence of the inhomogeneous effective dielectric surrounding a substrate within the protein environment of an enzymic reaction requires an explicit treatment. The EIE (K_H/K_D) for transfer of cyclopentyl, cyclohexyl, tetrahydrofuranyl and tetrahydropyranyl cations from water to cyclohexane is predicted by B3LYP/6-31+G(d) calculations with implicit solvation and confirmed by B3LYP/6-31+G(d)/OPLS-AA calculations with averaging over many explicit solvation configurations. Atomic Hessian analysis, whereby the full Hessian is reduced to the elements belonging to a single atom at the site of isotopic substitution, reveals a remarkable result for both implicit and explicit solvation: the influence of the solvent environment on these EIEs is essentially captured completely by only a 3 × 3 block of the Hessian, although these values must correctly reflect the influence of the whole environment. QM/MM simulation with ensemble averaging has an important role to play in assisting the meaningful interpretation of observed isotope effects for chemical reactions both in solution and catalyzed by enzymes.

δ-Deuterium isotope effects as probes for transition-state structures of isoprenoid substrates

The biosynthetic pathways to isoprenoid compounds involve transfer of the prenyl moiety in allylic diphosphates to electron-rich (nucleophilic) acceptors. The acceptors can be many types of nucleophiles, while the allylic diphosphates only differ in the number of isoprene units and stereochemistry of the double bonds in the hydrocarbon moieties. Because of the wide range of nucleophilicities of naturally occurring acceptors, the mechanism for prenyltransfer reactions may be dissociative or associative with early to late transition states. We have measured δ-secondary kinetic isotope effects operating through four bonds for substitution reactions with dimethylallyl derivatives bearing deuterated methyl groups at the distal (C3) carbon atom in the double bond under dissociative and associative conditions. Computational studies with density functional theory indicate that the magnitudes of the isotope effects correlate with the extent of bond formation between the allylic moiety and the electron-rich acceptor in the transition state for alkylation and provide insights into the structures of the transition states for associative and dissociative alkylation reactions.

Probing the catalytic mechanism of bovine CD38/NAD+ glycohydrolase by site directed mutagenesis of key active site residues

Bovine CD38/NAD(+) glycohydrolase catalyzes the hydrolysis of NAD(+) to nicotinamide and ADP-ribose and the formation of cyclic ADP-ribose via a stepwise reaction mechanism. Our recent crystallographic study of its Michaelis complex and covalently-trapped intermediates provided insights into the modalities of substrate binding and the molecular mechanism of bCD38. The aim of the present work was to determine the precise role of key conserved active site residues (Trp118, Glu138, Asp147, Trp181 and Glu218) by focusing mainly on the cleavage of the nicotinamide-ribosyl bond. We analyzed the kinetic parameters of mutants of these residues which reside within the bCD38 subdomain in the vicinity of the scissile bond of bound NAD(+). To address the reaction mechanism we also performed chemical rescue experiments with neutral (methanol) and ionic (azide, formate) nucleophiles. The crucial role of Glu218, which orients the substrate for cleavage by interacting with the N-ribosyl 2'-OH group of NAD(+), was highlighted. This contribution to catalysis accounts for almost half of the reaction energy barrier. Other contributions can be ascribed notably to Glu138 and Asp147 via ground-state destabilization and desolvation in the vicinity of the scissile bond. Key interactions with Trp118 and Trp181 were also proven to stabilize the ribooxocarbenium ion-like transition state. Altogether we propose that, as an alternative to a covalent acylal reaction intermediate with Glu218, catalysis by bCD38 proceeds through the formation of a discrete and transient ribooxocarbenium intermediate which is stabilized within the active site mostly by electrostatic interactions.

Experimental and computational analysis of the transition state for ribonuclease A-catalyzed RNA 2'-O-transphosphorylation

Enzymes function by stabilizing reaction transition states; therefore, comparison of the transition states of enzymatic and nonenzymatic model reactions can provide insight into biological catalysis. Catalysis of RNA 2'-O-transphosphorylation by ribonuclease A is proposed to involve electrostatic stabilization and acid/base catalysis, although the structure of the rate-limiting transition state is uncertain. Here, we describe coordinated kinetic isotope effect (KIE) analyses, molecular dynamics simulations, and quantum mechanical calculations to model the transition state and mechanism of RNase A. Comparison of the (18)O KIEs on the 2'O nucleophile, 5'O leaving group, and nonbridging phosphoryl oxygens for RNase A to values observed for hydronium- or hydroxide-catalyzed reactions indicate a late anionic transition state. Molecular dynamics simulations using an anionic phosphorane transition state mimic suggest that H-bonding by protonated His12 and Lys41 stabilizes the transition state by neutralizing the negative charge on the nonbridging phosphoryl oxygens. Quantum mechanical calculations consistent with the experimental KIEs indicate that expulsion of the 5'O remains an integral feature of the rate-limiting step both on and off the enzyme. Electrostatic interactions with positively charged amino acid site chains (His12/Lys41), together with proton transfer from His119, render departure of the 5'O less advanced compared with the solution reaction and stabilize charge buildup in the transition state. The ability to obtain a chemically detailed description of 2'-O-transphosphorylation transition states provides an opportunity to advance our understanding of biological catalysis significantly by determining how the catalytic modes and active site environments of phosphoryl transferases influence transition state structure.

Automated Transition State Searches without Evaluating the Hessian

Accurate and speedy determination of transition structures (TSs) is essential for computational studies on reaction pathways, particularly when the process involves expensive electronic structure calculations. Many search algorithms require a good initial guess of the TS geometry, as well as a Hessian input that possesses a structure consistent with the desired saddle point. Among the double-ended interpolation methods for generation of the guess for the TS, the freezing string method (FSM) is proven to be far less expensive compared to its predecessor, the growing string method (GSM). In this paper, it is demonstrated that the efficiency of this technique can be improved further by replacing the conjugate gradient optimization step (FSM-CG) with a quasi-Newton line search coupled with a BFGS Hessian update (FSM-BFGS). A second crucial factor that affects the speed with which convergence to the TS is achieved is the quality and cost of the Hessian of the energy for the guessed TS. For electronic structure calculations, the cost of calculating an exact Hessian increases more rapidly with system size than the energy and gradient. Therefore, to sidestep calculation of the exact Hessian, an approximate Hessian is constructed, using the tangent direction and local curvature at the TS guess. It is demonstrated that the partitioned-rational function optimization algorithm for locating TSs with this approximate Hessian input performs at least as well as with an exact Hessian input in most test cases. The two techniques, FSM and approximate Hessian construction, therefore can significantly reduce costs associated with finding TSs.

Transition state analysis of acid-catalyzed hydrolysis of an enol ether, enolpyruvylshikimate 3-phosphate (EPSP)

Proton transfer to carbon represents a significant catalytic challenge because of the large intrinsic energetic barrier and the frequently unfavorable thermodynamics. Multiple kinetic isotope effects (KIEs) were measured for acid-catalyzed hydrolysis of the enol ether functionality of enolpyruvylshikimate 3-phosphate (EPSP) as a nonenzymatic analog of the EPSP synthase (AroA) reaction. The large solvent deuterium KIE demonstrated that protonating C3 was the rate-limiting step, and the lack of solvent hydron exchange into EPSP demonstrated that protonation was irreversible. The reaction mechanism was stepwise, with C3, the methylene carbon, being protonated to form a discrete oxacarbenium ion intermediate before water attack at the cationic center, that is, an AH(‡)*AN (or AH(‡) + AN) mechanism. The calculated 3-(14)C and 3,3-(2)H2 KIEs varied as a function of the extent of proton transfer at the transition state, as reflected in the C3-H(+) bond order, nC3-H+. The calculated 3-(14)C KIE was a function primarily of C3 coupling with the movement of the transferring proton, as reflected in the reaction coordinate contribution ((light)ν(‡)/(heavy)ν(‡)), rather than of changes in bonding. Coupling was strongest in early and late transition states, where the reaction coordinate frequency was lower. The other calculated (14)C and (18)O KIEs were more sensitive to interactions with counterions and solvation in the model structures than nC3-H+. The KIEs revealed a moderately late transition state with significant oxacarbenium ion character and with a C3-H(+) bond order ≈0.6.

Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines

Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(‡) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.

Transition-state analysis of Trypanosoma cruzi uridine phosphorylase-catalyzed arsenolysis of uridine

Uridine phosphorylase catalyzes the reversible phosphorolysis of uridine and 2'-deoxyuridine to generate uracil and (2-deoxy)ribose 1-phosphate, an important step in the pyrimidine salvage pathway. The coding sequence annotated as a putative nucleoside phosphorylase in the Trypanosoma cruzi genome was overexpressed in Escherichia coli , purified to homogeneity, and shown to be a homodimeric uridine phosphorylase, with similar specificity for uridine and 2'-deoxyuridine and undetectable activity toward thymidine and purine nucleosides. Competitive kinetic isotope effects (KIEs) were measured and corrected for a forward commitment factor using arsenate as the nucleophile. The intrinsic KIEs are: 1'-(14)C = 1.103, 1,3-(15)N(2) = 1.034, 3-(15)N = 1.004, 1-(15)N = 1.030, 1'-(3)H = 1.132, 2'-(2)H = 1.086, and 5'-(3)H(2) = 1.041 for this reaction. Density functional theory was employed to quantitatively interpret the KIEs in terms of transition-state structure and geometry. Matching of experimental KIEs to proposed transition-state structures suggests an almost synchronous, S(N)2-like transition-state model, in which the ribosyl moiety possesses significant bond order to both nucleophile and leaving groups. Natural bond orbital analysis allowed a comparison of the charge distribution pattern between the ground-state and the transition-state models.

Transition state of ADP-ribosylation of acetyllysine catalyzed by Archaeoglobus fulgidus Sir2 determined by kinetic isotope effects and computational approaches

Sirtuins are protein-modifying enzymes distributed throughout all forms of life. These enzymes bind NAD(+), a universal metabolite, and react it with acetyllysine residues to effect deacetylation of protein side chains. This NAD(+)-dependent deacetylation reaction has been observed for sirtuin enzymes derived from archaeal, eubacterial, yeast, metazoan, and mammalian species, suggesting conserved chemical mechanisms for these enzymes. The first chemical step of deacetylation is the reaction of NAD(+) with an acetyllysine residue which forms an enzyme-bound ADPR-peptidylimidate intermediate and nicotinamide. In this manuscript, the transition state for the ADP-ribosylation of acetyllysine is solved for an Archaeoglobus fulgidus sirtuin (Af2Sir2). Kinetic isotope effects (KIEs) were obtained by the competitive substrate method and were [1(N)-(15)N] = 1.024(2), [1'(N)-(14)C] = 1.014(4), [1'(N)-(3)H] = 1.300(3), [2'(N)-(3)H] = 1.099(5), [4'(N)-(3)H] = 0.997(2), [5'(N)-(3)H] = 1.020(5), [4'(N)-(18)O] = 0.984(5). KIEs were calculated for candidate transition state structures using computational methods (Gaussian 03 and ISOEFF 98) in order to match computed and experimentally determined KIEs to solve the transition state. The results indicate that the enzyme stabilizes a highly dissociated oxocarbenium ionlike transition state with very low bond orders to the leaving group nicotinamide and the nucleophile acetyllysine. A concerted yet highly asynchronous substitution mechanism forms the ADPR-peptidylimidate intermediate of the sirtuin deacetylation reaction.

Factors controlling the mechanism of NAD(+) non-redox reactions

β-Nicotinamide adenine dinucleotide (NAD(+)) is an indispensable coenzyme or substrate for enzymes involved in catalyzing redox and non-redox reactions. ADP-ribosylating enzymes catalyze cleavage of the nicotinamide-glycosyl bond of NAD(+) and addition of a nucleophilic group from their substrate proteins to the N-ribose anomeric carbon of NAD(+). Although the role of the nicotinamide-ribose fragment in the mechanism of NAD(+) hydrolysis has been examined, the role of the doubly negatively charged, flexible, and chemically reactive NAD(+) diphosphate moiety in the reaction process has largely been neglected. Thus, the participation of the pyrophosphate group in stabilizing intra- and intermolecular interactions in the ground state and transition state has not been explored. Furthermore, the roles of other factors such as the type/nucleophilicity of the attacking nucleophile and the medium in influencing the reaction pathway have not been systematically evaluated. In this study, we endeavor to fill in these gaps and elucidate the role of these factors in controlling the NAD(+) nicotinamide-glycosyl bond cleavage. Using density functional theory combined with continuum dielectric methods, we modeled both S(N)1 and S(N)2 reaction pathways and assessed the role of the diphosphate group in stabilizing the (i) NAD(+) ground state, (ii) oxocarbocation intermediate, (iii) reaction product, and (iv) nucleophile. We also assessed the chemical nature of the attacking nucleophile and the role of the protein matrix in affecting the reaction mechanism. Our results reveal an intricate interplay among various factors in controlling the reaction pathway, which in turn suggests ways in which the enzyme can accelerate the reaction.

Sirtuin chemical mechanisms

Sirtuins are ancient proteins widely distributed in all lifeforms of earth. These proteins are universally able to bind NAD(+), and activate it to effect ADP-ribosylation of cellular nucleophiles. The most commonly observed sirtuin reaction is the ADP-ribosylation of acetyllysine, which leads to NAD(+)-dependent deacetylation. Other types of ADP-ribosylation have also been observed, including protein ADP-ribosylation, NAD(+) solvolysis and ADP-ribosyltransfer to 5,6-dimethylbenzimidazole, a reaction involved in eubacterial cobalamin biosynthesis. This review broadly surveys the chemistries and chemical mechanisms of these enzymes.

Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations

Sir2 enzymes catalyze the NAD+-dependent protein deacetylation and play critical roles in epigenetics, cell death, and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm-catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

Transition-state analysis of the DNA repair enzyme MutY

The transition state (TS) structure of MutY-catalyzed DNA hydrolysis was solved using multiple kinetic isotope effect (KIE) measurements. MutY is a base excision repair enzyme which cleaves adenine from 8-oxo-G:A mismatches in vivo, and also from G:A mismatches in vitro. TS analysis of G:A-DNA hydrolysis revealed a stepwise S(N)1 (D(N)*A(N)(double dagger)) mechanism proceeding through a highly reactive oxacarbenium ion intermediate which would have a lifetime in solution of <10(-10) s. C-N bond cleavage is reversible; the N-glycoside bond breaks and reforms repeatedly before irreversible water attack on the oxacarbenium ion. KIEs demonstrated that MutY uses general acid catalysis by protonating N7. It enforces a 3'-exo sugar ring conformation and other sugar ring distortions to stabilize the oxacarbenium ion. Combining the experimental TS structure with the previously reported crystal structure of an abortive Michaelis complex elucidates the step-by-step catalytic sequence.

Poly(ADP-ribose)-polymerase-catalyzed hydrolysis of NAD+: QM/MM simulation of the enzyme reaction

Poly(ADP-ribose) polymerase (PARP) is a nuclear enzyme which uses NAD+ as substrate and catalyzes the transfer of multiple units of ADP-ribose to target proteins. PARP is an attractive target for the discovery of novel therapeutic agents and PARP inhibitors are currently evaluated for the treatment of a variety of pathological conditions such as brain ischemia, inflammation, and cancer. Herein, we use the PARP-catalyzed reaction of NAD+ hydrolysis as a model for gaining insight into the molecular details of the catalytic mechanism of PARP. The reaction has been studied in both the gas-phase and in the enzyme environment through a QM/MM approach. Our results indicate that the cleavage reaction of the nicotinamide-ribosyl bond proceeds through an SN2 dissociative mechanism via an oxacarbenium transition structure. These results confirm the importance of the structural water molecule in the active site and may constitute the basis for the design of transition-state-based PARP inhibitors.

Transition state analysis of adenosine nucleosidase from yellow lupin (Lupinus luteus)

The transition state of adenosine nucleosidase (EC 3.2.2.7) isolated from yellow lupin (Lupinus luteus) was determined based upon a series of heavy atom kinetic isotope effects. Adenosine labeled with 13C, 2H, and 15N was analyzed by liquid chromatography/electrospray mass spectrometry to determine kinetic isotope effects. Values of 1.024+/-0.004, 1.121+/-0.005, 1.093+/-0.004, 0.993+/-0.006, and 1.028+/-0.005 were found for [1'-13C], [1'-2H], [2'-2H], [5'-2H], and [9-15N] adenosine, respectively. Using a bond order bond energy vibrational analysis, a transition state consisting of a significantly broken C-N bond, formation of an oxocarbenium ion in the ribose ring, a conformation of C3-exo for the ribose ring, and protonation of the heterocyclic base was proposed. This transition state was found to be very similar to the transition state for nucleoside hydrolase, another purine metabolizing enzyme, isolated from Crithidia fasciculata.

Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase

It has been suggested that the magnitudes of secondary kinetic isotope effects (2 degrees KIEs) of enzyme-catalyzed reactions are an indicator of the extent of reaction-center rehybridization at the transition state. A 2 degrees KIE value close to the corresponding secondary equilibrium isotope effects (2 degrees EIE) is conventionally interpreted as indicating a late transition state that resembles the final product. The reliability of using this criterion to infer the structure of the transition state is examined by carrying out a theoretical investigation of the hybridization states of the hydride donor and acceptor in the Escherichia coli dihydrofolate reductase (ecDHFR)-catalyzed reaction for which a 2 degrees KIE close to the 2 degrees EIE was reported. Our results show that the donor carbon at the hydride transfer transition state resembles the reactant state more than the product state, whereas the acceptor carbon is more productlike, which is a symptom of transition state imbalance. The conclusion that the isotopically substituted carbon is reactant-like disagrees with the conclusion that would have been derived from the criterion of 2 degrees KIEs and 2 degrees EIEs, but the breakdown of the correlation with the equilibrium isotope effect can be explained by considering the effect of tunneling.

Calculations of Substituent and Solvent Effects on the Kinetic Isotope Effects of Menshutkin Reactions

Nitrogen, deuterium, halogen, and carbon kinetic isotope effects have been modeled for the Menshutkin reaction between methyl halides and substituted N,N-dimethylaniline at the HF/6-31G(d) level of theory augmented by the C-PCM continuum solvent model for several solvents. Systematic changes in geometries of the transition states and Gibbs free energies of activation have been found with phenyl ring substituents, solvent, and the leaving group. Kinetic isotope effects also change systematically; however, these changes are predicted to be small, inside the usual precision of the experimental measurements. On the contrary, no correlation has been found between the kinetic isotope effects and the Hammett constants for para substituents. Thus opposite to previous assumptions, our results indicate that kinetic isotope effects on the Menshutkin reaction cannot be used to predict the position of the transition state on the reaction coordinate.

Experimental and theoretical multiple kinetic isotope effects for an SN2 reaction. An attempt to determine transition-state structure and the ability of theoretical methods to predict experimental kinetic isotope effects

The secondary alpha-deuterium, the secondary beta-deuterium, the chlorine leaving-group, the nucleophile secondary nitrogen, the nucleophile (12)C/(13)C carbon, and the (11)C/(14)C alpha-carbon kinetic isotope effects (KIEs) and activation parameters have been measured for the S(N)2 reaction between tetrabutylammonium cyanide and ethyl chloride in DMSO at 30 degrees C. Then, thirty-nine readily available different theoretical methods, both including and excluding solvent, were used to calculate the structure of the transition state, the activation energy, and the kinetic isotope effects for the reaction. A comparison of the experimental and theoretical results by using semiempirical, ab initio, and density functional theory methods has shown that the density functional methods are most successful in calculating the experimental isotope effects. With two exceptions, including solvent in the calculation does not improve the fit with the experimental KIEs. Finally, none of the transition states and force constants obtained from the theoretical methods was able to predict all six of the KIEs found by experiment. Moreover, none of the calculated transition structures, which are all early and loose, agree with the late (product-like) transition-state structure suggested by interpreting the experimental KIEs.

Exploring nucleoside hydrolase catalysis in silico: molecular dynamics study of enzyme-bound substrate and transition state

The mechanism of action of inosine-uridine nucleoside hydrolase has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Five MD studies have been performed with enzyme substrate complex (E.S), enzyme substrate complex with protonated His241 (EH.S), enzyme transition state complex (E.TS), enzyme transition state complex with protonated His241 (EH.TS), and His241Ala transition state complex E(H241A).TS. Special attention has been given to the role of His241, which has been considered as the general acid catalyst to assist departure of the leaving nucleobase on the basis of its location in the active site in the X-ray crystal structure (). Yet on the basis of the location in the active site, Tyr229 is closer to the aniline ring of pAPIR as compared to His241. On initiation of MD simulations, His241 does not approach the nucleobase in the structures of EH.S, E.S, EH.TS, and E.TS. In the solvated enzyme, Tyr229, which is a member of the hydrogen bonding network inosine O2'.Asp14.His241.Tyr229.inosine N7, serves as a proton source to the leaving nucleobase. The loss of significant activity of His241Ala mutant is shown to be related to the disruption of the above hydrogen bonded network and the distancing of Tyr229 from inosine N7. The structures of the enzyme complexes with substrate or TS are not visibly altered on protonation of His241, a most unusual outcome. The bell-shaped pH dependence upon pK(app)'s of 7.1 and 9.1 may be attributed to the necessity of the dissociation of Asp10 or Asp15 and the acid form of Tyr229, respectively. In TS, the residue Ile81 migrated closer, whereas Arg233 moved away from the nucleobase. The probability of ribooxocarbenium ion stabilization by Asn168 and Asp14 is discussed. The Asp14-CO(2)(-) is hydrogen bonded to the ribose 2'-OH for 96% of the MD simulation time. Nucleophilic addition of water138 to ribooxocarbenium ion is suggested to be assisted by the proton shuttle from water138 --> Asp10 --> Asp15 --> water pool. An anticorrelation motion between Tyr229-OH and Asn168-OD1 in EH.S and E.S is observed. The relationship of this anticorrelated motion to mechanism, if any, deserves further exploration, perhaps the formation of a near attack conformation.

Transition state variation in enzymatic reactions

Experimental analysis of enzymatic transition states by kinetic isotope effect methods has established geometric variation in related transition state structures. Differences are apparent in development of the reaction coordinate, in solvolytic transition states relative to those in enzymatic catalytic sites, in the stereochemistry of related substrates at the transition state, and in reactions catalyzed by related enzymes.

Determining transition states from kinetic isotope effects

BEBOVA-based TS determination has been very successful in elucidating enzyme mechanisms at a level of detail that would be otherwise inaccessible. The resulting TS structures have been used successfully as the basis for designing TS mimics as enzyme inhibitors with dissociation constants to 10(-11) M. The structure interpolation approach has systematized the process of finding a TS, increasing both the speed and the accuracy of TS determination. The combination of information from several TSs into a unified model increases the accuracy of the process significantly and results in an extremely sensitive probe of changes in TS with varying reaction conditions (i.e., enzymatic vs nonenzymatic reactions, different enzymes, or different nucleophiles). The TS determination process is summarized in Fig. 15.

Enzymatic transition states and transition state analog design

All chemical transformations pass through an unstable structure called the transition state, which is poised between the chemical structures of the substrates and products. The transition states for chemical reactions are proposed to have lifetimes near 10(-13) sec, the time for a single bond vibration. No physical or spectroscopic method is available to directly observe the structure of the transition state for enzymatic reactions. Yet transition state structure is central to understanding catalysis, because enzymes function by lowering activation energy. An accepted view of enzymatic catalysis is tight binding to the unstable transition state structure. Transition state mimics bind tightly to enzymes by capturing a fraction of the binding energy for the transition state species. The identification of numerous transition state inhibitors supports the transition state stabilization hypothesis for enzymatic catalysis. Advances in methods for measuring and interpreting kinetic isotope effects and advances in computational chemistry have provided an experimental route to understand transition state structure. Systematic analysis of intrinsic kinetic isotope effects provides geometric and electronic structure for enzyme-bound transition states. This information has been used to compare transition states for chemical and enzymatic reactions; determine whether enzymatic activators alter transition state structure; design transition state inhibitors; and provide the basis for predicting the affinity of enzymatic inhibitors. Enzymatic transition states provide an understanding of catalysis and permit the design of transition state inhibitors. This article reviews transition state theory for enzymatic reactions. Selected examples of enzymatic transition states are compared to the respective transition state inhibitors.