Role of a critical water in scytalone dehydratase-catalyzed reaction

Scytalone dehydratase (EC 4.2.1.94) catalyzes the dehydration of two important intermediates in the biosynthesis of melanin, and it functions without metal ions or any cofactors. Using molecular orbital theory, we have examined the role of a critical water molecule in the mechanism of scytalone dehydratase. The water, together with an internal hydrogen bonding, contributes significantly to the stabilization of the transition state (or the enolate intermediate). The role of two active site tyrosines (Tyr-50 and Tyr-30) is (i) to hold the critical water in place so that it may stabilize the transition state without much structural rearrangement during the catalytic reaction, and (ii) to polarize the water, making it a better general acid. The stereochemistry of the scytalone dehydratase-catalyzed dehydration is also discussed.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Use of intrinsic binding energy for catalysis by an RNA enzyme

The contribution of several individual ribozyme⋅substrate base pairs to binding and catalysis has been investigated using hammerhead ribozyme substrates that were truncated at their 3′ or 5′ ends. The base pairs at positions 1.1–2.1 and 15.2–16.2, which flank the conserved core, each contribute 104-fold in the chemical step, without affecting substrate binding. In contrast, base pairs distal to the core contribute to substrate binding but have no effect on the chemical step. These results suggest a “fraying model” in which each ribozyme⋅substrate helix can exist in either an unpaired (“open”) state or a helical (“closed”) state, with the closed state required for catalysis. The base pairs directly adjacent to the conserved core contribute to catalysis by allowing the closed state to form. Once the number of base pairs is sufficient to ensure that the closed helical state predominates, additional residues provide stabilization of the helix, and therefore increase binding, but have no further effect on the chemical step. Remarkably, the >5 kcal/mol free energy contribution to catalysis from each of the internal base pairs is considerably greater than the free energy expected for formation of a base pair. It is suggested that this unusually large energetic contribution arises because free energy that is typically lost in constraining residues within a base pair is expressed in the transition state, where it is used for positioning. This extends the concept of “intrinsic binding energy” from protein to RNA enzymes, suggesting that intrinsic binding energy is a fundamental feature of biological catalysis.

On the reaction mechanism of l-lactate oxidase: Quantitative structure-activity analysis of the reaction with para-substituted l-mandelates

The rate constants for reduction of the flavoenzyme, l-lactate oxidase, and a mutant (in which alanine 95 is replaced by glycine), by a series of para-substituted mandelates, in both the 2-1H- and 2-2H- forms, have been measured by rapid reaction spectrophotometry. In all cases, significant isotope effects (1H/2H = 3–7) on the rate constants of flavin reduction were found, indicating that flavin reduction is a direct measure of α-C-H bond breakage. The rate constants show only a small influence of the electronic characteristics of the substituents, but show a good correlation when combined with some substituent volume parameters. A surprisingly good correlation is found with the molecular mass of the substrate. The results are compatible with any mechanism in which there is little development of charge in the transition state. This could be a transfer of hydride to the flavin N(5) position or a synchronous mechanism in which the α-C-H is formally abstracted as a H+ while the resulting charge is simultaneously neutralized by another event.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Conformation of coenzyme pyrroloquinoline quinone and role of Ca2+ in the catalytic mechanism of quinoprotein methanol dehydrogenase

The ab initio structures of 2,7,9-tricarboxypyrroloquinoline quinone (PQQ), semiquinone (PQQH), and dihydroquinone (PQQH2) have been determined and compared with ab initio structures of the (PQQ)Ca2+, (PQQH)Ca2+, and (PQQH2)Ca2+ complexes as well as the x-ray structure of (PQQ)Ca2+ bound at the active site of the methanol dehydrogenase (MDH) of methyltropic bacteria. Plausible mechanisms for the MDH oxidation of methanol involving the (PQQ)Ca2+ complex are explored via ab initio computations and discussed. Considering the reaction of methanol with PQQ in the absence of Ca2+, nucleophilic addition of methanol to the PQQ C-5 carbonyl followed by a retro-ene elimination is deemed unlikely due to large energy barrier. A much more favorable disposition of the methanol C-5 adduct to provide formaldehyde involves proton ionization of the intermediate followed by elimination of methoxide concerted with hydride transfer to the oxygen of the C-4 carbonyl. Much the same transition state is reached if one searches for the transition state beginning with Asp-303–CO2−general-base removal of the methanol proton of the (PQQ)Ca2+O(H)CH3 complex concerted with hydride transfer to the oxygen at C-4. For such a mechanism the role of the Ca2+ moiety would be to (i) contribute to the formation of the ES complex (ii) provide a modest decrease in the pKa of methanol substrate,; and (iii) polarize the oxygen at C-5.

