20∆ 9 -tetrahydrocannabinol (∆ 9 -THC), the main active ingredient of Cannabis sativa 21 (marijuana), interacts with the human brain cannabinoid (CB1) receptor and mimics 22 pharmacological effects of endocannabinoids (eCBs) N-arachidonylethanolamide (AEA) and 2-23 arachidonoylglycerol (2-AG). Given recent intriguing findings that some allosteric modulators 24 can enhance selectively the AEA-activated CB1 receptor, it is imperative to determine the 25 structure of the AEA-bound CB1 receptor. However, due to its highly flexible nature of AEA, 26 establishing its binding mode to the CB1 receptor is elusive. The aim of the present study was to 27 explore many possible binding conformations of AEA within the binding pocket of the CB1 28 receptor that is confirmed in the recently available X-ray crystal structures of the agonist-bound 29 CB1 receptors and predict essential AEA binding domains. We performed long time molecular 30 dynamics stimulations of plausible AEA docking poses until its receptor binding interactions 31 became optimally established. Our simulation results revealed that AEA favors to bind to the 32 hydrophobic channel of the CB1 receptor, suggesting that the hydrophobic channel holds 33 essential significance in AEA binding to the CB1 receptor. Our results also suggest that the 34 H2/H3 region of the CB1 receptor is an AEA binding subsite privileged possibly over the H7 35 region. The results of the present study contribute to identifying the (hidden) allosteric site(s) of 36 the CB1 receptor in our immediate future study.37 38 39 40 41 3 42 4 65 adopt millions of conformations. Identification of the bioactive conformation of AEA to the CB1 66receptor can be quite elusive due to its potential to adopt many low-energy binding 67 conformations only a few of which would be responsible for receptor activation. AEA is a partial 68 agonist at CB1 receptors [16,17] just like Δ 9 -THC but somewhat more potent than Δ 9 -THC in 69 activating the CB1 receptor [1]. Without any known X-ray crystal structure of the AEA-bound 70 CB1 receptor, the nature of binding interactions of AEA with the CB1 receptor remains poorly 71 understood. 72 73 106 AM11542 bound to the CB1 receptor in the X-ray crystal structure of the AM11542-bound CB1 107 receptor [6]. A typical setting of docking parameters for performing AutoDock runs using a 108 hybrid global-local Lamarkian genetic algorithm (LGA) [24] were: the rate of gene mutation 109 (0.02), rate of crossover (0.8), GA window size (10), the number of individuals in population 110 (150), the maximum number of energy evaluations in each run (25,000,000), the maximum 6 111 number of generations (27,000) and the number of LGA docking runs (10). Only the ligand was 112 allowed to freely move inside the grid box while the protein was rigidly fixed in position. The 113 resulting docking poses were evaluated by the AutoDock4 scoring function [25]. AutoDock runs 114were performed more than one hundred times using the best scoring docking pose from the 115 previous run as the star...