In this work, we extend a mathematical model, which has been proposed for susceptible and infected compartments together with pathogen population, by including recovered subgroup. It is known that environmental pollution, such as contaminated drinking water and lack of an ordinary toilet, affects individuals and such negative impacts can be defined as "stressors." In order to include the influence of such stressors, susceptible subpopulation has been divided into two groups as the one affected by stress or not. Thus, spread of the disease is expressed in terms of a five-dimensional system. Moreover, we extend this model with the use of a time fractional derivative due to non-local effects of water pollution and we prove that the solution is non-negative and bounded from above. Then, we perform stability analysis for the disease-free equilibrium point. Afterward, the next step is to apply optimal control theory to optimize the decay rate of pathogens and the stress related parameters so that the number of infected individuals and the pathogen population can be minimized. Finally, we present some numerical results to find out the most appropriate control policy and the effect of the fractional order.