This paper is dedicated to providing new tools and methods for studying the trend to equilibrium of gradient flows in metric spaces (m, d) in the entropy and metric sense, to establish decay rates and finite time of extinction, and to characterize Lyapunov stable equilibrium points. More precisely, our main results are as follows: Introduction of a gradient inequality in the metric space framework, which in the Euclidean space R-N was obtained by Lojasiewicz [Editions du Centre National de la Recherche Scientifique, Paris, 1963, pp. 87-89], later improved by Kurdyka [Ann. Inst. Fourier 48 (1998), no. 3, 769-783], and generalized to the Hilbert space framework by Simon [Ann. of Math. (2) 118 (1983), no. 3, 525-571]. Obtainment of the trend to equilibrium in the entropy and metric sense of gradient flows generated by a function epsilon : M -> (-infinity, +infinity] satisfying a Kurdyka Lojasiewicz Simon inequality in a neighborhood of an equilibrium point of epsilon. Sufficient conditions are given implying decay rates and finite time of extinction of gradient flows. Construction of a talweg curve in m with an optimal growth function yielding the validity of a Kurdyka-Lojasiewicz-Simon inequality. Characterization of Lyapunov stable equilibrium points of epsilon satisfying a Kurdyka-Lojasiewicz-Simon inequality near such points. Characterization of the entropy-entropy production inequality with the Kurdyka-Lojasiewicz-Simon inequality. As an application of these results, the following properties are established. New upper bounds on the extinction time of gradient flows associated with the total variational flow. If the metric space m is the p-Wasserstein space P-p (R-N) < p < infinity, then new HWI-, Talagrand, and logarithmic Sobolev inequalities are obtained for functions epsilon associated with nonlinear diffusion problems modeling drift, potential, and interaction phenomena. It is shown that these inequalities are equivalent to the Kurdyka-Lojasiewicz-Simon inequality, and hence they imply a trend to equilibrium of the gradient flows of epsilon with decay rates or arrival in finite time.