RORγt, an isoform of the retinoic acid-related orphan receptor gamma (RORc, RORγ), has been identified as the master regulator of T-helper 17 (TH 17) cell function and development, making it an attractive target for the treatment of autoimmune diseases. Validation for this target comes from antibodies targeting interleukin-17 (IL-17), the signature cytokine produced by TH 17 cells, which have shown impressive results in clinical trials. Through focused screening of our compound collection, we identified a series of N-sulfonylated benzoxazepines, which displayed micromolar affinity for the RORγ ligand-binding domain (LBD) in a radioligand binding assay. Optimization of these initial hits resulted in potent binders, which dose-dependently decreased the ability of the RORγ-LBD to interact with a peptide derived from steroid receptor coactivator 1, and inhibited the release of IL-17 secretion from isolated and cultured human TH 17 cells with nanomolar potency. A cocrystal structure of inverse agonist 15 (2-chloro-6-fluoro-N-(4-{[3-(trifluoromethyl)phenyl]sulfonyl}-2,3,4,5-tetrahydro-1,4-benzoxazepin-7-yl)benzamide) bound to the RORγ-LBD illustrated that both hydrophobic interactions, leading to an induced fit around the substituted benzamide moiety of 15, as well as a hydrogen bond from the amide NH to His479 seemed to be important for the mechanism of action. This structure is compared with the structure of agonist 25 (N-(2-fluorophenyl)-4-[(4-fluorophenyl)sulfonyl]-2,3,4,5-tetrahydro-1,4-benzoxazepin-6-amine ) and structures of other known RORγ modulators.

In this study, we performed an extensive exploration of the ligand entry mechanism for members of the steroid nuclear hormone receptor family (androgen receptor, estrogen receptor α, glucocorticoid receptor, mineralocorticoid receptor, and progesterone receptor) and their endogenous ligands. The exploration revealed a shared entry path through the helix 3, 7, and 11 regions. Examination of the x-ray structures of the receptor-ligand complexes further showed two distinct folds of the helix 6-7 region, classified as "open" and "closed", which could potentially affect ligand binding. To improve sampling of the helix 6-7 loop, we incorporated motion modes based on principal component analysis of existing crystal structures of the receptors and applied them to the protein-ligand sampling. A detailed comparison with the anisotropic network model (an elastic network model) highlights the importance of flexibility in the entrance region. While the binding (interaction) energy of individual simulations can be used to score different ligands, extensive sampling further allows us to predict absolute binding free energies and analyze reaction kinetics using Markov state models and Perron-cluster cluster analysis, respectively. The predicted relative binding free energies for three ligands binding to the progesterone receptor are in very good agreement with experimental results and the Perron-cluster cluster analysis highlighted the importance of a peripheral binding site. Our analysis revealed that the flexibility of the helix 3, 7, and 11 regions represents the most important factor for ligand binding. Furthermore, the hydrophobicity of the ligand influences the transition between the peripheral and the active binding site.

A prerequisite for successful drugs is effective binding of the desired target protein in the complex environment of a living system. Drug-target engagement has typically been difficult to monitor in physiologically relevant models, and with current methods, especially, while maintaining spatial information. One recent technique for quantifying drug-target engagement is the cellular thermal shift assay (CETSA), in which ligand-induced protein stabilization is measured after a heat challenge. Here, we describe a CETSA protocol in live A431 cells for p38α (MAPK14), where remaining soluble protein is detected in situ, using high-content imaging in 384-well, microtiter plates. We validate this assay concept using a number of known p38α inhibitors and further demonstrate the potential of this technology for chemical probe and drug discovery purposes by performing a small pilot screen for novel p38α binders. Importantly, this protocol creates a workflow that is amenable to adherent cells in their native state and yields spatially resolved target engagement information measurable at the single-cell level.

No abstract

The bromodomain-containing proteins are a ligandable family of epigenetic readers, which play important roles in oncological, cardiovascular, and inflammatory diseases. Achieving selective inhibition of specific bromodomains is challenging, due to the limited understanding of compound and target selectivity features. In this study we build and benchmark proteochemometric (PCM) classification models on bioactivity data for 15,350 data points across 31 bromodomains, using both compound fingerprints and binding site protein descriptors as input variables, achieving a maximum performance as measured by the Matthew's Correlation Coefficient (MCC) of 0.83 on the external test set. We also find that histone peptide binding data can be used as a target descriptor to build a high performing PCM model (MCC 0.80), showing the transferability of peptide interaction information to modeling small-molecule bioactivity. 1,139 compounds were selected for prospective experimental testing by performing a virtual screen using model predictions and implementing conformal prediction, which resulted in 319 correctly predicted compound-target pair actives and the correct prediction for certain selectivity profile combinations of the four bromodomains tested against. We identify that conformal prediction can be used to fine-tune the balance between hit retrieval and hit structural diversity in a virtual screening setting. PCM can be applied to future virtual screening and compound design, including off-target prediction for bromodomains.