Jason P. Schwans

Structural inference of native and partially folded RNA by high-throughput contact mapping

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (·OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4–P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4–P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a compact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes.

Quantitative, directional measurement of electric field heterogeneity in the active site of ketosteroid isomerase

Understanding the electrostatic forces and features within highly heterogeneous, anisotropic, and chemically complex enzyme active sites and their connection to biological catalysis remains a longstanding challenge, in part due to the paucity of incisive experimental probes of electrostatic properties within proteins. To quantitatively assess the landscape of electrostatic fields at discrete locations and orientations within an enzyme active site, we have incorporated site-specific thiocyanate vibrational probes into multiple positions within bacterial ketosteroid isomerase. A battery of X-ray crystallographic, vibrational Stark spectroscopy, and NMR studies revealed electrostatic field heterogeneity of 8 MV/cm between active site probe locations and widely differing sensitivities of discrete probes to common electrostatic perturbations from mutation, ligand binding, and pH changes. Electrostatic calculations based on active site ionization states assigned by literature precedent and computational pKa prediction were unable to quantitatively account for the observed vibrational band shifts. However, electrostatic models of the D40N mutant gave qualitative agreement with the observed vibrational effects when an unusual ionization of an active site tyrosine with a pKa near 7 was included. UV-absorbance and 13C NMR experiments confirmed the presence of a tyrosinate in the active site, in agreement with electrostatic models. This work provides the most direct measure of the heterogeneous and anisotropic nature of the electrostatic environment within an enzyme active site, and these measurements provide incisive benchmarks for further developing accurate computational models and a foundation for future tests of electrostatics in enzymatic catalysis.

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.

Determining the catalytic role of remote substrate binding interactions in ketosteroid isomerase

A fundamental difference between enzymes and small chemical catalysts is the ability of enzymes to use binding interactions with nonreactive portions of substrates to accelerate chemical reactions. Remote binding interactions can localize substrates to the active site, position substrates relative to enzymatic functional groups and other substrates, trigger conformational changes, induce local destabilization, and modulate an active site environment by solvent exclusion. We investigated the role of remote substrate binding interactions in the reaction catalyzed by the enzyme ketosteroid isomerase (KSI), which catalyzes a double bond migration of steroid substrates through a dienolate intermediate that is stabilized in an oxyanion hole. Comparison of a single-ring and multiple-ring substrate allowed the catalytic contribution of binding interactions with the distal substrate rings to be determined. The value of kcat/KM for a single-ring substrate is reduced 27,000-fold relative to a multiple-ring steroid substrate, suggesting that remote binding interactions with the steroid substrate contribute substantially to the KSI reaction. Nevertheless, the reaction rates for KSI-bound single- and multiple-ring substrates (kcat) are within 2-fold. Further, oxyanion hole mutations have the same effect on reactions of the single- and multiple-ring substrates. These results suggest that remote binding interactions contribute >5 kcal/mol to catalysis by KSI but that local rather than remote interactions dictate the catalytic contributions from KSIs general base and oxyanion hole.

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

Significance Hydrogen bond networks play critical structural and functional roles in proteins but have been challenging to study within this complex environment. We incorporated spectroscopic probes into the active site of the bacterial enzyme ketosteroid isomerase to systematically dissect the proton transfer equilibrium within a key hydrogen bond network formed to bound transition state analogs. Our study provides direct insight into the physical and energetic properties of a hydrogen bond network within an enzyme and presents a simple computational model of electrostatic effects within this protein that succeeds due to detailed knowledge of ionization states and a tightly controlled experimental system. 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.

Uncovering the determinants of a highly perturbed tyrosine pKa in 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.

Evaluating the Catalytic Contribution from the Oxyanion Hole in Ketosteroid Isomerase

Prior site-directed mutagenesis studies in bacterial ketosteroid isomerase (KSI) reported that substitution of both oxyanion hole hydrogen bond donors gives a 10(5)- to 10(8)-fold rate reduction, suggesting that the oxyanion hole may provide the major contribution to KSI catalysis. But these seemingly conservative mutations replaced the oxyanion hole hydrogen bond donors with hydrophobic side chains that could lead to suboptimal solvation of the incipient oxyanion in the mutants, thereby potentially exaggerating the apparent energetic benefit of the hydrogen bonds relative to water-mediated hydrogen bonds in solution. We determined the functional and structural consequences of substituting the oxyanion hole hydrogen bond donors and several residues surrounding the oxyanion hole with smaller residues in an attempt to create a local site that would provide interactions more analogous to those in aqueous solution. These more drastic mutations created an active-site cavity estimated to be ~650 Å(3) and sufficient for occupancy by 15-17 water molecules and led to a rate decrease of only ~10(3)-fold for KSI from two different species, a much smaller effect than that observed from more traditional conservative mutations. The results underscore the strong context dependence of hydrogen bond energetics and suggest that the oxyanion hole provides an important, but moderate, catalytic contribution relative to the interactions in the corresponding solution reaction.

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 reactions 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.

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.

Experimental and computational mutagenesis to investigate the positioning of a general base within an enzyme active site.

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