Nuclear receptor interdomain communication is mediated by the hinge with ligand specificity

Nuclear receptors are ligand-induced transcription factors that bind directly to target genes and regulate their expression. Ligand binding initiates conformational changes that propagate to other domains, allosterically regulating their activity. The nature of this interdomain communication in nuclear receptors is poorly understood, largely owing to the difficulty of experimentally characterizing full-length structures. We have applied computational modeling approaches to describe and study the structure of the full length farnesoid X receptor (FXR), approximated by the DNA binding domain (DBD) and ligand binding domain (LBD) connected by the flexible hinge region. Using extended molecular dynamics simulations (> 10 microseconds) and enhanced sampling simulations, we provide evidence that ligands selectively induce domain rearrangement, leading to interdomain contact. We use protein-protein interaction assays to provide experimental evidence of these interactions, identifying a critical role of the hinge in mediating interdomain contact. Our results illuminate previously unknown aspects of interdomain communication in FXR and provide a framework to enable characterization of other full length nuclear receptors.

How nuclear receptors transition between active and inactive forms: An energetic perspective

Nuclear receptors regulate transcriptional programs in response to the binding of natural and synthetic ligands. These ligands modulate the receptor by inducing dynamic changes in the ligand binding domain that shift the C-terminal helix (H12) between active and inactive conformations. Despite decades of study, many questions persist regarding the nature of the inactive state and how ligands shift receptors between different states. Here, we use molecular dynamics (MD) simulations to investigate the timescale and energetic landscape of the conformational transition between inactive and active forms of progesterone receptor (PR) bound to a partial agonist. We observe that the microsecond timescale is insufficient to observe any transitions; only at millisecond timescales achieved via accelerated MD simulations do we find the inactive PR switches to the active state. Energetic analysis reveals that both active and inactive PR states represent energy minima separated by a barrier that can be traversed. In contrast, little or no transition is observed between active and inactive states when an agonist or antagonist is bound, confirming that ligand identity plays a key role in defining the energy landscape of nuclear receptor conformations.

Somatic estrogen receptor α mutations that induce dimerization promote receptor activity and breast cancer proliferation

Physiologic activation of estrogen receptor α (ERα) is mediated by estradiol (E2) binding in the ligand-binding pocket of the receptor, repositioning helix 12 (H12) to facilitate binding of coactivator proteins in the unoccupied coactivator binding groove. In breast cancer, activation of ERα is often observed through point mutations that lead to the same H12 repositioning in the absence of E2. Through expanded genetic sequencing of breast cancer patients, we identified a collection of mutations located far from H12 but nonetheless capable of promoting E2-independent transcription and breast cancer cell growth. Using machine learning and computational structure analyses, this set of mutants was inferred to act distinctly from the H12-repositioning mutants and instead was associated with conformational changes across the ERα dimer interface. Through both in vitro and in-cell assays of full-length ERα protein and isolated ligand-binding domain, we found that these mutants promoted ERα dimerization, stability, and nuclear localization. Point mutations that selectively disrupted dimerization abrogated E2-independent transcriptional activity of these dimer-promoting mutants. The results reveal a distinct mechanism for activation of ERα function through enforced receptor dimerization and suggest dimer disruption as a potential therapeutic strategy to treat ER-dependent cancers.

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.

A molecular dynamics perspective into estrogen receptor inhibition by selective flavonoids as alternative therapeutic options

Zearalenone is an estrogenic mycotoxin which is a common food contaminant and has been implicated in increasing the incidence of carcinogenesis and other reproductive health ailments through the estrogen receptor alpha (ERα) pathway. Competitive ERα blockers such as 4-Hydroxytamoxifen (OHT), are synthetic FDA approved drugs which, albeit being an effective anticancer agent, induces life altering side effects. For this reason, there is an increased interest in the use of naturally occurring medicinal plant products such as flavonoids. This study aimed to identity flavonoid ERα inhibitors and provide insights into the mechanism of inhibition using computational techniques. The Molecular Mechanics/Generalized Born Surface Area calculations revealed that quercetrin, hesperidin, epigallocatechin 3-gallate and kaempferol 7-O-glucoside out of 14 flavonoids had higher binding affinity for ERα than OHT. The structural analysis revealed that the binding of the compounds to the receptor lead to dynamic alterations, which induced conformational shift in the structure and orientation of the receptor resulting in stabilised, compact and low energy systems. The results of this study provide imperative information that supports the use of flavonoids in the inhibition of ERα to prevent or ameliorate the consequential adverse effects associated with zearalenone exposure.Communicated by Ramaswamy H. Sarma.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Elucidation of Agonist and Antagonist Dynamic Binding Patterns in ER-α by Integration of Molecular Docking, Molecular Dynamics Simulations and Quantum Mechanical Calculations

