Exploring Therapeutic Dynamics: Mathematical Modeling and Analysis of Type 2 Diabetes Incorporating Metformin Dynamics

Type 2 diabetes (T2D) is a chronic metabolic disorder requiring effective management to avoid complications. Metformin is a first-line drug agent and is routinely prescribed for the control of glycemia, but its underlying dynamics are complicated and not fully quantified. This paper formulates a control-oriented and interpretable mathematical model that integrates metformin dynamics into a classic beta-cell–insulin–glucose (BIG) regulation system. The paper’s applicability to theoretical and clinical settings is enhanced by rigorous mathematical analysis, which guarantees the model is globally bounded, well-posed, and biologically meaningful. One of the key features of the study is its global stability analysis using Lyapunov functions, which demonstrates the asymptotic stability of critical equilibrium points under realistic physiological constraints. These findings support the predictive reliability of the model in explaining long-term glycemic regulation. Bifurcation analysis also clarifies the dynamic interplay between glucose production and utilization by identifying parameter thresholds that signify transitions between homeostasis and pathological states. Residual analysis, which detects Gaussian-distributed errors, underlines the robustness of the fitting process and suggests possible refinements by including temporal effects. Sensitivity analysis highlights the predominant effect of the initial dose of metformin on long-term glucose regulation and provides practical guidance for optimizing individual treatment. Furthermore, changing the two considered metformin parameters from their optimal values—altering the dose by ±50% and the decay rate by ±20%—demonstrates the flexibility of the model in simulating glycemic responses, confirming its adaptability and its potential for optimizing personalized treatment strategies.