Allosteric antagonist binding sites in class B GPCRs: corticotropin receptor 1

The 41 amino acid neuropeptide, corticotropin-releasing factor (CRF) and its associated receptors CRF(1)-R and CRF(2)-R have been targeted for treating stress related disorders. Both CRF(1)-R and CRF(2)-R belong to the class B G-protein coupled receptors for which little information is known regarding the small molecule antagonist binding characteristics. However, it has been shown recently that different non-peptide allosteric ligands stabilize different receptor conformations for CRF(1)-R and hence an understanding of the ligand induced receptor conformational changes is important in the pharmacology of ligand binding. In this study, we modeled the receptor and identified the binding sites of representative small molecule allosteric antagonists for CRF(1)-R. The predicted binding sites of the investigated compounds are located within the transmembrane (TM) domain encompassing TM helices 3, 5 and 6. The docked compounds show strong interactions with H228 on TM3 and M305 on TM5 that have also been implicated in the binding by site directed mutation studies. H228 forms a hydrogen bond of varied strengths with all the antagonists in this study and this is in agreement with the decreased binding affinity of several compounds with H228F mutation. Also mutating M305 to Ile showed a sharp decrease in the calculated binding energy whereas the binding energy loss on M305 to Leu was less significant. These results are in qualitative agreement with the decrease in binding affinities observed experimentally. We further predicted the conformational changes in CRF(1)-R induced by the allosteric antagonist NBI-27914. Movement of TM helices 3 and 5 are dominant and generates three degenerate conformational states two of which are separated by an energy barrier from the third, when bound to NBI-27914. Binding of NBI-27914 was predicted to improve the interaction of the ligand with M305 and also enhanced the aromatic stacking between the ligand and F232 on TM3. A virtual ligand screening of ~13,000 compounds seeded with ~350 CRF(1)-R specific active antagonists performed on the NBI-27914 stabilized conformation of CRF(1)-R yielded a 44% increase in enrichment compared to the initially modeled receptor conformation at a 10% cutoff. The NBI-27914 stabilized conformation also shows a high enrichment for high affinity antagonists compared to the weaker ones. Thus, the conformational changes induced by NBI-27914 improved the ligand screening efficiency of the CRF(1)-R model and demonstrate a generalized application of the method in drug discovery.