A new Skyrme-like energy density suitable for studies of strongly elongated nuclei has been determined in the framework of the Hartree-Fock-Bogoliubov theory using the recently developed model-based, derivative-free optimization algorithm POUNDerS. A sensitivity analysis at the optimal solution has revealed the importance of states at large deformations in driving the parameterization of the functional. The good agreement with experimental data on masses and separation energies, achieved with the previous parameterization UNEDF0, is largely preserved. In addition, the new energy density UNEDF1 gives a much improved description of the fission barriers in 240Pu and neighboring nuclei.