A computational protocol to probe the role of solvation effects on the reduction potential of azurine mutants