Nvm, this was easier than I originally thought. I ended up expanding the characteristic PIDF equation illustrated in the picture attached in the original post.
For anyone that is interested here is the solution:
CPz = tf([P],[1],Ts);
CIz = tf([I*Ts I*Ts],[2 -2],Ts);
CDz = tf([2*D*N -2*D*N],[2+N*Ts N*Ts-2],Ts);
CPIDz = CPz+CIz+CDz;
I verified the discrete transfer function with the original PIDF block in Simulink. The model differences were in the factor of 10^-10.