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.     

Phosphoryl Transfer Enzymes and Hypervalent Phosphorus Chemistry: A Keynote Lecture of the 7th International Conference on Heteroatom Chemistry

The coordination tendencies of phosphorus to form a hexacoordinated state from a pentacoordinated state which might assist in describing mechanistic action of phosphoryl transfer enzymes are delineated. The factors discussed include substrate and transition or intermediate state anionicity; hydrogen bonding; packing effects, i.e., van der Waals forces; the ease of formation of hexacoordinate phosphorus from lower coordinate states; and the pseudorotation problem common to nonrigid pentacoordinate phosphorus. In view of the work reported in this account and recent work on enzyme promiscuity and moonlighting activities, it is suggested that donor action should play a role in determining active-site interactions in phosphoryl transfer enzyme mechanisms.     

Structure-based redesign of GST A2-2 for enhanced catalytic efficiency with azathioprine

Glutathione transferase (GST) A2-2 is the most efficient human enzyme in the biotransformation of the prodrug azathioprine (Aza). The activation of Aza has therapeutic potential for possible use of GSTs in targeted enzyme-prodrug treatment of diseases. Based on the assumed catalytic mechanism and computational docking of Aza to the active site of the enzyme, active-site residues were selected for construction of focused mutant libraries, which were thereafter screened for Aza activity. Mutants with elevated Aza activity were identified, DNA sequenced, and the proteins purified. The two most active mutants showed up to 70-fold higher catalytic efficiency than the parental GST A2-2. The structure of the most active triple mutant (L107G/L108D/F222H) enzyme was determined by X-ray crystallography demonstrating significant changes in the topography of the active site facilitating productive binding of Aza as a substrate.     

Effects of slow substrate binding and release in redox enzymes: theory and application to periplasmic nitrate reductase

For redox enzymes, the technique called protein film voltammetry makes it possible to determine the entire profile of activity against driving force by having the enzyme exchanging directly electrons with the rotating-disc electrode onto which it is adsorbed. Both the potential location of the catalytic response and its detailed shape report on the sequence of catalytic events, electron transfers and chemical steps, but the models that have been used so far to decipher this signal lack generality. For example, it was often proposed that substrate binding to multiple redox states of the active site may explain that turnover is greater in a certain window of electrode potential, but no fully analytical treatment has been given. Here, we derive (i) the general current equation for the case of reversible substrate binding to any redox states of a two-electron active site (as exemplified by flavins and Mo cofactors), (ii) the quantitative conditions for an extremum in activity to occur, and (iii) the expressions from which the substrate-concentration dependence of the catalytic potential can be interpreted to learn about the kinetics of substrate binding and how this affects the reduction potential of the active site. Not only does slow substrate binding and release make the catalytic wave shape highly complex, but we also show that it can have important consequences which will escape detection in traditional experiments: the position of the wave (this is the driving force that is required to elicit catalysis) departs from the reduction potential of the active site even at the lowest substrate concentration, and this deviation may be large if substrate binding is irreversible. This occurs in the reductive half-cycle of periplasmic nitrate reductase where irreversibility lowers the driving force required to reduce the active site under turnover conditions and favors intramolecular electron transfer from the proximal [4Fe4S]+ cluster to the active site Mo(V).     

Enzyme flexibility and the catalytic mechanism of farnesyltransferase: targeting the relation

Farnesyltransferase enzyme (FTase) is an interesting target for anticancer therapy that has been the subject of particular attention over the past decade. However, despite of the thrilling achievements in the development of farnesyltransferase inhibitors (FTIs) over the past few years, the farnesylation mechanism remains, to some degree, a mystery. This work describes the application of molecular dynamics simulations to the study of enzyme flexibility in the 4 key intermediate states formed during the FTase catalytic mechanism--FTase resting state, binary complex (FTase-FPP), ternary complex (FTase-FPP-Peptide), and product complex (FTase-Product)--thereby covering the main states in the mechanistic pathway of this mysterious enzyme, while relating, dissecting, and exploring, in minute detail, the set of structural and dynamical changes taking place with FPP binding, peptide coordination and product formation. This study reveals the existence of a series of variational patterns involving the mechanistic events taking place at the active site of the enzyme, increasing in magnitude away from the active-site, demonstrating that relatively small-scale events such as substrate binding or product formation cause minor changes at the neighboring residues and corresponding helices, but ultimately induce much more dramatic effects on the more external regions of the enzyme.     

