Next Article in Journal
Constituents from the Fruiting Bodies of Trametes cubensis and Trametes suaveolens in Vietnam and Their Anti-Inflammatory Bioactivity
Next Article in Special Issue
The Discovery of New Drug-Target Interactions for Breast Cancer Treatment
Previous Article in Journal
Facile Synthesis of 5-Aryl-N-(pyrazin-2-yl)thiophene-2-carboxamides via Suzuki Cross-Coupling Reactions, Their Electronic and Nonlinear Optical Properties through DFT Calculations
Previous Article in Special Issue
Theoretical Analysis on Absorption of Carbon Dioxide (CO2) into Solutions of Phenyl Glycidyl Ether (PGE) Using Nonlinear Autoregressive Exogenous Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Analysis of Reaction–Diffusion Equations Modeling the Michaelis–Menten Kinetics in a Micro-Disk Biosensor

by
Naveed Ahmad Khan
1,
Fahad Sameer Alshammari
2,
Carlos Andrés Tavera Romero
3,
Muhammad Sulaiman
1,* and
Ghaylen Laouini
4
1
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan
2
Department of Mathematics, College of Science and Humanities in Alkharj, Prince Sattam bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
3
COMBA R&D Laboratory, Faculty of Engineering, Universidad Santiago de Cali, Cali 76001, Colombia
4
College of Engineering and Technology, American University of the Middle East, Egaila 54200, Kuwait
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(23), 7310; https://doi.org/10.3390/molecules26237310
Submission received: 10 October 2021 / Revised: 22 November 2021 / Accepted: 24 November 2021 / Published: 2 December 2021
(This article belongs to the Special Issue Advances in the Theoretical and Computational Chemistry)

Abstract

:
In this study, we have investigated the mathematical model of an immobilized enzyme system that follows the Michaelis–Menten (MM) kinetics for a micro-disk biosensor. The film reaction model under steady state conditions is transformed into a couple differential equations which are based on dimensionless concentration of hydrogen peroxide with enzyme reaction ( H ) and substrate ( S ) within the biosensor. The model is based on a reaction–diffusion equation which contains highly non-linear terms related to MM kinetics of the enzymatic reaction. Further, to calculate the effect of variations in parameters on the dimensionless concentration of substrate and hydrogen peroxide, we have strengthened the computational ability of neural network (NN) architecture by using a backpropagated Levenberg–Marquardt training (LMT) algorithm. NNs–LMT algorithm is a supervised machine learning for which the initial data set is generated by using MATLAB built in function known as “pdex4”. Furthermore, the data set is validated by the processing of the NNs–LMT algorithm to find the approximate solutions for different scenarios and cases of mathematical model of micro-disk biosensors. Absolute errors, curve fitting, error histograms, regression and complexity analysis further validate the accuracy and robustness of the technique.

1. Introduction

A biosensor is a device that converts the results of biological processes into analytically accessible data. The amount of product created during the biological process affects the information analysis. Generally, the process comprises of two components named physicochemical transducer and biochemical recognition system [1,2,3]. Physical signals are generated by converting biochemical results with a specified sensitivity, which is generally considered analytic concentration, by a bioreceptor. Another component of a biosensor is transducer that produces electric signals from receptor output categorized by transducer and bioreceptor. Biosensor uses enzymes, cell structure, bioreceptors, antibodies, hormones, nucleic acid, and tissues. Auxiliary enzymes are also used in the development of biosensors. Transducers systems are categorized into thermometric, magnetic, electrochemical, and piezoelectric [4]. Biochemical reactions between an immobilized bimolecular and the target analyte form an electro-chemical biosensor. As a result, the theoretical foundation of a biosensor measures the electric current [5]. Electrochemical biosensors are categorized into potentiometric and amperometric, which are used for the mass production [6,7,8]. The functioning of an amperometric biosensor is based on calculating the Faraday current, which is determined when the current is constant at the electrode. As a result, the current is generated by the product’s oxidation or reduction [9,10,11]. Michaelis–Menten kinetic equations are commonly used to simulate the process.
In 1975, Mell [12] developed a mathematical model for amperometric two enzyme biosensors with multiple enzymes. Various numerical and analytical techniques have been used to calculate the analytical and approximate solution to the model. Che and Dong [13] reported analytical expressions for the steady-state concentration of current at the micro-disk chemical sensor. An approximate solution in integral form has been calculated by Phanthong [14] for micro-disk biosensor. Eswari derived the analytical solution in terms of enzyme kinetics and film thickness [15] for all parameters of the Michaelis–Menten constant. Manimozhi [16] study the steady-state concentration of substrate of nonlinear equations representing the action of the biosensor using homotopy perturbation method (HPM). Eswari in [15,17] 2010, derived an analytical solution for steady-state current on enzyme-modified microcylinder electrodes, micro-disk, and spherical biosensor. Loghambal [18] uses the asymptotic method (AM) to study the modeling of amperometric oxidase on enzyme membrane electrodes.
Recently, various techniques have been implemented on nonlinear ECE reactions to study the steady-state concentration of substrate and products on rotating disk electrodes [19]. M.C. Devi [20] uses hyperbolic function method (HFM) to find the analytical expression for the EC-catalytic mechanism of the first order. Variational iteration method (VIM), differential transformation method [21,21] and homotopy perturbation methods [22] has been widely used to calculate the approximate series solutions for the mathematical model of micro-disk biosensors. The implementation of these techniques has not been straightforward. For most problems, the above-discussed techniques fail to converge the solution into closed form and are time-consuming. In recent times, a stochastic numerical technique based on artificial intelligence has been developed to solve stiff nonlinear problems arising in various fields. Such stochastic computing techniques use artificial neural networks to model approximate solutions. These numerical solvers have wide applications in various fields including petroleum engineering [23], wireless communication [24], heat transfer [25,26,27], fuzzy systems [28], plasma system [29], civil engineering [30,31], wire coating dynamics [32] and Diabetic retinopathy classification [33]. The techniques mentioned earlier inspire the authors to explore and incorporate the soft computing architectures as an alternative, precise and feasible way for solving the mathematical model of micro-disk biosensors. Some highlighted features of the presented study are illustrated as follows:
  • A mathematical model for micro-disk biosensor has been presented to investigate the influence of variations in different parameters on the dimensionless concentration of substrate and hydrogen peroxide.
  • An artificial-neural-networks-based backpropagated Levenberg–Marquardt training (LMT) algorithm is developed to train the hidden neurons, calculate the validation of reference data-set generated by “pdex4” for different cases and scenarios of micro-disk biosensor.
  • Extensive graphical analysis has been conducted based on mean square error (MSE), absolute errors, regression, curve fitting, and error histograms that show the technique’s convergence, accuracy, and computational complexity. 3D plots of dimensionless concentration for substrate ( S ) and hydrogen peroxide ( H ) are plotted against dimensionless distance R and reaction–diffusion parameters to study the behavior and changes in the model.

