[step:Transfer total boundedness from $M$ to $\hat{M}$ using density]Let $\iota: M \hookrightarrow \hat{M}$ denote the isometric embedding of $M$ into its completion $\hat{M}$. By the definition of completion, $\iota(M)$ is dense in $\hat{M}$ and $\hat{d}(\iota(x), \iota(y)) = d(x, y)$ for all $x, y \in M$.
Fix $\varepsilon > 0$. Since $M$ is totally bounded, there exists a finite set $\{x_1, \ldots, x_N\} \subset M$ with $M \subset \bigcup_{i=1}^N B_d(x_i, \varepsilon/2)$.
We claim $\{\iota(x_1), \ldots, \iota(x_N)\}$ is an $\varepsilon$-net for $\hat{M}$. Let $\hat{x} \in \hat{M}$. Since $\iota(M)$ is dense in $\hat{M}$, there exists $y \in M$ with $\hat{d}(\hat{x}, \iota(y)) < \varepsilon/2$. Since $y \in M \subset \bigcup_{i=1}^N B_d(x_i, \varepsilon/2)$, there exists $i$ with $d(y, x_i) < \varepsilon/2$. By the isometry property, $\hat{d}(\iota(y), \iota(x_i)) = d(y, x_i) < \varepsilon/2$. The triangle inequality gives
\begin{align*}
\hat{d}(\hat{x}, \iota(x_i)) \le \hat{d}(\hat{x}, \iota(y)) + \hat{d}(\iota(y), \iota(x_i)) < \frac{\varepsilon}{2} + \frac{\varepsilon}{2} = \varepsilon.
\end{align*}
Therefore $\hat{M} \subset \bigcup_{i=1}^N B_{\hat{d}}(\iota(x_i), \varepsilon)$, and $\hat{M}$ is totally bounded.[/step]