Towards energy-stable and conservative discontinuous Galerkin spectral element methods for Einstein's equations of general relativity in second order form