Unraveling <scp>GPCRs</scp> Allosteric Modulation. Cannabinoid 1 Receptor as a Case Study

G-protein-coupled receptors (GPCRs) constitute one of the most prominent families of integral membrane receptor proteins that mediate most transmembrane signaling processes. Malfunction of these signal transduction processes is one of the underlying causes of many human pathologies (Parkinson's, Huntington's, heart diseases, etc), provoking that GPCRs are the largest family of druggable proteins. However, these receptors have been targeted traditionally by orthosteric ligands, which usually causes side effects due to the simultaneous targeting of homologous receptor subtypes. Allosteric modulation offers a promising alternative approach to circumvent this problematic and, thus, comprehending its details is a most important task. Here we use the Cannabinoid type-1 receptor (CB1R) in trying to shed light on this issue, focusing on positive allosteric modulation. This is done by using the protein-dipole Langevin-dipole (PDLD) within the linear response approximation (LRA) framework (PDLD/S-2000) along with our coarse-grained (CG) model of membrane proteins to evaluate the dissociation constants (K _Bs) and cooperativity factors (αs) for a diverse series of CB1R positive allosteric modulators belonging to the 2-phenylindole structural class, considering CP55940 as an agonist. The agreement with the experimental data evinces that significantly populated allosteric modulator:CB1R and allosteric modulator:CP55940:CB1R complexes have been identified and characterized successfully. Analyzing them, it has been determined that CB1R positive allosteric modulation lies in an outwards displacement of transmembrane α helix (TM) 4 extracellular end and in the regulation of the range of motion of a compound TM7 movement for binary and ternary complexes, respectively. In this respect, we achieved a better comprehension of the molecular architecture of CB1R positive allosteric site, identifying Lys192^3.28 and Gly194^3.30 as key residues regarding electrostatic interactions inside this cavity, and to rationalize (at both structural and molecular level) the exhibited stereoselectivity in relation to positive allosteric modulation activity by considered CB1R allosteric modulators. Additionally, putative/postulated allosteric binding sites have been screened successfully, identifying the real CB1R positive allosteric site, and most structure-activity relationship (SAR) studies of CB1R 2-phenylindole allosteric modulators have been rationalized. All these findings point out towards the predictive value of the methodology used in the current work, which can be applied to other biophysical systems of interest. The results presented in this study contribute significantly to understand GPCRs allosteric modulation and, hopefully, will encourage a more thorough exploration of the topic.