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.

Floral Scent Production in Clarkia breweri1 : III. Enzymatic Synthesis and Emission of Benzenoid Esters

The fragrance of Clarkia breweri (Onagraceae), a California annual plant, includes three benzenoid esters: benzylacetate, benzylbenzoate, and methylsalicylate. Here we report that petal tissue was responsible for the benzylacetate and methylsalicylate emission, whereas the pistil was the main source of benzylbenzoate. The activities of two novel enzymes, acetyl-coenzyme A:benzylalcohol acetyltransferase (BEAT), which catalyzes the acetyl esterification of benzylalcohol, and S-adenosyl-l-methionine:salicylic acid carboxyl methyltransferase, which catalyzes the methyl esterification of salicylic acid, were also highest in petal tissue and absent in leaves. In addition, the activity of both enzymes in the various floral organs was developmentally and differentially regulated. S-Adenosyl-l-methionine:salicylic acid carboxyl methyltransferase activity in petals peaked in mature buds and declined during the next few days after anthesis, and it showed a strong, positive correlation with the emission of methylsalicylate. The levels of BEAT activity and benzylacetate emission in petals also increased in parallel as the buds matured and the flowers opened, but as emission began to decline on the 2nd d after anthesis, BEAT activity continued to increase and remained high until the end of the lifespan of the flower.

Molecular and Enzymatic Characterization of Three Phosphoinositide-Specific Phospholipase C Isoforms from Potato1

Many cellular responses to stimulation of cell-surface receptors by extracellular signals are transmitted across the plasma membrane by hydrolysis of phosphatidylinositol-4,5-bisphosphate (PIP2), which is cleaved into diacylglycerol and inositol-1,4,5-tris-phosphate by phosphoinositide-specific phospholipase C (PI-PLC). We present structural, biochemical, and RNA expression data for three distinct PI-PLC isoforms, StPLC1, StPLC2, and StPLC3, which were cloned from a guard cell-enriched tissue preparation of potato (Solanum tuberosum) leaves. All three enzymes contain the catalytic X and Y domains, as well as C2-like domains also present in all PI-PLCs. Analysis of the reaction products obtained from PIP2 hydrolysis unequivocally identified these enzymes as genuine PI-PLC isoforms. Recombinant StPLCs showed an optimal PIP2-hydrolyzing activity at 10 μm Ca2+ and were inhibited by Al3+ in equimolar amounts. In contrast to PI-PLC activity in plant plasma membranes, however, recombinant enzymes could not be activated by Mg2+. All three stplc genes are expressed in various tissues of potato, including leaves, flowers, tubers, and roots, and are affected by drought stress in a gene-specific manner.

Energy metabolism, enzymatic flux capacities, and metabolic flux rates in flying honeybees.

Honeybees rely primarily on the oxidation of hexose sugars to provide the energy required for flight. Measurement of VCO2 (equal to VO2, because VCO2/VO2 = 1.0 during carbohydrate oxidation) during flight allowed estimation of steady-state flux rates through pathways of flight muscle energy metabolism. Comparison of Vmax values for flight muscle hexokinase, phosphofructokinase, citrate synthase, and cytochrome c oxidase with rates of carbon and O2 flux during flight reveal that these enzymes operate closer to Vmax in the flight muscles of flying honeybees than in other muscles previously studied. Possible mechanistic and evolutionary implications of these findings are discussed.

Enzymatic phosphorylation of muscle glycogen synthase: a mechanism for maintenance of metabolic homeostasis.

We recently analyzed experimental studies of mammalian muscle glycogen synthesis using metabolic control analysis and concluded that glycogen synthase (GSase) does not control the glycogenic flux but rather adapts to the flux which is controlled bv the activity of the proximal glucose transport and hexokinase steps. This model did not provide a role for the well established relationship between GSase fractional activity, determined by covalent phosphorylation, and the rate of glycogen synthesis. Here we propose that the phosphorylation of GSase, which alters the sensitivity to allosteric activation by glucose 6-phosphate (G6P), is a mechanism for controlling the concentration of G6P instead of controlling the flux. When the muscle cell is exposed to conditions which favor glycogen synthesis such as high plasma insulin and glucose concentrations the fractional activity of GSase is increased in coordination with increases in the activity of glucose transport and hexokinase. This increase in GSase fractional activity helps to maintain G6P homeostasis by reducing the G6P concentration required to activate GSase allosterically to match the flux determined by the proximal reactions. This role for covalent phosphorylation also provides a novel solution to the Kacser and Acarenza paradigm which requires coordinated activity changes of the enzymes proximal and distal to a shared intermediate, to avoid unwanted flux changes.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

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.

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.

