First, in a continuous-time LINEAR TIME-INVARIANT (LTI) system, the Laplace Transform manages to give you the frequency response of your system. First, you compute the Laplace Transform H(s) of its impulse response h(t) and for reasons that are shown below, this can be called Transfer Function (TF). Then, it can be shown that if you compute the Laplace Transform of the time-response of sin(wt), you get the product of the TF H(jw)=H(s=jw) with the Laplace Transform of sin(wt). This is why H(s) can also be called Frequency Response, because a simple product gives you the response of the linear systems for any input frequency. Moreover, next computation managed to show that, in linear systems, for any input u(t), the (transform of the) convolution integral that would give you the output y(t) in time-domain, becomes a simple product of H(s) with U(s) in Laplace domain, that can now be called Frequency domain. That’s why one one can separate the linear system H(s) from the input U(s) and get the output Y(s) as a simple product Y(s)=H(s)U(s).
In discrete linear systems, you get similar relations, only you use the z variable.
All that ends when you move away from LTI systems. While we still can, and may also want, to compute the transform of the input command and of the output signal, there is no simple relation with some transform of the nonlinear system.
That’s why, for example, stability analysis in linear systems is simply based on the Transfer Function (poles and zeros), while in nonlinear systems becomes some pretty complex analysis (Lyapunov, etc.)
Sorry if this is disappointing, yet nonlinear systems are complex and there is no use of s-transform or z-transform of the nonlinear SYSTEMS.