Role of catalytic residues in enzymatic mechanisms of homologous ketosteroid isomerases

Bond Bundle Analysis of Ketosteroid Isomerase

Bond bundle analysis is used to investigate enzymatic catalysis in the ketosteroid isomerase (KSI) active site. We identify the unique bonding regions in five KSI systems, including those exposed to applied oriented electric fields and those with amino acid mutations, and calculate the precise redistribution of electron density and other regional properties that accompanies either enhancement or inhibition of KSI catalytic activity. We find that catalytic enhancement results from promoting both inter- and intra-molecular electron density redistribution, between bond bundles and bond wedges within the KSI-docked substrate molecule, in the forward direction of the catalyzed reaction. Though the redistribution applies to both types of perturbed systems and is thus suggestive of a general catalytic role, we observe that bond properties (e.g., volume vs energy vs electron count) can respond independently and disproportionately depending on the type of perturbation. We conclude that the resulting catalytic enhancement/inhibition proceeds via different mechanisms, where some bond properties are utilized more by one type of perturbation than the other. Additionally, we find that the correlations between bond wedge properties and catalyzed reaction barrier energies are additive to predict those of bond bundles and atomic basins, providing a rigorous grounding for connecting changes in local charge density to resulting shifts in reaction barrier energy.

Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations

The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.

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.

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.

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.

On catalytic preorganization in oxyanion holes: highlighting the problems with the gas-phase modeling of oxyanion holes and illustrating the need for complete enzyme models

Oxyanion holes play a major role in catalyzing enzymatic reactions, yet the corresponding energetics is frequently misunderstood. The main problem may be associated with the nontrivial nature of the electrostatic preorganization effect, without following the relevant formulation. That is, although the energetics of oxyanion holes have been fully quantified in early studies (which include both the enzymatic and reference solution reactions), the findings of these studies are sometimes overlooked, and, in some cases, it is assumed that gas-phase calculations with a fixed model of an oxyanion hole are sufficient for assessing the corresponding effect in the protein. Herein, we present a systematic analysis of this issue, clarifying the problems associated with modeling oxyanions by means of two fixed water molecules (or related constructs). We then re-emphasize the point that the effect of the oxyanion hole is mainly due to the fact that the relevant dipoles are already set in an orientation that stabilizes the TS charges, whereas the corresponding dipoles in solution are randomly oriented, resulting in the need to pay a very large reorganization energy. Simply calculating interaction energies with relatively fixed species cannot capture this crucial point, and considering it may help in advancing rational enzyme design.

Impact of mutation on proton transfer reactions in ketosteroid isomerase: insights from molecular dynamics simulations

The two proton transfer reactions catalyzed by ketosteroid isomerase (KSI) involve a dienolate intermediate stabilized by hydrogen bonds with Tyr14 and Asp99. Molecular dynamics simulations based on an empirical valence bond model are used to examine the impact of mutating these residues on the hydrogen-bonding patterns, conformational changes, and van der Waals and electrostatic interactions during the proton transfer reactions. While the rate constants for the two proton transfer steps are similar for wild-type (WT) KSI, the simulations suggest that the rate constant for the first proton transfer step is smaller in the mutants due to the significantly higher free energy of the dienolate intermediate relative to the reactant. The calculated rate constants for the mutants D99L, Y14F, and Y14F/D99L relative to WT KSI are qualitatively consistent with the kinetic experiments indicating a significant reduction in the catalytic rates along the series of mutants. In the simulations, WT KSI retained two hydrogen-bonding interactions between the substrate and the active site, while the mutants typically retained only one hydrogen-bonding interaction. A new hydrogen-bonding interaction between the substrate and Tyr55 was observed in the double mutant, leading to the prediction that mutation of Tyr55 will have a greater impact on the proton transfer rate constants for the double mutant than for WT KSI. The electrostatic stabilization of the dienolate intermediate relative to the reactant was greater for WT KSI than for the mutants, providing a qualitative explanation for the significantly reduced rates of the mutants. The active site exhibited restricted motion during the proton transfer reactions, but small conformational changes occurred to facilitate the proton transfer reactions by strengthening the hydrogen-bonding interactions and by bringing the proton donor and acceptor closer to each other with the proper orientation for proton transfer. Thus, these calculations suggest that KSI forms a preorganized active site but that the structure of this preorganized active site is altered upon mutation. Moreover, small conformational changes due to stochastic thermal motions are required within this preorganized active site to facilitate the proton transfer reactions.

