How did an overall factor of theta_dot5 in the question later become sin(theta_5) in the comments ?!? Taking the original posting as correct, and ignoring theta_dot5, this is
sin(th1)*cos(th4) - cos(th1)*cos(th2)*cos(th3)*sin(th4) ...
+ cos(th1)*sin(th2)*sin(th3)*sin(th4)
I did the following calculuation by hand since it's straightforward AND because going to symbolic variables, Matlab seems incapable of simplifying this expression at all. Highly disappointing, although I don't think many people contend that Matlab symbolics is very adept at its job.
By inspection the expression becomes
ans1 = sin(th1)*cos(th4) - cos(th1)*cos(th2+th3)*sin(th4)
which already makes progress because it shows that th2,th3 only come in as their sum.
With d = cos(th2+th3), ans1 is of the form
a - b*d = (1/2)((a+b)*(1-d) + (a-b)*(1+d))
and throwing in some basic trig identities you can arrive at
ans2 = sin(th1+th4)*sin((th2+th3)/2)^2 ...
+ sin(th1-th4)*cos((th2+th3)/2)^2 .
ans2 has a symmetric kind of form but whether it is simpler than ans1 is a matter of opinion. I doubt that further simplification is in the cards.