Statistical physics modeling and interpretation of the adsorption of enantiomeric terpenes onto the human olfactory receptor OR1A1

Decoding Molecular Mechanism Underlying Human Olfactory Receptor OR8D1 Activation by Sotolone Enantiomers

Sotolone, a chiral compound, plays an important role in the food industry. Herein, (R)-/(S)-sotolone were separated to determine their odor characteristics and thresholds in air (R-form: smoky, burned, herb, and green aroma, 0.0514 μg/m3; S-form: sweet, milk, acid, and nutty aroma, 0.0048 μg/m3). OR8D1 responses to (R)-/(S)-sotolone were detected in a HEK293 cell-based luminescence assay. (S)-Sotolone was a more potent agonist than (R)-sotolone (EC50 values of 84.98 ± 1.05 and 167.20 ± 0.25 μmol/L, respectively). Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area analyses confirmed that the combination of (S)-sotolone and OR8D1 was more stable than that of (R)-sotolone. Odorant docking, multiple sequence alignments, site-directed mutagenesis, and functional studies with recombinant odorant receptors (ORs) in a cell-based luminescence assay identified 11 amino-acid residues that influence the enantioselectivity of OR8D1 toward sotolone significantly and that N2065.46 was indispensable to the activation of OR8D1 by (S)-sotolone.

Quantitative investigations of Zebrafish olfactory receptor ORA1 responsiveness to three pheromones: Microscopic and macroscopic characterizations via an advanced statistical physics treatment

In this work, an adsorption phenomenon putatively involved in the olfactory sense of phenylacetic acid, 4-chlorophenylacetic acid, and 4-methoxyphenylacetic acid pheromones in the Zebrafish olfactory receptor ORA1 was a helpful mechanism in interpreting and characterizing the olfaction process at a molecular level. Hence, the experimental dose-olfactory response curves were fitted by applying the one-layer adsorption model with a single energy (1LM1E). On one hand, the different parameters introduced in the selected model were used to microscopically study the three olfactory systems. Indeed, the fitting results showed that phenylacetic acid displayed the greatest maximum olfactory response at saturation, due to the effect of functional groups at the R4 position. The three pheromones were docked via a non-parallel orientation and the adsorption process was a multi-molecular mechanism. The sizes of different binding pockets of ORA1 were determined through the estimation of the olfactory receptor site size distributions (stereographic characterization). The estimated adsorption energies, ranging from 17.340 to 21.332 kJ/mol, can be used to describe the energetic interactions between the studied pheromones and the Zebrafish ORA1 binding pockets. The spectrums of the adsorption energy distributions of phenylacetic acid, 4-chlorophenylacetic acid, and 4-methoxyphenylacetic acid, which were spread out from 10 to 32.5 kJ/mol, 5 to 30 kJ/mol, and 10 to 32.5 kJ/mol, respectively, was determined to estimate the corresponding olfactory bands (energetic characterization). On the other hand, three thermodynamic functions were estimated in order to macroscopically study the three olfactory systems.

Quantitative characterizations of mOR-EG activated by vanilla odorants using advanced statistical physics modeling

An advanced monolayer adsorption model of an ideal gas was successfully employed to investigate the adsorption of vanillin, vanillin methyl ether, vanillin ethyl ether, and vanillin acetate odorants on mouse eugenol olfactory receptor mOR-EG. In order to understand the adsorption process putatively introduced in olfactory perception, model parameters were analyzed. Hence, fitting results showed that the studied vanilla odorants were linked in mOR-EG binding pockets with a non-parallel orientation, and their adsorption was a multi-molecular process (n > 1). The adsorption energy values that ranged from 14.021 to 19.193 kJ/mol suggested that the four vanilla odorants were physisorbed on mOR-EG (ΔE^a < 40 kJ/mol) and the adsorption mechanism may be considered as an exothermic mechanism (ΔE^a > 0). The estimated parameters may also be utilized for the quantitative characterization of the interactions of the studied odorants with mOR-EG to determine the corresponding olfactory bands ranging from 8 to 24.5 kJ/mol.

Advanced investigation of a putative adsorption process of nine non key food odorants (non-KFOs) on the broadly tuned human olfactory receptor OR2W1: Statistical physics modeling and molecular docking study

