Research Article  Open Access
Theoretical Analysis of an Amperometric Biosensor Based on Parallel Substrates Conversion
Abstract
The behaviour of an amperometric biosensor based on parallel substrates conversion for steadystate condition has been discussed. This analysis contains a nonlinear term related to enzyme kinetics. Simple and closed form of analytical expressions of concentrations and of biosensor current is derived. This model was originally reported by Vytautas Aseris and his team (2012). Concentrations of substrate and product are expressed in terms of single dimensionless parameter. A new approach to Homotopy perturbation method (HPM) is employed to solve the system of nonlinear reaction diffusion equations. Furthermore, in this work, the numerical solution of the problem is also reported using Matlab program. The analytical results are compared with the numerical results. The analytical result provided is reliable and efficient to understand the behavior of the system.
1. Introduction
A biosensor is an analytical device, used for the detection of an analyte, which combines a biological component with a physiochemical detector [1]. Analyzed substrate is biochemically converted to a product when the biosensor operation is being undertaken. In most of the cases, biosensor response is directly proportional to the concentration of the reaction product [2]. The amperometric biosensors have proved to be reliable in various systems with applications in the field of medicine, food technology, and the environmental industry [3, 4]. The understanding of the kinetic peculiarities of biosensors is of crucial importance for their design. The mathematical modeling is rather widely used to improve the efficiency of the biosensors design and to optimize their configuration [5–8]. Since 1970s, various mathematical models of biosensors have been developed and used successfully to study and optimize analytical characteristics of biosensors [9–14]. Mathematical modeling of twoenzyme biosensors has been started in 1980s with the modeling of an amperometric monolayer enzyme electrode with two coimmobilized enzymes [15, 16]. Later, nonlinear mathematical models have been developed for amperometric twoenzyme biosensors with different enzymes [17, 18].
The numerical method of solving partial differential equations is to make calculation at all intervals of substrates concentration and at different diffusion and enzymatic reaction rates. The nonstationary diffusion equations [19], containing a nonlinear term related to the enzymatic reaction, are carried out using the implicit difference scheme [20]. In recent years, analytical solutions are reported for various types of biosensors [21–27].
To the best of our knowledge, there is no simple analytical expression for steadystate concentrations and current of amperometricbased biosensor utilizing parallel substrates conversion for all values of parameters. In this paper, we present the simple and closed form of the analytical expressions for the concentrations of substrate and of product using new Homotopy perturbation method (HPM) [28]. The current density was determined corresponding to all possible values of the parameters.
2. Mathematical Formulation of the Problem and Analysis
The reaction scheme employing two parallel substrates ( and ), two enzymes ( and ), and two reaction products ( and ) in biosensor may be as follows:
First phase of the biosensor action contains hydrogen peroxide () only. The biosensor response reaches the steadystate at the end of the first phase. The second phase of the biosensor action begins at the second substrate () poured into the solution. The relative changes between the responses in both phases are measured as the final response of the biosensor. The schematic diagram of the system is provided in Figure 1. The above reaction in the enzymatic layer with the onedimensionalinspace diffusion, described by Fick’s law, leads to the following nonlinear transient reactiondiffusion equations [29]: where denotes space, , are the molar concentration of the substrates, and , are the molar concentration of products. , , , and are the corresponding diffusion coefficients and , , and are the reaction rate constants. Graphical representation of the products , and substrates , , with steadystate boundary conditions employed in the simulated multilayered biosensor is provided in Figure 2. The mass transport by diffusion of the substrates and products takes place outside the enzyme layer. The diffusion layer is treated as the Nernst diffusion layer. The boundary conditions for the above system of nonlinear equation are as follows: where and denote the thickness of the enzyme and diffusion layer. For the steadystate condition, the following sets of dimensionless variables are introduced: where . The governing nonlinear reaction diffusion equations (2) are expressed in the following steadystate nondimensional form: where and represent the dimensionless substrate concentrations, and represent the dimensionless product concentrations, and represents the dimensionless distance. The reaction/diffusion parameters of substrates and products are considered to be , , , , and . The boundary conditions are represented as follows: The dimensionless current density is given by
3. A New Approach to Homotopy Perturbation Method and Its Advantages
The Homotopy perturbation method is first introduced to solve nonlinear wave equations. The advantage of the method is that it does not need a small parameter in the system, leading to wide application in nonlinear wave equations [30]. Recently, many authors have used HPM to various problems and reported the efficiency of the HPM for handling nonlinear engineering problems [31–34]. This method is a combination of Homotopy in topology and classic perturbation techniques. He used the HPM to solve the Light hill equation [35]; this wide variety of applications shows the power of the HPM to solve functional equations. The HPM has uniqueness in its applicability, accuracy, and efficiency. Recently, a new approach to HPM is introduced to solve the nonlinear problem, in which one will get better simple approximate solution in the zeroth iteration [28].
4. Analytical Expression for the Concentration of Substrate and Product Using New Homotopy Perturbation Method
In this paper, a new approach to Homotopy perturbation method is used (Appendix A) to solve the nonlinear differential equations (5)–(8). Using this method, the analytical expression of the two substrate concentrations ( and ) and two product concentrations ( and ) can be obtained as follows: where Equations (11)–(14) satisfy the boundary conditions (9). These equations represent the simple and closed form of analytical expressions of concentrations of substrate and product for all possible values of the parameters. Using (6) and (8) with boundary conditions (9), the following relation is also obtained: The analytical expression of the current using (10) is as follows: When , the above expression can be simplified as .
4.1. Limiting Case
The mathematical as well as the corresponding computational models of the biosensor were evaluated. Assume that zero concentration of the second substrate () leads to reducing the governing equations (5)–(8) to the following equations: The steadystate solution of the above equations for the boundary condition (9) is given by The dimensionless steadystate current using (10) is given by Equation (20) is identical to (17), when (i.e., ) in (17).
5. Comparison of Analytical Expressions with Experimental Work and Numerical Simulation
To show the efficiency, the steadystate result (11)–(14) of the reaction diffusion equations is compared with the numerical solution of (5)–(8). The numerical solutions are found using pdepe (finite element method) in Matlab [36]. This Matlab program [37] is also given in Appendix B. The comparison is shown in Figures 3–6 and it gives satisfactory agreement.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(f)
6. Discussion
Equations (11)–(14) are the new simple and closed approximate analytical expressions for the dimensionless substrate concentrations ( and ) and product concentrations ( and ) in terms of or . The analytical results (11)–(14) are compared with numerical solution, which is provided in Figures 3–7. Figure 3 shows the plots of all the concentration versus dimensionless distance .
(a)
(b)
(c)
(d)
The kinetic response of a biosensor depends on the concentration species. The concentration depends on the following factors: , , , , and . The diffusion parameters and represent the ratio of the characteristic time of the enzymatic reaction to that of concentration diffusion. These parameters can be varied by changing the thickness of the enzyme or diffusional layer. These parameters describe the relative importance of diffusion and reaction in both layers. When or is small, the kinetics is dominant resistance; the uptake of concentration species is kinetically controlled. Under these conditions, the concentration species across both the layers is essentially uniform. When the diffusion parameter or is large, diffusion limitations are the principal determining factor.
Figures 4 and 5 show the dimensionless concentration of substrates ( and ) versus dimensionless distance . Figures 3 and 4(a) and 4(b) are plotted for various values of diffusion parameters and for some fixed values of and . From these figures, it is clear that the concentration of substrates and decreases when the value of and increases, respectively. The concentration is approximately equal to one when and . Figures 4(c) and 5(c) are plotted for various values of dimensionless parameter for some fixed values of other existing parameters. From these figures, it is clear that concentration increases for large values of . Figures 4(d) and 5(d) are plotted for various values of and , respectively. It is evident that the concentration is uniform, when (or ) and (or ).
Figures 6(a)–6(f) show the dimensionless product concentration profile versus dimensionless distance using (13). Figures 6(a) and 6(b) are plotted for all various values of for some fixed values of , , and . Thus, it is concluded that there is a simultaneous increase in the values of the concentration as well as in . Figure 6(c) is plotted for various values of and the fixed value of , , , and . It shows that the concentration is always minimum for all values of . Concentration versus distance is plotted in Figures 6(d) and 6(e) for various values of and for some fixed values of other existing parameters. Thus, it is concluded that the concentration is always minimum and also there is a simultaneous decrease in the values of the concentration as the value of increases. Figure 6(f) is plotted for various values of , and from these graphs, it is inferred that concentration is minimum and uniform when .
Figures 7(a)–7(d) are plotted for the dimensionless product concentration as ordinate versus dimensionless space as abscissa using (14). From Figures 7(a) and 7(b), it is evident that the concentration of increases as the diffusivity parameter increases. The concentration is approximately equal to zero, that is, () when for all possible values of the parameters , and . From Figure 7(c), it is inferred that concentration decreases, when the value of increases and when . Figure 7(d) is plotted for various values of ; it is evident that the concentration increases as the value of increases. The concentration increases, when or .
Figures 8(a)–8(e) are plotted using (17), which shows the dynamics of normalized steadystate current versus dimensionless parameters , , , , and , respectively. From Figure 8(a), it is noted that, for small values of and for any values of and , the value of the current is significantly high. The current grows notably faster at higher values of but in the case of and , the current does not increase. Figure 8(b) is plotted for various values of and for some fixed values of and . It is observed that the current increases rapidly and reaches the steadystate when . From Figures 8(c) and 8(d), it is evident that the current is uniformly decreasing for large values of , that is, (), and constant elsewhere. Figure 8(e) is plotted for various values of the reaction/diffusion parameter as abscissa and normalized steadystate current as ordinate. It is observed that the current increases as the value of increases. Also current decreases for all values of the reaction/diffusion parameter and reaches zero when .
(a)
(b)
(c)
(d)
(e)
7. Conclusion
Mathematical model, describing the steadystate response of an amperometric biosensor based on parallel substrate conversion (catalaseperoxidase biochemical reactions), is presented. A nonlinear time independent coupled system of partial differential equations has been solved analytically using a new Homotopy perturbation method. Analytical expressions for substrate and product concentrations as function of single dimensionless parameter or are derived. The optimum value of or is also reported. The analytical expressions provided are much simpler and more efficient to understand and to predict the behavior of the system. Our analytical solution for substrates and products are compared with numerical solution providing consistent results. Concentration of substrates, products, and current depends upon reaction/diffusion parameters , , , , and . The response of the steadystate current with respect to all parameters is discussed. For all different concentrations, the analytical results matched well with the simulated result. The close matching of simulated result and the analytical data shows that our proposed solution is able to simulate the dynamic performance of sensing system using the parameters. The analytical results provided here are useful in sensor design, optimization, and prediction of electrode response.
Appendices
A. Approximate Analytical Solution of the Nonlinear Equations (5)–(8) Using the Boundary Conditions (9) by New Homotopy Perturbation Method
In this appendix, we indicated how to find the solution of (5)–(8). We first construct the Homotopy as follows: The boundary conditions are as follows: The approximate analytical solutions of (A.1) to (A.4) are given by Substituting (A.6) into (A.1)–(A.4) and arranging the coefficients of powers , Solving (A.7) for the boundary conditions (A.5), we can find the following results: Since it is a new approach to HPM, we need not to find further iterations. We consider , , , and as the better approximate solutions. According to the HPM, we can conclude that After putting (A.8) into (A.9), respectively, the final results can be described in (11) to (14) in the text.
B. Matlab Program to Find the Numerical Solution of (5)–(8)
See Algorithm 1.

