Analysis of progress curves by simulations generated by numerical integration

Analysis of continuous enzyme kinetic data using ICEKAT

ICEKAT (Interactive Continuous Enzyme Analysis Tool) is an interactive web-based program for calculating initial rates and kinetic parameters (e.g., Vmax, kcat, KM, EC50, IC50) from continuous enzyme kinetic assay data that satisfy Michaelis-Menten and steady-state kinetic assumptions. ICEKAT is valuable in educational and research settings to consistently and accurately calculate initial rates and kinetic parameters, increasing assay veracity and reproducibility. Provided freely online to the scientific community, ICEKAT has been cited in at least 26 publications, and the initial journal article has been accessed nearly 9000 times since its debut in 2020 (Olp et al., 2020). Here, we provide in-depth instructions for software use, offer vital considerations for data analysis, and highlight updated software features for new and existing users. Through ICEKAT, we aim for the analysis of data from continuous enzyme kinetic studies worldwide to become more rapid, reliable, and repeatable. ICEKAT remains free of charge and available to all scientists at https://icekat.herokuapp.com/icekat; the source code for local use is found at https://github.com/SmithLabMCW/icekat.

pH dependencies of glycolytic enzymes of yeast under in vivo-like assay conditions

Under carbon source transitions, the intracellular pH of Saccharomyces cerevisiae is subject to change. Dynamics in pH modulate the activity of the glycolytic enzymes, resulting in a change in glycolytic flux and ultimately cell growth. To understand how pH affects the global behavior of glycolysis and ethanol fermentation, we measured the activity of the glycolytic and fermentative enzymes in S. cerevisiae under in vivo-like conditions at different pH. We demonstrate that glycolytic enzymes exhibit differential pH dependencies, and optima, in the pH range observed during carbon source transitions. The forward reaction of GAPDH shows the highest decrease in activity, 83%, during a simulated feast/famine regime upon glucose removal (cytosolic pH drop from 7.1 to 6.4). We complement our biochemical characterization of the glycolytic enzymes by fitting the V_max to the progression curves of product formation or decay over time. The fitting analysis shows that the observed changes in enzyme activities require changes in V_max , but changes in K_m cannot be excluded. Our study highlights the relevance of pH as a key player in metabolic regulation and provides a large set of quantitative data that can be explored to improve our understanding of metabolism in dynamic environments.

Kinetic Modeling of Time-Dependent Enzyme Inhibition by Pre-Steady-State Analysis of Progress Curves: The Case Study of the Anti-Alzheimer's Drug Galantamine

The Michaelis–Menten model of enzyme kinetic assumes the free ligand approximation, the steady-state approximation and the rapid equilibrium approximation. Analytical methods to model slow-binding inhibitors by the analysis of initial velocities have been developed but, due to their inherent complexity, they are seldom employed. In order to circumvent the complications that arise from the violation of the rapid equilibrium assumption, inhibition is commonly evaluated by pre-incubating the enzyme and the inhibitors so that, even for slow inhibitors, the binding equilibrium is established before the reaction is started. Here, we show that for long drug-target residence time inhibitors, the conventional analysis of initial velocities by the linear regression of double-reciprocal plots fails to provide a correct description of the inhibition mechanism. As a case study, the inhibition of acetylcholinesterase by galantamine, a drug approved for the symptomatic treatment of Alzheimer’s disease, is reported. For over 50 years, analysis based on the conventional steady-state model has overlooked the time-dependent nature of galantamine inhibition, leading to an erroneous assessment of the drug potency and, hence, to discrepancies between biochemical data and the pharmacological evidence. Re-examination of acetylcholinesterase inhibition by pre-steady state analysis of the reaction progress curves showed that the potency of galantamine has indeed been underestimated by a factor of ~100.

Monitoring protein conformational changes using fluorescent nanoantennas

Understanding the relationship between protein structural dynamics and function is crucial for both basic research and biotechnology. However, methods for studying the fast dynamics of structural changes are limited. Here, we introduce fluorescent nanoantennas as a spectroscopic technique to sense and report protein conformational changes through noncovalent dye-protein interactions. Using experiments and molecular simulations, we detect and characterize five distinct conformational states of intestinal alkaline phosphatase, including the transient enzyme-substrate complex. We also explored the universality of the nanoantenna strategy with another model protein, Protein G and its interaction with antibodies, and demonstrated a rapid screening strategy to identify efficient nanoantennas. These versatile nanoantennas can be used with diverse dyes to monitor small and large conformational changes, suggesting that they could be used to characterize diverse protein movements or in high-throughput screening applications.

interferENZY: A Web-Based Tool for Enzymatic Assay Validation and Standardized Kinetic Analysis

Enzymatic assays are widely employed to characterize important allosteric and enzyme modulation effects. The high sensitivity of these assays can represent a serious problem if the occurrence of experimental errors surreptitiously affects the reliability of enzyme kinetics results. We have addressed this problem and found that hidden assay interferences can be unveiled by the graphical representation of progress curves in modified reaction coordinates. To render this analysis accessible to users across all levels of expertise, we have developed a webserver, interferENZY, that allows (i) an unprecedented tight quality control of experimental data, (ii) the automated identification of small and major assay interferences, and (iii) the estimation of bias-free kinetic parameters. By eliminating the subjectivity factor in kinetic data reporting, interferENZY will contribute to solving the "reproducibility crisis" that currently challenges experimental molecular biology. The interferENZY webserver is freely available (no login required) at https://interferenzy.i3s.up.pt.

ICEKAT: an interactive online tool for calculating initial rates from continuous enzyme kinetic traces

BACKGROUND

Continuous enzyme kinetic assays are often used in high-throughput applications, as they allow rapid acquisition of large amounts of kinetic data and increased confidence compared to discontinuous assays. However, data analysis is often rate-limiting in high-throughput enzyme assays, as manual inspection and selection of a linear range from individual kinetic traces is cumbersome and prone to user error and bias. Currently available software programs are specialized and designed for the analysis of complex enzymatic models. Despite the widespread use of initial rate determination for processing kinetic data sets, no simple and automated program existed for rapid analysis of initial rates from continuous enzyme kinetic traces.

