Aflatoxin B₁ (AFB₁) exhibits the most potent mutagenic and carcinogenic activity among aflatoxins. For this reason, AFB₁ is recognized as a human group 1 carcinogen by the International Agency of Research on Cancer. Consequently, it is essential to determine its properties and behavior in different chemical systems. The chemical properties of AFB₁ can be explored using computational chemistry, which has been employed complementarily to experimental investigations. The present review includes in silico studies (semiempirical, Hartree–Fock, DFT, molecular docking, and molecular dynamics) conducted from the first computational study in 1974 to the present (2022). This work was performed, considering the following groups: (a) molecular properties of AFB₁ (structural, energy, solvent effects, ground and the excited state, atomic charges, among others); (b) theoretical investigations of AFB₁ (degradation, quantification, reactivity, among others); (c) molecular interactions with inorganic compounds (Ag⁺, Zn²⁺, and Mg²⁺); (d) molecular interactions with environmentally compounds (clays); and (e) molecular interactions with biological compounds (DNA, enzymes, cyclodextrins, glucans, among others). Accordingly, in this work, we provide to the stakeholder the knowledge of toxicity of types of AFB₁-derivatives, the structure–activity relationships manifested by the bonds between AFB₁ and DNA or proteins, and the types of strategies that have been employed to quantify, detect, and eliminate the AFB₁ molecule.

Theoretical studies of 8-chloro-9-hydroxy-aflatoxin B₁ (2) were carried out by Density Functional Theory (DFT). This molecule is the reaction product of the treatment of aflatoxin B₁ (1) with hypochlorous acid, from neutral electrolyzed water. Determination of the structural, electronic and spectroscopic properties of the reaction product allowed its theoretical characterization. In order to elucidate the formation process of 2, two reaction pathways were evaluated—the first one considering only ionic species (Cl⁺ and OH⁻) and the second one taking into account the entire hypochlorous acid molecule (HOCl). Both pathways were studied theoretically in gas and solution phases. In the first suggested pathway, the reaction involves the addition of chlorenium ion to 1 forming a non-classic carbocation assisted by anchimeric effect of the nearest aromatic system, and then a nucleophilic attack to the intermediate by the hydroxide ion. In the second studied pathway, as a first step, the attack of the double bond from the furanic moiety of 1 to the hypochlorous acid is considered, accomplishing the same non-classical carbocation, and again in the second step, a nucleophilic attack by the hydroxide ion. In order to validate both reaction pathways, the atomic charges, the highest occupied molecular orbital and the lowest unoccupied molecular orbital were obtained for both substrate and product. The corresponding data imply that the C₉ atom is the more suitable site of the substrate to interact with the hydroxide ion. It was demonstrated by theoretical calculations that a vicinal and anti chlorohydrin is produced in the terminal furan ring. Data of the studied compound indicate an important reduction in the cytotoxic and genotoxic potential of the target molecule, as demonstrated previously by our research group using different in vitro assays.

Several indolylquinone analogues of perezone, a natural sesquiterpene quinone, were characterized in this work by theoretical methods. In addition, some physicochemical, toxicological and metabolic properties were predicted using bioinformatics software. The predicted physicochemical properties are in agreement with the solubility and cLogP values, the penetration across the cell membrane, and absorption values, as well as with a possible apoptosis-activated mechanism of cytotoxic action. The toxicological predictions suggest no mutagenic, tumorigenic or reproductive effects of the four target molecules. Complementarily, the results of a performed docking study show high scoring values and hydrogen bonding values in agreement with the cytotoxicity IC₅₀ value ranking, i.e: indolylmenadione > indolylperezone > indolylplumbagine > indolylisoperezone. Consequently, it is possible to suggest an appropriate apoptotic pathway for each compound. Finally, potential metabolic pathways of the molecules were proposed.

The structure of the ortho-, meta- and para- hybrid diindolylmethane-phenylboronic acids and their interactions were optimized with by a quantum chemical method, using density functional theory at the (DFT) level. Thus, infrared bands were assigned based on the scaled theoretical wavenumbers by correlating the respective experimental data of the molecules. In addition, the corresponding ¹H-/¹³C-/¹¹B-NMR experimental and theoretical chemical shifts were correlated. The target molecules showed a poor treatment of the OH shifts in the GIAO method due to the absence of explicit solvent effects in these calculations; therefore, they were explicitly considered with acetone molecules. Moreover, the electron density at the hydrogen bond critical point increased, generating stabilization energy, from weak to moderate or weak to strong, serving as an indicator of the strength of the hydrogen bond between the different intermolecular interactions. Finally, some properties related to the reactive behavior of the target molecules associated with their cytotoxic effects and metabolic pathways were also calculated.