Thermal inactivation and conformational lock studies on glucose oxidase

In this study, the dissociative thermal inactivation and conformational lock theories are applied for the homodimeric enzyme glucose oxidase (GOD) in order to analyze its structure. For this purpose, the rate of activity reduction of glucose oxidase is studied at various temperatures using b-D-glucose as the substrate by incubation of enzyme at various temperatures in the wide range between 40 and 70 �C using UV–Vis spectrophotometry. It was observed that in the two ranges of temperatures, the enzyme has two different forms. In relatively low temperatures, the enzyme is in its dimeric state and has normal activity. In high temperatures, the activity almost disappears and it aggregates. The above achievements are confirmed by dynamic light scattering. The experimental parameter ‘‘n’’ as the obvious number of conformational locks at the dimer interface of glucose oxidase is obtained by kinetic data, and the value is near to two. To confirm the above results, the X-ray crystallography structure of the enzyme, GOD (pdb, 1gal), was also studied. The secondary and tertiary structures of the enzyme to track the thermal inactivation were studied by circular dichroism and fluorescence spectroscopy, respectively. We proposed a mechanism model for thermal inactivation of GOD based on the absence of the monomeric form of the enzyme by circular dichroism and fluorescence spectroscopy

Detailed Enzyme Kinetics in Terms of Biochemical Species: Study of Citrate Synthase

The compulsory-ordered ternary catalytic mechanism for two-substrate two-product enzymes is analyzed to account for binding of inhibitors to each of the four enzyme states and to maintain the relationship between the kinetic constants and the reaction equilibrium constant. The developed quasi-steady flux expression is applied to the analysis of data from citrate synthase to determine and parameterize a kinetic scheme in terms of biochemical species, in which the effects of pH, ionic strength, and cation binding to biochemical species are explicitly accounted for in the analysis of the data. This analysis provides a mechanistic model that is consistent with the data that have been used support competing hypotheses regarding the catalytic mechanism of this enzyme

Dynamic Disorder in Quasi-Equilibrium Enzymatic Systems

Conformations and catalytic rates of enzymes fluctuate over a wide range of timescales. Despite these fluctuations, there exist some limiting cases in which the enzymatic catalytic rate follows the macroscopic rate equation such as the Michaelis-Menten law. In this paper we investigate the applicability of macroscopic rate laws for fluctuating enzyme systems in which catalytic transitions are slower than ligand binding-dissociation reactions. In this quasi-equilibrium limit, for an arbitrary reaction scheme we show that the catalytic rate has the same dependence on ligand concentrations as obtained from mass-action kinetics even in the presence of slow conformational fluctuations. These results indicate that the timescale of conformational dynamics – no matter how slow – will not affect the enzymatic rate in quasi-equilibrium limit. Our numerical results for two enzyme-catalyzed reaction schemes involving multiple substrates and inhibitors further support our general theory

Computer Controlled Automated Assay for Comprehensive Studies of Enzyme Kinetic Parameters

Stability and biological activity of proteins is highly dependent on their physicochemical environment. The development of realistic models of biological systems necessitates quantitative information on the response to changes of external conditions like pH, salinity and concentrations of substrates and allosteric modulators. Changes in just a few variable parameters rapidly lead to large numbers of experimental conditions, which go beyond the experimental capacity of most research groups. We implemented a computer-aided experimenting framework (“robot lab assistant”) that allows us to parameterize abstract, human-readable descriptions of micro-plate based experiments with variable parameters and execute them on a conventional 8 channel liquid handling robot fitted with a sensitive plate reader. A set of newly developed R-packages translates the instructions into machine commands, executes them, collects the data and processes it without user-interaction. By combining script-driven experimental planning, execution and data-analysis, our system can react to experimental outcomes autonomously, allowing outcome-based iterative experimental strategies. The framework was applied in a response-surface model based iterative optimization of buffer conditions and investigation of substrate, allosteric effector, pH and salt dependent activity profiles of pyruvate kinase (PYK). A diprotic model of enzyme kinetics was used to model the combined effects of changing pH and substrate concentrations. The 8 parameters of the model could be estimated from a single two-hour experiment using nonlinear least-squares regression. The model with the estimated parameters successfully predicted pH and PEP dependence of initial reaction rates, while the PEP concentration dependent shift of optimal pH could only be reproduced with a set of manually tweaked parameters. Differences between model-predictions and experimental observations at low pH suggest additional protonation-sites at the enzyme or substrates critical for enzymatic activity. The developed framework is a powerful tool to investigate enzyme reaction specifics and explore biological system behaviour in a wide range of experimental conditions

Substrate Cooperativity in Marine Luciferases

