Existence of fractional order semianalytical results for enzyme kinetics model

*Correspondence: hussam89@yahoo.com 1College of Engineering, Al Ain University, Al Ain, United Arab Emirates 2Department of Mathematics, Tafila Technical University, Tafila, Jordan Full list of author information is available at the end of the article Abstract This research work investigates the existence of semianalytical solutions of a chemical kinematics model. We develop the conditions for the existence of the solutions for a proposed enzyme kinetics model, via tools of the fixed point theory. The semianalytical results were obtained with the help of Laplace transformation and Adomian decomposition method. The results established by the proposed techniques are in the form of infinite series. Furthermore, with extending homotopy perturbation method (HPM), we develop a series solutions for the considered model. By using Matlab, we present the approximate solution for both methods up to a few series terms.


Introduction
Recently, experimental evidence shows that dynamics problems in nature follow a fractional calculus analysis. The relevant field of research is a fast growing area, due to its numerous applications in diverse and widespread fields of engineering and science; such as chemical models, physics, signal and image processing, quantum mechanics, control theory, nonlinear dynamics, biological population models, optimization theory, and much more [1][2][3][4][5][6][7][8][9][10]. Instantly, it is evident that dealing with a dynamical system with memory effects is one of the biggest challenges for researchers. Fractional calculus has a direct link to dynamical systems (with memory effect). Therefore, fractional differential equations (FDEs) present a novel technique developed to model phenomena related to the dynamics of the aforesaid fields of science [11][12][13][14]. Fractional derivatives are global in nature and offer a greater degree of freedom compared to the conventional derivatives. Numerous researchers have investigated various features of FDEs concerning the existence, stability analysis, and approximate solutions. They utilized different techniques of fixed-point theory and numerical analysis to investigate the existence theory, stability analysis, and approximate solutions of FDEs [15][16][17][18].
An important aspect of the this concerned field is mathematical modeling, which is a tool that describes almost all dynamical phenomena in mathematical language and con-cepts. Mathematical modeling is an effective technique that produces some of the best results for optimization and quantitative productivity (see [19,20]). The mathematical approach for enzymes' kinetics has been discussed in different articles, and also different mathematical techniques have been used for analytical solution of proposed models of chemical reactions (see [19,[21][22][23][24][25]). In this regard, Sharpe and Lotka in 1923 studied a kinematic enzyme in mathematical form for ecology and epidemiology. Along the same lines, L. Michaelis and M. Menton [23] proposed a mathematical model for a basic reaction of an enzyme given by Studying the chemical reaction (1), one can infer that when enzyme E of one molecule is combined with substrate S of one molecule we will obtain a complex SE of one molecule.
Here η 1 denotes the rate of formation of enzyme, η 2 is the rate of creation of the product. Assume that the concentrations of S, E, SE are equal to S, E, C, respectively, i.e., where the concentration of a substance is represented by square brackets [ ]. The author in [26] used the law of mass action to derive the following four nonlinear differential equations:Ṡ In addition, one supplements the conditions S(0) = S 0 , C(0) = 0, E(0) = E 0 , and P(0) = 0.
Here η 1 , η 2 , η 3 , S 0 , and E 0 are positive constants. Alicea [24] solved the system of chemical reaction (2) by using multiple time scales method asymptotically. Alawneh [27] solved model (2) by a generalized differential transform method with multisteps and discussed the estimation analysis via fractional order derivative. Using the physical interpretation of the solution of system (2) with concentrations S, E, C, and P at any time τ obtained mathematically, one observes that most of the time it does not coincide with the results obtained experimentally. In order to overcome such deficiency, one way is to use another type of differential operator instead of integer order derivative, which is usually known as fractional derivative. The researchers paid considerable attention to the aforementioned class of derivatives because they are more flexible and accurate compared with the classical derivatives, we refer the interested readers to [4,28]. Keeping these applications of fractional differential operator in mind, we introduce the following noninteger order enzyme kinematic system: subject to the initial conditions S(0) = n 1 , C(0) = n 3 , E(0) = n 2 , P(0) = n 4 , where c D μ is noninteger order Caputo derivative and 0 < μ < 1. Further, n 1 , n 2 > 0 and n 3 ≥ 0, n 4 ≥ 0. In the proposed system (3), the supplied conditions are independent of each other and also satisfy the relation M(τ ) = S(τ ) + C(τ ) + E(τ ) + P(τ ). The value of M is used to represent terms presented in the concerned system of reactions.
Corresponding to model (3), we use the tools of fixed point theory to investigate some results that ensure the existence of the aforementioned model and its solution. We utilize Banach and Schauder's theorems from the proposed existence theory. We obtained the estimated solution of the model of noninteger order via Laplace transform combined with Adomian decomposition method which is known as Laplace Adomian Decomposition Method (LADM). LADM is an efficient technique by which we can find both explicit and analytic solutions for the system of differential equations. Those techniques are efficient and work outstandingly in both cases, i.e, in initial and boundary value problems. This method also works accurately for a system of stochastic differential equations. LADM does not need liberalization or perturbation like other existing computational and analytical schemes that need exploring the dynamical behavior of complex dynamical systems. The adopted techniques provide significant results for the solutions of FODEs, as well as for analytical solutions for the variety of problems of nonlinear differential equations.
In this paper, we utilize techniques of Adomian polynomials to decompose the nonlinearity and Laplace transform to convert the problem in hand into the form of algebraic equations, see [29]. Recently, the proposed techniques have been used to deal with nonsingular FODEs to obtain a very fruitful results (see [30]). To justify the results obtained via proposed techniques, we use Maple-13 and assign different values for the parameters and the supplied conditions. Furthermore, we remark that the obtained results via this method are in a form of convergent series, converging to the exact result uniformly. Thanks to the results of analysis obtained in [31][32][33], one can easily prove the convergence of the proposed method. Also, we construct the He homotopy perturbation method (HPM) to compute the concerned solutions for the proposed model. We compare both solutions and this HPM method worked very well (see [34,35]).

