Possible radical reaction products issuing from II-atom addition to cytosine have been characterized and analyzed by means of a comprehensive quantum mechanical approach including density functional computations (B3LYP), together with simulation of the solvent by the polarizable continuum model (PCM), and averaging of spectroscopic properties over the most important vibrational motions. The hyperfine couplings of the semirigid 5,6-dihydrocytos-6yl radical computed at the optimized geometry are in good agreement with their experimental counterparts. On the other hand, vibrational averaging is mandatory for obtaining an effectively planar structure for the 5,6-dihydrocytos-5yl radical with the consequent equivalence of beta -hydrogens. Finally, only proper consideration of environmental effects restores the agreement between computed and experimental couplings for the base anion protonated at N3.