The control of opportunistic infections among HIV infected individuals should be one of the major public health concerns in reducing mortality rate of individuals living with HIV/AIDS. In this study a deterministic co-infection mathematical model is developed to provide a quantification of treatment at each contagious stage against Pneumocystis Pneumonia (PCP) among HIV infected individuals on ART. The goal is to minimize the co-infection burden by putting the curable PCP under control. The disease-free equilibria for the HIV/AIDS sub-model, PCP sub-model and the co-infection model are shown to be locally asymptotically stable when their associated disease threshold parameter is less than a unity. By use of suitable Lyapunov functions, the endemic equilibria corresponding to HIV/AIDS and PCP sub-models are globally asymptotically stable whenever the HIV/AIDS related basic reproduction number <I>R</I><sub>0<I>H</I></sub> and the PCP related reproduction number <I>R</I><sub>0<I>P</I></sub> are respectively greater than a unity. The sensitivity analysis results implicate that the effective contact rates are the main mechanisms fueling the proliferation of the two diseases and on the other hand treatment efforts play an important role in reducing the incidence. The model numerical results reveal that PCP carriers have a considerable contribution in the transmission dynamics of PCP. Furthermore, treatment of PCP at all contagious phases significantly reduces the burden with HIV/AIDS and PCP co-infection.