The calculation of optical and spectroscopic properties of chromophores in solution or in complex environments requires methodologies and computational protocols that can account for the intrinsic electronic structure and vibrations of the excited molecule, as well as environmental effects. As mentioned previously, this necessitates reliable sampling of solute-environment configurations. To achieve this, sequential CL/QM methods with varying levels of accuracy have been developed. These approaches can be broadly classified as classical and mixed quantum-classical (MQC) schemes. CL methods account for the effects of solute and environment fluctuations by obtaining the spectrum as the distribution of vertical energies. These can be computed at a reasonable computational cost exploiting the time dependent density functional theory (TDDFT), for example within a QM/MM approach, over an ensemble of frames extracted from a reliable MD trajectory. Such a protocol, which relies on the classical Franck-Condon principle, is referred to in this work as the Classical Ensemble Average of Vertical Energies (CEA-VE) and is used to compute the optical properties of molecule V to VII. Computational details for the CEA-VE method are given in the Supporting Information in Section S4.