Swiveling-domain mechanism for enzymatic phosphotransfer between remote reaction sites.

The crystal structure of pyruvate phosphate dikinase, a histidyl multiphosphotransfer enzyme that synthesizes adenosine triphosphate, reveals a three-domain molecule in which the phosphohistidine domain is flanked by the nucleotide and the phosphoenolpyruvate/pyruvate domains, with the two substrate binding sites approximately 45 angstroms apart. The modes of substrate binding have been deduced by analogy to D-Ala-D-Ala ligase and to pyruvate kinase. Coupling between the two remote active sites is facilitated by two conformational states of the phosphohistidine domain. While the crystal structure represents the state of interaction with the nucleotide, the second state is achieved by swiveling around two flexible peptide linkers. This dramatic conformational transition brings the phosphocarrier residue in close proximity to phosphoenolpyruvate/pyruvate. The swiveling-domain paradigm provides an effective mechanism for communication in complex multidomain/multiactive site proteins.

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.

Mutagenesis of the potato ADPglucose pyrophosphorylase and characterization of an allosteric mutant defective in 3-phosphoglycerate activation.

ADPglucose pyrophosphorylase (glucose-1-phosphate adenylyltransferase; ADP:alpha-D-glucose-1-phosphate adenylyltransferase, EC 2.7.7.27) catalyzes a key regulatory step in alpha-glucan synthesis in bacteria and higher plants. We have previously shown that the expression of the cDNA sequences of the potato tuber large (LS) and small (SS) subunits yielded a functional heterotetrameric enzyme capable of complementing a mutation in the single AGP (glgC) structural gene of Escherichia coli. This heterologous complementation provides a powerful genetic approach to obtain biochemical information on the specific roles of LS and SS in enzyme function. By mutagenizing the LS cDNA with hydroxylamine and then coexpressing with wild-type SS in an E. coli glgC- strain, >350 mutant colonies were identified that were impaired in glycogen production. One mutant exhibited enzymatic and antigen levels comparable to the wild-type recombinant enzyme but required 45-fold greater levels of the activator 3-phosphoglycerate for maximum activity. Sequence analysis identified a single nucleotide change that resulted in the change of Pro-52 to Leu. This heterologous genetic system provides an efficient means to identify residues important for catalysis and allosteric functioning and should lead to novel approaches to increase plant productivity.

Selection of transduced CD34+ progenitors and enzymatic correction of cells from Gaucher patients, with bicistronic vectors.

The gene transfer efficiency of human hematopoietic stem cells is still inadequate for efficient gene therapy of most disorders. To overcome this problem, a selectable retroviral vector system for gene therapy has been developed for gene therapy of Gaucher disease. We constructed a bicistronic retroviral vector containing the human glucocerebrosidase (GC) cDNA and the human small cell surface antigen CD24 (243 bp). Expression of both cDNAs was controlled by the long terminal repeat enhancer/promoter of the Molony murine leukemia virus. The CD24 selectable marker was placed downstream of the GC cDNA and its translation was enhanced by inclusion of the long 5' untranslated region of encephalomyocarditis virus internal ribosomal entry site. Virus-producing GP+envAM12 cells were created by multiple supernatant transductions to create vector producer cells. The vector LGEC has a high titer and can drive expression of GC and the cell surface antigen CD24 simultaneously in transduced NIH 3T3 cells and Gaucher skin fibroblasts. These transduced cells have been successfully separated from untransduced cells by fluorescence-activated cell sorting, based on cell surface expression of CD24. Transduced and sorted NIH 3T3 cells showed higher GC enzyme activity than the unsorted population, demonstrating coordinated expression of both genes. Fibroblasts from Gaucher patients were transduced and sorted for CD24 expression, and GC enzyme activity was measured. The transduced sorted Gaucher fibroblasts had a marked increase in enzyme activity (149%) compared with virgin Gaucher fibroblasts (17% of normal GC enzyme activity). Efficient transduction of CD34+ hematopoietic progenitors (20-40%) was accomplished and fluorescence-activated cell sorted CD24(+)-expressing progenitors generated colonies, all of which (100%) were vector positive. The sorted, CD24-expressing progenitors generated erythroid burst-forming units, colony-forming units (CFU)-granulocyte, CFU-macrophage, CFU-granulocyte/macrophage, and CFU-mix hematopoietic colonies, demonstrating their ability to differentiate into these myeloid lineages in vitro. The transduced, sorted progenitors raised the GC enzyme levels in their progeny cells manyfold compared with untransduced CD34+ progenitors. Collectively, this demonstrates the development of high titer, selectable bicistronic vectors that allow isolation of transduced hematopoietic progenitors and cells that have been metabolically corrected.