OptZyme: Computational Enzyme Redesign Using Transition State Analogues

OptZyme is a new computational procedure for designing improved enzymatic activity (i.e., k_cat or k_cat/K_M) with a novel substrate. The key concept is to use transition state analogue compounds, which are known for many reactions, as proxies for the typically unknown transition state structures. Mutations that minimize the interaction energy of the enzyme with its transition state analogue, rather than with its substrate, are identified that lower the transition state formation energy barrier. Using Escherichia coli β-glucuronidase as a benchmark system, we confirm that K_M correlates (R² = 0.960) with the computed interaction energy between the enzyme and the para-nitrophenyl- β, D-glucuronide substrate, k_cat/K_M correlates (R² = 0.864) with the interaction energy of the transition state analogue, 1,5-glucarolactone, and k_cat correlates (R² = 0.854) with a weighted combination of interaction energies with the substrate and transition state analogue. OptZyme is subsequently used to identify mutants with improved K_M, k_cat, and k_cat/K_M for a new substrate, para-nitrophenyl- β, D-galactoside. Differences between the three libraries reveal structural differences that underpin improving K_M, k_cat, or k_cat/K_M. Mutants predicted to enhance the activity for para-nitrophenyl- β, D-galactoside directly or indirectly create hydrogen bonds with the altered sugar ring conformation or its substituents, namely H162S, L361G, W549R, and N550S.     

Tissue kallikrein processes small proenkephalin peptides

Tissue kallikrein may play a role in processing precursor polypeptide hormones. We investigated whether hydrolysis of natural enkephalin precursors, peptide F and bovine adrenal medulla docosapeptide (BAM-22P), by hog pancreatic kallikrein is consistent with this concept. Incubation of peptide F with this tissue kallikrein resulted in the release of Met⁵-enkephalin and Met⁵-Lys⁶-enkephalin. Met⁵-Lys⁶-enkephalin was the main peptide released, indicating that the major cleavage site was between two lysine residues. At 37°C and pH 8.5, the K_M values for formation of Met⁵-enkephalin and Met⁵-Lys⁶-enkephalin were 129 and 191 μM, respectively. Corresponding k_cat values were 0.001 and 0.03 s⁻¹ and k_cat/K_M ratios were 8 and 1.6·10² M⁻¹ · s⁻¹, respectively. Cleavage of peptide F at acidic pH (5.5) was negligible. When BAM-22P was used as a substrate, Met⁵-Arg⁶-enkephalin was released, thus indicating cleavage between two arginine residues. At pH 8.5, K_M was 64 μM, k_cat was 4.5 s⁻¹, and the k_cat/K_M ratio was 7 · 10⁴ M⁻¹ · s⁻¹. At 5.5, the pH of the secretory granules, K_M, k_cat and k_cat/K_M were 184 μM, 1.9 s⁻¹ and 10⁴ M⁻¹ · s⁻¹, respectively. It is unlikely that peptide F could be a substrate for kallikrein in vivo; however, tissue kallikrein could aid in processing proenkephalin precursors such as BAM-22P by cleaving Arg-Arg peptide bonds.     

Evaluation of kinetic data: What the numbers tell us about PRMTs

Protein arginine N-methyltransferase (PRMT) kinetic parameters have been catalogued over the past fifteen years for eight of the nine mammalian enzyme family members. Like the majority of methyltransferases, these enzymes employ the highly ubiquitous cofactor S-adenosyl-l-methionine as a co-substrate to methylate arginine residues in peptidic substrates with an approximately 4-μM median K_M. The median values for PRMT turnover number (k_cat) and catalytic efficiency (k_cat/K_M) are 0.0051 s⁻¹ and 708 M⁻¹ s⁻¹, respectively. When comparing PRMT metrics to entries found in the BRENDA database, we find that while PRMTs exhibit high substrate affinity relative to other enzyme-substrate pairs, PRMTs display largely lower k_cat and k_cat/K_M values. We observe that kinetic parameters for PRMTs and arginine demethylase activity from dual-functioning lysine demethylases are statistically similar, paralleling what the broader enzyme families in which they belong reveal, and adding to the evidence in support of arginine methylation reversibility.     

