Research Article  Open Access
Mathematical Modeling of Multienzyme Biosensor System
Abstract
A mathematical model of hybrid inhibitor biosensor system is discussed. This model consists of five nonlinear partial differential equations for bisubstrate sensitive amperometric system. Simple and closed form of analytical expressions for concentration of glucose6phosphate (substrate), potassium dihydrogen phosphate (inhibitor), oxygen (cosubstrate), glucose (product 1), and hydrogen peroxide (product 3) is obtained in terms of rate constant using modified Adomian decomposition method (MADM). In this study, behavior of biokinetic parameters is analyzed using this theoretical result. The obtained analytical results (concentrations) are compared with the numerical results and are found to be in satisfactory agreement.
1. Introduction
A biosensor is an analytical device, used for the detection of an analyte or substrate that combines a biological component with a physicochemical detector [1]. Biosensors are potentially useful as they detect components quickly and have been active in the research area for certain years. The development of these biosensing systems depends on a quantitative measurement of the enzyme activity before and after exposure to a target analyte. Inhibitionbased biosensors have been the topic of several recent reviews [2, 3]. Biosensors based on the principle of enzyme inhibition have been applied for a wide range of significant analytes. Some recent approaches and developments in enzyme inhibitionbased biosensors for food safety and environmental monitoring have been discussed by Amine et al. [4]. The application of wholecellbased biosensors to different environmental media, such as water, soil, and atmospheric monitoring has been reviewed by Gu et al. [5]. An algal biosensor for toxicity assessment of estuarine waters was tested at different concentrations [6]. On biosensor technology, initial studies for the construction of a diffusive sampling device was described for monitoring personal exposure to formic acid [7]. The determination of benzene in workplace air samples using a bacterial biosensor based on flow injection analysis (FIA) has been developed by Lanyon et al. [8]. Kok and Hasirci [9] reported that the total anticholinesterase activity of binary pesticide mixtures was lower than the sum of the individual inhibition values. Mazzei et al. [10] immobilized acid phosphatase and glucoseoxidase on a separate dialytic membrane. The bienzyme modified membrane was placed on platinum electrode using an “O” ring. On enzymatic biosensor, the steady state analytical expression of concentration of substrate is derived using variational iteration method [11]. A monolayer potentiometric biosensor response for a polymermodified electrode system has been discussed by Rahamathunissa et al. [12]. The biosensors can be used for biomedical applications [13] and electric circuits [14]. Some nonlinear problems in different types of biosensors are solved using various numerical methods [15] and analytical methods [16].
The multienzyme biosensor system was used for the detection of phosphorous and fluoride ions. Rangelova et al. [17] investigated the influence of concentration and reaction rate constant of inhibitor with two enzymes in the dynamic mode. In this work, we study the steady state performance of multienzyme inhibitor system performing glucose6phosphate hydrolysis at different conditions. The analytical expression for concentrations and current are derived using modified Adomian decomposition method. The effect of internal mass transfer parameters and kinetic parameters of inhibitor is discussed.
2. Mathematical Formulation of the Problem
In enzymatic reactions, the molecules at the beginning of the process, called substrates, are converted into different molecules, called products. Almost all chemical reactions in a biological cell need enzymes in order to occur at rates sufficient for life. The working mechanism of hybrid biosensor for phosphate is based on the following sequential biochemical reactions [8]: The biochemical reaction, which is biocatalyzed by acid phosphatase (AP), involves the hydrolysis of glucose6phosphate (C_{6}H_{13}O_{9}P) with the liberation of glucose (C_{6}H_{12}O_{6}) and inorganic phosphate (). In the second reaction with glucose oxidase (GOD), the oxygen present oxidizes the glucose into glucono deltalactone (C_{6}H_{10}O_{6}) and hydrogen peroxide (H_{2}O_{2}). The amount of hydrogen peroxide being produced is measured electrochemically. On the basis of the biochemical reactions, the following successive enzyme reactions with competitive inhibition have been proposed: where , and are the concentrations of the glucose6phosphate (substrate), potassium dihydrogen phosphate (inhibitor), oxygen (cosubstrate), glucose (product 1), and hydrogen peroxide (product 3) respectively. (acid phosphatase) is the first enzyme and (glucose oxidase) is the second enzyme. and represent the inorganic phosphate and glucono deltalactone. Figure 1 presents the overall reaction scheme for a multienzyme inhibitor system. We assume that diffusion is onedimensional in space and is described with second Fick’s law. The concentration profiles of five reagents of the system are described by the following nonlinear equations: where , and are diffusion coefficients for substrate, cosubstrate, product 1, and product 3. denote the reaction constant for concentration profiles (), respectively. is the reaction velocity of acid phosphatase and is the reaction velocity of glucose oxidase. The initial conditions are Equations (4) are subjected to the following boundary conditions: The output current is defined as where is the number of electrons, is the Faraday’s constant, and is the area of the electrode surface.
3. Governing Equations in Dimensionless Form
Equations (4)–(7) can be reduced to the corresponding dimensionless forms by introducing the following dimensionless parameters: Consequently, (4)–(7) can be written in dimensionless form for steady state as follows: The boundary conditions are Here we can assume that is very small () at .
The dimensionless current is as follows:
4. Analytical Expressions of the Concentrations and Current Using Modified Adomian Decomposition Method
In order to solve the boundary value problem, (9)–(13), we used the modified Adomian decomposition method [18, 19]. The basic idea of this method is described in Appendix A [20]. Detailed derivations of the dimensionless concentrations , and of glucose6phosphate, potassium dihydrogen phosphate, oxygen, glucose, and hydrogen peroxide are described in Appendix B. As a result, we have obtained where Equations (16)–(20) represent the analytical expression of the concentration profiles of five reagents for all values of the dimensionless parameters and . Also, these expressions satisfy the boundary conditions (14). The current response is given by
5. Numerical Simulation
The nonlinear differential equations (9)–(13) with boundary conditions (14) have been solved numerically using SCILAB/MATLAB software. A respective script pdex4 is provided in Algorithm 1. The numerical solution is compared with our analytical results in Figures 2–6. The comparison reveals that the analytical dimensionless concentrations and numerical results are in satisfactory agreement. The dimensionless concentration of substrate and product 3 obtained using (16) and (20) are reported in Tables 1 and 2. In these tables, our result is compared with simulation results. The average relative error between our analytical result for the concentration of substrate and simulation result are 1.487% when and 0.00% when and (Table 1) whereas the average relative error in the case for concentration of product 3 is 0.00% for all values of (Table 2).
 
