The present study aims to clarify the role of the fraction of patients under antiretroviral therapy (ART) achieving viral suppression (VS) (i.e. having plasma viral load below the detectability threshold) on the human immunodeficiency virus (HIV) epidemic in Italy. Based on the hypothesis that VS makes the virus untransmittable, we extend a previous model and we develop a time-varying ordinary differential equation model with immigration and treatment, where the naive and non-naive populations of infected are distinguished, and different compartments account for treated subjects virally suppressed and not suppressed. Moreover, naive and non-naive individuals with acquired immune deficiency syndrome (AIDS) are considered separately. Clinical data stored in the nationwide database Antiviral Response Cohort Analysis are used to reconstruct the history of the fraction of virally suppressed patients since highly active ART introduction, as well as to assess some model parameters. Other parameters are set according to the literature and the final model calibration is obtained by fitting epidemic data over the years 2003-2015. Predictions on the evolution of the HIV epidemic up to the end of 2035 are made assuming different future trends of the fraction of virally suppressed patients and different eligibility criteria for treatment. Increasing the VS fraction is found to reduce the incidence, the new cases of AIDS and the deaths from AIDS per year, especially in combination with early ART initiation. The asymptotic properties of a time-invariant formulation of the model are studied, and the existence and global asymptotic stability of a unique positive equilibrium are proved.
Keywords: HIV epidemic; antiretroviral therapy; epidemic ODE model; stability of equilibrium points; viral suppression.
© The Author(s) 2019. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.