Glucocorticoid Receptor Ligand Binding Domain Research Articles

Screening the binding affinity of bile acid derivatives for the glucocorticoid receptor ligand-binding domain

The necessity of anti-inflammatory drugs such as glucocorticoids has been evident during the COVID-19 pandemic. Glucocorticoids, are the standard therapy for the treatment of moderate and severe COVID-19 patients. However, serious side effects limit the use of these drugs, and anti-inflammatory drugs with better pharmacological properties are urgently required. Bile acids are of interest, because of their anti-inflammatory and immunomodulatory properties, facilitated through an unclear mechanism involving transmembrane and nuclear receptors. In this work, we screened the binding activity of a number of bile acid derivatives, for the ligand-binding domain of glucocorticoid receptor (GR-LBD), the most important receptor for anti-inflammatory processes. Tested compounds include oximes, lactones, lactams, tetrazoles, dienones, C-24 alcohols and cholic acid amides. Cholic acid oxime, deoxycholic acid dienone, 3-keto-24-cholic alcohol and cholic acid amide showed best binding affinities for GR-LBD among tested compounds. The in silico molecular docking explanation is provided. SAR analysis showed that expansion of B and C steroid rings or attachment of heterocycle to C ring is not beneficial for binding; side chain should contain hydrogen donor group; the GR-LBD tolerate well different functionalities on C-3 position. These results provide valuable information toward synthesis of the new glucocorticoids based on bile acids.

The Hippocampal Response to Acute Corticosterone Elevation Is Altered in a Mouse Model for Angelman Syndrome.

Angelman Syndrome (AS) is a severe neurodevelopmental disorder, caused by the neuronal absence of the ubiquitin protein ligase E3A (UBE3A). UBE3A promotes ubiquitin-mediated protein degradation and functions as a transcriptional coregulator of nuclear hormone receptors, including the glucocorticoid receptor (GR). Previous studies showed anxiety-like behavior and hippocampal-dependent memory disturbances in AS mouse models. Hippocampal GR is an important regulator of the stress response and memory formation, and we therefore investigated whether the absence of UBE3A in AS mice disrupted GR signaling in the hippocampus. We first established a strong cortisol-dependent interaction between the GR ligand binding domain and a UBE3A nuclear receptor box in a high-throughput interaction screen. In vivo, we found that UBE3A-deficient AS mice displayed significantly more variation in circulating corticosterone levels throughout the day compared to wildtypes (WT), with low to undetectable levels of corticosterone at the trough of the circadian cycle. Additionally, we observed an enhanced transcriptomic response in the AS hippocampus following acute corticosterone treatment. Surprisingly, chronic corticosterone treatment showed less contrast between AS and WT mice in the hippocampus and liver transcriptomic responses. This suggests that UBE3A limits the acute stimulation of GR signaling, likely as a member of the GR transcriptional complex. Altogether, these data indicate that AS mice are more sensitive to acute glucocorticoid exposure in the brain compared to WT mice. This suggests that stress responsiveness is altered in AS which could lead to anxiety symptoms.

The multivalency of the glucocorticoid receptor ligand-binding domain explains its manifold physiological activities.

The glucocorticoid receptor (GR) is a ubiquitously expressed transcription factor that controls metabolic and homeostatic processes essential for life. Although numerous crystal structures of the GR ligand-binding domain (GR-LBD) have been reported, the functional oligomeric state of the full-length receptor, which is essential for its transcriptional activity, remains disputed. Here we present five new crystal structures of agonist-bound GR-LBD, along with a thorough analysis of previous structural work. We identify four distinct homodimerization interfaces on the GR-LBD surface, which can associate into 20 topologically different homodimers. Biologically relevant homodimers were identified by studying a battery of GR point mutants including crosslinking assays in solution, quantitative fluorescence microscopy in living cells, and transcriptomic analyses. Our results highlight the relevance of non-canonical dimerization modes for GR, especially of contacts made by loop L1-3 residues such as Tyr545. Our work illustrates the unique flexibility of GR's LBD and suggests different dimeric conformations within cells. In addition, we unveil pathophysiologically relevant quaternary assemblies of the receptor with important implications for glucocorticoid action and drug design.

