Exploring the molecular basis of action of ring D aromatic steroidal antiestrogens

Salpichrolides are natural plant steroids that contain an unusual six‐membered aromatic ring D. We recently reported that some of these compounds, and certain analogs with a simplified side chain, exhibited antagonist effects toward the human estrogen receptor (ER), a nuclear receptor whose endogenous ligand has an aromatic A ring (estradiol). Drugs acting through the inhibition or modulation of ERs are frequently used as a hormonal therapy for ER(+) breast cancer. Previous results suggested that the aromatic D ring was a key structural motif for the observed activity; thus, this modified steroid nucleus may provide a new scaffold for the design of novel antiestrogens. Using molecular dynamics (MD) simulation we have modeled the binding mode of the natural salpichrolide A and a synthetic analog with an aromatic D ring within the ERα. These results taken together with the calculated energetic contributions associated to the different ligand‐binding modes are consistent with a preferred inverted orientation of the steroids in the ligand‐binding pocket with the aromatic ring D occupying a position similar to that observed for the A ring of estradiol. Major changes in both dynamical behavior and global positioning of H11 caused by the loss of the ligand–His524 interaction might explain, at least in part, the molecular basis of the antagonism exhibited by these compounds. Using steered MD we also found a putative unbinding pathway for the steroidal ligands through a cavity formed by residues in H3, H7, and H11, which requires only minor changes in the overall receptor conformation. Proteins 2015; 83:1297–1306. © 2015 Wiley Periodicals, Inc.     

Identification of a GPER/GPR30 antagonist with improved estrogen receptor counterselectivity

GPER/GPR30 is a seven-transmembrane G protein-coupled estrogen receptor that regulates many aspects of mammalian biology and physiology. We have previously described both a GPER-selective agonist G-1 and antagonist G15 based on a tetrahydro-3H-cyclopenta[c]quinoline scaffold. The antagonist lacks an ethanone moiety that likely forms important hydrogen bonds involved in receptor activation. Computational docking studies suggested that the lack of the ethanone substituent in G15 could minimize key steric conflicts, present in G-1, that limit binding within the ERα ligand binding pocket. In this report, we identify low-affinity cross-reactivity of the GPER antagonist G15 to the classical estrogen receptor ERα. To generate an antagonist with enhanced selectivity, we therefore synthesized an isosteric G-1 derivative, G36, containing an isopropyl moiety in place of the ethanone moiety. We demonstrate that G36 shows decreased binding and activation of ERα, while maintaining its antagonist profile towards GPER. G36 selectively inhibits estrogen-mediated activation of PI3K by GPER but not ERα. It also inhibits estrogen- and G-1-mediated calcium mobilization as well as ERK1/2 activation, with no effect on EGF-mediated ERK1/2 activation. Similar to G15, G36 inhibits estrogen- and G-1-stimulated proliferation of uterine epithelial cells in vivo. The identification of G36 as a GPER antagonist with improved ER counterselectivity represents a significant step towards the development of new highly selective therapeutics for cancer and other diseases.     

Unbinding pathways of an agonist and an antagonist from the 5-HT3 receptor

The binding sites of 5-HT3 and other Cys-loop receptors have been extensively studied, but there are no data on the entry and exit routes of ligands for these sites. Here we have used molecular dynamics simulations to predict the pathway for agonists and antagonists exiting from the 5-HT3 receptor binding site. The data suggest that the unbinding pathway follows a tunnel at the interface of two subunits, which is approximately 8 A long and terminates approximately 20 A above the membrane. The exit routes for an agonist (5-HT) and an antagonist (granisetron) were similar, with trajectories toward the membrane and outward from the ligand binding site. 5-HT appears to form many hydrogen bonds with residues in the unbinding pathway, and experiments show that mutating these residues significantly affects function. The location of the pathway is also supported by docking studies of granisetron, which show a potential binding site for granisetron on the unbinding route. We propose that leaving the binding pocket along this tunnel places the ligands close to the membrane and prevents their immediate reentry into the binding pocket. We anticipate similar exit pathways for other members of the Cys-loop receptor family.     

Validation of Alexa-647-ATP as a powerful tool to study P2X receptor ligand binding and desensitization

