Constraining the Skyrme energy density functional with quantum Monte Carlo calculations