RESULTS

An Interactive Continuous Enzyme Kinetics Analysis Tool (ICEKAT) was developed for semi-automated calculation of initial rates from continuous enzyme kinetic traces with particular application to the evaluation of Michaelis-Menten and EC₅₀/IC₅₀ kinetic parameters, as well as the results of high-throughput screening assays. ICEKAT allows users to interactively fit kinetic traces using convenient browser-based selection tools, ameliorating tedious steps involved in defining ranges to fit in general purpose programs like Microsoft Excel and Graphpad Prism, while still maintaining simplicity in determining initial rates. As a test case, we quickly analyzed over 500 continuous enzyme kinetic traces resulting from experimental data on the response of the protein lysine deacetylase SIRT1 to small-molecule activators.

CONCLUSIONS

ICEKAT allows simultaneous visualization of individual initial rate fits and the resulting Michaelis-Menten or EC₅₀/IC₅₀ kinetic model fits, as well as hits from high-throughput screening assays. In addition to serving as a convenient program for practicing enzymologists, ICEKAT is also a useful teaching aid to visually demonstrate in real-time how incorrect initial rate fits can affect calculated Michaelis-Menten or EC₅₀/IC₅₀ kinetic parameters. For the convenience of the research community, we have made ICEKAT freely available online at https://icekat.herokuapp.com/icekat.

Single particle tracking-based reaction progress kinetic analysis reveals a series of molecular mechanisms of cetuximab-induced EGFR processes in a single living cell

Cellular processes occur through the orchestration of multi-step molecular reactions. Reaction progress kinetic analysis (RPKA) can provide the mechanistic details to elucidate the multi-step molecular reactions. However, current tools have limited ability to simultaneously monitor dynamic variations in multiple complex states at the single molecule level to apply RPKA in living cells. In this research, a single particle tracking-based reaction progress kinetic analysis (sptRPKA) was developed to simultaneously determine the kinetics of multiple states of protein complexes in the membrane of a single living cell. The subpopulation ratios of different states were quantitatively (and statistically) reliably extracted from the diffusion coefficient distribution rapidly acquired by single particle tracking at constant and high density over a long period of time using super-resolution microscopy. Using sptRPKA, a series of molecular mechanisms of epidermal growth factor receptor (EGFR) cellular processing induced by cetuximab were investigated. By comprehensively measuring the rate constants and cooperativity of the molecular reactions involving four EGFR complex states, a previously unknown intermediate state was identified that represents the rate limiting step responsible for the selectivity of cetuximab-induced EGFR endocytosis to cancer cells.

Single-molecule enzymology of steroid transforming enzymes: Transient kinetic studies and what they tell us

Structure-function studies on steroid transforming enzymes often use site-directed mutagenesis to inform mechanisms of catalysis and effects on steroid binding, and data are reported in terms of changes in steady state kinetic parameters kcat, Km and kcat/Km. However, this dissection of function is limited since kcat is governed by the rate-determining step and Km is a complex macroscopic kinetic constant. Often site-directed mutagenesis can lead to a change in the rate-determining step which cannot be revealed by just reporting a decrease in kcat alone. These issues are made more complex when it is considered that many steroid transforming enzymes have more than one substrate and product. We present the case for using transient-kinetics performed with stopped-flow spectrometry to assign rate constants to discrete steps in these multi-substrate reactions and their use to interpret enzyme mechanism and the effects of disease and engineered mutations. We demonstrate that fluorescence kinetic transients can be used to measure ligand binding that may be accompanied by isomerization steps, revealing the existence of new enzyme intermediates. We also demonstrate that single-turnover reactions can provide a klim for the chemical step and Ks for steroid-substrate binding and that when coupled with kinetic isotope effect measurements can provide information on transition state intermediates. We also demonstrate how multiple turnover experiments can provide evidence for either "burst-phase" kinetics, which can reveal a slow product release step, or linear-phase kinetics, in which the chemical step can be rate-determining. With these assignments it becomes more straightforward to analyze the effects of mutations. We use examples from the hydroxysteroid dehydrogenases (AKR1Cs) and human steroid 5β-reductase (AKR1D1) to illustrate the utility of the approach, which are members of the aldo-keto reductase (AKR) superfamily.

Kinesin Motor Enzymology: Chemistry, Structure, and Physics of Nanoscale Molecular Machines

Molecular motors are enzymes that convert chemical potential energy into controlled kinetic energy for mechanical work inside cells. Understanding the biophysics of these motors is essential for appreciating life as well as apprehending diseases that arise from motor malfunction. This review focuses on kinesin motor enzymology with special emphasis on the literature that reports the chemistry, structure and physics of several different kinesin superfamily members.

Analyzing ATP utilization by DEAD-Box RNA helicases using kinetic and equilibrium methods

DEAD-box proteins (DBPs) couple ATP utilization to conformational rearrangement of RNA. In this chapter, we outline a combination of equilibrium and kinetic methods that have been developed and applied to the analysis of ATP utilization and linked RNA remodeling by DBPs, specifically Escherichia coli DbpA and Saccharomyces cerevisiae Mss116. Several important considerations are covered, including solution conditions, DBP assembly/aggregation, and RNA substrate properties. We discuss practical experimental methods for determination of DBP-RNA-nucleotide binding affinities and stoichiometries, steady-state ATPase activity, ATP binding, hydrolysis and product release rate constants, and RNA unwinding. We present general methods to integrate and analyze this combination of experimental data to identify the preferred kinetic pathway of ATP utilization and linked dsRNA unwinding.

Intermolecular migratory insertion of unactivated olefins into palladium-nitrogen bonds. Steric and electronic effects on the rate of migratory insertion

