Next Article in Journal
Inertial Projection Algorithm for Solving Split Best Proximity Point and Mixed Equilibrium Problems in Hilbert Spaces
Next Article in Special Issue
An MDL-Based Wavelet Scattering Features Selection for Signal Classification
Previous Article in Journal
On Newton’s First Law of Motion
Previous Article in Special Issue
An Adaptive Three-Dimensional Improved Virtual Force Coverage Algorithm for Nodes in WSN
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultimate Bounds for a Diabetes Mathematical Model Considering Glucose Homeostasis

Postgraduate Program in Engineering Sciences, BioMath Research Group, Tecnológico Nacional de México/IT Tijuana, Blvd. Alberto Limón Padilla s/n, Tijuana 22454, Mexico
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(7), 320; https://doi.org/10.3390/axioms11070320
Submission received: 27 May 2022 / Revised: 20 June 2022 / Accepted: 28 June 2022 / Published: 30 June 2022
(This article belongs to the Special Issue Mathematical and Computational Applications)

Abstract

:
This paper deals with a recently reported mathematical model formulated by five first-order ordinary differential equations that describe glucoregulatory dynamics. As main contributions, we found a localization domain with all compact invariant sets; we settled on sufficient conditions for the existence of a bounded positively-invariant domain. We applied the localization of compact invariant sets and Lyapunov’s direct methods to obtain these results. The localization results establish the maximum cell concentration for each variable. On the other hand, Lyapunov’s direct method provides sufficient conditions for the bounded positively-invariant domain to attract all trajectories with non-negative initial conditions. Further, we illustrate our analytical results with numerical simulations. Overall, our results are valuable information for a better understanding of this disease. Bounds and attractive domains are crucial tools to design practical applications such as insulin controllers or in silico experiments. In addition, the model can be used to understand the long-term dynamics of the system.
MSC:
34A34; 34D23; 34D45

1. Introduction

With high rates of mortality and morbidity, diabetes mellitus is considered a worldwide public health problem. Diabetes complications killed around 4.2 million adults between 20 and 79 years in 2019, which represents about 11% of deaths from all causes. On the other hand, 1.9 million of the diabetes mortalities (46%) happen in adults aged 60 or less [1,2]. There are two categories of diabetes mellitus: (i) Type-1 (T1DM) is considered an autoimmune disease that destroys insulin-producing pancreatic β -cells; it is normally juvenile-onset and commonly known as insulin-dependent diabetes. (ii) Type-2 (T2DM), develops usually in middle-aged or older people with insulin resistance; also called adult-onset diabetes.

1.1. Mathematical Modeling of Diabetes Mellitus

Diabetes mellitus is an open problem [3]. We have to understand the corresponding interactions and contributions in the parthenogenesis of this disease with respect to various defects of the glucose-insulin regulatory system associated with β -cells mass, the responsiveness level of β -cells to glucose, and the sensitivity of tissues to insulin. Regarding mathematical models based on T1DM and T2DM dynamical sense behavior, a realistic representation of the long-term physiological adaptation in developing insulin resistance is necessary for effectively design clinical trials and evaluate diabetes prevention or disease modification therapies [4,5]. During the past century, numerous mathematical models of the glucoregulatory system have been formulated or proposed in the field of diabetes with different purposes [6,7]. One of the main characteristics of these in silico models is dependence on the variability of their parameters. However, the main difficulty is dealing with many parameters, and complex metabolic studies are necessary to estimate them. These difficulties compel researchers to use approximations of the parameters’ values, decreasing the validity of models and their applicability. However, biomedical researchers are skeptical about utilizing mathematical modeling to investigate human diseases, so in silico experimentation is a powerful tool to reach this community [8]. A trustworthy understanding of diabetes progression is well known thanks to descriptive research by cross-sectional studies and rational extrapolation techniques. Nevertheless, the lack of a mathematical model interpretation by ordinary differential equations (ODEs) remains a research challenge [9]. In [10], authors provided a summary of how eight diabetes models match up to data from published clinical studies, to highlight differences between models, and to offer insights into the challenges facing diabetes models. They conclude that the used models predict well in typical cases, but fail with complicated scenarios, such as T2DM patients. In this direction, mathematical modeling gives us a reasonably precise means of approximating the phenomenon’s overall structure.

