Hybrid organic–inorganic semiconductor systems have important applications in both molecular electronics and photoresponsive materials. The characterizations of the interface and of the electronic excited-states of these hybrid systems remain a challenge for state-of-the-art computational methods, as the systems of interest are large. In the present investigation, we present for the first time a many-body Green’s function Bethe–Salpeter investigation of a series of photochromic molecules adsorbed onto TiO2 nanoclusters. On the basis of these studies, the performance of time-dependent density functional theory (TD-DFT) calculations is assessed. In addition, the photochromic properties of different hybrid systems are also evaluated. This work shows that qualitatively different conclusions can be reached with TD-DFT relying on various exchange–correlation functionals for such organic–inorganic interfaces and paves the way to more accurate simulation of many hybrid materials.