To shed light on the nature of the electronic states at play in N-doped TiO2 nanoparticles, we have performed detailed ground and excited state calculations on pure and N-doped TiO2 rutile using an embedding model. We have validated our model by comparing ground-state embedded results with those obtained from periodic DFT calculations. Our results are consistent with periodic calculations. Using this embedding model we have performed B3LYP based TDDFT calculations of the excited state spectrum. We have also studied the lowest excitations using high-level equation-of-motion coupled cluster (EOMCC) approaches involving all single and inter-band double excitations. We compare and contrast the nature of the excitations in detail for the pure and doped systems using these calculations. Our calculations indicate a lowering of the bandgap and confirm the role of the N3- states on the UV/Vis spectrum of N-doped TiO2 rutile supported by experimental findings.