hello Chiara
I tried to fix this issue and got to this solution (attached)
To me and according to the documentation of the 2D PID block, your input signals were flipped
I tried manually to tune the P and I terms (no D required IMHO) and the best results I could get was with
Kp = 0.1
Ki = 0.0003
the very small Ki gain is due to the huge time delay in your model

I don"t have the toolbox to do the automatic tuning, so I let this to you
I removed a few things to speed up the simulation (the Step Response Checker seems to slow it quite significantly)
also I introduced a -1 gain block after the PID just to compensate for the negative dc gain in your tf (my mental preferences is to have positive dc gain tf)
attached wy revised simulink model
maybe there are better solution for first order tf with delays than just a regular PID (lots of publications available on internet)
all the best
