The glucocorticoid receptor (GR) is a ligand-binding dependent transcription factor that ultimately regulates vital biological processes and inflammation response through specific gene expression control, thus representing a notable drug target to explore. Structurally, its ligand binding domain (LBD) harbors the region for the ligand-dependent transcriptional activation function 2 (AF-2), a majorly hydrophobic groove formed by residues from helices H3, H4, and H12, where the H12 position plays a critical role in AF-2 spatial conformation and GR function as a whole. However, the exact mechanisms underlying how regulatory ligands control the H12 structure and dynamics are yet to be elucidated. In this work, we have explored the correlation between ligand identity and GR LBD H12 behavior through different molecular dynamics (MD) simulations. After building diverse GR LBD systems in agonist and nonagonist states, we studied each system's response in the absence or the presence of an agonist ligand (dexamethasone) or an antagonist ligand (RU486) using classical MD simulations. We complemented them with steered MD for assessing the transition between those states and with the Umbrella Sampling method for free-energy evaluation. On the one hand, successfully obtaining fully folded nonagonist GR LBD states from the partially unfolded crystal GR LBD/RU486 underlines the role of the H1 in the GR LBD folding pathway. On the other hand, our results describe the H12 as a dynamic ensemble of conformations whose relative population is in the end determined by the interacting ligand: while dexamethasone privileges only a few poses (determined by a potential energy surface with a deep minimum), RU486 favors a wider H12 conformational amplitude, as indicated by a flatter potential landscape. By characterizing the H12 conformation in different conditions, we provide novel GR LBD models that represent potential targets for rational glucocorticoid drugs design, with the aim of accurately modulating GR activity.

The recognition of Cannabis as a source of new compounds suitable for medical use has attracted strong interest from the scientific community in its research, and substantial progress has accumulated regarding cannabinoids' activity; however, a thorough description of their molecular mechanisms of action remains a task to complete. Highlighting their complex pharmacology, the list of cannabinoids' interactors has vastly expanded beyond the canonical cannabinoid receptors. Among those, we have focused our study on the glycine receptor (GlyR), an ion channel involved in the modulation of nervous system responses, including, to our interest, sensitivity to peripheral pain. Here, we report the use of computational methods to investigate possible binding modes between the GlyR and Δ9‐tetrahydrocannabinol (THC). After obtaining a first pose for the THC binding from a biased molecular docking simulation and subsequently evaluating it by molecular dynamic simulations, we found a dynamic system with an identifiable representative binding mode characterized by the specific interaction with two transmembrane residues (Phe293 and Ser296). Complementarily, we assessed the role of membrane cholesterol in this interaction and positively established its relevance for THC binding to GlyR. Lastly, the use of restrained molecular dynamics simulations allowed us to refine the description of the binding mode and of the cholesterol effect. Altogether, our findings contribute to the current knowledge about the GlyR‐THC mode of binding and propose a new starting point for future research on how cannabinoids in general, and THC in particular, modulate pain perception in view of its possible clinical applications.