In the present paper, statistical physics formalism was used to understand the olfactory perception via the investigation of dose-olfactory response curves of a putative adsorption process of nine non key food odorants (non-KFOs) on the broadly tuned human olfactory receptor OR2W1, in order to quantitative characterize the interactions between the nine studied non-KFOs, i. e., furfuryl sulfide, furfuryl disulfide, benzyl methyl disulfide, furfuryl methyl disulfide, benzyl methyl sulfide, 1-phenylethanethiol, benzyl mercaptan, furfuryl methyl sulfide and 3-phenylpropanol molecules and OR2W1 binding sites at a molecular level. Two advanced adsorption models have been proposed: the advanced monolayer monoenergy model (monolayer model with identical and independent olfactory receptor binding sites) (Model 1) and the advanced monolayer model with two independent types of olfactory receptor binding sites (Model 2). It was concluded that the monolayer monoenergy model was selected as the most adequate model to fit the experimental dose-olfactory response curves tabulated in literature. Actually, the numerical values of the three fitted physico-chemical parameters (RM1, n and C1) were obtained by a non-linear regression. Indeed, modeling results suggested that the number of docked non-KFOs per OR2W1 binding site n values (1.24 < n < 1.94) was always superior to 1, which indicated the non-parallel orientation of the studied odorants on the olfactory receptor and the multi-molecular adsorption mechanism. The estimated molar adsorption energy ΔEa values (ranged from 6.07 to 12.16 kJ/mol) for the nine olfactory systems confirmed the physical the exothermic characters of the adsorption process since ΔEa values were lower than 40 kJ/mol and positive. Furthermore, these estimated parameters were applied to characterize stereographically and energetically the interaction between the nine non-KFOs and OR2W1 through the determination of the human receptor binding site size distributions (RSDs) and the adsorption energy distributions (AEDs), which were spread out from 0.25 to 6.50 nm and from 0 to 22.50 kJ/mol, respectively. The docking computation between these nine non-KFOs and OR2W1 proved that the estimated binding affinities were belonged to the adsorption energies spectrum in general and the specific adsorption energy band or the molecular vibration modes limited spectrum (between 2.50 kJ/mol and 17 kJ/mol) (approximate olfactory band).

New insights on the adsorption of floral odorants on Apis cerana cerana olfactory receptor AcerOr1: Theoretical modeling and thermodynamic study

Apis cerana cerana counted on its sensitive olfactory system to make survival activities in the surrounding environment and the olfactory receptors can be considered as a primary requirement of odorant detection, recognition and coding. Indeed, the exploitation of the olfactory system of insects in particular the Asian honeybee "Apis cerana cerana" can be the best experimental model to investigate the essentials of the chemosensitivity and may help to better understand the olfactory perception in insects. Hence, an advanced statistical physics modeling via the monolayer model with single energy (n ≠ 1) of the three dose-olfactory responses curves indicated that undecanoic acid, 1-octyl alcohol and 1-nonanol were docked with a mixed parallel and non-parallel orientation on AcerOr1. Furthermore, in the present work, the Apis cerana cerana olfactory receptor AcerOr1 showed high sensitivity and discrimination power to detect undecanoic acid, 1-octyl alcohol and 1-nonanol with concentrations at half saturations values of 10^-7 mol/L and the molar adsorption energy values obtained from data fitting results, which were ranged from 17.91 to 24.00 kJ/mol, confirmed the exothermic and the physisorption nature of the adsorption of the studied floral odorants on AcerOr1. The studied experimental dose-response curves of undecanoic acid, 1-octyl alcohol and 1-nonanol provided access to quantitative (i.e., stereographic and energetic) characterizations of AcerOr1 via the determination of the olfactory receptor site size distributions (RSDs) and the adsorption energy distributions (AEDs). The stereographic characterization showed RSDs spread out from 0.20 to 8 nm presenting average values corresponding to the maximum of the peaks at 1.50 nm, at 1.10 nm and at 1.04 nm for undecanoic acid, 1-octyl alcohol and 1-nonanol, respectively. The energetic characterization presented AEDs ranged from 0 to 40 kJ/mol showing an approximate adsorption energy bands defined between 7.50 and 27.50 kJ/mol, between 15 and 33 kJ/mol and between 13.50 and 34.50 kJ/mol for undecanoic acid, 1-octyl alcohol and 1-nonanol, respectively. The utilization of the analytical expression of the olfactory threshold allowed giving important and helpful informations about the occupation rate of AcerOr1 binding sites that fired a minimal olfactory response at a honeybee olfactory receptor. Hence, the olfactory response can be detected only when 1.97 %, 1.13 % and 2.00 % of AcerOr1 binding sites were occupied by undecanoic acid, 1-octyl alcohol and 1-nonanol, respectively. Lastly, by means of the selected model, the thermodynamic potentials, such as the adsorption entropy, the Gibbs free enthalpy and the internal energy could be calculated and interpreted.