We report a detailed examination of the effect of the steric and electronic properties of the ancillary ligand and the alkene reactant on the rate of migratory insertion of unactivated alkenes into the palladium-nitrogen bond of isolated palladium amido complexes. A series of THF-bound and THF-free amidopalladium complexes ligated by cyclometalated benzylphosphine ligands possessing varied steric and electronic properties were synthesized. The THF-free complexes react with ethylene at -50 °C to form olefin-amido complexes that were observed directly and that undergo migratory insertion, followed by β-hydride elimination to generate enamine products. The effect of the steric properties of the ancillary ligand on the binding of the alkene and the rate of migratory insertion were evaluated individually. The relative binding affinity of ethylene vs THF is larger for the less sterically hindered complex than for the more hindered complex, but the less hindered complex undergoes the insertion of ethylene more slowly than does the more hindered complex. These two changes in relative equilibrium and rate constants cause the rates of reaction of ethylene with the two THF-ligated species having different steric properties to be similar to each other. Reactions of the complexes containing electronically varied ancillary ligands showed that the more electron-poor complexes underwent the migratory insertion step faster than the more electron-rich complexes. Reactions of a THF-ligated palladium-amide with substituted vinylarenes showed that electron-poor vinylarenes reacted with the amido complex slightly faster than electron-rich vinylarenes. Separation of the energetics of binding and insertion indicate that the complex of an electron-rich vinylarene is more stable in this system than the complex of a more electron-poor vinylarene but that the insertion step of the bound, electron-rich vinylarene is slower than the insertion step with the bound, electron-poor vinylarene.

Exact and approximate solutions for the decades-old Michaelis-Menten equation: Progress-curve analysis through integrated rate equations

The Michaelis-Menten rate equation can be found in most general biochemistry textbooks, where the time derivative of the substrate is a hyperbolic function of two kinetic parameters (the limiting rate V, and the Michaelis constant K(M) ) and the amount of substrate. However, fundamental concepts of enzyme kinetics can be difficult to understand fully, or can even be misunderstood, by students when based only on the differential form of the Michaelis-Menten equation, and the variety of methods available to calculate the kinetic constants from rate versus substrate concentration "textbook data." Consequently, enzyme kinetics can be confusing if an analytical solution of the Michaelis-Menten equation is not available. Therefore, the still rarely known exact solution to the Michaelis-Menten equation is presented here through the explicit closed-form equation in terms of the Lambert W(x) function. Unfortunately, as the W(x) is not available in standard curve-fitting computer programs, the practical use of this direct solution is limited for most life-science students. Thus, the purpose of this article is to provide analytical approximations to the equation for modeling Michaelis-Menten kinetics. The elementary and explicit nature of these approximations can provide students with direct and simple estimations of kinetic parameters from raw experimental time-course data. The Michaelis-Menten kinetics studied in the latter context can provide an ideal alternative to the 100-year-old problems of data transformation, graphical visualization, and data analysis of enzyme-catalyzed reactions. Hence, the content of the course presented here could gradually become an important component of the modern biochemistry curriculum in the 21st century.

Kinetic analysis of the inhibition of matrix metalloproteinases: lessons from the study of tissue inhibitors of metalloproteinases

Tissue inhibitors of metalloproteinases (TIMPs) are a group of highly potent inhibitors of matrix metalloproteinases (MMPs) and disintegrin metalloproteinases (ADAMs). The high affinity and "tight-binding" nature of the inhibition of MMPs or ADAMs by TIMPs presents challenges for the determination of both equilibrium and dissociation rate constants of these inhibitory events. Methodologies that enable some of these challenges to be overcome are described in this chapter and represent valuable lessons for the in vitro assessment of MMP or ADAM inhibitors within a drug discovery context.

Enzyme catalyzed reactions: from experiment to computational mechanism reconstruction

The traditional experimental practice in enzyme kinetics involves the measurement of substrate or product concentrations as a function of time. Advances in computing have produced novel approaches for modeling enzyme catalyzed reactions from time course data. One example of such an approach is the selection of appropriate chemical reactions that best fit the data. A common limitation of this approach resides in the number of chemical species considered. The number of possible chemical reactions grows exponentially with the number of chemical species, which makes difficult to select reactions that uniquely describe the data and diminishes the efficiency of the methods. In addition, a method's performance is also dependent on several quantitative and qualitative properties of the time course data, of which we know very little. This information is important to experimentalists as it could allow them to setup their experiments in ways that optimize the network reconstruction. We have previously described a method for inferring reaction mechanisms and kinetic rate parameters from time course data. Here, we address the limitations in the number of chemical reactions by allowing the introduction of information about chemical interactions. We also address the unknown properties of the input data by determining experimental data properties that maximize our method's performance. We investigate the following properties: initial substrate-enzyme concentration ratios; initial substrate-enzyme concentration variation ranges; number of data points; number of different experiments (time courses); and noise. We test the method using data generated in silico from the Michaelis-Menten and the Hartley-Kilby reaction mechanisms. Our results demonstrate the importance of experimental design for time course assays that has not been considered in experimental protocols. These considerations can have far reaching implications for the computational mechanism reconstruction process.

Elucidating the kinetic mechanism of DNA polymerization catalyzed by Sulfolobus solfataricus P2 DNA polymerase B1

Transient-state kinetic techniques were used to resolve the kinetic mechanism of DNA polymerization catalyzed by an exonuclease-deficient mutant of Sulfolobus solfataricus P2 DNA polymerase B1 (PolB1 exo-). Here, we report the kinetic parameters of several elementary steps for the forward polymerization reaction. PolB1 exo- binds tightly to DNA (K(d)(DNA) = 1.8 nM) and a correct incoming nucleotide (apparent K(d)(dTTP) = 11 microM). Moreover, several lines of kinetic evidence suggested that correct nucleotide incorporation catalyzed by PolB1 exo- was limited by a protein conformational change which precedes the chemistry step. The utilization of an "induced fit" mechanism by PolB1 exo- was supported by the following: a small, alpha-thio elemental effect of 1.5, varying DNA dissociation rates for the binary complex (0.043 s(-1)) as well as ternary complexes before (0.18 s(-1)) and after (0.0071 s(-1)) a conformational change, a greater amplitude for the pulse-chase than the pulse-quench reaction, and an activation energy barrier of 38 kcal/mol which is greater than the predicted values of phosphodiester bond formation both in solution and within a polymerase active site. Lastly, PolB1 exo- exhibited a low processivity value of 15, thereby suggesting a protein cofactor confers this replicative DNA polymerase with higher processivity in vivo.

