Search

BackgroundDNA sequencing of archived dried blood spots (DBS) collected by newborn screening programs constitutes a potential health resource to study newborn disorders and understand genotype-phenotype relationships. However, its essential to verify that sequencing reads from DBS derived DNA are suitable for variant discovery.ResultsWe explored 16 metrics to comprehensively assess the quality of sequencing reads from 180 DBS and 35 whole blood (WB) samples. These metrics were used to assess a) mapping of reads to the reference genome, b) degree of DNA damage, and c) variant calling. Reads from both sets mapped with similar efficiencies, had similar overall DNA damage rates, measured by the mismatch rate with the reference genome, and produced variant calls sets with similar Transition-Transversion ratios. While evaluating single nucleotide changes that may have arisen from DNA damage, we observed that the A>T and T>A changes were more frequent in DNA from DBS than from WB. However, this did not affect the accuracy of variant calling, with DBS samples yielding a comparable count of high quality SNVs and indels in samples with at least 50x coverage.ConclusionsOverall, DBS DNA provided exome sequencing data of sufficient quality for clinical interpretation.

The Nudix homology clan encompasses over 80,000 protein domains from all three domains of life, defined by homology to each other. Proteins with a domain from this clan fall into four general functional classes: pyrophosphohydrolases, isopentenyl diphosphate isomerases (IDIs), adenine/guanine mismatch-specific adenine glycosylases (A/G-specific adenine glycosylases), and nonenzymatic activities such as protein/protein interaction and transcriptional regulation. The largest group, pyrophosphohydrolases, encompasses more than 100 distinct hydrolase specificities. To understand the evolution of this vast number of activities, we assembled and analyzed experimental and structural data for 205 Nudix proteins collected from the literature. We corrected erroneous functions or provided more appropriate descriptions for 53 annotations described in the Gene Ontology Annotation database in this family, and propose 275 new experimentally-based annotations. We manually constructed a structure-guided sequence alignment of 78 Nudix proteins. Using the structural alignment as a seed, we then made an alignment of 347 "select" Nudix homology domains, curated from structurally determined, functionally characterized, or phylogenetically important Nudix domains. Based on our review of Nudix pyrophosphohydrolase structures and specificities, we further analyzed a loop region downstream of the Nudix hydrolase motif previously shown to contact the substrate molecule and possess known functional motifs. This loop region provides a potential structural basis for the functional radiation and evolution of substrate specificity within the hydrolase family. Finally, phylogenetic analyses of the 347 select protein domains and of the complete Nudix homology clan revealed general monophyly with regard to function and a few instances of probable homoplasy. Proteins 2017; 85:775-811. © 2016 Wiley Periodicals, Inc.

Over the past thirty years, site-directed mutagenesis has become established as one of the most powerful techniques to probe enzyme reaction mechanisms1-3. Substitutions of active site residues are most likely to yield significant perturbations in kinetic parameters, but there are many examples of profound changes in these values elicited by remote mutations4-6. Ortholog comparisons of extant sequences show that many mutations do not have profound influence on enzyme function. As the number of potential single natural amino acid substitutions that can be introduced in a protein ofNamino acids in length by directed mutation is very large (19 * N), it would be useful to have a method to predict which amino acid substitutions are more likely to introduce significant changes in kinetic parameters in order to design meaningful probes into enzyme function. What is especially desirable is the identification of critical residues that do not contact the substrate directly, and may be remote from the active site.We collected literature data reflecting the effects of 2,804 mutations on kinetic properties for 12 enzymes. These data along with characteristic predictors were used in a machine-learning scheme to train a classifier to predict the effect of mutation. Use of this algorithm allows one to predict with a 2.5-fold increase in precision, if a given mutation, made anywhere in the enzyme, will cause a decrease in kcat/Kmvalue of ≥ 95%. The improved precision allows the experimentalist to reduce the number of mutations necessary to probe the enzyme reaction mechanism.

The subfamily Iα aminotransferases are typically categorized as having narrow specificity toward carboxylic amino acids (AATases), or broad specificity that includes aromatic amino acid substrates (TATases). Because of their general role in central metabolism and, more specifically, their association with liver-related diseases in humans, this subfamily is biologically interesting. The substrate specificities for only a few members of this subfamily have been reported, and the reliable prediction of substrate specificity from protein sequence has remained elusive. In this study, a diverse set of aminotransferases was chosen for characterization based on a scoring system that measures the sequence divergence of the active site. The enzymes that were experimentally characterized include both narrow-specificity AATases and broad-specificity TATases, as well as AATases with broader-specificity and TATases with narrower-specificity than the previously known family members. Molecular function and phylogenetic analyses underscored the complexity of this family's evolution as the TATase function does not follow a single evolutionary thread, but rather appears independently multiple times during the evolution of the subfamily. The additional functional characterizations described in this article, alongside a detailed sequence and phylogenetic analysis, provide some novel clues to understanding the evolutionary mechanisms at work in this family.

The nearly 50,000 known Nudix proteins have a diverse array of functions, of which the most extensively studied is the catalyzed hydrolysis of aberrant nucleotide triphosphates. The functions of 171 Nudix proteins have been characterized to some degree, although physiological relevance of the assayed activities has not always been conclusively demonstrated. We investigated substrate specificity for eight structurally characterized Nudix proteins, whose functions were unknown. These proteins were screened for hydrolase activity against a 74‐compound library of known Nudix enzyme substrates. We found substrates for four enzymes with k cat/K m values >10,000 M−1 s−1: Q92EH0_LISIN of Listeria innocua serovar 6a against ADP‐ribose, Q5LBB1_BACFN of Bacillus fragilis against 5‐Me‐CTP, and Q0TTC5_CLOP1 and Q0TS82_CLOP1 of Clostridium perfringens against 8‐oxo‐dATP and 3'‐dGTP, respectively. To ascertain whether these identified substrates were physiologically relevant, we surveyed all reported Nudix hydrolytic activities against NTPs. Twenty‐two Nudix enzymes are reported to have activity against canonical NTPs. With a single exception, we find that the reported k cat/K m values exhibited against these canonical substrates are well under 105 M−1 s−1. By contrast, several Nudix enzymes show much larger k cat/K m values (in the range of 105 to >107 M−1 s−1) against noncanonical NTPs. We therefore conclude that hydrolytic activities exhibited by these enzymes against canonical NTPs are not likely their physiological function, but rather the result of unavoidable collateral damage occasioned by the enzymes' inability to distinguish completely between similar substrate structures. Proteins 2016; 84:1810–1822. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

The existence of low levels of intersubunit communication in homooligomeric enzymes is often difficult to discover, as the identical active sites cannot be probed individually to dissect their interdependent contributions. The homodimeric paralogs, E. coli aspartate- (AATase) and tyrosine aminotransferase (TATase), have not been demonstrated to show allostery. To address this question, we engineered a hybrid aminotransferase containing two distinct catalytic pockets: an AATase and a TATase site. The TATase/AATase hybrid was constructed by grafting an engineered TATase active site into one of the catalytic pockets of E. coli AATase. Each active site conserves its specific catalytic and inhibitor binding properties, and the hybrid catalyzes simultaneously each aminotransferase reaction at the respective site. Importantly, association of a selective inhibitor into one of the catalytic pockets decreases the activity of the second active site by up to 25%, thus proving unequivocally the existence of allosteric communication between active sites. The procedure may be applicable to other homologous sets of enzymes.