Dissecting the paradoxical effects of hydrogen bond mutations in the ketosteroid isomerase oxyanion hole

The catalytic importance of enzyme active-site interactions is frequently assessed by mutating specific residues and measuring the resulting rate reductions. This approach has been used in bacterial ketosteroid isomerase to probe the energetic importance of active-site hydrogen bonds donated to the dienolate reaction intermediate. The conservative Tyr16Phe mutation impairs catalysis by 10(5)-fold, far larger than the effects of hydrogen bond mutations in other enzymes. However, the less-conservative Tyr16Ser mutation, which also perturbs the Tyr16 hydrogen bond, results in a less-severe 10(2)-fold rate reduction. To understand the paradoxical effects of these mutations and clarify the energetic importance of the Tyr16 hydrogen bond, we have determined the 1.6-A resolution x-ray structure of the intermediate analogue, equilenin, bound to the Tyr16Ser mutant and measured the rate effects of mutating Tyr16 to Ser, Thr, Ala, and Gly. The nearly identical 200-fold rate reductions of these mutations, together with the 6.4-A distance observed between the Ser16 hydroxyl and equilenin oxygens in the x-ray structure, strongly suggest that the more moderate rate effect of this mutant is not due to maintenance of a hydrogen bond from Ser at position 16. These results, additional spectroscopic observations, and prior structural studies suggest that the Tyr16Phe mutation results in unfavorable interactions with the dienolate intermediate beyond loss of a hydrogen bond, thereby exaggerating the apparent energetic benefit of the Tyr16 hydrogen bond relative to the solution reaction. These results underscore the complex energetics of hydrogen bonding interactions and site-directed mutagenesis experiments.

Hybrid quantum/classical molecular dynamics simulations of the proton transfer reactions catalyzed by ketosteroid isomerase: analysis of hydrogen bonding, conformational motions, and electrostatics

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.

Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole

Enzymes are classically proposed to accelerate reactions by binding substrates within active-site environments that are structurally preorganized to optimize binding interactions with reaction transition states rather than ground states. This is a remarkably formidable task considering the limited 0.1-1 A scale of most substrate rearrangements. The flexibility of active-site functional groups along the coordinate of substrate rearrangement, the distance scale on which enzymes can distinguish structural rearrangement, and the energetic significance of discrimination on that scale remain open questions that are fundamental to a basic physical understanding of enzyme active sites and catalysis. We bring together 1.2-1.5 A resolution X-ray crystallography, (1)H and (19)F NMR spectroscopy, quantum mechanical calculations, and transition-state analogue binding measurements to test the distance scale on which noncovalent forces can constrain the structural relaxation or translation of side chains and ligands along a specific coordinate and the energetic consequences of such geometric constraints within the active site of bacterial ketosteroid isomerase (KSI). Our results strongly suggest that packing and binding interactions within the KSI active site can constrain local side-chain reorientation and prevent hydrogen bond shortening by 0.1 A or less. Further, this constraint has substantial energetic effects on ligand binding and stabilization of negative charge within the oxyanion hole. These results provide evidence that subtle geometric effects, indistinguishable in most X-ray crystallographic structures, can have significant energetic consequences and highlight the importance of using synergistic experimental approaches to dissect enzyme function.

De novo design approach based on nanorecognition toward development of functional molecules/materials and nanosensors/nanodevices

