Understanding olfaction at the molecular level is challenging due to the lack of crystallographic models of odorant receptors (ORs). To better understand the molecular mechanism of OR activation, we focused on chiral (R)-muscone and other musk-smelling odorants due to their great importance and widespread use in perfumery and traditional medicine, as well as environmental concerns associated with bioaccumulation of musks with estrogenic/antiestrogenic properties. We experimentally and computationally examined the activation of human receptors OR5AN1 and OR1A1, recently identified as specifically responding to musk compounds. OR5AN1 responds at nanomolar concentrations to musk ketone and robustly to macrocyclic sulfoxides and fluorine-substituted macrocyclic ketones; OR1A1 responds only to nitromusks. Structural models of OR5AN1 and OR1A1 based on quantum mechanics/molecular mechanics (QM/MM) hybrid methods were validated through direct comparisons with activation profiles from site-directed mutagenesis experiments and analysis of binding energies for 35 musk-related odorants. The experimentally found chiral selectivity of OR5AN1 to (R)- over (S)-muscone was also computationally confirmed for muscone and fluorinated (R)-muscone analogs. Structural models show that OR5AN1, highly responsive to nitromusks over macrocyclic musks, stabilizes odorants by hydrogen bonding to Tyr260 of transmembrane α-helix 6 and hydrophobic interactions with surrounding aromatic residues Phe105, Phe194, and Phe207. The binding of OR1A1 to nitromusks is stabilized by hydrogen bonding to Tyr258 along with hydrophobic interactions with surrounding aromatic residues Tyr251 and Phe206. Hydrophobic/nonpolar and hydrogen bonding interactions contribute, respectively, 77% and 13% to the odorant binding affinities, as shown by an atom-based quantitative structure-activity relationship model.
Keywords: molecular dynamics; musk; odorant receptor; olfaction; quantum mechanics.