Open access peer-reviewed chapter

Mathematical Modeling and Simulation of Nonlinear Process in Enzyme Kinetics

By Lakshmanan Rajendran, Mohan Chitra Devi, Carlos Fernandez and Qiuming Peng

Submitted: March 27th 2017Reviewed: September 11th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.70914

Downloaded: 1169


A deep and analytical understanding of the enzyme kinetics has attracted a great attention of scientists from biology, medicine, chemistry, and pharmacy. Mathematical models of enzyme kinetics offer several advances for this deep and analytical understanding due to their in compensable potential in predicting kinetic processes and anticipating appropriate interventions when required. This chapter concerns mathematical modeling analysis and simulation of enzyme kinetics. Experimental data and available knowledge on enzyme mechanics are used in constituting a mathematical model. The models are either in the form of linear or nonlinear ordinary differential equations or partial differential equations. These equations are composed of kinetic parameters such as kinetic rate constants, initial rates, and concentrations of enzymes. The nonlinear nature of enzymatic reactions and a large number of parameters have caused major issues with regard to efficient simulation of those reactions. In this work, an enzymatic system that includes Michaelis-Menten and Ping Pong kinetics is modeled in the form of differential equations. These equations are solved numerically in which the system parameters are estimated. The numerical results are compared with the results from an existing work in literature.


  • mathematical modeling
  • enzyme kinetics
  • chemical kinetics
  • nonlinear reaction-diffusion equation
  • amperometric
  • cyclic voltammetry
  • chronoamperometric

1. Introduction

Enzyme kinetics is a challenging research field nowadays incorporating modern applied mathematics into biotechnology, engineering science, and pharmacy. Moreover, in medical studies, scientists work on human metabolism to improve the capabilities of some metabolites or enzymes in metabolic pathways. In industrial applications, kinetics methods are also widely used to develop certain methods for improving functionality of some molecules in a cell. Many problems in theoretical and experimental biology/chemistry involve the solution of the steady-state reaction diffusion equation with nonlinear chemical kinetics. Such problems also arise in the formulation of substrate and product material balances for enzymes immobilized within particles [1, 2], in the description of substrate transport into microbial cells [3, 4, 5], in membrane transport, in the transfer of oxygen to respiring tissue [6, 7], and in the analysis of any artificial kidney system [8].

To impose the functionality of some molecules in a cell, a mathematical model of such metabolic systems must be constructed and simulated. Most of the dynamical systems can be approximated by various types of differential and integral equations involving finite number of variables and parameters. Thus, the future behavior of the system can be predicted if model kinetics parameters and initial states of the variables are available. In particular, ordinary and partial differential equations (ODEs and PDEs) are popular in modeling of the metabolic pathways or enzyme kinetics.

Releasing enzyme-substrate reactions under single-molecule kinetics was reported by Shlomi et al. [9]. An integral equation method with Michaelis-Menten kinetics to solve nonlinear diffusion problems in spherical coordinates was stated by Tosaka and Miyale [10]. Maalmi et al. [11] reported numerical and semianalytical solutions of nonlinear equations, which covered diffusivity, size, bulk concentration of reactant, binding constant of Michaelis-Menten kinetics, and site reactivity values. Merchant [12] stated the M-M decay reaction terms and the Gray-Scott scheme along with the semianalytical method to nonlinear reaction-diffusion systems. Indira and Rajendran [13] described a homotopy perturbation method to obtain substrate and product concentrations within the enzymatic layers. Removal of substrate from Michaelis-Menten kinetics governed the extravascular partition in which the analytical solution for the steady-state condition was investigated by Bucolo and Tripathi [14]. Dang Do and Greenfield [15] utilized the finite integral transform method to elucidate the problem based on the nonlinear reaction diffusion coupled with the chemical kinetics of a general shape solid. Chapwanya et al. [16] conveyed an epidemiological model with the Michaelis-Menten contact rate formulation to investigate variations in the enzyme kinetics with a simple susceptible infected recovered (SIR) model. Napper [17] proposed the Michaelis-Menten kinetics model to investigate the oxygen transport to heart tissue. Regalbuto et al. [18] presented an analytical methodology for obtaining solutions based on the maximum principle to nonlinear reaction-diffusion boundary value problems.

