Statistics of A simple mass-action model for the eukaryotic heat shock response and its mathematical validation