Gradient-based dimension reduction for solving Bayesian inverse problems