In the paper, an universality theorem on the approximation of analytic functions by generalized discrete shifts of zeta functions of Hecke-eigen cusp forms is obtained. These shifts are defined by using the function having continuous derivative satisfying certain natural growth conditions and, on positive integers, uniformly distributed modulo 1.