1.2. Recent Contributions, ODE Model Description, and Proposed Analysis

The Localization of Compact Invariant Sets (LCIS) is an analytical method applied to mathematical models described by first-order ODEs to provide bounding conditions for the state variables of a system; see [11,12]. Furthermore, related to diabetes, the LCIS method provides a broad understanding of β -cells behavior in the presence of glucose [13] or with the immunological response [14]. In this paper, we are interested in analyzing through the LCIS method an ODE mathematical model obtained by in vitro and in silico estimation. Therefore, results give us the solution using the considered variables over time based on upper bounds computation. The mathematical model we use was presented in [15]. The numerical values of parameters were estimated by in vivo experimentation; thus, background values were validated in vivo, showing a good approximation to the time-series of the real measurements of insulin levels and plasma glucose. Although this model was initially introduced in [16], it was used as a reference to evaluate the time series of glucose homeostasis in other kinds of studies [17,18]. Modeling, analysis, and in silico experimentation for diabetes mellitus are valuable tools that allows better understanding of the complicated nature of this disease. In particular, recently reported work in this direction offers valuable information with respect to specific aspects of the development of this disease. A strictly analytical mathematical tool, such as the method of LCIS, provides a way to understand the dynamical properties of a biological problem, revealing the non-intuitive properties of the system. However, Lyapunov’s theory allows us to determine sufficient conditions for the existence of the bounded positively invariant domain (BPID), which is an asymptotically attractive region including all the significant dynamics of the system; see Lemma 4.1 by Khalil at Section 4.2 in [19].
The organization of this paper is the following. The next section describes the system formulated by five first-order ODEs. Following that, the LCIS method is briefly discussed. Section 3 is devoted to applying the LCIS method to the system. Further, in this section we show that Lyapunov’s direct method allows calculating the BPID. Then, numerical simulations illustrate our results. Finally, the last section presents the conclusions of this research.

2. Materials and Methods

The mathematical model proposed in [15] is used as a baseline to study the T1DM and a T2DM. A previous version of this model was published in [16], where the authors aimed to create a very good glucoregulatory system description without involving a large number of mathematical equations; see [20,21]. That model consists of three ODEs with eight parameters describing the dynamic concentrations of blood insulin, blood glucose concentration, and glucose in the intestine. However, specialized control algorithms such as HOMA Calculator [22] and UVA/Padova [23] are included in the mathematical model under study in this work, and it is given as follows:
d I ( t ) d t = k 6 I + k 1 G + k 7 Y ,
d Y ( t ) d t = k 8 Y ,
d D ( t ) d t = k a D ,
d G ( t ) d t = k 4 ( I I p i ) k 3 k 2 I + k 0 D k 5 ( G G u ) H ,
d U ( t ) d t = k 5 ( G G u ) H ,
where I ( t ) represents the plasma insulin variation (pmol/l), Y ( t ) is the amount of insulin in the subcutaneous compartment (UI), D ( t ) is the amount of glucose in the digestive system (mg), G ( t ) is the related plasma glucose level (mg/dl), and U ( t ) is associated with the amount of glucose in the urine (mg). Equations (4) and (5) include a term H that is a function of glycaemia given in terms of G, as H ( G ) , given by Equation (6).
H G = 0 i f G G u , 1 i f G > G u .
For the sake of the reader’s convenience, we include Figure 1 to interpret the contribution of each differential equation and its relationship with the parameters. The descriptions of the model and its parameters are presented in [15]; see Table 1.

The Localization of Compact Invariant Sets Method

The LCIS method proposed by Krishchenko and Starkov in [24,25] aims to understand the long-time dynamics of nonlinear systems formulated by autonomous and non-autonomous first-order ODEs. The later is achieved by computing both lower and upper bounds of the so-called localizing domain, which is a region in the state space R n where all compact invariant sets of the system are located. Equilibrium points, periodic orbits, limit cycles, and chaotic attractors are examples of compact invariant sets. Localizing functions are selected by a heuristic process. This means that one may need to analyze several functions in order to find a proper set that will allow one to fulfill Theorem A1. The LCIS method is deeply described in Appendix A.