A plant-derived glucocorticoid receptor modulator with potency to attenuate the side effects of glucocorticoid therapy.

Continuous efforts have been made to move towards maintaining the beneficial anti-inflammatory functions of glucocorticoids (GCs) while minimizing side effects. Here, we investigated the selective glucocorticoid receptor (GR) modulator-like properties of a plant-derived compound caesaldekarin e (CA-e). The therapeutic efficacy of CA-e was evaluated in several mouse models, including dextran sulfate sodium-induced colitis, ovalbumin-induced lung allergic inflammation, imiquimod-induced psoriasis-like skin inflammation and skin atrophy. The action of CA-e targeting the GR was analysed using molecular docking, cellular thermal shift assays and microscale thermophoresis. Other methods included DNA-protein pull-down assays and mass spectrometry. CA-e selectively inhibited positive GC response element ((+) GRE)-mediated direct transactivation while maintaining and even enhancing the anti-inflammatory effects of treatment with dexamethasone. CA-e, alone and in combination with dexamethasone, efficiently alleviated inflammation in several mouse models with milder side effects compared with dexamethasone alone. Mechanistically, CA-e inhibited the formation of dimers by binding to the dimerization interface located in the ligand-binding domain of GR and facilitated embryonic ectoderm development that is involved in the regulation of transcriptional repression to compete for binding to (+) GRE, eventually leading to the repression of (+) GRE-regulated genes. In addition, CA-e repressed NF-κB-dependent genes by enhancing the interaction between GR and p65. Our results reveal that CA-e is a novel GR modulator with strong potency to attenuate the side effects of GC therapy and can be used as a potential molecular tool for deciphering GR signalling.

Cortisol resistance in the degu (Octodon degus)

Cortisol resistance has also been reported in the degu, Octodon degus, a New World hystricomorph endemic to central Chile. The degu is used as a model for studies of stress and diurnal rhythms, parental behaviour and female masculinization. Another New World hystricomorph, the guinea pig, also exhibits glucocorticoid resistance, a result of amino acid sequences that differ from other mammalian glucocorticoid receptors (GR). Mutations in the ligand-binding domain of the human GR have been identified in familial or sporadic generalised cortisol resistance as have variants in the guinea pig. To address the possibility that the high levels of cortisol observed in the degu are a result of the same or similar sequence variations observed in the guinea pig GR, we have cloned, expressed and characterised the ligand-binding domain (LBD) of the degu GR. Somewhat unexpectedly, neither the amino acids nor the region involved in the resistance observed in the guinea pig GR are relevant in the degu GR. The relative resistance to cortisol observed in the degu GR is conferred by the substitution of two isoleucine residues, which are highly conserved in the GR across species, with a valine doublet. These amino acids lie in the region between helices 5 and 6 of the GR LBD, a region known to be important in determining the affinity of ligand-binding in steroid receptors.

Active unfolding of the glucocorticoid receptor by the Hsp70/Hsp40 chaperone system in single-molecule mechanical experiments

The glucocorticoid receptor (GR) is an important transcription factor and drug target linked to a variety of biological functions and diseases. It is one of the most stringent physiological clients of the Hsp90/Hsp70/Hsp40 chaperone system. In this study, we used single-molecule force spectroscopy by optical tweezers to observe the interaction of the GR’s ligand-binding domain (GR-LBD) with the Hsp70/Hsp40 chaperone system (Hsp70/40). We show in real time that Hsp70/40 can unfold the complete GR-LBD in a stepwise manner. Each unfolding step involves binding of an Hsp70 to the GR-LBD and subsequent adenosine triphosphate (ATP) hydrolysis, stimulated by Hsp40. The kinetics of chaperone-mediated unfolding depend on chaperone concentrations as well as the presence of the nucleotide exchange factor BAG1. We find that Hsp70/40 can stabilize new unfolding intermediates, showing that Hsp70/40 can directly interact with the folded core of the protein when working as an unfoldase. Our results support an unfolding mechanism where Hsp70 can directly bind to folded protein structures and unfold them upon ATP hydrolysis. These results provide important insights into the regulation of GR by Hsp70/40.

