Numerical approximation of fractional order HIV/AIDS model with treatment using Levenberg-Marquardt backpropagation

Muhammad Awais1, Tayyab Nawaz2, Abu Safyan Ali1

  • 1Department of Mathematics and Computer Science, University of Ferrara, Via Machiavelli 30, Ferrara 44121, Italy.

PubMed

HIV specifically targets the immune system's T cells, mainly the CD4+ T cells, leading to a chronic and severe illness characterized by a long incubation period. In this study, we propose a novel nonlinear fractional-order model for analyzing and controlling HIV proliferation dynamics. The model categorizes the HIV-infected population into four distinct compartments: susceptible individuals (S), acutely infected (I1), asymptomatically infected (I2), and individuals with AIDS (A). Mathematical properties like positivity and boundedness provide for model validation. The analytical strategy blends practical wisdom with accurate mathematical analysis. We calculated two equilibrium states - disease-free equilibrium (DFE) and endemic equilibrium (EE)- and established their local stability regions. The analysis demonstrates that DFE is stable when the basic reproduction number R0<1, whereas EE becomes stable when R0>1. The Normalized Sensitivity Index (NSI) technique evaluates the model's sensitivity to key biological parameters. The fractional derivative is defined in the Caputo sense, and numerical solutions are obtained via the Adams-Bashforth-Moulton (ABM) method. To enhance computational efficiency and accuracy, an artificial neural network (ANN) with 10 hidden layers is trained using Levenberg-Marquardt backpropagation (LMBNN). Three distinct cases of the fractional derivative order (θ=0.55,0.7,0.9) are investigated. The LMBNN framework minimizes the Mean Squared Error (MSE) across training (80%), validation (10%), and testing (10%) phases using the ABM-generated solutions as the reference dataset. The numerical results demonstrate that the LMBNN approach achieves high precision, with MSE values as low as 10-3, supported by regression analysis (R2≈1). Simulation results suggest that the proposed model, coupled with the LMBNN technique, effectively predicts the control of HIV progression to AIDS, thereby offering significant potential for optimizing treatment strategies. Reducing contact with infected individuals and accelerating disease management interventions can significantly lower the burden of infection. Ultimately, our study offers a bridge to improve the forecasting and controlling of HIV/AIDS disease challenges, effectively connecting surrogate mathematics formulation with actual disease treatment measures.

Related Concept Videos

Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving01:29

Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving

Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
277
Linear Approximation in Frequency Domain01:26

Linear Approximation in Frequency Domain

Linear systems are characterized by two main properties: superposition and homogeneity. Superposition allows the response to multiple inputs to be the sum of the responses to each individual input. Homogeneity ensures that scaling an input by a scalar results in the response being scaled by the same scalar.
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
341
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation01:24

One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation

This lesson introduces two critical methods in pharmacokinetics, the Wagner-Nelson and Loo-Riegelman methods, used for estimating the absorption rate constant (ka) for drugs administered via non-intravenous routes. The Wagner-Nelson method relates ka to the plasma concentration derived from the slope of a semilog percent unabsorbed time plot. However, it is limited to drugs with one-compartment kinetics and can be impacted by factors like gastrointestinal motility or enzymatic degradation.
On...
1.1K
Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model01:13

Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model

Drugs administered through various routes can lead to nonlinear elimination, resulting in complex pharmacokinetic behaviors crucial to understanding efficacious drug dosing.
When a drug is administered through a constant intravenous infusion and eliminated via nonlinear pharmacokinetics, it follows zero-order input. For example, oral drugs undergo first-order absorption upon administration and are eliminated through nonlinear pharmacokinetics.
In the case of subcutaneously administered drugs,...
281
Linear Approximation in Time Domain01:21

Linear Approximation in Time Domain

Nonlinear systems often require sophisticated approaches for accurate modeling and analysis, with state-space representation being particularly effective. This method is especially useful for systems where variables and parameters vary with time or operating conditions, such as in a simple pendulum or a translational mechanical system with nonlinear springs.
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
334
Model Approaches for Pharmacokinetic Data: Distributed Parameter Models01:06

Model Approaches for Pharmacokinetic Data: Distributed Parameter Models

Pharmacokinetic models are mathematical constructs that represent and predict the time course of drug concentrations in the body, providing meaningful pharmacokinetic parameters. These models are categorized into compartment, physiological, and distributed parameter models.
The distributed parameter models are specifically designed to account for variations and differences in some drug classes. This model is particularly useful for assessing regional concentrations of anticancer or...
231