Summary: HIV is one of the major causes of deaths, especially in Sub-Saharan Africa. In this paper, an in vivo deterministic model of differential equations is presented and analyzed for HIV dynamics. Optimal control theory is applied to investigate the key roles played by the various HIV treatment strategies. In particular, we establish the optimal strategies for controlling the infection using three treatment regimes as the system control variables. We have applied Pontryagin’s maximum principle in characterizing the optimality control, which then has been solved numerically by applying the Runge-Kutta forth-order scheme. The numerical results indicate that an optimal controlled treatment strategy would ensure significant reduction in viral load and also in HIV transmission. It is also evident from the results that protease inhibitor plays a key role in virus suppression; this is not to underscore the benefits accrued when all the three drug regimes are used in combination.