We propose a modification of collocation methods extending the `averaged vector field method' to high order. These new integrators exactly preserve energy for Hamiltonian systems, are of arbitrarily high order, and fall into the class of B-series integrators. We discuss their symmetry and conjugate-symplecticity, and we compare them to energy-preserving composition methods.