Marine luciferases are increasingly used as reporters to study gene regulation. These luciferases have utility in bioluminescent assay development, although little has been reported on their catalytic properties in response to substrate concentration. Here, we report that the two marine luciferases from the copepods, Gaussia princeps (GLuc) and Metridia longa (MLuc) were found, surprisingly, to produce light in a cooperative manner with respect to their luciferin substrate concentration; as the substrate concentration was decreased 10 fold the rate of light production decreased 1000 fold. This positive cooperative effect is likely a result of allostery between the two proposed catalytic domains found in Gaussia and Metridia. In contrast, the marine luciferases from Renilla reniformis (RLuc) and Cypridina noctiluca (CLuc) demonstrate a linear relationship between the concentration of their respective luciferin and the rate of light produced. The consequences of these enzyme responses are discussed

SBMLsqueezer: A CellDesigner plug-in to generate kinetic rate equations for biochemical networks

<p>Abstract</p> <p>Background</p> <p>The development of complex biochemical models has been facilitated through the standardization of machine-readable representations like SBML (Systems Biology Markup Language). This effort is accompanied by the ongoing development of the human-readable diagrammatic representation SBGN (Systems Biology Graphical Notation). The graphical SBML editor CellDesigner allows direct translation of SBGN into SBML, and vice versa. For the assignment of kinetic rate laws, however, this process is not straightforward, as it often requires manual assembly and specific knowledge of kinetic equations.</p> <p>Results</p> <p>SBMLsqueezer facilitates exactly this modeling step via automated equation generation, overcoming the highly error-prone and cumbersome process of manually assigning kinetic equations. For each reaction the kinetic equation is derived from the stoichiometry, the participating species (e.g., proteins, mRNA or simple molecules) as well as the regulatory relations (activation, inhibition or other modulations) of the SBGN diagram. Such information allows distinctions between, for example, translation, phosphorylation or state transitions. The types of kinetics considered are numerous, for instance generalized mass-action, Hill, convenience and several Michaelis-Menten-based kinetics, each including activation and inhibition. These kinetics allow SBMLsqueezer to cover metabolic, gene regulatory, signal transduction and mixed networks. Whenever multiple kinetics are applicable to one reaction, parameter settings allow for user-defined specifications. After invoking SBMLsqueezer, the kinetic formulas are generated and assigned to the model, which can then be simulated in CellDesigner or with external ODE solvers. Furthermore, the equations can be exported to SBML, LaTeX or plain text format.</p> <p>Conclusion</p> <p>SBMLsqueezer considers the annotation of all participating reactants, products and regulators when generating rate laws for reactions. Thus, for each reaction, only applicable kinetic formulas are considered. This modeling scheme creates kinetics in accordance with the diagrammatic representation. In contrast most previously published tools have relied on the stoichiometry and generic modulators of a reaction, thus ignoring and potentially conflicting with the information expressed through the process diagram. Additional material and the source code can be found at the project homepage (URL found in the Availability and requirements section).</p

Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug targets

BACKGROUND: Targeting persistent tubercule bacilli has become an important challenge in the development of anti-tuberculous drugs. As the glyoxylate bypass is essential for persistent bacilli, interference with it holds the potential for designing new antibacterial drugs. We have developed kinetic models of the tricarboxylic acid cycle and glyoxylate bypass in Escherichia coli and Mycobacterium tuberculosis, and studied the effects of inhibition of various enzymes in the M. tuberculosis model. RESULTS: We used E. coli to validate the pathway-modeling protocol and showed that changes in metabolic flux can be estimated from gene expression data. The M. tuberculosis model reproduced the observation that deletion of one of the two isocitrate lyase genes has little effect on bacterial growth in macrophages, but deletion of both genes leads to the elimination of the bacilli from the lungs. It also substantiated the inhibition of isocitrate lyases by 3-nitropropionate. On the basis of our simulation studies, we propose that: (i) fractional inactivation of both isocitrate dehydrogenase 1 and isocitrate dehydrogenase 2 is required for a flux through the glyoxylate bypass in persistent mycobacteria; and (ii) increasing the amount of active isocitrate dehydrogenases can stop the flux through the glyoxylate bypass, so the kinase that inactivates isocitrate dehydrogenase 1 and/or the proposed inactivator of isocitrate dehydrogenase 2 is a potential target for drugs against persistent mycobacteria. In addition, competitive inhibition of isocitrate lyases along with a reduction in the inactivation of isocitrate dehydrogenases appears to be a feasible strategy for targeting persistent mycobacteria. CONCLUSION: We used kinetic modeling of biochemical pathways to assess various potential anti-tuberculous drug targets that interfere with the glyoxylate bypass flux, and indicated the type of inhibition needed to eliminate the pathogen. The advantage of such an approach to the assessment of drug targets is that it facilitates the study of systemic effect(s) of the modulation of the target enzyme(s) in the cellular environment