Cloning,expression and some properties of α-carbonic anhydrase from Helicobacter pylori

The α-carbonic anhydrase gene from Helicobacter pylori strain 26695 has been cloned and sequenced. The full-length protein appears to be toxic to Escherichia coli, so we prepared a modified form of the gene lacking a part that presumably encodes a cleavable signal peptide. This truncated gene could be expressed in E. coli yielding an active enzyme comprising 229 amino acid residues. The amino acid sequence shows 36% identity with that of the enzyme from Neisseria gonorrhoeae and 28% with that of human carbonic anhydrase II. The H. pylori enzyme was purified by sulfonamide affinity chromatography and its circular dichroism spectrum and denaturation profile in guanidine hydrochloride have been measured. Kinetic parameters for CO₂ hydration catalyzed by the H. pylori enzyme at pH 8.9 and 25°C are k_cat=2.4×10⁵ s⁻¹, K_M=17 mM and k_cat/K_M=1.4×10⁷ M⁻¹ s⁻¹. The pH dependence of k_cat/K_M fits with a simple titration curve with pK_a=7.5. Thiocyanate yields an uncompetitive inhibition pattern at pH 9 indicating that the maximal rate of CO₂ hydration is limited by proton transfer between a zinc-bound water molecule and the reaction medium in analogy to other forms of the enzyme. The 4-nitrophenyl acetate hydrolase activity of the H. pylori enzyme is quite low with an apparent catalytic second-order rate constant, k_enz, of 24 M⁻¹ s⁻¹ at pH 8.8 and 25°C. However, with 2-nitrophenyl acetate as substrate a k_enz value of 665 M⁻¹ s⁻¹ was obtained under similar conditions.     

Rational design of a SHP-2 targeted,fluorogenic peptide substrate

Protein tyrosine phosphatase (PTP) targeted, peptide based chemical probes are valuable tools for studying this important family of enzymes, despite the inherent difficulty of developing peptides targeted towards an individual PTP. Here, we have taken a rational approach to designing a SHP-2 targeted, fluorogenic peptide substrate based on information about the potential biological substrates of SHP-2. The fluorogenic, phosphotyrosine mimetic phosphocoumaryl aminopropionic acid (pCAP) provides a facile readout for monitoring PTP activity. By optimizing the amino acids surrounding the pCAP residue, we obtained a substrate with the sequence Ac-DDPI-pCAP-DVLD-NH₂ and optimized kinetic parameters (k_cat = 0.059 ± 0.008 s⁻¹, K_m = 220 ± 50 µM, k_cat/K_m of 270 M⁻¹s⁻¹). In comparison, the phosphorylated coumarin moiety alone is an exceedingly poor substrate for SHP-2, with a k_cat value of 0.0038 ± 0.0003 s⁻¹, a K_m value of 1100 ± 100 µM and a k_cat/K_m of 3 M⁻¹s⁻¹. Furthermore, this optimized peptide has selectivity for SHP-2 over HePTP, MEG1 and PTPµ. The data presented here demonstrate that PTP-targeted peptide substrates can be obtained by optimizing the sequence of a pCAP containing peptide.     

The base exchange reaction of NAD+ glycohydrolase: identification of novel heterocyclic alternative substrates