Rajendran and Saravanakumar [19] discussed mediated bioelectrocatalysis in order to build bioreactors, bio fuel cells, and biosensors.

Due to the difficulties in solving nonlinear differential equations in enzyme kinetics, some recent advanced analytical and numerical simulation techniques are used to solve the problems in chemical kinetics. Thus, in this review, all analytical and numerical works in enzyme kinetics are summarized.


2. Reaction diffusion systems

Reaction diffusion system is a mathematical model based on how the concentration of substances/products is disseminated over space changes under the influence of diffusion and a local chemical reaction. The substances are transformed into each other in local chemical reaction, whereas the substances are spread out over a surface in space in diffusion. Reaction-diffusion (RD) systems arise in many branches of physics, chemistry, biology, ecology, etc. Reviews of the theory and applications of reaction-diffusion systems can be found in books and numerous articles (see, for example [20, 21, 22, 23]). These arise in a large variety of application areas, such as flow in porous media [24], heat conduction in plasma [25], combustion problems [26], liquid evaporation [27], and of more recent interest, image processing [28]. A great effort is being made in the development of the mathematical theory of nonlinear diffusion equations and to obtain exact solutions for special cases. Their significance not only relies on the huge number of their applications but also on the fact that they provide with a rather general class of linear and nonlinear differential operators. In mathematical analysis, it has shown to be a milestone for the development of applied, abstract, and numerical analysis as well as for algebra, geometry, and topology.

3. Nonlinear phenomena

The modern theory of the nonlinear reaction diffusion process is an important field in today’s science. The nonlinear system and coherent structures represent an interdisciplinary area with many nonlinear applications in various fields. Those applications can be divided into six disciplines: chemistry (autocatalytic chemical and enzyme reactions), physics (nonlinear optics and electric circuits, plasmas and states of solid, condensed atomic gases, hydrodynamics, galaxy dynamics and cosmology, fluid dynamics, and celestial mechanics), general relativity, biology (biofuel cell, bioreactor and biosensor, atmosphere and oceans, and animal dispersal), random media, and modern telecommunications. A great variety of phenomena in physics, chemistry, or biology can be described by nonlinear ODE/PDEs and particularly by reaction-diffusion equations. For these reasons, the theory of the analytical solutions of the reaction-diffusion equations is considered.

In reaction diffusion systems, nonlinear phenomena play a crucial role in applied mathematics and chemistry. Exact (closed-form) solution of nonlinear reaction diffusion equations plays an important role in the proper understanding of qualitative features of many phenomena and processes in various areas of natural science. The main result obtained from reaction and diffusion systems is that nonlinear phenomena include diversity of stationary and spatio-temporary dissipative patterns, oscillations, different types of waves, excitability, biostability, etc. But it is difficult for us to obtain the exact solution for these problems. The investigation of exact solution of nonlinear equation is interesting and important. In general, this results in the need to solve linear and nonlinear reaction diffusion equations with complex boundary conditions. The enzyme kinetics in biochemical systems have usually been modeled by differential equations, which are based only on reaction without spatial dependence of the various concentrations. The dimensionless nonlinear reaction diffusion equations are described below:


where Sand Prepresent the dimensionless concentrations of substrate and product, τrepresents the dimensionless time, and Ris the dimensionless radial co-ordinate of the particle. The first term on the right-hand side of the above equation accounts for active species (substrate or product) diffusion, whereas the second term f(R, τ, S, P) and g(R, τ, S, P) represents the homogeneous reaction term (nonlinear term), generally polynomial in the concentrations and time.

4. Common geometries and nonlinear reaction

