Matrix product simulations of non-equilibrium steady states of quantum spin chains