Reliable Method for Steady-State Concentrations and Current over the Diagnostic Biosensor Transducers ()
1. Introduction
A biosensor is an analytical device used for the detection of an analyte that combines a biological component with a physicochemical detector [1] [2] . The earlier biosensors were catalytic systems that integrated especially enzymes with transducers that converted the biological response into an electronic signal. The next generation of biosensors, took advantage of different biological elements, such as antibodies, receptors (natural or synthetic), or nucleic acids [3] [4] .
Biosensors for environmental application include the detection of harmful bacteria or pesticides in air, water, or food. New technologies are likely to encompass all-printed systems capitalising on the printed electronics revolution and systems with high compatibility with future mobile technology such as tablets and 4G phones [5] .
Rangelova et al. [6] [7] [8] described the models in biosensors and investigated the influence of the diffusion and kinetic parameters on the response of the biosensor. Tothil et al., [9] deals with the recent developments in biosensors and their potential use in the agricultural diagnostic market. Mishra et al. [10] reviewed various cancer biomarkers in saliva and compared the biomarkers efficacy with traditional diagnostics and state-of-the-art bioelectronics. Cortina et al. [11] presented the development and validation of a portable, robust and inexpensive electrochemical magnetic biosensor.
Lawal et al. [12] summarized the fabrication of carbon nanotubes-based electrochemical biosensors. They also discussed the synthesis, along with the application of carbon nanotubes to the assembly of carbon nanotube-based electrochemical sensors, its analytical performance and future expectations. Recently Gruhl et al. [13] described the latest applications of biosensors in diagnostic applications. In this paper the current state and future trends of biosensors are presented. Also Mascini et al. [3] reviewed the application of biosensor sin medical diagnostics, taking into account several crucial features.
The numerical method of solving the system of partial differential equations is to make calculation at all intervals of substrates concentration and at different diffusion and enzymatic reaction rates. The diffusion equations [14] , containing a nonlinear term related to the enzymatic reaction, are carried out using the implicit difference scheme [15] . In recent years, analytical solutions are reported for various types of biosensors [16] [17] [18] [19] . The analytical results of diagnostic biosensor are relevant because their solutions describe important applications such as bioreactors and biofuel cells, among others [20] [21] .
To the researcher’s knowledge no rigorous analytical solution of substrate concentration product with concentration profiles co-substrate concentration and corresponding current response has been derived for all possible values of parameters under steady-state conditions [22] . The purpose of this communication is to derive approximate analytical expressions for the steady-state concentrations and current over the diagnostic of biosensor transducers for first order, Michaelis-Menten and ping-pong kinetics using Homotopy perturbation method.
2. Mathematical Description of the Boundary Value Problem
Only biosensors systems will be investigated in the active membrane, because it is known that the concentrations of substrate, Product and co-substrate in other two membranes are changed linearly. Biosensors are function under diffusion control. It is assumed that the electrode has symmetrical geometry and the enzyme is homogeneously distributed in the active membrane. The diffusion is one dimensional in space and is described with the second Fick’s law. The two parameters diagnostic biosensor transducers are based on oxygen electrode. The steady-state reaction-diffusion equation for biosensor systems in the dynamic mode has the following form [22] :
(1)
The non-dimensional coordinates, variables and parameters are as follows:
(1a)
where are diffusion coefficients for substrate, co-substrate and product. denotes the reaction constant for concentration profiles (), respectively. is the enzyme rate and represents the coordinate distance. The diagnosis of the biosensor system depends on the enzyme kinetics and the enzyme reaction as well as on the basic transducer. The kinetics is distinguished in to three kinds:
First order kinetic:
(2)
Michaelis-Menten kinetic:
(3)
Two substrateping-pong kinetic:
(4)
The three types of biosensors can be described with the following system of differential equations:
Equations (5)-(13) are subjected to the following boundary conditions:
(14)
where represents the thickness of active membrane, is the Thiele Module, is diffusion coefficient of product, is the diffusion coefficient of co-substrate and is reaction rate constant for co-substrate. The initial current of the biosensor system is recorded normally in substrate, product and co-substrate concentrations at the electrode and are as follows:
(15)
(16)
(17)
where is the number of electrons taking part in electrochemical reaction, is the Faraday’s number, and is the area of the electrode surface [m2].
Analytical Solutions of Concentrations of Substrate, Product and Co-Substrate under Steady-State Condition Using the New Homotopy Perturbation Method
Recently, many authors have applied the HPM to solve the various non-linear problem in engineering sciences [23] - [28] . This method is a combination of Homotopy in topology and classic perturbation techniques. The HPM has uniqueness in its applicability, accuracy, and efficiency. Recently, a new approach of HPM with zeroth iteration has been applied to solve the nonlinear problem. In this work, a new approach to Homotopy perturbation method is used (Appen- dix A and Appendix B) to solve the nonlinear differential Equations (5)-(13). Using this method, the analytical expression of the concentration of substrate, Product and co-substrate can be obtained as follows:
The first order kinetic Equations (18)-(20) represent the simple and closed form of analytical expressions of concentrations of substrate, product and co- substrate for all possible values of the parameters. By using Equations (5)-(7) with boundary conditions (14), the following relation is also obtained:
First order kinetic
The Michaelis?Menten kinetic Equations (24)-(26) represent the simple and closed form of analytical expressions of concentrations of substrate, product and co-substrate for all possible values of the parameters. By using Equations (8)-(10) with boundary conditions (14), the following relation is also obtained:
Michaelis-Menten kinetic:
The Ping pong kinetic: Equations (30)-(32) represent the simple and closed forms of analytical expressions of concentrations of substrate, product and co- substrate for all possible values of the parameters. By using Equations (11)-(13) with boundary conditions (14), the following relation is also obtained:
Ping pong kinetic:
where, and
3. Numerical Simulation
The system of non-linear differential Equations (5)-(13) with boundary conditions (14) have been solved numerically using MATLAB software. A MATLAB script pdex4 is provided in Appendix C. In Figure 4, Tables 1-3 the numerical solutions are compared with the analytical results. The maximum average relative error between our analytical and numerical result is 0.94% for first order kinetic, 0.91% for Michaelis?Menten kinetic and 1.54% for Ping pong kinetic.
4. Results and Discussion
The dinensonless non-linear differential equations are solved using a new Homotopy perturbation method. Equations (18)-(20), (24)-(26) and (30)-(32) re- present the analytical expression of the concentrations of substrate, product and co-substrate for various values of Thiele modulus and the dimensionless parameters for first order, Michalies-Menten and Ping-Pong kinetics respectively. The analytical results are compared with the numerical results.
Table 1. Comparison of dimensionless concentrations (Equations (18)-(20)) and numerical simulation for first order kinetics when fixed value of .
Table 2. Comparison of dimensionless concentrations (Equation (24)-(26)) and numerical simulation for Michaelis-Menten kinetics when fixed value of.
Table 3. Comparison of dimensionless concentrations (Equations (30)-(32)) and numerical simulation for ping pong kinetics when fixed value of .
Figures 1-3 show the plots of all the concentrations versus dimensionless distance for various values of parameters. For, biosensors act in diffusion regime, and for, the biosensors act in rule of limiting kinetic. Reaction rate constant for substrate is dependable from enzyme concentration and characterized enzyme. There are different enzymes for various tissues, where as reaction rate constant is permanent for the given tissue. The strong affinity between enzyme and substrate shows low value of kinetics and poor affinity shows high value.
The substrate concentration approaches unity at. The substrate concentration increases with decreasing Thiele module. When the ratio of diffusion coefficient increases, the Thiele module increases the product. The concentration of co-substrate increases, when Thiele module decreases. If the ratio of diffusion coefficient and ratio of reaction rate constant increases, the concentration of co-substrate decreases.
Figure 4 focuses the concentration on the first-order kinetics, Michaelis- Menten kinetics, and ping-pong kinetic mechanism of the substrate, product and co-substrate for particular values of parameters. The analytical results are compared with the numerical results as given in Tables 1-3 for fixed values of parameters and satisfactory agreement is noted. In all the cases, the average relative error is less than 1.54%.
Figure 5 represents the concentration of the substrate, product and co-sub- strate vursus distance for the first-order kinetics, Michaelis?Menten kinetics, and ping-pong kinetic mechanism for particular values of parameters. From this figure, it is inferred that the concentration of the ping-pong kinetics largely corresponds to the other two mechanisms. But for the co substrate, the concentration does not show much difference. From these Figure, it is concluded that the dimensionless concentration of substrate and co-substrateis greater for the Ping-Pong than the first order and M-M kinetics.
Figure 6 and Figure 7 represent the dimensionless current profiles of product and co-substrate for various values of dimensionless parameters. The current depends on the product and co-substrate gradient at the electrode surface. has no influenced over biosensor response, but and has been increasing the value of the diffusion constant and leads to small values of response time.
4. Conclusion
In this work, a mathematical model that describes the steady-state response of a two parameters diagnostic of biosensor is discussed. Anew Homotopy perturbation method is employed to solve the system of steady-state non-linear differential equations for three types of kinetics. Analytical expressions corresponding to substrate, product and co-substrate concentrations are derived as the function of dimensionless parameters. For all different concentrations, the analytical results match well with the simulated results. The analytical results provided in this work are useful to understand the behaviour of the system. The extension
of the procedure to other systems nonlinear equation in various type of biosensor seems to be possible.
Appendix A
Approximate analytical solutions for Equations (8)-(10) (Michaelis-Men- ten kinetic) using HPM:
In order to solve Equation (8-10) by means of the new HPM, first the Zeroth order deformation equation is constructed.
(A1)
(A2)
(A3)
The approximate solutions of Equations (A1), (A2) and (A3) are as follows
(A4)
Substituting (A4) in Equation (A1) and equating the like powers of p, we get
(A5)
(A6)
Substituting (A4) in Equation (A2) and equating the like powers of p, we get
(A7)
(A8)
Substituting (A4) in Equation (A3) and equating the like powers of p, we get
(A9)
(A10)
The boundary conditions in Equation (14) becomes
(A11)
and
(A12)
Now by applying the boundary conditions (A11) in (A5), (A7) and (A9), we get
(A13)
(A14)
(A15)
Substituting the values of and in Equation (A6), (A8) and (A10), and respectively by solving the equations using the boundary conditions (A12), the following results are obtained:
(A16)
(A17)
(A18)
Adding Equations (A13) and (A16), we get Equation (24) in the text. Similarly, Equation (25) and Equation (26) are got in the text.
Appendix B
Approximate analytical solutions for Equations (11-13) (ping-pong kinetic) using HPM:
In order to solve Equations (11-13) by means of the new HPM, first the Zeroth order deformation equation is constructed.
(B1)
(B2)
(B3)
The approximate solutions of Equations (B1), (B2) and (B3) are as follows
(B4)
Substituting (B4) in Equations (B1) and equating the like powers of p, we get
(B5)
(B6)
Substituting (B4) in Equation (B2) and equating the like powers of p, we get
(B7)
(B8)
Substituting (B4) in Equation (B3) and equating the like powers of p, we get
(B9)
(B10)
The boundary conditions in Equation (14) becomes
(B11)
and
(B12)
Now by applying the boundary conditions (B11) in (B5), (B7) and (B9), we get
(B13)
(B14)
(B15)
Substituting the values of and in Equation (B6), (B8) and (B10), and solving the equations using the boundary conditions (B12), the following results are obtained:
(B16)
(B17)
(B18)
Adding Equations (B13) and (B16), we get Equation (30) in the text. Similarly, Equation (31) and Equation (32) are got in the text.
Appendix C
Scilab/Matlab program for the numerical solution of the system of non- linear Equations (5-13).
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);
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)')
%--------------------------------------------------------------
function [c,f,s] = pdex4pde(x,t,u,DuDx)
c = [1; 1; 1];
f = [1; 1; 1] .* DuDx;
l=0.1;mu=0.5;q=1;p=5;
F=-q^2*u(1);
F1=l*q^2*u(1);
F2=-mu*p*q^2*u(1);
s=[F; F1; F2];
% --------------------------------------------------------------
function u0 = pdex4ic(x);
u0 = [1; 1; 1];
% --------------------------------------------------------------
function [pl,ql,pr,qr]=pdex4bc(xl,ul,xr,ur,t)
pl = [ul(1)-1;ul(2)-0;ul(3)-1];
ql = [0; 0; 0];
pr = [0; ur(2)-0;ur(3)-0];
qr = [1;0; 0];
Nomenclature
Symbol
Dimensionless parameters
Concentration of substrate (mmol)
Concentration of substrate (mmol)
Concentration of substrate (mmol)
Dimensionless concentration of substrate (mmol)
Dimensionless concentration of product (mmol)
Dimensionless concentration of co-substrate (mmol)
Initial concentration of substrate (mmol)
Initial concentration of product (mmol)
Initial concentration of co-substrate (mmol)
Reaction rate constant (mmol)
Diffusion coefficient (m2/s)
Enzymatic rate (mmol/s)
Output current
Number of electrons
Faraday’s number (C/mol)
Area of the electrode surface (m2)
Dimensionless distance
Greek symbols
Diffusion coefficient of substrate and product (Dimensionless)
Diffusion coefficient of substrate and co-substrate (Dimensionless)
Reaction rate constant for substrate and co-substrate (Dimensionless)
Thiele Module (Dimensionless)
dimension distance (m)
Current (Dimensionless)
Subscripts
Substrate
Product
Co-Substrate