The spectroscopic properties of the organic chromophore 4-naphthoyloxy-1-methoxy-2,2,6,6-tetramethylpiperidine (NfO-TEMPO-Me) in toluene solution are explored through an integrated computational strategy combining a classical dynamic sampling with a quantum mechanical description within the framework of the time-dependent density functional theory (TDDFT) approach. The atomistic simulations are based on an accurately parametrized force field, specifically designed to represent the conformational behavior of the molecule in its ground and bright excited states, whereas TDDFT calculations are performed through a selected combination of hybrid functionals and basis sets to obtain optical spectra closely matching the experimental findings. Solvent effects, crucial to obtain good accuracy, are taken into account through explicit molecules and polarizable continuum descriptions. Although, in the case of toluene, specific solvation is not fundamental, the detailed conformational sampling in solution has confirmed the importance of a dynamic description of the molecular geometry for a reliable description of the photophysical properties of the dye. The agreement between theoretical and experimental data is established and a robust protocol for the prediction of the optical behaviour of flexible fluorophores in solution is set.