Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Analysis of Agonist and Antagonist Effects on Thyroid Hormone Receptor Conformation by Hydrogen/Deuterium Exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH(3)). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH(3), increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH(3.) We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011)

On the Denaturation Mechanisms of the Ligand Binding Domain of Thyroid Hormone Receptors

The ligand binding domain (LBD) of nuclear hormone receptors adopts a very compact, mostly alpha-helical structure that binds specific ligands with very high affinity. We use circular dichroism spectroscopy and high-temperature molecular dynamics Simulations to investigate unfolding of the LBDs of thyroid hormone receptors (TRs). A molecular description of the denaturation mechanisms is obtained by molecular dynamics Simulations of the TR alpha and TR beta LBDs in the absence and in the presence of the natural ligand Triac. The Simulations Show that the thermal unfolding of the LBD starts with the loss of native contacts and secondary Structure elements, while the Structure remains essentially compact, resembling a molten globule state. This differs From most protein denaturation simulations reported to date and suggests that the folding mechanism may start with the hydrophobic collapse of the TR LBDs. Our results reveal that the stabilities of the LBDs of the TR alpha and TR beta Subtypes are affected to different degrees by the binding of the isoform selective ligand Triac and that ligand binding confers protection against thermal denaturation and unfolding in a subtype specific manner. Our Simulations indicate two mechanisms by which the ligand stabilizes the LBD: (1) by enhancing the interactions between H8 and H 11, and the interaction of the region between H I and the Omega-loop with the core of the LBD, and (2) by shielding the hydrophobic H6 from hydration.

Adenosine Modulatesa alpha(2)-Adrenergic Receptors within Specific Subnuclei of the Nucleus Tractus Solitarius in Normotensive and Spontaneously Hypertensive Rats

Adenosine Is known to modulate neuronal activity within the nucleus tractus solitarius (NTS). The modulatory effect of adenosine A, receptors (A(1R)) on alpha(2)-adrenoceptors (Adr(2R)) was evaluated using quantitative radioautography within NTS subnuclei and using neuronal culture of normotensive (WKY) and spontaneously hypertensive rats (SHR). Radioautography was used in a saturation experiment to measure Adr2R binding parameters (B(max), K(d)) In the presence of 3 different concentrations of N(6)-cyclopentyladenosine (CPA), an A(1R) agonist. Neuronal culture confirmed our radioautographic results. [(3)H]RX821002, an Adr(2R) antagonist, was used as a ligand for both approaches. The dorsomedial/dorsolateral subnucleus of WKY showed an increase in B(max) values (21%) Induced by 10 nmol/L of CPA. However, the subpostremal subnucleus showed a decrease in Kd values (24%) induced by 10 nmol/L of CPA. SHR showed the same pattern of changes as WKY within the same subnuclei; however, the modulatory effect of CPA was induced by I nmol/L (increased B(max), 17%; decreased K(d), 26%). Cell culture confirmed these results, because 10(-5) and 10(-7) mol/L of CPA promoted an Increase in [3H]RX821002 binding of WKY (53%) and SHR cells (48%), respectively. DPCPX, an AIR antagonist, was used to block the modulatory effect promoted by CPA with respect to Adr2R binding. In conclusion, our study shows for the first time an interaction between A(1R) that increases the binding of Adr2R within specific subnuclei of the NTS. This may be important In understanding the complex autonomic response induced by adenosine within the NTS. In addition, changes in interactions between receptors might be relevant to understanding the development of hypertension. (Hypertens Res 2008; 31: 2177-2186)

Actions of octocoral and tobacco cembranoids on nicotinic receptors

Nicotinic acetylcholine receptors (AChRs) are pentameric proteins that form agonist-gated cation channels through the plasma membrane. AChR agonists and antagonists are potential candidates for the treatment of neurodegenerative diseases. Cembranoids are naturally occurring diterpenoids that contain a 14-carbon ring. These diterpenoids interact with AChRs in complex ways: as irreversible inhibitors at the agonist sites, as noncompetitive inhibitors, or as positive modulators, but no cembranoid was ever shown to have agonistic activity on AChRs. The cembranoid eupalmerin acetate displays positive modulation of agonist-induced currents in the muscle-type AChR and in the related gamma-aminobutyric acid (GABA) type A receptor. Moreover, cembranoids display important biological effects, many of them mediated by nicotinic receptors. Cembranoids from tobacco are neuroprotective through a nicotinic anti-apoptotic mechanism preventing excitotoxic neuronal death which in part could result from anti-inflammatory properties of cembranoids. Moreover, tobacco cembranoids also have anti-inflammatory properties which could enhance their neuroprotective properties. Cembranoids from tobacco affect nicotine-related behavior: they increase the transient initial ataxia caused by first nicotine injection into naive rats and inhibit the expression of locomotor sensitization to repeated injections of nicotine. In addition, cembranoids are known to act as anti-tumor compounds. In conclusion, cembranoids provide a promising source of lead drugs for many clinical areas, including neuroprotection, smoking-cessation, and anti-cancer therapies. (C) 2009 Elsevier Ltd. All rights reserved.

Estrogen receptor: Structural differences and potential implications on selectivity examined by the GRID/CPCA approach

Selective Estrogen Receptor Modulators ( SERMs) have been developed, but the selectivity towards the subtypes ( ER or ER is not well understood. Based on three-dimensional structural properties of ligand binding domains, a model that takes into account this aspect was developed via molecular interaction fields and consensus principal component analysis (GRID/CPCA).