Pre-tRNA turnover catalyzed by the yeast nuclear RNase P holoenzyme is limited by product release

Ribonuclease P (RNase P) is a ribonucleoprotein that catalyzes the 5' maturation of precursor transfer RNA in the presence of magnesium ions. The bacterial RNase P holoenzyme consists of one catalytically active RNA component and a single essential but catalytically inactive protein. In contrast, yeast nuclear RNase P is more complex with one RNA subunit and nine protein subunits. We have devised an affinity purification protocol to gently and rapidly purify intact yeast nuclear RNase P holoenzyme for transient kinetic studies. In pre-steady-state kinetic studies under saturating substrate concentrations, we observed an initial burst of tRNA formation followed by a slower, linear, steady-state turnover, with the burst amplitude equal to the concentration of the holoenzyme used in the reaction. These data indicate that the rate-limiting step in turnover occurs after pre-tRNA cleavage, such as mature tRNA release. Additionally, the steady-state rate constants demonstrate a large dependence on temperature that results in nonlinear Arrhenius plots, suggesting that a kinetically important conformational change occurs during catalysis. Finally, deletion of the 3' trailer in pre-tRNA has little or no effect on the steady-state kinetic rate constants. These data suggest that, despite marked differences in subunit composition, the minimal kinetic mechanism for cleavage of pre-tRNA catalyzed by yeast nuclear RNase P holoenzyme is similar to that of the bacterial RNase P holoenzyme.

Global kinetic explorer: a new computer program for dynamic simulation and fitting of kinetic data

We describe a new dynamic kinetic simulation program that allows multiple data sets to be fit simultaneously to a single model based on numerical integration of the rate equations describing the reaction mechanism. Unlike other programs that allow fitting based on numerical integration of rate equations, in the dynamic simulation rate constants, output factors, and starting concentrations of reactants can be scrolled while observing the change in the shape of the simulated reaction curves. Fast dynamic simulation facilitates the exploration of initial parameters that serve as the starting point for nonlinear regression in fitting data and facilitates exploration of the relationships between individual constants and observable reactions. The exploration of parameter space by dynamic simulation provides a powerful tool for learning kinetics and for evaluating the extent to which parameters are constrained by the data. This feature is critical to avoid overly complex models that are not supported by the data.

FitSpace explorer: an algorithm to evaluate multidimensional parameter space in fitting kinetic data

Fitting several sets of kinetic data directly to a model based on numerical integration provides the best method to extract kinetic parameters without relying on the simplifying assumptions required to achieve analytical solutions of rate equations. However, modern computer programs make it too easy to enter an overly complex model, and standard error analysis grossly underestimates errors when a system is underconstrained and fails to reveal the full degree to which multiple parameters are linked through the complex relationships common in kinetic data. Here we describe the application of confidence contour analysis obtained by measuring the dependence of the sum square error on each pair of parameters while allowing all remaining parameters to be adjusted in seeking the best fit. The confidence contours reveal complex relationships between parameters and clearly outline the space over which parameters can vary (the "FitSpace"). The utility of the method is illustrated by examples of well-constrained fits to published data on tryptophan synthase and the kinetics of oligonucleotide binding to a ribozyme. In contrast, analysis of alanine racemase clearly refutes claims that global analysis of progress curves can be used to extract the free energy profiles of enzyme-catalyzed reactions.

Uses and misuses of progress curve analysis in enzyme kinetics

Progress curve analysis is a convenient tool for the characterization of enzyme action: a single reaction mixture provides multiple experimental measured points for continuously varying amounts of substrates and products with exactly the same enzyme and modulator concentrations. The determination of kinetic parameters from the progress curves, however, requires complex mathematical evaluation of the time-course data. Some freely available programs (e.g. FITSIM, DYNAFIT) are widely applied to fit kinetic parameters to user-defined enzymatic mechanisms, but users often overlook the stringent requirements of the analytic procedures for appropriate design of the input experiments. Flaws in the experimental setup result in unreliable parameters with consequent misinterpretation of the biological phenomenon under study. The present commentary suggests some helpful mathematical tools to improve the analytic procedure in order to diagnose major errors in concept and design of kinetic experiments.

Use of fluorescence spectroscopy to elucidate RNA folding pathways

This overview unit discusses fluorescence spectroscopy as a tool for studying RNA folding. Ribozymes and oligonucleotides can be labeled with a fluorescent probe and analyzed to give information about both slow and fast kinetic processes with real-time data acquisition. The unit discusses the advantages and disadvantages of various pendant probes and nucleotide analogs, the analytical methods that can be used, instrument setup, control experiments, and a variety of kinetic experiments that can be performed, such as determination of rate constants.

Locating the rate-determining step(s) for three-step hydrolase-catalyzed reactions with DYNAFIT

Hydrolytic reactions of oligopeptide 4-nitroanilides catalyzed by human-alpha-thrombin, human activated protein C and human factor Xa were studied at pH 8.0-8.4 and 25.0+/-0.1 degrees C by the progress curve method and individual rate constants were calculated mostly within 10% internal error using DYNAFITV. A systematic strategy has been developed for fitting a three-step consecutive mechanism to eighteen hundred to six thousand time-course data points polled from two to four independent kinetic experiments. Enzyme and substrate concentrations were also calculated. Individual rate constants well reproduce published values obtained under comparable conditions and the Michaelis-Menten kinetic parameters calculated from these elementary rate constants are also within reasonable limits of published values. For comparison, the integrated Michaelis-Menten equation was also fitted to data from twelve sets. Both the k(cat) and k(cat)/K(m) values are within 15% agreement with those calculated using the elementary rate constants obtained with DYNAFITV. Rate constants for the second and third consecutive steps are within 3-4 fold indicating that both determine the overall rate. The Factor Xa-catalyzed hydrolysis of N-alpha-Z-D-Arg-Gly-Arg-pNA.2HCl at pH 8.4 in a series of buffers containing increasing fractions of deuterium at 25.0+/-0.1 degrees C shows a very strong dependence of k(3) and a moderate dependence of k(2) on D content in the buffer: the fractionation factors are: 0.49+/-0.03 for K(1,) 0.70+/-0.05 for k(2), and (0.32+/-0.03)(2) for k(3).

