You have shown that $e^{-Ht}\colon L^2\to L^1$ is bounded with the given operator norm.  It follows that the dual operator
$$\left(e^{-Ht}\right)^*\colon (L^1)^* \to (L^2)^*$$
is bounded with the same operator norm.  Since $(L^1)^* = L^\infty$ and $(L^2)^* = L^2$ (these are isomorphisms, actually), it remains to check that $\left(e^{-Ht}\right)^* = e^{-Ht}$.

For this note that $$\int f(x)(e^{-Ht}g)(x)dx = \iint f(x)e(x,y;t)g(y)dxdy$$ where $e(x,y;t)$ is the heat kernel, which is real and symmetric in the variables $x,y$.  It follows that $\left(e^{-Ht}\right)^*$ is again integration against the heat kernel as claimed.