I reformulate you question as follows: suppose we a given a Riemannian metric. How can one decide whether there exists a complex structure $J$ such that the metric is Kähler w.r.t. this complex structure?
Below I will explain how to answer this question in practice, i.e., if say a Kähler metric is explicitly given. Actually, if the metric is given by explicit formulas, sometimes one can relatively easy give a negative answer, and in small dimensions as the rule one can also find the complex structure $J$
or prove its existence.
Easy version of the answer: Kähler metrics have the property that the curvatures commute with the complex structure: for any two vectors $X,Y$ we have
$$
[R(X,Y), J]=0. \ \ \ (\ast) $$
Let us view this as as a system of linear equations on $J$ (the vectors $X$ and $Y$ are the one we choose freely, and we know the curvature since we know the metric).
It is a very overdetermined system, in the case it has no nontrivial (i.e., $\ne Id$) solution (which is a generic case)
the metric is definitely not Kähler. In the case the space of solution is two-dimensional, one need to understant whether the solution has only two eigevalues and both eigenvalues are not real valued. This is again an algorithmically doable procedure. If not, then the metric is not Kähler. If yes, we have the only candidate for the complex structure and one simply need to check whether Nijenhuis is zero. In the case we have many solutions, one could try to enlarge the system $(\ast)$ to get only one solution by adding to it the equations $[(\nabla_XR)(Y,Z), J]=0$ and so on.
Formally correct but in practice not applyable version of the aswer.
One need to consider the holonomy group of the metric, and look whether it preserves a complex structure.
Relation between formally correct and easy versions of the answer. By Ambrose-Singer, the holonomy group is generated by $R(X,Y)$ and also by
$
(\tau^*R)(X,Y)
$
over all parallel translations $\tau$ along all possible curves. Since $
(\tau^*R)(X,Y)
$ is almost impossible to calculate in practice, in the easy version of the answer it is replaced by its derviative at the initial point, which is $(\nabla_XR)(Y,Z)$.