The switch from client holding to folding in the Hsp70/Hsp90 chaperone machineries is regulated by a direct interplay between co-chaperones

Folding of stringent clients requires transfer from Hsp70 to Hsp90. The co-chaperone Hop physically connects the chaperone machineries. Here, we define its role from the remodeling of Hsp70/40-client complexes to the mechanism of client transfer and the conformational switching from stalled to active client-processing states of Hsp90. We show that Hsp70 together with Hsp40 completely unfold a stringent client, the glucocorticoid receptor ligand-binding domain (GR-LBD) in large assemblies. Hop remodels these for efficient transfer onto Hsp90. As p23 enters, Hsp70 leaves the complex via switching between binding sites in Hop. Current concepts assume that to proceed to client folding, Hop dissociates and the co-chaperone p23 stabilizes the Hsp90 closed state. In contrast, we show that p23 functionally interacts with Hop, relieves the stalling Hsp90-Hop interaction, and closes Hsp90. This reaction allows folding of the client and is thus the key regulatory step for the progression of the chaperone cycle.

Structure of Hsp90-p23-GR reveals the Hsp90 client-remodelling mechanism.

Hsp90 is a conserved and essential molecular chaperone responsible for the folding and activation of hundreds of 'client' proteins1-3. The glucocorticoid receptor (GR) is a model client that constantly depends on Hsp90 for activity4-9. GR ligand binding was previously shown to nr inhibited by Hsp70 and restored by Hsp90, aided by the co-chaperone p2310. However, a molecular understanding of the chaperone-mediated remodelling that occurs between the inactive Hsp70-Hsp90 'client-loading complex' and an activated Hsp90-p23 'client-maturation complex' is lacking for any client, including GR. Here we present a cryo-electron microscopy (cryo-EM) structure of the human GR-maturation complex (GR-Hsp90-p23), revealing that the GR ligand-binding domain is restored to a folded, ligand-bound conformation, while being simultaneously threaded through the Hsp90 lumen. In addition, p23 directly stabilizes native GR using a C-terminal helix, resulting in enhanced ligand binding. This structure of a client bound to Hsp90 in a native conformation contrasts sharply with the unfolded kinase-Hsp90 structure11. Thus, aided by direct co-chaperone-client interactions, Hsp90 can directly dictate client-specific folding outcomes. Together with the GR-loading complex structure12, we present the molecular mechanism of chaperone-mediated GR remodelling, establishing the first, to our knowledge,complete chaperone cycle for anyHsp90 client.

Discovery of Novel GR Ligands toward Druggable GR Antagonist Conformations Identified by MD Simulations and Markov State Model Analysis.

Binding of different ligands to glucocorticoid receptor (GR) may induce different conformational changes and even trigger completely opposite biological functions. To understand the allosteric communication within the GR ligand binding domain, the folding pathway of helix 12 (H12) induced by the binding of the agonist dexamethasone (DEX), antagonist RU486, and modulator AZD9567 are explored by molecular dynamics simulations and Markov state model analysis. The ligands can regulate the volume of the activation function‐2 through the residues Phe737 and Gln738. Without ligand or with agonist binding, H12 swings from inward to outward to visit different folding positions. However, the binding of RU486 or AZD9567 perturbs the structural state, and the passive antagonist state appears more stable. Structure‐based virtual screening and in vitro bioassays are used to discover novel GR ligands that bias the conformation equilibria toward the passive antagonist state. HP‐19 exhibits the best anti‐inflammatory activity (IC50 = 0.041 ± 0.011 µm) in nuclear factor‐kappa B signaling pathway, which is comparable to that of DEX. HP‐19 also does not induce adverse effect‐related transactivation functions of GR. The novel ligands discovered here may serve as promising starting points for the development of GR modulators.

