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

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

Extreme Catalytic Power of Ketosteroid Isomerase Related to the Reversal of Proton Dislocations in Hydrogen-Bond Network

Dynamic electrostatic catalytic field (DECF) vectors derived from transition state and reactant wavefunctions for the two-step reaction occurring within ketosteroid isomerase (KSI) have been calculated using MP2/aug-cc-pVTZ and lower theory levels to determine the magnitude of the catalytic effect and the optimal directions of proton transfers in the KSI hydrogen-bond network. The most surprising and meaningful finding is that the KSI catalytic activity is enhanced by proton dislocations proceeding in opposite directions for each of the two consecutive reaction steps in the same hydrogen network. Such a mechanism allows an ultrafast switching of the catalytic proton wire environment, possibly related to the exceptionally high KSI catalytic power.

Influence of local microenvironment on the double hydrogen transfer in porphycene

We performed time-resolved transient absorption and fluorescence anisotropy measurements in order to study tautomerization of porphycene in rigid polymer matrices at cryogenic temperatures. Studies were carried out in poly(methyl methacrylate) (PMMA), poly(vinyl butyral) (PVB), and poly(vinyl alcohol) (PVA). The results prove that in all studied media hydrogen tunnelling plays a significant role in the double hydrogen transfer which becomes very sensitive to properties of the environment below approx. 150 K. We also demonstrate that there exist two populations of porphycene molecules in rigid media: "hydrogen-transferring" molecules, in which tautomerization occurs on time scales below 1 ns and "frozen" molecules in which double hydrogen transfer is too slow to be monitored with nanosecond techniques. The number of "frozen" molecules increases when the sample is cooled. We explain this effect by interactions of guest molecules with a rigid host matrix which disturbs symmetry of porphycene and hinders tunnelling. Temperature dependence of the number of hydrogen-transferring molecules suggests that the factor which restores the symmetry of the double-minimum potential well in porphycene are intermolecular vibrations localized in separated regions of the amorphous polymer.

Characterization of Cercospora nicotianae Hypothetical Proteins in Cercosporin Resistance

The photoactivated toxin, cercosporin, produced by Cercospora species, plays an important role in pathogenesis of this fungus to host plants. Cercosporin has almost universal toxicity to cells due to its production of reactive oxygen species including singlet oxygen. For that reason, Cercospora species, which are highly resistant to their own toxin, are good candidates to identify genes for resistance to cercosporin and to the reactive oxygen species it produces. In previous research, the zinc cluster transcription factor CRG1 (cercosporin resistance gene 1) was found to be crucial for Cercospora species' resistance against cercosporin, and subtractive hybridization analysis identified 185 genes differentially expressed between Cercospora nicotianae wild type (wt) and a crg1 mutant. The focus of this work was to identify and characterize the hypothetical proteins that were identified in the Cercospora nicotianae subtractive library as potential resistance factors. Quantitative RT-PCR analysis of the 20 genes encoding hypothetical proteins showed that two, 24cF and 71cR, were induced under conditions of cercosporin toxicity, suggesting a role in resistance. Transformation and expression of 24cF and 71cR in the cercosporin-sensitive fungus, Neurospora crassa, showed that 71cR provided increased resistance to cercosporin toxicity, whereas no significant increase was observed in 24cF transformants. Gene disruption was used to generate C. nicotianae 71cR mutants; these mutants did not differ from wt C. nicotianae in cercosporin resistance or production. Quantitative RT-PCR analysis showed induction of other resistance genes in the 71cR mutant that may compensate for the loss of 71cR. Analysis of 71cR conserved domains and secondary and tertiary structure identify the protein as having an NTF2-like superfamily DUF1348 domain with unknown function, to be intracellular and localized in the cytosol, and to have similarities to proteins in the steroid delta-isomerase family.

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.

Challenges in computational studies of enzyme structure, function and dynamics

In this review we give an overview of the field of Computational enzymology. We start by describing the birth of the field, with emphasis on the work of the 2013 chemistry Nobel Laureates. We then present key features of the state-of-the-art in the field, showing what theory, accompanied by experiments, has taught us so far about enzymes. We also briefly describe computational methods, such as quantum mechanics-molecular mechanics approaches, reaction coordinate treatment, and free energy simulation approaches. We finalize by discussing open questions and challenges.

Studying allosteric regulation in metal sensor proteins using computational methods