3. Results

3.1. Bounding the Diabetes Mathematical Model including Glucose Homeostasis Using LCIS Method

In this work, we apply the LCIS method to find the ultimate bounds of the system (1)–(6). Recently, the behavior of that model was studied in [15]. Moreover, the authors of that work claim that a deep understanding of the mathematical model can show how to reduce the flaws in strategies based on glucose control that use an insulin infusion pump. We now present the LCIS analysis to find the ultimate bounds of the system, which will allow us to have better background knowledge about the qualitative properties of the model. All parameters of the system are positive; their descriptions and values are given in Table 1; one can consult [15] for additional information about the model (1)–(5). Due to the biological meaning of the system, trajectories are confined to the non-negative orthant:
R + , 0 5 = I ( t ) 0 , Y ( t ) 0 , G ( t ) 0 , D ( t ) 0 , U ( t ) 0 ,
which is possible if the next condition
k 4 I p i > k 3
is fulfilled, according to [16]. Clearly, H ( G ) is a discontinuous function, which prevents analyzing the system with standard tools. We avoid that restriction by using instead of the Heaviside function (6), a sigmoid function given by Equation (9). Sigmoid monotonically increases, is continuous everywhere, and usually is interchangeable with the Heaviside function in biological modeling applications using the Hausdorff metric. Here we propose
H 1 G = ε + 1 1 + e c 1 G G u ,
where c 1 is a positive constant that allows changing the slope of the function (Hausdorff metric), and ε is an infinitesimal positive constant. Then,
ε H 1 G 1 + ε .
Now we define H 1 min : = ε and H 1 max : = 1 + ε . These bounds of the H 1 G function allow us to take the constant H 1 max to find the upper bounds by applying the LCIS method. From this point we define H 1 : = H 1 G H G . Now we apply the first localizing function by considering a free positive parameter q 1 :
h 1 = G + q 1 D .
The Lie derivative of (11) with respect to the system (1)–(5) is
L f h 1 = k 0 D k 2 I k 3 k 4 ( I I p i ) k 5 ( G G u ) H 1 G q 1 k a D
and the set S ( h 1 ) gives us the following mathematical expression:
h 1 S ( h 1 ) = k 5 1 H 1 max 1 k 5 G u + k 4 I p i k 3 q 1 k a k 0 q 1 k 5 H 1 max D ( k 2 + k 4 ) I .
Condition (8) must be fulfilled in order to take the following inequality:
k a > k 5 ,
and our free parameter q 1 is given by
q 1 k 0 k a k 5 H 1 max .
Therefore, we can ensure an upper bound for the function h 1 if conditions (8), (12), and (13) are satisfied, hence,
h 1 | S ( h 1 ) k 4 I p i k 3 + k 5 G u k 5 H 1 max ,
implying that set K 1 is defined as follows:
K 1 = G + q 1 D η 1 : = k 4 I p i k 3 + k 5 G u k 5 H 1 max .
Thus, for the sake of notations, we define G max : = η 1 and D max : = η 1 q 1 .
On the other hand, the bound of U ( t ) was mentioned in [16] with the following description: “when the plasma glucose level G G u , then H ( G ) = 0 , and there is not renal glucose excretion; i.e., U ( t ) = 0 . On the other hand, when the plasma glucose level G > G u , then H ( G ) = 1 , and the renal glucose excretion U ( t ) depends on ( G G u ) ." Then, taking our bound G max and considering the sigmoid function (9), we have the bound given by U max { k 5 G max ( 1 + ε ) } .
Now, let us propose a second localizing function with a free positive parameter q 2 :
h 2 = I + q 2 Y ,
where the Lie derivative of (11) with respect to the system (1)–(5) is
L f h 2 = k 1 G k 6 I + k 7 Y q 2 k 8 Y ,
leading to
h 2 | S ( h 2 ) = k 6 1 k 1 G + k 7 Y + q 2 ( k 6 k 8 ) Y ,
when fulfilling the following condition:
k 6 < k 8 .
The free parameter of q 2 must fulfill the next condition:
q 2 k 7 k 8 k 6 .
When applying the Iterative theorem to Equation (18) by considering the set K 1 through the upper bound of G m a x , we have
h 2 | S ( h 2 ) K 1 k 1 k 6 G max ,
implying that
K 2 K 1 = I + q 2 Y h 2 | S ( h 2 ) K 1 : = η 2 = k 1 k 6 G max .
Then,
K 2 = I + q 2 Y η 2 ,
with I max = η 2 and Y max = η 2 q 2 . Hence, we derive the following theorem.
Theorem 1.
Suppose that conditions (8), (12), (13), (19), and (20) are fulfilled for a set of the parameters of the system (1)–(5). Let K 1 and K 2 be given by Equations (15) and (23), respectively. The set K ( h ) defined by Equation (24) contains all compact invariant sets for the model (1)–(5), with
K h = K 1 K 2 .
For the parameters given in Table 1, all the conditions (8), (12), (13), (19), and (20) are fulfilled.

