Active Site Of Ketosteroid Isomerase Research Articles

An electrostatic energy-based charge model for molecular dynamics simulation.

The interactions of the polar chemical bonds such as C=O and N-H with an external electric field were investigated, and a linear relationship between the QM/MM interaction energies and the electric field along the chemical bond is established in the range of weak to intermediate electrical fields. The linear relationship indicates that the electrostatic interactions of a polar group with its surroundings can be described by a simple model of a dipole with constant moment under the action of an electric field. This relationship is employed to develop a general approach to generating an electrostatic energy-based charge (EEC) model for molecules containing single or multiple polar chemical bonds. Benchmark test studies of this model were carried out for (CH3)2-CO and N-methyl acetamide in explicit water, and the result shows that the EEC model gives more accurate electrostatic energies than those given by the widely used charge model based on fitting to the electrostatic potential (ESP) in direct comparison to the energies computed by the QM/MM method. The MD simulations of the electric field at the active site of ketosteroid isomerase based on EEC demonstrated that EEC gave a better representation of the electrostatic interaction in the hydrogen-bonding environment than the Amber14SB force field by comparison with experiment. The current study suggests that EEC should be better suited for molecular dynamics study of molecular systems with polar chemical bonds such as biomolecules than the widely used ESP or RESP (restrained ESP) charge models.

Charge Density in Enzyme Active Site as a Descriptor of Electrostatic Preorganization

Large protein macromolecules in enzymatic catalysis have been shown to exert a specific electric field that reduces the reorganization energy upon barrier crossing and thus reduces the reaction free energy barrier. In this work we suggest that the charge density in the active site of an enzyme investigated using formalisms embodied by the quantum theory of atoms in molecules (QTAIM) provides a sensitive and quantum mechanically rigorous probe of electrostatic preorganization. We focus on the active site of ketosteroid isomerase, a well-studied enzyme for which electrostatic preorganization has been modeled theoretically and studied experimentally. We study the charge density in the active site and the reaction mechanism in the presence of small external electric fields of various directions and magnitudes. We show that the geometry of the full charge density is a sensitive reporter on the external field experienced by the active site. Changes are observed in the relative positions of critical points and amount of charge at critical points as a function of the field. At the same time, a subset of these features correlates linearly with the barrier of the first reaction step in catalysis. Small changes in the barrier, within 1-2 kcal/mol, are reflected in the charge density, suggesting the existence of a field - reactant state charge density - reaction barrier correlation. Hence, QTAIM can be used for the analysis of electric field in enzyme active sites, and further investigations and exploitations of the found correlations may prove useful in enzyme design where preorganization is optimized.

An Ab Initio QM/MM Study of the Electrostatic Contribution to Catalysis in the Active Site of Ketosteroid Isomerase.

The electric field in the hydrogen-bond network of the active site of ketosteroid isomerase (KSI) has been experimentally measured using vibrational Stark effect (VSE) spectroscopy, and utilized to study the electrostatic contribution to catalysis. A large gap was found in the electric field between the computational simulation based on the Amber force field and the experimental measurement. In this work, quantum mechanical (QM) calculations of the electric field were performed using an ab initio QM/MM molecular dynamics (MD) simulation and electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method. Our results demonstrate that the QM-derived electric field based on the snapshots from QM/MM MD simulation could give quantitative agreement with the experiment. The accurate calculation of the electric field inside the protein requires both the rigorous sampling of configurations, and a QM description of the electrostatic field. Based on the direct QM calculation of the electric field, we theoretically confirmed that there is a linear correlation relationship between the activation free energy and the electric field in the active site of wild-type KSI and its mutants (namely, D103N, Y16S, and D103L). Our study presents a computational protocol for the accurate simulation of the electric field in the active site of the protein, and provides a theoretical foundation that supports the link between electric fields and enzyme catalysis.

Correction to Uncovering the Determinants of a Highly Perturbed Tyrosine p Ka in the Active Site of Ketosteroid Isomerase.

Correction to Uncovering the Determinants of a Highly Perturbed Tyrosine p Ka in the Active Site of Ketosteroid Isomerase.

Proton Network Flexibility Enables Robustness and Large Electric Fields in the Ketosteroid Isomerase Active Site.