Establishing a quantitative understanding of the determinants of affinity in protein-protein interactions remains challenging. For example, TEM-1/β-lactamase inhibitor protein (BLIP) and SHV-1/BLIP are homologous β-lactamase/β-lactamase inhibitor protein complexes with disparate K(d) values (3 nM and 2 μM, respectively), and a single substitution, D104E in SHV-1, results in a 1000-fold enhancement in binding affinity. In TEM-1, E104 participates in a salt bridge with BLIP K74, whereas the corresponding SHV-1 D104 does not in the wild type SHV-1/BLIP co-structure. Here, we present a 1.6 A crystal structure of the SHV-1 D104E/BLIP complex that demonstrates that this point mutation restores this salt bridge. Additionally, mutation of a neighboring residue, BLIP E73M, results in salt bridge formation between SHV-1 D104 and BLIP K74 and a 400-fold increase in binding affinity. To understand how this salt bridge contributes to complex affinity, the cooperativity between the E/K or D/K salt bridge pair and a neighboring hot spot residue (BLIP F142) was investigated using double mutant cycle analyses in the background of the E73M mutation. We find that BLIP F142 cooperatively stabilizes both interactions, illustrating how a single mutation at a hot spot position can drive large perturbations in interface stability and specificity through a cooperative interaction network.

The field of enzymology was profoundly affected by the death of Bill Jencks on 3 January 2007, after a 17-year struggle against Alzheimer’s disease. Jencks was a major innovator, a superb teacher, a devoted husband and father, and a helpful friend and colleague. Those of us who had the privilege of working with Bill, or of knowing him professionally, are poorer for his passing.

β-lactamase/β-lactamase Inhibitor Protein (BLIP) complexes are emerging as a well characterized experimental model system for studying protein-protein interactions. β-lactamases are enzymes that catalyze the hydrolysis of β-lactam antibiotics. BLIP is a 165 amino acid protein that inhibits several class A β-lactamases with a wide range of affinities: pM affinity for K1; nM affinity for TEM-1, SME-1, and BlaI but only µM affinity for SHV-1 β-lactamase. The large differences in affinity coupled with the availability of extensive mutagenesis data and high resolution crystal structures for the TEM-1/BLIP and SHV-1/BLIP complexes make them attractive systems for the further development of protein engineering and computational design methodologies. We used EGAD, a physics-based computational design program, to redesign BLIP with the goal of increasing affinity for SHV-1. The resulting designed sequences are highly similar to wildtype, with the exception of BLIP residues surrounding β-lactamase position 104. Interestingly, this residue is a known specificity determinant between TEM-1 and SHV-1. Characterization of several of the designs and point mutants revealed that in all cases, the mutations stabilize the interface by 10 to 1000 fold relative to wildtype BLIP. The calculated changes in binding affinity for the mutants were within a mean absolute error of 0.87 kcal/mol from the experimental values, and comparison of calculated and experimental values for a set of 30 SHV-1/BLIP complexes yielded a correlation coefficient of 0.77. Although binding specificity for SHV-1 versus TEM-1 was not explicitly considered in the design process, two of the BLIP variants exhibit a small specificity switch. Structures of the two highest affinity complexes, SHV-1/BLIP (E73M) and SHV-1/BLIP (E73M, S130K, S146M), are presented at 1.7 Å resolution. While the predicted structures have much in common with the experimentally determined structures, they do not coincide perfectly; in particular a salt bridge between SHV-1 D104 and BLIP K74 is observed in the experimental structures, but not in the predicted design conformations. This discrepancy highlights the difficulty of modeling salt bridge interactions with a protein design algorithm that approximates sidechains as discrete rotamers. Nevertheless, while local structural features of the interface were sometimes miscalculated, EGAD is globally successful in designing complexes with increased affinity.

Aminotransferases are essential enzymes involved in the central metabolism of all organisms. The Iα subfamily of aspartate and tyrosine aminotransferases (AATases and TATases) is the best-characterized grouping, but only eight enzymes from this subfamily, representing relatively little sequence diversity, have been experimentally characterized for substrate specificity ( i.e ., AATase vs . TATase). Genome annotation, based on this limited dataset, provides tentative assignments for all sequenced members of this subfamily. This procedure is, however, subject to error, particularly when the experimental basis set is limited. To address this problem we cloned twelve additional subfamily Iα enzymes from an evolutionarily divergent set of organisms. Nine were purified to homogeneity after heterologous expression in Escherichia coli in native, intein-tagged or His 6 -tagged forms. The two Saccharomyces cerevisiae isoforms were recombinantly produced in yeast. The effects of the C-terminal tags on expression, purification and enzyme activity are discussed.

Five positions in the Escherichia coli malate dehydrogenase (eMDH) sequence, which distinguish MDH from lactate dehydrogenase (LDH) activity, were identified through a combination of Venn diagrams constructed from whole genomic data and from unbiased representative sequences from terminal clades. Incorporation of the five changes in eMDH sufficed to convert the enzyme from one with ( k cat / K m pyruvate )/( k cat / K m oxaloacetate ) = 6.1 × 10 −9 to one with that ratio = 28. The substrate specificity was thus changed by a factor of 4.6 × 10 9 . The k cat / K m pyruvate value for the pentamutant (eMDH I12V/R81Q/M85E/G210A/V214I) is 3,500 M −1 ·s −1 , which is ≈1/1,000 of the values found for typical wild-type LDHs. The procedure isolates an intersection of “strong forcing sets” that should prove to be of general use in switching paralog function.

