In the present work a general systems biology approach has been used to study the complex regulatory network controlling the transcription of the spa gene, encoding protein A, a major surface protein and an important virulence factor of Staphylococcus aureus. A valid mathematical model could be formulated using parameter values, which were fitted to quantitative Northern blot data from various S. aureus regulatory mutants using a gradient search method. The model could correctly predict spa expression levels in 4 different regulatory mutants not included in the parameter value search, and in 2 other S. aureus strains, SH1000 and UAMS-1. The mathematical model revealed that sarA and sarS seem to balance each other in a way that when the activating impact of sarS is small, e.g. in the wild-type, the repressive impact of sarA is small, while in an agr-deficient background, when the impact of sarS is maximal, the repressive impact of sarA is close to its maximum. Furthermore, the model revealed that Rot and SarS act synergistically to stimulate spa expression, something that was not obvious from experimental data. We believe that this mathematical model can be used to evaluate the significance of other putative interactions in the regulatory network governing spa transcription.