A Computational Approach to the Study of the Binding Mode of S1PR Agonists Based on the Active-Like Receptor Model.

Sphingosine-1-phosphate receptor 1 (S1PR), a member of the G protein-coupled receptor (GPCR) family, is an attractive protein target for the treatment of autoimmune diseases, and a diverse array of S1PR agonists have been developed. Rational drug design based on S1PR remains challenging due to the limited information available on the binding mode between S1PR and its agonists. In this work, the active-like state of S1PR was modeled via Gaussian accelerated molecular dynamics (GaMD) based on its inactive form, which was further validated by docking studies with two representative S1PR agonists. Moreover, with the usage of the induced active-like state, the binding mode between S1PR and its agonists was studied through molecular dynamics simulations and MM-GBSA calculations. The results of those studies indicated that four groups of binding site residues were the major contributors to the ligand and receptor interactions. In addition, this model was verified by five chemically similar compounds synthesized in-house and 1145 known S1PR agonists collected from the BindingDB database. The elucidation of the key binding characteristics will further complete the cognition of S1PR, which can guide the rational design of novel S1PR agonists.