Nomenclature
:  Substrate concentrations () 
:  Enzyme concentrations () 
:  Product concentrations () 
:  Diffusion coefficients () 
:  Reaction rate constant () 
:  Thickness of the enzyme and diffusion layer () 
:  Bulk concentrations () 
:  Space () 
:  Current (). 
:  Dimensionless substrate concentrations 
:  Dimensionless product concentrations 
:  Dimensionless reaction/diffusion parameters 
:  Dimensionless parameters 
:  Dimensionless space 
:  Dimensionless parameter 
:  Dimensionless parameter. 
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the CSIR (no. 01 (2442)/10/EMRII), Government of India.
References
 A. Turner, G. Wilson, and I. Kaube, Biosensors: Fundamentals and Applications, Oxford University Press, Oxford, UK, 1987.
 F. W. Scheller, F. Schubert, and J. Fedrowitz, Frontiers in Biosensorics, Birkhauser, Basel, Switzerland, 1997.
 A. Sadana and N. Sadana, Handbook of Biosensors and Biosensor Kinetics, Elsevier, Amsterdam, The Netherlands, 2010.
 D. M. Fraser, Biosensors in the Body, Continuous in Vivo Monitoring, John Wiley and Sons, New York, NY, USA, 1997.
 L. S. Ferreira, M. B. De Souza Jr., J. O. Trierweiler, O. Broxtermann, R. O. M. Folly, and B. Hitzmann, “Aspects concerning the use of biosensors for process control: experimental and simulation investigations,” Computers and Chemical Engineering, vol. 27, no. 89, pp. 1165–1173, 2003. View at: Publisher Site  Google Scholar
 J. R. D. Corcuera, R. Cavalieri, J. Powers, and J. Tang, “Amperometric enzyme biosensor optimization using mathematical modeling,” in Proceeding of the/Case 18th Annual International Meeting (ASAE '04), page papers No. 047030, American Society of Agricultural Engineers, Ottawa, Canada, August 2004. View at: Google Scholar
 C. Amatore, A. I. Oleinick, and I. Svir, “Construction of optimal quasiconformal mappings for the 2Dnumerical simulation of diffusion at microelectrodes. Part 1: principle of the method and its application to the inlaid disk microelectrode,” Journal of Electroanalytical Chemistry, vol. 597, no. 1, pp. 69–76, 2006. View at: Publisher Site  Google Scholar
 I. Stamatin, C. Berlic, and A. Vaseashta, “On the computeraided modelling of analytereceptor interactions for an efficient sensor design,” Thin Solid Films, vol. 495, no. 12, pp. 312–315, 2006. View at: Publisher Site  Google Scholar
 L. D. Mell and J. T. Maloy, “A model for the amperometric enzyme electrode obtained through digital simulation and applied to the immobilized glucose oxidase system,” Analytical Chemistry, vol. 47, no. 2, pp. 299–307, 1975. View at: Google Scholar
 J. P. Kernevez, Enzyme Mathematics. Studies in Mathematics and Its Applications, Elsevier Science, Amsterdam, The Netherlands, 1980.
 J. Kulys, “The development of new analytical systems based on biocatalysts,” Analytical Letters, vol. 14, article 377, 1981. View at: Google Scholar
 P. N. Bartlett and K. F. E. Pratt, “Modelling of processes in enzyme electrodes,” Biosensors and Bioelectronics, vol. 8, no. 910, pp. 451–462, 1993. View at: Publisher Site  Google Scholar
 M. E. G. Lyons, J. C. Greer, C. A. Fitzgerald, T. Bannon, and P. N. Barlett, “Reaction/diffusion with MichaelisMenten kinetics in electroactive polymer films. Part 1. The steadystate amperometric response,” Analyst, vol. 121, no. 6, pp. 715–731, 1996. View at: Publisher Site  Google Scholar
 P. Pérusse and D. Leech, “Dual electrode cyclic voltammetry under computer control using graphical programming of a bipotentiostat,” Instrumentation Science and Technology, vol. 28, no. 1, pp. 59–70, 2000. View at: Publisher Site  Google Scholar
 J. J. Kulys, V. V. Sorochinskii, and R. A. Vidziunaite, “Transient response of bienzyme electrodes,” Biosensors, vol. 2, no. 3, pp. 135–146, 1986. View at: Google Scholar
 V. V. Sorochinskii and B. I. Kurganov, “Steadystate kinetics of cyclic conversions of substrate in amperometric bienzyme sensors,” Biosensors and Bioelectronics, vol. 11, no. 3, pp. 225–238, 1996. View at: Publisher Site  Google Scholar
 T. Schulmeister, J. Rose, and F. Scheller, “Mathematical modelling of exponential amplification in membranebased enzyme sensors,” Biosensors and Bioelectronics, vol. 12, no. 910, pp. 1021–1030, 1997. View at: Publisher Site  Google Scholar
 A. Neykov and V. Rangelova, “Mathematical modeling of the biosensor systems,” Biotechnology and Biotechnological Equipment, vol. 1998, no. 2, pp. 100–109, 1998. View at: Google Scholar
 R. Baronas, F. Ivanauskas, and J. Kulys, “Modelling a biosensor based on the heterogeneous microreactor,” Journal of Mathematical Chemistry, vol. 25, no. 23, pp. 245–252, 1999. View at: Google Scholar
 W. F. Ames, Numerical Methods for Partial Differential Equations, Academic Press, New York, NY, USA, 2nd edition, 1977.
 A. Eswari and L. Rajendran, “Analytical solution of steady state current at a microdisk biosensor,” Journal of Electroanalytical Chemistry, vol. 641, no. 12, pp. 35–44, 2010. View at: Publisher Site  Google Scholar
 P. Manimozhi, A. Subbiah, and L. Rajendran, “Solution of steadystate substrate concentration in the action of biosensor response at mixed enzyme kinetics,” Sensors and Actuators, B: Chemical, vol. 147, no. 1, pp. 290–297, 2010. View at: Publisher Site  Google Scholar
 A. Meena and L. Rajendran, “Mathematical modeling of amperometric and potentiometric biosensors and system of nonlinear equations—homotopy perturbation approach,” Journal of Electroanalytical Chemistry, vol. 644, no. 1, pp. 50–59, 2010. View at: Publisher Site  Google Scholar
 A. Eswari and L. Rajendran, “Analytical solution of steadystate current an enzymemodified microcylinder electrodes,” Journal of Electroanalytical Chemistry, vol. 648, no. 1, pp. 36–46, 2010. View at: Publisher Site  Google Scholar
 S. Loghambal and L. Rajendran, “Mathematical modeling in amperometric oxidase enzymemembrane electrodes,” Journal of Membrane Science, vol. 373, no. 12, pp. 20–28, 2011. View at: Publisher Site  Google Scholar
 A. Eswari and L. Rajendran, “Analytical expressions pertaining to the concentration of catechol, oquinone and current at PPOmodified microcylinder biosensor for diffusionkinetic model,” Journal of Electroanalytical Chemistry, vol. 660, no. 1, pp. 200–208, 2011. View at: Publisher Site  Google Scholar
 K. Venugopal, A. Eswari, and L. Rajendran, “Mathematical model for steady state current at PPOmodified microcylinder biosensors,” Journal of Biomedical Science and Engineering, vol. 4, pp. 631–641, 2011. View at: Google Scholar
 L. Rajendran and S. Anitha, “Reply to ‘ Comments on analytical solution of amperometric enzymatic reactions based on HPM ’,” Electrochimica Acta, vol. 102, pp. 474–476, 2013. View at: Google Scholar
 V. Aseris, R. Baronas, and J. Kulys, “Modeling the biosensor utilizing parallel substrate conversion,” Journal of Electroanalytical Chemistry, vol. 685, pp. 63–71, 2012. View at: Google Scholar
 J.H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar
 Q. K. Ghori, M. Ahmed, and A. M. Siddiqui, “Application of homotopy perturbation method to squeezing flow of a newtonian fluid,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 179–184, 2007. View at: Google Scholar
 T. Öziş and A. Yildirim, “A comparative study of He's homotopy perturbation method for determining frequencyamplitude relation of a nonlinear oscillator with discontinuities,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 243–248, 2007. View at: Google Scholar
 S.J. Li and Y.X. Liu, “An improved approach to nonlinear dynamical system identification using PID neural networks,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 2, pp. 177–182, 2006. View at: Google Scholar
 M. M. Mousa and S. F. Ragab, “Application of the homotopy perturbation method to linear and nonlinear schrödinger equations,” Zeitschrift fur Naturforschung— A Journal of Physical Sciences, vol. 63, no. 34, pp. 140–144, 2008. View at: Google Scholar
 J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, pp. 257–262, 1999. View at: Google Scholar
 R. D. Skeel and M. Berzins, “A method for the spatial discretization of parabolic equations in one space variable,” SIAM Journal on Scientific & Statistical Computing, vol. 11, pp. 1–32, 1990. View at: Google Scholar
 “MATLAB 6. 1, The Math Works Inc., Natick, MA,” 2000, http://www.scilabenterprises.com. View at: Google Scholar
Copyright
Copyright © 2014 T. Praveen and L. Rajendran. 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.