Estrogen receptor alpha (ERα) is a ligand-dependent transcriptional factor in the nuclear receptor superfamily. Many structures of ERα bound with agonists and antagonists have been determined. However, the dynamic binding patterns of agonists and antagonists in the binding site of ERα remains unclear. Therefore, we performed molecular docking, molecular dynamics (MD) simulations, and quantum mechanical calculations to elucidate agonist and antagonist dynamic binding patterns in ERα. 17β-estradiol (E2) and 4-hydroxytamoxifen (OHT) were docked in the ligand binding pockets of the agonist and antagonist bound ERα. The best complex conformations from molecular docking were subjected to 100 nanosecond MD simulations. Hierarchical clustering was conducted to group the structures in the trajectory from MD simulations. The representative structure from each cluster was selected to calculate the binding interaction energy value for elucidation of the dynamic binding patterns of agonists and antagonists in the binding site of ERα. The binding interaction energy analysis revealed that OHT binds ERα more tightly in the antagonist conformer, while E2 prefers the agonist conformer. The results may help identify ERα antagonists as drug candidates and facilitate risk assessment of chemicals through ER-mediated responses.

Bisphenol-S and Bisphenol-F alter mouse pancreatic β-cell ion channel expression and activity and insulin release through an estrogen receptor ERβ mediated pathway

Bisphenol-S (BPS) and Bisphenol-F (BPF) are current Bisphenol-A (BPA) substitutes. Here we used pancreatic β-cells from wild type (WT) and estrogen receptor β (ERβ) knockout (BERKO) mice to investigate the effects of BPS and BPF on insulin secretion, and the expression and activity of ion channels involved in β-cell function. BPS or BPF rapidly increased insulin release and diminished ATP-sensitive K⁺ (K_ATP) channel activity. Similarly, 48 h treatment with BPS or BPF enhanced insulin release and decreased the expression of several ion channel subunits in β-cells from WT mice, yet no effects were observed in cells from BERKO mice. PaPE-1, a ligand designed to preferentially trigger extranuclear-initiated ER pathways, mimicked the effects of bisphenols, suggesting the involvement of extranuclear-initiated ERβ pathways. Molecular dynamics simulations indicated differences in ERβ ligand-binding domain dimer stabilization and solvation free energy among different bisphenols and PaPE-1. Our data suggest a mode of action involving ERβ whose activation alters three key cellular events in β-cell, namely ion channel expression and activity, and insulin release. These results may help to improve the hazard identification of bisphenols.

Monitoring ligand-mediated helix 12 transitions within the human estrogen receptor α using bipartite tetracysteine display

Estrogen receptor α ligand-binding domains (ERα-LBD) expressing tetracysteine motifs bind FlAsH-EDT₂ upon transition of helix 12 (H12) to a folded state. Changes in fluorescence intensity allowed surveillance of ligand-mediated H12 transitions and facilitated the determination of FlAsH association rates (k_on) and apparent equilibrium dissociation constants (K_app) to ERα-LBDs in the presence of estrogenic ligands.

Computational investigations of the binding mechanism of novel benzophenone imine inhibitors for the treatment of breast cancer