Most commonly used electrodes/microelectrodes consist of a conducting metal/glassy carbon or semiconducting surface embedded in an insulating wall. When the conducting surface is a rectangle or disc of a few millimeters, this is known as a “planar” electrode. Diffusion to this surface is effectively planar (the effects of the edges are negligible), hence the nonlinear one-dimensional reaction diffusion equation is given by:


Two other electrode geometries where diffusion occurs in only one spatial dimension are the hemispherical and hemicylindrical electrodes. The nonlinear two-dimensional (hemispherical or spherical) reaction diffusion equation is:


and for the latter is:


The hemisphere can be achieved experimentally via a small drop of mercury positioned over a smaller conducting disc. A soft polymer, rubber, or other similar materials are usually employed to fabricate a hemicylinder. The electrodes are usually employed in theoretical studies due to the low dimensionality of the mass-transport equation. Additional terms such as diffusion and nonlinear reaction allow the equation to be solved analytically. Furthermore, the electrodes are not accurately or easily fabricated for practical geometries.

The corresponding nonlinear reaction-diffusion issues in enzyme kinetics are focused on the mathematical resolution. Table 1 shows the response of particular electrodes with special emphasis on earlier theoretical works in the field.

AuthorReferenceExperimental techniqueEnzymatic schemeModeling method
Analytical solutions
G. Rahamathunissa et al.Journal of theoretical and Computational Chemistry, 7(1)(2008)113–138AmperometricS+CKMSCPC'kcP+CC'kE'CDanckwort’s expression
R. Senthamarai et al.Electrochemical Acta 53(2008)3566–3578ChronoamperometricA + eB
G. Rahamathunissa
L. Rajendran
Journal Mathematical Chemistry 44(2008)849–801AmperometricE+SKMESK2E+PVariation iteration method (VIM)
A. Meena et al.Journal Mathematical Chemistry, 48(2010)179–186AmperometricE+SK1K1ESE+PHe’s variation iteration method
A. Eswari, L. RajendranJournal of Electroanalytical Chemistry 641(2010)35–44AmperometricS+E1K2K1E1SKcatP+E2Homotopy perturbation method (HPM)
P. Manimozhi et al.Sensors and Actuators B 147(2010)290–297AmperometricE+Sk1k1ESkcE+P
Variational iteration and homotopy perturbation method (VIM & HPM)
S. Logambal, L. RajendranElectrochemical Acta 55(2010)5230–5238AmperometricA+E2KAB+E1
Homotopy perturbation method (HPM)
A. Meena, L. RajendranJournal of Electroanalytical Chemistry, 6411(2010)50–59Amperometric and PotentiometricE + S↔[ES]→E + PHomotopy perturbation method (HPM)
S. Anitha, L. RajendranJournal of Physical Chemistry 114(2010)7030–7037AmperometricBDEB+SKA+ZABReduction of order method
P. Manimozhi, L. RajendranJournal of Electroanalytical Chemistry 647(2010)87–92AmperometricS+EKαKaES
A. Eswari, L. RajendranJournal of Electroanalytical Chemistry 648(2010)36–46AmperometricS+E1K1K1E1SP+E2Homotopy perturbation method (HPM)
A. Eswari, L. RajendranRussian Journal of Electroanalytical Chemistry 47(2011)195–204Cyclic voltammetryEA+eBCBk1ProductsLaplace Transformation
A. Eswari, L. RajendranRussian Journal of Electroanalytical Chemistry 47(2011)205–212Cyclic voltammetryEA+eBC2B+Bk1ProductsHomotopy perturbation method (HPM)
A. Eswari, L. RajendranJournal of Electroanalytical Chemistry 651(2011)173–184ChronoamperometricO+neRR+ZkO+ProductsHomotopy perturbation method (HPM)
G. Rahamathunissa et al.Journal of Mathematical Chemistry 9(2011)457–474ChronoamperometricS+EkMESk2E+PVIM
S. Logambal, L. RajendranJournal of Membrane Sciences 373(2011)20–28AmperometricEOX+SkMk1ESk2Ered+PEred+O2k3EOX+H2O2Homotopy perturbation method (HPM)
S. Anitha et al.Electrochimica Acta 56(2011)3345–3352AmperometricS+E1KME1SkcatP+E2ABHomotopy perturbation method (HPM)
K. Indra, L. RajendranElectrochimica Acta 56(2011)6411–6419ChronoamperometricS1+O2PPOP2+H2OV1P2+2e+2H+krk0S2E0S2+12O2PPOP2+H2OV2Homotopy perturbation method (HPM)
S. Thiagarajan et al.Journal of Mathematical Chemistry DOI: 10.1007/s10919-011-9854-zChronoamperometricS+MoxkMSMoxkcatP+MredHomotopy perturbation method (HPM)
M. Uma Maheswari, L. RajendranJournal of Mathematical Chemistry DOI: 10.1007/s10919-011-9853-0ChronoamperometricE+SK1k1ESk2E+PHomotopy perturbation method (HPM)
P. Rijiravanich et al.Electroanalytical Chemistry 589(2006)249AmperometricO2+2catechol2oquinone+2H2Ooquinone+2H++2ecatecholTheory and experiment
A. Eswari, L. RajendranJournal of Electroanalytical Chemistry 660(2011)200–208AmperometricO2+2catechol2oquinone+2H2Ooquinone+2H++2ecatecholVIM
G. Varatharajan, L. RajendranApplied Mathematics 2(2011)1140–1147AmperometricS+EK1k1CkcatP+EEk3EiHomotopy perturbation method (HPM)
K. Venugopal et al.Journal of Biomedical Science and Engineering 4(2011)631–641ChronoamperometricO2+2catechol2oquinone+2H2Ooquinone+2H++2ecatecholHomotopy perturbation method (HPM)
K. Indra, L. RajendranJournal of Mathematical Chemistry DOI: 10.1007/s10910-011-9968-3ChronoamperometricAB+CB±eproductsHomotopy perturbation method (HPM)
V. Margret Ponrani, L. RajendranJournal of Mathematical Chemistry DOI: 10.1007/s10910-011-9973-6AmperometricG+EK1k1Xk2k2F+EHomotopy perturbation method (HPM)
S. Sevukaperumal et al.Applied Mathematics 3(2012)373–381ChronoamperometricGlucose+O2Glucoseoxidasegluconicacid+H2O2H2O2CatalaseH2O2+12O2Homotopy analysis method (HPM)
Numerical solution
R. Baronas et al.Biosensors and Bioelectronics 19(2004)915–922AmperometricSPESFinite-difference technique
R. BaronasElectrochimica Acta 240(2017)399–407Amperometric biosensorS+Ek1k1ESP+ESEPNumerical simulation and analytical solution
V. Ašerisa et al.Journal of Electroanalytical Chemistry 685(2012)63–71Amperometricparallel substrates conversionS1E112P1S1+S2E2P2Digital simulation-finite-difference technique
V. Flexer et al.Bioelectrochemistry 74(2008)201–209Cyclic voltammetryS+EOXk1k1ESkcatP+EredNumerical simulation
R. Baronas et al.Chemometrics and Intelligent Laboratory Systems 126(2013)108–116AmperometricE+Sik1iESik2iE+Pi,i=1,,kNumerical
R. BaronasNonlinear Analysis: Modeling and Control 9(3)(2004)203–218AmperometricSEPDigital simulation-finite-difference technique.
R. Baronas et al.Sensors 12(2012)9146–9160AmperometricEOX+Sk1Ered+PEredk2EOX+neeFinite-difference
R. Baronas et al.J. Mathematical Chemistry 32 (2)(2002)225–237AmperometricSEPNumerical simulation
R. Baronas et al.Mathematical Modeling of Biosensors, Springer Series on chemical sensors and biosensors (2009)AmperometricAll enzyme reactionsAnalytical and numerical methods
L. RajendranBiosensor: Modeling and Simulation of Diffusion-Limited Process, Chemical Sensors: Simulation and Modeling, GhenadiiKorotcenkov (Ed.), Electrochemical Sensors, Vol. 5, Momentum Press, LLC, New York (2013)AmperometricAll enzyme reactionsAnalytical, HPM&HAM, VIM,ADM, etc.