A designed synthetic analogue of 4-OT is specific for a non-natural substrate

The substrate specificity of 4-oxalocrotonate tautomerase (4-OT) is characterized by electrostatic interactions between positively charged arginine (Arg) side chains on the enzyme and the dianionic substrate, 4-oxalocrotonate. To generate specific hydrogen-bonding interactions with a monoanionic substrate analogue, we have introduced a urea functional group into the active site by replacing arginine side chains with isosteric citrulline (Cit) residues. This design was based on the complementarity between the urea functionality of citrulline and the uncharged amide function of the substrate, as opposed to the guanidinium-carboxylate electrostatic interaction between the wild-type enzyme and the natural substrate. Indeed, the synthetic (Arg39Cit)4-OT analogue catalyzed the tautomerization of the non-natural monoamide-monoacid substrate while it was a poor catalyst for the natural diacid substrate. The specificity of (Arg39Cit)4-OT for the monoamide-monoacid substrate relative to that of the diacid substrate was found to be 740-fold greater than that of the wild-type enzyme for tautomerization of the non-natural substrate as compared with the natural one. The role of electrostatic interactions in the tautomerization of the monoamide-monoacid substrate was probed in detail with several other Arg to Cit analogues of this enzyme. This study has demonstrated that chemical manipulation of the functional groups within the active site of an enzyme can modify its catalytic activity and substrate specificity in a predictable way, suggesting that the incorporation of noncoded amino acids into proteins has great promise for the development of new enzymatic mechanisms and new binding interactions.     

Understanding the mechanism of action of B12-dependent ethanolamine ammonia-lyase: synergistic interactions at play

Ab initio molecular orbital calculations are used to examine the mechanism of action of B(12)-dependent ethanolamine ammonia-lyase involving the conversion of 2-aminoethanol to acetaldehyde plus ammonia. We attempt to elucidate the mechanism by which the enzyme facilitates this reaction through interactions between active-site residues and the substrate. Our calculations suggest a preferred pathway involving a 1,2-shift in the associated radical and also suggest that interactions between the enzyme and the migrating group of the substrate that afford an almost fully protonated migrating group will lead to the most efficient catalysis. However, this criterion on its own is insufficient to fully understand the rearrangement. Additional synergistic interactions between the spectator hydroxyl group in the substrate and active-site residues on the enzyme are required to lower the barrier height to a value consistent with experimental observations.     

Wide-open flaps are key to urease activity

Substrate ingress and product egress from the active site of urease is tightly controlled by an active-site flap. Molecular dynamics simulations of urease have revealed a previously unobserved wide-open flap state that, unlike the well-characterized closed and open states, allows ready access to the metal cluster in the active site. This state is easily reached from the open state via low free energy barriers. Additionally, we have found that even when the flap is closed, a region of the binding pocket is solvent-exposed, leading to the hypothesis that it may act as a substrate/product reservoir. The newly identified wide-open state offers further opportunities for small-molecule drug discovery by defining a more extensive active-site pocket than has been previously described.     

DFT calculations for intermediate and active states of the diiron center with a tryptophan or tyrosine radical in Escherichia coli ribonucleotide reductase

Class Ia ribonucleotide reductase subunit R2 contains a diiron active site. In this paper, active-site models for the intermediate X-Trp48(?+) and X-Tyr122(?), the active Fe(III)Fe(III)-Tyr122(?), and the met Fe(III)Fe(III) states of Escherichia coli R2 are studied, using broken-symmetry density functional theory incorporated with the conductor-like screening solvation model. Different structural isomers and different protonation states have been explored. Calculated geometric, energetic, Mo?ssbauer, hyperfine, and redox properties are compared with available experimental data. Feasible detailed structures of these intermediate and active states are proposed. Asp84 and Trp48 are most likely the main contributing residues to the result that the transient Fe(IV)Fe(IV) state is not observed in wild-type class Ia E. coli R2. Asp84 is proposed to serve as a proton-transfer conduit between the diiron cluster and Tyr122 in both the tyrosine radical activation pathway and the first steps of the catalytic proton-coupled electron-transfer pathway. Proton-coupled and simple redox potential calculations show that the kinetic control of proton transfer to Tyr122(?) plays a critical role in preventing reduction from the active Fe(III)Fe(III)-Tyr122(?) state to the met state, which is potentially the reason why Tyr122(?) in the active state can be stable over a very long period.     

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.     

