Graphic rules in steady and non-steady state enzyme kinetics

Graphic methods, when applied to enzyme kinetics, can provide a visually intuitive relation between calculations and reaction graphs. This will not only greatly raise the efficiency of calculations but also significantly help the analysis of enzyme kinetic mechanisms. In this paper, four graphic rules are presented. Rules 1-3 are established for steady state enzyme-catalyzed reaction systems and Rule 4 is for non-steady state ones. In comparison with conventional graphic methods which can only be applied to steady state systems, the present rules have the following merits. 1) Complicated and tedious calculations can be greatly simplified; for example, in calculating the concentrations of enzyme species for the bi-bi random mechanism, the calculation work can be reduced 8-fold compared with the King-Altman's method. 2) A great deal of wasted labor can be avoided; for example, in calculating the rate of product formation for the same mechanism, the operation of finding and removing the 96 reciprocally canceled terms is no longer needed because they automatically disappear during the derivation. 3) Final results can be easily and safely checked by a formula provided in each of the graphic rules. 4) Non-steady state systems can also be treated by the present graphic method; for example, applying Rule 4, one can directly write out the solution for a non-steady state enzyme-catalyzed system, without the need to follow more difficult and complicated operations to solve differential equations. The mathematical proofs of Rules 1-4 are given in Appendices A-D (in the Miniprint), respectively.     

ENZO: a web tool for derivation and evaluation of kinetic models of enzyme catalyzed reactions

We describe a web tool ENZO (Enzyme Kinetics), a graphical interface for building kinetic models of enzyme catalyzed reactions. ENZO automatically generates the corresponding differential equations from a stipulated enzyme reaction scheme. These differential equations are processed by a numerical solver and a regression algorithm which fits the coefficients of differential equations to experimentally observed time course curves. ENZO allows rapid evaluation of rival reaction schemes and can be used for routine tests in enzyme kinetics. It is freely available as a web tool, at http://enzo.cmm.ki.si.     

Graphical rules for enzyme-catalysed rate laws.

Two graphical rules for calculating the reaction rates of steady-state enzyme-catalysed systems are presented. According to Rule 1, the number of non-zero terms can be easily calculated, which is useful to check the results and avoid mistakes. Rule 2 provides a new graphical method, which will be much more effective than the existing methods in dealing with complicated problems. In addition, during the process of deriving Rule 1, the mathematical principles of the Wong-Hanes structural rule have been naturally obtained too.     

Systems approach to the study of drug transport across membranes using suspension cultures of mammalian cells V. Simultaneous passive transport and biosynthesis

A physical model is described for the simultaneous enzymatic bioconversion of a nonelectrolyte solute and the passive transport of both the solute and product of the enzymatic reaction out of cells in culture suspension. The plasma membrane is assumed to be the rate-determining transport barrier. This model provides the basis for the experimental design and analysis of the Michaelis-Menten kinetic parameters of simple enzymatic reactions in situ, the phenomenological transport parameters and other factors. The primary set of differential equations describing the quasisteady state rate of change in the concentration of the solute and product within the cell due to enzyme reaction and transport are given. These are nonlinear and must be solved by numerical methods. However, analytical mathematical expressions have been derived for various cases in the limit when the rate of enzymatic reaction is first or zero order.     

Determination the rate constants of some biexponential reactions

Reactions that are described by biexponential functions are typical for many biological processes. The kinetics of these reactions is described by transcendental irrational equations interconnecting the reagent concentrations, time and rate constants. Meantime, their graphical representation in the semi-logarithmic coordinates can be decomposed into two straight lines that intercept at some angle. New simple methods for asymptotic numerical solution of the equations describing these reactions are suggested. These methods permit determining the rate constants using the kinetic data of initial substance concentration, which transform into final product according to a two-component model, a sequential model or a competitive model.     

Mathematical modelling of glycolysis and adenine nucleotide metabolism of human erythrocytes. I. Reaction-kinetic statements, analysis of in vivo state and determination of starting conditions for in vitro experiments