3.2. Bounded Positively Invariant Domain

Now, let us consider a Lyapunov candidate function to define a BPID, ensuring that all trajectories of the system (1)–(5) that go inside of BPID will never leave this domain. In Equation (25) is proposed a function given the mathematical domain of the variables R + , 0 5 . Free positive parameters q 4 and q 5 are associated.
V = I + q 4 Y + D + q 5 G + U ,
Then, the Lyapunov derivative is the following expression when β 1 = k 6 + q 5 ( k 4 + k 2 ) and β 2 = q 5 ( k 4 I p i k 3 + k 5 G u H 1 m a x ) k 5 G u H 1 m a x .
V ˙ = β 1 I ( q 5 k 5 H 1 m a x k 5 H 1 m a x k 1 ) G ( q 4 k 8 k 7 ) Y ( k a q 5 k 0 ) D + β 2 .
Given function V ˙ in (26), and if the following conditions are fulfilled
q 4 > k 6 k 8 ,
k 5 H 1 m a x + k 1 k 5 H 1 m a x < q 5 < k a k 0 ,
k 4 I p i + k 5 G u > k 3 ,
then with V ˙ = 0 , the polytope (30)
β 1 I + ( k 5 H 1 m a x ( q 5 1 ) k 1 ) G + ( q 4 k 8 k 7 ) Y + ( k a q 5 k 0 ) D β 2
is obtained. From the expression (30), we can get that V ˙ < 0 outside this polytope. Therefore, for any ρ 1 > 0 , the solid polytope (31)
Ξ 1 ( ρ 1 ) = β 1 I + ( k 5 H 1 m a x ( q 5 1 ) k 1 ) G + ( q 4 k 8 k 7 ) Y + ( k a q 5 k 0 ) D β 2 + ρ 1
is a positive invariant domain, i.e., each trajectory of the system (1)–(5) entering into Ξ 1 ( ρ 1 ) remains there; see Lemma 4.1 by Khalil in [19]
Furthermore, every compact invariant set is contained in Ξ 1 ( ρ 1 ) . This result can be combined with K ( h ) (24) to sharpen the BPID.

3.3. Analysis for the Model When a T2DM Case Is Considered

The model can be considered for a T2DM case if Y ( t ) = 0 (where insulin into the subcutaneous compartment is not considered). Localization K ( h ) of Equation (24) remains valid, so I max = η 2 , G max = η 1 , and D max = η 1 q 1 ; and condition (20) must be eliminated.
On the other hand, the BPID can be calculated in a similar way, and is given by
Ξ 2 ( ρ 2 ) = β 1 I + ( k 5 H 1 m a x ( q 5 1 ) k 1 ) G + ( k a q 5 k 0 ) D β 2 + ρ 2 ,
with β 1 and β 2 defined before. That means that for any ρ 2 > 0 , the solid polytope (32) is a positive invariant domain; i.e., each trajectory of the system (1)–(5) entering into Ξ 2 ( ρ 2 ) remains there; see Lemma 4.1 by Khalil in [19].

