The isomers of the nitro-en-substituted fullerenes (azafullerenes) C19N, C59N, C69N, and C75N are examined using all-electron Gaussian atomic orbital basis density functional theory, to determine the doublet radical geometries and hyperfine coupling constants. We find that the inaccuracy of previously calculated hyperfine coupling constants of C59N resulted from a poor treatment of the geometry optimization. We find that UB3LYP minimization of the radical geometry in the 6-31G* basis, followed by single-point evaluation of the hyperfine constants in which an expanded basis is used on the atomic sites of interest, forms an efficient compromise between computational cost and accuracy with respect to experimental hyperfine constants. Using this approach, we assign the hyperfine signals observed in experiments on the C69N radical by calculating the hyperfine coupling constants for all five of the isomers and examine the electron spin density distribution. Finally, we present predicted hyperfine coupling constants for the isomers of C19N and C75N for use in the interpretation of future experiments.