Network models are a well established tool for studying the robustness of complex systems, including modelling the effect of loss of function mutations in protein interaction networks. Past work has concentrated on average damage caused by random node removal, with little attention to the shape of the damage distribution. In this work, we use fission yeast co-expression networks before and after exposure to stress to model the effect of stress on mutational robustness. We find that exposure to stress decreases the average damage from node removal, suggesting stress induces greater tolerance to loss of function mutations. The shape of the damage distribution is also changed upon stress, with a greater incidence of extreme damage after exposure to stress. We demonstrate that the change in shape of the damage distribution can have considerable functional consequences, highlighting the need to consider the damage distribution in addition to average behaviour. 