Proton affinity of the oxyanion hole in the active site of ketosteroid isomerase

Negative influence of pKa on activation energy barrier: A case study for double proton transfer reaction in inorganic acid dimers

Strength of acid can be determined by means of pK_a value. Attempts have been made to find a relationship between pK_a and activation energy barrier for a double proton transfer (DPT) reaction in inorganic acid dimers. Negative influence of pK_a is observed on activation energy (E_a ) which is contrary to the general convention of pK_a . Four different levels of theories with two different basis sets have been used to calculate the activation energy barrier of the DPT reaction in inorganic acid dimers. A model based on first and second order polynomial has been created to find the relationship between activation energy for DPT reaction. © 2018 Wiley Periodicals, Inc.

De novo endo-functionalized organic cages as cooperative multi-hydrogen-bond-donating catalysts

Two endo-functionalized organic cages as oxyanion hole mimics were achieved via dynamic covalent chemistry, which exhibit good size selectivity, catalytic activity and broad substrate scopes for Friedel-Crafts reactions. The modular character of the synthesis provides an easy way to modify the size, shape and inner function of the molecular cavity.

A Critical Test of the Electrostatic Contribution to Catalysis with Noncanonical Amino Acids in Ketosteroid Isomerase

The vibrational Stark effect (VSE) has been used to measure the electric field in the active site of ketosteroid isomerase (KSI). These measured fields correlate with ΔG(⧧) in a series of conventional mutants, yielding an estimate for the electrostatic contribution to catalysis (Fried et al. Science 2014, 346, 1510-1513). In this work we test this result with much more conservative variants in which individual Tyr residues in the active site are replaced by 3-chlorotyrosine via amber suppression. The electric fields sensed at the position of the carbonyl bond involved in charge displacement during catalysis were characterized using the VSE, where the field sensitivity has been calibrated by vibrational Stark spectroscopy, solvatochromism, and MD simulations. A linear relationship is observed between the electric field and ΔG(⧧) that interpolates between wild-type and more drastic conventional mutations, reinforcing the evaluation of the electrostatic contribution to catalysis in KSI. A simplified model and calculation are developed to estimate changes in the electric field accompanying changes in the extended hydrogen-bond network in the active site. The results are consistent with a model in which the O-H group of a key active site tyrosine functions by imposing a static electrostatic potential onto the carbonyl bond. The model suggests that the contribution to catalysis from the active site hydrogen bonds is of similar weight to the distal interactions from the rest of the protein. A similar linear correlation was also observed between the proton affinity of KSI's active site and the catalytic rate, suggesting a direct connection between the strength of the H-bond and the electric field it exerts.

Short Hydrogen Bonds and Proton Delocalization in Green Fluorescent Protein (GFP)

Short hydrogen bonds and specifically low-barrier hydrogen bonds (LBHBs) have been the focus of much attention and controversy for their possible role in enzymatic catalysis. The green fluorescent protein (GFP) mutant S65T, H148D has been found to form a very short hydrogen bond between Asp148 and the chromophore resulting in significant spectral perturbations. Leveraging the unique autocatalytically formed chromophore and its sensitivity to this interaction we explore the consequences of proton affinity matching across this putative LBHB. Through the use of noncanonical amino acids introduced through nonsense suppression or global incorporation, we systematically modify the acidity of the GFP chromophore with halogen substituents. X-ray crystal structures indicated that the length of the interaction with Asp148 is unchanged at ∼2.45 Å while the absorbance spectra demonstrate an unprecedented degree of color tuning with increasing acidity. We utilized spectral isotope effects, isotope fractionation factors, and a simple 1D model of the hydrogen bond coordinate in order to gain insight into the potential energy surface and particularly the role that proton delocalization may play in this putative short hydrogen bond. The data and model suggest that even with the short donor-acceptor distance (∼2.45 Å) and near perfect affinity matching there is not a LBHB, that is, the barrier to proton transfer exceeds the H zero-point energy.

