Intramolecular coevolution of amino acid sites has repeatedly been studied to improve predictions on protein structure and function. Thereby, the focus was on bacterial proteins with available crystallographic data. However, intramolecular coevolution has not yet been compared between protein sets along a gradient of functional proximity to fertilization. This is especially true for the potential effect of external selective forces on intraprotein coevolution. In this study, we investigated both aspects in equally sized sets of mammalian proteins representing spermatozoa, testis, entire body, and liver. For coevolutionary analyses, we derived the proportion of covarying sites per protein from amino acid alignments of 10 mammalian orthologues each. In confirmation of the validity of our coevolution proxy, we found positive associations with the nonsynonymous or amino acid substitution rate in all protein sets. However, our coevolution proxy negatively correlated with the number of protein interactants (node degree) in male reproductive protein