In spectroscopy simulations using time-dependent DFT, calculations have two distinct parts where the first involves the SCF optimization of the electron density of the unperturbed reference state and the second involves the application of response theory to determine the time-dependent oscillations in the electron density as due to the external electromagnetic fields. The overarching goal of our software development is to minimize the total wall time required to produce final spectra of realistic model systems. For biochemical systems in particular, that requires attention to be paid to both the environment (often aqueous solutions and/or protein embeddings) and the effects of molecular dynamics (typically handled by configuration sampling). The former of these two aspects is here taken into account with use of QM/MM modeling whereas the latter aspect in not addressed. Configuration space sampling is, however, embarrassingly parallel in nature and can therefore be left out from the discussion of computational performance and scaling. From a physical standpoint, it is also well motivated to leave out this issue in the case of our specific application as the G-quadruplex has a very low degree of conformational flexibility.