For the design of functional molecules and nanodevices, it is very useful to utilize nanorecognition (which is governed mainly by interaction forces such as hydrogen bonding, ionic interaction, π-H/π-π interactions, and metallic interactions) and nanodynamics (involving capture, transport, and release of electrons, photons, or protons). The manifestation of these interaction forces has led us to the design and realization of diverse ionophores/receptors, organic nanotubes, nanowires, molecular mechanical devices, molecular switches, enzyme mimetics, protein folding/unfolding, etc. In this review, we begin with a brief discussion of the interaction forces, followed by some of our representative applications. We discuss ionophores with chemo-sensing capability for biologically important cations and anions and explain how the understanding of hydrogen bonding and π-interactions has led to the design of self-assembled nanotubes from calix[4]hydroquinone (CHQ). The binding study of neutral and cationic transition metals with the redox system of hydroquinone (HQ) and quinone (Q) predicts what kind of nanostructures would form. Finally, we look into the conformational changes between stacked and edge-to-face conformers in π-benzoquinone-benzene complexes controlled by alternating electrochemical potential. The resulting flapping motion illustrates a promising pathway toward the design of mobile nanomechanical devices.

Dissolution nature of the lithium hydroxide by water molecules

The structures, stabilities, thermodynamic quantities, dissociation energies, infrared spectra, and electronic properties of LiOH hydrated by up to seven water molecules are investigated by using the density-functional theory and the Møller-Plesset second-order perturbation theory (MP2). Further accurate analysis based on the coupled-cluster theory with singles, doubles, and perturbative triples excitations agrees with the MP2 results. The Li-OH stretch mode significantly shifts with the increase of water molecules, and it eventually disappears upon dissociation. It is revealed that seven water molecules are needed for the stable dissociation of LiOH (as a completely dissociated conformation), in contrast to the cases of RbOH and CsOH which require four and three water molecules, respectively.

Enzymatic mechanisms for catalysis of enolization: ketosteroid isomerase

Breaking a carbon-hydrogen bond adjacent to a carbonyl is a slow step in a large number of chemical reactions. However, many enzymes are capable of catalyzing this reaction with great efficiency. One of the most proficient of these enzymes is 3-oxo-Delta5-steroid isomerase (KSI), which catalyzes the isomerization of a wide variety of 3-oxo-Delta5-steroids to their Delta4-conjugated isomers. In this review, the mechanism of KSI is discussed, with particular emphasis on energetic considerations. Both experimental and theoretical approaches are considered to explain the mechanistic details of the reaction.

Theoretical Studies on the Mechanism of Acid-Promoted Hydrolysis of N-Formylaziridine in Comparison with Formamide

[reaction: see text] We present an ab initio study of the acid-promoted hydrolysis reaction mechanism of N-formylaziridine in comparison with formamide. Since the rate of amide hydrolysis reactions depends on the formation of the tetrahedral intermediate, we focused our attention mainly on the reactant complex, the tetrahedral intermediate, and the transition state connecting these two stationary points. Geometries were optimized using the density functional theory, and the energetics were refined using ab initio theory including electron correlation. Solvent effects were investigated by using polarizable continuum method calculations. The proton-transfer reaction between the O-protonated and N-protonated amides was investigated. In acidic media, despite that the N-protonated species is more stable than the O-protonated one, it is predicted that both N-protonated and O-protonated pathways compete in the hydrolysis reaction of N-formylaziridine.

Theoretical insights into the mechanism of acetylcholinesterase-catalyzed acylation of acetylcholine

Acylation of acetylcholine (ACh) catalyzed by acetylcholinesterase (AChE) has been studied using high-level theoretical calculations on a model system that mimics the reaction center of the enzyme, and compared with uncatalyzed acylation reaction. The geometries of all the intermediates and transition states, activation energies, and solvent effects have been calculated. The calculations predict simultaneous formation of two short-strong hydrogen bonds (SSHB) in the rate-determining transition state structures [the first SSHB involves the hydrogen atom of Ser-200 (H(s)) and another involves the hydrogen atom of His-440 (H(h))]. In the intermediate states, the H-bond corresponding to H(h) involves SSHB, whereas the one corresponding to H(s) does not.

