Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

Efficient and accurate experimental design for enzyme kinetics: Bayesian studies reveal a systematic approach

In areas such as drug development, clinical diagnosis and biotechnology research, acquiring details about the kinetic parameters of enzymes is crucial. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. We demonstrate that a Bayesian approach (the use of prior knowledge) can produce major gains quantifiable in terms of information, productivity and accuracy of each experiment. Developing the use of Bayesian Utility functions, we have used a systematic method to identify the optimum experimental designs for a number of kinetic model data sets. This has enabled the identification of trends between kinetic model types, sets of design rules and the key conclusion that such designs should be based on some prior knowledge of K-M and/or the kinetic model. We suggest an optimal and iterative method for selecting features of the design such as the substrate range, number of measurements and choice of intermediate points. The final design collects data suitable for accurate modelling and analysis and minimises the error in the parameters estimated. (C) 2003 Elsevier Science B.V. All rights reserved.

Prediction of Intrinsic kinetics in diffusion disguised immobilized enzyme systems

When immobilized enzyme kinetics are disguised by inter- and intraparticle diffusion effects, an approximate mathematical procedure is indicated whereby experimental data obtained in the limiting ranges of first- and zeroth-order Michaelis-Menten kinetics could be used for the prediction of the kinetic constants.

Stochastic chemical kinetics and the quasi-steady-state assumption : application to the stochastic simulation algorithm and chemical master equation

Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis–Menten enzyme kinetics, double phosphorylation, the Goldbeter–Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.

Studies on active site mutants of P-falciparum adenylosuccinate synthetase: Insights into enzyme catalysis and activation

Adenylosuccinate synthetase catalyzes a reversible reaction utilizing IMP, GTP and aspartate in the presence of Mg2+ to form adenylosuccinate, GDP and inorganic phosphate. Comparison of similarly liganded complexes of Plasmodium falciparum, mouse and Escherichia coil AdSS reveals H-bonding interactions involving nonconserved catalytic loop residues (Asn429, Lys62 and Thr307) that are unique to the parasite enzyme. Site-directed mutagenesis has been used to examine the role of these interactions in catalysis and structural organization of P. falciparum adenylosuccinate synthetase (PfAdSS). Mutation of Asn429 to Val, Lys62 to Leu and Thr307 to Val resulted in an increase in K-m values for IMP, GTP and aspartate, respectively along with a 5 fold drop in the k(cat) value for N429V mutant suggesting the role of these residues in ligand binding and/or catalysis. We have earlier shown that the glycolytic intermediate, fructose 1,6 bisphosphate, which is an inhibitor of mammalian AdSS is an activator of the parasite enzyme. Enzyme kinetics along with molecular docking suggests a mechanism for activation wherein F16BP seems to be binding to the Asp loop and inducing a conformation that facilitates aspartate binding to the enzyme active site. Like in other AdSS, a conserved arginine residue (Arg155) is involved in dimer crosstalk and interacts with IMP in the active site of the symmetry related subunit of PfAdSS. We also report on the iochemical characterization of the arginine mutants (R155L, R155K and R155A) which suggests that unlike in E. coil AdSS, Arg155 in PfAdSS influences both ligand binding and catalysis. (C) 2010 Elsevier B.V. All rights reserved.

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

Exploring the enzymatic landscape: distribution and kinetics of hydrolytic enzymes in soil particle-size fractions

The location of extracellular enzymes within the soil architecture and their association with the various soil components affects their catalytic potential. A soil fractionation study was carried out to investigate: (a) the distribution of a range of hydrolytic enzymes involved in C, N and P transformations, (b) the effect of the location on their respective kinetics, (c) the effect of long-term N fertilizer management on enzyme distribution and kinetic parameters. Soil (silty clay loam) from grassland which had received 0 or 200 kg N ha(-1) yr(-1) was fractionated, and four particle-size fractions (> 200, 200-63, 63-2 and 0. 1-2 mum) were obtained by a combination of wet-sieving and centrifugation, after low-energy ultrasonication. All fractions were assayed for four carbohydrases (beta-cellobiohydrolase, N-acetyl-beta-glucosammidase, beta-glucosidase and beta-xylosidase), acid phosphatase and leucine-aminopeptidase using a microplate fluorimetric assay based on MUB-substrates. Enzyme kinetics (V-max and K-m) were estimated in three particle-size fractions and the unfractionated soil. The results showed that not all particle-size fractions were equally enzymatically active and that the distribution of enzymes between fractions depended on the enzyme. Carbohydrases predominated in the coarser fractions while phosphatase and leucine-aminopeptidase were predominant in the clay-size fraction. The Michaelis constant (K.) varied among fractions, indicating that the association of the same enzyme with different particle-size fractions affected its substrate affinity. The same values of Km were found in the same fractions from the soil under two contrasting fertilizer management regimes, indicating that the Michaelis constant was unaffected by soil changes caused by N fertilizer management. (C) 2004 Elsevier Ltd. All rights reserved.