NS: numerical simulation. 
 
NS: numerical simulation. 

(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
6. Results and Discussion
The normalized nonlinear differential equations were solved using modified Adomian decomposition method. Equations (16)–(20) represent the analytical expressions for the dimensionless concentrations of glucose6phosphate , potassium dihydrogen phosphate , oxygen , glucose , and hydrogen peroxide . Figure 2 demonstrates quantitatively the effects of dimensionless substrate inhibition modulus, dimensionless substrate Thiele modulus, and dimensionless inhibitor inhibition modulus on the dimensionless concentration profile of substrate as a function of dimensionless distance using (16). It shows that the value of approaches unity when for all values of , and . It is evident that the substrate concentration increases as a result of increasing dimensionless inhibition modulus and dimensionless substrate modulus . Further, the substrate concentration decreases when dimensionless Thiele modulus or thickness of the membrane increases.
Figure 3 represents the dimensionless concentration profile of the inhibitor versus dimensionless distance using (17). It is clear that for the increasing values of dimensionless inhibitor Thiele modulus and dimensionless substrate inhibition modulus , the concentration profile deviates from the linear form. Also, the inhibitor concentration increases when dimensionless inhibitor inhibition modulus increases (see Figure 3(f)). The dimensionless concentration profile of cosubstrate against dimensionless distance is plotted in Figure 4. From these Figures, we conclude that the value of decreases when the dimensionless cosubstrate Thiele modulus increases. Furthermore, the cosubstrate profile increases when dimensionless MichaelisMenten constant of glucose and oxygen increases (see Figures 4(e) and 4(f)).
From Figure 5, it is clear that the dimensionless concentration profile of product 1 increases when product 1 dimensionless Thiele modulus or thickness of the membrane increases. It also decreases with decreasing , approaching zero at . It shows that the value of decreases as a result of increasing dimensionless inhibition modulus and dimensionless substrate inhibition modulus . Figure 6 illustrates the dimensionless concentration profile of the product 3 versus dimensionless distance using (20). From these Figures, it is inferred that the product 3 concentration is in the form of parabola with as the axis of parabola. The value of the concentration of product increases when dimensionless MichaelisMenten constants and decrease. It also shows that the value of decreases when product 3 dimensionless Thiele modulus decreases.
Figures 7(a)7(b) show the effect of dimensionless MichaelisMenten constant on the dimensionless current as a function of product 3 dimensionless Thiele modulus with . The dimensionless current approaches zero and it is independent of at low dimensionless Thiele modulus. Also, the dimensionless current increases with decreasing at high dimensionless Thiele modulus. Figures 7(c)7(d) represent the effect of Thiele modulus on the dimensionless current as a function of product 3 dimensionless MichaelisMenten constant with . The dimensionless current decreases with decreasing dimensionless Thiele modulus . From the Figures 7(e)7(f), it is inferred that the dimensionless current increases when dimensionless MichaelisMenten constant of glucose decreases for some fixed values of .
(a)
(b)
(c)
(d)
(e)
(f)
7. Conclusions
The system of nonlinear, steadystate reaction diffusion equations in multienzyme biosensor has been solved analytically using modified Adomian decomposition method. The accuracy of the approximate analytical solutions has been verified by comparison with numerical solutions. Satisfactory agreement is noted. This method is very effective, simple, and time saving. It also gives the better results. This method can be easily extensive to all kinds of nonlinear differential equations. Moreover, the behavior of inhibitor starting concentration over biosensor can be analyzed through the current expression.
Appendices
A. Basic Concepts of the Modified Adomian Decomposition Method
The modified Adomian decomposition method consists of decomposing the nonlinear differential equation into two components where and are the linear and the nonlinear parts of respectively. The operator is assumed to be an invertible operator. Solving for leads to Applying the inverse operator to both sides of (A.3) yields where is the function that satisfies the condition . Now suppose that the solution can be represented as an infinite series of the form The modified Adomian decomposition method assumes that the nonlinear term can be written as an infinite series in terms of the Adomian polynomials : where the Adomian polynomials of are evaluated using the formula where is a hypothetical parameter.
Substituting (A.5) and (A.6) in (A.4) gives By equating the terms in the linear system of (A.8) one obtains the recurrence formula However, in practice all terms of the series (A.6) cannot be determined, and the solution is approximated by the truncated series .
B. Analytical Solutions for the Concentration Profiles of Five Reagents
The analytical solution is necessary to predict the concentration profiles of glucose6phosphate, potassium dihydrogen phosphate, oxygen, glucose, and hydrogen peroxide solution of (9)–(13). Modified Adomian decomposition method is employed [21, 22] to construct approximate nonperturbative solutions of five reagents described by (9)–(13). In this appendix the expressions are derived for the concentrations of glucose6phosphate and potassium dihydrogen phosphate. The operator form is where . Applying the inverse operator to both sides of (B.1) yields, in agreement with (A.4), (A.6), and (A.8) where , and are the integration constant and the Adomian polynomial coefficients can be obtained using (A.7). By equating the terms of (B.3) and using the boundary conditions (14), we get The first Adomian polynomial coefficients are From the above equations we obtain , as follows: where Adding (B.4) and (B.6) we obtain (16) and (17) in the text. Similarly we can find the analytical solution for (11)–(13).
C. Matlab Code
The Matlab program to find the numerical solution of (9)–(13) is as shown in Algorithm 1.
Notations
Concentration of glucose6phosphate (mol/cm^{3})  
:  Concentration of potassium dihydrogen phosphate (mol/cm^{3}) 
:  Concentration of oxygen (mol/cm^{3}) 
:  Concentration of glucose (mol/cm^{3}) 
:  Concentration of inorganic phosphate, (mol/cm^{3}) 
:  Concentration of hydrogen peroxide (mol/cm^{3}) 
:  Concentration of glucono deltalactone, C_{6}H_{10}O_{6} (mol/cm^{3}) 
:  Acid phosphatase (AP) 
:  Glucose oxidase (GOD) 
:  Reaction velocity of acid phosphatase (mol/s cm^{3}) 
:  Reaction velocity of glucose oxidase (mol/s cm^{3}) 
:  Reaction constant for glucose6phosphate (mol/cm^{3}) 
:  Reaction constant for potassium dihydrogen phosphate (mol/cm^{3}) 
:  Reaction constant for oxygen (mol/cm^{3}) 
:  Reaction constant for glucose (mol/cm^{3}) 
:  Diffusion coefficient of glucose6phosphate (cm^{2}/s) 
:  Diffusion coefficient of potassium dihydrogen phosphate (cm^{2}/s) 
:  Diffusion coefficient of oxygen (cm^{2}/s) 
:  Diffusion coefficient of glucose (cm^{2}/s) 
:  Diffusion coefficient of hydrogen peroxide (cm^{2}/s) 
:  Distance (cm) 
:  Concentration of glucose6phosphate in bulk solution (mol/cm^{3}) 
:  Concentration of potassium dihydrogen phosphate in bulk solution (mol/cm^{3}) 
:  Concentration of oxygen in bulk solution (mol/cm^{3}) 
:  Thickness of the membrane (cm) 
:  Output current at the electrode surface (ampere). 
:  Dimensionless concentration of glucose6phosphate 
:  Dimensionless concentration of potassium dihydrogen phosphate 
:  Dimensionless concentration of oxygen 
:  Dimensionless concentration of glucose 
:  Dimensionless concentration of hydrogen peroxide 
:  Dimensionless distance 
:  Dimensionless current. 
:  Inhibitor inhibition modulus (dimensionless) 
:  Substrate inhibition modulus (dimensionless) 
:  MichaelisMenten constants (dimensionless) 
:  Substrate Thiele modulus (dimensionless) 
:  Inhibitor Thiele modulus (dimensionless) 
:  Cosubstrate Thiele modulus (dimensionless) 
:  Product 1 Thiele modulus (dimensionless) 
:  Product 3 Thiele modulus (dimensionless). 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the Department of Science and Technology (DST) (no. SB/SI/PC50/2012), New Delhi, India. The authors are thankful to The Principal, The Madura College, Madurai, and The Secretary, The Madura College Board, for their encouragement. The authors are grateful to the reviewer’s valuable comments that improved the paper.
References
 A. Sadana and N. Sadana, Handbook of Biosensors and Biosensor Kinetics, Elsevier, New York, NY, USA, 2010.
 S. V. Dzyadevych, A. P. Soldatkin, Y. I. Korpan et al., “Biosensors based on enzyme fieldeffect transistors for determination of some substrates and inhibitors,” Analytical and Bioanalytical Chemistry, vol. 377, no. 3, pp. 496–506, 2003. View at: Publisher Site  Google Scholar
 M. D. L. de Castro and M. C. Herrera, “Enzyme inhibitionbased biosensors and biosensing systems: questionable analytical devices,” Biosensors and Bioelectronics, vol. 18, no. 23, pp. 279–294, 2003. View at: Publisher Site  Google Scholar
 A. Amine, H. Mohammadi, I. Bourais, and G. Palleschi, “Enzyme inhibitionbased biosensors for food safety and environmental monitoring,” Biosensors and Bioelectronics, vol. 21, no. 8, pp. 1405–1423, 2006. View at: Publisher Site  Google Scholar
 M. B. Gu, R. J. Mitchell, and B. C. Kim, “Wholecellbased biosensors for environmental biomonitoring and application,” Advances in Biochemical Engineering/Biotechnology, vol. 87, pp. 269–305, 2004. View at: Publisher Site  Google Scholar
 L. Campanella, F. Cubadda, M. P. Sammartino, and A. Saoncella, “An algal biosensor for the monitoring of water toxicity in estuarine environments,” Water Research, vol. 35, no. 1, pp. 69–76, 2000. View at: Publisher Site  Google Scholar
 K. J. M. Sandström, J. Newman, A.L. Sunesson, J.O. Levin, and A. P. F. Turner, “Amperometric biosensor for formic acid in air,” Sensors and Actuators B: Chemical, vol. 70, no. 1–3, pp. 182–187, 2000. View at: Publisher Site  Google Scholar
 Y. H. Lanyon, G. Marrazza, I. E. Tothill, and M. Mascini, “Benzene analysis in workplace air using an FIAbased bacterial biosensor,” Biosensors and Bioelectronics, vol. 20, no. 10, pp. 2089–2096, 2005. View at: Publisher Site  Google Scholar
 F. N. Kok and V. Hasirci, “Determination of binary pesticide mixtures by an acetylcholinesterase choline oxidase biosensor,” Biosensors and Bioelectronics, vol. 19, no. 7, pp. 661–665, 2004. View at: Publisher Site  Google Scholar
 F. Mazzei, F. Botrè, and C. Botrè, “Acid phosphatase/glucose oxidasebased biosensors for the determination of pesticides,” Analytica Chimica Acta, vol. 336, no. 1–3, pp. 67–75, 1996. View at: Publisher Site  Google Scholar
 G. Rahamathunissa and L. Rajendran, “Application of He's variational iteration method in nonlinear boundary value problems in enzymesubstrate reaction diffusion processes: part 1. The steadystate amperometric response,” Journal of Mathematical Chemistry, vol. 44, no. 3, pp. 849–861, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 G. Rahamathunissa, P. Manisankar, L. Rajendran, and K. Venugopal, “Modeling of nonlinear boundary value problems in enzymecatalyzed reaction diffusion processes,” Journal of Mathematical Chemistry, vol. 49, no. 2, pp. 457–474, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. Zhang, H. Ju, and J. Wang, Eds., Electrochemical Sensors, Biosensors and their Biomedical Applications, Academic Press, New York, NY, USA, 2007.
 J.Y. Yoon, Introduction to Biosensors: From Electric Circuits to Immunosensors, Springer, 2013.
 R. Baronas, F. Ivanauskas, and J. Kulys, Mathematical Modeling of Biosensors: An Introduction for Chemists and Mathematicians, Springer Series on Chemical Sensors and Biosensors, Springer, Dordrecht, The Netherlands, 2010. View at: Publisher Site  MathSciNet
 L. Rajendran, Chemical Sensors: Simulation and Modeling Volume 5: Electrochemical Sensors, Edited by: G. Korotcenkov, Momentum press, New York, NY, USA, 2013.
 V. Rangelova, A. Pandelova, and N. Stoiyanov, “Inhibitor multienzyme biosensor system in dynamic modephosphate measurement,” Annals of Faculty Engineering HunedoaraInternational Journal of Engineering, vol. 9, no. 2, pp. 83–86, 2011. View at: Google Scholar
 G. Adomian, “Convergent series solution of nonlinear equations,” Journal of Computational and Applied Mathematics, vol. 11, no. 2, pp. 225–230, 1984. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Mathematical and Computer Modelling, vol. 13, no. 7, pp. 17–43, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Q. Hasan and L. M. Zhu, “Modified adomian decomposition method for singular initial value problems in the secondorder ordinary differential equations,” Surveys in Mathematics and Its Applications, vol. 3, pp. 183–193, 2008. View at: Google Scholar  MathSciNet
 M. M. Hosseini, “Adomian decomposition method with Chebyshev polynomials,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1685–1693, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A.M. Wazwaz, “A new algorithm for calculating Adomian polynomials for nonlinear operators,” Applied Mathematics and Computation, vol. 111, no. 1, pp. 53–69, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2014 SP. Ganesan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.