ADP-ribosyl cyclase and NAD⁺ glycohydrolase (CD38, E.C.3.2.2.5) efficiently catalyze the exchange of the nicotinamidyl moiety of NAD⁺, nicotinamide adenine dinucleotide phosphate (NADP⁺) or nicotinamide mononucleotide (NMN⁺) with an alternative base. 4′-Pyridinyl drugs (amrinone, milrinone, dismerinone and pinacidil) were efficient alternative substrates (k_cat/K_M = 0.9-10 μM⁻¹ s⁻¹) in the exchange reaction with ADP-ribosyl cyclase. When CD38 was used as a catalyst the k_cat/K_M values for the exchange reaction were reduced two or more orders of magnitude (0.015-0.15 μM⁻¹ s⁻¹). The products of this reaction were novel dinucleotides. The values of the equilibrium constants for dinucleotide formation were determined for several drugs. These enzymes also efficiently catalyze the formation of novel mononucleotides in an exchange reaction with NMN⁺, k_cat/K_M = 0.05-0.4 μM⁻¹ s⁻¹. The k_cat/K_M values for the exchange reaction with NMN⁺ were generally similar (0.04-0.12 μM⁻¹ s⁻¹) with CD38 and ADP-ribosyl cyclase as catalysts. Several novel heterocyclic alternative substrates were identified as 2-isoquinolines, 1,6-naphthyridines and tricyclic bases. The k_cat/K_M values for the exchange reaction with these substrates varied over five orders of magnitude and approached the limit of diffusion with 1,6-naphthyridines. The exchange reaction could be used to synthesize novel mononucleotides or to identify novel reversible inhibitors of CD38.     

Kinetic studies of Gly28:Ser mutant form of Bacillus pumilus lipase: Changes in kcat and thermal dependence

Lipases are useful catalysts for a wide variety of industrial purposes. Herein we report the stability and thermal dependence of the activity of wild-type Bacillus pumilus lipase (BplA) and four site-directed mutants designed to improve its thermal stability. The Gly28:Ser mutation produces a dramatic four-fold increase in its k_cat and a remarkable increase in its stability. While the increase in k_cat is temperature-independent, the increase in stability shows that the resultant interactions of this mutation have a strong enthalpic component. Thermal dependence of stability, k_cat, K_M and k_cat/K_M were analysed to gain insight on the structural effects of mutations on BplA. Our results are consistent with a gain in enzyme mobility for those mutants displaying enhanced catalytic properties; the analysis of thermal dependence of kinetic parameters indicates that the mutations did not change either the catalytic mechanism or the rate-limiting step of catalysis.     

A lysine conserved in the monoamine oxidase family is involved in oxidation of the reduced flavin in mouse polyamine oxidase

Lysine 315 of mouse polyamine amine oxidase corresponds to a lysine residue that is conserved in the flavoprotein amine oxidases of the monoamine oxidase structural family. In several structures, this lysine residue forms a hydrogen bond to a water molecule that is hydrogen-bonded to the flavin N(5). Mutation of Lys315 in polyamine oxidase to methionine was previously shown to have no effect on the kinetics of the reductive half-reaction of the enzyme (M. Henderson Pozzi, V. Gawandi, P.F. Fitzpatrick, Biochemistry 48 (2009) 1508-1516). In contrast, the mutation does affect steps in the oxidative half-reaction. The k_cat value is unaffected by the mutation; this kinetic parameter likely reflects product release. At pH 10, the k_cat/K_m value for oxygen is 25-fold lower in the mutant enzyme. The k_cat/K_O2 value is pH-dependent for the wild-type enzyme, decreasing below a pK_a of 7.0, while this kinetic parameter for the mutant enzyme is pH-independent. This is consistent with the neutral form of Lys315 being required for more rapid flavin oxidation. The solvent isotope effect on the k_cat/K_O2 value increases from 1.4 in the wild-type enzyme to 1.9 in the mutant protein, and the solvent inventory changes from linear to bowed. The effects of the mutation can be explained by the lysine orienting the bridging water so that it can accept the proton from the flavin N(5) during flavin oxidation. In the mutant enzyme the lysine amine would be replaced by a water chain.     

Solvent isotope effects on α-glucosidase

