Consider the parabolic equation $p: \mathbb R^2\to\mathbb R$

$$\partial_t p + b(t)\partial_x p + D(t,x)\partial^2_{xx}p=0,$$

where $b, D$ are functions nice enough. I look for the continuity of the derivatives $\partial p, \partial_x p$ of the solution. It is known by Nash's paper that, under very reasonable conditions on $b, D$, we have the Holder type continuity of $p$. Is there any work concerning such continuity analysis of $\partial_t p, \partial_x p$? 

PS : My idea is to consider $q:=\partial_x p$. Then   

$$\partial_t q + b(t)\partial_xq + D(t,x)\partial^2_{xx}q=-\partial_x D(t,x)\partial_{xx}p$$

is a similar parabolic equation for $q$ with an additional source. While the term $\partial_x D(t,x)\partial_{xx}p$ containt $\partial_{xx}p$, which makes the estimation even harder....