  A description of the NN scattering is constructed starting from the $\pi,
\eta, \eta'$ pseudoscalar-, the $\rho, \phi, \omega$ vector-, and the
$\varepsilon(600), a_0, f_0(1400)$ scalar - meson-nucleon coupling constants,
which we obtain within a relativistic quark model. Working within the
Blankenbecler-Sugar-Logunov-Tavkhelidze quasipotential dynamics we thus
describe the NN phase shifts in a relativistically invariant way. In this
procedure we use the phenomenological form factor cutoff masses and the
effective $\varepsilon$ and $\omega$ meson-nucleon coupling constants, only.
The comparison of our NN phase shifts to the both empirical data and the Bonn
OBEP fit shows good agreement -- the ratio of the $\chi^2$ for the present
results to the $\chi^2$ for Bonn OBEP description is 1.2.
