We present a new approach to the definition of Grand and Small discrete Lebesgue spaces. Previously, we developed such an approach in the continuous case. The fundamental difference of our studies is that we base the definition of norms on extrapolation theory, and our approach permits including the extreme cases p = 1 and p = infinity into consideration, which is the main result of the study presented here. The Small space for p = infinity is realized as the sum l(s)(xi) = Sigma(p is an element of[2,infinity))xi(p)l(p) of spaces with a certain fundamental function., and the Grand space l(g)(zeta) for p = 1 is given by the product l(g)(zeta) = boolean AND (p is an element of( 1,2])zeta(p)l(p) of spaces with a certain fundamental function zeta. As one of the main results, we show that if the function xi satisfies the Delta(2) (doubling) condition, then the space l(s)(xi) coincides, up to norm equivalence, with the discrete Lorentz space lambda(psi), where psi(k) similar or equal to xi(1/ ln k). We also show that if the function. satisfies the Delta(2) condition, then the space l(g)(zeta) coincides, up to norm equivalence, with the discrete Marcinkiewicz space m(psi), where psi(k) similar or equal to k/(zeta(1/ ln k)). We expect that our general new construction of the norms in Grand and Small discrete Lebesgue spaces will imply further studies of the spaces and operators in these spaces in such a general setting.