The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Functional and physiological consequences of genetic variation at phosphoglucose isomerase: Heat shock protein expression is related to enzyme genotype in a montane beetle

Allele frequency variation at the phosphoglucose isomerase (PGI) locus in Californian populations of the beetle Chrysomela aeneicollis suggests that PGI may be undergoing natural selection. We quantified (i) apparent Michaelis-Menten constant (Km) of fructose 6-phosphate at different temperatures and (ii) thermal stability for three common PGI genotypes (1–1, 1–4, and 4–4). We also measured air temperature (Ta) and beetle body temperature (Tb) in three montane drainages in the Sierra Nevada, California. Finally, we measured 70-kDa heat shock protein (Hsp70) expression in field-collected and laboratory-acclimated beetles. We found that PGI allele 1 predominated in the northernmost drainage, Rock Creek (RC), which was also significantly cooler than the southernmost drainage, Big Pine Creek (BPC), where PGI allele 4 predominated. Allele frequencies and air temperatures were intermediate in the middle drainage, Bishop Creek (BC). Differences among genotypes in Km (1–1 > 1–4 > 4–4) and thermal stability (4–4 > 1–4 > 1–1) followed a pattern consistent with temperature adaptation. In nature, Tb was closely related to Ta. Hsp70 expression in adult beetles decreased with elevation and differed among drainages (BPC > BC > RC). After laboratory acclimation (8 days, 20°C day, 4°C night) and heat shock (4 h, 28–36°C), Hsp70 expression was greater for RC than BPC beetles. In RC, field-collected beetles homozygous for PGI 1–1 had higher Hsp70 levels than heterozygotes or a 4–4 homozygote. These results reveal functional and physiological differences among PGI genotypes, which suggest that montane populations of this beetle are locally adapted to temperature.

Human sulfite oxidase R160Q: Identification of the mutation in a sulfite oxidase-deficient patient and expression and characterization of the mutant enzyme

Sulfite oxidase catalyzes the terminal reaction in the degradation of sulfur amino acids. Genetic deficiency of sulfite oxidase results in neurological abnormalities and often leads to death at an early age. The mutation in the sulfite oxidase gene responsible for sulfite oxidase deficiency in a 5-year-old girl was identified by sequence analysis of cDNA obtained from fibroblast mRNA to be a guanine to adenine transition at nucleotide 479 resulting in the amino acid substitution of Arg-160 to Gln. Recombinant protein containing the R160Q mutation was expressed in Escherichia coli, purified, and characterized. The mutant protein contained its full complement of molybdenum and heme, but exhibited 2% of native activity under standard assay conditions. Absorption spectroscopy of the isolated molybdenum domains of native sulfite oxidase and of the R160Q mutant showed significant differences in the 480- and 350-nm absorption bands, suggestive of altered geometry at the molybdenum center. Kinetic analysis of the R160Q protein showed an increase in Km for sulfite combined with a decrease in kcat resulting in a decrease of nearly 1,000-fold in the apparent second-order rate constant kcat/Km. Kinetic parameters for the in vitro generated R160K mutant were found to be intermediate in value between those of the native protein and the R160Q mutant. Native sulfite oxidase was rapidly inactivated by phenylglyoxal, yielding a modified protein with kinetic parameters mimicking those of the R160Q mutant. It is proposed that Arg-160 attracts the anionic substrate sulfite to the binding site near the molybdenum.

Diversity and phylogeny of gephyrin: Tissue-specific splice variants, gene structure, and sequence similarities to molybdenum cofactor-synthesizing and cytoskeleton-associated proteins

Gephyrin is essential for both the postsynaptic localization of inhibitory neurotransmitter receptors in the central nervous system and the biosynthesis of the molybdenum cofactor (Moco) in different peripheral organs. Several alternatively spliced gephyrin transcripts have been identified in rat brain that differ in their 5′ coding regions. Here, we describe gephyrin splice variants that are differentially expressed in non-neuronal tissues and different regions of the adult mouse brain. Analysis of the murine gephyrin gene indicates a highly mosaic organization, with eight of its 29 exons corresponding to the alternatively spliced regions identified by cDNA sequencing. The N- and C-terminal domains of gephyrin encoded by exons 3–7 and 16–29, respectively, display sequence similarities to bacterial, invertebrate, and plant proteins involved in Moco biosynthesis, whereas the central exons 8, 13, and 14 encode motifs that may mediate oligomerization and tubulin binding. Our data are consistent with gephyrin having evolved from a Moco biosynthetic protein by insertion of protein interaction sequences.