TAZ inhibits glucocorticoid receptor and coordinates hepatic glucose homeostasis in normal physiological states.

The elucidation of the mechanisms whereby the liver maintains glucose homeostasis is crucial for the understanding of physiological and pathological states. Here, we show a novel role of hepatic transcriptional co-activator with PDZ-binding motif (TAZ) in the inhibition of glucocorticoid receptor (GR). TAZ is abundantly expressed in pericentral hepatocytes and its expression is markedly reduced by fasting. TAZ interacts via its WW domain with the ligand-binding domain of GR to limit the binding of GR to the GR response element in gluconeogenic gene promoters. Therefore, liver-specific TAZ knockout mice show increases in glucose production and blood glucose concentration. Conversely, the overexpression of TAZ in mouse liver reduces the binding of GR to gluconeogenic gene promoters and glucose production. Thus, our findings demonstrate that hepatic TAZ inhibits GR transactivation of gluconeogenic genes and coordinates gluconeogenesis in response to physiological fasting and feeding.

Glucocorticoid receptor condensates link DNA-dependent receptor dimerization and transcriptional transactivation

The glucocorticoid receptor (GR) is a ligand-regulated transcription factor (TF) that controls the tissue- and gene-specific transactivation and transrepression of thousands of target genes. Distinct GR DNA-binding sequences with activating or repressive activities have been identified, but how they modulate transcription in opposite ways is not known. We show that GR forms phase-separated condensates that specifically concentrate known coregulators via their intrinsically disordered regions (IDRs) in vitro. A combination of dynamic, multivalent (between IDRs) and specific, stable interactions (between LxxLL motifs and the GR ligand-binding domain) control the degree of recruitment. Importantly, GR DNA binding directs the selective partitioning of coregulators within GR condensates such that activating DNAs cause enhanced recruitment of coactivators. Our work shows that condensation controls GR function by modulating coregulator recruitment and provides a mechanism for the up- and down-regulation of GR target genes controlled by distinct DNA recognition elements.

Unraveling the Allosteric Cross-Talk between the Coactivator Peptide and the Ligand-Binding Site in the Glucocorticoid Receptor.

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known allosteric cross-talk between the ligand and coregulator binding sites within the GR ligand-binding domain that is crucial for the control of the functional response. However, the molecular mechanisms underlying such an allosteric control remain elusive. Here, molecular dynamics (MD) simulations, bioinformatic analysis, and biophysical measurements are integrated to capture the structural and dynamic features of the allosteric cross-talk within the GR. We identified a network of evolutionarily conserved residues that enables the allosteric signal transduction, in agreement with experimental data. MD simulations clarify how such a network is dynamically interconnected and offer a mechanistic explanation of how different peptides affect the intensity of the allosteric signal. This study provides useful insights to elucidate the GR allosteric regulation, ultimately providing a foundation for designing novel drugs.

Glucocorticoid Receptor Condensates Link DNA-Dependent Receptor Dimerization and Transcriptional Transactivation

The glucocorticoid receptor (GR) is a ligand-regulated transcription factor (TF) that controls the tissue- and gene-specific transactivation and transrepression of thousands of target genes. Distinct GR DNA binding sequences with activating or repressive activities have been identified, but how they modulate transcription in opposite ways is not known. We show that GR forms phase-separated condensates that specifically concentrate known co-regulators via their intrinsically disordered regions (IDRs) in vitro. A combination of dynamic, multivalent (between IDRs) and specific, stable interactions (between LxxLL motifs and the GR ligand binding domain) control the degree of recruitment. Importantly, GR DNA-binding directs the selective partitioning of co-regulators within GR condensates such that activating DNAs cause enhanced recruitment of co-activators. Our work shows that condensation controls GR function by modulating co-regulator recruitment and provides a mechanism for the up- and down-regulation of GR target genes controlled by distinct DNA recognition elements.

Unravelling Allosteric Cross-Talk between Co-Activator Peptide and Ligand Binding Site in Glucocorticoid Receptor

Glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known interconnection between the ligand and co-regulator binding sites of the GR ligand binding domain that is crucial for the functional differentiation of the downstream gene transcription. However, the molecular mechanisms underlying such an allosteric control remain elusive. The integration of molecular dynamics simulations, bioinformatic analysis and biophysical measurements led us to detect a network of dynamically interconnected evolutionary conserved residues that mediates the allosteric signal propagation between the co-regulator docking site and the ligand binding pocket. Co-regulator peptides alter the functionality of this network, explaining why the allosteric response depends on the peptide sequence. These findings provide useful insights for the comprehension of GR allosteric regulation that should help the design of novel drugs.

Dynamic allosteric communication pathway directing differential activation of the glucocorticoid receptor.

Allosteric communication within proteins is a hallmark of biochemical signaling, but the dynamic transmission pathways remain poorly characterized. We combined NMR spectroscopy and surface plasmon resonance to reveal these pathways and quantify their energetics in the glucocorticoid receptor, a transcriptional regulator controlling development, metabolism, and immune response. Our results delineate a dynamic communication network of residues linking the ligand-binding pocket to the activation function-2 interface, where helix 12, a switch for transcriptional activation, exhibits ligand- and coregulator-dependent dynamics coupled to graded activation. The allosteric free energy responds to variations in ligand structure: subtle changes gradually tune allostery while preserving the transmission pathway, whereas substitution of the entire pharmacophore leads to divergent allosteric control by apparently rewiring the communication network. Our results provide key insights that should aid in the design of mechanistically differentiated ligands.

Structural Insights into the Ligand Binding Domain of the Glucocorticoid Receptor: A Molecular Dynamics Study.

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.

First High-Resolution Crystal Structures of the Glucocorticoid Receptor Ligand-Binding Domain-Peroxisome Proliferator-Activated γ Coactivator 1-α Complex with Endogenous and Synthetic Glucocorticoids.

Both synthetic and endogenous glucocorticoids are important pharmaceutic drugs known to bind to the ligand-binding domain (LBD) of glucocorticoid receptor (GR), a member of the nuclear receptor (NR) superfamily. Ligand binding induces conformational changes within GR, resulting in subsequent DNA binding and differential coregulator recruitment, ultimately activating or repressing target gene expression. One of the most crucial coregulators is peroxisome proliferator-activated γ coactivator 1-α (PGC1α), which acts to regulate energy metabolism by directly interacting with GR to modulate gene expression. However, the mechanisms through which PGC1α senses GR conformation to drive transcription are not completely known. Here, an ancestral variant of the GR (AncGR2) LBD was used as a tool to produce stable protein for biochemical and structural studies. PGC1α is found to interact more tightly and form a more stable complex with AncGR2 LBD than nuclear receptor coactivator 2. We report the first high-resolution X-ray crystal structures of AncGR2 LBD in complex with PGC1α and dexamethasone (DEX) or hydrocortisone (HCY). Structural analyses reveal how distinct steroid drugs bind to GR with different affinities by unique hydrogen bonds and hydrophobic interactions. Important charge clamps are formed between the activation function 2 and PGC1α to mediate their specific interactions. These interactions lead to a high level of protection from hydrogen-deuterium exchange at the coregulator interaction site and strong intramolecular allosteric communication to ligand binding site. This is the first structure detailing the GR-PGC1α interaction providing a foundation for future design of specific therapeutic agents targeting these critical metabolic regulators. SIGNIFICANCE STATEMENT: High-resolution structures of AncGR2 LBD bound to DEX and HCY in complex with PGC1α reveal the molecular mechanism of PGC1α binding to AncGR2 LBD as well as the distinct affinities between DEX and HCY binding. Identifying the structural mechanisms that drive drug affinity is of pharmacologic interest to the glucocorticoid receptor field as an avenue to guide future drug design targeting GR-PGC1α signaling, which plays a crucial role in controlling hepatic glucose output.

Deciphering the structural basis for glucocorticoid resistance caused by missense mutations in the ligand binding domain of glucocorticoid receptor

