I would just approximate $\sin^2$ by 1/2, it oscillates rapidly for any $p\neq 0$ and the average over one period should be a sensible estimate, then 
$$S(p,q,s)\approx \tfrac{1}{2}H_q^{(s)},$$
the Harmonic number. For large $q$ this gives
$$S(p,q,s)\rightarrow q^{-s} \left(\frac{q}{2 (1-s)}+\tfrac{1}{4}+{\cal O}(q^{-2})\right)+\tfrac{1}{2}\zeta (s).$$