Mathematical Study of Medicine Propagation in Biological Tissue and Some of Its Applications

The paper deals with the problem of the distribution of the medicine (en-zyme) in the damaged biological tissue where the reaction enzyme—substrat takes place. The biological problem is reduced to a singular degenerate ini-tial-boundary value problem for two coupled ordinary differential equations. Analytical solution of the singular degenerated IBV-problem was obtained by power series. The solution demonstrates the real situation and found suitable to depict the degeneration of singular system, caused by low concentration of the enzyme.


Introduction and Auxiliary Results
There are many different works dealing with medicine (enzyme) transportation to the damaged zone. The mathematical modeling of biological processes often leads to highly complex systems involving many state-variables and reactions. For a wide range of analysis model complexity can present an insurmountable barrier [1]. Methods of model reduction therefore remain a vital topic and a widely applicable tool in the analysis and modeling of biochemical system. To describe a biochemical reaction of medicine, injected into the tissue, we present here mathematical model [2]. The model fully reflects the formation of the enzyme-substrate complex in the damaged tissue [3] [4]. It makes possible to demonstrate the real situation and found suitable to depict the degeneration of singular system, caused by low concentration of the enzyme. The same mathemati-cal model was considered also in the previous work [5]. The analytical solution of the corresponding Cauchy problem using the power series method and numerical solutions using the Runge-Kutta-Felberg method [6] were obtained for a large time scales. Some other authors [7] discus changes in substrate concentration (decreases with time) and enzyme-extract complex (increases with time). On an example the natural rubber bio synthesis process using the Michaelis-Menten model the quasi stationary states with constant rates of intermediate reactions in a chain were obtained. Conditions of convergence are developed for different situations in the case of a biochemical reaction were examined [8]. In other work "Extensions to Michaelis-Menten Kinetics for Single Parameters" [9] they defined the equation for substrate and product and rate of reaction based on rate and time perspectives. Two quadratic equations were developed. The first, represents a body entity that gave a useful relationship of enzyme and the second of nutrients. By combining Enzymes and Enzyme-substrate complexes values, they derive a Michaelis-Menten hyperbolic function and the proportionate rate of reaction and Enzymes values of the quadratics resulting in another Michaelis-Menten hyperbolic. It was obtained that between these two hyperbolic functions, in-competitive inhibitions exist, indicating metabolic activities and growth in terms of energy levels. The effect of decomposition to the kinetics of the enzyme reaction based on a microscopic pharmacokinetic model and the application of the Monte Carlo method to analyze that explains possible reaction scenarios in disordered media using the kinetics of the Michaelis-Menten model in pharmacology [10]. In [11] complex variables method combined with Monte Carlo method.
It this paper we provide the overview of the main methods proposed for an analytical solution. The solution of the Cauchy problem with initial conditions is given for two functions ( ) x t and ( ) y t , which is characterized by a singular degeneracy with respect to the parameter  (Section 2). An analytical method for representing the solution by power series was proposed. The functions in power series up to the fifth order were calculated (in Section 3). On this basis, solutions for changes in the concentration of the substrate and enzyme-substrate complex were obtained (in Section 4). Accordingly, when using the Runge-Kutta-Felberg method, calculations were carried out up to the 9 th order. In the above analysis, in contrast to previous studies, it is not required to take into account the limitations on the very small amount of the introduced enzyme  .

The Problem Definition
The mathematical model can be represented as the Cauchy problem for a system of coupled ordinary differential equations for two functions ( ) with the following initial conditions

Analytical Solution for Arbitrary Functions
Equations (2.1) and (2.2) are strongly coupled with nonlinear terms. So, to construct an analytical solutions for the Cauchy problem is not trivial. In this work, we use the Michaelis-Menten method for a general case with no restrictions on the initial enzyme concentration. For solution of the corresponding Cauchy problem we use a power series method, applying to the time coordinate. The method is valid for a short interactions, wihout any restrictions on the enzime concentration. Now, more relative terms were held in the expansion, than in the in series in earlier works [5]. The second approach is the numerical solution of the problem by the Runge-Kutta-Felberg method, which is valid for long time intervals and large values of parameters and concentrations, i.e., for strong interactions when the values of the enzyme and enzyme-substrate complex are of the same order. We present the results in dimensionless form. The results are suitable for interactions in the case of any real enzyme and substrate structures. Usually, the medicine transported for a sufficiently large distance, the concentration of the medicine is significantly decreases. Consequently, the parameter is characterized by the initial concentration of the enzyme for which the initial reaction is very weak. In this case, problem (2.1) and (2.3) is singularly degenerate, since in this case the coefficient of derivative is very small, and the concentration of the enzyme is low, and an approximate solution to the singularly degenerate problem can be obtained by the method of coalescence asymptotic expansions [4]. In the case when the medicine (enzyme) is distributed in the immediate Near the point where the concentration is large and no singular degeneration appears, we consider the problem (2.1) and (2.3) for arbitrary parameters. Then the problem contains strongly nonlinear systems of differential Equations (2.1) and (2.2) and becomes significantly more complicated. In this case, the concentrations of the substrate and enzyme are of the same order 0 0 s e ≅ . The order ( ) and initial reaction rate can be calculated using the corrected formula:  ( ) ( ) 2 According to (3.4), (3.5) we preserve more terms than (3.7) and (3.8) (k = 3.5... 9). Expressions are also obtained for the coefficients 3 4 5 , , a a a and

Analysis of Results
In Figure 1, the results of the numerical Runge-Kutta-Felberg method and the power series approach are presented and compared. According to the results, we conclude that for relatively short time periods, the series method gives very good approximation, and may be used for descriptions of the expected medicine distribution, characterize the features of its transportation, subsequent localization and its possible impact. From Figure 2, it can be seen that for smaller  value provides maximum concentration of the enzyme substrate complex. The maximum values increase with decreasing of  .