3.4. Numerical Simulations: In Silico Experimentation

This section includes numerical simulations as in silico experiments. Figures provide a visual way to interpret our results. We took the parameters of Table 1 to solve the system (1)–(5) with the Runge–Kutta method, using initial conditions outside of the BPID to show the attractiveness of that region. Here, Figure 2 and Figure 3 represent the simulation for the T1DM case. Moreover, Figure 4 and Figure 5 depict the curve when the T2DM case is considered. The trajectory goes inside the upper bounds G max and I max defined in (24), meaning that trajectories are going to the BPID. If the trajectories of Y ( t ) and D ( t ) are trivial (exponential type), then we do not show them here.

4. Discussion and Conclusions

This work focused on examining the ultimate bounds for a diabetes mathematical model, including glucose homeostasis. We applied the so-called method localization of compact invariant sets, which allowed us to determine the upper bounds for the system (1)–(5). On the other hand, we used a Lyapunov function to find a bounded positive invariant domain where the trajectories of the system are attracted. The original system contains a Heaviside function that limits the application of the LCIS method. We proposed a smooth function given by a sigmoid, which allowed us to solve the LCIS and BPID problems. The I m a x , Y m a x , D m a x , and G m a x bounds were computed, and there are valid domains for both the T1DM and T2DM cases. These bounds were constructed with inequalities containing parameters of the system. Application of the localization theorem gave us some inequalities given by (8), (12), (13), (19), and (20), which are necessary assumptions for the existence of localization sets K 1 and K 2 given by Equations (15) and (23), respectively.
With respect to BPID, we obtained a solid polytope Ξ 1 ( ρ 1 ) , which attracts trajectories of the system with conditions (27)–(29). In both cases, (LCIS and BPID) conditions were fulfilled for the parameters of Table 1, and that could be used to apply our bounds in specific control problems. We performed numerical simulations to show the most relevant bounds I m a x and G m a x .
Regarding biological implications, we found a BPID that depends directly on the parameters k 5 and k 6 , which are related to insulin, plasma glucose level, and glucose in urine. If we can change at least one parameter externally, it is possible to have a smaller BPID. That implies in some cases, better health conditions for a patient. This kind of application is possible with a controller or observer design but is necessary to evaluate hardware and medical implications that are out of the reach of this paper.
Future work will be focused on the development of observers and controllers based on our bounds, which allows solving several problems, i.e., robustness of the closed-loop system, Lyapunov based controllers or observers.

Author Contributions

Conceptualization, L.N.C. and D.G.; methodology, P.A.V. and L.N.C.; software, D.G. and L.N.C.; validation, P.A.V. and L.N.C.; formal analysis, L.N.C.; investigation, D.G. and L.N.C.; resources, P.A.V. and D.G.; writing—original draft preparation, L.N.C. and D.G.; writing—review and editing, L.N.C. and P.A.V.; visualization, P.A.V. and L.N.C.; project administration, D.G. and P.A.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by Tecnologico Nacional de Mexico/IT Tijuana grants 13886.22-P Análisis y estudio matemático no lineal aplicado en modelos biológicos asociados a Diabetes Mellitus Tipo 1; and 13747.22-P Modelizado computacional y experimentos in silico aplicados al análisis y control de sistemas biológicos.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

We thank the anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Appendix A