4-Hydroxytamoxifen (4-OHT), the most common hormone used for the treatment of breast cancer, is a selective estrogen receptor modulator (SERM) inhibitor that acts as an antagonist in breast tissue and a partial agonist in the endometrium. However, the detailed molecular mechanism of 4-OHT structure modification has not been well investigated to date. Herein, molecular docking, molecular dynamics simulations and free energy calculations were performed to explore the mechanisms of the molecular interactions between newly designed benzophenone imines (BIs) and the three forms apo, antagonist and agonist of the human estrogen receptor hERα. The proposed inhibitors were designed by replacing the triarylethylene estrogenic scaffold found in 4-OHT with Schiff base triarylimine derivatives. The antiestrogen scaffold i.e. the O-alkyl side chain in 4-OHT was developed by incorporating an alanine amino acid side chain functionality into the triarylimine scaffold. Docking results reveal that the newly designed BIs bind to the hydrophobic open pocket of the apo and antagonist hERα conformations with higher affinity as compared to the natural and synthetic estrogen estradiol (E2) and 4-OHT. The analysis of the molecular dynamics simulation results based on six different systems of the best docked BI (5c) with hERα receptors demonstrates stable interactions, and the complex undergoes fewer conformational fluctuations in the open apo/antagonist hERα receptors as compared to the case of the closed agonist. In addition, the calculated binding free energies indicate that the main factor that contributes to the stabilization of the receptor-inhibitor complexes is hydrophobic interactions. This study suggests that the development of these Schiff base derivatives may be worth exploring for the preparation of new 4-OHT analogues.

Computational insights on agonist and antagonist mechanisms of estrogen receptor α induced by bisphenol A analogues

Structural analogues of bisphenol A (BPA) have become widely used as alternatives in BPA-free products. Most toxicological investigations have focused on the estrogenic activities of these analogues, which have been considered as potential environmental estrogens. However, recent studies revealed that certain BPA analogues could dramatically inhibit the proliferation of breast cancer cells, and exhibited strong anti-estrogenic effects compared with the antagonist 4-hydroxytamoxifen (OHT). Thus, we adopted computational models combining molecular dynamics simulations and binding free energy calculations to explore the underlying molecular basis of BPA analogues binding to estrogen receptor α (ERα). We also evaluated ligand-induced structural rearrangements of ERα at the atomic level. Conformational analyses showed that induced-fit H-bonding recognition by Thr347 was an important factor distinguishing antagonist from agonist BPA analogues. Moreover, antagonists of BPA analogues could indirectly induce the structural reposition of key helix 12 and produce an antagonistic conformation of ERα. Compared with OHT, the binding affinity of BPA analogues is stronger for antagonists than agonists. Taken together, we therefore propose computational indicators for screening of anti-estrogenic activities of BPA analogues, which may be beneficial for predicting the estrogenic or anti-estrogenic effects of BPA alternatives.

Prediction of Accurate Binding Modes Using Combination of Classical and Accelerated Molecular Dynamics and Free-Energy Perturbation Calculations: An Application to Toxicity Studies

Estimating the correct binding modes of ligands in protein-ligand complexes is crucial not only in the drug discovery process but also for elucidating potential toxicity mechanisms. In the current paper, we propose a computational modeling workflow using the combination of docking, classical molecular dynamics (cMD), accelerated molecular dynamics (aMD) and free-energy perturbation (FEP+ protocol) for identification of possible ligand binding modes. It was applied for investigation of selected perfluorocarboxyl acids (PFCAs) in the PPARγ nuclear receptor. Although both regular and induced fit docking failed to reproduce the experimentally determined binding mode of the ligands when docked into a non-native X-ray structure, cMD and aMD simulations successfully identified the most probable binding conformations. Moreover, multiple binding modes were identified for all of these compounds and the shorter-chain PFCAs continuously moved between a few energetically favorable binding conformations. On the basis of MD predictions of binding conformations, we applied the default and also redesigned FEP+ sampling protocols, which accurately reproduced experimental differences in the binding energies. Thus, the preliminary MD simulations can also provide helpful information about correct setup of the FEP+ calculations. These results show that the PFCA binding modes were accurately predicted and that the FEP+ protocol can be used to estimate free energies of binding of flexible ligands that are not typical druglike compounds. Our in silico workflow revealed the specific ligand-residue interactions within the ligand binding domain and the main characteristics of the PFCAs, and it was concluded that these compounds are week PPARγ partial agonists. This work also suggests a common pipeline for identification of ligand binding modes, ligand-protein dynamics description, and relative free-energy calculations.

