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

Abstract

Since estrogen-related receptor alpha (ERRa), one of three estrogen-related receptors, displays constitutively active transcriptional activities and important implications in both physiological and pathological processes of breast cancers, ERRa was recently recognized as a new target to fight breast cancers, and regulating the activity of ERRa with inverse agonists has thus become a promising new therapeutic strategy. A few inverse agonists cyclohexylmethyl-(1-p-tolyl-1H-indol-3-ylmethyl)-amine (compound1), thiadiazoacrylamide (XCT790), and 1-(2,5-diethoxy-benzyl)-3-phenyl-area analogues (compounds2and3) were reported to be capable of targeting ERRa. However, the detailed mechanism by which the inverse agonists deactivate ERRa 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 ERRa ligand binding domain, all-atom molecular dynamics (MD) simulations were firstly carried out for the complexes of inverse agonists and ERRa. The binding free energies were then calculated with MM-PBSA method to quantitatively discuss the binding of the inverse agonists with ERRa. 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 ERRa 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-1a binds to ERRa. Binding energy analysis indicates that compound3andXCT790bind more tightly to ERRa than compounds1and2, 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 compound1XCT790, and compound3have similar binding spectra to ERRa, primarily interacting with the residues of H3, H5, H6/H7 loop, and H11 helix, while compound2lacks 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.