Inhibition of regulator of G protein signaling function by two mutant RGS4 proteins

Regulators of G protein signaling (RGS) proteins limit the lifetime of activated (GTP-bound) heterotrimeric G protein α subunits by acting as GTPase-activating proteins (GAPs). Mutation of two residues in RGS4, which, based on the crystal structure of RGS4 complexed with Giα1-GDP-AlF4−, directly contact Giα1 (N88 and L159), essentially abolished RGS4 binding and GAP activity. Mutation of another contact residue (S164) partially inhibited both binding and GAP activity. Two other mutations, one of a contact residue (R167M/A) and the other an adjacent residue (F168A), also significantly reduced RGS4 binding to Giα1-GDP-AlF4−, but in addition redirected RGS4 binding toward the GTPγS-bound form. These two mutant proteins had severely impaired GAP activity, but in contrast to the others behaved as RGS antagonists in GAP and in vivo signaling assays. Overall, these results are consistent with the hypothesis that the predominant role of RGS proteins is to stabilize the transition state for GTP hydrolysis. In addition, mutant RGS proteins can be created with an altered binding preference for the Giα-GTP conformation, suggesting that efficient RGS antagonists can be developed.

Structural basis for efficient phosphorylation of 3′-azidothymidine monophosphate by Escherichia coli thymidylate kinase

The crystal structures of Escherichia coli thymidylate kinase (TmpK) in complex with P1-(5′-adenosyl)-P5-(5′-thymidyl)pentaphosphate and P1-(5′-adenosyl)P5-[5′-(3′-azido-3′-deoxythymidine)] pentaphosphate have been solved to 2.0-Å and 2.2-Å resolution, respectively. The overall structure of the bacterial TmpK is very similar to that of yeast TmpK. In contrast to the human and yeast TmpKs, which phosphorylate 3′-azido-3′-deoxythymidine 5′-monophosphate (AZT-MP) at a 200-fold reduced turnover number (kcat) in comparison to the physiological substrate dTMP, reduction of kcat is only 2-fold for the bacterial enzyme. The different kinetic properties toward AZT-MP between the eukaryotic TmpKs and E. coli TmpK can be rationalized by the different ways in which these enzymes stabilize the presumed transition state and the different manner in which a carboxylic acid side chain in the P loop interacts with the deoxyribose of the monophosphate. Yeast TmpK interacts with the 3′-hydroxyl of dTMP through Asp-14 of the P loop in a bidentate manner: binding of AZT-MP results in a shift of the P loop to accommodate the larger substituent. In E. coli TmpK, the corresponding residue is Glu-12, and it interacts in a side-on fashion with the 3′-hydroxyl of dTMP. This different mode of interaction between the P loop carboxylic acid with the 3′ substituent of the monophosphate deoxyribose allows the accommodation of an azido group in the case of the E. coli enzyme without significant P loop movement. In addition, although the yeast enzyme uses Arg-15 (a glycine in E. coli) to stabilize the transition state, E. coli seems to use Arg-153 from a region termed Lid instead. Thus, the binding of AZT-MP to the yeast TmpK results in the shift of a catalytic residue, which is not the case for the bacterial kinase.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

How does a β-hairpin fold/unfold? Competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a β-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

Interaction of two arginine residues in lactate oxidase with the enzyme flavin: Conversion of FMN to 8-formyl-FMN

Two arginine residues, Arg-181 and Arg-268, are conserved throughout the known family of FMN-containing enzymes that catalyze the oxidation of α-hydroxyacids. In the lactate oxidase from Aerococcus viridans, these residues have been changed to lysine in two single mutations and in a double mutant form. In addition, Arg-181 has been replaced by methionine to determine the effect of removing the positive charge on the residue. The effects of these replacements on the kinetic and thermodynamic properties are reported. With all mutant forms, there are only small effects on the reactivity of the reduced flavin with oxygen. On the other hand, the efficiency of reduction of the oxidized flavin by l-lactate is greatly reduced, particularly with the R268K mutant forms. The results demonstrate the importance of the two arginine residues in the binding of substrate and its interaction with the flavin, and are consistent with a previous hypothesis that they also play a role of charge neutralization in the transition state of substrate dehydrogenation. The replacement of Arg-268 by lysine also results in a slow conversion of the 8-CH3- substituent of FMN to yield 8-formyl-FMN, still tightly bound to the enzyme, and with significantly different physical and chemical properties from those of the FMN-enzyme.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.