A Computational Assay of Estrogen Receptor α Antagonists Reveals the Key Common Structural Traits of Drugs Effectively Fighting Refractory Breast Cancers

Somatic mutations of the Estrogen Receptor α (ERα) occur with an up to 40% incidence in ER sensitive breast cancer (BC) patients undergoing prolonged endocrine treatments. These polymorphisms are implicated in acquired resistance, disease relapse, and increased mortality rates, hence representing a current major clinical challenge. Here, multi-microseconds (12.5 µs) molecular dynamics simulations revealed that recurrent ERα polymorphisms (i. e. L536Q, Y537S, Y537N, D538G) (mERα) are constitutively active in their apo form and that they prompt the selection of an agonist (active)-like conformation even upon antagonists binding. Interestingly, our simulations rationalize, for the first time, the efficacy profile of (pre)clinically used Selective Estrogen Receptor Modulators/Downregulators (SERMs/SERDs) against these variants, enlightening, at atomistic level of detail, the key common structural traits needed by drugs able to effectively fight refractory BC types. This knowledge represents a key advancement for mechanism-based therapeutics targeting resistant ERα isoforms, potentially allowing the community to move a step closer to ‘precision medicine’ calibrated on patients’ genetic profiles and disease progression.

An alternative conformation of ERβ bound to estradiol reveals H12 in a stable antagonist position

The natural ligand 17β-estradiol (E2) is so far believed to induce a unique agonist-bound active conformation in the ligand binding domain (LBD) of the estrogen receptors (ERs). Both subtypes, ERα and ERβ, are transcriptionally activated in the presence of E2 with ERβ being somewhat less active than ERα under similar conditions. The molecular bases for this intriguing behavior are mainly attributed to subtype differences in the amino-terminal domain of these receptors. However, structural details that confer differences in the molecular response of ER LBDs to E2 still remain elusive. In this study, we present a new crystallographic structure of the ERβ LBD bound to E2 in which H12 assumes an alternative conformation that resembles antagonist ERs structures. Structural observations and molecular dynamics simulations jointly provide evidence that alternative ERβ H12 position could correspond to a stable conformation of the receptor under physiological pH conditions. Our findings shed light on the unexpected role of LBD in the lower functional response of ERβ subtype.

PPARγ helix 12 exhibits an antagonist conformation

PPARγ activation helix 12 can exist in an antagonist form: evidence from high-throughput accelerated molecular dynamics and metadynamics.

Computational investigation of the interaction mechanism between the estrogen related receptor α and its agonists

The agonists may promote the binding of coactivator PGC-1α to ERRα by stabilizing the conformation and the site of H12.

Structural and Dynamical Insight into PPARγ Antagonism: In Silico Study of the Ligand-Receptor Interactions of Non-Covalent Antagonists

The structural and dynamical properties of the peroxisome proliferator-activated receptor γ (PPARγ) nuclear receptor have been broadly studied in its agonist state but little is known about the key features required for the receptor antagonistic activity. Here we report a series of molecular dynamics (MD) simulations in combination with free energy estimation of the recently discovered class of non-covalent PPARγ antagonists. Their binding modes and dynamical behavior are described in details. Two key interactions have been detected within the cavity between helices H3, H11 and the activation helix H12, as well as with H12. The strength of the ligand-amino acid residues interactions has been analyzed in relation to the specificity of the ligand dynamical and antagonistic features. According to our results, the PPARγ activation helix does not undergo dramatic conformational changes, as seen in other nuclear receptors, but rather perturbations that occur through a significant ligand-induced reshaping of the ligand-receptor and the receptor-coactivator binding pockets. The H12 residue Tyr473 and the charge clamp residue Glu471 play a central role for the receptor transformations. Our results also demonstrate that MD can be a helpful tool for the compound phenotype characterization (full agonists, partial agonists or antagonists) when insufficient experimental data are available.