Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Binding and disruption of phospholipid bilayers by supramolecular RNA complexes

In an RNA world, RNAs would have regulated traffic through normally impermeable bilayer membranes. Using selection-amplification we previously found RNAs that bind stably and increase the ionic conductance of phospholipid membranes at high Mg2+ and Ca2+ concentrations. Now selection in reduced divalents yields RNAs that bind phosphatidylcholine liposomes under conditions closer to physiological. Such affinity for phospholipid membranes requires interactions between RNAs. In fact, we detected no functional monomeric membrane-binding RNAs. A membrane-active end-to-end heterotrimer consisting of 2 RNA 9 and 1 RNA 10 is defined by nucleotide protection, oligonucleotide competition, and mutant analysis. Oligomers of the heterotrimer bind stably, cause release of liposome-encapsulated solutes, and disrupt model black membranes. Individual RNA molecules do not show any of these activities. This novel mechanism of RNA binding to lipid membranes may not only regulate membrane permeability, but suggests that arrays of catalytic or structural RNAs on membranes are plausible. Finally, a selection met only by RNA complexes evokes new possibilities for selection-amplification itself.

Sir3-dependent assembly of supramolecular chromatin structures in vitro

Baculovirus-expressed recombinant Sir3p (rSir3p) has been purified to near homogeneity, and its binding to naked DNA, mononucleosomes, and nucleosomal arrays has been characterized in vitro. At stoichiometric levels rSir3p interacts with intact nucleosomal arrays, mononucleosomes, and naked DNA, as evidenced by formation of supershifted species on native agarose gels. Proteolytic removal of the core histone tail domains inhibits but does not completely abolish rSir3p binding to nucleosomal arrays. The linker DNA in the supershifted complexes remains freely accessible to restriction endonuclease digestion, suggesting that both the tail domains and nucleosomal DNA contribute to rSir3p–chromatin interactions. Together these data indicate that rSir3p cross-links individual nucleosomal arrays into supramolecular assemblies whose physical properties transcend those of typical 10-nm and 30-nm fibers. Based on these data we hypothesize that Sir3p functions, at least in part, by mediating reorganization of the canonical chromatin fiber into functionally specialized higher order chromosomal domains.

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.

The hard-soft acid-base principle in enzymatic catalysis: dual reactivity of phosphoenolpyruvate.

In this paper, the chemical reactivity of C3 of phosphoenolpyruvate (PEP) has been analyzed in terms of density functional theory quantified through quantum chemistry calculations. PEP is involved in a number of important enzymatic reactions, in which its C3 atom behaves like a base. In three different enzymatic reactions analyzed here, C3 sometimes behaves like a soft base and sometimes behaves like a hard base in terms of the hard-soft acid-base principle. This dual nature of C3 of PEP was found to be related to the conformational change of the molecule. This leads to a testable hypothesis: that PEP adopts particular conformations in the enzyme-substrate complexes of different PEP-using enzymes, and that the enzymes control the reactivity through controlling the dihedral angle between the carboxylate and the C==C double bond of PEP.

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.

A general assay for antibody catalysis using acridone as a fluorescent tag.

A simple and highly sensitive catalysis assay is demonstrated based on analyzing reactions with acridonetagged compounds by thin-layer chromatography. As little as 1 pmol of product is readily visualized by its blue fluorescence under UV illumination and identified by its retention factor (Rf). Each assay requires only 10 microliters of solution. The method is reliable, inexpensive, versatile, and immediately applicable in repetitive format for screening catalytic antibody libraries. Three examples are presented: (i) the epoxidation of acridone labeled (S)-citronellol. The pair of stereoisomeric epoxides formed is resolved on the plate, which provides a direct selection method for enantioselective epoxidation catalysts. (ii) Oxidation of acridone-labeled 1-hexanol to 1-hexanal. The activity of horse liver alcohol dehydrogenase is detected. (iii) Indirect product labeling of released aldehyde groups by hydrazone formation with an acridone-labeled hydrazide. Activity of catalytic antibodies for hydrolysis of enol ethers is detected.

Visualization of intermediate and transition-state structures in protein-tyrosine phosphatase catalysis.

Engineering site-specific amino acid substitutions into the protein-tyrosine phosphatase (PTPase) PTP1 and the dual-specific vaccinia H1-related phosphatase (VHR), has kinetically isolated the two chemical steps of the reaction and provided a rare opportunity for examining transition states and directly observing the phosphoenzyme intermediate. Changing serine to alanine in the active-site sequence motif HCXXGXXRS shifted the rate-limiting step from intermediate formation to intermediate hydrolysis. Using phosphorus 31P NMR, the covalent thiol-phosphate intermediate was directly observed during catalytic turnover. The importance of the conserved aspartic acid (D92 in VHR and D181 in PTP1) in both chemical steps was established. Kinetic analysis of D92N and D181N mutants indicated that aspartic acid acts as a general acid by protonating the leaving-group phenolic oxygen. Structure-reactivity experiments with native and aspartate mutant enzymes established that proton transfer is concomitant with P-O cleavage, such that no charge develops on the phenolic oxygen. Steady- and presteady-state kinetics, as well as NMR analysis of the double mutant D92N/S131A (VHR), suggested that the conserved aspartic acid functions as a general base during intermediate hydrolysis. As a general base, aspartate would activate a water molecule to facilitate nucleophilic attack. The amino acids involved in transition-state stabilization for cysteinylphosphate hydrolysis were confirmed by the x-ray structure of the Yersinia PTPase complexed with vanadate, a transition-state mimic that binds covalently to the active-site cysteine. Consistent with the NMR, x-ray, biochemical, and kinetic data, a unifying mechanism for catalysis is proposed.