Table 1.

Contributions to the theoretical modeling of enzymatic electrodes.

Example 1: Michaelis-Menten kinetics and microcylinder electrodes

The model is written for an enzyme reaction to generate an electro-active product (e.g., hydrogen peroxide from an oxidase enzyme) that reacts at an immobilization matrix, which is metallically conducting sites/particles. The reaction within the film under the Michaelis-Menten kinetics may be written as follows:


The consumption rate of Sis given by k1cScE − k−1cES, where cidenotes the concentration of species i. The rate is equivalent to (kcat/KM) cScE, where KMis the Michaelis constant, defined as KM = (k−1 + kcat)/k1. The consumption rate of Sin the film is compensated by diffusion. If the solution is stirred uniformly, so that Sis constantly supplied to the film, the mass balance for Scan be written in cylindrical coordinates:


where cSis the concentration profile of substrate, cEis the concentration profile of enzyme, DSis its diffusion coefficient, and KMis the Michaelis constant. The rate of consumption will be v(r) = k cH, where kis the rate constant for the hydrogen peroxide reaction and cHis the peroxide concentration. Then, the equation of continuum for hydrogen peroxide is generally expressed in the steady-state by


At the electrode surface (r0) and at the film surface (r1), the boundary conditions are [29]:


where cSis the bulk concentration of Sscaled by the partition coefficient of the film. The current is provided by the consumption rate at each site. Thus, the total current at an electrode of length Lis expressed by [29]