Preliminaries
To make things easier, in this section we provide some well-known definitions, theorems, and lemmas. Their details are available in [4,29,30].

Definition 1
The fractional integral in the sense of Riemann-Liouville of an arbitrary order μ for f ∈ L 1 ([0, ∞), R) is given as such that the integral on the right-hand side is defined on (0, ∞).

Definition 2 The Caputo fractional order derivative of a function f is given by
the integral part of μ is represented by [μ] and n = [μ] + 1. The concerned derivative is used through this work.

Lemma 1
In case of fractional differential equations, the following result holds: where ξ k is any real number, k is a positive integer up to n -1 and the integral part of μ is denoted by [μ], while n = [μ] + 1.

Definition 3
The Laplace transform in the sense of proposed derivative (Caputo) as defined by

Existence theory of enzyme kinetics model (3)
This section is devoted to determining existence results for the proposed problem (3). In this connection, we define the following functions: Further, we set B = C[0, T] to be a Banach space with Keeping in mind the above notation, system (3) may be expressed as In view of I μ on both sides of (5), we have Let us define A : B → B to be an operator by using (6) as Further, we assume that the following hypothesis holds: (T 2 ) There exists a constant L H > 0 such that, for every X, X ∈ B, We established the results for at least one solution by utilizing Leray-Schauder fixed point theory.
This means that A(X) B ≤ ρ. Hence A(D) ⊂ D. Also A is continuous. Letting τ 1 < τ 2 ∈ [0, T], we show that A is a completely continuous operator. For this, we have Now in (8), the right-hand side approaches zero when τ 2 → τ 1 . Since A is bounded, Thus by Arzelá-Ascoli theorem, A is relatively compact. Hence A is completely continuous. Thus via Schauder theorem, system (3) has a solution.
To ensure uniqueness of solution of (3), we have the next result.
Proof Here we state that the first equation in this proof is given in (6). Further, if X, X ∈ B, A : B → B is the operator defined above. Consider Hence, the operator A is a contraction. Thus the system has a unique solution.
In the same way, the other terms can be computed too. We assigned random values to μ, in order to observe the mathematical dynamics of the above results (S(τ ), E(τ ), C(τ ), p(τ )).

Conclusions
The presented work is devoted to a fractional-order model of enzyme kinetics, which plays a role of a catalyst in various reactions taking place in living organisms. By means of nonlinear analysis, we have derived some results about the existence of at least one solution for the model under consideration. With the help of LADM, we have developed the semianalytic results for the proposed noninteger order model (3). The presented approximate results have been obtained by using famous He's HPM. Both methods have resulted the same solution. We have plotted the proposed solutions for both methods up to the first three terms, which completely agree with each other. Both methods can be used as a powerful mathematical tools to deal with many nonlinear problems of FDEs.