A deep cavitand catalyzes the Diels-Alder reaction of bound maleimides

A deep cavitand catalyzes Diels-Alder reactions of bound maleimides via activation of the dienophile by interaction with the organized hydrogen bonding network at the cavitand rim. Rapid in-out exchange of reactant and product allows efficient turnover. The increase in steric bulk of the reaction product lessens its binding affinity, reducing (and in some cases completely eliminating) product inhibition.

Assembly dynamics of Mycobacterium tuberculosis FtsZ

We have investigated the assembly of FtsZ from Mycobacterium tuberculosis (MtbFtsZ). Electron microscopy confirmed the previous observation that MtbFtsZ assembled into long, two-stranded filaments at pH 6.5. However, we found that assembly at pH 7.2 or 7.7 produced predominantly short, one-stranded protofilaments, similar to those of Escherichia coli FtsZ (EcFtsZ). Near pH 7, which is close to the pH of M. tuberculosis cytoplasm, MtbFtsZ formed a mixture of single- and two-stranded filaments. We developed a fluorescence resonance energy transfer assay to measure the kinetics of initial assembly and the dynamic properties at steady state. Assembly of MtbFtsZ reached a plateau after 60-100 s, about 10 times slower than EcFtsZ. The initial assembly kinetics were similar at pH 6.5 and 7.7, despite the striking difference in the polymer structures. Both were fit with a cooperative assembly mechanism involving a weak dimer nucleus, similar to EcFtsZ but with slower kinetics. Subunit turnover and GTPase at steady state were also about 10 times slower for MtbFtsZ than for EcFtsZ. Specifically, the half-time for subunit turnover in vitro at pH 7.7 was 42 s for MtbFtsZ compared with 5.5 s for EcFtsZ. Photobleaching studies in vivo showed a range of turnover half-times with an average of 25 s for MtbFtsZ as compared with 9 s for EcFtsZ.

Kinetics of binding the mRNA cap analogues to the translation initiation factor eIF4E under second-order reaction conditions

The kinetics of binding of five analogues of the 5'-mRNA cap, differing in size and electric charge, to the eukaryotic initiation factor eIF4E, at 20 degrees C, pH 7.2, and ionic strength of 150 mM, were measured, after mixing solutions of comparable concentrations of the reagents, in a stopped-flow spectrofluorimeter. The registered stopped-flow signals were fitted using an efficient software package, called Dyna Fit, based on a numerical solution of the kinetic rate equations for assumed reaction mechanisms. One-, two-, and three-step binding models were considered. The quality of fits for these models were compared using two statistical criteria: Akaike's Information Criterion and Bayesian Information Criterion. Based on resulting probabilities of the models, it was concluded that for all investigated ligands a one-step binding model has essentially no support in the experimental observations. Our conclusions were also analysed from the perspective of kinetic transients obtained for cap-eIF4E systems under the so called pseudo-first order reaction condition, which result in the linear correlation of the observed association rate constant with ligand concentration. The existence of such a linear correlation is usually considered as proof of a one-step binding mechanism. The kinetic and optical parameters, derived from fitting a two-step cap-binding model with the DynaFit, were used to simulate kinetic transients under pseudo-first order reaction conditions. It appeared that the observed association rate constants derived from these simulated transients are also linearly correlated with the ligand concentration. This indicated that these linear dependencies are not sufficient to conclude a one-step binding.

Transition state analysis of acid-catalyzed dAMP hydrolysis

Multiple kinetic isotope effects (KIEs) on deoxyadenosine monophosphate (dAMP) hydrolysis in 0.1 M HCl were used to determine the transition state (TS) structure and probe its intrinsic reactivity. The experimental KIEs revealed a stepwise (SN1) mechanism, with a discrete oxacarbenium ion intermediate. This is the first direct evidence for the deoxyribosyl oxacarbenium ion in solution. In 50% methanol/0.1 M HCl the products were deoxyribose 5-phosphate (dRMP) and alpha- and beta-methyl dRMP. The alpha-Me-dRMP/beta-Me-dRMP ratio was 8.5:1. Assuming that a free oxacarbenium ion is equally susceptible to nucleophilic attack on either face, this indicated that approximately 20% proceeded through a solvent-separated ion pair complex, or free oxacarbenium ion, a DN+AN mechanism, while approximately 80% of the reaction proceeded through a contact ion pair complex. The oxacarbenium ion lifetime was estimated at 10(-11)-10(-10) s. Computational transition states were found for ANDN, DN*AN, DN*AN, and DN+AN mechanisms using hybrid density functional theory calculations. After taking into account 20% of DN+AN, there was an excellent match of calculated to experimental KIEs for 80% of the reaction having a DN*AN mechanism. That is, C-N bond cleavage is reversible, with dAMP and the {oxacarbenium ion*adenine} complex in equilibrium. The first irreversible step is water attack on the oxacarbenium ion. The calculated 1'-14C KIE for a stepwise mechanism with irreversible C-N bond cleavage (DN*AN) was 1.052, in the range previously associated only with ANDN transition states, and close to the calculated ANDN value, 1.059. The 1'-14C KIE was strongly dependent on the adenine protonation state.

DEAD-box-protein-assisted RNA structure conversion towards and against thermodynamic equilibrium values

RNAs in biological processes often interconvert between defined structures. These RNA structure conversions are assisted by proteins and are frequently coupled to ATP hydrolysis. It is not well understood how proteins coordinate RNA structure conversions and which role ATP hydrolysis has in these processes. Here, we have investigated in vitro how the DEAD-box ATPase Ded1 facilitates RNA structure conversions in a simple model system. We find that Ded1 assists RNA structure conversions via two distinct pathways. One pathway requires ATP hydrolysis and involves the complete disassembly of the RNA strands. This pathway represents a kinetically controlled steady state between the RNA structures, which allows formation of less stable from more stable RNA conformations and thus RNA structure conversion against thermodynamic equilibrium values. The other pathway is ATP-independent and proceeds via multipartite intermediates that are stabilized by Ded1. Our results provide a basic mechanistic framework for protein-assisted RNA structure conversions that illuminates the role of ATP hydrolysis and reveal an unexpected diversity of pathways.

