As remarked above by Yuri Vyatkin, the Hodge * being an isometry is more or less the same as the double Hodge being $(-1)^{k(n-k)}$ times the identity.
This is my attempt to prove it is an isometry (well,"isometry like": the sign of the innerproduct/nondegenerate symmetric bilinear form is allowed to flip sign) without using a basis. Only for simple (fully decomposable) $k$-vectors however. By (multi)linearity it should work on all of $\Lambda^k(V)$. I guess this is how many properties on $\Lambda(V)$ are defined.
Consider the Grassmann complement. Given a special $n$-covector $\omega$ and also a privileged $n$-vector $\Omega$, the Grassmann complement is a map $*:\Lambda^k(V)\to\Lambda^{n-k}(V^*)$ defined by $$i_{a_1\wedge\dots\wedge a_k}\omega(u_1,\dots,u_{n-k})=\omega(a_1,\dots,a_k,u_1,\dots,u_{n-k}). $$ Similarly $*:\Lambda^k(V^*)\to\Lambda^{n-k}(V)$ - with $\Omega$ instead of $\omega$. (Both directions being called * as well as the Hodge dual denoted * is confusing. Sorry.)
The action of a simple $k$-covector on a simple $k$-vector is given by the determinant
$$ \langle a_1\wedge\dots\wedge a_k,\alpha_1\wedge\dots\wedge \alpha_k\rangle =\left|\begin{array}{ccc}
\langle a_{1},\alpha_{1}\rangle & \cdots & \langle a_{1},\alpha_{k}\rangle\\
\vdots & \ddots\\
\langle a_{k},\alpha_{1}\rangle & & \langle a_{k},\alpha_{k}\rangle
\end{array}\right|.$$
To find $\langle *(\alpha_{1}\wedge\dots\wedge\alpha_{k}), *(a_{1}\wedge \dots\wedge a_k)\rangle$, we can use that the adjoint of contraction with $\alpha\in V^*$ is wedging with $\alpha$.
$$
\langle *(\alpha_{1}\wedge\dots\wedge\alpha_{k}), *(a_{1}\wedge \dots\wedge a_k)\rangle=\langle i_{\alpha_{1}\wedge\dots\wedge\alpha_{k}}\Omega, i_{a_{1}\wedge \dots\wedge a_k} \omega\rangle\\
=\langle \Omega,\alpha_{1}\wedge\dots\wedge\alpha_{k}\wedge i_{a_{1}\wedge \dots\wedge a_k} \omega\rangle
$$
We also need $i_a (\sigma \wedge \tau)=(i_a \sigma)\wedge \tau + (-1)^{\deg \sigma}\sigma\wedge i_a \tau $ ("$i_a$ is an anti-derivation").
For $k=1$, we simply get
$$
\alpha\wedge i_a\omega = -i_a(\alpha\wedge\omega)+\langle a,\alpha\rangle \omega=\langle a,\alpha\rangle \omega.
$$
Next, $k>1$. By bringing the $i_{a_k}$ to the front in the LHS below, using the anti-derivational character of $i_a$ you get an expansion of the determinant above. (When working with orthonormal/dual bases this step is particularly easy). We get
$$
\alpha_{1}\wedge\dots\wedge\alpha_{k}\wedge i_{a_{k}}\dots i_{a_{1}}\omega=\langle a_1\wedge\dots\wedge a_k,\alpha_1\wedge\dots\wedge \alpha_k\rangle\omega
$$
So
$$
\langle a_1\wedge\dots\wedge a_k,\alpha_1\wedge\dots\wedge \alpha_k\rangle\langle \Omega,\omega\rangle =\langle *(\alpha_{1}\wedge\dots\wedge\alpha_{k}), *(a_{1}\wedge \dots\wedge a_k)\rangle.
$$
In the case of the Hodge * dual we have the volume form $\omega=\sqrt{|g|}dx^1\wedge\dots\wedge dx^n$ and also a special $\Omega=\varphi(\omega)$, where $\varphi$ is the Riesz isomorphism (is it called that on Minkowski space too?). We should have $\langle\Omega,\omega\rangle=\sqrt{|g|}\cdot\sqrt{|g|}/\det{g}=\text{sign}(\text{metric})$
(The Clifford answer by Nicholas Todoroff looks pretty neat to me btw. I don't fully understand it.) First post, can't comment.