The Water Effect on the Kinetics of the Bovine Liver Catalase

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Kinetics and crystal structure of human purine nucleoside phosphorylase in complex with 7-methyl-6-thio-guanosine

Purine nucleoside phosphorylase (PNP) catalyzes the reversible phosphorolysis of nucleosides and deoxynucleosides, generating ribose 1-phosphate and the purine base, which is an important step of purine catabolism pathway. The lack of such an activity in humans, owing to a genetic disorder, causes T-cell impairment, and drugs that inhibit this enzyme may have the potential of being utilized as modulators of the immunological system to treat leukemia, autoimmune diseases, and rejection in organ transplantation. Here, we describe kinetics and crystal structure of human PNP in complex with 7-methyl-6-thio-guanosine, a synthetic substrate, which is largely used in activity assays. Analysis of the structure identifies different protein conformational changes upon ligand binding, and comparison of kinetic and structural data permits an understanding of the effects of atomic substitution on key positions of the synthetic substrate and their consequences to enzyme binding and catalysis. Such knowledge may be helpful in designing new PNP inhibitors. © 2005 Elsevier Inc. All rights reserved.

The Mode of Action of Recombinant Mycobacterium tuberculosis Shikimate Kinase: Kinetics and Thermodynamics Analyses

Tuberculosis remains as one of the main cause of mortality worldwide due to a single infectious agent, Mycobacterium tuberculosis. The aroK-encoded M. tuberculosis Shikimate Kinase (MtSK), shown to be essential for survival of bacilli, catalyzes the phosphoryl transfer from ATP to the carbon-3 hydroxyl group of shikimate (SKH), yielding shikimate-3-phosphate and ADP. Here we present purification to homogeneity, and oligomeric state determination of recombinant MtSK. Biochemical and biophysical data suggest that the chemical reaction catalyzed by monomeric MtSK follows a rapid-equilibrium random order of substrate binding, and ordered product release. Isothermal titration calorimetry (ITC) for binding of ligands to MtSK provided thermodynamic signatures of non-covalent interactions to each process. A comparison of steady-state kinetics parameters and equilibrium dissociation constant value determined by ITC showed that ATP binding does not increase the affinity of MtSK for SKH. We suggest that MtSK would more appropriately be described as an aroL-encoded type II shikimate kinase. Our manuscript also gives thermodynamic description of SKH binding to MtSK and data for the number of protons exchanged during this bimolecular interaction. The negative value for the change in constant pressure heat capacity (ΔCp) and molecular homology model building suggest a pronounced contribution of desolvation of non-polar groups upon binary complex formation. Thermodynamic parameters were deconvoluted into hydrophobic and vibrational contributions upon MtSK:SKH binary complex formation. Data for the number of protons exchanged during this bimolecular interaction are interpreted in light of a structural model to try to propose the likely amino acid side chains that are the proton donors to bulk solvent following MtSK:SKH complex formation. © 2013 Rosado et al.

Cyclodextrin Production by Bacillus lehensis Isolated from Cassava Starch: Characterisation of a Novel Enzyme

The properties of a previously unknown enzyme, denominated cyclodextrin glycosyltransferase, produced from Bacillus lehensis, were evaluated using affinity chromatography for protein purification. Enzyme characteristics (optimum pH and temperature; pH and temperature stability), the influence of substances on the enzyme activity, enzyme kinetics, and cyclodextrin production were analysed. Cyclodextrin glycosyltransferase was purified up to 320.74-fold by affinity chromatography using beta-cyclodextrin as the binder and it exhibited 8.71% activity recovery. This enzyme is a monomer with a molecular weight of 81.27 kDa, as estimated by SDS-PAGE. Optimum temperature and pH for cydodextrin glycosyltransferase were 55 degrees C and 8.0, respectively. The Michaelis-Menten constant was 8.62 g/l during maximum velocity of 0.858 g/l.h.

