Search

Redox-active residues, such as tyrosine and tryptophan, play important roles in a wide range of biological processes. The α 3 Y de novo protein, which is composed of three α helices and a tyrosine residue Y32, provides a platform for investigating the redox properties of tyrosine in a well-defined protein environment. Herein, the proton-coupled electron transfer (PCET) reaction that occurs upon oxidation of tyrosine in this model protein by a ruthenium photosensitizer is studied by using a vibronically nonadiabatic PCET theory that includes hydrogen tunneling and excited vibronic states. The input quantities to the analytical nonadiabatic rate constant expression, such as the diabatic proton potential energy curves and associated proton vibrational wave functions, reorganization energy, and proton donor-acceptor distribution functions, are obtained from density functional theory calculations on model systems and molecular dynamics simulations of the solvated α 3 Y protein. Two possible proton acceptors, namely, water or a glutamate residue in the protein scaffold, are explored. The PCET rate constant is greater when glutamate is the proton acceptor, mainly due to the more favorable driving force and shorter equilibrium proton donor-acceptor distance, although contributions from excited vibronic states mitigate these effects. Nevertheless, water could be the dominant proton acceptor if its equilibrium constant associated with hydrogen bond formation is significantly greater than that for glutamate. Although these calculations do not definitively identify the proton acceptor for this PCET reaction, they elucidate the conditions under which each proton acceptor can be favored. These insights have implications for tyrosine-based PCET in a wide variety of biochemical processes.

The projected decline of available phosphorus necessitates alternative methods to derive usable phosphate for fertilizer and other applications. Phosphite dehydrogenase oxidizes phosphite to phosphate with the cofactor NAD + serving as the hydride acceptor. In addition to producing phosphate, this enzyme plays an important role in NADH cofactor regeneration processes. Mixed quantum mechanical/molecular mechanical free energy simulations were performed to elucidate the mechanism of this enzyme and to identify the protonation states of the substrate and product. Specifically, the finite temperature string method with umbrella sampling was used to generate the free energy surfaces and determine the minimum free energy paths for six different initial conditions that varied in the protonation state of the substrate and the position of the nucleophilic water molecule. In contrast to previous studies, the mechanism predicted by all six independent strings is a concerted but asynchronous dissociative mechanism in which hydride transfer from the phosphite substrate to NAD + occurs prior to attack by the nucleophilic water molecule. His292 is identified as the most likely general base that deprotonates the attacking water molecule. However, Arg237 could also serve as this base if it were deprotonated and His292 were protonated prior to the main chemical transformation, although this scenario is less probable. The simulations indicate that the phosphite substrate is monoanionic in its active form and that the most likely product is dihydrogen phosphate. These mechanistic insights may be helpful for designing mutant enzymes or artificial constructs that convert phosphite to phosphate and NAD + to NADH more effectively.

The recently discovered twister ribozyme is thought to utilize general acid-base catalysis in its self-cleavage mechanism, but the roles of nucleobases and metal ions in the mechanism are unclear. Herein, molecular dynamics simulations of the env22 twister ribozyme are performed to elucidate the structural and equilibrium dynamical properties, as well as to examine the role of Mg(2+) ions and possible candidates for the general base and acid in the self-cleavage mechanism. The active site region and the ends of the pseudoknots were found to be less mobile than other regions of the ribozyme, most likely providing structural stability and possibly facilitating catalysis. A purported catalytic Mg(2+) ion and the closest neighboring Mg(2+) ion remained chelated and relatively immobile throughout the microsecond trajectories, although removal of these Mg(2+) ions did not lead to any significant changes in the structure or equilibrium motions of the ribozyme on the microsecond time scale. In addition, a third metal ion, a Na(+) ion remained close to A1(O5'), the leaving group atom, during the majority of the microsecond trajectories, suggesting that it might stabilize the negative charge on A1(O5') during self-cleavage. The locations of these cations and their interactions with key nucleotides in the active site suggest that they may be catalytically relevant. The P1 stem is partially melted at its top and bottom in the crystal structure and further unwinds in the trajectories. The simulations also revealed an interconnected network comprised of hydrogen-bonding and π-stacking interactions that create a relatively rigid network around the self-cleavage site. The nucleotides involved in this network are among the highly conserved nucleotides in twister ribozymes, suggesting that this interaction network may be important to structure and function.

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates. The RP substrate gives large normal thio effects in the presence of all monovalent ions. The SP substrate also gives normal or no thio effects, but only for smaller monovalent and divalent cations, such as Li(+), Mg(2+), Ca(2+), and Sr(2+); in contrast, sizable inverse thio effects are found for larger monovalent and divalent cations, including Na(+), K(+), NH4(+), and Ba(2+). Proton inventories are found to be unity in the presence of the larger monovalent and divalent ions, but two in the presence of Mg(2+). Additionally, rate-pH profiles are inverted for the low charge density ions, and only imidazole plus ammonium ions rescue an inactive C75Δ variant in the absence of Mg(2+). Results from the thio effect experiments, rate-pH profiles, proton inventories, and ammonium/imidazole rescue experiments, combined with QM/MM free energy simulations, support a change in the mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized to stepwise and proton transfer-stabilized as the charge density of the metal ion decreases.