Quantitative assessment of enzyme specificity in vivo: P2 recognition by Kex2 protease defined in a genetic system

The specificity of the yeast proprotein-processing Kex2 protease was examined in vivo by using a sensitive, quantitative assay. A truncated prepro-α-factor gene encoding an α-factor precursor with a single α-factor repeat was constructed with restriction sites for cassette mutagenesis flanking the single Kex2 cleavage site (-SLDKR↓EAEA-). All of the 19 substitutions for the Lys (P2) residue in the cleavage site were made. The wild-type and mutant precursors were expressed in a yeast strain lacking the chromosomal genes encoding Kex2 and prepro-α-factor. Cleavage of the 20 sites by Kex2, expressed at the wild-type level, was assessed by using a quantitative-mating assay with an effective range greater than six orders of magnitude. All substitutions for Lys at P2 decreased mating, from 2-fold for Arg to >106-fold for Trp. Eviction of the Kex2-encoding plasmid indicated that cleavage of mutant sites by other cellular proteases was not a complicating factor. Mating efficiencies of strains expressing the mutant precursors correlated well with the specificity (kcat/KM) of purified Kex2 for comparable model peptide substrates, validating the in vivo approach as a quantitative method. The results support the conclusion that KM, which is heavily influenced by the nature of the P2 residue, is a major determinant of cleavage efficiency in vivo. P2 preference followed the rank order: Lys > Arg > Thr > Pro > Glu > Ile > Ser > Ala > Asn > Val > Cys > AsP > Gln > Gly > His > Met > Leu > Tyr > Phe > Trp.

Evolution of an enzyme active site: The structure of a new crystal form of muconate lactonizing enzyme compared with mandelate racemase and enolase

Muconate lactonizing enzyme (MLE), a component of the β-ketoadipate pathway of Pseudomonas putida, is a member of a family of related enzymes (the “enolase superfamily”) that catalyze the abstraction of the α-proton of a carboxylic acid in the context of different overall reactions. New untwinned crystal forms of MLE were obtained, one of which diffracts to better than 2.0-Å resolution. The packing of the octameric enzyme in this crystal form is unusual, because the asymmetric unit contains three subunits. The structure of MLE presented here contains no bound metal ion, but is very similar to a recently determined Mn2+-bound structure. Thus, absence of the metal ion does not perturb the structure of the active site. The structures of enolase, mandelate racemase, and MLE were superimposed. A comparison of metal ligands suggests that enolase may retain some characteristics of the ancestor of this enzyme family. Comparison of other residues involved in catalysis indicates two unusual patterns of conservation: (i) that the position of catalytic atoms remains constant, although the residues that contain them are located at different points in the protein fold; and (ii) that the positions of catalytic residues in the protein scaffold are conserved, whereas their identities and roles in catalysis vary.

The domain structure of protoporphyrinogen oxidase, the molecular target of diphenyl ether-type herbicides