Ion channel opening and desensitization is a fundamental process in neurotransmission. The ATP-gated P2X1 receptor (P2X1R) shows rapid and long-lasting desensitization upon agonist binding. This makes the electrophysiological investigation of its desensitization process, agonist unbinding, and recovery from desensitization a challenging task. Here, we show that the fluorescent agonist Alexa-647-ATP is a potent agonist at the P2X1R and a versatile tool to directly visualize agonist binding and unbinding. We demonstrate that the long-lasting desensitization of the P2X1R is due to both slow unbinding of agonist from the desensitized receptor and agonist mediated receptor internalization. Furthermore, the unbinding of the agonist Alexa-647-ATP from the desensitized receptor is accelerated in the continuous presence of competitive ligand. Modeling of our data indicates that three agonist molecules are required to drive the receptor into desensitization. Direct visualization of ligand unbinding from the desensitized receptor demonstrates the cooperativity of this process.     

Microscopic kinetics and energetics distinguish GABA(A) receptor agonists from antagonists

Although agonists and competitive antagonists presumably occupy overlapping binding sites on ligand-gated channels, these interactions cannot be identical because agonists cause channel opening whereas antagonists do not. One explanation is that only agonist binding performs enough work on the receptor to cause the conformational changes that lead to gating. This idea is supported by agonist binding rates at GABA(A) and nicotinic acetylcholine receptors that are slower than expected for a diffusion-limited process, suggesting that agonist binding involves an energy-requiring event. This hypothesis predicts that competitive antagonist binding should require less activation energy than agonist binding. To test this idea, we developed a novel deconvolution-based method to compare binding and unbinding kinetics of GABA(A) receptor agonists and antagonists in outside-out patches from rat hippocampal neurons. Agonist and antagonist unbinding rates were steeply correlated with affinity. Unlike the agonists, three of the four antagonists tested had binding rates that were fast, independent of affinity, and could be accounted for by diffusion- and dehydration-limited processes. In contrast, agonist binding involved additional energy-requiring steps, consistent with the idea that channel gating is initiated by agonist-triggered movements within the ligand binding site. Antagonist binding does not appear to produce such movements, and may in fact prevent them.     

How maltose influences structural changes to bind to maltose‐binding protein: Results from umbrella sampling simulation

A well‐studied periplasmic‐binding protein involved in the abstraction of maltose is maltose‐binding protein (MBP), which undergoes a ligand‐induced conformational transition from an open (ligand‐free) to a closed (ligand‐bound) state. Umbrella sampling simulations have been us to estimate the free energy of binding of maltose to MBP and to trace the potential of mean force of the unbinding event using the center‐of‐mass distance between the protein and ligand as the reaction coordinate. The free energy thus obtained compares nicely with the experimentally measured value justifying our theoretical basis. Measurement of the domain angle (N‐terminal‐domain – hinge – C‐terminal‐domain) along the unbinding pathway established the existence of three different states. Starting from a closed state, the protein shifts to an open conformation during the initial unbinding event of the ligand then resides in a semi‐open conformation and later resides predominantly in an open‐state. These transitions along the ligand unbinding pathway have been captured in greater depth using principal component analysis. It is proposed that in mixed‐model, both conformational selection and an induced‐fit mechanism combine to the ligand recognition process in MBP. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Phytoestrogens genistein and daidzin enhance the acetylcholinesterase activity of the rat pheochromocytoma cell line PC12 by binding to the estrogen receptor

Some compounds derived from plants have been known to possess estrogenic properties and can thus alter the physiology of higher organisms. Genistein and daidzin are examples of these phytoestrogens, which have recently been the subject of extensive research. In this study, genistein and daidzin were found to enhance the acetylcholinesterase (AChE) activity of the rat neuronal cell line PC12 at concentrations as low as 0.08 μM by binding to the estrogen receptor (ER). Results have shown that this enhancement was effectively blocked by the known estrogen receptor antagonist tamoxifen, indicating the involvement of the ER in AChE induction. That genistein and daidzin are estrogenic were confirmed in a cell proliferation assay using the human breast cancer cell line MCF7. This proliferation was also blocked by tamoxifen, again indicating the involvement of the ER. On the other hand, incubating the PC12 cells in increasing concentrations of 17 β-estradiol (E2) did not lead to enhanced AChE activity, even in the presence of genistein or daidzin. This suggests that mere binding of an estrogenic compound to the ER does not necessarily lead to enhanced AChE activity. Moreover, the effect of the phytoestrogens on AChE activity cannot be expressed in the presence of E2 since they either could not compete with the natural ligand in binding to the ER or that E2 down-regulates its own receptor. This study clearly suggests that genistein and daidzin enhance AChE activityin PC12 cells by binding to the ER; however, the actual mechanism of enhancement is not known. This revised version was published online in August 2006 with corrections to the Cover Date.     

Differential DNA-binding abilities of estrogen receptor occupied with two classes of antiestrogens: studies using human estrogen receptor overexpressed in mammalian cells.

