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

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

Glucocorticoid receptor signaling: intricacies and therapeutic opportunities

The glucocorticoid receptor (GR) is a major nuclear receptor (NR) drug target for the treatment of inflammatory disorders and several cancers. Despite the effectiveness of GR ligands, their systemic action triggers a plethora of side effects, limiting long-term use. Here, we discuss new concepts of and insights into GR mechanisms of action to assist in the identification of routes toward enhanced therapeutic benefits. We zoom in on the communication between different GR domains and how this is influenced by different ligands. We detail findings on the interaction between GR and chromatin, and highlight how condensate formation and coregulator confinement can perturb GR transcriptional responses. Last, we discuss the potential of novel ligands and the therapeutic exploitation of crosstalk with other NRs.

Influence of dexamethasone on the interaction between glucocorticoid receptor and SOX9: A molecular dynamics study

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. GR transcriptional activity is modulated by a series of ligands and coenzymes, where a ligand can act as an agonist or antagonist. GR agonists, such as the glucocorticoids dexamethasone (DEX) and prednisolone, are widely prescribed to patients with inflammatory and autoimmune diseases. DEX is also used to induce osteogenic differentiation in vitro. Recently, it has been highlighted that DEX induces changes in the osteogenic differentiation of human mesenchymal stromal cells by downregulating the transcription factor SRY-box transcription factor 9 (SOX9) and upregulating the peroxisome proliferator-activated receptor γ (PPARG). SOX9 is fundamental in the control of chondrogenesis, but also in osteogenesis by acting as a dominant-negative of RUNX2. Many processes remain to be clarified during cell fate determination, such as the interplay between the key transcription factors. The main objective pursued by this work is to shed light on the interaction between GR and SOX9 in the presence and absence of DEX at an atomic level of resolution using molecular dynamics simulations. The outcome of this research could help the understanding of possible molecular interactions between GR and SOX9 and their role in the determination of cell fate. The results highlight the key residues at the interface between GR and SOX9 involved in the complexation process and shed light on the mechanism through which DEX modulates GR-SOX9 binding and exerts its biological activity.

Combining structural and coevolution information to unveil allosteric sites

Understanding allosteric regulation in biomolecules is of great interest to pharmaceutical research and computational methods emerged during the last decades to characterize allosteric coupling. However, the prediction of allosteric sites in a protein structure remains a challenging task. Here, we integrate local binding site information, coevolutionary information, and information on dynamic allostery into a structure-based three-parameter model to identify potentially hidden allosteric sites in ensembles of protein structures with orthosteric ligands. When tested on five allosteric proteins (LFA-1, p38-α, GR, MAT2A, and BCKDK), the model successfully ranked all known allosteric pockets in the top three positions. Finally, we identified a novel druggable site in MAT2A confirmed by X-ray crystallography and SPR and a hitherto unknown druggable allosteric site in BCKDK validated by biochemical and X-ray crystallography analyses. Our model can be applied in drug discovery to identify allosteric pockets.

Deciphering the Shared and Specific Drug Resistance Mechanisms of Anaplastic Lymphoma Kinase via Binding Free Energy Computation

Anaplastic lymphoma kinase (ALK), a tyrosine receptor kinase, has been proven to be associated with the occurrence of numerous malignancies. Although there have been already at least 3 generations of ALK inhibitors approved by FDA or in clinical trials, the occurrence of various mutations seriously attenuates the effectiveness of the drugs. Unfortunately, most of the drug resistance mechanisms still remain obscure. Therefore, it is necessary to reveal the bottom reasons of the drug resistance mechanisms caused by the mutations. In this work, on the basis of verifying the accuracy of 2 main kinds of binding free energy calculation methodologies [end-point method of Molecular Mechanics with Poisson-Boltzmann/Generalized Born and Surface Area (MM/PB(GB)SA) and alchemical method of Thermodynamic Integration (TI)], we performed a systematic analysis on the ALK systems to explore the underlying shared and specific drug resistance mechanisms, covering the one-drug-multiple-mutation and multiple-drug-one-mutation cases. Through conventional molecular dynamics (cMD) simulation in conjunction with MM/PB(GB)SA and umbrella sampling (US) in conjunction with contact network analysis (CNA), the resistance mechanisms of the in-pocket, out-pocket, and multiple-site mutations were revealed. Especially for the out-pocket mutation, a possible transfer chain of the mutation effect was revealed, and the reason why different drugs exhibited various sensitivities to the same mutation was also uncovered. The proposed mechanisms may be prevalent in various drug resistance cases.

Deciphering the Structural Determinants Critical in Attaining the FXR Partial Agonism