While the urea-mediated unfolding pathway of the Escherichia coli aspartate aminotransferase (eAATase) homodimer proceeds through a reversible three-state process with a partially folded dimeric intermediate, D D* 2U (E. Deu and J. F. Kirsch, accompanying paper), that of a cofactor-stabilized form differs. Pyridoxal phosphate, which binds at the intersubunit active sites, stabilizes the native form by 6 kcal mol-1 and dissociates during the D==D* transition. Reductive trapping of the cofactor to a nondissociable derivative (PPL-eAATase) precludes the formation of D*. A novel monomeric intermediate (M'-PPL) with 70% of the native secondary structure (circular dichroism) was identified in the unfolding pathway of PPL-eAATase: D-PPL2==2M'-PPL==2U-PPL. The combined results define two structural regions with distinct stabilities: the active site region (ASR) and the generally more stable, dimerization region (DMR). The DMR includes the key intersubunit contacts. It is responsible for the multimeric nature of D*, and its disorder leads to dimer dissociation. Selective strengthening of the ASR-cofactor interactions by cofactor trapping reverses the relative stabilities of the two regions (from DMRASR in the apoenzyme to ASRDMR in PPL-eAATase) and results in a reordering of the eAATase denaturation pathway.

The advantages of electrospray ionization mass spectrometry (ESIMS) to measure relative solution-phase affinities of tightly bound protein-protein complexes are demonstrated with selected variants of the Bacillus amyloliquefaciens protein barstar (b*) and the RNAase barnase (bn), which form protein-protein complexes with a range of picomolar to nanomolar dissociation constants. A novel chemical annealing procedure rapidly establishes equilibrium in solutions containing competing b* variants with limiting bn. The relative ion abundances of the complexes and those of the competing unbound monomers are shown to reflect the relative solution-phase concentrations of those respective species. No measurable dissociation of the complexes occurs either during ESI or mass detection, nor is there any evidence for nonspecific binding at protein concentrations25 microM. Differences in DeltaDeltaG of dissociation between variants were determined with precisions0.1 kcal/mol. The DeltaDeltaG values obtained deviate on average by 0.26 kcal/mol from those measured with a solution-phase enzyme assay. It is demonstrated that information about the protein conformation and covalent modifications can be obtained from differences in mass and charge state distributions. This method serves as a rapid and precise means to interrogate protein-protein-binding surfaces for complexes that have affinities in the picomolar to nanomolar range.

Human tyrosine aminotransferase (hTATase) is the pyridoxal phosphate-dependent enzyme that catalyzes the reversible transamination of tyrosine to p-hydrophenylpyruvate, an important step in tyrosine metabolism. hTATase deficiency is implicated in the rare metabolic disorder, tyrosinemia type II. This enzyme is a member of the poorly characterized Iγ subfamily of the family I aminotransferases. The full length and truncated forms of recombinant hTATase were expressed in Escherichia coli, and purified to homogeneity. The pH-dependent titration of wild-type reveals a spectrum characteristic of family I aminotransferases with an aldimine pKa of 7.22. I249A mutant hTATase exhibits an unusual spectrum with a similar aldimine pKa (6.85). hTATase has very narrow substrate specificity with the highest enzymatic activity for the Tyr/α-ketoglutarate substrate pair, which gives a steady state kcat value of 83 s−1. In contrast there is no detectable transamination of aspartate or other cosubstrates. The present findings show that hTATase is the only known aminotransferase that discriminates significantly between Tyr and Phe: the kcat/Km value for Tyr is about four orders of magnitude greater than that for Phe. A comparison of substrate specificities of representative Iα and Iγ aminotransferases is described along with the physiological significance of the discrimination between Tyr and Phe by hTATase as applied to the understanding of the molecular basis of phenylketonuria.

OmpT is a surface protease of gram-negative bacteria that has been shown to cleave antimicrobial peptides, activate human plasminogen, and degrade some recombinant heterologous proteins. We have analyzed the substrate specificity of OmpT by two complementary substrate filamentous phage display methods: (i) in situ cleavage of phage that display protease-susceptible peptides by Escherichia coli expressing OmpT and (ii) in vitro cleavage of phage-displayed peptides using purified enzyme. Consistent with previous reports, OmpT was found to exhibit a virtual requirement for Arg in the P1 position and a slightly less stringent preference for this residue in the P1′ position (P1 and P1′ are the residues immediately prior to and following the scissile bond). Lys, Gly, and Val were also found in the P1′ position. The most common residues in the P2′ position were Val or Ala, and the P3 and P4 positions exhibited a preference for Trp or Arg. Synthetic peptides based upon sequences selected by bacteriophage display were cleaved very efficiently, with k cat / K m values up to 7.3 × 10 6 M −1 s −1 . In contrast, a peptide corresponding to the cleavage site of human plasminogen was hydrolyzed with a k cat / K m almost 10 6 -fold lower. Overall, the results presented in this work indicate that in addition to the P1 and P1′ positions, additional amino acids within a six-residue window (between P4 and P2′) contribute to the binding of substrate polypeptides to the OmpT binding site.

The crystal structure of 1-aminocyclopropane-1-carboxylate (ACC) synthase in complex with the substrate analogue [2-(amino-oxy)ethyl](5′-deoxyadenosin-5′-yl)(methyl)sulfonium (AMA) was determined at 2.01-A resolution. The crystallographic results show that a covalent adduct (oxime) is formed between AMA (an amino-oxy analogue of the natural substrate S -adenosyl- l -methionine (SAM)) and the pyridoxal 5′-phosphate (PLP) cofactor of ACC synthase. The oxime formation is supported by spectroscopic data. The ACC synthase–AMA structure provides reliable and detailed information on the binding mode of the natural substrate of ACC synthase and complements previous structural and functional work on this enzyme.

Cystathionine beta-synthase (CBS) is a pyridoxal-phosphate-dependent enzyme that catalyzes a beta-replacement reaction in which the hydroxyl group of serine (L-Ser) is displaced by the thiol of homocysteine (L-Hcys) to form cystathionine (L-Cth) in the first step of the trans-sulfuration pathway. A new continuous assay for the forward reaction, employing cystathionine beta-lyase and L-lactate dehydrogenase as coupling enzymes, is described. It alleviates product inhibition by L-Cth and revealed that the values for (1.2 mM) and for substrate inhibition by L-Hcys ( = 2.0 mM) are lower than those previously reported. A continuous, 5,5'-dithio-bis-(2-nitrobenzoic acid) (DTNB)-based assay for the CBS-catalyzed hydrolysis of L-Cth to L-Ser and L-Hcys provides a tool for investigation of the reverse reaction (k(catR) = 0.56 s(-)(1), = 0.083 mM). The (catR)/ versus pH profile of ytCBS is bell-shaped with a pH optimum of 8.3, and the pK(a) values for the acidic and basic limbs are 8.05 and 8.63, respectively. The latter is assigned to the alpha-amino group of L-Cth (pK(a) = 8.54). The internal aldimine of ytCBS remains protonated at pH11; therefore, the acidic pK(a) is assigned to an enzyme functionality that is not associated with the internal aldimine. K(eq) was determined directly and from the kinetic parameters, and the values are 0.61 and 1.2 microM, respectively.

The 1.6-A crystal structure of the covalent ketimine complex of apple 1-aminocyclopropane-1-carboxylate (ACC) synthase with the potent inhibitor l-aminoethoxyvinylglycine (AVG) is described. ACC synthase catalyzes the committed step in the biosynthesis of ethylene, a plant hormone that is responsible for the initiation of fruit ripening and for regulating many other developmental processes. AVG is widely used in plant physiology studies to inhibit the activity of ACC synthase. The structural assignment is supported by the fact that the complex absorbs maximally at 341 nm. These results are not in accord with the recently reported crystal structure of the tomato ACC synthase AVG complex, which claims that the inhibitor only associates noncovalently. The rate constant for the association of AVG with apple ACC synthase was determined by stopped-flow spectrophotometry (2.1 x 10(5) m(-1) s(-1)) and by the rate of loss of enzyme activity (1.1 x 10(5) m(-1) s(-1)). The dissociation rate constant determined by activity recovery is 2.4 x 10(-6) s(-1). Thus, the calculated K(d) value is 10-20 pm.

The role of the Lys68*:Glu265 intersubunit salt bridge that is conserved (Csb) in all known aspartate aminotransferases (AATases), except those of animal cytosolic, Ac (His68*:Glu265), and plant mitochondrial, Pm (Met68*:Gln265), origins, was evaluated in the Escherichia coli AATase. Two double-mutant cycles, to K68M/E265Q and the charge reversed K68E/E265K, were characterized with the context dependence (C) and impact (I) formalism, previously defined for functional chimeric analysis. Mutations of Lys68* with Glu265 fixed are generally more deleterious than the converse mutations of Glu265 with Lys68* fixed, showing that buried negative charges have greater effects than buried positive charges in this context. Replacement of the charged Lys68*:Glu265 with the K68M/E265Q neutral pair introduces relatively small effects on the kinetic parameters. The differential sensitivity of k(cat)/K(M, L-Asp) and k(cat)/K(M, alpha-KG) to salt bridge mutagenic replacements is shown by a linear-free energy relationship, in which the logarithms of the latter second order rate constants are generally decreased by a factor of two more than are those of the former. Thus, k(cat)/K(M, L-Asp) and k(cat)/K(M, alpha-KG) are 133 and 442 mM(-1)s(-1) for the wild-type (WT) enzyme, respectively, but their relative order is reversed in the more severely compromised mutants (14.8 and 5.3 mM(-1)s(-1) for K68E). A Venn diagram illustrates apparent forced covariances of groups of amino acids that accompany the naturally occurring salt bridge replacements in the Pm and Ac classes. The more deeply rooted tree indicates that the Csb variant was the ancestral specie.

The active sites of the homologous pyridoxal phosphate- (PLP-) dependent enzymes 1 -aminocyclopropane- I -carboxylate (ACC) synthase and aspartate aminotransferase (AATase) are almost entirely conserved, yet the pK a 's ofthe two internal aldimines are 9.3 and 7.0, respectively, to complement the substrate pK a 's (S-adenosylmethionine pK a = 7.8 and aspartate pK a = 9.9). This complementation is required for maximum enzymatic activity in the physiological pH range. The most prominent structural difference in the active site is that Ile232 of ACC synthase is replaced by alanine in AATase. The I232A mutation was introduced into ACC synthase with a resulting 1.1 unit decrease (from 9.3 to 8.2) in the aldimine pK a , thus identifying Ile232 as a major determinant of the high pK a of ACC synthase. The mutation also resulted in reduced k c a t (0.5 vs 11 s - 1 ) and k c a t /K m values (5.0 x 10 4 vs 1.2 x 10 6 M - 1 s - 1 ). The effect of the mutation is interpreted as the result of shortening of the Tyr233-PLP hydrogen bond. Addition of the Y233F mutation to the I232A ACC synthase to generate the double mutant I232A/Y233F raised the pK a from 8.2 to 8.8, because the Y233F mutation eliminates the hydrogen bond between that residue and PLP. The introduction of the retro mutation A224I into AATase raised the aldimine pK a of that enzyme from 6.96 to 7.16 and resulted in a decrease in single-turnover k m a x (108 vs 900 s - 1 for aspartate) and k m a x /K m a p p (7.5 x 10 4 vs 3.8 x 10 5 M - 1 s - 1 ) values. The distance from the pyridine nitrogen of the cofactor to a conserved aspartate residue is 2.6 A in AATase and 3.8 A in ACC synthase. The D230E mutation introduced into ACC synthase to close this distance increases the aldimine pK a from 9.3 to 10.0, as would be predicted from a shortened hydrogen bond.

Glutamate 47 is conserved in 1-aminocyclopropane-1-carboxylate (ACC) synthases and is positioned near the sulfonium pole of (S,S)-S-adenosyl-L-methionine (SAM) in the modeled pyridoxal phosphate quinonoid complex with SAM. E47Q and E47D constructs of ACC synthase were made to investigate a putative ionic interaction between Glu47 and SAM. The k(cat)/K(m) values for the conversion of (S,S)-SAM to ACC and methylthioadenosine (MTA) are depressed 630- and 25-fold for the E47Q and E47D enzymes, respectively. The decreases in the specificity constants are due to reductions in k(cat) for both mutant enzymes, and a 5-fold increase in K(m) for the E47Q enzyme. Importantly, much smaller effects were observed for the kinetic parameters of reactions with the alternate substrates L-vinylglycine (L-VG) (deamination to form alpha-ketobutyrate and ammonia) and L-alanine (transamination to form pyruvate), which have uncharged side chains. L-VG is both a substrate and a mechanism-based inactivator of the enzyme [Feng, L., and Kirsch, J. F. (2000) Biochemistry 39, 2436-2444], but the partition ratio, k(cat)/k(inact), is unaffected by the Glu47 mutations. ACC synthase primarily catalyzes the beta,gamma-elimination of MTA from the (R,S) diastereomer of SAM to produce L-VG [Satoh, S., and Yang, S. F. (1989) Arch.Biochem. Biophys. 271, 107-112], but catalyzes the formation of ACC to a lesser extent via alpha,gamma-elimination of MTA. The partition ratios for (alpha,gamma/beta,gamma)-elimination on (R,S)-SAM are 0.4, < or =0.014, and < or =0.08 for the wild-type, E47Q, and E47D enzymes, respectively. The results of these experiments strongly support a role for Glu47 as an anchor for the sulfonium pole of (S,S)-SAM, and consequently a role as an active site determinant of reaction specificity.

The low-barrier hydrogen bond (LBHB) between the Asp and His residues of the catalytic triad in a serine protease was perturbed via the D32C mutation in subtilisin BPN' (Bacillus protease N'). This mutant enzyme catalyzes the hydrolysis of N-Suc-Ala-Ala-Pro-Phe-SBzl with a k(cat)/K(m) value that is only 8-fold reduced from that of the wild-type (WT) enzyme. The value of k(cat)/K(m) for the corresponding p-nitroanilide (pNA) substrate is only 50-fold lower than that of the WT enzyme (DeltaDeltaG++ = 2.2 kcal/mol). The pK(a) controlling the ascending limb of the pH versus k(cat)/K(m) profile is lowered from 7.01 (WT) to 6.53 (D32C), implying that any hydrogen bond replacing that between Asp32 and His64 of the WT enzyme most likely involves the neutral thiol rather than the thiolate form of Cys32. It is shown by viscosity variation that the reaction of WT subtilisin with N-Suc-Ala-Ala-Pro-Phe-SBzl is 50% (sucrose) to 100% (glycerol) diffusion-controlled, while that of the D32C construct is 29% (sucrose) to 76% (glycerol) diffusion-controlled. The low-field NMR resonance of 18 ppm that has been assigned to a proton shared by Asp32 and His64, and is considered diagnostic of a LBHB in the WT enzyme, is not present in D32C subtilisin. Thus, the LBHB is not an inherent requirement for substantial rate enhancement for subtilisin.