2. Problem Formulation

In this section, a mathematical derivation of the microdisk biosensor is presented. A generalized form of the polymer solution of drop coating is considered the special case for micro-disk electrodes. It has been observed that microdisk is insulated through droplets of enzymes or polymers, taking the form of hemisphere on isolating plane. Phanthong and Somasundrum [14] describe the mathematical formulation of a micro-disk biosensor. Figure 1 shows the schematic view of a micro-disk biosensor. It can be observed that the radius of the film is greater than the size of the disk. Measurements can be simplified by moving the micro-disk sensor to the micro-hemisphere. In [14] the micro-disk electrode is modified by using redox polymer. In such case, the enzymatic reaction favors Michaelis–Menten kinetics, and film reaction is given as
S + E 1 k 1 k 2 E 1 S k c a t P + E 2 ,
if the solution is agitated evenly such that S is continuously applied to the film, the mass balance at steady state for S will be given by
D S r 2 d d r r 2 d C S d r k c a t C E C S C S + K M = 0 ,
where C S , C E denotes the concentration profiles of substrate and enzymes, D S represents the diffusion coefficient of reaction, K M is Michaelis constant which is defined as
K M = ( k 1 + k c a t ) k 1 ,
steady state mass balance for H is given as
D H r 2 d d r r 2 d C H d r + k c a t C E C S C S + K M = 0 ,
concentration profile of steady state hydrogen peroxide is denoted by C H . Boundary conditions at surfaces of electrode r 0 and film r 1 are defined as
r = r 0 : d C S d r = 0 , C H = 0 ,
r = r 1 : C S = C S * , C H = 0 ,
where C S * represents the bulk concentration of S.
Now introducing the following set of dimensionless variables
S = C S C S * , H = C H C H * , R = r 1 r 0 , α = C S * K M , γ E = k c a t C E r 0 2 D S K M , γ S = k c a t C E r 0 2 D H K M , D H D S = γ E γ S ,
where S and H represent the dimensionless concentration profiles of substrate and hydrogen peroxide. R is dimensionless distance, γ E , γ S and α are reaction–diffusion and saturation parameters respectively. Thus, dimensionless form of Equations (2) and (3) along with boundary conditions can be written as
d 2 S d R 2 + 2 R d S d R γ E S 1 + α S = 0 ,
d 2 H d R 2 + 2 R d H d R + γ S S 1 + α S = 0 ,
with
d S d R = 0 , H = 0 w h e n R = 1 ,
S = 1 , H = 0 w h e n R = r 1 / r 0 .

3. Reference Solutions

Various analytical and numerical techniques have been previously developed to solve the mathematical model of a micro-disk biosensor in the literature. A number of these analytical methods focus on obtaining an analytical expression for the model. These methods include variational iteration method (VIM) [34], Li-He’s variational principle methods [35,36,37], Akbari-Ganji Method (AGM Method) [38], homotopy perturbation method (HPM) [39], Modified Adomian decomposition method (MADM) [40,41], exp-function method [42], Green’s function iteration method [20] and Taylor series method [43]. An analytical expression for concentration of substrate and hydrogen peroxide obtained by HPM [15] is given as
C S P = r r 0 0.5 r 2 r 0 2 + 0.5 r 1 2 r 0 2 r 1 r 0 ,
C H P = 1 2 r 1 r r 0 2 r 2 r 0 2 r 1 r 0 + r r 0 ,
where C S P = c S c S * c S * γ r 0 2 and C H P = c H D H x r 0 2 D S c s * . Analytical solution by MADM [40,41] are given as
S ( R ) = 1 + γ E 3 ( 1 + α ) r 1 r 0 γ E 6 ( 1 + α ) r 1 r 0 2 γ E 3 ( 1 + α ) R + γ E 6 ( 1 + α ) R 2 ,
H ( R ) = γ S 6 ( 1 + α ) r 1 r 0 + γ S 6 ( 1 + α ) 1 + r 1 / r 0 R γ S 6 ( 1 + α ) R 2 .
Analytical solution obtained by HAM [44] are given as
S ( R ) = 1 + h γ E 2 ( 1 + α ) r 1 r 0 2 h γ E ( 1 + α ) r 1 r 0 + h γ E ( 1 + α ) R h γ E 2 ( 1 + α ) R 2 ,
H ( R ) = h γ S 2 ( 1 + α ) r 1 r 0 h γ S 2 ( 1 + α ) r 1 / r 0 + 1 R + h γ S 2 ( 1 + α ) R 2 .

4. Design Methodology