Elucidation of structural determinants is pivotal for structure-based drug discovery. The Farnesoid X receptor (FXR) is a proven target for NASH; however, its full agonism causes certain clinical complications. Therefore, partial agonism (PA) appears as a viable alternative for improved therapeutics. Since the agonist and PA both share the same binding site, i.e., ligand-binding pocket (LBP), which is highly dynamic and has synergy with the substrate binding site, the selective designing of PA is challenging. The identification of structural and conformational determinants is critical for PA compared with an agonist. Furthermore, the mechanism by which PA modulates the structural dynamics of FXR at the residue level, a prerequisite for PA designing, is still elusive. Here, by using ∼4.5 μs of MD simulations and residue-wise communication network analysis, we identified the structural regions which are flexible with PA but frozen with an agonist. Also, the network analysis identified the considerable changes between an agonist and PA in biologically essential zones of FXR such as helix H10/H11 and loop L:H11/H12, which lead to the modulation of synergy between LBP and the substrate binding site. Furthermore, the thermodynamic profiling suggested the methionine residues, mainly M³²⁸, M³⁶⁵, and M⁴⁵⁰, seem to be responsible for the recruitment of PA. The other residues I³⁵⁷, Y³⁶¹, L⁴⁶⁵, F³⁰⁸, Q³¹⁶, and K³²¹ are also identified, exclusively interacting with PA. This study offers novel structural and mechanistic insights that are critical for FXR targeted drug discovery for PA designing.

Cooperation of structural motifs controls drug selectivity in cyclin-dependent kinases: an advanced theoretical analysis

Understanding drug selectivity mechanism is a long-standing issue for helping design drugs with high specificity. Designing drugs targeting cyclin-dependent kinases (CDKs) with high selectivity is challenging because of their highly conserved binding pockets. To reveal the underlying general selectivity mechanism, we carried out comprehensive analyses from both the thermodynamics and kinetics points of view on a representative CDK12 inhibitor. To fully capture the binding features of the drug-target recognition process, we proposed to use kinetic residue energy analysis (KREA) in conjunction with the community network analysis (CNA) to reveal the underlying cooperation effect between individual residues/protein motifs to the binding/dissociating process of the ligand. The general mechanism of drug selectivity in CDKs can be summarized as that the difference of structural cooperation between the ligand and the protein motifs leads to the difference of the energetic contribution of the key residues to the ligand. The proposed mechanisms may be prevalent in drug selectivity issues, and the insights may help design new strategies to overcome/attenuate the drug selectivity associated problems.

Ligand-induced shifts in conformational ensembles that describe transcriptional activation

Nuclear receptors function as ligand-regulated transcription factors whose ability to regulate diverse physiological processes is closely linked with conformational changes induced upon ligand binding. Understanding how conformational populations of nuclear receptors are shifted by various ligands could illuminate strategies for the design of synthetic modulators to regulate specific transcriptional programs. Here, we investigate ligand-induced conformational changes using a reconstructed, ancestral nuclear receptor. By making substitutions at a key position, we engineer receptor variants with altered ligand specificities. We combine cellular and biophysical experiments to characterize transcriptional activity, as well as elucidate mechanisms underlying altered transcription in receptor variants. We then use atomistic molecular dynamics (MD) simulations with enhanced sampling to generate ensembles of wildtype and engineered receptors in combination with multiple ligands, followed by conformational analysis and correlation of MD-based predictions with functional ligand profiles. We determine that conformational ensembles accurately describe ligand responses based on observed population shifts. These studies provide a platform which will allow structural characterization of physiologically-relevant conformational ensembles, as well as provide the ability to design and predict transcriptional responses in novel ligands.

Allosteric Antagonism of the Pregnane X Receptor (PXR): Current-State-of-the-Art and Prediction of Novel Allosteric Sites

The pregnane X receptor (PXR, NR1I2) is a xenobiotic-activated transcription factor with high levels of expression in the liver. It not only plays a key role in drug metabolism and elimination, but also promotes tumor growth, drug resistance, and metabolic diseases. It has been proposed as a therapeutic target for type II diabetes, metabolic syndrome, and inflammatory bowel disease, and PXR antagonists have recently been considered as a therapy for colon cancer. There are currently no PXR antagonists that can be used in a clinical setting. Nevertheless, due to the large and complex ligand-binding pocket (LBP) of the PXR, it is challenging to discover PXR antagonists at the orthosteric site. Alternative ligand binding sites of the PXR have also been proposed and are currently being studied. Recently, the AF-2 allosteric binding site of the PXR has been identified, with several compounds modulating the site discovered. Herein, we aimed to summarize our current knowledge of allosteric modulation of the PXR as well as our attempt to unlock novel allosteric sites. We describe the novel binding function 3 (BF-3) site of PXR, which is also common for other nuclear receptors. In addition, we also mention a novel allosteric site III based on in silico prediction. The identified allosteric sites of the PXR provide new insights into the development of safe and efficient allosteric modulators of the PXR receptor. We therefore propose that novel PXR allosteric sites might be promising targets for treating chronic metabolic diseases and some cancers.

Wandering beyond small molecules: peptides as allosteric protein modulators

Recent years have seen the rise of allosteric modulation as an innovative approach for drug design and discovery, efforts which culminated in the development of several clinical candidates. Allosteric modulation of many drug targets, including mainly membrane-embedded receptors, have been vastly explored through small molecule screening campaigns, but much less attention has been paid to peptide-based allosteric modulators. However, peptides have a significant impact on the pharmaceutical industry due to the typically higher potency and selectivity for their targets, as compared with small molecule therapeutics. Therefore, peptides represent one of the most promising classes of molecules that can modulate key biological pathways. Here, we report on the allosteric modulation of proteins (ranging from G protein-coupled receptors to specific protein-protein interactions) by peptides for applications in drug discovery.