Hydrogen-bond networks play vital roles in biological functions ranging from protein folding to enzyme catalysis. Here we combine electronic structure calculations and ab initio path integral molecular dynamics simulations, which incorporate both nuclear and electronic quantum effects, to show why the network of short hydrogen bonds in the active site of ketosteroid isomerase is remarkably robust to mutations along the network and how this gives rise to large local electric fields. We demonstrate that these properties arise from the network's ability to respond to a perturbation by shifting proton positions and redistributing electronic charge density. This flexibility leads to small changes in properties such as the partial ionization of residues and pKa isotope effects upon mutation of the residues, consistent with recent experiments. This proton flexibility is further enhanced when an extended hydrogen-bond network forms in the presence of an intermediate analogue, which allows us to explain the chemical origins of the large electric fields in the enzyme's active site observed in recent experiments.

An Activator-Blocker Pair Provides a Controllable On-Off Switch for a Ketosteroid Isomerase Active Site Mutant.

Control of enzyme activity is fundamental to biology and represents a long-term goal in bioengineering and precision therapeutics. While several powerful molecular strategies have been developed, limitations remain in their generalizability and dynamic range. We demonstrate a control mechanism via separate small molecules that turn on the enzyme (activator) and turn off the activation (blocker). We show that a pocket created near the active site base of the enzyme ketosteriod isomerase (KSI) allows efficient and saturable base rescue when the enzyme's natural general base is removed. Binding a small molecule with similar properties but lacking general-base capability in this pocket shuts off rescue. The ability of small molecules to directly participate in and directly block catalysis may afford a broad controllable dynamic range. This approach may be amenable to numerous enzymes and to engineering and screening approaches to identify activators and blockers with strong, specific binding for engineering and therapeutic applications.

Kemp Eliminase Activity of Ketosteroid Isomerase.

Kemp eliminases represent the most successful class of computationally designed enzymes, with rate accelerations of up to 109-fold relative to the rate of the same reaction in aqueous solution. Nevertheless, several other systems such as micelles, catalytic antibodies, and cavitands are known to accelerate the Kemp elimination by several orders of magnitude. We found that the naturally occurring enzyme ketosteroid isomerase (KSI) also catalyzes the Kemp elimination. Surprisingly, mutations of D38, the residue that acts as a general base for its natural substrate, produced variants that catalyze the Kemp elimination up to 7000-fold better than wild-type KSI does, and some of these variants accelerate the Kemp elimination more than the computationally designed Kemp eliminases. Analysis of the D38N general base KSI variant suggests that a different active site carboxylate residue, D99, performs the proton abstraction. Docking simulations and analysis of inhibition by active site binders suggest that the Kemp elimination takes place in the active site of KSI and that KSI uses the same catalytic strategies of the computationally designed enzymes. In agreement with prior observations, our results strengthen the conclusion that significant rate accelerations of the Kemp elimination can be achieved with very few, nonspecific interactions with the substrate if a suitable catalytic base is present in a hydrophobic environment. Computational design can fulfill these requirements, and the design of more complex and precise environments represents the next level of challenges for protein design.

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.

Dissecting Proton Delocalization in an Enzyme's Hydrogen Bond Network with Unnatural Amino Acids.

Extended hydrogen bond networks are a common structural motif of enzymes. A recent analysis proposed quantum delocalization of protons as a feature present in the hydrogen bond network spanning a triad of tyrosines (Y(16), Y(32), and Y(57)) in the active site of ketosteroid isomerase (KSI), contributing to its unusual acidity and large isotope shift. In this study, we utilized amber suppression to substitute each tyrosine residue with 3-chlorotyrosine to test the delocalization model and the proton affinity balance in the triad. X-ray crystal structures of each variant demonstrated that the structure, notably the O-O distances within the triad, was unaffected by 3-chlorotyrosine substitutions. The changes in the cluster's acidity and the acidity's isotope dependence in these variants were assessed via UV-vis spectroscopy and the proton sharing pattern among individual residues with (13)C nuclear magnetic resonance. Our data show pKa detuning at each triad residue alters the proton delocalization behavior in the H-bond network. The extra stabilization energy necessary for the unusual acidity mainly comes from the strong interactions between Y(57) and Y(16). This is further enabled by Y(32), which maintains the right geometry and matched proton affinity in the triad. This study provides a rich picture of the energetics of the hydrogen bond network in enzymes for further model refinement.

BIOPHYSICS. Response to Comments on "Extreme electric fields power catalysis in the active site of ketosteroid isomerase".

Natarajan et al. and Chen and Savidge comment that comparing the electric field in ketosteroid isomerase's (KSI's) active site to zero overestimates the catalytic effect of KSI's electric field because the reference reaction occurs in water, which itself exerts a sizable electrostatic field. To compensate, Natarajan et al. argue that additional catalytic weight arises from positioning of the general base, whereas Chen and Savidge propose a separate contribution from desolvation of the general base. We note that the former claim is not well supported by published results, and the latter claim is intriguing but lacks experimental basis. We also take the opportunity to clarify some of the more conceptually subtle aspects of electrostatic catalysis.

