Scalable parameter estimation for genome-scale biochemical reaction networks