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

Effects and Influence of External Electric Fields on the Equilibrium Properties of Tautomeric Molecules

In this review, we have attempted to briefly summarize the influence of an external electric field on an assembly of tautomeric molecules and to what experimentally observable effects this interaction can lead to. We have focused more extensively on the influence of an oriented external electric field (OEEF) on excited-state intramolecular proton transfer (ESIPT) from the studies available to date. The possibilities provided by OEEF for regulating several processes and studying physicochemical processes in tautomers have turned this direction into an attractive area of research due to its numerous applications.

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.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

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.

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.

Measuring electrostatic fields in both hydrogen-bonding and non-hydrogen-bonding environments using carbonyl vibrational probes

Vibrational probes can provide a direct readout of the local electrostatic field in complex molecular environments, such as protein binding sites and enzyme active sites. This information provides an experimental method to explore the underlying physical causes of important biomolecular processes such as binding and catalysis. However, specific chemical interactions such as hydrogen bonds can have complicated effects on vibrational probes and confound simple electrostatic interpretations of their frequency shifts. We employ vibrational Stark spectroscopy along with infrared spectroscopy of carbonyl probes in different solvent environments and in ribonuclease S to understand the sensitivity of carbonyl frequencies to electrostatic fields, including those due to hydrogen bonds. Additionally, we carried out molecular dynamics simulations to calculate ensemble-averaged electric fields in solvents and in ribonuclease S and found excellent correlation between calculated fields and vibrational frequencies. These data enabled the construction of a robust field-frequency calibration curve for the C═O vibration. The present results suggest that carbonyl probes are capable of quantitatively assessing the electrostatics of hydrogen bonding, making them promising for future study of protein function.

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.