Spatio-temporal dynamics of a variety of proteins is, among other things, regulated by posttranslational modifications of these proteins. Such modifications can thus influence stability and biochemical activities of the proteins, activity and stability of their upstream targets within specific signalling pathways. Commonly used mathematical tools for such protein-protein (and/or protein-mRNA) interactions in single cells, namely, Michaelis-Menten and Hill kinetics, yielding a system of ordinary differential equations, are extended here into (non-linear) partial differential equations by taking into account a more realistic spatial representation of the environment where these reactions occur. In the modelling framework under consideration, all interactions occur in a cell divided into two compartments, the nucleus and the cytoplasm, connected by the semipermeable nuclear membrane and bounded by the impermeable cell membrane. Passive transport mechanism, modelled by the so-called Kedem-Katchalsky boundary conditions, is used here to represent migration of species throughout the nuclear membrane. Nonlinear systems of partial differential equations are solved by the semi-implicit Rothe method. Examples of two spatial oscillators are shown. Namely, these are the circadian rhythm for concentration of the FRQ protein in Neurospora crassa and oscillatory dynamics observed in the activation and regulation of the p53 protein following DNA damage.