A mathematical model of energy metabolism of human red cells is presented, which includes besides the glycolytic reactions the adenine nucleotide metabolism. The model is based on the network of chemical reactions, the thermodynamic equilibrium constants of fast reversible reactions and on the kinetic equations for irreversible enzyme reactions. The model consists of a system of 16 differential equations and allows the mathematical evaluation of metabolic levels in the steady state of energy metabolism corresponding to the in vivo state erythrocytes with the kinetic data for the enzymes derived from in vitro experiments. The dependence of the levels of metabolites in the steady state on the activity of some enzymes is analysed to characterize the regulatory properties of the system. The comparison of the steady state levels of the model with experimental data makes it possible to estimate values of some controversial enzyme parameters. Estimates of the kinetic parameters of the following intracellular processes are presented: 1) rate constant of AMP-phosphatase, 2) maximum rate of adenylate deaminase, 3) activity of adenine phosphoribosylpyrophosphate transferase and 4) adenosine transport through the cell membrane. The simulation of the preparatory phase before incubation of erythrocytes indicates, that the model also permits to compute the time course of changes of levels of metabolites. To solve the initial problem the stiff differential equation system is integrated numerically by an efficient program without the application of the quasi-steady-state approximation.     

Stochastic simulations on a model of circadian rhythm generation

Biological phenomena are often modeled by differential equations, where states of a model system are described by continuous real values. When we consider concentrations of molecules as dynamical variables for a set of biochemical reactions, we implicitly assume that numbers of the molecules are large enough so that their changes can be regarded as continuous and they are described deterministically. However, for a system with small numbers of molecules, changes in their numbers are apparently discrete and molecular noises become significant. In such cases, models with deterministic differential equations may be inappropriate, and the reactions must be described by stochastic equations. In this study, we focus a clock gene expression for a circadian rhythm generation, which is known as a system involving small numbers of molecules. Thus it is appropriate for the system to be modeled by stochastic equations and analyzed by methodologies of stochastic simulations. The interlocked feedback model proposed by Ueda et al. as a set of deterministic ordinary differential equations provides a basis of our analyses. We apply two stochastic simulation methods, namely Gillespie's direct method and the stochastic differential equation method also by Gillespie, to the interlocked feedback model. To this end, we first reformulated the original differential equations back to elementary chemical reactions. With those reactions, we simulate and analyze the dynamics of the model using two methods in order to compare them with the dynamics obtained from the original deterministic model and to characterize dynamics how they depend on the simulation methodologies.     

Effectiveness factor calculations for immobilized enzyme catalysts

The steady state, nonlinear diffusion equations which describe reactions in constrained enzyme solutions are of great interest in many biological and engineering applications. As in other types of nonlinear differential equations, exact analytical solutions do not exist except in some simplified cases. In this paper, a general procedure is presented for solving numerically for the substrate concentration profile and effectiveness factor utilizing the transformation method suggested by Na and Na. Design correlations for enzyme solutions constrained within spherical membranes are included. The use of a unique definition of the Thiele Modulus in these charts permits the clear illustration of the effects of substrate concentration and external mass transfer resistances on the overall effectiveness factor for the catalyst particle.     

Enzyme kinetics for a two-step enzymic reaction with comparable initial enzyme-substrate ratios

We extend the validity of the quasi-steady state assumption for a model double intermediate enzyme-substrate reaction to include the case where the ratio of initial enzyme to substrate concentration is not necessarily small. Simple analytical solutions are obtained when the reaction rates and the initial substrate concentration satisfy a certain condition. These analytical solutions compare favourably with numerical solutions of the full system of differential equations describing the reaction. Experimental methods are suggested which might permit the application of the quasi-steady state assumption to reactions where it may not have been obviously applicable before.     

Stochastic approach to molecular interactions and computational theory of metabolic and genetic regulations

The underlying molecular mechanisms of metabolic and genetic regulations are computationally identical and can be described by a finite state Markov process. We establish a common computational model for both regulations based on the stationary distribution of the Markov process with the aim of establishing a unified, quantitative model of general biological regulations. Various existing results regarding intracellular regulations are derived including the classical Michaelis-Menten equation and its generalization to more complex allosteric enzymes in a systematic way. The notion of probability flow is introduced to distinguish the equilibrium stationary distribution from the non-equilibrium one; it plays a crucial role in the analysis of stationary state equations. A graphical criterion to guarantee the existence of an equilibrium stationary distribution is derived, which turns out to be identical to the classical Wegscheider condition. Simple graphical methods to compute the equilibrium and non-equilibrium stationary distributions are derived based crucially on the probability flow, which dramatically simplifies the classical methods still used in enzymology.     

A common active site model for catalysis by chorismate mutase--prephenate dehydrogenase.

