Ligand pathways in estrogen-related receptors

The three subtypes of estrogen-related receptors ERRα, ERRβ, and ERRγ are nuclear receptors mediating metabolic processes in various tissues such as the skeletal muscle, fat tissue, bone, and liver. Although the knowledge on their physiological ligands is limited, they have been implicated as drug targets for important indications including diabetes, cardiovascular diseases, and osteoporosis. As in other nuclear receptors, their ligand binding pocket is buried within the core of the receptor and connected to its surrounding by ligand pathways. Here, we investigated these pathways with conventional molecular dynamics as well as metadynamics simulations to reveal their distribution and their capability to facilitate ligand translocation. Dependent on the ERR subtype and the conformational state of the receptor, we could detect different pathways to be favored. Overall, the results suggested pathways IIIa and IIIb to be favored in the agonistic conformation, while antagonists preferred pathways I, II, and V. Along the pathways, the ligands passed different gating mechanisms of the receptor, including groups of protein residues as well as whole secondary structure elements, to leave the binding site. Even though these pathways are suggested to influence ligand specificity of the receptors and their elucidation might advance rational drug design, they have not yet been studied in ERRs.Communicated by Ramaswamy H. Sarma.

Nuclear receptors in the kidney during health and disease

Over the last 30 years, nuclear receptors (NRs) have been increasingly recognized as key modulators of systemic homeostasis and as contributing factors in many diseases. In the kidney, NRs play numerous important roles in maintaining homeostasis-many of which continue to be unraveled. As "master regulators", these important transcription factors integrate and coordinate many renal processes such as circadian responses, lipid metabolism, fatty acid oxidation, glucose handling, and inflammatory responses. The use of recently-developed genetic tools and small molecule modulators have allowed for detailed studies of how renal NRs contribute to kidney homeostasis. Importantly, while NRs are intimately involved in proper kidney function, they are also implicated in a variety of renal diseases such as diabetes, acute kidney injury, and other conditions such as aging. In the last 10 years, our understanding of renal disease etiology and progression has been greatly shaped by knowledge regarding how NRs are dysregulated in these conditions. Importantly, NRs have also become attractive therapeutic targets for attenuation of renal diseases, and their modulation for this purpose has been the subject of intense investigation. Here, we review the role in health and disease of six key renal NRs including the peroxisome proliferator-activated receptors (PPAR), estrogen-related receptors (ERR), the farnesoid X receptors (FXR), estrogen receptors (ER), liver X receptors (LXR), and vitamin D receptors (VDR) with an emphasis on recent findings over the last decade. These NRs have generated a wealth of data over the last 10 years that demonstrate their crucial role in maintaining normal renal homeostasis as well as their capacity to modulate disease progression.

Insights into the interaction mechanisms of estrogen-related receptor alpha (ERRα) with ligands via molecular dynamics simulations

Estrogen-related receptor alpha (ERRα), a member of nuclear receptors (NRs), participates in energy metabolism. Recent experiments identified that several agonists to increase the activity of ERRα, which have a therapeutic effect in improving insulin sensitivity and lowering blood glucose levels. However, the detailed molecular mechanism about how the ligands affect the structure of ERRα remains elusive. To better understand the conformational change of ERRα complexed with agonists and inverse agonists, unbiased molecular dynamics (MD) simulations were performed on the ligand binding domain of ERRα (ERRα-LBD) bound with different ligands. According to the results, the ERRα-agonist interactions were more stable in the presence of the peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α). More importantly, we observed that the binding of inverse agonists would decrease the stability of helix 12 (H12) of ERRα. Moreover, we suggested that Phe232 and Phe414 should be key residues in the interaction pathway from ligands to H12, which provided a possible explanation about how ligands impact the structure of ERRα. These results would provide insights into the design of novel and efficient agonists of ERRα to treat metabolic diabetes.Communicated by Ramaswamy H. Sarma.

Subresidue-Resolution Footprinting of Ligand-Protein Interactions by Carbene Chemistry and Ion Mobility-Mass Spectrometry