In this chapter, we describe advances made in understanding the mechanism of allosteric regulation of DNA operator binding in the ArsR/SmtB family of metal-sensing proteins using computational methods. The paradigm, zinc-sensing transcriptional repressor Staphylococcus aureus CzrA represents an excellent model system to understand how metal sensor proteins maintain cellular metal homeostasis. Here, we discuss studies that helped to characterize a metal ion-mediated hydrogen-bonding pathway (HBP) that plays a dominant role in the allosteric mechanism of DNA operator binding in these proteins. The chapter discusses computational methods used to provide a molecular basis for the large conformational motions and allosteric coupling free energy (~6kcal/mol) associated with Zn(II) binding in CzrA. We present an accurate and convenient means by which to include metal ions in the nuclear magnetic resonance (NMR) structure determination process using molecular dynamics (MD) constrained by NMR-derived data. The method provides a realistic and physically viable description of the metal-binding site(s) and has potentially broad applicability in the structure determination of metal ion-bound proteins, protein folding, and metal template protein-design studies. Finally, our simulations provide strong support for a proposed HBP that physically connects the metal-binding residue, His97, to the DNA-binding interface through the αR helix that is present only in the Zn(II)-bound state. We find the interprotomer hydrogen bond interaction to be significantly stronger (~8kcal/mol) at functional allosteric metal-binding sites compared to the apo proteins. This interaction works to overcome the considerable disorder at these hydrogen-bonding sites in apo protein and functions as a "switch" to lock in a weak DNA-binding conformation once metal is bound. This interaction is found to be considerably weaker in nonresponsive metal-binding sites. These findings suggest a conserved functional role of metal-mediated second-shell coordination hydrogen bonds at allosterically responsive sites in zinc-sensing transcription regulators.

Characterizing solution surface loop conformational flexibility of the GM2 activator protein

GM2AP has a β-cup topology with numerous X-ray structures showing multiple conformations for some of the surface loops, revealing conformational flexibility that may be related to function, where function is defined as either membrane binding associated with ligand binding and extraction or interaction with other proteins. Here, site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and molecular dynamic (MD) simulations are used to characterize the mobility and conformational flexibility of various structural regions of GM2AP. A series of 10 single cysteine amino acid substitutions were generated, and the constructs were chemically modified with the methanethiosulfonate spin label. Continuous wave (CW) EPR line shapes were obtained and subsequently simulated using the microscopic order macroscopic disorder (MOMD) program. Line shapes for sites that have multiple conformations in the X-ray structures required two spectral components, whereas spectra of the remaining sites were adequately fit with single-component parameters. For spin labeled sites L126C and I66C, spectra were acquired as a function of temperature, and simulations provided for the determination of thermodynamic parameters associated with conformational change. Binding to GM2 ligand did not alter the conformational flexibility of the loops, as evaluated by EPR and NMR spectroscopies. These results confirm that the conformational flexibility observed in the surface loops of GM2AP crystals is present in solution and that the exchange is slow on the EPR time scale (>ns). Furthermore, MD simulation results are presented and agree well with the conformational heterogeneity revealed by SDSL.

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 10³-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.

Solution NMR refinement of a metal ion bound protein using metal ion inclusive restrained molecular dynamics methods

Correctly calculating the structure of metal coordination sites in a protein during the process of nuclear magnetic resonance (NMR) structure determination and refinement continues to be a challenging task. In this study, we present an accurate and convenient means by which to include metal ions in the NMR structure determination process using molecular dynamics (MD) simulations constrained by NMR-derived data to obtain a realistic and physically viable description of the metal binding site(s). This method provides the framework to accurately portray the metal ions and its binding residues in a pseudo-bond or dummy-cation like approach, and is validated by quantum mechanical/molecular mechanical (QM/MM) MD calculations constrained by NMR-derived data. To illustrate this approach, we refine the zinc coordination complex structure of the zinc sensing transcriptional repressor protein Staphylococcus aureus CzrA, generating over 130 ns of MD and QM/MM MD NMR-data compliant sampling. In addition to refining the first coordination shell structure of the Zn(II) ion, this protocol benefits from being performed in a periodically replicated solvation environment including long-range electrostatics. We determine that unrestrained (not based on NMR data) MD simulations correlated to the NMR data in a time-averaged ensemble. The accurate solution structure ensemble of the metal-bound protein accurately describes the role of conformational sampling in allosteric regulation of DNA binding by zinc and serves to validate our previous unrestrained MD simulations of CzrA. This methodology has potentially broad applicability in the structure determination of metal ion bound proteins, protein folding and metal template protein-design studies.

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.

Catalytic efficiency of enzymes: a theoretical analysis

This brief review analyzes the underlying physical principles of enzyme catalysis, with an emphasis on the role of equilibrium enzyme motions and conformational sampling. The concepts are developed in the context of three representative systems, namely, dihydrofolate reductase, ketosteroid isomerase, and soybean lipoxygenase. All of these reactions involve hydrogen transfer, but many of the concepts discussed are more generally applicable. The factors that are analyzed in this review include hydrogen tunneling, proton donor-acceptor motion, hydrogen bonding, pKa shifting, electrostatics, preorganization, reorganization, and conformational motions. The rate constant for the chemical step is determined primarily by the free energy barrier, which is related to the probability of sampling configurations conducive to the chemical reaction. According to this perspective, stochastic thermal motions lead to equilibrium conformational changes in the enzyme and ligands that result in configurations favorable for the breaking and forming of chemical bonds. For proton, hydride, and proton-coupled electron transfer reactions, typically the donor and acceptor become closer to facilitate the transfer. The impact of mutations on the catalytic rate constants can be explained in terms of the factors enumerated above. In particular, distal mutations can alter the conformational motions of the enzyme and therefore the probability of sampling configurations conducive to the chemical reaction. Methods such as vibrational Stark spectroscopy, in which environmentally sensitive probes are introduced site-specifically into the enzyme, provide further insight into these aspects of enzyme catalysis through a combination of experiments and theoretical calculations.

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 pK(a) 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 pK(a) near 7 was included. UV-absorbance and (13)C 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.