Comparison of the calculated structures for the transition states of the two reactions catalysed by chorismate mutase prephenate dehydrogenase suggests that both reactions could be catalysed at a common active site. Kinetic data for the enzyme from Aerobacteraerogenes are consistent with this possibility. On the basis of these theoretical and experimental data a model for a common active site is developed. In the model, the transition state for each reaction is bound to the enzyme via both of the two substrate carboxyl groups, and can also interact with the coenzyme nicotinamide adenine dinucleotide through a hydrogen bond between the amide moiety of the nicotinamide ring and the hydroxyl group of the substrate. Chorismate, prephenate and 4-hydroxyphenylpyruvate in their ground states form the same hydrogen bond to the coenzyme, but are bound to the enzyme via a single carboxyl group only. The additional bond formed between the enzyme and the transition state structures thus provides the transition state stabilization required for catalysis of both reactions.     

Introducing total substrates simplifies theoretical analysis at non-negligible enzyme concentrations: pseudo first-order kinetics and the loss of zero-order ultrasensitivity

The Briggs–Haldane standard quasi-steady state approximation and the resulting rate expressions for enzyme driven biochemical reactions provide crucial theoretical insight compared to the full set of equations describing the reactions, mainly because it reduces the number of variables and equations. When the enzyme is in excess of the substrate, a significant amount of substrate can be bound in intermediate complexes, so-called substrate sequestration. The standard quasi-steady state approximation is known to fail under such conditions, a main reason being that it neglects these intermediate complexes. Introducing total substrates, i.e., the sums of substrates and intermediate complexes, provides a similar reduction of the number of variables to consider but without neglecting the contribution from intermediate complexes. The present theoretical study illustrates the usefulness of such simplifications for the understanding of biochemical reaction schemes. We show how introducing the total substrates allows a simple analytical treatment of the relevance of significant enzyme concentrations for pseudo first-order kinetics and reconciles two proposed criteria for the validity of the pseudo first-order approximation. In addition, we show how the loss of zero-order ultrasensitivity in covalent modification cycles can be analyzed, in particular that approaches such as metabolic control analysis are immediately applicable to scenarios described by the total substrates with enzyme concentrations higher than or comparable to the substrate concentrations. A simple criterion which excludes the possibility of zero-order ultrasensitivity is presented.     

Steady state kinetics of consecutive enzyme catalysed reactions involving single substrates: procedures for the interpretation of coupled assays

Sets of differential rate equations are written describing a linear sequence of reactions occurring in solution each catalysed by a control enzyme or one of the Michaelis-Menten type. It is shown that the solutions of these equations may be formulated as a set of Maclaurin polynomials, expressing the concentration of each reactant and of final product as a function of time. From arrays of such polynomials, general expressions are induced for the first non-zero term of the series. These are used to formulate a procedure (illustrated with an example simulated by numerical integration) by which results of coupled enzymic assays may be analysed in terms of maximal velocities and apparent Michaelis constants: correlation is made with other established methods for conducting coupled assays. The present procedure assumes a steady state of enzyme-substrate complexes but not of intermediate reactants.     

A mathematical model for the branched chain amino acid biosynthetic pathways of Escherichia coli K12

As a first step toward the elucidation of the systems biology of the model organism Escherichia coli, it was our goal to mathematically model a metabolic system of intermediate complexity, namely the well studied end product-regulated pathways for the biosynthesis of the branched chain amino acids L-isoleucine, L-valine, and L-leucine. This has been accomplished with the use of kMech (Yang, C.-R., Shapiro, B. E., Mjolsness, E. D., and Hatfield, G. W. (2005) Bioinformatics 21, in press), a Cellerator (Shapiro, B. E., Levchenko, A., Meyerowitz, E. M., Wold, B. J., and Mjolsness, E. D. (2003) Bioinformatics 19, 677-678) language extension that describes a suite of enzyme reaction mechanisms. Each enzyme mechanism is parsed by kMech into a set of fundamental association-dissociation reactions that are translated by Cellerator into ordinary differential equations. These ordinary differential equations are numerically solved by Mathematica. Any metabolic pathway can be simulated by stringing together appropriate kMech models and providing the physical and kinetic parameters for each enzyme in the pathway. Writing differential equations is not required. The mathematical model of branched chain amino acid biosynthesis in E. coli K12 presented here incorporates all of the forward and reverse enzyme reactions and regulatory circuits of the branched chain amino acid biosynthetic pathways, including single and multiple substrate (Ping Pong and Bi Bi) enzyme kinetic reactions, feedback inhibition (allosteric, competitive, and non-competitive) mechanisms, the channeling of metabolic flow through isozymes, the channeling of metabolic flow via transamination reactions, and active transport mechanisms. This model simulates the results of experimental measurements.     