The sequences of genes encoding homologues of 1-aminocyclopropane-1-carboxylate (ACC) synthase, the first enzyme in the two-step biosynthetic pathway of the important plant hormone ethylene, have recently been found in Fugu rubripes and Homo sapiens (Peixoto et al., Gene 246 (2000) 275). ACC synthase (ACS) catalyzes the formation of ACC from S-adenosyl-L-methionine. ACC is oxidized to ethylene in the second and final step of ethylene biosynthesis. Profound physiological questions would be raised if it could be demonstrated that ACC is formed in animals, because there is no known function for ethylene in these organisms. We describe the cloning of the putative human ACS (PHACS) cDNA that encodes a 501 amino acid protein that exhibits 58% sequence identity to the putative Fugu ACS and approximately 30% sequence identity to plant ACSs. Purified recombinant PHACS, expressed in Pichia pastoris, contains bound pyridoxal-5'-phosphate (PLP), but does not catalyze the synthesis of ACC. PHACS does, however, catalyze the deamination of L-vinylglycine, a known side-reaction of apple ACS. Bioinformatic analysis indicates that PHACS is a member of the alpha-family of PLP-dependent enzymes. Molecular modeling data illustrate that the conservation of residues between PHACS and the plant ACSs is dispersed throughout its structure and that two active site residues that are important for ACS activity in plants are not conserved in PHACS.

