Variational scheme to compute protein reaction pathways using atomistic force fields with explicit solvent