Protoporphyrinogen oxidase (EC 1–3-3–4), the 60-kDa membrane-bound flavoenzyme that catalyzes the final reaction of the common branch of the heme and chlorophyll biosynthesis pathways in plants, is the molecular target of diphenyl ether-type herbicides. It is highly resistant to proteases (trypsin, endoproteinase Glu-C, or carboxypeptidases A, B, and Y), because the protein is folded into an extremely compact form. Trypsin maps of the native purified and membrane-bound yeast protoporphyrinogen oxidase show that this basic enzyme (pI > 8.5) was cleaved at a single site under nondenaturing conditions, generating two peptides with relative molecular masses of 30,000 and 35,000. The endoproteinase Glu-C also cleaved the protein into two peptides with similar masses, and there was no additional cleavage site under mild denaturing conditions. N-terminal peptide sequence analysis of the proteolytic (trypsin and endoproteinase Glu-C) peptides showed that both cleavage sites were located in putative connecting loop between the N-terminal domain (25 kDa) with the βαβ ADP-binding fold and the C-terminal domain (35 kDa), which possibly is involved in the binding of the isoalloxazine moiety of the FAD cofactor. The peptides remained strongly associated and fully active with the Km for protoporphyrinogen and the Ki for various inhibitors, diphenyl-ethers, or diphenyleneiodonium derivatives, identical to those measured for the native enzyme. However, the enzyme activity of the peptides was much more susceptible to thermal denaturation than that of the native protein. Only the C-terminal domain of protoporphyrinogen oxidase was labeled specifically in active site-directed photoaffinity-labeling experiments. Trypsin may have caused intramolecular transfer of the labeled group to reactive components of the N-terminal domain, resulting in nonspecific labeling. We suggest that the active site of protoporphyrinogen oxidase is in the C-terminal domain of the protein, at the interface between the C- and N-terminal domains.

The reductive enzyme thioredoxin 1 acts as an oxidant when it is exported to the Escherichia coli periplasm

Thioredoxin 1 is a major thiol-disulfide oxidoreductase in the cytoplasm of Escherichia coli. One of its functions is presumed to be the reduction of the disulfide bond in the active site of the essential enzyme ribonucleotide reductase. Thioredoxin 1 is kept in a reduced state by thioredoxin reductase. In a thioredoxin reductase null mutant however, most of thioredoxin 1 is in the oxidized form; recent reports have suggested that this oxidized form might promote disulfide bond formation in vivo. In the Escherichia coli periplasm, the protein disulfide isomerase DsbC is maintained in the reduced and active state by the membrane protein DsbD. In a dsbD null mutant, DsbC accumulates in the oxidized form. This oxidized form is then able to promote disulfide bond formation. In both these cases, the inversion of the function of these thiol oxidoreductases appears to be due to an altered redox balance of the environment in which they find themselves. Here, we show that thioredoxin 1 attached to the alkaline phosphatase signal sequence can be exported into the E. coli periplasm. In this new environment for thioredoxin 1, we show that thioredoxin 1 can promote disulfide bond formation and, therefore, partially complement a dsbA strain defective for disulfide bond formation. Thus, we provide evidence that by changing the location of thioredoxin 1 from cytoplasm to periplasm, we change its function from a reductant to an oxidant. We conclude that the in vivo redox function of thioredoxin 1 depends on the redox environment in which it is localized.

Studies on the enzymatic and transcriptional activity of the dimerization cofactor for hepatocyte nuclear factor 1

The relationship between the enzymatic and the transcriptional activity of the bifunctional protein pterin-4a-carbinolamine dehydratase/dimerization cofactor for hepatocyte nuclear factor 1 (DCoH) has been elucidated by site-directed mutagenesis. DCoH dimers harbor a binding site for hepatocyte nuclear factor 1 (HNF1), two active centers that bind pterins, and a saddle-shaped surface that resembles nucleic acid binding domains. Two domains of the protein have been selectively targeted to determine if a change in one activity affects the other. No strong correlation has been found, supporting the idea that carbinolamine dehydratase activity is not required for HNF1 binding in vitro or transcriptional coactivation in vivo. Double mutations in the active center, however, influence the in vivo transcriptional activity but not HNF1 binding. This finding suggests that some active center residues also are used during transcription, possibly for binding of another (macro)molecule. Several mutations in the saddle led to a surprising increase in transcription, therefore linking this domain to transcriptional regulation as well. The transcriptional function of DCoH therefore is composed of two parts, HNF1 binding and another contributing effect that involves the active site and, indirectly, the saddle.

B12-dependent ribonucleotide reductases from deeply rooted eubacteria are structurally related to the aerobic enzyme from Escherichia coli

