Molecular Mechanism of Action of RORγt Agonists and Inverse Agonists: Insights from Molecular Dynamics Simulation

Development and therapeutic potential of allosteric retinoic acid receptor-related orphan receptor γt (RORγt) inverse agonists for autoimmune diseases

The transcription factor retinoic acid receptor-related orphan receptor γt (RORγt) is an attractive drug target for some autoimmune diseases owing to its roles in the differentiation of human T helper 17 (Th17) cells which produce pro-inflammatory cytokine interleukin (IL)-17. RORγt agonists and inverse agonists are classically targeted to the hydrophobic and highly conserved orthosteric binding pocket of RORγt ligand binding domain (LBD). Although successful, this approach also brings some challenges, including off-target effects due to lack of selectivity over other nuclear receptors (NRs). Allosteric regulation of RORγt by synthetic small molecules has recently emerged as novel research interests for its interesting modes of action (MOA), satisfying bioactivity profile and improved selectivity. In this review, we delineated the discovery and identification of the allosteric pocket of RORγt. Subsequently, we focused on examples of small molecules that allosterically inhibit RORγt, with a central attention on structural-activity-relationship (SAR) information, biological activity, pharmacokinetic (PK) property, and the ligand binding mode of these compounds. We also discussed the potential role of RORγt allosteric inverse agonists as small molecule therapeutics for autoimmune diseases.

Discovery, Synthesis, and In Vitro Characterization of 2,3 Derivatives of 4,5,6,7-Tetrahydro-Benzothiophene as Potent Modulators of Retinoic Acid Receptor-Related Orphan Receptor γt

Retinoic acid receptor-related orphan receptor γt (RORγt) is a nuclear receptor that is expressed in a variety of tissues and is a potential drug target for the treatment of inflammatory and auto-immune diseases, metabolic diseases, and resistant cancer types. We herein report the discovery of 2,3 derivatives of 4,5,6,7-tetrahydro-benzothiophene modulators of RORγt. We also report the solubility in acidic/neutral pH, mouse/human/dog/rat microsomal stability, Caco-2, and MDR1-MDCKII permeabilities of a set of these derivatives. For this group of modulators, inverse agonism by steric clashes and push-pull mechanisms induce greater instability to protein conformation compared to agonist lock hydration. Independent of the two mechanisms, we observed a basal modulatory activity of the tested 2,3 derivatives of 4,5,6,7-tetrahydro-benzothiophene toward RORγt due to the interactions with the Cys320-Glu326 and Arg364-Phe377 hydrophilic regions. The drug discovery approach reported in the current study can be employed to discover modulators of nuclear receptors and other globular protein targets.

Virtual Screening Strategy to Identify Retinoic Acid-Related Orphan Receptor γt Modulators

Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target's binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%.

Discovery of a Novel Class of Dual GPBAR1 Agonists-RORγt Inverse Agonists for the Treatment of IL-17-Mediated Disorders

Retinoic acid receptor-related orphan receptor γ-t (RORγt) and GPBAR1, a transmembrane G-protein-coupled receptor for bile acids, are attractive drug targets to develop clinically relevant small modulators as potent therapeutics for autoimmune diseases. Herein, we designed, synthesized, and evaluated several new bile acid-derived ligands with potent dual activity. Furthermore, we performed molecular docking and MD calculations of the best dual modulators in the two targets to identify the binding modes as well as to better understand the molecular basis of the inverse agonism of RORγt by bile acid derivatives. Among these compounds, 7 was identified as a GPBAR1 agonist (EC50 5.9 μM) and RORγt inverse agonist (IC50 0.107 μM), with excellent pharmacokinetic properties. Finally, the most promising ligand displayed robust anti-inflammatory activity in vitro and in vivo in a mouse model of 2,4,6-trinitrobenzene sulfonic acid (TNBS)-induced colitis.

Modes of action insights from the crystallographic structures of retinoic acid receptor-related orphan receptor-γt (RORγt)

RORγt plays an important role in mediating IL-17 production and some tumor cells. It has four functional domains, of which the ligand-binding domain (LBD) is responsible for binding agonists to recruit co-activators or inverse agonists to prevent co-activator recruiting the agonists. Thus, potent ligands targeting the LBD of this protein could provide novel treatments for cancer and autoimmune diseases. In this perspective, we summarized and discussed various modes of action (MOA) of RORγt-ligand binding structures. The ligands can bind with RORγt at either orthosteric site or the allosteric site, and the binding modes at these two sites are different for agonists and inverse agonist. At the orthosteric site, the binding of agonist is to stabilize the H479-Y502-F506 triplet interaction network of RORγt. The binding of inverse agonist features as these four apparent ways: (1) blocking the entrance of the agonist pocket in RORγt; (2) directly breaking the H479-Y502 pair interactions; (3) destabilizing the triplet H479-Y502-F506 interaction network through perturbing the conformation of the side chain in M358 at the bottom of the binding pocket; (4) and destabilizing the triplet H479-Y502-F506 through changing the conformation of the side chain of residue W317 side chain. At the allosteric site of RORγt, the binding of inverse agonist was found recently to inhibit the activation of protein by interacting directly with H12, which results in unfolding of helix 11' and orientation of H12 to directly block cofactor peptide binding. This overview of recent advances in the RORγt structures is expected to provide a guidance of designing more potent drugs to treat RORγt-related diseases.