In this section, a novel machine learning technique based on supervised learning of neurons in artificial neural networks (ANNs) is utilized to study the mathematical model of an immobilized enzyme system that follows the Michaelis–Menten (MM) kinetics for micro-disk biosensor. An Artificial Neural Networks (ANNs) is a collection of interconnected, basic components known as neurons with multiple inputs and a single output, each neuron represents a mapping. The output of a neuron is a function of the sum of its inputs which is generated with the help of activation function. In this paper, multilayer perceptron (MLP) is considered, with an objective to perform the optimization of the hidden units number in the hidden layer. Additionally, the optimization of the connection weights and biases has been conducted. The typical structure of MLP with one hidden layer is given as
S j = i = 1 n ω i j R + β j ,
where, R denotes the input, β is a biased vector and ω i j represents the connection weights. Log-sigmoid is utilized as an activation function for the Feed-forword neural network model which is given as
f j ( x ) = 1 1 + e S j ,
Further, the working procedure of the design soft computing technique has been discussed. Implementation of the proposed technique is based on two steps. Initially, a data set is generated for the mathematical model using an efficient numerical technique of Matlab builtin function known as “pdex4”. The date set of 1143 and 1251 points are generated from 0 to 5 and 0 to 1.5, respectively. Furthermore, an intelligent strength of neural networks with 60 hidden neurons is utilized by using the Levenberg–Marquardt technique to find approximate solutions for different scenarios and cases of the problem. The proposed NN’s-LMT algorithm in the form of a single neuron model is shown through Figure 2. The supervised learning of the Levenberg–Marquardt technique uses the data set generated in the first step by using the “nftool” package of MATLAB. The working procedure of processing data for validation and testing is shown through the flow chart given in Figure 3. Moreover, performance measures are defined in terms of mean square error of the objective function, regression study, error histograms, and absolute errors to study the accuracy and convergence of the design scheme.

5. Experimentation Setup and Discussion