The ribonucleotide reductases from three ancient eubacteria, the hyperthermophilic Thermotoga maritima (TM), the radioresistant Deinococcus radiodurans (DR), and the thermophilic photosynthetic Chloroflexus aurantiacus, were found to be coenzyme-B12 (class II) enzymes, similar to the earlier described reductases from the archaebacteria Thermoplasma acidophila and Pyrococcus furiosus. Reduction of CDP by the purified TM and DR enzymes requires adenosylcobalamin and DTT. dATP is a positive allosteric effector, but stimulation of the TM enzyme only occurs close to the temperature optimum of 80–90°C. The TM and DR genes were cloned by PCR from peptide sequence information. The TM gene was sequenced completely and expressed in Escherichia coli. The deduced amino acid sequences of the two eubacterial enzymes are homologous to those of the archaebacteria. They can also be aligned to the sequence of the large protein of the aerobic E. coli ribonucleotide reductase that belongs to a different class (class I), which is not dependent on B12. Structure determinations of the E. coli reductase complexed with substrate and allosteric effectors earlier demonstrated a 10-stranded β/α-barrel in the active site. From the conservation of substrate- and effector-binding residues we propose that the B12-dependent class II enzymes contain a similar barrel.

Epilepsy in mice deficient in the 65-kDa isoform of glutamic acid decarboxylase

γ-Aminobutyric acid (GABA), the major inhibitory neurotransmitter in the mammalian brain, is synthesized by two glutamate decarboxylase isoforms, GAD65 and GAD67. The separate role of the two isoforms is unknown, but differences in saturation with cofactor and subcellular localization suggest that GAD65 may provide reserve pools of GABA for regulation of inhibitory neurotransmission. We have disrupted the gene encoding GAD65 and backcrossed the mutation into the C57BL/6 strain of mice. In contrast to GAD67−/− animals, which are born with developmental abnormalities and die shortly after birth, GAD65−/− mice appear normal at birth. Basal GABA levels and holo-GAD activity are normal, but the pyridoxal 5′ phosphate-inducible apo-enzyme reservoir is significantly decreased. GAD65−/− mice develop spontaneous seizures that result in increased mortality. Seizures can be precipitated by fear or mild stress. Seizure susceptibility is dramatically increased in GAD65−/− mice backcrossed into a second genetic background, the nonobese diabetic (NOD/LtJ) strain of mice enabling electroencephalogram analysis of the seizures. The generally higher basal brain GABA levels in this backcross are significantly decreased by the GAD65−/− mutation, suggesting that the relative contribution of GABA synthesized by GAD65 to total brain GABA levels is genetically determined. Seizure-associated c-fos-like immunoreactivity reveals the involvement of limbic regions of the brain. These data suggest that GABA synthesized by GAD65 is important in the dynamic regulation of neural network excitability, implicate at least one modifier locus in the NOD/LtJ strain, and present GAD65−/− animals as a model of epilepsy involving GABA-ergic pathways.

Interaction of the herbicide glyphosate with its target enzyme 5-enolpyruvylshikimate 3-phosphate synthase in atomic detail

Biosynthesis of aromatic amino acids in plants, many bacteria, and microbes relies on the enzyme 5-enolpyruvylshikimate 3-phosphate (EPSP) synthase, a prime target for drugs and herbicides. We have identified the interaction of EPSP synthase with one of its two substrates (shikimate 3-phosphate) and with the widely used herbicide glyphosate by x-ray crystallography. The two-domain enzyme closes on ligand binding, thereby forming the active site in the interdomain cleft. Glyphosate appears to occupy the binding site of the second substrate of EPSP synthase (phosphoenol pyruvate), mimicking an intermediate state of the ternary enzyme⋅substrates complex. The elucidation of the active site of EPSP synthase and especially of the binding pattern of glyphosate provides a valuable roadmap for engineering new herbicides and herbicide-resistant crops, as well as new antibiotic and antiparasitic drugs.

Telomeric repeat amplification, without shortening or lengthening of the telomerase products: a method to analyze the processivity of telomerase enzyme

The Telomeric Repeat Amplification Protocol (TRAP) and its modified versions (including ours, TP-TRAP) change the size and/or the ratio of the telomerase products in the amplification stage of the assay. Based on our recently published method we developed a new TRAP. This method ensures that the number of telomeric repeats present in the original telomerase products does not change on PCR amplification. The usefulness of the method was proved with amplification of chemically synthesized telomerase products and a newly designed telomerase substrate oligonucleotide. This is the first report in which the PCR products directly reflect the size distribution of telomerase products generated by the enzyme.