Advanced investigation of the olfactory perception of semiochemical TMT on OR5K1 and Olfr175 by statistical physics approach

The adsorption of the trimethylthiazoline (TMT) on the human olfactory receptor OR5K1 and the mouse olfactory receptor Olfr175 was the object of the present paper. The main contribution of this work was to characterize stereographically and energetically OR5K1 and Olfr175 activated by trimethylthiazoline molecules docked on the human and the mouse olfactory binding pockets using the grand canonical ensemble in statistical physics. The experimental data and the advanced statistical physics models revealed that the adsorption of the trimethylthiazoline on the human olfactory receptor OR5K1 can be interpreted using the monolayer model with single energy, while the monolayer model with two energies described the interaction between the trimethylthiazoline molecules and the mouse olfactory receptor Olfr175. In fact, the investigated odorant was shown to be docked by a multi-docking process and non parallel orientation on OR5K1 and Olfr175 since the values of the number of TMT molecules per binding site n were superior to 1. The proposed models were applied to calculate the human and the mouse olfactory receptor binding site size distributions relative to TMT, which were spread out from 0.30 to 20 nm with a maximum at about 1.75 nm for OR5K1 and from 1 to 25 nm with a peak at about 4.25 nm for Olfr175. Furthermore, it was found from the calculated molar adsorption energies, which were lower than 11 kJ/mol, that physical adsorption process was occurred in the two olfactory systems. The adsorption energy distributions relative to TMT can be also calculated in order to understand of olfaction process in general through the determination of olfactory bands (i. e., adsorption energy distribution bands), which were situated between 0 and 10.50 kJ/mol and between 3 and 12.50 kJ/mol for OR5K1 and Olfr175, respectively. Referring to the investigation of thermodynamic functions governing the adsorption process such as the adsorption entropy, the Gibbs free enthalpy and the internal energy, it may be noted that the disorder peak of the two olfactory systems was reached when the equilibrium concentration was equal to the concentration at half saturation. In addition, the Gibbs free enthalpy and the internal energy were calculated and their negative values indicated that the adsorption process involved in the olfactory mechanism was exothermic and spontaneous nature.

Theoretical study of the olfactory perception of floral odorant on OR10J5 and Olfr16 using the grand canonical ensemble in statistical physics approach

In this work, two experimental dose-response curves of lyral molecules on the OR10J5 and the Olfr16 were employed in order to examine the evolution of physico-chemical parameters involved in the selected statistical physics model(s) to investigate the human and the mouse smelling of a floral scent. Indeed, one layer adsorption model on one type of sites with one energy (1LAM1T1E) and one layer adsorption model on two types of sites with two energies (1LAM2T2E), considered as appropriate models for the adsorption of lyral molecules on the OR10J5 and Olfr16, respectively, have been applied to fit the experimental data. Stereographic and energetic physico-chemical parameters, namely: the maximum response(s) at saturation, the number of docked molecules per olfactory receptor binding site and the concentration(s) at half saturation, were investigated to retrieve helpful information to describe the adsorption process putatively introduced in the olfaction perception. Thus, the advanced modeling results indicated that the studied molecules were docked with a non-parallel orientation (n > 1). Furthermore, for the two olfactory systems, the molar adsorption energies estimated from curves modeling were inferior to 11 kJ/mol, which showed the physisorption process of the adsorption of lyral molecules on OR10J5 and Olfr16. The 1LAM2T2E and the 1LAM1T1E were applied to estimate the OR10J5 and the Olfr175 RSDs, respectively. Hence, lyral RSDs were spread out from 0.7 to 20 nm with maximums at about 4 nm for OR10J5 and at about 3.65 nm for Olfr16. In addition, by using the two advanced models, the olfactory responses of lyral on OR10J5 and Olfr16 can be used for the energetic characterization of the lyral-OR10J5/Olfr16 binding sites interactions and allowed access to the adsorption energy distributions (AEDs). Then, two approximate olfactory bands can be determined for lyral molecules docked on OR10J5 and Olfr16, which are defined between 3 and 15.5 kJ/mol and between 3.5 and 13.5 kJ/mol, respectively. Lastly, thanks to the proposed models the adsorption entropy of the studied systems can be calculated to describe the disorder and the order on OR10J5 and Olfr16 surfaces (disorder peak of the two olfactory systems was attained when the equilibrium concentration was equal to the concentration at half saturation). Furthermore, the Gibbs free enthalpy and the internal energy were estimated and their negative values indicated that the adsorption phenomenon involved in the olfactory perception was spontaneous and exothermic nature.