Enolpyruvyl Activation by Enolpyruvylshikimate-3-phosphate Synthase

Enolpyruvylshikimate-3-phosphate synthase (AroA, also called EPSP synthase) is a carboxyvinyl transferase involved in aromatic amino acid biosynthesis, forming EPSP from shikimate 3-phosphate and phosphoenolpyruvate. Upon extended incubation, EPSP ketal, a side product, forms by intramolecular nucleophilic addition of O4 to C2' of the enolpyruvyl group. The catalytic significance of this reaction was unclear, as it was initially proposed to arise from nonenzymatic breakdown of tetrahedral intermediate that had dissociated from AroA. This study shows that EPSP ketal formed in AroA's active site, not nonenzymatically, by demonstrating its formation in the presence of excess AroA. It formed both in the normal reaction and during AroA-catalyzed EPSP hydrolysis. In addition, nonenzymatic EPSP hydrolysis was studied to elucidate the catalytic imperative for enolpyruvyl reactions. Hydrolysis was acid-catalyzed, with a rate enhancement of >5 x 10(8)-fold. There was no detectable EPSP breakdown after 16 days at 90 degrees C in 1 M KOH, a solution that is 1000-fold more nucleophilic than neutral aqueous solutions. Thus, an unactivated enolpyruvyl group is not susceptible to nucleophilic attack. Enzymatic EPSP ketal formation therefore requires enolpyruvyl activation through protonation of C3' to form either a cationic intermediate or a highly cation-like transition state. Forming an EPSP cation requires the investment of considerable catalytic power by AroA. Such an intermediate is a potential target motif for inhibitor design.

Diastereoselective Oxidative Addition of Dihydrogen to IrI(CO)((R)-BINAP) and [Ir(CO)2((R)-BINAP)][SbF6]

Two chiral iridium(I) (R)-BINAP complexes, IrI(CO)((R)-BINAP) (1) and [Ir(CO)2((R)-BINAP)][SbF6] (2), have been synthesized and characterized, and their reactivity with dihydrogen has been studied. Complex 1 is formed on the addition of (R)-BINAP to [Bu4N][IrI2(CO)2] in toluene, and 2 is generated by the addition of AgSbF6 to a solution of 1 in dichloromethane under CO. A structure determination of complex 2 confirms a square planar coordination geometry, while that of 1 reveals a significant tetrahedral distortion from the expected planar coordination. Additionally, the structure of 1 shows a disorder between iodide and CO ligands. The reaction of 1 with H2 proceeds under kinetic control and shows a high degree of kinetic and thermodynamic selectivity; the kinetic product is formed by H2 addition across the P-Ir-CO axis of IrI(CO)((R)-BINAP) and yields two diastereomers which then convert over time to two more stable diastereomers which correspond to oxidative addition across the P-Ir-I axis. The kinetically favored diastereomers are formed in an initial ratio of 8.6:1, corresponding to a DeltaDeltaG* of 1.27 kcal/mol. The reaction of H2 with the C2-symmetric complex 2 also leads to the formation of two diastereomers, with one favored over the other kinetically by a 9.9:1 ratio on extrapolation to t = 0. When these reactions are followed using parahydrogen NMR methods, only one of the initially formed diastereomers in each case is found to exhibit substantial parahydrogen-induced polarization in the hydride resonances at room temperature.

Transient kinetic experiments demonstrate the existence of a unique catalytic enzyme form in the peptide-stimulated ATPase mechanism of Escherichia coli Lon protease

Lon is an oligomeric serine protease whose proteolytic activity is mediated by ATP hydrolysis. Although each monomeric subunit has an identical sequence, Lon contains two types of ATPase sites that hydrolyze ATP at drastically different rates. The catalytic low-affinity sites display pre-steady-state burst kinetics and hydrolyze ATP prior to peptide cleavage. The high-affinity sites are able to hydrolyze ATP at a very slow rate. By utilizing the differing Kd's, the high-affinity site can be blocked with unlabeled nucleotide while the activity at the low-affinity site is monitored. Little kinetic data are available that describe microscopic events along the reaction pathway of Lon. In this study we utilize MANT-ATP, a fluorescent analogue of ATP, to monitor the rate constants for binding of ATP as well as the release of ADP from Escherichia coli Lon protease. All of the adenine nucleotides tested bound to Lon on the order of 10(5) M(-1) s(-1), and the previously proposed conformational change associated with nucleotide binding was also detected. On the basis of the data obtained in this study we propose that the rate of ADP release is slightly different for the two ATPase sites. As the model peptide substrate [S2; YRGITCSGRQK(Bz)] [Thomas-Wohlever, J., and Lee, I. (2002) Biochemistry 41, 9418-9425] or the protein substrate casein affects only the steady-state ATPase activity of the low-affinity sites, we propose that Lon adopts a different form after its first turnover as an ATP-dependent protease. Based on the obtained rate constants, a revised kinetic model is presented for ATPase activity in Lon protease in both the absence and presence of the model peptide substrate (S2).

Kinetics of ATP-stimulated nuclease activity of the Escherichia coli RecBCD enzyme

