The potentiality of the time dependent density functional theory (TDDFT) for the description of core excitation spectra (XAS) in transition metal oxides is analyzed, considering the rutile form of TiO2 as a test case. Cluster models are adopted to mimic the bulk, embedded within an array of point charges to simulate the Madelung potential. All of the edges, titanium and oxygen K and titanium L edges, are considered, and the TDDFT results are compared with the experimental data in order to assess the performance of the theoretical approach in dealing with this complex class of compounds. Satisfactory results have been obtained for the Ti and O K edges, while in the case of the Ti L edge some discrepancies with the experiment are still present. The configuration mixing explicitly included in the TDDFT model strongly influences the distribution of the 2p metal oscillator strength. The origin of the spectral features is investigated with the help of the partial density of the virtual states (PDOS) calculated for each core hole considered, which can be qualitatively compared with the theoretical spectra calculated in the Kohn-Sham one-electron approach.