This work investigates the co-dynamics of Hepatitis E and HIV. Initially, we formulate a co-infection dynamics model of Hepatitis E and HIV. Then, we analyze each model and discuss their mathematical results. After that, we investigate the full model and present their basic mathematical results. A bifurcation analysis for full model is investigated. Further, we formulate a mathematical model with five controls. Optimal control model is formulated and the necessary results of the optimal control characterization are presented. Moreover, numerical results with different control strategies are presented. It is shown that each strategy has its own importance but for the disease elimination the combination of all the five controls at the same time can best decrease the disease burden from the community