Rotation of subunits during catalysis by Escherichia coli F1-ATPase.

During oxidative and photo-phosphorylation, F0F1-ATP synthases couple the movement of protons down an electrochemical gradient to the synthesis of ATP. One proposed mechanistic feature that has remained speculative is that this coupling process requires the rotation of subunits within F0F1. Guided by a recent, high-resolution structure for bovine F1 [Abrahams, J. P., Leslie, A. G., Lutter, R. & Walker, J. E. (1994) Nature (London) 370, 621-628], we have developed a critical test for rotation of the central gamma subunit relative to the three catalytic beta subunits in soluble F1 from Escherichia coli. In the bovine F1 structure, a specific point of contact between the gamma subunit and one of the three catalytic beta subunits includes positioning of the homolog of E. coli gamma-subunit C87 (gamma C87) close to the beta-subunit 380DELSEED386 sequence. A beta D380C mutation allowed us to induce formation of a specific disulfide bond between beta and gamma C87 in soluble E. coli F1. Formation of the crosslink inactivated beta D380C-F1, and reduction restored full activity. Using a dissociation/reassembly approach with crosslinked beta D380C-F1, we incorporated radiolabeled beta subunits into the two noncrosslinked beta-subunit positions of F1. After reduction of the initial nonradioactive beta-gamma crosslink, only exposure to conditions for catalytic turnover results in similar reactivities of unlabeled and radiolabeled beta subunits with gamma C87 upon reoxidation. The results demonstrate that gamma subunit rotates relative to the beta subunits during catalysis.

Carbohydrate gluing, an architectural mechanism in the supramolecular structure of an annelid giant hemoglobin.

We report a carbohydrate-dependent supramolecular architecture in the extracellular giant hemoglobin (Hb) from the marine worm Perinereis aibuhitensis; we call this architectural mechanism carbohydrate gluing. This study is an extension of our accidental discovery of deterioration in the form of the Hb caused by a high concentration of glucose. The giant Hbs of annelids are natural supramolecules consisting of about 200 polypeptide chains that associate to form a double-layered hexagonal structure. This Hb has 0.5% (wt) carbohydrates, including mannose, xylose, fucose, galactose, glucose, N-acetylglucosamine (GlcNAc), and N-acetylgalactosamine (GalNAc). Using carbohydrate-staining assays, in conjunction with two-dimensional polyacrylamide gel electrophoresis, we found that two types of linker chains (L1 and L2; the nomenclature of the Hb subunits followed that for another marine worm, Tylorrhynchus heterochaetus) contained carbohydrates with both GlcNAc and GalNAc. Furthermore, two types of globins (a and A) have only GlcNAc-containing carbohydrates, whereas the other types of globins (b and B) had no carbohydrates. Monosaccharides including mannose, fucose, glucose, galactose, GlcNAc, and GalNAc reversibly dissociated the intact form of the Hb, but the removal of carbohydrate with N-glycanase resulted in irreversible dissociation. These results show that carbohydrate acts noncovalently to glue together the components to yield the complete quaternary supramolecular structure of the giant Hb. We suggest that this carbohydrate gluing may be mediated through lectin-like carbohydrate-binding by the associated structural chains ("linkers").

Use of binding energy by an RNA enzyme for catalysis by positioning and substrate destabilization.

A fundamental catalytic principle for protein enzymes in the use of binding interactions away from the site of chemical transformation for catalysis. We have compared the binding and reactivity of a series of oligonucleotide substrates and products of the Tetrahymena ribozyme, which catalyzes a site-specific phosphodiester cleavage reaction: CCCUCUpA+G<-->CCCUCU-OH+GpA. The results suggest that this RNA enzyme, like protein enzymes, can utilize binding interactions to achieve substantial catalysis via entropic fixation and substrate destabilization. The stronger binding of the all-ribose oligonucleotide product compared to an analog with a terminal 3' deoxyribose residue gives an effective concentration of 2200 M for the 3' hydroxyl group, a value approaching those obtained with protein enzymes and suggesting the presence of a structurally well defined active site capable of precise positioning. The stabilization from tertiary binding interactions is 40-fold less for the oligonucleotide substrate than the oligonucleotide product, despite the presence of the reactive phosphoryl group in the substrate. This destabilization is accounted for by a model in which tertiary interactions away from the site of bond cleavage position the electron-deficient 3' bridging phosphoryl oxygen of the oligonucleotide substrate next to an electropositive Mg ion. As the phosphodiester bond breaks and this 3' oxygen atom develops a negative charge in the transition state, the weak interaction of the substrate with Mg2+ becomes strong. These strategies of "substrate destabilization" and "transition state stabilization" provide estimated rate enhancements of approximately 280- and approximately 60-fold, respectively. Analogous substrate destabilization by a metal ion or hydrogen bond donor may be used more generally by RNA and protein enzymes catalyzing reactions of phosphate esters.