Gradient-based optimization for variational empirical Bayes multiple regression