Computer-aided drug design supporting sunscreen research: a showcase study using previously synthesized hybrid UV filter-antioxidant compounds

Abstract

Context: Although quantum mechanical calculations have proven effective in accurately predicting UV absorption and assessing the antioxidant potential of compounds, the utilization of computer-aided drug design (CADD) to support sustainable synthesis research of new sunscreen active ingredients remains an area with limited exploration. Furthermore, there are ongoing concerns about the safety and effectiveness of existing sunscreens. Therefore, it remains crucial to investigate photoprotection mechanisms and develop enhanced strategies for mitigating the harmful effects of UVR exposure, improving both the safety and efficacy of sunscreen products. A previous study conducted synthesis research on eight novel hybrid compounds (I–VIII) for use in sunscreen products by molecular hybridization of trans-resveratrol (RESV), avobenzone (AVO), and octinoxate (OMC). Herein, time-dependent density functional theory (TD-DFT) calculations performed in the gas phase on the isolated hybrid compounds (I–VIII) proved to reproduce the experimental UV absorption. Resveratrol-avobenzone structure-based hybrids (I–IV) present absorption maxima in the UVB range with slight differences between them, while resveratrol-OMC structure-based hybrids (V–VIII) showed main absorption in the UVA range. Among RESV-OMC hybrids, compounds V and VI exhibited higher UV absorption intensity, and compound VIII stood out for its broad-spectrum coverage in our simulations. Furthermore, both in silico and in vitro analyses revealed that compounds VII and VIII exhibited the highest antioxidant activity, with compound I emerging as the most reactive antioxidant within RESV-AVO hybrids. The study suggests a preference for the hydrogen atom transfer (HAT) mechanism over single-electron transfer followed by proton transfer (SET-PT) in the gas phase. With a strong focus on sustainability, this approach reduces costs and minimizes effluent production in synthesis research, promoting the eco-friendly development of new sunscreen active ingredients. Methods: The SPARTAN'20 program was utilized for the geometry optimization and energy calculations of all compounds. Conformer distribution analysis was performed using the Merck molecular force field 94 (MMFF94), and geometry optimization was carried out using the parametric method 6 (PM6) followed by density functional theory (DFT/B3LYP/6-31G(d)). The antioxidant behavior of the hybrid compounds (I–VIII) was determined using the highest occupied molecular orbital (εHOMO) and the lowest unoccupied molecular orbital (εLUMO) energies, as well as the bond dissociation enthalpy (BDE), ionization potential (IP), and proton dissociation enthalpy (PDE) values, all calculated at the same level of structural optimization. TD-DFT study is carried out to calculate the excitation energy using the B3LYP functional with the 6-31G(d) basis set. The calculated transitions were convoluted with a Gaussian profile using the Gabedit program.