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.