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.

Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage

Ribozymes of hepatitis delta virus have been proposed to use an active-site cytosine as an acid-base catalyst in the self-cleavage reaction. In this study, we have examined the role of cytosine in more detail with the antigenomic ribozyme. Evidence that proton transfer in the rate-determining step involved cytosine 76 (C76) was obtained from examining cleavage activity of the wild-type and imidazole buffer-rescued C76-deleted (C76Δ) ribozymes in D2O and H2O. In both reactions, a similar kinetic isotope effect and shift in the apparent pKa indicate that the buffer is functionally substituting for the side chain in proton transfer. Proton inventory of the wild-type reaction supported a mechanism of a single proton transfer at the transition state. This proton transfer step was further characterized by exogenous base rescue of a C76Δ mutant with cytosine and imidazole analogues. For the imidazole analogues that rescued activity, the apparent pKa of the rescue reaction, measured under kcat/KM conditions, correlated with the pKa of the base. From these data a Brønsted coefficient (β) of 0.51 was determined for the base-rescued reaction of C76Δ. This value is consistent with that expected for proton transfer in the transition state. Together, these data provide strong support for a mechanism where an RNA side chain participates directly in general acid or general base catalysis of the wild-type ribozyme to facilitate RNA cleavage.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

Femtochemistry

The topic of femtochemistry is surveyed from both theoretical and experimental points of view. A time-dependent wave packet description of the photodissociation of the O—C—S molecule reveals vibrational motion in the transition-state region and suggests targets for direct experimental observation. Theoretical approaches for treating femtosecond chemical phenomena in condensed phases are featured along with prospects for laser-controlled chemical reactions by using tailored ultrashort chirped pulses. An experimental study of the photoisomerization of retinal in the protein bacteriorhodopsin is discussed with an aim to gain insight into the potential energy surfaces on which this remarkably efficient and selective reactions proceeds. Finally, a prospective view of new frontiers in femtochemistry is given.

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.

Monitoring the GAP catalyzed H-Ras GTPase reaction at atomic resolution in real time

The molecular reaction mechanism of the GTPase-activating protein (GAP)-catalyzed GTP hydrolysis by Ras was investigated by time resolved Fourier transform infrared (FTIR) difference spectroscopy using caged GTP (P3-1-(2-nitro)phenylethyl guanosine 5′-O-triphosphate) as photolabile trigger. This approach provides the complete GTPase reaction pathway with time resolution of milliseconds at the atomic level. Up to now, one structural model of the GAP⋅Ras⋅GDP⋅AlFx transition state analog is known, which represents a “snap shot” along the reaction-pathway. As now revealed, binding of GAP to Ras⋅GTP shifts negative charge from the γ to β phosphate. Such a shift was already identified by FTIR in GTP because of Ras binding and is now shown to be enhanced by GAP binding. Because the charge distribution of the GAP⋅Ras⋅GTP complex thus resembles a more dissociative-like transition state and is more like that in GDP, the activation free energy is reduced. An intermediate is observed on the reaction pathway that appears when the bond between β and γ phosphate is cleaved. In the intermediate, the released Pi is strongly bound to the protein and surprisingly shows bands typical of those seen for phosphorylated enzyme intermediates. All these results provide a mechanistic picture that is different from the intrinsic GTPase reaction of Ras. FTIR analysis reveals the release of Pi from the protein complex as the rate-limiting step for the GAP-catalyzed reaction. The approach presented allows the study not only of single proteins but of protein–protein interactions without intrinsic chromophores, in the non-crystalline state, in real time at the atomic level.

SURVEY AND SUMMARY: Recent advances in the elucidation of the mechanisms of action of ribozymes

The cleavage of RNA can be accelerated by a number of factors. These factors include an acidic group (Lewis acid) or a basic group that aids in the deprotonation of the attacking nucleophile, in effect enhancing the nucleophilicity of the nucleophile; an acidic group that can neutralize and stabilize the leaving group; and any environment that can stabilize the pentavalent species that is either a transition state or a short-lived intermediate. The catalytic properties of ribozymes are due to factors that are derived from the complicated and specific structure of the ribozyme–substrate complex. It was postulated initially that nature had adopted a rather narrowly defined mechanism for the cleavage of RNA. However, recent findings have clearly demonstrated the diversity of the mechanisms of ribozyme-catalyzed reactions. Such mechanisms include the metal-independent cleavage that occurs in reactions catalyzed by hairpin ribozymes and the general double-metal-ion mechanism of catalysis in reactions catalyzed by the Tetrahymena group I ribozyme. Furthermore, the architecture of the complex between the substrate and the hepatitis delta virus ribozyme allows perturbation of the pKa of ring nitrogens of cytosine and adenine. The resultant perturbed ring nitrogens appear to be directly involved in acid/base catalysis. Moreover, while high concentrations of monovalent metal ions or polyamines can facilitate cleavage by hammerhead ribozymes, divalent metal ions are the most effective acid/base catalysts under physiological conditions.

Distortion in the spacer region of Pm during activation of middle transcription of phage Mu.

Transcription from the middle promoter, Pm, of phage Mu is initiated by Escherichia coli RNA polymerase holoenzyme (E sigma 70; RNAP) and the phage-encoded activator, Mor. Point mutations in the spacer region between the -10 hexamer and the Mor binding site result in changes of promoter activity in vivo. These mutations are located at the junction between a rigid T-tract and adjacent, potentially deformable G + C-rich DNA segment, suggesting that deformation of the spacer region may play a role in the transcriptional activation of Pm. This prediction was tested by using dimethyl sulfate and potassium permanganate footprinting analyses. Helical distortion involving strand separation was detected at positions -32 to -34, close to the predicted interface between Mor and RNAP. Promoter mutants in which this distortion was not detected exhibited a lack of melting in the -12 to -1 region and reduced promoter activity in vivo. We propose that complexes containing the distortion represent stressed intermediates rather than stable open complexes and thus can be envisaged as a transition state in the kinetic pathway of Pm activation in which stored torsional energy could be used to facilitate melting around the transcription start point.

Ribonuclease P (RNase P) RNA is converted to a Cd(2+)-ribozyme by a single Rp-phosphorothioate modification in the precursor tRNA at the RNase P cleavage site.

To study the cleavage mechanism of bacterial Nase P RNA, we have synthesized precursor tRNA substrates carrying a single Rp- or Sp-phosphorothioate modification at the RNase P cleavage site. Both the Sp- and the Rp-diastereomer reduced the rate of processing by Escherichia coli RNase P RNA at least 1000-fold under conditions where the chemical step is rate-limiting. The Rp-modification had no effect and the Sp-modification had a moderate effect on precursor tRNA ground state binding to RNase P RNA. Processing of the Rp-diastereomeric substrate was largely restored in the presence of the "thiophilic" Cd2+ as the only divalent metal ion, demonstrating direct metal ion coordination to the (pro)-Rp substituent at the cleavage site and arguing against a specific role for Mg(2+)-ions at the pro-Sp oxygen. For the Rp-diastereomeric substrate, Hill plot analysis revealed a cooperative dependence upon [Cd2+] of nH = 1.8, consistent with a two-metal ion mechanism. In the presence of the Sp-modification, neither Mn2+ nor Cd2+ was able to restore detectable cleavage at the canonical site. Instead, the ribozyme promotes cleavage at the neighboring unmodified phosphodiester with low efficiency. Dramatic inhibition of the chemical step by both the Rp- and Sp-phosphorothioate modification is unprecedented among known ribozymes and points to unique features of transition state geometry in the RNase P RNA-catalyzed reaction.