Localization of the Compact Invariant Sets Method The LCIS method is formulated as follows. Take a nonlinear system with the form x ˙ = f x , where f x is a C differentiable vector function and x R n is the state vector. Let h ( x ) : R n R be a C differentiable function; h S denotes the restriction of h x on a set S R n . The function h x used in this statement is called localizing. S ( h ) denotes the set x R n L f h x = 0 , where L f h ( x ) represents the Lie derivative of f x and is given by L f h x = h / x f x . From the latter, one can define the lower and upper bounds h inf = inf h x x S h and h sup = sup h x x S h . Hence, the general theorem concerning the localization of all compact invariant sets of a dynamical system establishes the following:
Theorem A1.
General theorem. Section 2 by Krishchenko and Starkov in [25]. Each compact invariant set Γ of x ˙ = f x is contained in the localizing domain:
K ( h ) = h inf h x h sup .
If the location of all compact invariant sets is inside the domain Λ R n , then the set K h Λ is formulated. Further, if all compact invariant sets are located in the sets K h i and K h j , with K h i , K h j R n , then they are also located in the set K h i K h j . Hence, a refinement of Theorem A1 is realized with the help of the iterative theorem as follows:
Theorem A2.
Iterative theorem. See Section 2 in [25]. Let h m x , m = 0 , 1 , 2 , be a sequence of C differentiable functions. Sets
K 0 = K h 0 , K m = K m 1 K m 1 , m , m > 0 ,
with
K m 1 , m = x : h m , inf h m x h m , sup , h m , sup = sup S ( h m ) K m 1 h m x , h m , inf = inf S ( h m ) K m 1 h m x ,
contain any compact invariant set of the system x ˙ = f x and
K 0 K 1 K m ·

