An extension of the spin density functional theory simultaneously accounting for dielectric mismatch between neighboring materials and nonparabolicity corrections originating from interactions between conduction and valence bands is presented. This method is employed to calculate ground state and addition energy spectra of homogeneous and multishell spherical quantum dots. Our calculations reveal that corrections become especially relevant when they come into play simultaneously in strong regimes of spatial confinement.