A variational approach to parameter estimation in ordinary differential equations

ABSTRACT: BACKGROUND: Ordinary differential equations are widely-used in the field of systems biology andchemical engineering to model chemical reaction networks. Numerous techniques havebeen developed to estimate parameters like rate constants, initial conditions or steady stateconcentrations from time-resolved data. In contrast to this countable set of parameters, theestimation of entire courses of network components corresponds to an innumerable set ofparameters. RESULTS: The approach presented in this work is able to deal with course estimation for extrinsicsystem inputs or intrinsic reactants, both not being constrained by the reaction networkitself. Our method is based on variational calculus which is carried out analytically toderive an augmented system of differential equations including the unconstrainedcomponents as ordinary state variables. Finally, conventional parameter estimation isapplied to the augmented system resulting in a combined estimation of courses andparameters. CONCLUSIONS: The combined estimation approach takes the uncertainty in input courses correctly intoaccount. This leads to precise parameter estimates and correct confidence intervals. Inparticular this implies that small motifs of large reaction networks can be analysedindependently of the rest. By the use of variational methods, elements from control theoryand statistics are combined allowing for future transfer of methods between the two fields.     

Description of the ratios between metabolite concentrations in a polyenzyme system

It was demonstrated that the relations between substrate and product concentrations for a reaction catalyzed by michaelian enzyme incorporated in a multienzyme system can be graphically represented by a diverging set of straight lines intersecting in one point, the flux velocity being treated as a parameter. A competitive inhibitor shifts the intersection point along the line of equilibrium state. The relations between the concentrations of more than two reagents are represented by a set of equivelocity surfaces. The relations between substrate and product concentrations for a kinetically cooperative reaction conforming to the graphical representation by the second--order curves were analyzed. The stability criterion was obtained for a multienzyme system with the first enzyme allosterically regulated by products of subsequent reactions.     

Pitfalls in the study of steady state kinetics of enzymes: spurious inhibition patterns due to stray light errors

When reaction velocity measurements of enzyme reactions are carried out with single beam, single monochromator spectrophotometers, stray light in the spectrophotometer can produce systematic errors in the apparent velocities when highly absorbing solutions (optical density >2.0) are used. These errors can give rise to spurious “inhibition” patterns of the steady state kinetics. Because of a suspected error of this kind, this laboratory has recently reinvestigated the kinetics of glucose 6-phosphate dehydrogenase from Escherichia coli and found that the reported noncompetitive inhibition of the enzyme by DPNH is explained more readily by an unnoticed effect of stray light on the apparent reaction velocity than by a true enzyme inhibition. Methods for estimating and correcting such errors in spectrophotometers are presented in detail.     

pH Feedback control of enzyme membranes

pH feedback on immobilized enzymes is theoretically examined with respect to substrate and pH levels, strength of acids produced by the reaction, buffering and asymmetry of the system. All the productions of proton by the different reactions are taken into account by using a ‘symbolic species’ H^*. The system of differential diffusion-reaction equations is then integrated using numerical methods. The local ‘effective enzyme activity’ modulated by an acidity factor enables us to predict and quantify evolutions of the systems: NonMichacIian behavior of an immobilized MichaeIis-Mentcn-type enzyme is shown, even when pH back-actions are excluded: the analysis of intramembranc pH profiles shows that the shift of the optimal pH is a complex function of the substrate and pH levels, the intrinsic pH dependence of the enzyme, and the membrane characteristics. This study may easily be transposed to other types of effector such as divalent cations and used in examining self-regulations of multienzyme systems where pH-active reactions are involved.     

Membrane transport models with fast and slow reactions: General analytical solution for a single relaxation

Summary Membrane transport models are usually expressed on the basis of chemical kinetics. The states of a transporter are related by rate constants, and the time-dependent changes of these states are given by linear differential equations of first order. To calculate the time-dependent transport equation, it is necessary to solve a system of differential equations which does not have a general analytical solution if there are more than five states. Since transport measurements in a complex system rarely provide all the time constants because some of them are too rapid, it is more appropriate to obtain approximate analytical solutions, assuming that there are fast and slow reaction steps. The states of the fast steps are related by equilibrium constants, thus permitting the elimination of their differential equations and leaving only those for the slow steps. With a system having only two slow steps, a single differential equation is obtained and the state equations have a single relaxation. Initial conditions for the slow reactions are determined after the perturbation which redistribute the states related by fast reactions. Current and zero-trans uptake equations are calculated. Curve fitting programs can be used to implement the general procedure and obtain the model parameters.