This paper develops a mathematical model to investigate the Human Immunodeficiency Virus (HIV) infection dynamics. The model includes two transmission modes (cell-to-cell and cell-free), two adaptive immune responses (cytotoxic T-lymphocyte (CTL) and antibody), a saturated CTL immune response, and latent HIV infection. The existence and local stability of equilibria are fully characterized by four reproduction numbers. Through sensitivity analyses, we assess the partial rank correlation coefficients of these reproduction numbers and identify that the infection rate via cell-to-cell transmission, the number of new viruses produced by each infected cell during its life cycle, the clearance rate of free virions, and immune parameters have the greatest impact on the reproduction numbers. Additionally, we compare the effects of immune stimulation and cell-to-cell spread on the model’s dynamics. The findings highlight the significance of adaptive immune responses in increasing the population of uninfected cells and reducing the numbers of latent cells, infected cells, and viruses. Furthermore, cell-to-cell transmission is identified as a facilitator of HIV transmission. The analytical and numerical results presented in this study contribute to a better understanding of HIV dynamics and can potentially aid in improving HIV management strategies.