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.

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.

Ras-catalyzed hydrolysis of GTP: a new perspective from model studies.

Despite the biological and medical importance of signal transduction via Ras proteins and despite considerable kinetic and structural studies of wild-type and mutant Ras proteins, the mechanism of Ras-catalyzed GTP hydrolysis remains controversial. We take a different approach to this problem: the uncatalyzed hydrolysis of GTP is analyzed, and the understanding derived is applied to the Ras-catalyzed reaction. Evaluation of previous mechanistic proposals from this chemical perspective suggests that proton abstraction from the attacking water by a general base and stabilization of charge development on the gamma-phosphoryl oxygen atoms would not be catalytic. Rather, this analysis focuses attention on the GDP leaving group, including the beta-gamma bridge oxygen of GTP, the atom that undergoes the largest change in charge in going from the ground state to the transition state. This leads to a new catalytic proposal in which a hydrogen bond from the backbone amide of Gly-13 to this bridge oxygen is strengthened in the transition state relative to the ground state, within an active site that provides a template complementary to the transition state. Strengthened transition state interactions of the active site lysine, Lys-16, with the beta-nonbridging phosphoryl oxygens and a network of interactions that positions the nucleophilic water molecule and gamma-phosphoryl group with respect to one another may also contribute to catalysis. It is speculated that a significant fraction of the GAP-activated GTPase activity of Ras arises from an additional interaction of the beta-gamma bridge oxygen with an Arg side chain that is provided in trans by GAP. The conclusions for Ras and related G proteins are expected to apply more widely to other enzymes that catalyze phosphoryl (-PO(3)2-) transfer, including kinases and phosphatases.

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Insights into antibody catalysis: structure of an oxygenation catalyst at 1.9-angstrom resolution.

The x-ray crystal structures of the sulfide oxidase antibody 28B4 and of antibody 28B4 complexed with hapten have been solved at 2.2-angstrom and 1.9-angstrom resolution, respectively. To our knowledge, these structures are the highest resolution catalytic antibody structures to date and provide insight into the molecular mechanism of this antibody-catalyzed monooxygenation reaction. Specifically, the data suggest that entropic restriction plays a fundamental role in catalysis through the precise alignment of the thioether substrate and oxidant. The antibody active site also stabilizes developing charge on both sulfur and periodate in the transition state via cation-pi and electrostatic interactions, respectively. In addition to demonstrating that the active site of antibody 28B4 does indeed reflect the mechanistic information programmed in the aminophosphonic acid hapten, these high-resolution structures provide a basis for enhancing turnover rates through mutagenesis and improved hapten design.

Exploring the active site of chorismate mutase by combinatorial mutagenesis and selection: the importance of electrostatic catalysis.

Chorismate mutase (EC 5.4.99.5) catalyzes the intramolecular rearrangement of chorismate to prephenate. Arg-90 in the active site of the enzyme from Bacillus subtilis is in close proximity to the substrate's ether oxygen and may contribute to efficient catalysis by stabilizing the presumed dipolar transition state that would result upon scission of the C--O bond. To test this idea, we have developed a novel complementation system for chorismate mutase activity in Escherichia coli by reengineering parts of the aromatic amino acid biosynthetic pathway. The codon for Arg-90 was randomized, alone and in combination with that for Cys-88, and active clones were selected. The results show that a positively charged residue either at position 88 (Lys) or 90 (Arg or Lys) is essential. Our data provide strong support for the hypothesis that the positive charge is required for stabilization of the transition state of the enzymatic chorismate rearrangement. The new selection system, in conjunction with combinatorial mutagenesis, renders the mechanism of the natural enzyme(s) accessible to further exploration and opens avenues for the improvement of first generation catalytic antibodies with chorismate mutase activity.