Model reduction of biochemical reaction networks