The solvent kinetic isotope effects (SKIE) on the yeast α-glucosidase-catalyzed hydrolysis of p-nitrophenyl and methyl-d-glucopyranoside were measured at 25 °C. With p-nitrophenyl-d-glucopyranoside (pNPG), the dependence of k_cat/K_m on pH (pD) revealed an unusually large (for glycohydrolases) solvent isotope effect on the pL-independent second-order rate constant, ^DOD(k_cat/K_m), of 1.9 (±0.3). The two pK_as characterizing the pH profile were increased in D₂O. The shift in pK_a2 of 0.6 units is typical of acids of comparable acidity (pK_a=6.5), but the increase in pK_a1 (=5.7) of 0.1 unit in going from H₂O to D₂O is unusually small. The initial velocities show substrate inhibition (K_is/K_m～200) with a small solvent isotope effect on the inhibition constant [^DODK_is=1.1 (±0.2)]. The solvent equilibrium isotope effects on the K_is for the competitive inhibitors d-glucose and α-methyl d-glucoside are somewhat higher [^DODK_i=1.5 (±0.1)]. Methyl glucoside is much less reactive than pNPG, with k_cat 230 times lower and k_cat/K_m 5×10⁴ times lower. The solvent isotope effect on k_cat for this substrate [=1.11 (±0. 02)] is lower than that for pNPG [=1.67 (±0.07)], consistent with more extensive proton transfer in the transition state for the deglucosylation step than for the glucosylation step.     

Steady-state kinetic mechanism of LodA,a novel cysteine tryptophylquinone-dependent oxidase

LodA is a novel lysine-ε-oxidase which possesses a cysteine tryptophylquinone cofactor. It is the first tryptophylquinone enzyme known to function as an oxidase. A steady-state kinetic analysis shows that LodA obeys a ping-pong kinetic mechanism with values of k_cat of 0.22 ± 0.04 s⁻¹, K_lysine of 3.2 ± 0.5 μM and K_O2 of 37.2 ± 6.1 μM. The k_cat exhibited a pH optimum at 7.5 while k_cat/K_lysine peaked at 7.0 and remained constant to pH 8.5. Alternative electron acceptors could not effectively substitute for O₂ in the reaction. A mechanism for the reductive half reaction of LodA is proposed that is consistent with the ping-pong kinetics.     

Investigation of cobalt(II) substituted carboxypeptidase a interacting with azide and cyanate ions

Cobalt(II)-substituted carboxypeptidase A has been found to reversibly bind N3^? and NCO^?, but not NCS^?, in the pH range 5–10, thus including the pH range of activity of the enzyme. The pH dependence of the anion binding constant is affected by two ionizations, which are assigned as those regulating k_cat and K_M. The electronic and ₁H NMR spectra are consistent with a substantially pseudotetrahedral geometry of the anion derivatives.     

Substrate specificity and kinetic framework of a DNAzyme with an expanded chemical repertoire: a putative RNaseA mimic that catalyzes RNA hydrolysis independent of a divalent metal cation

This work addresses the binding, cleavage and dissociation rates for the substrate and products of a synthetic RNaseA mimic that was combinatorially selected using chemically modified nucleoside triphosphates. This trans-cleaving DNAzyme, 9₂₅-11t, catalyzes sequence-specific ribophosphodiester hydrolysis in the total absence of a divalent metal cation, and in low ionic strength at pH 7.5 and in the presence of EDTA. It is the first such sequence capable of multiple turnover. 9₂₅-11t consists of 31 bases, 18 of which form a catalytic domain containing 4 imidazole and 6 allylamino modified nucleotides. This sequence cleaves the 15 nt long substrate, S1, at one embedded ribocytosine at the eighth position to give a 5′-product terminating in a 2′,3′-phosphodiester and a 3′-product terminating in a 5′-OH. Under single turnover conditions at 24°C, 9₂₅-11t displays a maximum first-order rate constant, k_cat, of 0.037 min⁻¹ and a catalytic efficiency, k_cat/K_m, of 5.3 × 10⁵ M⁻¹ min⁻¹. The measured value of k_cat under catalyst excess conditions agrees with the value of k_cat observed for steady-state multiple turnover, implying that slow product release is not rate limiting with respect to multiple turnover. The substrate specificity of 9₂₅-11t was gauged in terms of k_cat values for substrate sequence variants. Base substitutions on the scissile ribose and at the two bases immediately downstream decrease k_cat values by a factor of 4 to 250, indicating that 9₂₅-11t displays significant sequence specificity despite the lack of an apparent Watson–Crick base-pairing scheme for recognition.     

