In this article we consider the KPZ fixed point starting from a two-sided Brownian motion with an arbitrary diffusion coefficient. We apply the integration by parts formula from Malliavin calculus to establish a key relation between the two-point (covariance) function of the spatial derivative process and the location of the maximum of an Airy process plus Brownian motion minus a parabola. Integration by parts also allows us to deduce the density of this location in terms of the second derivative of the variance of the KPZ fixed point. In the stationary regime, we find the same density related to limit fluctuations of a second-class particle. We further develop an adaptation of Stein's method that implies asymptotic independence of the spatial derivative process from the initial data.