The RecBCD enzyme is an ATP-dependent nuclease on both single-stranded and double-stranded DNA substrates. We have investigated the kinetics of the RecBCD-catalyzed reaction with small, single-stranded oligodeoxyribonucleotide substrates under single-turnover conditions using rapid-quench flow techniques. RecBCD-DNA complexes were allowed to form in pre-incubation mixtures. The nuclease reactions were initiated by mixing with ATP. The reaction time-courses were fit to several possible reaction mechanisms and quantitative estimates were obtained for rate constants for individual reaction steps. The relative rates of forward reaction versus dissociation from the DNA, and the fact that inclusion of excess non-radiolabeled single-stranded DNA to trap free RecBCD has no effect on the nuclease reaction, indicates that the reaction is processive. The reaction products show that the reaction begins near the 3'-end of the [5'-32P]DNA substrates and the major cleavage sites are two to four phosphodiester bonds apart. The product distribution is unchanged as the ATP concentration varies from 10 microM to 100 microM ATP, while the overall reaction rate varies by about tenfold. These observations suggest that DNA cleavage is tightly coordinated with movement of the enzyme along the DNA. The reaction time-courses at low concentrations of ATP (10 microM and 25 microM) have a significant lag before cleavage products appear. We propose that the lag represents ATP-dependent movement of the DNA from an initial binding site in the helicase domain of the RecB subunit to the nuclease active site in a separate domain of RecB. The extent of reaction of the substrate is limited (approximately 50%) under all conditions. This may indicate the formation of a non-productive RecBCD-DNA complex that does not dissociate in the 1-2 s time-scale of our experiments.

Role of the PH domain in regulating in vitro autophosphorylation events required for reconstitution of PDK1 catalytic activity

In addition to its catalytic domain, phosphoinsositide-dependent protein kinase-1 (PDK1) contains a C-terminal pleckstrin homology (PH) domain, which binds the membrane-bound phosphatidylinositol (3,4,5)-triphosphate [PI(3,4,5)P3] second messenger. Here, we report in vitro kinetic, phosphopeptide mapping, and oligomerization studies that address the role of the PH domain in regulating specific autophosphorylation events, which are required for PDK1 catalytic activation. First, 'inactive' unphosphorylated forms of N-terminal His6 tagged full length (His6-PDK1) and catalytic domain constructs [His6-PDK1(Delta PH)] were generated by treatment with Lambda protein phosphatase (lambda PP). Reconstitution of lambda PP-treated His6-PDK1(Delta PH) catalytic activity required activation loop Ser-241 phosphorylation, which occurred only upon trans-addition of 'active' PDK1 with an apparent bimolecular rate constant of (app)k1(S241) = 374+/-29 M(-1) s(-1). In contrast, full length lambda PP-treated His6-PDK1 catalyzed Ser-241 cis-autophosphorylation with an apparent first-order rate constant of (app)k1(S241) = (5.0+/-1.5) x 10(-4) s(-1) but remained 'inactive'. Reconstitution of lambda PP-treated His(6)-PDK1 catalytic activity occurred only when autophosphorylated in the presence of PI(3,4,5)P3 containing vesicles. PI(3,4,5)P3 binding to the PH domain activated apparent first-order Ser-241 autophosphorylation by 20-fold [(app)k1(S241) = (1.1+/-0.1) x 10(-2) s(-1)] and also promoted biphasic Thr-513 trans-autophosphorylation [(app)k2(T513) = (4.9+/-1.1) x 10(2) M(-1) s(-1) and(app)k3(T513) = (1.5+/-0.2) x 10(3) M(-1) s(-1)]. The results of mutagenesis studies suggest that Thr-513 phosphorylation may cause dissociation of autoinhibitory contacts formed between the contiguous regulatory PH and catalytic kinase domains.

Human mesotrypsin exhibits restricted S1' subsite specificity with a strong preference for small polar side chains

Mesotrypsin, an inhibitor-resistant human trypsin isoform, does not activate or degrade pancreatic protease zymogens at a significant rate. These observations led to the proposal that mesotrypsin is a defective digestive protease on protein substrates. Surprisingly, the studies reported here with alpha1-antitrypsin (alpha1AT) revealed that, even though mesotrypsin was completely resistant to this serpin-type inhibitor, it selectively cleaved the Lys10-Thr11 peptide bond at the N-terminus. Analyzing a library of alpha1AT mutants in which Thr11 was mutated to various amino acids, we found that mesotrypsin hydrolyzed lysyl peptide bonds containing Thr or Ser at the P1' position with relatively high specificity (kcat/KM approximately 10(5) m(-1) x s(-1)). Compared with Thr or Ser, P1' Gly or Met inhibited cleavage 13- and 25-fold, respectively, whereas P1' Asn, Asp, Ile, Phe or Tyr resulted in 100-200-fold diminished rates of proteolysis, and Pro abolished cleavage completely. Consistent with the Ser/Thr P1' preference, mesotrypsin cleaved the Arg358-Ser359 reactive-site peptide bond of alpha1AT Pittsburgh and was rapidly inactivated by the serpin mechanism (ka approximately 10(6) m(-1) s(-1)). Taken together, the results indicate that mesotrypsin is not a defective protease on polypeptide substrates in general, but exhibits a relatively high specificity for Lys/Arg-Ser/Thr peptide bonds. This restricted, thrombin-like subsite specificity explains why mesotrypsin cannot activate pancreatic zymogens, but might activate certain proteinase-activated receptors. The observations also identify alpha1AT Pittsburgh as an effective mesotrypsin inhibitor and the serpin mechanism as a viable stratagem to overcome the inhibitor-resistance of mesotrypsin.

Effects of Mutations and Truncations on the Kinetic Behavior of IIAGlc, a Phosphocarrier and Regulatory Protein of the Phosphoenolpyruvate Phosphotransferase System of Escherichia coli

IIAGlc, a component of the glucose-specific phosphoenolpyruvate:phosphotransferase system (PTS) of Escherichia coli, is important in regulating carbohydrate metabolism. In Glc uptake, the phosphotransfer sequence is: phosphoenolpyruvate --> Enzyme I --> HPr --> IIAGlc --> IICBGlc --> Glc. (HPr is the first phosphocarrier protein of the PTS.) We previously reported two classes of IIAGlc mutations that substantially decrease the P-transfer rate constants to/from IIAGlc. A mutant of His75 which adjoins the active site (His90), (H75Q), was 0.5% as active as wild-type IIAGlc in the reversible P-transfer to HPr. Two possible explanations were offered for this result: (a) the imidazole ring of His75 is required for charge delocalization and (b) H75Q disrupts the hydrogen bond network: Thr73, His75, phospho-His90. The present studies directly test the H-bond network hypothesis. Thr73 was replaced by Ser, Ala, or Val to eliminate the network. Because the rate constants for phosphotransfer to/from HPr were largely unaffected, we conclude that the H-bond network hypothesis is not correct. In the second class of mutants, proteolytic truncation of seven residues of the IIAGlc N terminus caused a 20-fold reduction in phosphotransfer to membrane-bound IICBGlc from Salmonella typhimurium. Here, we report the phosphotransfer rates between two genetically constructed N-terminal truncations of IIAGlc (Delta7 and Delta16) and the proteins IICBGlc and IIBGlc (the soluble cytoplasmic domain of IICBGlc). The truncations did not significantly affect reversible P-transfer to IIBGlc but substantially decreased the rate constants to IICBGlc in E. coli and S. typhimurium membranes. The results support the hypothesis (Wang, G., Peterkofsky, A., and Clore, G. M. (2000) J. Biol. Chem. 275, 39811-39814) that the N-terminal 18-residue domain "docks" IIAGlc to the lipid bilayer of membranes containing IICBGlc.