Quantitative analysis of the effect of salt concentration on enzymatic catalysis.

Like pH, salt concentration can have a dramatic effect on enzymatic catalysis. Here, a general equation is derived for the quantitative analysis of salt-rate profiles: k(cat)/K(M) = (k(cat)/K(M))(MAX)/[1+([Na+]/K[Na+])(n')], where (k(cat)/K(M))(MAX) is the physical limit of k(cat)/K(M), K(Na+) is the salt concentration at which k(cat)/K(M) = (k(cat)/K(M))(MAX)/2, and -n' is the slope of the linear region in a plot of log(k(cat)/K(M)) versus log [Na+]. The value of n' is of special utility, as it reflects the contribution of Coulombic interactions to the uniform binding of the bound states. This equation was used to analyze salt effects on catalysis by ribonuclease A (RNase A), which is a cationic enzyme that catalyzes the cleavage of an anionic substrate, RNA, with k(cat)/K(M) values that can exceed 10(9) M(-1) s(-1). Lys7, Arg10, and Lys66 comprise enzymic subsites that are remote from the active site. Replacing Lys7, Arg10, and Lys66 with alanine decreases the charge on the enzyme as well as the value of n'. Likewise, decreasing the number of phosphoryl groups in the substrate decreases the value of n'. Replacing Lys41, a key active-site residue, with arginine creates a catalyst that is limited by the chemical conversion of substrate to product. This change increases the value of n', as expected for a catalyst that is more sensitive to changes in the binding of the chemical transition state. Hence, the quantitative analysis of salt-rate profiles can provide valuable insight into the role of Coulombic interactions in enzymatic catalysis.     

Fluoride-mediated capture of a noncovalent bound state of a reversible covalent enzyme inhibitor: X-ray crystallographic analysis of an exceptionally potent α-ketoheterocycle inhibitor of fatty acid amide hydrolase

Two cocrystal X-ray structures of the exceptionally potent α-ketoheterocycle inhibitor 1 (K(i) = 290 pM) bound to a humanized variant of rat fatty acid amide hydrolase (FAAH) are disclosed, representing noncovalently and covalently bound states of the same inhibitor with the enzyme. Key to securing the structure of the noncovalently bound state of the inhibitor was the inclusion of fluoride ion in the crystallization conditions that is proposed to bind the oxyanion hole precluding inhibitor covalent adduct formation with stabilization of the tetrahedral hemiketal. This permitted the opportunity to detect important noncovalent interactions stabilizing the binding of the inhibitor within the FAAH active site independent of the covalent reaction. Remarkably, noncovalently bound 1 in the presence of fluoride appears to capture the active site in the same "in action" state with the three catalytic residues Ser241-Ser217-Lys142 occupying essentially identical positions observed in the covalently bound structure of 1, suggesting that this technique of introducing fluoride may have important applications in structural studies beyond inhibiting substrate or inhibitor oxyanion hole binding. Key insights to emerge from the studies include the observations that noncovalently bound 1 binds in its ketone (not gem diol) form, that the terminal phenyl group in the acyl side chain of the inhibitor serves as the key anchoring interaction overriding the intricate polar interactions in the cytosolic port, and that the role of the central activating heterocycle is dominated by its intrinsic electron-withdrawing properties. These two structures are also briefly compared with five X-ray structures of α-ketoheterocycle-based inhibitors bound to FAAH recently disclosed.     

Hidden Conformations in Aspergillus niger Monoamine Oxidase are Key for Catalytic Efficiency

Enzymes exist as an ensemble of conformational states, whose populations can be shifted by substrate binding, allosteric interactions, but also by introducing mutations to their sequence. Tuning the populations of the enzyme conformational states through mutation enables evolution towards novel activity. Herein, Markov state models are used to unveil hidden conformational states of monoamine oxidase from Aspergillus niger (MAO‐N). These hidden conformations, not previously observed by any other technique, play a crucial role in substrate binding and enzyme activity. This reveals how distal mutations regulate MAO‐N activity by stabilizing these hidden, catalytically important conformational states, but also by modulating the communication pathway between both MAO‐N subunits.     

A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.     

Tuning the regio- and stereoselectivity of C-H activation in n-octanes by cytochrome P450 BM-3 with fluorine substituents: evidence for interactions between a C-F bond and aromatic π systems

We employed the water-soluble cytochrome P450 BM-3 to study the activity and regiospecificity of oxidation of fluorinated n-octanes. Three mutations, A74G, F87V, and L188Q, were introduced into P450 BM-3 to allow the system to undergo n-octane oxidation. In addition, the alanine at residue 328 was replaced with a phenylalanine to introduce an aromatic residue into the hydrophobic pocket to examine whether or not van der Waals interactions between a C-F substituent in the substrate and the polarizable π system of the phenylalanine may be used to steer the positioning of the substrate within the active-site pocket of the enzyme and control the regioselectivity and stereoselectivity of hydroxylation. Interestingly, not only was the regioselectivity controlled when the fluorine substituent was judiciously positioned in the substrate, but the electron input into the iron-heme group became tightly coupled to the formation of product, essentially without abortive side reactions. Remarkable enhancement of the coupling efficiency between electron input and product formation was observed for a range of fluorinated octanes in the enzyme even without the A328F mutation, presumably because of interactions of the C-F substituent with the π system of the porphyrin macrocycle within the active-site pocket. Evidently, tightening the protein domain containing the heme pocket tunes the distribution of accessible enzyme conformations and the associated protein dynamics that activate the iron porphyrin for substrate hydroxylation to allow the reactions mediated by the high-valent Fe(IV)=O to become kinetically more commensurate with electron transfer from the flavin adenine dinucleotide (FAD)/flavin mononucleotide (FMN) reductase. These observations lend compelling evidence to support significant van der Waals interactions between the CF(2) group and aromatic π systems within the heme pocket when the fluorinated octane substrate is bound.     

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.     

Silver metallation of hen egg white lysozyme: X-ray crystal structure and NMR studies

The X-ray crystal structure, NMR binding studies, and enzyme activity of silver(I) metallated hen egg white lysozyme are presented. Primary bonding of silver is observed through His15 with secondary bonding interactions coming from nearby Arg14 and Asp87. A covalently bound nitrate completes a four coordinate binding pocket.     

Hybrid QM/MM and DFT investigations of the catalytic mechanism and inhibition of the dinuclear zinc metallo-beta-lactamase CcrA from Bacteroides fragilis

Based on hybrid QM/MM molecular dynamics simulation and density functional theoretical (DFT) calculations, we investigate the mechanistic and energetic features of the catalytic action of dizinc metallo-beta-lactamase CcrA from Bacteroides fragilis. The 200 ps QM/MM simulation of the CcrA enzyme in complex with nitrocefin shows that the substrate beta-lactam moiety is directed toward the active site dizinc center through the interactions of aminocarbonyl and carboxylate groups with the two active site zinc ions and the two conserved residues, Lys167 and Asn176. From the determination of the potential energy profile of a relevant enzymatic reaction model, it is found that the nucleophilic displacement reaction step proceeds with a low-barrier height, leading to the formation of an energetically favored reaction intermediate. The results also show that the high catalytic activity of the CcrA enzyme stems from a simultaneous operation of three catalytic components: activation of the bridging hydroxide nucleophile by zinc-coordinated Asp86; polarization of the substrate aminocarbonyl group by the first zinc ion; stabilization of the negative charge developed on the departing amide nitrogen by the second zinc ion. Consistent with the previous experimental finding that the proton-transfer reaction step is rate-limiting, the activation energy of the second step is found to be 1.6 kcal/mol higher than that of the first step. Finally, through an examination of the structural and energetic features of binding of a thiazolidinecarboxylic acid inhibitor to the active site dizinc center, a two-step inhibition mechanism involving a protonation-induced ligand exchange reaction is proposed for the inhibitory action of a tight-binding inhibitor possessing a thiol group.