In this section, to check the validity of the proposed technique and study the effect of variations in reaction–diffusion parameters on substrate and hydrogen peroxide concentration profiles, we have considered certain scenarios. The scenarios are based on different values of reaction–diffusion and saturation parameters with fixed values of thin-film r 1 / r 0 . The details of different scenarios and cases are dictated in the flow chart given in Figure 4.
The data set values generated by the numerical technique are used to test, train and validate the data with a probability of 75 % , 15 % , and 15 % , respectively. Figure 5 and Figure 6 represent the convergence of mean square error (MSE) function for scenarios I and II, respectively. The performance values for each case of different scenario lies around 3.9474 × 10 13 , 6.4304 × 10 12 , 2.5348 × 10 13 , 4.56 × 10 11 , 1.3821 × 10 13 , 6.161 × 10 12 , 1.8469 × 10 14 , 1.6659 × 10 12 , 2.9481 × 10 14 , 7.3806 × 10 12 , 1.4037 × 10 14 and 4.4535 × 10 14 . Figure 7 and Figure 8 are plotted to investigate the influence of variations in saturation and reaction–diffusion parameters with fixed values of film thickness. These parameters describe the importance of reaction and diffusion in the enzyme layer. It is evident from the figures that the concentration profile of substrate significantly increases when the values of reaction–diffusion ( γ E ) and saturation parameter α are increased. Additionally, a simultaneous increase is observed in the concentration of hydrogen peroxide and substrate as γ S and α decreases. The approximate solutions for concentration of substrate obtained by design algorithm are compared with the homotopy analysis method (HAM) [44], modified Adomain decomposition method [40], hyperbolic function method [45,46] and numerical solutions by Pdex4 as shown in Table 1. The fitting of approximate solution and targeted date overlap each other with minimum absolute errors as shown in the Figure 9 and Figure 10. The values of absolute errors in solutions for S ( R ) lies around 10 6 to 10 8 , 10 6 to 10 9 , 10 6 to 10 8 , 10 7 to 10 8 , 10 6 to 10 8 and 10 7 to 10 8 respectively. Additionally, AE for the solution concentration of hydrogen peroxide lies around 10 6 to 10 7 , 10 5 to 10 7 , 10 6 to 10 7 , 10 6 to 10 8 , 10 6 to 10 8 and 10 7 to 10 8 . The statistical performance in term of gradient, mu, validation failures and regression analysis are plotted in Figure 11, Figure 12, Figure 13 and Figure 14 and dictated in Table 2 and Table 3 respectively. The value gradient for all the cases of micro-disk biosensor are 9.92 × 10 13 , 4.11 × 10 7 , 9.76 × 10 8 , 9.65 × 10 8 , 9.94 × 10 8 , 4.53 × 10 9 , 9.91 × 10 8 , 8.39 × 10 8 , 9.76 × 10 8 , 8.98 × 10 8 , 2.59 × 10 9 and 9.74 × 10 8 .
Three-dimensional plots of substrate and hydrogen peroxide are plotted against dimensionless distance R, reaction–diffusion parameters, and saturation parameters as shown in Figure 15 and Figure 16 respectively. From the figures, the influence of variations in parameters can be observed. From Figure 15a,b, the increase in reaction–diffusion parameters ( γ E , γ S ) with fixed α = 50 causes a decrease in the concentration profile of substrate. In contrast, a significant increase is observed in the concentration profile of hydrogen peroxide. Figure 16 represents the effect of variations in α . It is evident that slight increase in S ( R ) is observed till α = 40 but increases sharply for α > 40 . Simultaneously, a decrease in the concentration of hydrogen peroxide is observed for R 1 . Normal probability curves are plotted to study the computational complexity. The time taken by the system for obtaining solutions for each scenario lies around 0.45 and 0.50 s, as shown in Figure 17.
To check the performance and validity of the proposed ANN–GNDO–SQP algorithm, we defined different statistical operators along with their global form. The performance operators are fitness functions, Theil’s inequality coefficient (TIC), mean absolute deviations (MAD), Nash Sutcliffe efficiency (NSE), and error in Nash Sutcliffe efficiency (ENSE). A mathematical formulation of these indices is defined as
MAD = 1 N m = 1 N S m ( R ) S ^ m ( R ) ,
TIC = 1 N n = 1 N S m ( R ) S ^ m ( R ) 2 ( 1 N m = 1 N ( S m ( R ) ) 2 + 1 N m = 1 N ( S ^ m ( R ) ) 2 ) ,
NSE = 1 m = 1 N S m ( H ) S ^ m ( H ) 2 m = 1 N ( S m ( H ) S ^ ̲ m ( H ) 2 , S ^ ̲ m ( H ) = 1 N m = 1 N S ^ ( H ) ,
ENSE = 1 NSE .
where, S m is analytical solution and S ^ m represents the approximate solution by proposed algorithm. N denotes the grid points.
The stability, efficiency and accuracy of the proposed algorithm is established by executing the design algorithms for twenty multiple runs. Global values of performance function in terms of mean square error for each scenario lies around 3.8110 10 to 3.7510 12 and 1.1410 11 to 2.7510 13 with standard deviations 5.8910 9 to 1.0710 12 and 1.9810 11 to 3.1610 13 respectively. In addition, from Table 4 and Table 5 the values of MAD, TIC and ENSE are also approaching zero, which shows the perfect modeling of the approximate solutions by the design algorithm.

6. Conclusions

This paper investigated the mathematical model of an immobilized enzyme system that follows the Michaelis–Menten (MM) kinetics for micro-disk biosensors. The model was based on reaction–diffusion phenomena, which was given by a system of nonlinear differential equations. Furthermore, a soft computing technique based on supervised learning of Levenberg–Marquardt backpropagation neural networks is used to calculate substrate and hydrogen peroxide concentration under the influence of variations in several parameters, including reaction diffusion and film thickness, and saturation parameter. The results illustrated in the figures conclude that increase in reaction–diffusion and saturation parameters causes a decrease in substrate concentration while an increase in the concentration profile of hydrogen peroxide. It is also concluded that an increase in saturation parameter directly relates to substrate concentration, while inversely relating to the concentration of hydrogen peroxide. An extensive graphical analysis based on MSE, error histogram, absolute errors, regressions, and computational complexity are conducted, showing the robustness, accuracy, and efficiency of the designed scheme.

Author Contributions

Data curation, N.A.K.; Formal analysis, N.A.K.; Funding acquisition, F.S.A., G.L. and C.A.T.R.; Investigation, N.A.K. and M.S.; Methodology, N.A.K. and M.S.; Project administration, M.S.; Resources, F.S.A., G.L. and C.A.T.R., M.S.; Software, M.S.; Supervision, M.S.; Visualization, N.A.K.; Writing—original draft, N.A.K.; Writing—review and editing, F.S.A., G.L., C.A.T.R., M.S. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Dirección General de Investigaciones of Universidad Santiago de Cali under call No. 01-2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This research has been funded by Dirección General de Investigaciones of Universidad Santiago de Cali under call No. 01-2021.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Turner, A.; Karube, I.; Wilson, G.S. Biosensors: Fundamentals and Applications; Oxford University Press: New York, NY, USA, 1987; p. 770. [Google Scholar]
  2. Monosik, R.; Stredansky, M.; Sturdik, E. Biosensors—Classification, characterization and new trends. Acta Chim. Slovaca 2012, 5, 109–120. [Google Scholar] [CrossRef] [Green Version]
  3. Salomi, R.J.; Sylvia, S.V.; Rajendran, L.; Abukhaled, M. Electric potential and surface oxygen ion density for planar, spherical and cylindrical metal oxide grains. Sens. Actuators B Chem. 2020, 321, 128576. [Google Scholar] [CrossRef]
  4. McNaught, A.D.; Wilkinson, A. Compendium of Chemical Terminology; Blackwell Science: Oxford, UK, 1997; Volume 1669. [Google Scholar]
  5. Thevenot, D.R.; Toth, K.; Durst, R.A.; Wilson, G.S. Electrochemical biosensors: Recommended definitions and classification. Pure Appl. Chem. 1999, 71, 2333–2348. [Google Scholar] [CrossRef] [Green Version]
  6. Ghindilis, A.L.; Atanasov, P.; Wilkins, M.; Wilkins, E. Immunosensors: Electrochemical sensing and other engineering approaches. Biosens. Bioelectron. 1998, 13, 113–131. [Google Scholar] [CrossRef]
  7. Yulaev, M.; Sitdikov, R.; Dmitrieva, N.; Yazynina, E.; Zherdev, A.; Dzantiev, B. Development of a potentiometric immunosensor for herbicide simazine and its application for food testing. Sens. Actuators B Chem. 2001, 75, 129–135. [Google Scholar] [CrossRef]
  8. Wilson, M.S.; Rauh, R.D. Novel amperometric immunosensors based on iridium oxide matrices. Biosens. Bioelectron. 2004, 19, 693–699. [Google Scholar] [CrossRef]
  9. Rogers, K.R. Biosensors for environmental applications. Biosens. Bioelectron. 1995, 10, 533–541. [Google Scholar] [CrossRef]
  10. Krishnamoorthy, S.; Makhijani, V.; Lei, M.; Giridharan, M.; Tisone, T. Computational studies of membrane-based test formats. In Proceedings of the Technical International Conference on Modeling and Simulation of Microsystems, San Diego, CA, USA, 27–29 March 2000; pp. 590–593. [Google Scholar]
  11. Park, S.J.; Taton, T.A.; Mirkin, C.A. Array-based electrical detection of DNA with nanoparticle probes. Science 2002, 295, 1503–1506. [Google Scholar] [CrossRef] [PubMed]
  12. Mell, L.D.; Maloy, J. Model for the amperometric enzyme electrode obtained through digital simulation and applied to the immobilized glucose oxidase system. Anal. Chem. 1975, 47, 299–307. [Google Scholar] [CrossRef]
  13. Dong, S.; Che, G. Electrocatalysis at a microdisk electrode modified with a redox species. J. Electroanal. Chem. Interfacial Electrochem. 1991, 309, 103–114. [Google Scholar] [CrossRef]
  14. Phanthong, C.; Somasundrum, M. The steady state current at a microdisk biosensor. J. Electroanal. Chem. 2003, 558, 1–8. [Google Scholar] [CrossRef]
  15. Eswari, A.; Rajendran, L. Analytical solution of steady state current at a microdisk biosensor. J. Electroanal. Chem. 2010, 641, 35–44. [Google Scholar] [CrossRef]
  16. Manimozhi, P.; Subbiah, A.; Rajendran, L. Solution of steady-state substrate concentration in the action of biosensor response at mixed enzyme kinetics. Sens. Actuators B Chem. 2010, 147, 290–297. [Google Scholar] [CrossRef]
  17. Eswari, A.; Rajendran, L. Analytical expressions of concentration and current in homogeneous catalytic reactions at spherical microelectrodes: Homotopy perturbation approach. J. Electroanal. Chem. 2011, 651, 173–184. [Google Scholar] [CrossRef]
  18. Loghambal, S.; Rajendran, L. Mathematical modeling in amperometric oxidase enzyme–membrane electrodes. J. Membr. Sci. 2011, 373, 20–28. [Google Scholar] [CrossRef]
  19. Visuvasam, J.; Meena, A.; Rajendran, L. New analytical method for solving nonlinear equation in rotating disk electrodes for second-order ECE reactions. J. Electroanal. Chem. 2020, 869, 114106. [Google Scholar] [CrossRef]
  20. Devi, M.C.; Pirabaharan, P.; Abukhaled, M.; Rajendran, L. Analysis of the steady-state behavior of pseudo-first-order EC-catalytic mechanism at a rotating disk electrode. Electrochim. Acta 2020, 345, 136175. [Google Scholar] [CrossRef]
  21. Odibat, Z.M.; Bertelle, C.; Aziz-Alaoui, M.; Duchamp, G.H. A multi-step differential transform method and application to non-chaotic or chaotic systems. Comput. Math. Appl. 2010, 59, 1462–1472. [Google Scholar] [CrossRef] [Green Version]
  22. Babolian, E.; Azizi, A.; Saeidian, J. Some notes on using the homotopy perturbation method for solving time-dependent differential equations. Math. Comput. Model. 2009, 50, 213–224. [Google Scholar] [CrossRef]
  23. Khan, N.A.; Sulaiman, M.; Aljohani, A.J.; Kumam, P.; Alrabaiah, H. Analysis of Multi-Phase Flow Through Porous Media for Imbibition Phenomena by Using the LeNN-WOA-NM Algorithm. IEEE Access 2020, 8, 196425–196458. [Google Scholar] [CrossRef]
  24. Khan, N.A.; Sulaiman, M.; Aljohani, A.J.; Bakar, M.A. Mathematical models of CBSC over wireless channels and their analysis by using the LeNN-WOA-NM algorithm. Eng. Appl. Artif. Intell. 2022, 107, 104537. [Google Scholar] [CrossRef]
  25. Ahmad, A.; Sulaiman, M.; Alhindi, A.; Aljohani, A.J. Analysis of temperature profiles in longitudinal fin designs by a novel neuroevolutionary approach. IEEE Access 2020, 8, 113285–113308. [Google Scholar] [CrossRef]
  26. Khan, N.A.; Sulaiman, M.; Kumam, P.; Bakar, M.A. Thermal analysis of conductive-convective-radiative heat exchangers with temperature dependent thermal conductivity. IEEE Access 2021, 9, 138876–138902. [Google Scholar] [CrossRef]
  27. Khan, N.A.; Khalaf, O.I.; Romero, C.A.T.; Sulaiman, M.; Bakar, M.A. Application of Euler Neural Networks with Soft Computing Paradigm to Solve Nonlinear Problems Arising in Heat Transfer. Entropy 2021, 23, 1053. [Google Scholar] [CrossRef] [PubMed]
  28. Ahmad, S.; Sulaiman, M.; Kumam, P.; Hussain, Z.; Asif Jan, M.; Mashwani, W.K.; Ullah, M. A novel population initialization strategy for accelerating levy flights based multi-verse optimizer. J. Intell. Fuzzy Syst. 2020, 39, 1–17. [Google Scholar] [CrossRef]
  29. Bukhari, A.H.; Sulaiman, M.; Raja, M.A.Z.; Islam, S.; Shoaib, M.; Kumam, P. Design of a hybrid NAR-RBFs neural network for nonlinear dusty plasma system. Alex. Eng. J. 2020, 59, 3325–3345. [Google Scholar] [CrossRef]
  30. Huang, W.; Jiang, T.; Zhang, X.; Khan, N.A.; Sulaiman, M. Analysis of Beam-Column Designs by Varying Axial Load with Internal Forces and Bending Rigidity Using a New Soft Computing Technique. Complexity 2021, 2021, 6639032. [Google Scholar] [CrossRef]
  31. Zhang, Y.; Lin, J.; Hu, Z.; Khan, N.A.; Sulaiman, M. Analysis of Third-Order Nonlinear Multi-Singular Emden–Fowler Equation by Using the LeNN-WOA-NM Algorithm. IEEE Access 2021, 9, 72111–72138. [Google Scholar] [CrossRef]
  32. Khan, N.A.; Sulaiman, M.; Kumam, P.; Aljohani, A.J. A new soft computing approach for studying the wire coating dynamics with Oldroyd 8-constant fluid. Phys. Fluids 2021, 33, 036117. [Google Scholar] [CrossRef]
  33. Ali, A.; Qadri, S.; Khan Mashwani, W.; Kumam, W.; Kumam, P.; Naeem, S.; Goktas, A.; Jamal, F.; Chesneau, C.; Anam, S.; et al. Machine learning based automated segmentation and hybrid feature analysis for diabetic retinopathy classification using fundus image. Entropy 2020, 22, 567. [Google Scholar] [CrossRef]
  34. He, J.H. Homotopy perturbation method: A new nonlinear analytical technique. Appl. Math. Comput. 2003, 135, 73–79. [Google Scholar] [CrossRef]
  35. He, J.H. Lagrange crisis and generalized variational principle for 3D unsteady flow. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 1189–1196. [Google Scholar] [CrossRef]
  36. Li, X.; He, J.H. Variational multi-scale finite element method for the two-phase flow of polymer melt filling process. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 1407–1426. [Google Scholar] [CrossRef]
  37. He, J.H.; Sun, C. A variational principle for a thin film equation. J. Math. Chem. 2019, 57, 2075–2081. [Google Scholar] [CrossRef]
  38. Akbari, M.; Ganji, D.; Majidian, A.; Ahmadi, A. Solving nonlinear differential equations of Vanderpol, Rayleigh and Duffing by AGM. Front. Mech. Eng. 2014, 9, 177–190. [Google Scholar] [CrossRef]
  39. Yu, D.N.; He, J.H.; Garcıa, A.G. Homotopy perturbation method with an auxiliary parameter for nonlinear oscillators. J. Low Freq. Noise Vib. Act. Control 2019, 38, 1540–1554. [Google Scholar] [CrossRef]
  40. Hasan, Y.Q.; Zhu, L.M. Modified Adomian decomposition method for singular initial value problems in the second-order ordinary differential equations. Surv. Math. Appl. 2008, 3, 183–193. [Google Scholar]
  41. Wazwaz, A.M. A reliable modification of Adomian decomposition method. Appl. Math. Comput. 1999, 102, 77–86. [Google Scholar] [CrossRef]
  42. He, J.H.; Ji, F.Y. Taylor series solution for Lane–Emden equation. J. Math. Chem. 2019, 57, 1932–1934. [Google Scholar] [CrossRef]
  43. He, J.H. Taylor series solution for a third order boundary value problem arising in architectural engineering. Ain Shams Eng. J. 2020, 11, 1411–1414. [Google Scholar] [CrossRef]
  44. Liao, S. Beyond Perturbation: Introduction to the Homotopy Analysis Method; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  45. Eswari, A.; Saravanakumar, S. New Mathematical Analysis for Nonlinear Simultaneous Differential Equation in Micro-Disk Biosensor Using Hyperbolic Function Method. J. Electroanal. Chem. 2021, 880, 114827. [Google Scholar] [CrossRef]
  46. Aljohani, J.L.; Alaidarous, E.S.; Raja, M.A.Z.; Shoaib, M.; Alhothuali, M.S. Intelligent Computing through Neural Networks for Numerical Treatment of Non-Newtonian Wire Coating Analysis Model. Sci. Rep. 2021, 11, 9072. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic of micro-disk biosensor.
Figure 1. Schematic of micro-disk biosensor.
Molecules 26 07310 g001
Figure 2. Architecture of single neuron model.
Figure 2. Architecture of single neuron model.
Molecules 26 07310 g002
Figure 3. Working mechanism of design NNs–LM technique for numerical solution of mathematical model of micro-disk biosensor.
Figure 3. Working mechanism of design NNs–LM technique for numerical solution of mathematical model of micro-disk biosensor.
Molecules 26 07310 g003
Figure 4. Different scenarios and cases of micro-disk biosensor.
Figure 4. Different scenarios and cases of micro-disk biosensor.
Molecules 26 07310 g004
Figure 5. Performance of NNs–LM technique in terms of mean square error for (a,c,e) dimensionless concentration of substrate and (b,d,f) hydrogen peroxide of scenario I.
Figure 5. Performance of NNs–LM technique in terms of mean square error for (a,c,e) dimensionless concentration of substrate and (b,d,f) hydrogen peroxide of scenario I.
Molecules 26 07310 g005
Figure 6. Performance of NNs–LM technique in terms of mean square error for (a,c,e) dimensionless concentration of substrate and (b,d,f) hydrogen peroxide of scenario II.
Figure 6. Performance of NNs–LM technique in terms of mean square error for (a,c,e) dimensionless concentration of substrate and (b,d,f) hydrogen peroxide of scenario II.
Molecules 26 07310 g006
Figure 7. Comparison of approximate solutions obtained by NNs–LMT with numerical solution for different cases of Scenario-I.
Figure 7. Comparison of approximate solutions obtained by NNs–LMT with numerical solution for different cases of Scenario-I.
Molecules 26 07310 g007
Figure 8. Comparison of approximate solutions obtained by NNs–LMT with numerical solution for different cases of Scenario-II.
Figure 8. Comparison of approximate solutions obtained by NNs–LMT with numerical solution for different cases of Scenario-II.
Molecules 26 07310 g008
Figure 9. Error histogram analysis between targeted data and approximate solutions for scenario I.
Figure 9. Error histogram analysis between targeted data and approximate solutions for scenario I.
Molecules 26 07310 g009
Figure 10. Error histogram analysis between targeted data and approximate solutions for scenario II.
Figure 10. Error histogram analysis between targeted data and approximate solutions for scenario II.
Molecules 26 07310 g010
Figure 11. Performance of NNs–LM based on gradient, mu and validations failure during the process of optimization for different cases of Scenario I.
Figure 11. Performance of NNs–LM based on gradient, mu and validations failure during the process of optimization for different cases of Scenario I.
Molecules 26 07310 g011
Figure 12. Performance of NNs–LM based on gradient, mu and validations failure during the process of optimization for different cases of Scenario II.
Figure 12. Performance of NNs–LM based on gradient, mu and validations failure during the process of optimization for different cases of Scenario II.
Molecules 26 07310 g012
Figure 13. (af) shows the regressionanalysis for dimensionless concentration of substrate of micro-disk biosensor model.
Figure 13. (af) shows the regressionanalysis for dimensionless concentration of substrate of micro-disk biosensor model.
Molecules 26 07310 g013
Figure 14. (af) shows the regression analysis for dimensionless concentration of hydrogen peroxide of micro-disk biosensor model.
Figure 14. (af) shows the regression analysis for dimensionless concentration of hydrogen peroxide of micro-disk biosensor model.
Molecules 26 07310 g014
Figure 15. Plotof three dimensional (a) concentration profiles of substrate and (b) hydrogen peroxide against dimensionless distance R and reaction–diffusion parameters. Saturation parameter α is equal to 50.
Figure 15. Plotof three dimensional (a) concentration profiles of substrate and (b) hydrogen peroxide against dimensionless distance R and reaction–diffusion parameters. Saturation parameter α is equal to 50.
Molecules 26 07310 g015
Figure 16. Plot of three dimensional (a) concentration profiles of substrate and (b) hydrogen peroxide against dimensionless distance R and Saturation parameter α . Reaction–diffusion parameters are equal to 50.
Figure 16. Plot of three dimensional (a) concentration profiles of substrate and (b) hydrogen peroxide against dimensionless distance R and Saturation parameter α . Reaction–diffusion parameters are equal to 50.
Molecules 26 07310 g016
Figure 17. Normal probability curves for the computational complexity of different scenarios for (a) concentration profiles of substrate and (b) hydrogen peroxide of micro-disk biosensor.
Figure 17. Normal probability curves for the computational complexity of different scenarios for (a) concentration profiles of substrate and (b) hydrogen peroxide of micro-disk biosensor.
Molecules 26 07310 g017
Table 1. Comparison of approximate solutions obtained by NNs–LM algorithm with homotopy analysis method, modified Adomain decomposition method, hyperbolic function method and numerical solver Pdex4.
Table 1. Comparison of approximate solutions obtained by NNs–LM algorithm with homotopy analysis method, modified Adomain decomposition method, hyperbolic function method and numerical solver Pdex4.
α = 100 , γ E = 30 α = 50 , γ E = 20 α = 10 , γ E = 10
RHAMHFMMADMNumericalNNs–LMTHAMHFMMADMNumericalNNs–LMTHAMHFMMADMNumericalNNs–LMT
1.00.70380.70130.70250.70030.70030.61020.60000.61100.60100.60100.53160.53100.53200.53000.5300
1.50.71090.70740.70990.70750.70750.60140.60700.61210.60720.60720.54320.53970.53900.53980.5398
2.00.71120.71030.71080.71050.71050.61230.60810.61290.60800.60800.54730.54180.54540.54200.5420
2.50.71960.71590.71770.71580.71580.61580.61320.61450.61250.61250.54670.54740.54980.54760.5476
3.00.73320.73260.73290.73270.73270.62970.62820.62980.62800.62800.56780.56500.56980.56510.5651
3.50.75650.75900.76620.75950.75950.66710.66210.67010.66210.66210.60890.60350.60940.60370.6037
4.00.81320.81310.81320.81310.81310.72190.72890.73240.72890.72890.67800.67940.67680.67950.6795
4.50.88960.88940.88950.88940.88940.81900.82600.82700.82600.82600.79360.79050.78970.79050.7905
5.0111110.99670.99650.9963110.99990.99930.99980.99920.9992
Table 2. Statistical analysis of performance measures including MSE, Gradient, mu, number of iterations and time taken by the system for obtaining the results of Scenario I.
Table 2. Statistical analysis of performance measures including MSE, Gradient, mu, number of iterations and time taken by the system for obtaining the results of Scenario I.
Case ICase IICase III
S ( R ) H ( R ) S ( R ) H ( R ) S ( R ) H ( R )
Hidden Neurons606060606060
Training 3.02 × 10 13 5.44 × 10 12 1.91 × 10 13 2.76 × 10 11 9.21 × 10 14 8.53 × 10 12
Validation 3.95 × 10 13 6.43 × 10 12 2.53 × 10 13 4.56 × 10 11 1.38 × 10 13 6.16 × 10 12
Testing 5.29 × 10 13 7.67 × 10 12 2.68 × 10 13 4.76 × 10 11 3.11 × 10 13 5.12 × 10 13
Gradient 9.92 × 10 13 4.11 × 10 7 9.76 × 10 8 9.65 × 10 8 9.94 × 10 8 4.53 × 10 9
Mu 1.00 × 10 15 1.00 × 10 15 1.00 × 10 15 1.00 × 10 13 1.00 × 10 13 1.00 × 10 12
Epochs58166596378138
Regression111111
Time (s)<1 s<1 s<1 s<1 s<1 s<1 s
Table 3. Statistical analysis of performance measures including MSE, Gradient, mu, number of iterations and time taken by the system for obtaining the results of Scenario II.
Table 3. Statistical analysis of performance measures including MSE, Gradient, mu, number of iterations and time taken by the system for obtaining the results of Scenario II.
Case ICase IICase III
S ( R ) H ( R ) S ( R ) H ( R ) S ( R ) H ( R )
Hidden Neurons606060606060
Training 2.03 × 10 14 1.55 × 10 12 2.47 × 10 14 4.98 × 10 12 1.61 × 10 14 3.02 × 10 12
Validation 1.85 × 10 14 1.67 × 10 12 2.95 × 10 14 7.38 × 10 12 1.40 × 10 14 4.45 × 10 12
Testing 3.31 × 10 14 1.75 × 10 12 3.84 × 10 14 6.51 × 10 12 1.67 × 10 13 3.58 × 10 12
Gradient 9.91 × 10 8 8.39 × 10 8 9.76 × 10 8 8.98 × 10 8 2.59 × 10 9 9.74 × 10 8
Mu 1.00 × 10 16 1.00 × 10 17 1.00 × 10 16 1.00 × 10 15 1.00 × 10 15 1.00 × 10 15
Epochs442139184121
Time (s)<1 s<1 s<1 s<1 s<1 s<1 s
Table 4. Stability analysis on the values of performance function and performance indicators in terms of minimum, mean and standard deviations for different cases of scenario-I.
Table 4. Stability analysis on the values of performance function and performance indicators in terms of minimum, mean and standard deviations for different cases of scenario-I.
MSE MAD TIC ENSE
Min.MeanStd.Min.MeanStd.Min.MeanStd.Min.MeanStd.
Case I S ( R ) 3.95 × 10 13 3.75 × 10 12 4.86 × 10 12 6.54 × 10 6 1.23 × 10 5 4.42 × 10 5 1.96 × 10 7 3.78 × 10 7 1.76 × 10 7 7.23 × 10 14 2.81 × 10 13 1.72 × 10 13
H ( R ) 6.43 × 10 12 4.27 × 10 11 3.71 × 10 11 2.55 × 10 6 2.22 × 10 5 6.74 × 10 5 2.35 × 10 8 5.35 × 10 6 2.88 × 10 7 6.35 × 10 13 9.17 × 10 12 5.32 × 10 13
Case II S ( R ) 2.53 × 10 13 1.13 × 10 12 4.60 × 10 12 1.06 × 10 6 6.75 × 10 5 3.06 × 10 6 4.64 × 10 8 4.51 × 10 6 4.95 × 10 6 5.05 × 10 14 4.97 × 10 13 2.41 × 10 13
H ( R ) 4.56 × 10 11 3.81 × 10 10 5.89 × 10 9 6.92 × 10 6 4.40 × 10 5 3.18 × 10 5 6.23 × 10 7 1.36 × 10 8 6.25 × 10 7 1.18 × 10 13 1.67 × 10 12 2.10 × 10 13
Case III S ( R ) 1.38 × 10 13 6.51 × 10 12 1.07 × 10 12 2.08 × 10 7 6.17 × 10 6 2.20 × 10 6 4.25 × 10 7 4.51 × 10 6 3.35 × 10 6 2.09 × 10 15 7.95 × 10 14 1.59 × 10 13
H ( R ) 6.16 × 10 12 8.96 × 10 11 1.23 × 10 11 7.56 × 10 6 8.64 × 10 5 3.38 × 10 6 1.49 × 10 7 2.19 × 10 6 2.49 × 10 7 2.75 × 10 13 2.95 × 10 12 9.41 × 10 12
Table 5. Stability analysis on the values of performance function and performance indicators in terms of minimum, mean and standard deviations for different cases of scenario-II.
Table 5. Stability analysis on the values of performance function and performance indicators in terms of minimum, mean and standard deviations for different cases of scenario-II.
MSE MAD TIC ENSE
Min.MeanStd.Min.MeanStd.Min.MeanStd.Min.MeanStd.
Case I S ( R ) 1.84 × 10 14 3.40 × 10 12 5.87 × 10 11 1.37 × 10 6 9.08 × 10 6 2.12 × 10 6 1.97 × 10 7 1.49 × 10 6 7.33 × 10 7 3.10 × 10 14 2.76 × 10 13 2.13 × 10 14
H ( R ) 1.66 × 10 12 1.14 × 10 11 1.98 × 10 11 2.58 × 10 6 2.83 × 10 5 2.68 × 10 6 2.15 × 10 7 1.28 × 10 6 2.06 × 10 7 9.32 × 10 14 1.97 × 10 13 2.75 × 10 13
Case II S ( R ) 2.95 × 10 14 2.75 × 10 13 3.16 × 10 13 3.13 × 10 6 2.66 × 10 5 1.82 × 10 6 7.44 × 10 8 4.11 × 10 6 1.46 × 10 6 2.24 × 10 15 2.63 × 10 14 2.56 × 10 14
H ( R ) 7.38 × 10 12 9.12 × 10 11 1.80 × 10 11 2.34 × 10 6 2.44 × 10 5 2.08 × 10 6 1.30 × 10 7 7.59 × 10 7 8.68 × 10 7 1.87 × 10 13 2.55 × 10 12 2.79 × 10 13
Case III S ( R ) 1.40 × 10 14 7.43 × 10 13 4.49 × 10 13 1.56 × 10 6 2.37 × 10 5 2.72 × 10 6 2.56 × 10 7 1.90 × 10 6 8.85 × 10 7 2.85 × 10 14 1.37 × 10 13 2.14 × 10 13
H ( R ) 4.45 × 10 12 2.69 × 10 11 1.03 × 10 11 2.23 × 10 6 3.10 × 10 5 3.44 × 10 6 3.19 × 10 7 2.38 × 10 6 1.15 × 10 6 2.73 × 10 13 1.61 × 10 12 2.61 × 10 12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khan, N.A.; Alshammari, F.S.; Romero, C.A.T.; Sulaiman, M.; Laouini, G. Mathematical Analysis of Reaction–Diffusion Equations Modeling the Michaelis–Menten Kinetics in a Micro-Disk Biosensor. Molecules 2021, 26, 7310. https://doi.org/10.3390/molecules26237310

AMA Style

Khan NA, Alshammari FS, Romero CAT, Sulaiman M, Laouini G. Mathematical Analysis of Reaction–Diffusion Equations Modeling the Michaelis–Menten Kinetics in a Micro-Disk Biosensor. Molecules. 2021; 26(23):7310. https://doi.org/10.3390/molecules26237310

Chicago/Turabian Style

Khan, Naveed Ahmad, Fahad Sameer Alshammari, Carlos Andrés Tavera Romero, Muhammad Sulaiman, and Ghaylen Laouini. 2021. "Mathematical Analysis of Reaction–Diffusion Equations Modeling the Michaelis–Menten Kinetics in a Micro-Disk Biosensor" Molecules 26, no. 23: 7310. https://doi.org/10.3390/molecules26237310

Article Metrics

Back to TopTop