The hydrogen bond strength of the phenol–phenolate anionic complex: a computational and photoelectron spectroscopic study

The phenol–phenolate anionic complex was studied in vacuo by negative ion photoelectron spectroscopy using 193 nm photons and by density functional theory (DFT) computations at the ωB97XD/6-311+G(2d,p) level.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site

Hydrogen bonds are ubiquitous in enzyme active sites, providing binding interactions and stabilizing charge rearrangements on substrate groups over the course of a reaction. But understanding the origin and magnitude of their catalytic contributions relative to hydrogen bonds made in aqueous solution remains difficult, in part because of complexities encountered in energetic interpretation of traditional site-directed mutagenesis experiments. It has been proposed for ketosteroid isomerase and other enzymes that active site hydrogen bonding groups provide energetic stabilization via "short, strong" or "low-barrier" hydrogen bonds that are formed due to matching of their pKa or proton affinity to that of the transition state. It has also been proposed that the ketosteroid isomerase and other enzyme active sites provide electrostatic environments that result in larger energetic responses (i.e., greater "sensitivity") to ground-state to transition-state charge rearrangement, relative to aqueous solution, thereby providing catalysis relative to the corresponding reaction in water. To test these models, we substituted tyrosine with fluorotyrosines (F-Tyr's) in the ketosteroid isomerase (KSI) oxyanion hole to systematically vary the proton affinity of an active site hydrogen bond donor while minimizing steric or structural effects. We found that a 40-fold increase in intrinsic F-Tyr acidity caused no significant change in activity for reactions with three different substrates. F-Tyr substitution did not change the solvent or primary kinetic isotope effect for proton abstraction, consistent with no change in mechanism arising from these substitutions. The observed shallow dependence of activity on the pKa of the substituted Tyr residues suggests that the KSI oxyanion hole does not provide catalysis by forming an energetically exceptional pKa-matched hydrogen bond. In addition, the shallow dependence provides no indication of an active site electrostatic environment that greatly enhances the energetic response to charge accumulation, consistent with prior experimental results.

Thermodynamic framework for identifying free energy inventories of enzyme catalytic cycles

Pauling's suggestion that enzymes are complementary in structure to the activated complexes of the reactions they catalyze has provided the conceptual basis to explain how enzymes obtain their fantastic catalytic prowess, and has served as a guiding principle in drug design for over 50 y. However, this model by itself fails to predict the magnitude of enzymes' rate accelerations. We construct a thermodynamic framework that begins with the classic concept of differential binding but invokes additional terms that are needed to account for subtle effects in the catalytic cycle's proton inventory. Although the model presented can be applied generally, this analysis focuses on ketosteroid isomerase (KSI) as an example, where recent experiments along with a large body of kinetic and thermodynamic data have provided strong support for the noncanonical thermodynamic contribution described. The resulting analysis precisely predicts the free energy barrier of KSI's reaction as determined from transition-state theory using only empirical thermodynamic data. This agreement is suggestive that a complete free energy inventory of the KSI catalytic cycle has been identified.

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.

Calculation of vibrational shifts of nitrile probes in the active site of ketosteroid isomerase upon ligand binding

The vibrational Stark effect provides insight into the roles of hydrogen bonding, electrostatics, and conformational motions in enzyme catalysis. In a recent application of this approach to the enzyme ketosteroid isomerase (KSI), thiocyanate probes were introduced in site-specific positions throughout the active site. This paper implements a quantum mechanical/molecular mechanical (QM/MM) approach for calculating the vibrational shifts of nitrile (CN) probes in proteins. This methodology is shown to reproduce the experimentally measured vibrational shifts upon binding of the intermediate analogue equilinen to KSI for two different nitrile probe positions. Analysis of the molecular dynamics simulations provides atomistic insight into the roles that key residues play in determining the electrostatic environment and hydrogen-bonding interactions experienced by the nitrile probe. For the M116C-CN probe, equilinen binding reorients an active-site water molecule that is directly hydrogen-bonded to the nitrile probe, resulting in a more linear C≡N--H angle and increasing the CN frequency upon binding. For the F86C-CN probe, equilinen binding orients the Asp103 residue, decreasing the hydrogen-bonding distance between the Asp103 backbone and the nitrile probe and slightly increasing the CN frequency. This QM/MM methodology is applicable to a wide range of biological systems and has the potential to assist in the elucidation of the fundamental principles underlying enzyme catalysis.

Three hydrogen bond donor catalysts: oxyanion hole mimics and transition state analogues

Enzymes and their mimics use hydrogen bonds to catalyze chemical transformations. Small-molecule transition state analogues of oxyanion holes have been characterized by computations, gas-phase IR and photoelectron spectroscopy, and determination of their binding constants in acetonitrile. A new class of hydrogen bond catalysts is proposed (donors that can contribute three hydrogen bonds to a single functional group) and demonstrated in a Friedel-Crafts reaction. The employed catalyst was observed to react 100 times faster than its rotamer that can employ only two hydrogen bonds. The former compound also binds anions more tightly and was found to have a thermodynamic advantage.

Site-specific measurement of water dynamics in the substrate pocket of ketosteroid isomerase using time-resolved vibrational spectroscopy

Little is known about the reorganization capacity of water molecules at the active sites of enzymes and how this couples to the catalytic reaction. Here, we study the dynamics of water molecules at the active site of a highly proficient enzyme, Δ(5)-3-ketosteroid isomerase (KSI), during a light-activated mimic of its catalytic cycle. Photoexcitation of a nitrile-containing photoacid, coumarin183 (C183), mimics the change in charge density that occurs at the active site of KSI during the first step of the catalytic reaction. The nitrile of C183 is exposed to water when bound to the KSI active site, and we used time-resolved vibrational spectroscopy as a site-specific probe to study the solvation dynamics of water molecules in the vicinity of the nitrile. We observed that water molecules at the active site of KSI are highly rigid, during the light-activated catalytic cycle, compared to the solvation dynamics observed in bulk water. On the basis of this result, we hypothesize that rigid water dipoles at the active site might help in the maintenance of the preorganized electrostatic environment required for efficient catalysis. The results also demonstrate the utility of nitrile probes in measuring the dynamics of local (H-bonded) water molecules in contrast to the commonly used fluorescence methods which measure the average behavior of primary and subsequent spheres of solvation.

Evaluation of the energetics of the concerted acid-base mechanism in enzymatic catalysis: the case of ketosteroid isomerase

Structures of enzymes invariably reveal the proximity of acidic and basic residues to reactive sites on the substrate, so it is natural and common to suggest that enzymes employ concerted mechanisms to catalyze their difficult reactions. Ketosteroid isomerase (KSI) has served as a paradigm of enzymatic proton transfer chemistry, and its catalytic effect has previously been attributed to concerted proton transfer. We employ a specific inhibitor that contains an IR probe that reports directly and quantitatively on the ionization state of the ligand when bound in the active site of KSI. Measurement of the fractional ionization provides a missing link in a thermodynamic cycle that can discriminate the free energy advantage of a concerted versus nonconcerted mechanism. It is found that the maximum thermodynamic advantage that KSI could capture from a concerted mechanism (ΔΔG° = 0.5 kcal mol(-1)) is quite small.

Direct measurement of the protein response to an electrostatic perturbation that mimics the catalytic cycle in ketosteroid isomerase

Understanding how electric fields and their fluctuations in the active site of enzymes affect efficient catalysis represents a critical objective of biochemical research. We have directly measured the dynamics of the electric field in the active site of a highly proficient enzyme, Δ(5)-3-ketosteroid isomerase (KSI), in response to a sudden electrostatic perturbation that simulates the charge displacement that occurs along the KSI catalytic reaction coordinate. Photoexcitation of a fluorescent analog (coumarin 183) of the reaction intermediate mimics the change in charge distribution that occurs between the reactant and intermediate state in the steroid substrate of KSI. We measured the electrostatic response and angular dynamics of four probe dipoles in the enzyme active site by monitoring the time-resolved changes in the vibrational absorbance (IR) spectrum of a spectator thiocyanate moiety (a quantitative sensor of changes in electric field) placed at four different locations in and around the active site, using polarization-dependent transient vibrational Stark spectroscopy. The four different dipoles in the active site remain immobile and do not align to the changes in the substrate electric field. These results indicate that the active site of KSI is preorganized with respect to functionally relevant changes in electric fields.

Stabilization of the second oxyanion intermediate by 1,4-dihydroxy-2-naphthoyl-coenzyme A synthase of the menaquinone pathway: spectroscopic evidence of the involvement of a conserved aspartic acid

1,4-Dihydroxy-2-naphthoyl-coenzyme A (DHNA-CoA) synthase, or MenB, catalyzes an intramolecular Claisen condensation involving two oxyanion intermediates in the biosynthetic pathway of menaquinone, an essential respiration electron transporter in many microorganisms. Here we report the finding that the DHNA-CoA product and its analogues bind and inhibit the synthase from Escherichia coli with significant ultraviolet--visible spectral changes, which are similar to the changes induced by deprotonation of the free inhibitors in a basic solution. Dissection of the structure--affinity relationships of the inhibitors identifies the hydroxyl groups at positions 1 (C1-OH) and 4 (C4-OH) of DHNA-CoA or their equivalents as the dominant and minor sites, respectively, for the enzyme--ligand interaction that polarizes or deprotonates the bound ligands to cause the observed spectral changes. In the meantime, spectroscopic studies with active site mutants indicate that C4-OH of the enzyme-bound DHNA-CoA interacts with conserved polar residues Arg-91, Tyr-97, and Tyr-258 likely through a hydrogen bonding network that also includes Ser-161. In addition, site-directed mutation of the conserved Asp-163 to alanine causes a complete loss of the ligand binding ability of the protein, suggesting that the Asp-163 side chain is most likely hydrogen-bonded to C1-OH of DHNA-CoA to provide the dominant polarizing effect. Moreover, this mutation also completely eliminates the enzyme activity, strongly supporting the possibility that the Asp-163 side chain provides a strong stabilizing hydrogen bond to the tetrahedral oxyanion, which takes a position similar to that of C1-OH of the enzyme-bound DHNA-CoA and is the second high-energy intermediate in the intracellular Claisen condensation reaction. Interestingly, both Arg-91 and Tyr-97 are located in a disordered loop forming part of the active site of all available DHNA-CoA synthase structures. Their involvement in the interaction with the small molecule ligands suggests that the disordered loop is folded in interaction with the substrates or reaction intermediates, supporting an induced-fit catalytic mechanism for the enzyme.

Hydrogen bonding in the active site of ketosteroid isomerase: electronic inductive effects and hydrogen bond coupling

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.

Solvation response along the reaction coordinate in the active site of ketosteroid isomerase

A light-activated reaction analog has been developed to mimic the catalytic reaction cycle of Delta(5)-3-ketosteroid isomerase to probe the functionally relevant protein solvation response to the catalytic charge transfer. Delta(5)-3-ketosteroid isomerase from Pseudomonas putida catalyzes a C-H bond cleavage and formation through an enolate intermediate. Conversion of the ketone substrate to the enolate intermediate is simulated by a photoacid bound to the active site oxyanion hole. In the ground state, the photoacid electrostatically resembles the enolate intermediate while the low pK(a) excited state resembles the ketone starting material. Time-resolved fluorescence experiments with photoacids coumarin 183 and equilenin show the active site of Delta(5)-3-ketosteroid isomerase to be largely unperturbed by the light-activated reaction. The small solvation response for the photoacid at the active site as compared with a simple solvent suggests the active site does not significantly change its electrostatic environment during the catalytic cycle. Instead, the reaction takes place in an electrostatically preorganized environment.