Structures, energetics, and spectra of hydrated hydroxide anion clusters

The structures, energetics, electronic properties, and spectra of hydrated hydroxide anions are studied using density functional and high level ab initio calculations. The overall structures and binding energies are similar to the hydrated anion clusters, in particular, to the hydrated fluoride anion clusters except for the tetrahydrated clusters and hexahydrated clusters. In tetrahydrated system, tricoordinated structures and tetracoordinated structures are compatible, while in pentahydrated systems and hexahydrated systems, tetracoordinated structures are stable. The hexahydrated system is similar in structure to the hydrated chloride cluster. The thermodynamic quantities (enthalpies and free energies) of the clusters are in good agreement with the experimental values. The electronic properties induced by hydration are similar to hydrated chloride anions. The charge-transfer-to-solvent energies of these hydrated-hydroxide anions are discussed, and the predicted ir spectra are used to explain the experimental data in terms of the cluster structures. The low-energy barriers between the conformations along potential energy surfaces are reported.

Efficient Fluoride-Selective Fluorescent Host: Experiment and Theory

A new naphthalene derivative containing a urea group at the 1,8-position of naphthalene was synthesized and showed a unique absorption and fluorescence peak with fluoride ion. Calculations suggested that a new peak was attributed to the increased anionic character of urea nitrogen due to the strongly charged hydrogen bonding between fluoride and amide protons of the urea. The fluoride selectivity among halides (F(-), Cl(-), Br(-)) comes from the fact that the fluoride approaches much closer to the amide protons than other halides and resides in the cavity with fast dynamics. The nature of electronic transitions that were analyzed from the calculations by the collective electronic oscillator method also supports the anionic nature of the complex between host and fluoride.

Many chemistries could be used to build living systems

It has been widely suggested that life based around carbon, hydrogen, oxygen, and nitrogen is the only plausible biochemistry, and specifically that terrestrial biochemistry of nucleic acids, proteins, and sugars is likely to be "universal." This is not an inevitable conclusion from our knowledge of chemistry. I argue that it is the nature of the liquid in which life evolves that defines the most appropriate chemistry. Fluids other than water could be abundant on a cosmic scale and could therefore be an environment in which non-terrestrial biochemistry could evolve. The chemical nature of these liquids could lead to quite different biochemistries, a hypothesis discussed in the context of the proposed "ammonochemistry" of the internal oceans of the Galilean satellites and a more speculative "silicon biochemistry" in liquid nitrogen. These different chemistries satisfy the thermodynamic drive for life through different mechanisms, and so will have different chemical signatures than terrestrial biochemistry.

Molecular recognition of fluoride anion: benzene-based tripodal imidazolium receptor

A benzene-based tripodal imidazolium receptor utilizing the strong (C-H)(+)...X(-) hydrogen bonding interaction between imidazolium moieties and halide anions is extensively investigated both theoretically and experimentally. Ab initio calculations predict that this receptor has a very high affinity for fluoride ion (F(-)). The association constant and free energy gain of the N-butyl receptor 2 for F(-) in acetonitrile were measured to be 2.1 x 10(5) M(-1) and -7.25 kcal/mol, respectively, showing that the receptor has a high affinity for F(-) in highly polar organic solvents.

Tripodal nitro-imidazolium receptor for anion binding driven by (C-H)+- - -X- hydrogen bonds

[structure: see text] A positively charged tripodal receptor with nitro groups in the imidazolium rings was designed, synthesized, and characterized for its anion binding strength. The receptor shows strong affinity and high selectivity for Cl- through (C-H)+- - -X(-) hydrogen bonds wherein charge-charge and charge-dipole electrostatic interactions dominate. The association constant with chloride anion in a 9:1 mixture of acetonitrile-d3 and DMSO-d6 is measured to be 1.1 x 10(6) M(-1). The receptor also shows reasonably high affinity toward H2PO4-.