The problem of calculating the mass spectrum of the two-body Bethe-Salpeter equation is studied with no reduction to the three-dimensional (''quasipotential'') equation. The method to find the ground state and excited states for a channel with any quantum numbers is presented. The problem of the confining interaction for the Bethe-Salpeter equation is discussed from the point of view of formal properties of the bound state spectrum, but with only inspiration from QCD. We study the kernel that is nonvanishing at large Euclidean intervals, i.e. R(E) --> infinity, which is constructed as a special limiting case of the sum of the covariant one-boson-exchange kernels. In the coordinate space this kernel is just a positive constant and corresponds to the kernel proportional to delta(k(E)) in the momentum space. When the usual attractive interaction is added, it is found that this kernel is similar in its effect to the nonrelativistic potential in coordinate space, V(r), with V(r --> infinity). The positive real constant V-infinity gives the scale that defines the limit of the bound state spectrum compared to the sum of the constituent masses, M < 2m + V-infinity. At the same time, the self-energy corrections remove the singularities from the propagators of the constituents, i.e. constituents do not propagate as free particles. The combination of these features of the solutions allows an interpretation of this type of interaction as a confining one. The illustrative analytical and numerical calculations are presented for a model of massive scalar particles with scalar interaction, i.e. the ''massive Wick model.''