The Cholesky factorization of the moment matrix is applied to discrete orthogonal polynomials on the homogeneous lattice. In particular, semiclassical discrete orthogonal polynomials, which are built in terms of a discrete Pearson equation, are studied. The Laguerre-Freud structure semiinfinite matrix that models the shifts by ±1 in the independent variable of the set of orthogonal polynomials is introduced. In the semiclassical case it is proven that this Laguerre-Freud matrix is banded. From the well-known fact that moments of the semiclassical weights are logarithmic derivatives of generalized hypergeometric functions, it is shown how the contiguous relations for these hypergeometric functions translate as symmetries for the corresponding moment matrix. It is found that the 3D Nijhoff-Capel discrete Toda lattice describes the corresponding contiguous shifts for the squared norms of the orthogonal polynomials. The continuous 1D Toda equation for these semiclassical discrete orthogonal polynomials is discussed and the compatibility equations are derived. It is also shown that the Kadomtesev-Petviashvilii This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.