The cubature on Wiener space method, a high-order weak approximation scheme, is established for SPDEs in the case of unbounded characteristics and unbounded test functions. We first describe a recently introduced flexible functional analytic framework, so called weighted spaces, where Feller-like properties hold. A refined analysis of vector fields on weighted spaces then yields optimal convergence rates of cubature methods for stochastic partial differential equations of Da PratoZabczyk type. The ubiquitous stability for the local approximation operator within the functional analytic setting is proved for SPDEs, however, in the infinite dimensional case we need a newly introduced technical assumption on weak symmetry of the cubature formula. Computational results for a cubature discretization of a spatially extended stochastic FitzHugh-Nagumo model, an SPDE model from mathematical biology, are shown, illustrating the applicability of our theory.