References

  1. Saeedi, P.; Salpea, P.; Karuranga, S.; Petersohn, I.; Malanda, B.; Gregg, E.W.; Unwin, N.; Wild, S.H.; Williams, R. Mortality Attributable to Diabetes in 20–79 Years Old Adults, 2019 Estimates: Results from the International Diabetes Federation Diabetes Atlas. Diabetes Res. Clin. Pract. 2020, 162, 108086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Williams, R.; Karuranga, S.; Malanda, B.; Saeedi, P.; Basit, A.; Besançon, S.; Bommer, C.; Esteghamati, A.; Ogurtsova, K.; Zhang, P.; et al. Global and regional estimates and projections of diabetes-related health expenditure: Results from the international diabetes federation diabetes atlas. Diabetes Res. Clin. Pract. 2020, 162, 108072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Makroglou, A.; Li, J.; Kuang, Y. Mathematical models and software tools for the glucose-insulin regulatory system and diabetes: An overview. Appl. Numer. Math. 2006, 56, 559–573. [Google Scholar] [CrossRef]
  4. Gaetano, A.D.; Hardy, T.; Beck, B.; Abu-Raddad, E.; Palumbo, P.; Bue-Valleskey, J.; Pørksen, N. Mathematical models of diabetes progression. Am. J. Physiol.-Endocrinol. Metab. 2008, 295, E1462–E1479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Nath, A.; Biradar, S.; Balan, A.; Dey, R.; Padhi, R. Physiological models and control for type 1 diabetes mellitus: A brief review. IFAC-PapersOnLine 2018, 51, 289–294. [Google Scholar] [CrossRef]
  6. Kovatchev, B. A century of diabetes technology: Signals, models, and artificial pancreas control. Trends Endocrinol. Metab. 2019, 30, 432–444. [Google Scholar] [CrossRef] [PubMed]
  7. Tejedor, M.; Woldaregay, A.Z.; Godtliebsen, F. Reinforcement learning application in diabetes blood glucose control: A systematic review. Artif. Intell. Med. 2020, 104, 101836. [Google Scholar] [CrossRef] [PubMed]
  8. Vera, J.; Lischer, C.; Nenov, M.; Nikolov, S.; Lai, X.; Eberhardt, M. Mathematical modelling in biomedicine: A primer for the curious and the skeptic. Int. J. Mol. Sci. 2021, 22, 547. [Google Scholar] [CrossRef] [PubMed]
  9. Rathee, S. ODE models for the management of diabetes: A review. Int. J. Diabetes Dev. Ctries 2016, 37, 4–15. [Google Scholar] [CrossRef]
  10. Palmer, A.J. Computer modeling of diabetes and its complications: A report on the fifth mount hood challenge meeting. Value Health 2013, 16, 670–685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Valle, P.A.; Coria, L.N.; Salazar, Y. Tumor clearance analysis on a cancer chemo-immunotherapy mathematical model. Bull. Math. 2019, 81, 4144–4173. [Google Scholar] [CrossRef]
  12. Starkov, K.E. A cancer model for the angiogenic switch and immunotherapy: Tumor eradication in analysis of ultimate dynamics. Int. J. Bifurcation Chaos 2020, 30, 2050150. [Google Scholar] [CrossRef]
  13. Gamboa, D.; Coria, L.N.; Cárdenas, J.R.; Ramírez, R.; Valle, P.A. Hardware implementation of a non-linear observer for a diabetes mellitus type 1 mathematical model. Comput. Sist. 2019, 23, 1475–1486. [Google Scholar] [CrossRef]
  14. Gamboa, D.; Vázquez, C.E.; Campos, P.J. Nonlinear analysis for a type-1 diabetes model with focus on t-cells and pancreatic β-cells behavior. Math. Comput. Appl. 2020, 25, 23. [Google Scholar] [CrossRef]
  15. Lombarte, M.; Lupo, M.; Brenda, L.F.; Campetelli, G.; Marilia, A.B.; Basualdo, M.; Rigalli, A. In vivo measurement of the rate constant of liver handling of glucose and glucose uptake by insulin-dependent tissues, using a mathematical model for glucose homeostasis in diabetic rats. J. Theor. Biol. 2018, 439, 205–215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Lombarte, M.; Lupo, M.; Campetelli, G.; Basualdo, M.; Rigalli, A. Mathematical model of glucose–insulin homeostasis in healthy rats. Math. Biosci. 2013, 245, 269–277. [Google Scholar] [CrossRef]
  17. Brenner, M.; Kwon, G.; Lee, H.F.; Johns, M.; Malik, N. Simulating and predicting glucose response in lean and obese mice. In Proceedings of the 2014 40th Annual Northeast Bioengineering Conference (NEBEC), Boston, MA, USA, 25–27 April 2014. [Google Scholar]
  18. Lombarte, M.; Fina, B.; Lupion, P.; Lupo, M.; Rigalli, A. In vivo measurementof fluoride effects on glucose homeostasis: An explanation for the decrease in intelligence quotient andinsulin resistance induced by fluoride. Fluoride 2016, 3 Pt 1, 204–210. [Google Scholar]
  19. Khalil, H. Nonlinear Systems, 3rd ed.; Prentice-Hall: Hoboken, NJ, USA, 2002. [Google Scholar]
  20. Ha, J.; Sherman, A. Type 2 diabetes: One disease, many pathways. Am. J. Physiol.-Endocrinol. Metab. 2020, 319, E410–E426. [Google Scholar] [CrossRef] [PubMed]
  21. Sánchez, O.D.; Ruiz-Velázquez, E.; Alanís, A.Y.; Quiroz, G.; Torres-Treviño, L. Parameter estimation of a meal glucose–insulin model for TIDM patients from therapy historical data. IET Syst. Biol. 2019, 13, 8–15. [Google Scholar] [CrossRef]
  22. Wallace, T.M.; Levy, J.C.; Matthews, D.R. Use and abuse of HOMA modeling. Diabetes Care 2004, 27, 1487–1495. [Google Scholar] [CrossRef] [Green Version]
  23. Man, C.D.; Micheletto, F.; Lv, D.; Breton, M.; Kovatchev, B.; Cobelli, C. The UVA/PADOVA type 1 diabetes simulator. J. Diabetes Sci. Technol. 2014, 8, 26–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Krishchenko, A.P. Localization of invariant compact sets of dynamical systems. Differ. Equ. 2005, 41, 1669–1676. [Google Scholar] [CrossRef]
  25. Krishchenko, A.P.; Starkov, K.E. Localization of compact invariant sets of the lorenz system. Phys. Lett. A 2006, 353, 383–388. [Google Scholar] [CrossRef]