BIOPHYSICS. Comment on "Extreme electric fields power catalysis in the active site of ketosteroid isomerase".

Fried et al. (Reports, 19 December 2014, p. 1510) demonstrate electric field-dependent acceleration of biological catalysis using ketosteroid isomerase as a prototypic example. These findings were not extended to aqueous solution because water by itself has field fluctuations that are too large and fast to provide a catalytic effect. Given physiological context, when water electrostatic interactions are considered, electric fields play a less important role in the catalysis.

BIOPHYSICS. Comment on "Extreme electric fields power catalysis in the active site of ketosteroid isomerase".

Fried et al. (Reports, 19 December 2014, p. 1510) demonstrated a strong correlation between reaction rate and the carbonyl stretching frequency of a product analog bound to ketosteroid isomerase oxyanion hole mutants and concluded that the active-site electric field provides 70% of catalysis. Alternative comparisons suggest a smaller contribution, relative to the corresponding solution reaction, and highlight the importance of atomic-level descriptions.

Extreme electric fields power catalysis in the active site of ketosteroid isomerase.

Enzymes use protein architecture to impose specific electrostatic fields onto their bound substrates, but the magnitude and catalytic effect of these electric fields have proven difficult to quantify with standard experimental approaches. Using vibrational Stark effect spectroscopy, we found that the active site of the enzyme ketosteroid isomerase (KSI) exerts an extremely large electric field onto the C=O chemical bond that undergoes a charge rearrangement in KSI's rate-determining step. Moreover, we found that the magnitude of the electric field exerted by the active site strongly correlates with the enzyme's catalytic rate enhancement, enabling us to quantify the fraction of the catalytic effect that is electrostatic in origin. The measurements described here may help explain the role of electrostatics in many other enzymes and biomolecular systems.

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.

Quantum Delocalization of Protons in the Ketosteroid Isomerase Active Site

Ketosteroid isomerase (KSI) catalyzes steroid isomerization with extremely high efficiency and has become a paradigm of enzymatic proton transfer chemistry. In this poster I will show our recent research which has identified how the hydrogen bond network formed by this enzyme facilitate proton delocalization and sharing in the active site. This quantum delocalization greatly stabilizes a deprotonated tyrosine residue in the active site and leads to a 10,000 fold increase in its acid dissociation constant. This study is made possible by a series of technical breakthroughs developed in our group that greatly enhance the efficiency of including quantum mechanical fluctuations in simulating condensed phase systems and extracting isotope effects. These advances allow us to investigate the quantum delocalization of protons in the active site of KSI, which might shed light on the origin of the enzyme's remarkable catalytic efficiency and guide the development of synthetic catalysts.

Uncovering the Determinants of a Highly Perturbed Tyrosine pKain the Active Site of Ketosteroid Isomerase

Within the idiosyncratic enzyme active-site environment, side chain and ligand pKa values can be profoundly perturbed relative to their values in aqueous solution. Whereas structural inspection of systems has often attributed perturbed pKa values to dominant contributions from placement near charged groups or within hydrophobic pockets, Tyr57 of a Pseudomonas putida ketosteroid isomerase (KSI) mutant, suggested to have a pKa perturbed by nearly 4 units to 6.3, is situated within a solvent-exposed active site devoid of cationic side chains, metal ions, or cofactors. Extensive comparisons among 45 variants with mutations in and around the KSI active site, along with protein semisynthesis, (13)C NMR spectroscopy, absorbance spectroscopy, and X-ray crystallography, was used to unravel the basis for this perturbed Tyr pKa. The results suggest that the origin of large energetic perturbations are more complex than suggested by visual inspection. For example, the introduction of positively charged residues near Tyr57 raises its pKa rather than lowers it; this effect, and part of the increase in the Tyr pKa from the introduction of nearby anionic groups, arises from accompanying active-site structural rearrangements. Other mutations with large effects also cause structural perturbations or appear to displace a structured water molecule that is part of a stabilizing hydrogen-bond network. Our results lead to a model in which three hydrogen bonds are donated to the stabilized ionized Tyr, with these hydrogen-bond donors, two Tyr side chains, and a water molecule positioned by other side chains and by a water-mediated hydrogen-bond network. These results support the notion that large energetic effects are often the consequence of multiple stabilizing interactions rather than a single dominant interaction. Most generally, this work provides a case study for how extensive and comprehensive comparisons via site-directed mutagenesis in a tight feedback loop with structural analysis can greatly facilitate our understanding of enzyme active-site energetics. The extensive data set provided may also be a valuable resource for those wishing to extensively test computational approaches for determining enzymatic pKa values and energetic effects.