The positioning of catalytic groups within proteins plays an important role in enzyme catalysis, and here we investigate the positioning of the general base in the enzyme ketosteroid isomerase (KSI). The oxygen atoms of Asp38, the general base in KSI, were previously shown to be involved in anion-aromatic interactions with two neighboring Phe residues. Here we ask whether those interactions are sufficient, within the overall protein architecture, to position Asp38 for catalysis or whether the side chains that pack against Asp38 and/or the residues of the structured loop that is capped by Asp38 are necessary to achieve optimal positioning for catalysis. To test positioning, we mutated each of the aforementioned residues, alone and in combinations, in a background with the native Asp general base and in a D38E mutant background, as Glu at position 38 was previously shown to be mispositioned for general base catalysis. These double-mutant cycles reveal positioning effects as large as 10(3)-fold, indicating that structural features in addition to the overall protein architecture and the Phe residues neighboring the carboxylate oxygen atoms play roles in positioning. X-ray crystallography and molecular dynamics simulations suggest that the functional effects arise from both restricting dynamic fluctuations and disfavoring potential mispositioned states. Whereas it may have been anticipated that multiple interactions would be necessary for optimal general base positioning, the energetic contributions from positioning and the nonadditive nature of these interactions are not revealed by structural inspection and require functional dissection. Recognizing the extent, type, and energetic interconnectivity of interactions that contribute to positioning catalytic groups has implications for enzyme evolution and may help reveal the nature and extent of interactions required to design enzymes that rival those found in biology.

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

This brief review analyzes the underlying physical principles of enzyme catalysis, with an emphasis on the role of equilibrium enzyme motions and conformational sampling. The concepts are developed in the context of three representative systems, namely, dihydrofolate reductase, ketosteroid isomerase, and soybean lipoxygenase. All of these reactions involve hydrogen transfer, but many of the concepts discussed are more generally applicable. The factors that are analyzed in this review include hydrogen tunneling, proton donor-acceptor motion, hydrogen bonding, pKa shifting, electrostatics, preorganization, reorganization, and conformational motions. The rate constant for the chemical step is determined primarily by the free energy barrier, which is related to the probability of sampling configurations conducive to the chemical reaction. According to this perspective, stochastic thermal motions lead to equilibrium conformational changes in the enzyme and ligands that result in configurations favorable for the breaking and forming of chemical bonds. For proton, hydride, and proton-coupled electron transfer reactions, typically the donor and acceptor become closer to facilitate the transfer. The impact of mutations on the catalytic rate constants can be explained in terms of the factors enumerated above. In particular, distal mutations can alter the conformational motions of the enzyme and therefore the probability of sampling configurations conducive to the chemical reaction. Methods such as vibrational Stark spectroscopy, in which environmentally sensitive probes are introduced site-specifically into the enzyme, provide further insight into these aspects of enzyme catalysis through a combination of experiments and theoretical calculations.

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pK(a) shifted toward neutrality. Less is known about the role of metal ions in the mechanism. A recent crystal structure of the precleavage ribozyme identified a Mg²⁺ ion that interacts through its partial hydration sphere with the G25·U20 reverse wobble. In addition, this Mg²⁺ ion is in position to directly coordinate the nucleophile, the 2'-hydroxyl of U(-1), suggesting it can serve as a Lewis acid to facilitate deprotonation of the 2'-hydroxyl. To test the role of the active site Mg²⁺ ion, we replaced the G25·U20 reverse wobble with an isosteric A25·C20 reverse wobble. This change was found to significantly reduce the negative potential at the active site, as supported by electrostatics calculations, suggesting that active site Mg²⁺ binding could be adversely affected by the mutation. The kinetic analysis and molecular dynamics of the A25·C20 double mutant suggest that this variant stably folds into an active structure. However, pH-rate profiles of the double mutant in the presence of Mg²⁺ are inverted relative to the profiles for the wild-type ribozyme, suggesting that the A25·C20 double mutant has lost the active site metal ion. Overall, these studies support a model in which the partially hydrated Mg²⁺ positioned at the G25·U20 reverse wobble is catalytic and could serve as a Lewis acid, a Brønsted base, or both to facilitate deprotonation of the nucleophile.

