After simplifying the expression, what does mean?
If you want to simulate the system in transfer functions, you can try the following.
s = tf('s');
% Plant
Gp = 410.39/(s^2 + 4.376)
step(Gp, 20)
% PID controller
kp = - 0.0026657569;
ki = 0.0053315139;
kd = 0.0025512318;
Tf = 0.5;
Gc = pid(kp, ki, kd, Tf)
% Closed-loop transfer function
Gcl = minreal(feedback(Gc*Gp, 1))
step(Gcl, 20)