Heterodimeric capping protein from Arabidopsis is regulated by phosphatidic acid

The cytoskeleton is a key regulator of morphogenesis, sexual reproduction, and cellular responses to extracellular stimuli. Changes in the cellular architecture are often assumed to require actin-binding proteins as stimulus-response modulators, because many of these proteins are regulated directly by binding to intracellular second messengers or signaling phospholipids. Phosphatidic acid (PA) is gaining widespread acceptance as a major, abundant phospholipid in plants that is required for pollen tube tip growth and mediates responses to osmotic stress, wounding, and phytohormones; however, the number of identified effectors of PA is rather limited. Here we demonstrate that exogenous PA application leads to significant increases in filamentous actin levels in Arabidopsis suspension cells and poppy pollen grains. To investigate further these lipid-induced changes in polymer levels, we analyzed the properties of a key regulator of actin filament polymerization, the heterodimeric capping protein from Arabidopsis thaliana (AtCP). AtCP binds to PA with a K(d) value of 17 muM and stoichiometry of approximately 1:2. It also binds well to PtdIns(4,5)P(2), but not to several other phosphoinositide or acidic phospholipids. The interaction with PA inhibited the actin-binding activity of CP. In the presence of PA, CP is unable to block the barbed or rapidly growing and shrinking end of actin filaments. Precapped filament barbed ends can also be uncapped by addition of PA, allowing rapid filament assembly from an actin monomer pool that is buffered with profilin. The findings support a model in which the inhibition of CP activity in cells by elevated PA results in the stimulation of actin polymerization from a large pool of profilin-actin. Such regulation may be important for the response of plant cells to extracellular stimuli as well as for the normal process of pollen tube tip growth.

Chymotrypsin C (caldecrin) stimulates autoactivation of human cationic trypsinogen

Trypsin-mediated trypsinogen activation (autoactivation) facilitates digestive zymogen activation in the duodenum but may precipitate pancreatitis if it occurs prematurely in the pancreas. Autoactivation of human cationic trypsinogen is inhibited by a repulsive electrostatic interaction between the unique Asp218 on the surface of cationic trypsin and the conserved tetra-aspartate (Asp19-22) motif in the trypsinogen activation peptide (Nemoda, Z., and Sahin-Tóth, M. (2005) J. Biol. Chem. 280, 29645-29652). Here we describe that this interaction is regulated by chymotrypsin C (caldecrin), which can specifically cleave the Phe18-Asp19 peptide bond in the trypsinogen activation peptide and remove the N-terminal tripeptide. In contrast, chymotrypsin B, elastase 2A, or elastase 3A (proteinase E) are ineffective. Autoactivation of N-terminally truncated cationic trypsinogen is stimulated approximately 3-fold, and this effect is dependent on the presence of Asp218. Because chymotrypsinogen C is activated by trypsin, and chymotrypsin C stimulates trypsinogen activation, these reactions establish a positive feedback mechanism in the digestive enzyme cascade of humans. Furthermore, inappropriate activation of chymotrypsinogen C in the pancreas may contribute to the development of pancreatitis. Consistent with this notion, the pancreatitis-associated mutation A16V in cationic trypsinogen increases the rate of chymotrypsin C-mediated processing of the activation peptide 4-fold and causes accelerated trypsinogen activation in vitro.

Transient State Kinetics of Enzyme IICBGlc, a Glucose Transporter of the Phosphoenolpyruvate Phosphotransferase System of Escherichia coli

During translocation across the cytoplasmic membrane of Escherichia coli, glucose is phosphorylated by phospho-IIA(Glc) and Enzyme IICB(Glc), the last two proteins in the phosphotransfer sequence of the phosphoenolpyruvate:glucose phosphotransferase system. Transient state (rapid quench) methods were used to determine the second order rate constants that describe the phosphotransfer reactions (phospho-IIA(Glc) to IICB(Glc) to Glc) and also the second order rate constants for the transfer from phospho-IIA(Glc) to molecularly cloned IIB(Glc), the soluble, cytoplasmic domain of IICB(Glc). The rate constants for the forward and reverse phosphotransfer reactions between IIA(Glc) and IICB(Glc) were 3.9 x 10(6) and 0.31 x 10(6) m(-1) s(-1), respectively, and the rate constant for the physiologically irreversible reaction between [P]IICB(Glc) and Glc was 3.2 x 10(6) m(-1) s(-1). From the rate constants, the equilibrium constants for the transfer of the phospho-group from His90 of [P]IIA(Glc) to the phosphorylation site Cys of IIB(Glc) or IICB(Glc) were found to be 3.5 and 12, respectively. These equilibrium constants signify that the thiophospho-group in these proteins has a high phosphotransfer potential, similar to that of the phosphohistidinyl phosphotransferase system proteins. In these studies, preparations of IICB(Glc) were invariably found to contain endogenous, firmly bound Glc (estimated K'(D) approximately 10(-7) m). The bound Glc was kinetically competent and was rapidly phosphorylated, indicating that IICB(Glc) has a random order, Bi Bi, substituted enzyme mechanism. The equilibrium constant for the binding of Glc was deduced from differences in the statistical goodness of fit of the phosphotransfer data to the kinetic model.