We construct stochastically inhomogeneous epithelial cell models via simulation of Gaussian random fields; the extent and correlation length of subnuclear refractive index fluctuations are based on values quantified from high-resolution images of cervical tissue. We then employ the finite-difference time-domain method to simulate azimuth-resolved light scattering patterns of the constructed models. We process these two-dimensional patterns and calculate a series of Haralick features with the ultimate goal of identifying signatures that directly point to changes in subnuclear refractive index profile. Our results show that azimuthal contrast calculated over specific angular ranges is highly sensitive to the extent and correlation length of refractive index fluctuations. This metric is insensitive to changes in the overall size and mean refractive index of the constructed models, thereby allowing for selective tracking of changes in subnuclear refractive index variations.