This brief review discusses our current understanding of the molecular basis of enzyme catalysis. A historical development is presented, beginning with steady state kinetics and progressing through modern fast reaction methods, nuclear magnetic resonance, and single-molecule fluorescence techniques. Experimental results are summarized for ribonuclease, aspartate aminotransferase, and especially dihydrofolate reductase (DHFR). Multiple intermediates, multiple conformations, and cooperative conformational changes are shown to be an essential part of virtually all enzyme mechanisms. In the case of DHFR, theoretical investigations have provided detailed information about the movement of atoms within the enzyme-substrate complex as the reaction proceeds along the collective reaction coordinate for hydride transfer. A general mechanism is presented for enzyme catalysis that includes multiple intermediates and a complex, multidimensional standard free energy surface. Protein flexibility, diverse protein conformations, and cooperative conformational changes are important features of this model.

Classical molecular dynamics simulations were utilized to investigate the structural and dynamical properties of water in the active site of ketosteroid isomerase (KSI) to provide insight into the role of these water molecules in the enzyme-catalyzed reaction. This reaction is thought to proceed via a dienolate intermediate that is stabilized by hydrogen bonding with residues Tyr16 and Asp103. A comparative study was performed for the wild-type (WT) KSI and the Y16F, Y16S, and Y16F/Y32F/Y57F (FFF) mutants. These systems were studied with three different bound ligands: equilenin, which is an intermediate analog, and the intermediate states of two steroid substrates. Several distinct water occupation sites were identified in the active site of KSI for the WT and mutant systems. Three additional sites were identified in the Y16S mutant that were not occupied in WT KSI or the other mutants studied. The number of water molecules directly hydrogen bonded to the ligand oxygen was approximately two in the Y16S mutant and one in the Y16F and FFF mutants, with intermittent hydrogen bonding of one water molecule in WT KSI. The molecular dynamics trajectories of the Y16F and FFF mutants reproduced the small conformational changes of residue 16 observed in the crystal structures of these two mutants. Quantum mechanical/molecular mechanical calculations of (1)H NMR chemical shifts of the protons in the active site hydrogen-bonding network suggest that the presence of water in the active site does not prevent the formation of short hydrogen bonds with far-downfield chemical shifts. The molecular dynamics simulations indicate that the active site water molecules exchange much more frequently for WT KSI and the FFF mutant than for the Y16F and Y16S mutants. This difference is most likely due to the hydrogen-bonding interaction between Tyr57 and an active site water molecule that is persistent in the Y16F and Y16S mutants but absent in the FFF mutant and significantly less probable in WT KSI., (© 2011 American Chemical Society)

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Computational studies are performed to analyze the physical properties of hydrogen bonds donated by Tyr16 and Asp103 to a series of substituted phenolate inhibitors bound in the active site of ketosteroid isomerase (KSI). As the solution pK(a) of the phenolate increases, these hydrogen bond distances decrease, the associated nuclear magnetic resonance (NMR) chemical shifts increase, and the fraction of protonated inhibitor increases, in agreement with prior experiments. The quantum mechanical/molecular mechanical calculations provide insight into the electronic inductive effects along the hydrogen bonding network that includes Tyr16, Tyr57, and Tyr32, as well as insight into hydrogen bond coupling in the active site. The calculations predict that the most-downfield NMR chemical shift observed experimentally corresponds to the Tyr16-phenolate hydrogen bond and that Tyr16 is the proton donor when a bound naphtholate inhibitor is observed to be protonated in electronic absorption experiments. According to these calculations, the electronic inductive effects along the hydrogen bonding network of tyrosines cause the Tyr16 hydroxyl to be more acidic than the Asp103 carboxylic acid moiety, which is immersed in a relatively nonpolar environment. When one of the distal tyrosine residues in the network is mutated to phenylalanine, thereby diminishing this inductive effect, the Tyr16-phenolate hydrogen bond becomes longer and the Asp103-phenolate hydrogen bond shorter, as observed in NMR experiments. Furthermore, the calculations suggest that the differences in the experimental NMR data and electronic absorption spectra for pKSI and tKSI, two homologous bacterial forms of the enzyme, are due predominantly to the third tyrosine that is present in the hydrogen bonding network of pKSI but not tKSI. These studies also provide experimentally testable predictions about the impact of mutating the distal tyrosine residues in this hydrogen bonding network on the NMR chemical shifts and electronic absorption spectra.