We have developed a transient transfection system using the Cytomegalovirus (CMV) promoter to express the human estrogen receptor (ER) at very high levels in COS-1 cells and have used it to study the interaction of agonist and antagonist receptor complexes with estrogen response element (ERE) DNA. ER can be expressed to levels of 20-40 pmol/mg or 0.2-0.3% of total soluble protein and all of the soluble receptor is capable of binding hormone. The ER binds estradiol with high affinity (Kd 0.2 nM), and is indistinguishable from native ER in that the receptor is capable of recognizing its cognate DNA response element with high affinity, and of transactivating a transgene in an estradiol-dependent manner. Gel mobility shift assays reveal interesting ligand-dependent differences in the binding of receptor complexes to ERE DNA. Receptors occupied by estradiol or the type I antiestrogen transhydroxytamoxifen bind to DNA response elements when exposed to the ligand in vitro or in vivo. Likewise, receptors exposed to the type II antiestrogen ICI 164,384 in vitro bind to ERE DNA. However, when receptor exposure to ICI 164,384 is carried out in vivo, the ER-ICI 164,384 complexes do not bind to ERE DNA, or do so only weakly. This effect is not reversed by subsequent incubation with estradiol in vitro, but is rapidly reversible by in vivo estradiol exposure of intact COS-1 cells. This suggests there may be some cellular process involved in the mechanism of antagonism by the pure antiestrogen ICI 164,384, which is not observed in cell-free extracts.     

Conformational selection in protein binding and function

Protein binding and function often involves conformational changes. Advanced nuclear magnetic resonance (NMR) experiments indicate that these conformational changes can occur in the absence of ligand molecules (or with bound ligands), and that the ligands may “select” protein conformations for binding (or unbinding). In this review, we argue that this conformational selection requires transition times for ligand binding and unbinding that are small compared to the dwell times of proteins in different conformations, which is plausible for small ligand molecules. Such a separation of timescales leads to a decoupling and temporal ordering of binding/unbinding events and conformational changes. We propose that conformational‐selection and induced‐change processes (such as induced fit) are two sides of the same coin, because the temporal ordering is reversed in binding and unbinding direction. Conformational‐selection processes can be characterized by a conformational excitation that occurs prior to a binding or unbinding event, while induced‐change processes exhibit a characteristic conformational relaxation that occurs after a binding or unbinding event. We discuss how the ordering of events can be determined from relaxation rates and effective on‐ and off‐rates determined in mixing experiments, and from the conformational exchange rates measured in advanced NMR or single‐molecule fluorescence resonance energy transfer experiments. For larger ligand molecules such as peptides, conformational changes and binding events can be intricately coupled and exhibit aspects of conformational‐selection and induced‐change processes in both binding and unbinding direction.     

Unraveling Subunit Cooperativity in Homotetrameric HCN2 Channels

In a multimeric receptor protein, the binding of a ligand can modulate the binding of a succeeding ligand. This phenomenon, called cooperativity, is caused by the interaction of the receptor subunits. By using a complex Markovian model and a set of parameters determined previously, we analyzed how the successive binding of four ligands leads to a complex cooperative interaction of the subunits in homotetrameric HCN2 pacemaker channels. The individual steps in the model were characterized by Gibbs free energies for the equilibria and activation energies, specifying the affinity of the binding sites and the transition rates, respectively. Moreover, cooperative free energies were calculated for each binding step in both the closed and the open channel. We show that the cooperativity sequence positive-negative-positive determined for the binding affinity is generated by the combined effect of very different cooperativity sequences determined for the binding and unbinding rates, which are negative-negative-positive and no-negative-no, respectively. It is concluded that in the ligand-induced activation of HCN2 channels, the sequence of cooperativity based on the binding affinity is caused by two even qualitatively different sequences of cooperativity that are based on the rates of ligand binding and unbinding.     

Forced unbinding of GPR17 ligands from wild type and R255I mutant receptor models through a computational approach

Background  

GPR17 is a hybrid G-protein-coupled receptor (GPCR) activated by two unrelated ligand families, extracellular nucleotides and cysteinyl-leukotrienes (cysteinyl-LTs), and involved in brain damage and repair. Its exploitment as a target for novel neuro-reparative strategies depends on the elucidation of the molecular determinants driving binding of purinergic and leukotrienic ligands. Here, we applied docking and molecular dynamics simulations (MD) to analyse the binding and the forced unbinding of two GPR17 ligands (the endogenous purinergic agonist UDP and the leukotriene receptor antagonist pranlukast from both the wild-type (WT) receptor and a mutant model, where a basic residue hypothesized to be crucial for nucleotide binding had been mutated (R255I) to Ile.