You could define 2 (or 3) building blocks such as the Cysteine
dipeptide & thio-phosphate derivative(s) and use inter-molecular
charge constraint(s) (INTER-MCC) during the fit to reconstruct your
phosphorylated cysteine.
R-SH + Me-S-PO3(2-) -> R-S-PO3(2-)
<-------------------->
INTER-MCC sets to zero
+ Me-S-PO3H(-) -> R-S-PO3H(-)
<-------------------->
INTER-MCC sets to zero
To study the impact of the INTER-MCC(s) used R.E.D. Server includes by
now a statistics module.