Abstract
Sphingolipids have been identified as important signaling compounds in stress responses. However, it is not always clear how different sphingolipid profiles are achieved in a particular stress situation. Here we propose a detailed mass action model, containing 42 dependent variables and 137 reactions, that offers explanations of the roles of variant ceramides species, which differ in the lengths of their fatty acyl chains and their saturation state, in the response to hydroxyurea stress. The simulations demonstrate that the cells manage to achieve hydroxyurea tolerance through a well-coordinated, differential usage of the variant ceramide species. Moreover, the results suggest that key enzymes have different affinities toward saturated and unsaturated fatty acyl chains, which implies that the saturation state affords the cells with an additional mode of regulation that had not been recognized so far. These conclusions from our computational analysis are yet to be validated experimentally.