We introduce a systematic construction of higher-order matrix product operator (MPO) approximations of the time evolution operator for generic (short and long range) one-dimensional Hamiltonians. We demonstrate the utility of our construction, by showing an order of magnitude speedup in simulation cost compared to conventional first-order MPO time evolution schemes.