Aims. We investigate the gravitational redshift in clusters of galaxies inthe symmetron and Hu-Sawicky $f(R)$ models. The characteristic feature of bothmodels is the screening mechanism that hides the fifth force in denseenvironments recovering general relativity. Methods. We use $N$-bodysimulations that were run using the code \texttt{Isis}, which includes scalarfields, to analyse the deviation of observables in modified gravity models withrespect to $\Lambda$CDM. Results. We find that due to the presence of thescreening, the deviation is highly dependent on the halo mass. For instance,the $f(R)$ parameters $|f_{R0}|=10^{-5}$, $n=1$ cause an enhancement of thegravitational signal by up to 50% for halos with masses between$10^{13}\,M_\odot h^{-1}$ and $10^{14}\,M_\odot h^{-1}$. The characteristicmass range where the fifth force is most active varies with the modelparameters. The usual assumption is that the presence of a fifth force leads toa deeper potential well and thus, a stronger gravitational redshift. However,we find that in cases in which only the central regions of the halos arescreened, there could also be a weaker gravitational redshift.