In this paper, a mathematical model is developed to simulate the activation of regulatory T lymphocytes dynamics. The model considers the adaptive immune response and consists of epithelial cells, infected cells, free virus particles, helper and cytotoxic T lymphocytes, B lymphocytes, and regulatory T lymphocytes. A mathematical analysis was carried out to discuss the conditions of existence and stability of equilibrium solutions in terms of the basic reproductive number. In addition, the definitions and properties necessary to preserve the positivity and stability of the model are shown. The precision of these mathematical models can be affected by numerous sources of uncertainty, partly due to the balance between the complexity of the model and its predictive capacity to depict the biological process accurately. Nevertheless, these models can provide remarkably perspectives on the dynamics of infection and assist in identification specific immunological traits that improve our comprehension of immune mechanisms. The theoretical results are validated by numerical simulations using data reported in the literature. The construction, analysis, and simulation of the developed models demonstrate that the increased induced regulatory T lymphocytes effectively suppress the inflammatory response in contrast to similar cells at lower contents, playing a key role in maintaining self-tolerance and immune homeostasis.