We prove that nonsmooth quasilinear parabolic systems admit a local solution in Lp strongly differentiable with respect to time over a bounded three-dimensional polyhedral space domain. The proof rests essentially on new elliptic regularity results for polyhedral Laplace interface problems for anisotropic materials. These results are based on sharp pointwise estimates for Green's function, which are also of independent interest. To treat the nonlinear problem, we then apply a classical theorem of Sobolevskii for abstract parabolic equations and recently obtained resolvent estimates for elliptic operators and interpolation results. As applications we have in mind primarily reaction-diffusion systems. The treatment of such equations in an Lp context seems to be new and allows (by Gauss' theorem) the proper definition of the normal component of currents across the boundary.