The analytical results of the problem are discussed by Eswari and Rajendran [30].

Example 2: enzyme catalysis reaction

The reactions without spatial dependence on various concentrations have modeled the enzyme kinetics in biochemical systems. Nonlinear systems of ordinary differential equations are solely based on that. Michaelis and Menten were pioneers in explaining the enzyme reaction model. In addition, they also reported the free enzyme binding to the reactant, which produced an enzyme-reactant complex. Eq. (11) illustrates the Michaelis-Menten kinetics, in which the enzyme-substrate complex is formed after the enzyme is combined with the substrate.


As can be seen from Eq. (11), the product Pis released by the binding of substrate Swith enzyme E. The product released is not reversible; however, the substrate binding is reversible. The reactants’ concentrations in Eq. (11) are represented by the following letters:


The law of mass action leads to the system of following nonlinear reaction equations [31],


where k1 is the forward rate of ES complex formation and k−1 is the backward rate constant. The above problem is discussed theoretically by Meena et al. [32].

Example 3: Michaelis-Menten mechanism for co-substrate and substrate

Figure 1 illustrates Michaelis-Menten reaction kinetics scheme for co-substrate and substrate. Limoges et al. [33] reported for a redox enzymatic homogenous system along with one-dimensional mass transport equation a concise discussion and derivation.

Figure 1.

Reaction scheme for substrate and co-substrate.

When the enzyme is being solubilized, the electrochemical signal that is produced during the reaction is governed by the following set of nonlinear partial differential equations.


where DP , DSare the diffusion coefficients of co-substrate and substrate, respectively; Q , Sare the concentrations of co-substrate and substrate, respectively; xis the distance from the electrode surface; CS0is the bulk concentration of substrate; CE0is the total concentration of enzyme; k1,  k2 , 2, and  k2 are the reaction rate constants; and tis the time. The initial and boundary conditions for Eqs. (14) and (15) are given by:


The analytical expressions corresponding to the concentration of co-substrate for steady and nonsteady state conditions have been obtained by solving the above nonlinear equation using a new approach to homotopy perturbation method (HPM). Analytical expressions of the plateau current are also presented for steady and nonsteady state conditions:


where Eis the electrode potential, EPQ0is the standard potential of the P/Q couple, Fis the Faraday constant, and Sis the surface area of the electrode. The above problem is discussed theoretically by Rasi et al. [34].

5. Analytical solutions

To study many of the physical phenomena, the exact solutions of nonlinear partial or ordinary differential equations play an important role. In order to understand the mechanism of complicated dynamical processes and physical phenomena modeled by nonlinear differential equations, the existence of approximate analytical and exact solutions is very important. In addition, nonlinear differential equations can also assist to investigate the stability of these solutions as well as checking the simulation analysis. Nonlinear partial differential equations govern a significant variety of phenomena including physical, chemical, and biological. The development of techniques aimed at exact solutions of nonlinear differential equations with nonsteady and steady state [35] has been one of the most exciting advances of nonlinear science and theoretical physics/chemistry. An important role in nonlinear science is played by exact solutions of differential equations. Furthermore, this can be especially observed in nonlinear physical chemistry science. This can be attributed to the provision of physical information as well as more insight into the physical aspects of the problem, which could lead to further applications. Over the past few decades, different methods have been reported to solve analytical solutions such as Tanh-sech [36], extended tanh [37], Jacobi elliptic function expansion [39], hyperbolic function [38], F-expansion [40], and the First integral [41]. To solve different types of nonlinear systems of PDEs, the sine-cosine method [42] has been employed. A variety of powerful analytical methods such as homotopy perturbation method [43, 44, 45], homotopy analysis method [46, 47], Adomian decomposition method [48, 49], wavelet transform method [50], etc. are applied to solve the nonlinear problems (e.g., Eqs. (8) and (13)(15)) in chemical kinetics [51].


6. Numerical solutions

Many differential equations cannot be solved analytically. For practical purpose, however, such as in physical engineering sciences, a numerical approximation to the solution is often sufficient. The numerical method is mainly to solve complex problem physically or geometrically. It is also used to validate the experimental results. Some of the nonlinear equations in chemical kinetics were solved using numerical methods [52, 53, 54, 55, 56].

7. Summary

Most mathematical models of enzyme kinetics are based on reaction diffusion equations or rate equations containing nonlinear terms related to the kinetics of the enzyme reaction. Powerful and accurate analytical (HPM, HAM, ADM, etc.) and numerical mathematical methods have been employed for their resolution under steady and nonsteady state conditions. The theoretical results provide very useful insight into the effects on the performance of the thickness and structure of the enzymatic film, the loading of the different species, the diffusivity of the mediator, etc. Also, the theoretical modeling and simulation of these systems enable us to characterize the enzymatic reactions (i.e., rate constant, turnover rate, and Michaelis-Menten constants).

In spite of the above-mentioned benefits, there are only limited theoretical studies addressing kinetics of enzyme reaction and most of them include a number of simplifying assumptions mainly related to the mass and charge transport inside and outside the biocatalyst film, the enzymatic kinetic scheme, and the electrode morphology. Experimental validation of proposed models is even more seldom. Therefore, more effort in the future research is needed in this direction in order to develop more detailed models and accurate simulations that can assist the rational development and optimization of enzyme electrodes.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Lakshmanan Rajendran, Mohan Chitra Devi, Carlos Fernandez and Qiuming Peng (December 20th 2017). Mathematical Modeling and Simulation of Nonlinear Process in Enzyme Kinetics, Advanced Chemical Kinetics, Muhammad Akhyar Farrukh, IntechOpen, DOI: 10.5772/intechopen.70914. Available from:

chapter statistics

1169total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Autoignition and Chemical-Kinetic Mechanisms of Homogeneous Charge Compression Ignition Combustion for the Fuels with Various Autoignition Reactivity

By Dongwon Jung

Related Book

First chapter

Electronic (Absorption) Spectra of 3d Transition Metal Complexes

By S. Lakshmi Reddy, Tamio Endo and G. Siva Reddy

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us