Cloning, Expression, and Enzymatic Characterization of Isocitrate Dehydrogenase from Helicobacter pylori

Isocitrate dehydrogenase (IDH) catalyzes the oxidative decarboxylation of isocitrate to α-ketoglutarate with NAD(P) as a cofactor in the tricarboxylic acid cycle. As a housekeeping protein in Helicobacter pylori, IDH was considered as a possible candidate for serological diagnostics and detection. Here, we identified a new icd gene encoding IDH from H. pylori strain SS1. The recombinant H. pylori isocitrate dehydrogenase (HpIDH) was cloned, expressed, and purified in E. coli system. The enzymatic characterization of HpIDH demonstrates its activity with k _cat of 87 s^?1, K _m of 124 μM and k _cat/K _m of 7 × 10⁵ M^?1s^?1 toward isocitrate, k _cat of 80 s^?1, K _m of 176 μM and k _cat/K _m of 4.5 × 10⁵ M^?1s^?1 toward NADP. The optimum pH of the enzyme activity is around 9.0, and the optimum temperature is around 50 °C. This current work is expected to help better understand the features of HpIDH and provide useful information for H. pylori serological diagnostics and detection.     

Novel internally quenched substrate of the trypsin-like subunit of 20S eukaryotic proteasome

This article describes the synthesis, using combinatorial chemistry, of internally quenched substrates of the trypsin-like subunit of human 20S proteasome. Such substrates were optimized in both the nonprime and prime regions of the peptide chain. Two were selected as the most susceptible for proteasomal proteolysis with excellent kinetic parameters: (i) ABZ-Val-Val-Ser-Arg-Ser-Leu-Gly-Tyr(3-NO₂)-NH₂ (k_cat/K_M = 934,000 M⁻¹ s⁻¹) and (ii) ABZ-Val-Val-Ser-GNF-Ala-Met-Gly-Tyr(3-NO₂)-NH₂ (k_cat/K_M = 1,980,000 M⁻¹ s⁻¹). Both compounds were efficiently hydrolyzed by the 20S proteasome at picomolar concentrations, demonstrating significant selectivity over other proteasome entities.     

Methanothermobacter thermautotrophicus tRNA Gln confines the amidotransferase GatCAB to asparaginyl-tRNA Asn formation

Many prokaryotes form the amide aminoacyl-tRNAs glutaminyl-tRNA and asparaginyl-tRNA by tRNA-dependent amidation of the mischarged tRNA species, glutamyl-tRNA^Gln or aspartyl-tRNA^Asn. Archaea employ two such amidotransferases, GatCAB and GatDE, while bacteria possess only one, GatCAB. The Methanothermobacter thermautotrophicus GatDE is slightly more efficient using Asn as an amide donor than Gln (k_cat/K_M of 5.4 s⁻¹/mM and 1.2 s⁻¹/mM, respectively). Unlike the bacterial GatCAB enzymes studied to date, the M. thermautotrophicus GatCAB uses Asn almost as well as Gln as an amide donor (k_cat/K_M of 5.7 s⁻¹/mM and 16.7 s⁻¹/mM, respectively). In contrast to the initial characterization of the M. thermautotrophicus GatCAB as being able to form Asn-tRNA^Asn and Gln-tRNA^Gln, our data demonstrate that while the enzyme is able to transamidate Asp-tRNA^Asn (k_cat/K_M of 125 s⁻¹/mM) it is unable to transamidate M. thermautotrophicus Glu-tRNA^Gln. However, M. thermautotrophicus GatCAB is capable of transamidating Glu-tRNA^Gln from H. pylori or B. subtilis, and M. thermautotrophicus Glu-tRNA^Asn. Thus, M. thermautotrophicus encodes two amidotransferases, each with its own activity, GatDE for Gln-tRNA and GatCAB for Asn-tRNA synthesis.     