Figure 1. Graphical explanation of the contribution of each differential equation and its relationships with the parameters.
Figure 1. Graphical explanation of the contribution of each differential equation and its relationships with the parameters.
Axioms 11 00320 g001
Figure 2. Numerical simulation for type-1 diabetes mellitus system (1)–(5): Solution for G ( t ) (plasma glucose level), taking initial conditions I ( t ) = 10 , Y ( t ) = 10 , D ( t ) = 800 , G ( t ) = 800 , and U ( t ) = 550 , and parameter values of Table 1. G * ( t ) is the equilibrium point, and G m a x is the upper bound.
Figure 2. Numerical simulation for type-1 diabetes mellitus system (1)–(5): Solution for G ( t ) (plasma glucose level), taking initial conditions I ( t ) = 10 , Y ( t ) = 10 , D ( t ) = 800 , G ( t ) = 800 , and U ( t ) = 550 , and parameter values of Table 1. G * ( t ) is the equilibrium point, and G m a x is the upper bound.
Axioms 11 00320 g002
Figure 3. Numerical simulation for type-1 diabetes mellitus system (1)–(5): Solution for I ( t ) (plasma insulin variation), taking initial conditions I ( t ) = 10 , Y ( t ) = 10 , D ( t ) = 800 , G ( t ) = 800 , and U ( t ) = 550 , and parameter values of Table 1. I * ( t ) is the equilibrium point, and I m a x is the upper bound.
Figure 3. Numerical simulation for type-1 diabetes mellitus system (1)–(5): Solution for I ( t ) (plasma insulin variation), taking initial conditions I ( t ) = 10 , Y ( t ) = 10 , D ( t ) = 800 , G ( t ) = 800 , and U ( t ) = 550 , and parameter values of Table 1. I * ( t ) is the equilibrium point, and I m a x is the upper bound.
Axioms 11 00320 g003
Figure 4. Numerical simulation for type-2 diabetes mellitus system (1)–(5): Solution for G ( t ) (plasma glucose level), taking initial conditions I ( t ) = 1800 , D ( t ) = 3000 , G ( t ) = 300 , and U ( t ) = 100 , and parameter values of Table 1. G * ( t ) is the equilibrium point, and G m a x is the upper bound.
Figure 4. Numerical simulation for type-2 diabetes mellitus system (1)–(5): Solution for G ( t ) (plasma glucose level), taking initial conditions I ( t ) = 1800 , D ( t ) = 3000 , G ( t ) = 300 , and U ( t ) = 100 , and parameter values of Table 1. G * ( t ) is the equilibrium point, and G m a x is the upper bound.
Axioms 11 00320 g004
Figure 5. Numerical simulation for type-2 diabetes mellitus system (1)–(5): Solution for I ( t ) (plasma insulin variation), taking initial conditions I ( t ) = 1800 , D ( t ) = 3000 , G ( t ) = 300 , and U ( t ) = 100 , and parameter values of Table 1. I * ( t ) is the equilibrium point, and I m a x is the upper bound.
Figure 5. Numerical simulation for type-2 diabetes mellitus system (1)–(5): Solution for I ( t ) (plasma insulin variation), taking initial conditions I ( t ) = 1800 , D ( t ) = 3000 , G ( t ) = 300 , and U ( t ) = 100 , and parameter values of Table 1. I * ( t ) is the equilibrium point, and I m a x is the upper bound.
Axioms 11 00320 g005
Table 1. Parameters reported in [15].
Table 1. Parameters reported in [15].
ParameterT1DMT2DM
k a 0.206 0.060
k 0 0.160 0.026
k 1 0.002 0.330
k 2 0.0067 0.000235
k 3 2.244 2.342
k 4 0.017 0.008
k 5 0.009 0.258
k 6 0.170 0.052
k 7 2.397 -
k 8 0.934 -
G u 352.4 232.0
I p i 156.3 1103.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gamboa, D.; Coria, L.N.; Valle, P.A. Ultimate Bounds for a Diabetes Mathematical Model Considering Glucose Homeostasis. Axioms 2022, 11, 320. https://doi.org/10.3390/axioms11070320

AMA Style

Gamboa D, Coria LN, Valle PA. Ultimate Bounds for a Diabetes Mathematical Model Considering Glucose Homeostasis. Axioms. 2022; 11(7):320. https://doi.org/10.3390/axioms11070320

Chicago/Turabian Style

Gamboa, Diana, Luis N. Coria, and Paul A. Valle. 2022. "Ultimate Bounds for a Diabetes Mathematical Model Considering Glucose Homeostasis" Axioms 11, no. 7: 320. https://doi.org/10.3390/axioms11070320

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop