The Hydrophobic Ligands Entry and Exit from the GPCR Binding Site-SMD and SuMD Simulations

Most G protein-coupled receptors that bind the hydrophobic ligands (lipid receptors and steroid receptors) belong to the most populated class A (rhodopsin-like) of these receptors. Typical examples of lipid receptors are: rhodopsin, cannabinoid (CB), sphingosine-1-phosphate (S1P) and lysophosphatidic (LPA) receptors. The hydrophobic ligands access the receptor binding site from the lipid bilayer not only because of their low solubility in water but also because of a large N-terminal domain plug preventing access to the orthosteric binding site from the extracellular milieu. In order to identify the most probable ligand exit pathway from lipid receptors CB1, S1P1 and LPA1 orthosteric binding sites we performed at least three repeats of steered molecular dynamics simulations in which ligands were pulled in various directions. For specific ligands being agonists, the supervised molecular dynamics approach was used to simulate the ligand entry events to the inactive receptor structures. For all investigated receptors the ligand entry/exit pathway goes through the gate between transmembrane helices TM1 and TM7, however, in some cases it combined with a direction toward water milieu.