Three mathematical models were developed to analyze the dynamics of fluorescent dyes penetration into 3D cellular spheroids. Two fluorescent dyes were chosen to verify mathematical models: rhodamine 6G (R6G) as a small molecule, which can freely penetrate through the cells, and wheat germ agglutinin (WGA) conjugated with Alexa488 fluorescent label, which reacts with the cells plasma membrane, and its cellular penetration is significantly lower. Dye penetration and binding to cells were modeled with nonlinear diffusion–reaction equations. System of differential equations was solved using numerical methods, and good correspondence with physical experiment was shown. Diffusion coefficients in extracellular matrix were determined for both fluorescent dyes, and the influence of reactions parameters to WGA penetration was analyzed. Dynamics of dyes accumulation into cell spheroids were also determined.