Molecular Interactions Governing the Rat Aryl Hydrocarbon Receptor Activities of Polycyclic Aromatic Compounds and Predictive Model Development.

Polycyclic aromatic compounds (PACs) exhibit rat aryl hydrocarbon receptor (rAhR) activities, leading to diverse biological or toxic effects. In this study, the key amino residues and molecular interactions that govern the rAhR activity of PACs were investigated using in silico strategies. The homology model of rAhR was first docked with 90 PACs to yield complexes, and the results of the molecular dynamics simulations of 16 typical complexes showed that the binding energies of the complexes range from -7.37 to -26.39 kcal/mol. The major contribution to the molecular interaction comes from van der Waals forces, and Pro295 and Arg316 become the key residues involved in most complexes. Two QSAR models were further developed to predict the rAhR activity of PACs (in terms of log IEQ for PACs without halogen substitutions and log%-TCDD-max for halogenated PACs). Both models have good predictive ability, robustness, and extrapolation ability. Molecular polarizability, electronegativity, size, and nucleophilicity are identified as the important factors affecting the rAhR activity of PACs. The developed models could be employed to predict the rAhR activity of other reactive PACs. This work provides insight into the mechanisms and interactions of the rAhR activity of PACs and assists in the assessment of their fate and risk in organisms.