The present paper provides a mathematical analysis of the model of intracellular protein dynamics proposed in [14]. The model describes protein and mRNA transport inside a cell and takes into account diffusion in the nucleus and cytoplasm as well as active transport of protein molecules along microtubules in the cytoplasm. The model is a complex system of nonlinear PDEs with appropriate boundary conditions. The model reproduces, at least in numerical simulations, the oscillatory changes in protein concentration observed in the experimental data. To our knowledge this is the first paper that, in the multidimensional case, deals with a rigorous mathematical analysis of a model of intracellular dynamics with active transport on microtubules. In particular, in the present paper, we prove the existence and uniqueness result for the model in arbitrary space dimension. The model may be adapted to other signaling pathways.