Note that
$$X_t := - \int_0^t H_s \, dB_s - \frac{1}{2} \int_0^t H_s^2 \, ds$$
defines an Itô process. Therefore we may apply Itô's formula to $f(x) := e^x$:
$$\begin{align*} f(X_t)-f(X_0) &= \int_0^t f'(X_s) \, dX_s + \frac{1}{2} \int_0^t f''(X_s) \, d\langle X \rangle_s \\ \Leftrightarrow M_t - 1 &= - \int_0^t M_s H_s \, dB_s - \frac{1}{2} \int_0^t M_s H_s^2 \, ds + \frac{1}{2} \int_0^t M_s H_s^2 \, ds \\ &= -\int_0^t M_s H_s \, dB_s \tag{1} \end{align*}$$
Here we have used that $M_t = f(X_t) = f'(X_t)=f''(X_t)$ and that the quadratic variation process $\langle X \rangle_t$ is given by $$\langle X \rangle_t = \int_0^t H_s^2 \, ds.$$
Since $H$ is bounded, the right-hand side of $(1)$ is a martingale. Consequently, $(M_t)_{t \geq 0}$ is a martingale. On the other hand, if we set $g(x) := e^{2x}$, then $g(X_t)= M_t^2$ and a similar calculation yields
$$M_t^2 - 1 = -2 \int_0^t M_s^2 H_s \, dB_s + \int_0^t M_s^2 H_s^2 \, ds.$$
Therefore, it follows from the definition of quadratic variation that
$$\langle M \rangle_t = \int_0^t M_s^2 H_s^2 \, ds.$$
Moreover, $(1)$ shows
$$\langle M,B \rangle_t = -\int_0^t M_s H_s \, ds.$$