Molecular and structural basis of interactions of vitamin D3 hydroxyderivatives with aryl hydrocarbon receptor (AhR): An integrated experimental and computational study

To better understand the molecular and structural basis underlying the interaction of vitamin D3 hydroxyderivatives with AhR, molecular simulation was used to probe the binding of 1,20(OH)₂D3, 1,25(OH)₂D3, 20,23(OH)₂D3 and 20(OH)D3 to AhR. qPCR showed that vitamin D3 derivatives stimulate expression of cyp1A1 and cyp1B1 genes that are downstream targets of AhR signaling. These secosteroids stimulated the translocation of the AhR to the nucleus, as measured by flow cytometry and western blotting. Molecular dynamics simulations were used to model the binding of vitamin D3 derivatives to AhR to examine their influence on the structure, conformation and dynamics of the AhR ligand binding domain (LBD). Binding thermodynamics, conformation, secondary structure, dynamical motion and electrostatic potential of AhR were analyzed. The molecular docking scores and binding free energy were all favorable for the binding of D3 derivatives to the AhR. These established ligands and the D3 derivatives are predicted to have different patterns of hydrogen bond formation with the AhR, and varied residue conformational fluctuations and dynamical motion for the LBD. These changes could alter the shape, size and electrostatic potential distribution of the ligand binding pocket, contributing to the different binding affinities of AhR for the natural ligands and D3 derivatives.