A general solution of the exchange problem in the high-nuclearity spin clusters (HNSC) containing arbitrary number of exchange-coupled centers and topology is developed. All constituent magnetic centers are supposed to possess well-isolated orbitally non-degenerate ground states so that the isotropic Heisenberg-Dirac-Van Vleck (HDVV) term is the leading part of the exchange spin Hamiltonian. Along with the HDVV term, we consider higher-order isotropic exchange terms (biquadratic exchange), as well as the anisotropic terms (anisotropic and antisymmetric exchange interactions and local single-ion anisotropies). All these terms are expressed as irreducible tensor operators (ITO). This allows us to take full advantage of the spin symmetry of the system. At the same time, we have also benefitted by taking into account the point group symmetry of the cluster, which allows us to work with symmetrized spin functions. This results in an additional reduction of the matrices to diagonalize. The approach developed here is accompanied by an efficient computational procedure that allows us to calculate the bulk magnetic properties (magnetic susceptibility, magnetization, and magnetic specific heat) as well as the spectroscopic properties of HNSC. Special attention is paid to calculate the magnetic excitations observed by inelastic neutron scattering (INS), their intensities, and their Q and temperature dependencies. This spectroscopic technique provides direct access to the energies and wave functions of the different spin states of the cluster; thus, it can be applied to spin clusters in order to obtain deep and detailed information on the nature of the magnetic exchange phenomenon. The general expression for the INS cross-section of spin clusters interacting by all kinds of exchange interactions, including also the single-ion zero-field splitting term, is derived for the first time. A closed-form expression is also derived for the particular case in which only the isotropic exchange interactions are involved. Finally this approach has been used to model the magnetic properties as well as the INS spectra of the polyoxometalate anion [Ni₉(OH)₃(H₂O)₆(HPO ₄)₂(PW₉O₃₄)₃] ^16-, which contains a central magnetic cluster formed by nine exchange-coupled Ni(II) ions surrounded by diamagnetic phosphotungstate ligands (PW₉O₃₄)^9-.