Modeling physiological processes in plankton on enzyme kinetic principles

[EN] Many ecologically important chemical transformations in the ocean are controlled by biochemical enzyme reactions in plankton. Nitrogenase regulates the transformation of N2 to ammonium in some cyanobacteria and serves as the entryway for N2 into the ocean biosphere. Nitrate reductase controls the reduction of NO3 to NO2 and hence new production in phytoplankton. The respiratory electron transfer system in all organisms links the carbon oxidation reactions of intermediary metabolism with the reduction of oxygen in respiration. Rubisco controls the fixation of CO2 into organic matter in phytoplankton and thus is the major entry point of carbon into the oceanic biosphere. In addition to these, there are the enzymes that control CO2 production, NH4 excretion and the fluxes of phosphate. Some of these enzymes have been recognized and researched by marine scientists in the last thirty years. However, until recently the kinetic principles of enzyme control have not been exploited to formulate accurate mathematical equations of the controlling physiological expressions. Were such expressions available they would increase our power to predict the rates of chemical transformations in the extracellular environment of microbial populations whether this extracellular environment is culture media or the ocean. Here we formulate from the principles of bisubstrate enzyme kinetics, mathematical expressions for the processes of NO3 reduction, O2 consumption, N2 fixation, total nitrogen uptake.

A windows program for the derivation of steady-state equations in enzyme systems

Despite the number of computer-assisted methods described for the derivation of steady-state equations of enzyme systems, most of them are focused on strict steady-state conditions or are not able to solve complex reaction mechanisms. Moreover, many of them are based on computer programs that are either not readily available or have limitations. We present here a computer program called WinStes, which derives equations for both strict steady-state systems and those with the assumption of rapid equilibrium, for branched or unbranched mechanisms, containing both reversible and irreversible conversion steps. It solves reaction mechanisms involving up to 255 enzyme species, connected by up to 255 conversion steps. The program provides all the advantages of the Windows programs, such as a user-friendly graphical interface, and has a short computation time. WinStes is available free of charge on request from the authors. (c) 2006 Elsevier Inc. All rights reserved.

Sequence-related structural DNA anomalies affect restriction enzyme activity and DNA polymerase I binding

DNA serves as a target molecule for several types of enzymes and may assume a wide variety of structural motifs depending upon the local sequence. The BssHII restriction site (GC)3 resides in a 9bp region of alternating pyrimidine and purine residues within the &phis;X174 genome. Such sequences are known to demonstrate non-canonical helical behavior under the appropriate conditions. The kinetics of BssHII cleavage was investigated in supercoiled and linear plasmid DNA, and in a 323bp DNA fragment obtained via amplification of &phis;X174. The rate of enzyme cleavage was enhanced in the supercoiled form and in the presence of 50μM cobalt hexamine. Similarly, cobalt hexamine was also found to enhance TaqI activity directly adjacent to the (GC)3 region. ^ Initial DNA polymerase I binding studies (including a gel mobility shift assay and a protection assay) indicated a notable interaction between DNA polymerase I and the BssHII site. An in-depth study revealed that equilibrium binding of DNA polymerase I to the T7 RNA polymerase promoter was comparable to that of the (GC)3 site, however the strongest interaction was observed with a cruciform containing region. Increasing the ionic strength of the solution environment, including the addition of DNA polymerase I reaction buffer significantly decreased the equilibrium dissociation constant values. ^ It is suggested that the region within or around the BssHII site experiences a conformational change generating a novel structure under the influence of supercoiled tension or 50μM cobalt hexamine. It is proposed that this transition may enhance enzyme activity and binding by providing an initial enzyme-docking site—the rate-limiting step in restriction enzyme kinetics. The high binding potential of DNA polymerase I for each of the motifs described, is hypothesized to be due to recognition of the structural DNA anomalies by the 3′–5′ exonuclease domain. ^