We establish well-posedness and maximal regularity estimates for linear parabolic SPDE in divergence form involving random coefficients that are merely bounded and measurable in the time, space, and probability variables. To reach this level of generality, and avoid any of the smoothness assumptions used in the literature, we introduce a notion of pathwise weak solution and develop a new harmonic analysis toolkit. The latter includes techniques to prove the boundedness of various maximal regularity operators on relevant spaces of square functions, the parabolic tent spaces Tp. Applied to deterministic parabolic PDE in divergence form with real coefficients, our results also give the first extension of Lions maximal regularity theorem on L-2(R+ x R-n) = T-2 to T-p, for all 1 - epsilon < p <= infinity in this generality. (c) 2024 The Author(s). Published by Elsevier Inc. This is an open access article under the CC BY license (http://creativecommons .org/licenses/by /4 .0/).