Let $(M,g)$ be a Riemannian manifold and let $f \in C^{\infty}(M)$. Let $X$ be a smooth vector field on $M$. In smooth local coordinates $(x^i)$ on $M$, we can write $g = g_{ij} dx^i \otimes dx^j$ as well as $X = X^i \partial_{x^i}$. Now I have computed that $\langle\operatorname{grad} f, X \rangle_g = X^i\frac{\partial f}{\partial x^i}$ as follows. If we let $g^{ij}$ be the inverse matrix of $g_{ij}$ We have $$\begin{eqnarray*} \langle\operatorname{grad} f, X \rangle_g &=& \left\langle g^{ij} \frac{\partial f}{\partial x^i}\partial_{x^j} , X^k\partial_{x^k} \right\rangle_g\\ &=& g^{ij} \frac{\partial f}{\partial x^i}X^k\left\langle \partial_{x^j},\partial_{x^k} \right\rangle_g\\ &=& g_{jk}g^{ij} \frac{\partial f}{\partial x^i}X^k \\ &=& g_{kj}g^{ji}\frac{\partial f}{\partial x^i}X^k \hspace{1cm} (\text{by using that $g_{ij}$ is a symmetric matrix})\\ &=& X^i\frac{\partial f}{\partial x^i}. \end{eqnarray*}$$
My question is: Is there any geometric reason as to why $\langle \operatorname{grad} f, X \rangle_g$ should be independent of the metric $g$?
Absolutely. The exterior differential of a function, $df$, is a one-form which is defined independently of metric: $df(X) = X(f)$. In coordinates, $df = \partial_if dx^i$.
Choice of a nondegenerate metric induces an isomorphism between vector fields and $1$-forms by pairing: $g(X)(V) = \langle X,V\rangle$. In coordinates, this pairing corresponds to multiplying by the inverse of the metric, or "raising indices." This isomorphism is also known as a musical isomorphism: if $X$ is a vector field, then the one-form $X^\sharp$ is defined by $X^\sharp(V) = \langle X,V\rangle$.
The gradient $\operatorname{grad}(f)$ is defined to be the metric dual of the one-form $df$, i.e., $\operatorname{grad}(f) = (df)^\sharp$, so if you pair $\operatorname{grad}(f)$ with an arbitrary vector field $X$, you get $df(X) = X(f)$, which does not depend at all on the metric. The gradient is constructed so that, in pairing, all mentions of the metric cancel.