Model reductions in biochemical reaction networks