Kinetic discretization of the multidimensional PIDE model of gene regulatory networks