The glucocorticoid resistance hereditary condition may emerge from the occurrence of point mutations in the glucocorticoid receptor (GR), which could impair its functionality. Because the main feature of such pathology is the resistance of the hypothalamic-pituitary-adrenal axis to the hormone cortisol, we used the GR ligand binding domain three-dimensional structure to perform computational analysis for eight variants known to cause this clinical condition (I559 N, V571A, D641V, G679S, F737L, I747 M, L753F and L773P), aiming to understand, on the atom scale, how they cause glucocorticoid resistance. We observed that the mutations generated a reduced affinity to cortisol and they alter some loop conformations, which could be a consequence from changes in protein motion, which in turn could result from the reduced stability of mutant GR structures. Therefore, the analyzed mutations compromise the GR ligand binding domain structure and cortisol binding, which could characterize the glucocorticoid resistance phenotype.

Structural and Dynamics of Glucocorticoid Receptor Ligand Binding Domain in Complex with PGC1α and Synthetic Glucocorticoids

Synthetic glucocorticoids are important pharmaceutical drugs that are known to bind to the ligand binding domain (LBD) of glucocorticoid receptor (GR), a member of the nuclear receptor superfamily. These binding will induce conformational changes within GR for the subsequent DNA binding and different co‐regulators recruitments, ultimately lead to activate or repress target gene expression. One of the crucial co‐regulators is peroxisome proliferator‐activated gamma coactivator 1‐α (PGC1α), known to regulates energy metabolism by direct interacting with GR to modulate gene expression. However, the mechanisms through which PGC1α changes GR conformation to drive transcription are unknown. Here, an ancestral variant of GR (GR2) was utilized as a tool to produce stable protein for biochemical and structural studies. We report the first high resolution X‐ray crystal structures of GR2 LBD in complexes with PGC1α and dexamethasone and hydrocortisone. These structural analyses reveal how distinct steroid drugs bind to GR with different affinities via unique hydrogen bonding and hydrophobic interactions. Hydrogen‐deuterium exchange mass spectrometry showed different effects on GR2 LBD allosteric communication after glucocorticoids binding. Molecular dynamics simulations and network analysis identified the key residues involved in the intramolecular communication between ligand binding site and activation function surface, which is the PGC1α binding site.Support or Funding InformationThis work was supported by W.M. Keck Foundation Medical Research Grant to E. A. O. X. L. was supported by an American Heart Association postdoctoral fellowship [17POST33660110].This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Discovery of New non-steroidal selective glucocorticoid receptor agonists

Glucocorticoids (GCs) are widely used as potent anti-inflammatory drugs; however, GC therapy is often accompanied by adverse side effects. The anti-inflammatory action of GCs is exerted through the glucocorticoid receptor (GR) in part by antagonizing the pro-inflammatory nuclear factor k B (NF-kB) whereas the majority of side effects are assumed to be mediated by transactivation of GR target genes. We set out to identify novel non-steroidal selective GR agonists (SEGRA) favoring transrepression of NF-kB target genes over transactivation of genes associated with undesirable effects. Our virtual screening protocol was driven by a pharmacophore model based on a pyrrolidinone amide analogue (named as ‘compound 12′ in Biggadike et al 2009, PNAS USA 106, 18,114) bound to the extended binding pocket of the GR ligand binding domain (GR-LBD). Ambinter library (7.8 million compounds) was queried by our validated pharmacophore hypothesis and the prioritized compounds were biologically evaluated using a series of well-established screening assays. Two structurally similar hits (1 and 13) were identified that bind to GR, induce its translocation to the nucleus, do not mediate transactivation of GR target genes whereas partially repress a number of pro-inflammatory NF-kB target genes, in a GR-dependent manner. Explanatory molecular dynamics (MD) calculations could detail the per-residue interactions accounting for the binding of 1 and 13 to the extended binding pocket of GR. The discovered 1,3-benzothiazole analogs introduce a new class of genuine SEGRA paving the way for hit-to-lead optimization.