Selectively Fluorinated Citronellol Analogues Support a Hydrogen Bonding Donor Interaction with the Human OR1A1 Olfactory Receptor

C-2 fluorinated and methylated stereoisomers of the fragrance citronellol 1 and its oxalate esters were prepared from (R)-pulegone 11 and explored as agonists of the human olfactory receptor OR1A1 and assayed also against site-specific mutants. There were clear isomer preferences and C-2 difluorination as in 18 led to the most active compound suggesting an important hydrogen bond donor role for citronellol 1. C-2 methylation and the corresponding oxalate ester analogues were less active.

Physico-chemical investigations of human olfactory receptors OR10G4 and OR2B11 activated by vanillin, ethyl vanillin, coumarin and quinoline molecules using statistical physics method

This research work is a contribution to understand the olfaction mechanism at a molecular level of vanillin, ethyl vanillin, coumarin and quinoline molecules using a modeling of a putative adsorption process by analytical model established by statistical physics formalism. A statistical physics modeling using the monolayer model with identical and independent binding sites of the responses of the two human olfactory receptors OR10G4 and OR2B11 showed that vanillin and quinoline were adsorbed with a mixed non-parallel and parallel orientation on OR10G4 and on OR2B11, respectively. However, ethyl vanillin and coumarin were anchored with a total non-parallel orientation. The adsorption energy values collected from data analysis, which were ranged from 12.51 to 20.91 kJ/mol, confirmed that the adsorption of vanillin and ethyl vanillin on OR10G4 and the adsorption of coumarin and quinoline on OR2B11were exothermic and were based on physical interactions. Furthermore, the dose-olfactory response curves of vanillin, ethyl vanillin, coumarin and quinoline provided access to OR10G4 and OR2B11 steric characterization via the calculation of the studied olfactory receptors site size distributions (RSDs). Indeed, vanillin, ethyl vanillin, coumarin and quinoline RSDs are spread from 0.3 to 12 nm, from 0.5 to 12 nm, from 0.40 to 12 nm and from 0.14 to 12 nm, respectively, with a maximum at 1.55 nm, 2.11 nm, 2.50 and 1.13 nm, respectively. Lastly, the physico-chemical model parameters can be used for the energetic characterization to confirm the physical nature of the vanillin/ethyl vanillin-OR10G4 and the coumarin/quinoline-OR2B11 interactions and to determine an olfactory band of order of 12 kJ/mol [11-23 kJ/mol], 10 kJ/mol [14-24 kJ/mol], 7 kJ/mol [9-16 kJ/mol], 15 kJ/mol [13-28 kJ/mol] for vanillin, ethyl vanillin, coumarin and quinoline, respectively, through the determination of the adsorption energy distributions (AEDs).

Steric and energetic characterizations of mouse and human musk receptors activated by nitro musk smelling compounds at molecular level: Statistical physics treatment and molecular docking analysis

Understanding olfaction process at a microscopic or molecular level needs more elucidation of the multiple stages involved in the olfaction mechanism. A worth full elucidation and a better understanding of this molecular mechanism, a necessary preamble should be achieved. The content of this work is a preamble for that. A study of the mouse and human olfactory receptors activation in response to two nitro musks stimuli, which are the musk xylol and the musk ketone, are considered here, first, for their wide expanded use in perfumery, but also to show some particular aspects of this process in the case of these two stimuli, which could help to deduce more details and more general aspects in the global olfactory mechanism. A statistical physics modeling using the monolayer model with two independent types of receptor binding sites of the response of the mouse olfactory receptor MOR215-1 and the human olfactory receptor OR5AN1, which are identified as specifically responding to musk compounds, is used to characterize the interaction between the two nitro musk molecules, the mouse and the human olfactory receptors and to determine the olfactory band of these two odorants through the determination of the molar adsorption energies and the adsorption energy distributions. The physico-chemical model parameters can be used for the steric characterization via the calculation of the receptor site size distributions. The docking computation between these two nitro musks and the human olfactory receptor OR5AN1 is performed demonstrating a large similarity in receptor-ligand detection process. Thus, docking finding results prove that the calculated binding affinities were belonging to the spectrum of adsorption energies.