A Nyström method for linear second kind Volterra integral equations on unbounded intervals, with sufficiently smooth kernels, is described. The procedure is based on the use of a truncated Lagrange interpolation process and of a truncated Gaussian quadrature formula. The stability and the convergence of the method in suitable weighted spaces of functions are studied and some numerical examples showing its reliability are presented. In particular, the proposed method has been tested for the numerical resolution of some Volterra integral equations arising from the reformulation of differential models describing metastatic tumor growth whose unknown solutions represent biological observables as the metastatic mass or the number of metastases.