Kinetic stability as a mechanism for protease longevity

The folding of the extracellular serine protease, α-lytic protease (αLP; EC 3.4.21.12) reveals a novel mechanism for stability that appears to lead to a longer functional lifetime for the protease. For αLP, stability is based not on thermodynamics, but on kinetics. Whereas this has required the coevolution of a pro region to facilitate folding, the result has been the optimization of native-state properties independent of their consequences on thermodynamic stability. Structural and mutational data lead to a model for catalysis of folding in which the pro region binds to a conserved β-hairpin in the αLP C-terminal domain, stabilizing the folding transition state and the native state. The pro region is then proteolytically degraded, leaving the active αLP trapped in a metastable conformation. This metastability appears to be a consequence of pressure to evolve properties of the native state, including a large, highly cooperative barrier to unfolding, and extreme rigidity, that reduce susceptibility to proteolytic degradation. In a test of survival under highly proteolytic conditions, homologous mammalian proteases that have not evolved kinetic stability are much more rapidly degraded than αLP. Kinetic stability as a means to longevity is likely to be a mechanism conserved among the majority of extracellular bacterial pro-proteases and may emerge as a general strategy for intracellular eukaryotic proteases subject to harsh conditions as well.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Mechanical and chemical unfolding of a single protein: A comparison

Is the mechanical unraveling of protein domains by atomic force microscopy (AFM) just a technological feat or a true measurement of their unfolding? By engineering a protein made of tandem repeats of identical Ig modules, we were able to get explicit AFM data on the unfolding rate of a single protein domain that can be accurately extrapolated to zero force. We compare this with chemical unfolding rates for untethered modules extrapolated to 0 M denaturant. The unfolding rates obtained by the two methods are the same. Furthermore, the transition state for unfolding appears at the same position on the folding pathway when assessed by either method. These results indicate that mechanical unfolding of a single protein by AFM does indeed reflect the same event that is observed in traditional unfolding experiments. The way is now open for the extensive use of AFM to measure folding reactions at the single-molecule level. Single-molecule AFM recordings have the added advantage that they define the reaction coordinate and expose rare unfolding events that cannot be observed in the absence of chemical denaturants.

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.