The Caenorhabditis elegans DAF-12 nuclear receptor: structure, dynamics, and interaction with ligands

A structure for the ligand binding domain (LBD) of the DAF-12 receptor from Caenorhabditis elegans was obtained from the X-ray crystal structure of the receptor LBD from Strongyloides stercoralis bound to (25R)-Δ(7)-dafachronic acid (DA) (pdb:3GYU). The model was constructed in the presence of the ligand using a combination of Modeller, Autodock, and molecular dynamics (MD) programs, and then its dynamical behavior was studied by MD. A strong ligand binding mode (LBM) was found, with the three arginines in the ligand binding pocket (LBP) contacting the C-26 carboxylate group of the DA. The quality of the ceDAF-12 model was then evaluated by constructing several ligand systems for which the experimental activity is known. Thus, the dynamical behavior of the ceDAF-12 complex with the more active (25S)-Δ(7)-DA showed two distinct binding modes, one of them being energetically more favorable compared with the 25R isomer. Then the effect of the Arg564Cys and Arg598Met mutations on the (25R)-Δ(7)-DA binding was analyzed. The MD simulations showed that in the first case the complex was unstable, consistent with the lack of transactivation activity of (25R)-Δ(7)-DA in this mutant. Instead, in the case of the Arg598Met mutant, known to produce a partial loss of activity, our model predicted smaller effects on the LBM with a more stable MD trajectory. The model also showed that removal of the C-25 methyl does not impede the simultaneous strong interaction of the carboxylate with the three arginines, predicting that 27-nor-DAs are putative ceDAF-12 ligands.     

Molecular mechanism of ligand recognition by NR3 subtype glutamate receptors

NR3 subtype glutamate receptors have a unique developmental expression profile, but are the least well-characterized members of the NMDA receptor gene family, which have key roles in synaptic plasticity and brain development. Using ligand binding assays, crystallographic analysis, and all atom MD simulations, we investigate mechanisms underlying the binding by NR3A and NR3B of glycine and D-serine, which are candidate neurotransmitters for NMDA receptors containing NR3 subunits. The ligand binding domains of both NR3 subunits adopt a similar extent of domain closure as found in the corresponding NR1 complexes, but have a unique loop 1 structure distinct from that in all other glutamate receptor ion channels. Within their ligand binding pockets, NR3A and NR3B have strikingly different hydrogen bonding networks and solvent structures from those found in NR1, and fail to undergo a conformational rearrangement observed in NR1 upon binding the partial agonist ACPC. MD simulations revealed numerous interdomain contacts, which stabilize the agonist-bound closed-cleft conformation, and a novel twisting motion for the loop 1 helix that is unique in NR3 subunits.     

Investigation of the influence of external factors on the conformational dynamics of rhodopsin-like receptors by means of molecular dynamics simulation

Abstract

The study reports about the influence of binding of orthosteric ligands on the conformational dynamics of β-2-adrenoreceptor. Using molecular dynamics (MD) simulation, we found that there was a little fraction of active states of the receptor in its apo (ligand-free) ensemble. Analysis of MD trajectories indicated that such spontaneous activation of the receptor is accompanied by the motion in intracellular part of its alpha-helices. Thus, receptor’s constitutive activity directly results from its conformational dynamics. On the other hand, the binding of a full agonist resulted in a significant shift of the initial equilibrium towards its active state. Finally, the binding of the inverse agonist stabilized the receptor in its inactive state. It is likely that the binding of inverse agonists might be a universal way of constitutive activity inhibition in vivo. Our results indicate that ligand binding redistribute pre-existing conformational degrees of freedom (in accordance to the Monod–Wyman–Changeux Model) of the receptor rather than cause induced fit in it. Therefore, the ensemble of biologically relevant receptor conformations is encoded in its spatial structure, and individual conformations from that ensemble might be used by the cell in conformity with the physiological behavior.     

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.     

Molecular dynamics simulations of apo and holo forms of fatty acid binding protein 5 and cellular retinoic acid binding protein II reveal highly mobile protein,retinoic acid ligand,and water molecules

Structural and dynamic properties from a series of 300 ns molecular dynamics, MD, simulations of two intracellular lipid binding proteins, iLBPs, (Fatty Acid Binding Protein 5, FABP5, and Cellular Retinoic Acid Binding Protein II, CRABP-II) in both the apo form and when bound with retinoic acid reveal a high degree of protein and ligand flexibility. The ratio of FABP5 to CRABP-II in a cell may determine whether it undergoes natural apoptosis or unrestricted cell growth in the presence of retinoic acid. As a result, FABP5 is a promising target for cancer therapy. The MD simulations presented here reveal distinct differences in the two proteins and provide insight into the binding mechanism. CRABP-II is a much larger, more flexible protein that closes upon ligand binding, where FABP5 transitions to an open state in the holo form. The traditional understanding obtained from crystal structures of the gap between two β-sheets of the β-barrel common to iLBPs and the α-helix cap that forms the portal to the binding pocket is insufficient for describing protein conformation (open vs. closed) or ligand entry and exit. When the high degree of mobility between multiple conformations of both the ligand and protein are examined via MD simulation, a new mode of ligand motion that improves understanding of binding dynamics is revealed.     

Study of the impact of the T877A mutation on ligand‐induced helix‐12 positioning of the androgen receptor resulted in design and synthesis of novel antiandrogens

