Identifying biochemical reaction networks from heterogeneous datasets