I will present some results on the shape estimates of the density $\rho(t,x; y)$ of the solution $u(t,x)$ to stochastic partial differential equation:$\frac{\partial }{\partial t} u(t,x)=\frac{1}{2} \Delta u(t,x)+u\diamond\dot W(t,x)$,where $\dot W$ is a general Gaussian noise and $\diamond$ denotes the Wick product. I mainly concern with the asymptotic behavior of $\rho(t,x; y)$ when $y\rightarrow \infty$ or when $t\to 0+$. Both upper and lower bounds are obtained and these two bounds match each other modulo some multiplicative constants. If the initial data is positive, then $\rho(t,x;y)$ is supported on the positive half line $y\in [0, \infty)$ and in this case $\rho(t,x; 0+)=0$ and I will give n an upper bound for $\rho(t,x; y)$ when $y\rightarrow 0+$.