Crystallography-guided discovery of carbazole-based retinoic acid-related orphan receptor gamma-t (RORγt) modulators: insights into different protein behaviors with "short" and "long" inverse agonists

A series of 6-substituted carbazole-based retinoic acid-related orphan receptor gamma-t (RORγt) modulators were discovered through 6-position modification guided by insights from the crystallographic profiles of the "short" inverse agonist 6. With the increase in the size of the 6-position substituents, the "short" inverse agonist 6 first reversed its function to agonists and then to "long" inverse agonists. The cocrystal structures of RORγt complexed with the representative "short" inverse agonist 6 (PDB: 6LOB), the agonist 7d (PDB: 6LOA) and the "long" inverse agonist 7h (PDB: 6LO9) were revealed by X-ray analysis. However, minor differences were found in the binding modes of "short" inverse agonist 6 and "long" inverse agonist 7h. To further reveal the molecular mechanisms of different RORγt inverse agonists, we performed molecular dynamics simulations and found that "short" or "long" inverse agonists led to different behaviors of helixes H11, H11', and H12 of RORγt. The "short" inverse agonist 6 destabilizes H11' and dislocates H12, while the "long" inverse agonist 7h separates H11 and unwinds H12. The results indicate that the two types of inverse agonists may behave differently in downstream signaling, which may help identify novel inverse agonists with different regulatory mechanisms.

Molecular Dynamics Simulations in Designing DARPins as Phosphorylation-Specific Protein Binders of ERK2

Extracellular signal-regulated kinases 1 and 2 (ERK1/2) play key roles in promoting cell survival and proliferation through the phosphorylation of various substrates. Remarkable antitumour activity is found in many inhibitors that act upstream of the ERK pathway. However, drug-resistant tumour cells invariably emerge after their use due to the reactivation of ERK1/2 signalling. ERK1/2 inhibitors have shown clinical efficacy as a therapeutic strategy for the treatment of tumours with mitogen-activated protein kinase (MAPK) upstream target mutations. These inhibitors may be used as a possible strategy to overcome acquired resistance to MAPK inhibitors. Here, we report a class of repeat proteins-designed ankyrin repeat protein (DARPin) macromolecules targeting ERK2 as inhibitors. The structural basis of ERK2-DARPin interactions based on molecular dynamics (MD) simulations was studied. The information was then used to predict stabilizing mutations employing a web-based algorithm, MAESTRO. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations. Two mutations, Ala → Asp and Ser → Leu, were found to perform better than the original sequence (DARPin E40) based on the associated energy and key residues involved in protein-protein interaction. MD simulations and analysis of the data obtained on these mutations supported our predictions.

Discovery of tetrahydroquinolines and benzomorpholines as novel potent RORγt agonists

The retinoic acid receptor-related orphan receptor γt (RORγt) is an important nuclear receptor that regulates the differentiation of Th17 cells and production of interleukin 17(IL-17). RORγt agonists increase basal activity of RORγt and could provide a potential approach to cancer immunotherapy. Herein, hit compound 1 was identified as a weak RORγt agonist during in-house library screening. Changes in LHS core of 1 led to the identification of tetrahydroquinoline compound 6 as a partial RORγt agonist (max. act. = 39.3%). Detailed structure-activity relationship on substituent of the LHS core, amide linker and RHS arylsulfonyl moiety was explored and a novel series of tetrahydroquinolines and benzomorpholines was discovered as potent RORγt agonists. Tetrahydroquinoline compound 8g (EC₅₀ = 8.9 ± 0.4 nM, max. act. = 104.5%) and benzomorpholine compound 9g (EC₅₀ = 7.5 ± 0.6 nM, max. act. = 105.8%) were representative compounds with high RORγt agonistic activity in dual FRET assay, and they showed good activity in cell-based Gal4 reporter gene assay and Th17 cell differentiation assay (104.5% activation at 300 nM of 8g; 59.4% activation at 300 nM of 9g). The binding modes of 8g and 9g as well as the two RORγt inverse agonists accidentally discovered were also discussed.