Use of anion–aromatic interactions to position the general base in the ketosteroid isomerase active site

Although the cation-pi pair, formed between a side chain or substrate cation and the negative electrostatic potential of a pi system on the face of an aromatic ring, has been widely discussed and has been shown to be important in protein structure and protein-ligand interactions, there has been little discussion of the potential structural and functional importance in proteins of the related anion-aromatic pair (i.e., interaction of a negatively charged group with the positive electrostatic potential on the ring edge of an aromatic group). We posited, based on prior structural information, that anion-aromatic interactions between the anionic Asp general base and Phe54 and Phe116 might be used instead of a hydrogen-bond network to position the general base in the active site of ketosteroid isomerase from Comamonas testosteroni as there are no neighboring hydrogen-bonding groups. We have tested the role of the Phe residues using site-directed mutagenesis, double-mutant cycles, and high-resolution X-ray crystallography. These results indicate a catalytic role of these Phe residues. Extensive analysis of the Protein Data Bank provides strong support for a catalytic role of these and other Phe residues in providing anion-aromatic interactions that position anionic general bases within enzyme active sites. Our results further reveal a potential selective advantage of Phe in certain situations, relative to more traditional hydrogen-bonding groups, because it can simultaneously aid in the binding of hydrophobic substrates and positioning of a neighboring general base.

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.

Ground State Destabilization from a Positioned General Base in the Ketosteroid Isomerase Active Site

We compared the binding affinities of ground state analogues for bacterial ketosteroid isomerase (KSI) with a wild-type anionic Asp general base and with uncharged Asn and Ala in the general base position to provide a measure of potential ground state destabilization that could arise from the close juxtaposition of the anionic Asp and hydrophobic steroid in the reaction's Michaelis complex. The analogue binding affinity increased ~1 order of magnitude for the Asp38Asn mutation and ~2 orders of magnitude for the Asp38Ala mutation, relative to the affinity with Asp38, for KSI from two sources. The increased level of binding suggests that the abutment of a charged general base and a hydrophobic steroid is modestly destabilizing, relative to a standard state in water, and that this destabilization is relieved in the transition state and intermediate in which the charge on the general base has been neutralized because of proton abstraction. Stronger binding also arose from mutation of Pro39, the residue adjacent to the Asp general base, consistent with an ability of the Asp general base to now reorient to avoid the destabilizing interaction. Consistent with this model, the Pro mutants reduced or eliminated the increased level of binding upon replacement of Asp38 with Asn or Ala. These results, supported by additional structural observations, suggest that ground state destabilization from the negatively charged Asp38 general base provides a modest contribution to KSI catalysis. They also provide a clear illustration of the well-recognized concept that enzymes evolve for catalytic function and not, in general, to maximize ground state binding. This ground state destabilization mechanism may be common to the many enzymes with anionic side chains that deprotonate carbon acids.

Vibrational Stark Effects in the Active Site of Ketosteroid Isomerase Point to Large Electric Fields Driving Chemical Catalysis

Enzymes are extraordinary catalysts that actuate nearly all biomolecular processes with speed and specificity. Nevertheless, the physical origins of enzymes’ catalytic power remain elusive despite investigations of many enzymes’ mechanisms over the last half-century. Ketosteroid isomerase (KSI) is a small and proficient enzyme - accelerating an isomerization reaction ∼1 trillion-fold over its intrinsic rate in water - that has been employed as a test system to examine the catalytic strategies at Nature's disposal. Electrostatic interactions are broadly purported to play an essential role in catalysis, but this proposal has yet to be experimentally tested in a quantitative fashion. Toward this end, vibrational Stark effect spectroscopy provides a toolkit to measure electric fields in biomolecular environments. By measuring the vibrational Stark effect of an intrinsic probe that mimics the reactive species of KSI's catalytic cycle, we found that KSI's active site focuses an extremely large electric field onto the scissile bond, potentially enabling its speedy chemical conversion through electrostatic interactions. Moreover, we observed a strong correlation between active site electric field and catalytic power across several KSI mutants. These studies are building toward a highly reductionist picture of KSI's catalysis and possibly enzyme catalysis in general.