Insight into the cation-π interaction at the metal binding site of the copper metallochaperone CusF

The periplasmic Cu(+)/Ag(+) chaperone CusF features a novel cation-π interaction between a Cu(+)/Ag(+) ion and Trp44 at the metal binding site. The nature and strength of the Cu(+)/Ag(+)-Trp44 interactions were investigated using computational methodologies. Quantum-mechanical (QM) calculations showed that the Cu(+) and Ag(+) interactions with Trp44 are of similar strength (~14 kcal/mol) and bond order. Quantum-mechanical/molecular-mechanical (QM/MM) calculations showed that Cu(+) binds in a distorted tetrahedral coordination environment in the Trp44Met mutant, which lacks the cation-π interaction. Molecular dynamics (MD) simulations of CusF in the apo and Cu(+)-bound states emphasized the importance of the Cu(+)-Trp44 interaction in protecting Cu(+) from water oxidation. The protein structure does not change over the time scale of hundreds of nanoseconds in the metal-bound state. The metal recognition site exhibits small motions in the apo state but remains largely preorganized toward metal binding. Trp44 remains oriented to form the cation-π interaction in the apo state and faces an energetic penalty to move away from the metal ion. Cu(+) binding quenches the protein's internal motions in regions linked to binding CusB, suggesting that protein motions play an essential role in Cu(+) transfer to CusB.

Reengineering rate-limiting, millisecond enzyme motions by introduction of an unnatural amino acid

Rate-limiting millisecond motions in wild-type (WT) Ribonuclease A (RNase A) are modulated by histidine 48. Here, we incorporate an unnatural amino acid, thia-methylimidazole, at this site (H48C-4MI) to investigate the effects of a single residue on protein motions over multiple timescales and on enzyme catalytic turnover. Molecular dynamics simulations reveal that H48C-4MI retains some crucial WT-like hydrogen bonding interactions but the extent of protein-wide correlated motions in the nanosecond regime is decreased relative to WT. NMR Carr-Purcell-Meiboom-Gill relaxation dispersion experiments demonstrate that millisecond conformational motions in H48C-4MI are present over a similar pH range compared to WT. Furthermore, incorporation of this nonnatural amino acid allows retention of WT-like catalytic activity over the full pH range. These studies demonstrate that the complexity of the protein energy landscape during the catalytic cycle can be maintained using unnatural amino acids, which may prove useful in enzyme design efforts.

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.

Water in the active site of ketosteroid isomerase

Classical molecular dynamics simulations were utilized to investigate the structural and dynamical properties of water in the active site of ketosteroid isomerase (KSI) to provide insight into the role of these water molecules in the enzyme-catalyzed reaction. This reaction is thought to proceed via a dienolate intermediate that is stabilized by hydrogen bonding with residues Tyr16 and Asp103. A comparative study was performed for the wild-type (WT) KSI and the Y16F, Y16S, and Y16F/Y32F/Y57F (FFF) mutants. These systems were studied with three different bound ligands: equilenin, which is an intermediate analog, and the intermediate states of two steroid substrates. Several distinct water occupation sites were identified in the active site of KSI for the WT and mutant systems. Three additional sites were identified in the Y16S mutant that were not occupied in WT KSI or the other mutants studied. The number of water molecules directly hydrogen bonded to the ligand oxygen was approximately two in the Y16S mutant and one in the Y16F and FFF mutants, with intermittent hydrogen bonding of one water molecule in WT KSI. The molecular dynamics trajectories of the Y16F and FFF mutants reproduced the small conformational changes of residue 16 observed in the crystal structures of these two mutants. Quantum mechanical/molecular mechanical calculations of (1)H NMR chemical shifts of the protons in the active site hydrogen-bonding network suggest that the presence of water in the active site does not prevent the formation of short hydrogen bonds with far-downfield chemical shifts. The molecular dynamics simulations indicate that the active site water molecules exchange much more frequently for WT KSI and the FFF mutant than for the Y16F and Y16S mutants. This difference is most likely due to the hydrogen-bonding interaction between Tyr57 and an active site water molecule that is persistent in the Y16F and Y16S mutants but absent in the FFF mutant and significantly less probable in WT KSI.

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.