Antiandrogen flutamide, an antagonist of the wild‐type androgen receptor (AR), is used in the clinics for treating metastatic prostate cancer. However, the T877A mutated AR is paradoxically activated by hydroxyflutamide, an active form of flutamide. Despite of crystallographic studies, how the T877A mutation results in antagonist‐agonist conversion of hydroxyflutamide remains a puzzle. Here, started from a structural model of the apo form of AR ligand‐binding domain (AR‐LBD), we have investigated the impact of the T877A mutation on ligand‐induced helix‐12 positioning by replica‐exchange molecular dynamics (REMD) simulations with an unique protocol, which is capable of simulating the H12 dynamics and keeping the main body of AR‐LBD unchanged. Specifically, (i) we have computationally demonstrated that on the binding of hydroxyflutamide, the apo form of H12 rearranges into the agonistic form in the T877A mutant, but into the antagonistic forms in the wild‐type receptor, shedding light on hydroxyflutamide agonism/antagonism; (ii) By REMD simulations, we have predicted antiandrogen SC184 is a non‐agonist of the T877A mutant. This was confirmed by luciferase assays; and (iii) on the basis of the binding modes of hydroxyflutamide and SC184 from the simulations, we designed a novel flutamide derivative called SC333, which was subsequently predicted to be a pure antagonist of the T877A mutant. We then synthesized and experimentally confirmed SC333 is a pan‐antiandrogen effective against the wild‐type and the T877A and W741C mutated ARs, showing low micromolar cytotoxicity in LNCaP cells. Importantly, we demonstrated that distribution of the H12 conformations from REMD simulations is correlated with ligand agonist/antagonist activity. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

核雌激素受体配体结合域的晶体结构特性与功能

雌激素或类雌激素活性物质通过细胞核雌激素受体（nuclear estrogen receptor, nER）通路发挥相应的生理性作用。当这些配体被nER的配体结合域（ligand binding domain, LBD）识别后进入疏水性配体结合空腔内并引起受体构象发生改变,使得原先处于高度活动性的helix 12(H12)被固定从而进一步稳定空腔结构|同时nER也能通过招募一系列辅助调节因子及其他共调节蛋白质,最终调控基因转录。但是,由于不同的配体和受体结合形成的晶体结构并不完全相同,导致这些复合体具有不同的性质,从而影响基因的转录活性。本文综述了nER配体结合域及结合配体后形成的相应晶体结构与活性以及不同配体对受体结构和基因转录的影响。     

Subunit-specific modulation of [(3)H]MK-801 binding to NMDA receptors mediated by dopamine receptor ligands in rodent brain

Dopamine D1 receptor (D1R) ligands may directly interact with the NMDA receptor (NMDAR), but detailed knowledge about this effect is lacking. Here we identify D1R ligands that directly modulate NMDARs and examine the contributions of NR2A and NR2B subunits to these interactions. Binding of the open channel blocker [(3)H]MK-801 in membrane preparations from rat- and mouse brain was used as a biochemical measure of the functional state of the NMDAR channel. We show that both D1R agonist A-68930 and dopamine receptor D2 antagonist haloperidol can decrease [(3)H]MK-801 binding with increased potency in membranes from the NR2A(-/-) mice (i.e. in membranes containing NR2B only), as compared to the inhibition obtained in wild-type membranes. Further, a wide range of D1R agonists such as A-68930, SKF-83959, SKF-83822, SKF-38393 and dihydrexidine were able to decrease [(3)H]MK-801 binding, all showing half maximal inhibitory concentrations ~20 μM, and with significant effects occurring at or above 1 μM. With membranes from D1R(-/-) mice, we demonstrate that these effects occurred through a D1R-independent mechanism. Our results demonstrate that dopamine receptor ligands can selectively influence NR2B containing NMDARs, and we characterize direct inhibitory NMDAR effects by different D1R ligands.     

Molecular dynamics simulations of human LRH-1: the impact of ligand binding in a constitutively active nuclear receptor

The liver receptor homologue 1 (LRH-1 (NR5A2)) belongs to the orphan nuclear receptor family, indicating that initially no ligand was known. Although recent studies have shown that ligand binding can be obtained, the biological relevance remains elusive. Here, we modify the observed X-ray ligand into a biologically more significant phospholipid (phosphatidylserine, PS) present in human, to study, by molecular dynamics (MD) simulations, the impact of the ligand on the receptor and the interaction with different cofactor peptides. Furthermore, we characterize the interactions between receptor and the cofactor peptides of DAX-1 (NR0B1), Prox1 and SHP LXXLL box 1 and 2 (NR0B2) in terms of specificity. Our MD simulation results show different interaction patterns for the SHP box2 compared to DAX-1, PROX1 and SHP box1. SHP box2 shows specific interactions at its more C-terminal end while the other investigated peptides show specific interactions at several positions but particularly at the +2 site. The peptide +2 side chain interacts with a charged amino acid of the receptor, in hLRH-1 Asp372. Together with the charge clamp residues Arg361 and Glu534, Asp372 forms a triangle shaped charge clamp responsible for peptide orientation and increased affinity. The binding of the PS ligand causes no overall structural changes of the receptor but affects the interactions with cofactor peptides. The cofactor peptides from SHP decrease its interaction with the receptor upon ligand binding while DAX-1 and PROX1 are unchanged or increase. The diverse ligand binding response of the cofactor provides an opportunity for drug design with the possibility to create agonist ligands to modify cofactor interaction.