RORγ Structural Plasticity and Druggability

Retinoic acid receptor-related orphan receptor γ (RORγ) is a transcription factor regulating the expression of the pro-inflammatory cytokine IL-17 in human T helper 17 (Th17) cells. Activating RORγ can induce multiple IL-17-mediated autoimmune diseases but may also be useful for anticancer therapy. Its deep immunological functions make RORɣ an attractive drug target. Over 100 crystal structures have been published describing atomic interactions between RORɣ and agonists and inverse agonists. In this review, we focus on the role of dynamic properties and plasticity of the RORɣ orthosteric and allosteric binding sites by examining structural information from crystal structures and simulated models. We discuss the possible influences of allosteric ligands on the orthosteric binding site. We find that high structural plasticity favors the druggability of RORɣ, especially for allosteric ligands.

Delineation of the molecular determinants of the unique allosteric binding site of the orphan nuclear receptor RORγt

Nuclear receptors (NRs) are high-interest targets in drug discovery because of their involvement in numerous biological processes and diseases. Classically, NRs are targeted via their hydrophobic, orthosteric pocket. Although successful, this approach comes with challenges, including off-target effects due to lack of selectivity. Allosteric modulation of NR activity constitutes a promising pharmacological strategy. The retinoic acid receptor-related orphan receptor-γt (RORγt) is a constitutively active NR that positively regulates the expression of interleukin-17 in T helper 17 cells. Inhibiting this process is an emerging strategy for managing autoimmune diseases. Recently, an allosteric binding pocket in the C-terminal region of the ligand-binding domain (LBD) of RORγt was discovered that is amenable to small-molecule drug discovery. Compounds binding this pocket induce a reorientation of helix 12, thereby preventing coactivator recruitment. Therefore, inverse agonists binding this site with high affinity are actively being pursued. To elucidate the pocket formation mechanism, verify the uniqueness of this pocket, and substantiate the relevance of targeting this site, here we identified the key characteristics of the RORγt allosteric region. We evaluated the effects of substitutions in the LBD on coactivator, orthosteric, and allosteric ligand binding. We found that two molecular elements unique to RORγt, the length of helix 11' and a Gln-487 residue, are crucial for the formation of the allosteric pocket. The unique combination of elements present in RORγt suggests a high potential for subtype-selective targeting of this NR to more effectively treat patients with autoimmune diseases.

Structural change of retinoic-acid receptor-related orphan receptor induced by binding of inverse-agonist: Molecular dynamics and ab initio molecular orbital simulations

To elucidate structural changes in the retinoic acid receptor-related orphan receptor gamma (RORγt) induced by the binding of an agonist or an inverse agonist, we conducted molecular dynamics (MD) simulations in explicit water. In addition, ab initio fragment molecular orbital calculations were carried out for certain characteristic structures obtained from the MD simulations to reveal important interactions between the amino acid residues of RORγt, and to distinguish the different effects in the binding of an agonist and an inverse agonist on the structure of RORγt. The results elucidate that the hydrogen bond between His479 of helix11 (H11) and Tyr502 of helix12 (H12) is important to keep the H12 conformation in the agonist-bound RORγt. In contrast, in the inverse-agonist-bound RORγt, the side chain of His479 rotates, significantly weakening the interaction between His479 and Tyr502, leading to a conformational change in H12. Therefore, the present molecular simulations clearly indicate that the conformational change in the side chain of His479 in the inverse-agonist-bound RORγt is the main reason for the H12 destabilization induced by the binding of the inverse agonist. Such a conformational change does not occur on the binding of the agonist in RORγt, owing to the strong hydrogen bond between His479 and Tyr502.

HDX-MS reveals structural determinants for RORγ hyperactivation by synthetic agonists

Members of the nuclear receptor (NR) superfamily regulate both physiological and pathophysiological processes ranging from development and metabolism to inflammation and cancer. Synthetic small molecules targeting NRs are often deployed as therapeutics to correct aberrant NR signaling or as chemical probes to explore the role of the receptor in physiology. Nearly half of NRs do not have specific cognate ligands (termed orphan NRs) and it’s unclear if they possess ligand dependent activities. Here we demonstrate that ligand-dependent action of the orphan RORγ can be defined by selectively disrupting putative endogenous—but not synthetic—ligand binding. Furthermore, the characterization of a library of RORγ modulators reveals that structural dynamics of the receptor assessed by HDX-MS correlate with activity in biochemical and cell-based assays. These findings, corroborated with X-ray co-crystallography and site-directed mutagenesis, collectively reveal the structural determinants of RORγ activation, which is critical for designing RORγ agonists for cancer immunotherapy.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.