Onset dynamics of type A botulinum neurotoxin-induced paralysis

Experimental studies have demonstrated that botulinum neurotoxin serotype A (BoNT/A) causes flaccid paralysis by a multi-step mechanism. Following its binding to specific receptors at peripheral cholinergic nerve endings, BoNT/A is internalized by receptor-mediated endocytosis. Subsequently its zinc-dependent catalytic domain translocates into the neuroplasm where it cleaves a vesicle-docking protein, SNAP-25, to block neurally evoked cholinergic neurotransmission. We tested the hypothesis that mathematical models having a minimal number of reactions and reactants can simulate published data concerning the onset of paralysis of skeletal muscles induced by BoNT/A at the isolated rat neuromuscular junction (NMJ) and in other systems. Experimental data from several laboratories were simulated with two different models that were represented by sets of coupled, first-order differential equations. In this study, the 3-step sequential model developed by Simpson (J Pharmacol Exp Ther 212:16–21,1980) was used to estimate upper limits of the times during which anti-toxins and other impermeable inhibitors of BoNT/A can exert an effect. The experimentally determined binding reaction rate was verified to be consistent with published estimates for the rate constants for BoNT/A binding to and dissociating from its receptors. Because this 3-step model was not designed to reproduce temporal changes in paralysis with different toxin concentrations, a new BoNT/A species and rate (kS) were added at the beginning of the reaction sequence to create a 4-step scheme. This unbound initial species is transformed at a rate determined by kS to a free species that is capable of binding. By systematically adjusting the values of kS, the 4-step model simulated the rapid decline in NMJ function (kS ≥0.01), the less rapid onset of paralysis in mice following i.m. injections (kS = 0.001), and the slow onset of the therapeutic effects of BoNT/A (kS < 0.001) in man. This minimal modeling approach was not only verified by simulating experimental results, it helped to quantitatively define the time available for an inhibitor to have some effect (tinhib) and the relation between this time and the rate of paralysis onset. The 4-step model predicted that as the rate of paralysis becomes slower, the estimated upper limits of (tinhib) for impermeable inhibitors become longer. More generally, this modeling approach may be useful in studying the kinetics of other toxins or viruses that invade host cells by similar mechanisms, e.g., receptor-mediated endocytosis

A Linear Framework for Time-Scale Separation in Nonlinear Biochemical Systems

Cellular physiology is implemented by formidably complex biochemical systems with highly nonlinear dynamics, presenting a challenge for both experiment and theory. Time-scale separation has been one of the few theoretical methods for distilling general principles from such complexity. It has provided essential insights in areas such as enzyme kinetics, allosteric enzymes, G-protein coupled receptors, ion channels, gene regulation and post-translational modification. In each case, internal molecular complexity has been eliminated, leading to rational algebraic expressions among the remaining components. This has yielded familiar formulas such as those of Michaelis-Menten in enzyme kinetics, Monod-Wyman-Changeux in allostery and Ackers-Johnson-Shea in gene regulation. Here we show that these calculations are all instances of a single graph-theoretic framework. Despite the biochemical nonlinearity to which it is applied, this framework is entirely linear, yet requires no approximation. We show that elimination of internal complexity is feasible when the relevant graph is strongly connected. The framework provides a new methodology with the potential to subdue combinatorial explosion at the molecular level

Cloning and endogenous expression of a Eucalyptus grandis UDP-glucose dehydrogenase cDNA

UDP-glucose dehydrogenase (UGDH) catalyzes the oxidation of UDP-glucose (UDP-Glc) to UDP-glucuronate (UDP-GlcA), a key sugar nucleotide involved in the biosynthesis of plant cell wall polysaccharides. A full-length cDNA fragment coding for UGDH was cloned from the cambial region of 6-month-old E. grandis saplings by RT-PCR. The 1443-bp-ORF encodes a protein of 480 amino acids with a predicted molecular weight of 53 kDa. The recombinant protein expressed in Escherichia coli catalyzed the conversion of UDP-Glc to UDP-GlcA, confirming that the cloned cDNA encodes UGDH. The deduced amino acid sequence of the cDNA showed a high degree of identity with UGDH from several plant species. The Southern blot assay indicated that more than one copy of UGDH is present in Eucalyptus. These results were also confirmed by the proteomic analysis of the cambial region of 3- and 22-year-old E. grandis trees by 2-DE and LC-MS/MS, showing that at least two isoforms are present. The cloned gene is mainly expressed in roots, stem and bark of 6-month-old saplings, with a lower expression in leaves. High expression levels were also observed in the cambial region of 3- and 22-year-old trees. The results described in this paper provide a further view of the hemicellulose biosynthesis during wood formation in E. grandis