The knowledge of ligand-protein interactions is essential for understanding fundamental biological processes and for the rational design of drugs that target such processes. Carbene footprinting efficiently labels proteinaceous residues and has been used with mass spectrometry (MS) to map ligand-protein interactions. Nevertheless, previous footprinting studies are typically performed at the residue level, and therefore, the resolution may not be high enough to couple with conventional crystallography techniques. Herein we developed a subresidue footprinting strategy based on the discovery that carbene labeling produces subresidue peptide isomers and the intensity changes of these isomers in response to ligand binding can be exploited to delineate ligand-protein topography at the subresidue level. The established workflow combines carbene footprinting, extended liquid chromatographic separation, and ion mobility (IM)-MS for efficient separation and identification of subresidue isomers. Analysis of representative subresidue isomers located within the binding cleft of lysozyme and those produced from an amyloid-β segment have both uncovered structural information heretofore unavailable by residue-level footprinting. Lastly, a "real-world" application shows that the reactivity changes of subresidue isomers at Phe399 can identify the interactive nuances between estrogen-related receptor α, a potential drug target for cancer and metabolic diseases, with its three ligands. These findings have significant implications for drug design. Taken together, we envision the subresidue-level resolution enabled by IM-MS-coupled carbene footprinting can bridge the gap between structural MS and the more-established biophysical tools and ultimately facilitate diverse applications for fundamental research and pharmaceutical development.

Exploring the binding mode and thermodynamics of inverse agonists against estrogen-related receptor alpha

Since estrogen-related receptor alpha (ERRα), one of three estrogen-related receptors, displays constitutively active transcriptional activities and important implications in both physiological and pathological processes of breast cancers, ERRα was recently recognized as a new target to fight breast cancers, and regulating the activity of ERRα with inverse agonists has thus become a promising new therapeutic strategy. A few inverse agonists cyclohexylmethyl-(1-p-tolyl-1H-indol-3-ylmethyl)-amine (compound 1), thiadiazoacrylamide (XCT790), and 1-(2,5-diethoxy-benzyl)-3-phenyl-area analogues (compounds 2 and 3) were reported to be capable of targeting ERRα. However, the detailed mechanism by which the inverse agonists deactivate ERRα remains unclear, especially in the aspects of quantitative binding and hot spot residues. Therefore, to gain insights into the interaction modes between inverse agonists and ERRα ligand binding domain, all-atom molecular dynamics (MD) simulations were firstly carried out for the complexes of inverse agonists and ERRα. The binding free energies were then calculated with MM-PBSA method to quantitatively discuss the binding of the inverse agonists with ERRα. The binding affinities were finally decomposed to per-residue contributions to identify the hot spot residues as well as assess their role in the binding mechanism. MD simulations show that the inverse agonists stretch downwards into the ERRα ligand binding pocket (LBP) formed by H3 and H11 helices, and upon the binding H12 adopts a well-defined position in the coactivator groove, where PGC-1α binds to ERRα. Binding energy analysis indicates that compound 3 and XCT790 bind more tightly to ERRα than compounds 1 and 2, and the energy difference mainly results from the contribution of van der Waals interaction. Both binding mode analysis and affinity decomposition per-residue indicate that compound 1, XCT790, and compound 3 have similar binding spectra to ERRα, primarily interacting with the residues of H3, H5, H6/H7 loop, and H11 helix, while compound 2 lacks a significant interaction with the H5 region. The hot spot residues significantly binding to the three inverse agonists in common include Leu324, Phe328, Phe382, Leu398, Phe495, and Leu500. It is essential for an effective inverse agonist to strongly bind with the aromatic ring cluster consisting of Phe328(H3), Phe495(H11), and Phe382(H5/H6 loop) as well as Leu500.

All-atom MD simulations were for the first time carried out for the complexes of inverse agonists and ERRα, and their binding free energies were also calculated with MM-PBSA to quantitatively discuss the binding of the inverse agonists with ERRα.

Computational insights into the interaction mechanisms of estrogen-related receptor alpha with endogenous ligand cholesterol

Estrogen-related receptor alpha (ERRα) has attracted increasing concerns. ERRα, orphan nuclear receptor, plays important roles in energy metabolism. Therefore, small molecule agonists of ERRα could be a potential therapeutic strategy in the treatment of metabolic diseases such as diabetes. Recently, Wei et al. identified cholesterol as the endogenous agonist of ERRα. However, the detailed molecular mechanism of cholesterol bound with ERRα remains ambiguous. Thus, in this study molecular docking and molecular dynamics (MD) simulations were performed to characterize how cholesterol affects the behavior of ERRα. Based on the results, we found that a proven residue Phe232 and others including Leu228, Glu235, Arg276, and Phe399 were key residues to ligand binding. A hydrogen-bonding interaction between cholesterol and Glu235 ensured the orientation of the ligand in the binding pocket, while hydrophobic interactions between cholesterol and the above-mentioned residues promoted the stability of ERRα-cholesterol complex. In the presence of the proliferator-activated receptor γ coactivator 1α (PGC-1α), the cholesterol-ERRα interaction became more stable. Interestingly, we observed that cholesterol facilitated the binding of ERRα with its coactivator PGC-1α via stabilizing the conformation of helix 12 and the interaction surface of ERRα/PGC-1α. Overall, these findings would be valuable for the future rational design of novel ERRα agonists.