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 glucose-6-phosphate (substrate), potassium dihydrogen phosphate (inhibitor), oxygen (co-substrate), 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. Inhibition-based 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 inhibition-based biosensors for food safety and environmental monitoring have been discussed by Amine et al. [4]. The application of whole-cell-based 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 polymer-modified 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 glucose-6-phosphate 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 glucose-6-phosphate (C6H13O9P) with the liberation of glucose (C6H12O6) and inorganic phosphate (). In the second reaction with glucose oxidase (GOD), the oxygen present oxidizes the glucose into glucono delta-lactone (C6H10O6) and hydrogen peroxide (H2O2). 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 glucose-6-phosphate (substrate), potassium dihydrogen phosphate (inhibitor), oxygen (co-substrate), 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 delta-lactone. Figure 1 presents the overall reaction scheme for a multienzyme inhibitor system. We assume that diffusion is one-dimensional 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, co-substrate, 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 glucose-6-phosphate, 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 26. 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).

function pdex4
m = 0;
x = linspace(0,1);
t=linspace(0,100000);
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
u2 = sol(:,:,2);
u3 = sol(:,:,3);
u4 = sol(:,:,4);
u5 = sol(:,:,5);
figure
plot(x,u1(end,:))
title(‘u1(x,t)’)
xlabel(‘Distance x’)
ylabel(‘u1(x,2)’)
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
figure
plot(x,u2(end,:))
title(‘u2(x,t)’)
xlabel(‘Distance x’)
ylabel(‘u2(x,2)’)
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
figure
plot(x,u3(end,:))
title(‘u3(x,t)’)
xlabel(‘Distance x’)
ylabel(‘u3(x,2)’)
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
figure
plot(x,u4(end,:))
title(‘u4(x,t)’)
xlabel(‘Distance x’)
ylabel(‘u4(x,2)’)
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
figure
plot(x,u5(end,:))
title(‘u5(x,t)’)
xlabel(‘Distance x’)
ylabel(‘u5(x,2)’)
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
function c,f,s = pdex4pde(x,t,u,DuDx)
c = 1; 1; 1; 1; 1;
f = 1; 1; 1; 1; 1.* DuDx;
a=0.01;b=0.01;p=0.01;e=0.01;
rp=0.01;kp=0.01;
rs=0.01;ki=0.01;kc=0.01;kb=5;
F1 =-(rs2*u(1))/(1+a*u(2)+b*u(1));
F2 =-(ki2*b*u(1))/(1+a*u(2)+b*u(1));
F3 =-(kc2)/(1+p/u(4)+e/u(3));
F4 =-(rp2)/(1+p/u(4)+e/u(3))+ (kp2*u(1))/(1+a*u(2)+b*u(1));
F5 =(kb2)/(1+p/u(4)+e/u(3));
s=F1; F2; F3; F4; F5;
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
function u0 = pdex4ic(x) %create a initial conditions
u0 = 1; 1; 1; 1; 1;
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
function pl,ql,pr,qr=pdex4bc(xl,ul,xr,ur,t) %create a boundary conditions
pl = ul(1)-1; ul(2)-1; ul(3)-1; ul(4)-0.01; ul(5)-0;
ql = 0; 0; 0; 0; 0;
pr = 0; ur(2)-0; ur(3)-0; 0; ur(5)-0;
qr = 1; 0; 0; 1; 0;

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 glucose-6-phosphate , 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 co-substrate against dimensionless distance is plotted in Figure 4. From these Figures, we conclude that the value of decreases when the dimensionless co-substrate Thiele modulus increases. Furthermore, the co-substrate profile increases when dimensionless Michaelis-Menten 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 Michaelis-Menten 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 Michaelis-Menten 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 Michaelis-Menten 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 Michaelis-Menten constant of glucose decreases for some fixed values of .

7. Conclusions

The system of nonlinear, steady-state 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 glucose-6-phosphate, 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 glucose-6-phosphate 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 glucose-6-phosphate (mol/cm3)
:Concentration of potassium dihydrogen phosphate (mol/cm3)
:Concentration of oxygen (mol/cm3)
:Concentration of glucose (mol/cm3)
: Concentration of inorganic phosphate, (mol/cm3)
: Concentration of hydrogen peroxide (mol/cm3)
: Concentration of glucono delta-lactone, C6H10O6 (mol/cm3)
:Acid phosphatase (AP)
: Glucose oxidase (GOD)
: Reaction velocity of acid phosphatase (mol/s cm3)
: Reaction velocity of glucose oxidase (mol/s cm3)
:Reaction constant for glucose-6-phosphate (mol/cm3)
:Reaction constant for potassium dihydrogen phosphate (mol/cm3)
:Reaction constant for oxygen (mol/cm3)
:Reaction constant for glucose (mol/cm3)
:Diffusion coefficient of glucose-6-phosphate (cm2/s)
:Diffusion coefficient of potassium dihydrogen phosphate (cm2/s)
:Diffusion coefficient of oxygen (cm2/s)
:Diffusion coefficient of glucose (cm2/s)
:Diffusion coefficient of hydrogen peroxide (cm2/s)
:Distance (cm)
:Concentration of glucose-6-phosphate in bulk solution (mol/cm3)
:Concentration of potassium dihydrogen phosphate in bulk solution (mol/cm3)
:Concentration of oxygen in bulk solution (mol/cm3)
:Thickness of the membrane (cm)
:Output current at the electrode surface (ampere).
Dimensionless Variables
:Dimensionless concentration of glucose-6-phosphate
:Dimensionless concentration of potassium dihydrogen phosphate
:Dimensionless concentration of oxygen
:Dimensionless concentration of glucose
:Dimensionless concentration of hydrogen peroxide
:Dimensionless distance
:Dimensionless current.
Greek Symbols
:Inhibitor inhibition modulus (dimensionless)
:Substrate inhibition modulus (dimensionless)
:Michaelis-Menten constants (dimensionless)
:Substrate Thiele modulus (dimensionless)
:Inhibitor Thiele modulus (dimensionless)
:Co-substrate 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/PC-50/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.