Seven of the 13 non-glycine contact amino acids in the hen (chicken) egg white lysozyme (HEWL) epitope for antibody Fab-10 each contribute ≤0.3 kcal/mol to the change in free energy (ΔΔGD) from wild type (WT) when replaced by alanine (nullspots), and three others each give (0.7 4.0 kcal/mol on alanine substitution. Proteins 2000;40:49–57. © 2000 Wiley-Liss, Inc.

L-Vinylglycine (L-VG) has been shown to be a mechanism-based inhibitor of 1-aminocyclopropane-1-carboxylate (ACC) synthase [Satoh, S., and Yang, S. F. (1989) Plant Physiol. 91, 1036-1039] as well as of other pyridoxal phosphate-dependent enzymes. This report demonstrates that L-VG is primarily an alternative substrate for the enzyme. The L-VG deaminase activity of ACC synthase yields the products alpha-ketobutyrate and ammonia with a k(cat) value of 1.8 s(-1) and a K(m) value of 1.4 mM. The k(cat)/K(m) of 1300 M(-1) s(-1) is 0.17% that of the diffusion-controlled reaction with the preferred substrate, S-adenosyl-L-methionine. The enzyme-L-VG complex partitions to products 500 times for every inactivation event. The catalytic mechanism proceeds through a spectrophotometrically detected quinonoid with lambda(max) of 530 nm, which must rearrange to a 2-aminocrotonate aldimine to yield final products. Alternative mechanisms for the inactivation reaction are presented, and the observed kinetics for the full reaction course are satisfactorily modeled by kinetic simulation. The inactive enzyme is an aldimine with lambda(max) of 432 nm. It is resistant to NaBH(3)CN but is reduced by NaBH(4). ACC synthase is now expressed in Pichia pastoris with an improved yield of 10 mg/L.

The 2.4 A crystal structure of the vitamin B6-dependent enzyme 1-aminocyclopropane-1-carboxylate (ACC) synthase is described. This enzyme catalyses the committed step in the biosynthesis of ethylene, a plant hormone that is responsible for the initiation of fruit ripening and for regulating many other developmental processes. ACC synthase has 15 % sequence identity with the well-studied aspartate aminotransferase, and a completely different catalytic activity yet the overall folds and the active sites are very similar. The new structure together with available biochemical data enables a comparative mechanistic analysis that largely explains the catalytic roles of the conserved and non-conserved active site residues. An external aldimine reaction intermediate (external aldimine with ACC, i.e. with the product) has been modeled. The new structure provides a basis for the rational design of inhibitors with broad agricultural applications.

Alanine scanning mutagenesis of the HyHEL-10 paratope of the HyHEL-10/HEWL complex demonstrates that the energetically important side chains (hot spots) of both partners are in contact. A plot of deltadeltaG(HyHEL-10_mutant) vs. deltadeltaG(HEWL_mutant) for the five of six interacting side-chain hydrogen bonds is linear (Slope = 1). Only 3 of the 13 residues in the HEWL epitope contribute >4 kcal/mol to the free energy of formation of the complex when replaced by alanine, but 6 of the 12 HyHEL-10 paratope amino acids do. Double mutant cycle analysis of the single crystallographically identified salt bridge, D32H/K97, shows that there is a significant energetic penalty when either partner is replaced with a neutral side-chain amino acid, but the D32(H)N/K97M complex is as stable as the WT. The role of the disproportionately high number of Tyr residues in the CDR was evaluated by comparing the deltadeltaG values of the Tyr --> Phe vs. the corresponding Tyr --> Ala mutations. The nonpolar contacts in the light chain contribute only about one-half of the total deltadeltaG observed for the Tyr --> Ala mutation, while they are significantly more important in the heavy chain. Replacement of the N31L/K96 hydrogen bond with a salt bridge, N31D(L)/K96, destabilizes the complex by 1.4 kcal/mol. The free energy of interaction, deltadeltaG(int), obtained from double mutant cycle analysis showed that deltadeltaG(int) for any complex for which the HEWL residue probed is a major immunodeterminant is very close to the loss of free energy observed for the HyHEL-10 single mutant. Error propagation analysis of double mutant cycles shows that data of atypically high precision are required to use this method meaningfully, except where large deltadeltaG values are analyzed.

The rate constants, k(on), for the formation of hen (chicken) lysozyme (HEWL). Fab-10 complexes have been determined for wild-type (WT) and epitope-mutated lysozymes by a homogeneous solution method based on the 95 percent reduced enzymatic activity of the complex. The values fall within a narrow 10-fold range [(0.18 to 1.92) x 10(6)M(-1)s(-l)]. The affinity constants, K(D), cover a broader, 440-fold, range from 0.075 to 33 nM. Values of K(D) as high as 7 microM were obtained for the complexes prepared from some mutations at HEWL positions 96 and 97, but the associated kinetic constants could not be determined. The values of k(on) are negatively correlated with side-chain volume at position 101HEWL, but are essentially independent of this parameter for position 21HEWL substitutions. The multiple mutations made at positions 21HEWL and 101HEWL provide sufficient experimental data on complex formation to evaluate phi values [phi = (deltadeltaGon)/(deltadeltaG(D))] at these two positions to begin to define trajectories for protein-protein association. The data, when interpreted within the concept of a two-step association sequence embracing a metastable encounter complex intermediate, argue that the rate determining step at position 21HEWL (phiavg = 0.2) is encounter complex formation, but the larger phi(avg) value of 0.36 experienced formore » most position 101HEWL mutations indicates a larger contribution from the post-encounter annealing process at this site for these replacements.« less

The hen (chicken) egg-white lysozyme (HEWL) epitope for the monoclonal antibody HyHEL-10 Fab (Fab-10) was investigated by alanine scan mutagenesis. The association rate constants (k(on)) for the HEWL Fab-10 complexes were obtained from the homogenous solution method described in the preceding paper (Taylor et al., 1998). A new method for determining the dissociation rate constant (k(off)) for the complex, by trapping nascent free antibody with an inactive HEWL mutant is described. The values of k(on) fall within a factor of 2 of the wild-type (WT) HEWL value (1.43+/-0.13 X 10(6)M(-1)s(-1)), while the increases in k(off)more nearly reflect the total change in free energies of the complex (deltadeltaG(D)). The dissociation constants (K(D)) were measured directly in those cases where satisfactory kinetic data could not be obtained. The Y20A, K96A, and K97A HEWL.Fab-10 complexes are destabilized by more than 4 kcal/mol compared to the WT complex. The R21A, L75A, and D101A antibody complexes are moderately destabilized (0.7 < deltadeltaG(D)< or = 1.0 kcal/mol). Additional mutations of the "hotspot" residues (Tyr20, Lys96, Lys97) were constructed to probe, more precisely, the nature of their contributions to complex formation. The results show that the entire hydrocarbon side chains of Tyr20 and Lys97, and only the epsilon-amino group of Lys96, contribute to the stability of the complex. The value of deltadeltaG(D) for the R21A mutant complex is a distinct outlier in the Arg21 replacement series demonstrating the importance of supplementing alanine scan mutagenesis with additional mutations.

Aspartate aminotransferase (AATase) and aminocyclopropane carboxylate synthase (ACC synthase) are pyridoxal phosphate (PLP)-dependent enzymes whose common junction of mechanistic divergence is after the formation of a Cα carbanion from the amino acid substrate bound to PLP as a Schiff base (aldimine). AATase catalyzes the reversible interconversion of α-amino acids and α-keto acids, while ACC synthase effects the irreversible decomposition of S-adenosylmethionine (SAM) to 1-aminocyclopropane-1-carboxylate (ACC) and 5′-methylthioadenosine (MTA). ACC is subsequently converted to ethylene, the plant ripening and senescence hormone, by ACC oxidase, the next enzyme in the pathway. AATase and ACC synthase exhibit many similar phenomenological characteristics that result from different detailed mechanistic origins. The kcat/KMversus pH profiles for both enzymes are similar (AATase, acidic pKa = 6.9, basic pKa = 9.6; ACC synthase, acidic pKa = 7.5, basic pKa = 8.9); however the acidic pKa of AATase reflects the ionization of an enzyme proton from the internal Schiff base, and the basic one is that of the α-amino group of the substrate, while the opposite situation obtains for ACC synthase, i.e. the apparent pKa of 7.4 is due to the α-amino group of SAM, whereas that of 9 reflects the Schiff base pKa. The mechanistic imperative underlying this reversal is dictated by the reaction mechanism and the low pKa of the α-amino group of SAM. The low pKa of SAM requires that the enzyme pKa be moved upward in order to have sufficient quantities of the reacting species at neutral pH. It is shown by viscosity variation experiments with wild-type and active site mutant controls of both enzymes that the reaction of SAM with ACC synthase is 100% diffusion controlled (kcat/KM = 1.2 × 106 l mol−1 s−1) while the corresponding reaction for the combination of L-aspartate with AATase is insensitive to viscosity, and is therefore chemically not diffusion limited. Tyr225 (AATase) or Tyr233 (ACC synthase) forms a hydrogen bond with the PLP in both enzymes, but that formed with the former enzyme is stronger and accounts for the lower pKa of the Schiff base. © 1998 John Wiley & Sons, Ltd.

The primary role of Tyr225 in the aspartate aminotransferase mechanism is to provide a hydrogen bond to stabilize the 3'O- functionality of bound pyridoxal phosphate. The strength of this hydrogen bond is perturbed by replacement of Tyr225 with 3-fluoro-L-tyrosine (FlTyr) by in vitro transcription/translation. This mutant enzyme exhibits kcat/values that are near to those of wild type enzyme; however, the kcat/vs pH profile is much sharper with similar pKas of approximately 7.5 for both the ascending and descending limbs. The pKas are assigned to the endocyclic proton of the internal aldimine and to the bridging hydrogen bond, respectively. The pKas in the kcat vs pH profile of 7.2 and 8.7 are assigned to the epsilon-NH3+ of lysine 258 and to the endocyclic protons of the ketimine complex, respectively. Arginine 292 forms a salt bridge with the beta-COOH of the substrate, aspartate. An improvement on the earlier attempt to invert the substrate charge specificity via R292D mutation-induced arginine transaminase activity [Cronin, C. N., & Kirsch, J. F. (1988) Biochemistry 27, 4572-4579] is described. Here Arg292 is replaced with homoglutamate (R292hoGlu). This construct exhibits 6.8 x 10(4)-fold greater activity for the cationic substrate D,L-[Calpha-3H]-alpha-amino-beta-guanidinopropionic acid (D,L-[Calpha-3H]AGPA) than does wild type enzyme. The gain in selectivity for this substrate is at least 4500-fold greater than that achieved in the 1988 experiment, i.e., [(kcat/KM)R292hoGlu/(kcat/KM)WT (D,L-[Calpha-3H]AGPA)] >/= 4500 x [(kcat/KM)R292D/(kcat/KM)WT (L-arginine)]. The value of (kcat/KM)R292D is 0.43 M-1 s-1 with L-Arg while (kcat/KM)R292hoGlu is 29 M-1 s-1 with D,L-[Calpha-3H]AGPA (it is assumed that the D-enantiomer is unreactive). The latter value is the lower limit because of the uncertain value of 3H kinetic isotope effect.

The common substrate structure for the functionally diverse Nudix protein superfamily is nucleotide-diphosphate-X, where X is a large variety of leaving groups. The substrate specificity is known for less than 1% of the 29,400 known members. Most activities result in the release of an inorganic phosphate ion or of a product bearing a terminal phosphate moiety. Reactions have typically been monitored by a modification of the discontinuous Fiske-SubbaRow assay, which is relatively insensitive and slow. We report here the development of a continuous fluorescence assay that enables the rapid and accurate determination of substrate specificities in a 96-well format. We used this novel assay to confirm the reported substrate characterizations of MutT and NudD of Escherichia coli and to characterize DR-1025 of Deinococcus radiodurans and MM-0920 of Methanosarcina mazei. Novel findings enabled by the new assay include the following. First, in addition to the well-characterized hydrolysis of 8-oxo-dGTP at the α-β position, MutT cleaves at the β-γ phosphate bond at a rate of 3% of that recorded for hydrolysis at the α-β position. Second, MutT also catalyzes the hydrolysis of 5-methyl-dCTP. Third, 8-oxo-dGTP was observed to be the best substrate for DR-1025 of the 41 compounds screened. © 2013 Elsevier Inc. All rights reserved.

The periplasmic histidine permease of Salmonella typhimurium is composed of a membrane-bound complex and a soluble histidine-binding protein (the periplasmic receptor), HisJ. Liganded receptor interacts with the membrane-bound complex, inducing ATP hydrolysis and substrate translocation. Preliminary evidence had shown a lack of direct correlation between the affinity of HisJ for a ligand and translocation efficiency, suggesting that the precise form of the receptor is important in determining its interaction with the membrane-bound complex. We have investigated the nature of the conformations assumed by HisJ upon binding a variety of ligands by tryptophan fluorescence enhancement, reaction with a closed form-specific monoclonal antibody, and changes in UV absorption spectra. It is demonstrated that although HisJ binds all the ligands and undergoes a conformational change, it assumes measurably different conformations. We also show that the interaction between HisJ and the membrane-bound complex depends on the nature of the ligand. Transport specificity appears to be defined, at least in part, by the conformation of the bound receptor, manifested either by the effect of a given ligand on the closed structure per se, or by the effect of ligand association on the equilibrium constant relating the open and the closed liganded forms.

It was previously suggested that the conserved Cys-191 of aspartate aminotransferases (AATases) is conserved, not because it is essential, but because it is frozen in the sequence, with no neutral corridor to traverse to the similar phenotype of Ala-191 (Gloss et al., Biochemistry 31:32–39, 1992). This hypothesis has now been tested by additional mutations. All possible one-base mutations from Cys were made at position 191. All of these variants display kinetic parameters (kcat and kcat/KM values) that differ from the wild-type enzyme by 30% or more. The non-conserved cysteines that are predominantly Ala in other AATase sequences (Cys-82, Cys-192, and Cys-401) were mutated to Ser to test the corollary that a neutral Cys → Ala corridor does exist for these positions. These Cys → Ser mutations yielded enzymes with wild-type-like kinetic parameters. The pKa values of the internal aldimines of the mutants, Cys-191 → Ser, Phe, Tyr, and Trp are higher than that of wild type by 0.6–0.8 pH units. The stabilities to urea denaturation of the Cys-191 mutants are similar to that of wild type, while those of the non-conserved cysteines show greater variation. Examination of the three-dimensional environment of the five cysteines showed that the van der Waals contacts of Cys-191 are more conserved than are those of the non-conserved cysteines. These data provide further support for the above hypothesis.

The X-ray structure of a chicken egg white lysozyme (ChEWL) complex with a peptidoglycan-derived inhibitor suggests that interactions of Asn46 and Asp52 with the D-subsite N-acetylmuramic acid residue help to distort that pyranose ring into the reactive half-chair conformation and that a hydrogen bond is formed between Asn46 and Asp52 [Strynadka, N. C. J., & James, M. N. G. (1991) J. Mol. Biol. 220, 401-424]. These hypotheses were investigated through the D52A, N46A, and D52A/N46A mutants of ChEWL. The Michaelis constants of the D52A and D52A/N46A ChEWL complexes with Micrococcus luteus cells are 3- and 4-fold higher, respectively, than the wild-type KM; the corresponding kcat values are 25- and 50-fold lower, respectively, than the wild-type kcat. These results support the proposal of Strynadka and James. The velocities of reactions catalyzed by the N46A and D52A mutants are approximately equal to each other for all classes of substrate, suggesting that the respective roles of Asn46 and Asp52 in transition state stabilization do not vary. The mutation of either Asn46 or Asp52 to Ala apparently disrupts the interactions of the other (nonmutated) residue with the substrate, supporting the crystallographic evidence of a hydrogen-bond interaction between the two residues. The mutations do not change the values of the dissociation constants of complexes with (carboxymethyl)chitin complexes, suggesting that ground state complexes of ChEWL with chitin-derived substrates differ in conformation from complexes with bacterial peptidoglycans.

A series of 24 mutants was made in the buried core of chicken lysozyme at positions 40, 55, and 91. The midpoint temperature of thermal denaturation transition (Tm) values of these core constructs range from 60.9 to 77.3 degrees C, extending an earlier, more limited investigation on thermostability. The Tm values of variants containing conservative replacements for the wild type (WT) (Thr 40-Ile 55-Ser 91) triplet are linearly correlated with hydrophobicity (r = 0.81) and, to a lesser degree, with combined side-chain volume (r = 0.75). The X-ray structures of the S91A (1.9 A) and I55L/S91T/D101S (1.7 A) mutants are presented. The former amino acid change is found in duck and mammalian lysozymes, and the latter contains the most thermostable core triplet. A network of four conserved, buried water molecules is associated with the core. It is postulated that these water molecules significantly influence the mutational tolerance at the individual triplet positions. The pH dependence of Tm for the S91D mutant was compared with that of WT enzyme. The pKa of S91D is 1.2 units higher in the native than in the denatured state, corresponding to delta delta G298 = 1.7 kcal/mol. This is a low value for charge burial and likely reflects the moderating influence of the buried water molecules or a conformational change. Thermal and chemical denaturation and far UV CD spectroscopy were used to characterize the in vitro properties of I55T. This variant, which buries a hydroxyl group, has similar properties to those of the human amyloidogenic variant I56T.

Although several high-resolution X-ray crystallographic structures have been determined for Escherichia coli aspartate aminotransferase (eAATase), efforts to crystallize E. coli tyrosine aminotransferase (eTATase) have been unsuccessful. Sequence alignment analyses of eTATase and eAATase show 43% sequence identity and 72% sequence similarity, allowing for conservative substitutions. The high similarity of the two sequences indicates that both enzymes must have similar secondary and tertiary structures. Six active site residues of eAATase were targeted by homology modeling as being important for aromatic amino acid reactivity with eTATase. Two of these positions (Thr 109 and Asn 297) are invariant in all known aspartate aminotransferase enzymes, but differ in eTATase (Ser 109 and Ser 297). The other four positions (Val 39, Lys 41, Thr 47, and Asn 69) line the active site pocket of eAATase and are replaced by amino acids with more hydrophobic side chains in eTATase (Leu 39, Tyr 41, Ile 47, and Leu 69). These six positions in eAATase were mutated by site-directed mutagenesis to the corresponding amino acids found in eTATase in an attempt to redesign the substrate specificity of eAATase to that of eTATase. Five combinations of the individual mutations were obtained from mutagenesis reactions. The redesigned eAATase mutant containing all six mutations (Hex) displays second-order rate constants for the transamination of aspartate and phenylalanine that are within an order of magnitude of those observed for eTATase. Thus, the reactivity of eAATase with phenylalanine was increased by over three orders of magnitude without sacrificing the high transamination activity with aspartate observed for both enzymes.(ABSTRACT TRUNCATED AT 250 WORDS)

Mutation of six residues of Escherichia coli aspartate aminotransferase results in substantial acquisition of the transamination properties of tyrosine aminotransferase without loss of aspartate transaminase activity. X-ray crystallographic analysis of key inhibitor complexes of the hexamutant reveals the structural basis for this substrate selectivity. It appears that tyrosine aminotransferase achieves nearly equal affinities for a wide range of amino acids by an unusual conformational switch. An active-site arginine residue either shifts its position to electrostatically interact with charged substrates or moves aside to allow access of aromatic ligands.

The pH dependence of Escherichia coli aspartate aminotransferase (AATase) has been investigated by the use of site-directed mutants and alternative substrates. Inhibition of the enzyme by CHES and variations in ionic strength are proposed to explain some of the qualitative differences in the published pH dependence of pig cytosolic AATase kinetics [Velick, S. F., & Vavra, J. (1962) J. Biol. Chem. 237, 2109-2122; Kiick, D.M., & Cook, P.F. (1983) Biochemistry 22, 375-382]. The pKa values of the basic limbs in the kcat/KM profiles for the amino acids, L-Asp and L-cysteinesulfinate (L-CS), are identical, within error, to those of free substrates, (L-Asp, pKa = 9.6; L-CS, pKa = 9.0). This pKa therefore is assigned to the alpha-amino group of the substrate. Replacement of the active site base, Lys-258, with the weaker base, gamma-thia-Lys, does not alter the intrinsic pKa for the profiles of the Ki values for the maleate-E.PMP complexes or the kcat/K alpha-KGM values. The mutation Y225F results in an alkaline shift of the pKa in the kcat/K alph-KGM profile. This pKa is assigned to the C4' amino group of PMP. E. coli AATase, unlike pig cytosolic AATase, shows a pH dependence on kcat between pH 5 and 10 that arises from a change in the rate-determining step at pH extremes. C alpha proton abstraction is partially rate-determining at neutral pH values, but not at pH extremes.(ABSTRACT TRUNCATED AT 250 WORDS)

The pyridoxal phosphate-dependent enzyme 1-aminocyclopropane-1-carboxylate synthase (ACC synthase; S-adenosyl-L-methionine methylthioadenosine-lyase, EC 4.4.1.14) catalyzes the conversion of S-adenosylmethionine (AdoMet) to ACC and 5'-methylthioadenosine, the committed step in ethylene biosynthesis in plants. Apple ACC synthase was overexpressed in Escherichia coli (3 mg/liter) and purified to near homogeneity. A continuous assay was developed by coupling the ACC synthase reaction to the deamination of 5'-methylthioadenosine by adenosine deaminase (adenosine aminohydrolase, EC 3.5.4.4) from Aspergillus oryzae. The enzyme is dimeric, with kcat = 9s-1 per monomer and Km = 12 microM for AdoMet. The pyridoxal phosphate-binding site of ACC synthase appears to be highly homologous to that of aspartate aminotransferase, suggesting similar roles for corresponding residues. Site-directed mutagenesis of Lys-273, Arg-407, and Tyr-233 (corresponding to residues 258, 386, and 225 in aspartate aminotransferase) and kinetic analyses of the mutants confirms their importance in the ACC synthase mechanism. The Lys-273 to Ala mutant has no detectable activity, supporting the identification of this residue as the base catalyzing C alpha proton abstraction. Mutation of Arg-407 to Lys results in a precipitous drop in kcat/Km and an increase in Km for AdoMet of at least 20-fold, in accordance with its proposed role as principal ligand for the substrate alpha-carboxylate group. Replacement of Tyr-233 with Phe causes a 24-fold increase in the Km for AdoMet and no change in kcat, suggesting that this residue plays a role in orienting the pyridoxal phosphate cofactor in the active site.

The strictly conserved active site residue, Asp222, which forms a hydrogen-bonded salt bridge with the pyridine nitrogen atom of the pyridoxal 5' phosphate (PLP) co-factor of aspartate aminotransferase (AATase), was replaced with alanine (D222A) in the Escherichia coli enzyme. The D222A mutant exhibits non-hyberbolic saturation behavior with amino acid substrates which appear as apparent negative cooperativity in steady-state kinetic analyses. Single turnover progress curves for D222A are well described by the sum of two exponentials, contrasting with the monophasic kinetics of the wild-type enzyme. An active/inactive heterodimer containing the D222A mutation retains this biphasic kinetic response, proving that the observed cooperativity is not the result of induced allostery. The anomalous behavior is explained by a hysteretic kinetic model involving two slowly interconverting enzyme forms, only one of which is catalytically competent. The slow functional transition between the two forms has a half-life of approximately 10 mins. Preincubation of the mutant with the dicarboxylic inhibitor maleate shifts the equilibrium population of the enzyme towards the catalytically active form, suggesting that the slow transition is related to the domain closure known to occur upon association of this inhibitor with the wild-type enzyme. The importance of Asp222 in the chemical steps of transamination is confirmed by the approximately 10(5)-fold decrease in catalytic competence in the D222A mutant, and by the large primary C alpha-deuterium kinetic isotope effect (6.7 versus 2.2 for the wild-type). The transamination activity of the D222A mutant is enhanced 4- to 20-fold by reconstitution with the co-factor analog N-methylpyridoxal-5'-phosphate (N-MePLP), and the C alpha-proton abstraction step is less rate determining, as evidenced by the decrease in the primary kinetic isotope effect from 6.7 to 2.3. These results suggest that the conserved interaction between the protonated pyridine nitrogen of PLP and the negatively charged carboxylate of Asp222 is important not only for efficient C alpha-proton abstraction, but also for conformational transitions concomitant with the transamination process.

The complex formed between hen egg white lysozyme (HEL) and the monoclonal antibody HyHEL-10 Fab fragment has an interface composed of van der Waals interactions, hydrogen bonds, and a single ion pair. The antibody overlaps part of the active site cleft. Putative critical residues within the epitope region of HEL, identified from the x-ray crystallographic structure of the complex, were replaced by site-directed mutagenesis to probe their relative importance in determining affinity of the antibody for HEL. Twenty single mutations of HEL at three contact residues (Arg-21HEL, Asp-101HEL, and Gly-102HEL) and at a partially buried residue (Asn-19HEL) in the epitope were made, and the effects on the free energies of dissociation were measured. A correlation between increased amino acid side-chain volume and reduced affinity for HELs with mutations at position 101 was observed. The D101GHEL mutant is bound to HyHEL-10 as tightly as wild-type enzyme, but the delta delta Gdissoc is increased by about 2.2 kcal (9.2 kJ)/mol for the larger residues in this position. HEL variants with lysine or histidine replacements for arginine at position 21 are bound 1.4-2.7 times more tightly than those with neutral or negatively charged amino acids in this position. These exhibit 1/40 the affinity for HyHEL-10 Fab compared with wild type. There is no side-chain volume correlation with delta delta Gdissoc at position 21. Although Gly-102HEL and Asn-19HEL are in the epitope, replacements at these positions have no effect on the affinity of HEL for the antibody.

Motivation: The identification of catalytic residues is a key step in understanding the function of enzymes. While a variety of computational methods have been developed for this task, accuracies have remained fairly low. The best existing method exploits information from sequence and structure to achieve a precision (the fraction of predicted catalytic residues that are catalytic) of 18.5% at a corresponding recall (the fraction of catalytic residues identified) of 57% on a standard benchmark. Here we present a new method, Discern, which provides a significant improvement over the state-of-the-art through the use of statistical techniques to derive a model with a small set of features that are jointly predictive of enzyme active sites. Results: In cross-validation experiments on two benchmark datasets from the Catalytic Site Atlas and CATRES resources containing a total of 437 manually curated enzymes spanning 487 SCOP families, Discern increases catalytic site recall between 12% and 20% over methods that combine information from both sequence and structure, and by ≥50% over methods that make use of sequence conservation signal only. Controlled experiments show that Discern's improvement in catalytic residue prediction is derived from the combination of three ingredients: the use of the INTREPID phylogenomic method to extract conservation information; the use of 3D structure data, including features computed for residues that are proximal in the structure; and a statistical regularization procedure to prevent overfitting. Contact: kimmen@berkeley.edu Supplementary information: Supplementary data are available at Bioinformatics online.

The five cysteines, at positions 82, 191, 192, 270, and 401, of Escherichia coli aspartate aminotransferase (AATase) were, individually and in some combinations, converted to alanine by site-directed mutagenesis (C82A, C191A, C192A, C270A, C401A). Cys-191, which is conserved in all AATase isozymes, was mutated to serine as well (C191S). A quintuple mutant, with all cysteines converted to alanines (Quint), was also constructed. The effects of these single and multiple mutations were examined by steady-state kinetics and urea denaturation. The thermal stabilities of Quint and of the wild-type enzyme (WT) were determined by differential scanning calorimetry. The mutants had kcat values up to 50% greater than that of WT and KMAsp and KM alpha-KG values up to 1.5- and 3.3-fold higher than that of WT. The mutants C82A and C191A exhibit nearly the same CM in urea denaturation experiments as WT, while the other single mutants and Quint are less stable, with CM differences of up to 0.7 M urea. Quint is also less thermostable than WT, with a delta TM of 3.3-4.4 degrees C. Thus the five cysteine replacements yield small, but significant, changes in catalytic and denaturation parameters, but none of the cysteines was found to be essential. The changes manifested in the mutation of the conserved Cys-191 to alanine are no greater than those observed with the four nonconserved cysteines. We consider the evolutionary implications of these findings.

Primary amines functionally replace lysine 258 by catalyzing both the 1,3-prototropic shift and external aldimine hydrolysis reactions with the inactive aspartate aminotransferase mutant K258A. This finding allows classical Bronsted analyses of proton transfer reactions to be applied to enzyme-catalyzed reactions. An earlier study of the reaction of K258A with cysteine sulfinate (Toney, M.D. & Kirsch, J.F., 1989, Science 243, 1485) provided a beta value of 0.4 for the 1,3-prototropic shift. The beta value reported here for the transamination of oxalacetate to aspartate is 0.6. The catalytic efficacy of primary amines is largely determined by basicity and molecular volume. The dependence of the rate constants for the reactions of K258A and K258M on amine molecular volume is nearly identical. This observation argues that the alkyl groups of the added amines do not occupy the position of the lysine 258 side chain in the wild type enzyme. Large primary C alpha and insignificant solvent deuterium kinetic isotope effects with amino acid substrates demonstrate that the amine nitrogen of the exogenous catalysts directly abstracts the labile proton in the rate-determining step.

Lys-258 of aspartate aminotransferase forms a Schiff base with pyridoxal phosphate and is responsible for catalysis of the 1,3-prototropic shift central to the transamination reaction sequence. Substitution of arginine for Lys-258 stabilizes the otherwise elusive quinonoid intermediate, as assessed by the long wavelength absorption bands observed in the reactions of this mutant with several amino acid substrates. The external aldimine intermediate is not detectable during reactions of this mutant with amino acids, although the inhibitor alpha-methylaspartate does slowly and stably form this species. These results suggest that external aldimine formation is one of the rate-determining steps of the reaction. The pyridoxamine-5'-phosphate-like enzyme form (330-nm absorption maximum) is unreactive toward keto acid substrates, and the coenzyme bound to this species is not dissociable from the protein.

The active-site essential catalytic residue of aspartate aminotransferase, Lys 258, has been converted to Cys (K258C) by site-directed mutagenesis. This mutant retains less than 10(-6) of the wild-type activity with L-aspartate. The deleted general base was functionally replaced by selective (with respect to the other five cysteines in wild type) aminoethylation of the introduced Cys 258 with (2-bromoethyl)amine following reversible protection of the nontarget sulfhydryl groups at different stages of unfolding. The chemically elaborated mutant (K258C-EA) is 10(5) times more reactive than is K258C and has a kcat value of approximately 7% of that of wild type (WT). Km and KI values are similar to those for WT. The acidic pKa controlling V/KAsp is shifted from 7.3 (WT) to 6.0 (mutant). V/K values for amino acids are approximately 3% of those found for WT, whereas they are approximately 20% for keto acids. The value of DV increases from 1.6 for WT to 3.4 for the mutant, indicating that C alpha proton abstraction constitutes a more significant kinetic barrier for the latter enzyme. A smaller, but still significant, increase in D(V/KAsp) from 1.9 in WT to 3.0 in the mutant shows that the forward and reverse commitment factors are inverted by the mutation. The acidic limb of the V/KAsp versus pH profile, is lowered by 1.3 pH units, probably reflecting the similar difference in the basicity of the epsilon-NH2 group in gamma-thialysine versus that in lysine.(ABSTRACT TRUNCATED AT 250 WORDS)