Origin of the Proton-transfer Step in the Cofactor-free (1H)-3-Hydroxy-4-oxoquinaldine 2,4-Dioxygenase: EFFECT OF THE BASICITY OF AN ACTIVE SITE HIS RESIDUE*

Dioxygenases catalyze a diverse range of chemical reactions that involve the incorporation of oxygen into a substrate and typically use a transition metal or organic cofactor for reaction. Bacterial (1H)-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (HOD) belongs to a class of oxygenases able to catalyze this energetically unfavorable reaction without any cofactor. In the quinaldine metabolic pathway, HOD breaks down its natural N-heteroaromatic substrate using a mechanism that is still incompletely understood. Experimental and computational approaches were combined to study the initial step of the catalytic cycle. We have investigated the role of the active site His-251/Asp-126 dyad, proposed to be involved in substrate hydroxyl group deprotonation, a critical requirement for subsequent oxygen reaction. The pH profiles obtained under steady-state conditions for the H251A and D126A variants show a strong pH effect on their k_cat and k_cat/K_m constants, with a decrease in k_cat/K_m of 5500- and 9-fold at pH 10.5, respectively. Substrate deprotonation studies under transient-state conditions show that this step is not rate-limiting and yield a pK_a value of ∼7.2 for WT HOD. A large solvent isotope effect was found, and the pK_a value was shifted to ∼8.3 in D₂O. Crystallographic and computational studies reveal that the mutations have a minor effect on substrate positioning. Computational work shows that both His-251 and Asp-126 are essential for the proton transfer driving force of the initial reaction. This multidisciplinary study offers unambiguous support to the view that substrate deprotonation, driven by the His/Asp dyad, is an essential requirement for its activation.     

Purification and Characterization of Aspartic Protease Derived from Sf9 Insect Cells

An aspartic protease that is significantly produced by baculovirus-infected Spodoptera frugiperda Sf9 insect cells was purified to homogeneity from a growth medium. To monitor aspartic protease activity, an internally quenched fluoresce (IQF) substrate specific to cathepsin D was used. The purified aspartic protease showed a single protein band on SDS–PAGE with an apparent molecular mass of 40 kDa. The N-terminal amino acid sequence of the enzyme had a high homology to a Bombyx mori aspartic protease. The enzyme showed greatest affinity for the IQF substrate at pH 3.0 with a K _m of 0.85 μM. The k _cat and k _cat?K _m values were 13 s^?1 and 15 s^?1 μM^?1 respectively. Pepstatin A proved to be a potent competitive inhibitor with inhibitor constant, K _i, of 25 pM.     

Molecular Engineering of Fungal GH5 and GH26 Beta-(1,4)-Mannanases toward Improvement of Enzyme Activity