Hybrid quantum/classical molecular dynamics simulations of the two proton transfer reactions catalyzed by ketosteroid isomerase are presented. The potential energy surfaces for the proton transfer reactions are described with the empirical valence bond method. Nuclear quantum effects of the transferring hydrogen increase the rates by a factor of approximately 8, and dynamical barrier recrossings decrease the rates by a factor of 3-4. For both proton transfer reactions, the donor-acceptor distance decreases substantially at the transition state. The carboxylate group of the Asp38 side chain, which serves as the proton acceptor and donor in the first and second steps, respectively, rotates significantly between the two proton transfer reactions. The hydrogen-bonding interactions within the active site are consistent with the hydrogen bonding of both Asp99 and Tyr14 to the substrate. The simulations suggest that a hydrogen bond between Asp99 and the substrate is present from the beginning of the first proton transfer step, whereas the hydrogen bond between Tyr14 and the substrate is virtually absent in the first part of this step but forms nearly concurrently with the formation of the transition state. Both hydrogen bonds are present throughout the second proton transfer step until partial dissociation of the product. The hydrogen bond between Tyr14 and Tyr55 is present throughout both proton transfer steps. The active site residues are more mobile during the first step than during the second step. The van der Waals interaction energy between the substrate and the enzyme remains virtually constant along the reaction pathway, but the electrostatic interaction energy is significantly stronger for the dienolate intermediate than for the reactant and product. Mobile loop regions distal to the active site exhibit significant structural rearrangements and, in some cases, qualitative changes in the electrostatic potential during the catalytic reaction. These results suggest that relatively small conformational changes of the enzyme active site and substrate strengthen the hydrogen bonds that stabilize the intermediate, thereby facilitating the proton transfer reactions. Moreover, the conformational and electrostatic changes associated with these reactions are not limited to the active site but rather extend throughout the entire enzyme.

The concept is developed that enzyme mechanisms should be viewed as "catalytic networks" with multiple conformations that occur serially and in parallel in the mechanism. These coupled ensembles of conformations require a multi-dimensional standard free-energy surface that is very "rugged", containing multiple minima and transition states. Experimental and theoretical evidence is presented to support this concept.

Experimental and theoretical data imply that enzyme motion plays an important role in enzymatic reactions. Enzyme motion can influence both the activation free energy barrier and the degree of barrier recrossing. A hybrid theoretical approach has been developed for the investigation of the relation between enzyme motion and activity. This approach includes both electronic and nuclear quantum effects. It distinguishes between thermally averaged promoting motions that influence the activation free energy barrier and dynamical motions that influence the barrier recrossings. Applications to hydride transfer in liver alcohol dehydrogenase and dihydrofolate reductase resulted in the identification and characterization of important enzyme motions. These applications have also led to the proposal of a network of coupled promoting motions in enzymatic reactions. These concepts have important implications for protein engineering and drug design.