Microbial mannanases are biotechnologically important enzymes since they target the hydrolysis of hemicellulosic polysaccharides of softwood biomass into simple molecules like manno-oligosaccharides and mannose. In this study, we have implemented a strategy of molecular engineering in the yeast Yarrowia lipolytica to improve the specific activity of two fungal endo-mannanases, PaMan5A and PaMan26A, which belong to the glycoside hydrolase (GH) families GH5 and GH26, respectively. Following random mutagenesis and two steps of high-throughput enzymatic screening, we identified several PaMan5A and PaMan26A mutants that displayed improved kinetic constants for the hydrolysis of galactomannan. Examination of the three-dimensional structures of PaMan5A and PaMan26A revealed which of the mutated residues are potentially important for enzyme function. Among them, the PaMan5A-G311S single mutant, which displayed an impressive 8.2-fold increase in k_cat/K_M due to a significant decrease of K_M, is located within the core of the enzyme. The PaMan5A-K139R/Y223H double mutant revealed modification of hydrolysis products probably in relation to an amino-acid substitution located nearby one of the positive subsites. The PaMan26A-P140L/D416G double mutant yielded a 30% increase in k_cat/K_M compared to the parental enzyme. It displayed a mutation in the linker region (P140L) that may confer more flexibility to the linker and another mutation (D416G) located at the entrance of the catalytic cleft that may promote the entrance of the substrate into the active site. Taken together, these results show that the directed evolution strategy implemented in this study was very pertinent since a straightforward round of random mutagenesis yielded significantly improved variants, in terms of catalytic efiiciency (k_cat/K_M).     

Exploring the active site of tripeptidyl-peptidase II through studies of pH dependence of reaction kinetics

Tripeptidyl-peptidase II (TPP II) is a subtilisin-like serine protease which forms a large enzyme complex (> 4 MDa). It is considered a potential drug target due to its involvement in specific physiological processes. However, information is scarce concerning the kinetic characteristics of TPP II and its active site features, which are important for design of efficient inhibitors. To amend this, we probed the active site by determining the pH dependence of TPP II catalysis. Access to pure enzyme is a prerequisite for kinetic investigations and herein we introduce the first efficient purification system for heterologously expressed mammalian TPP II. The pH dependence of kinetic parameters for hydrolysis of two different chromogenic substrates, Ala-Ala-Phe-pNA and Ala-Ala-Ala-pNA, was determined for murine, human and Drosophila melanogaster TPP II as well as mutant variants thereof. The investigation demonstrated that TPP II, in contrast to subtilisin, has a bell-shaped pH dependence of k_cat^app/K_M probably due to deprotonation of the N-terminal amino group of the substrate at higher pH. Since both the K_M and k_cat^app are lower for cleavage of AAA-pNA than for AAF-pNA we propose that the former can bind non-productively to the active site of the enzyme, a phenomenon previously observed with some substrates for subtilisin. Two mutant variants, H267A and D387G, showed bell-shaped pH-dependence of k_cat^app, possibly due to an impaired protonation of the leaving group. This work reveals previously unknown differences between TPP II orthologues and subtilisin as well as features that might be conserved within the entire family of subtilisin-like serine peptidases.     

An algebraic model to determine substrate kinetic parameters by global nonlinear fit of progress curves

We propose that the time course of an enzyme reaction following the Michaelis-Menten reaction mechanism can be conveniently described by a newly derived algebraic equation, which includes the Lambert Omega function. Following Northrop's ideas [Anal. Biochem.321, 457–461, 1983], the integrated rate equation contains the Michaelis constant (K_M) and the specificity number (k_S≡k_cat/K_M$k_{\text{S}}\equiv k_{\text{cat}}/K_{\text{M}}$) as adjustable parameters, but not the turnover number k_cat. A modification of the usual global-fit approach involves a combinatorial treatment of nominal substrate concentrations being treated as fixed or alternately optimized model parameters. The newly proposed method is compared with the standard approach based on the “initial linear region” of the reaction progress curves, followed by nonlinear fit of initial rates to the hyperbolic Michaelis-Menten equation. A representative set of three chelation-enhanced fluorescence EGFR kinase substrates is used for experimental illustration. In one case, both data analysis methods (linear and nonlinear) produced identical results. However, in another test case, the standard method incorrectly reported a finite (50–70 μM) K_M value, whereas the more rigorous global nonlinear fit shows that the K_M is immeasurably high.