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.

Molecular Basis of the Thermostability and Thermophilicity of Laminarinases: X-ray Structure of the Hyperthermostable Laminarinase from Rhodothermus marinus and Molecular Dynamics Simulations

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 angstrom resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis, for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting beta-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Prigogine-Defay Ratio for an Ionic Glass-Former: Molecular Dynamics Simulations

The pressure dependence of the glass-transition temperature, T(g)(P), of the ionic glass-former 2Ca(NO(3))(2) center dot 3KNO(3), CKN, has been obtained by molecular dynamics (MD) simulations The liquid-glass difference of thermal expansivity, Delta alpha, heat capacity, Delta C(p), and isothermal compressibility, Delta kappa, have been calculated as a function of pressure. It has been found that the Ehrenfest relation dT(g)/dP = TV Delta alpha/Delta C(p) predicts the pressure dependence of T, but the other Ehrenfest relation, dT(g)/dP = Delta kappa/Delta alpha, does not. Consequently, the Prigogine-Defay ratio, Pi = Delta C(p)Delta kappa/TV Delta alpha(2), is Pi similar to 1.2 at low pressures, but increases 1 order of magnitude at high pressures. The pressure dependence of the Prigogine-Defay ratio is interpreted in light of recent explanations for the finding Pi > 1.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

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.

Selectivity in the mechanism of action of antimicrobial mastoparan peptide Polybia-MP1

Many potent antimicrobial peptides also present hemolytic activity, an undesired collateral effect for the therapeutic application. Unlike other mastoparan peptides, Polybia-MP1 (IDWKKLLDAAKQIL), obtained from the venom of the social wasp Polybia paulista, is highly selective of bacterial cells. The study of its mechanism of action demonstrated that it permeates vesicles at a greater rate of leakage on the anionic over the zwitterionic, impaired by the presence of cholesterol or cardiolipin; its lytic activity is characterized by a threshold peptide to lipid molar ratio that depends on the phospholipid composition of the vesicles. At these particular threshold concentrations, the apparent average pore number is distinctive between anionic and zwitterionic vesicles, suggesting that pores are similarly formed depending on the ionic character of the bilayer. To prospect the molecular reasons for the strengthened selectivity in Polybia-MP1 and its absence in Mastoparan-X, MD simulations were carried out. Both peptides presented amphipathic alpha-helical structures, as previously observed in Circular Dichroism spectra, with important differences in the extension and stability of the helix; their backbone solvation analysis also indicate a different profile, suggesting that the selectivity of Polybia-MP1 is a consequence of the distribution of the charged and polar residues along the peptide helix, and on how the solvent molecules orient themselves according to these electrostatic interactions. We suggest that the lack of hemolytic activity of Polybia-MP1 is due to the presence and position of Asp residues that enable the equilibrium of electrostatic interactions and favor the preference for the more hydrophilic environment.

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide`s primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.

Theoretical Analysis of Aggregation in Block-Copolymer Films: The Optical Signature

We focus this work on the theoretical investigation of the block-copolymer poly [oxyoctyleneoxy-(2,6-dimethoxy-1,4phenylene-1,2-ethinylene-phenanthrene-2,4diyl) named as LaPPS19, recently proposed for optoelectronic applications. We used for that a variety of methods, from molecular mechanics to quantum semiempirical techniques (AMI, ZINDO/S-CIS). Our results show that as expected isolated LaPPS19 chains present relevant electron localization over the phenanthrene group. We found, however, that LaPPS19 could assemble in a pi-stacked form, leading to impressive interchain interaction; the stacking induces electronic delocalization between neighbor chains and introduces new states below the phenanthrene-related absorption; these results allowed us to associate the red-shift of the absorption edge, seen in the experimental results, to spontaneous pi-stack aggregation of the chains. (C) 2009 Wiley Periodicals, Inc. Int J Quantum Chem 110: 885-892, 2010

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Selective Wrapping and Supramolecular Structures of Polyfluorene-Carbon Nanotube Hybrids

We report on the photophysical properties of single-walled carbon nanotube (SWNT) suspensions In toluene solutions of poly[9,9-dioctylfluorenyl-2,7-diyl](PFO). Steady-state and time-resolved photoluminescence spectroscopy in the near-infrared and visible spectral regions are used to study the interaction of the dispersed SWNTs with the wrapped polymer. Molecular dynamics simulations of the PFO-SWNT hybrids in toluene were carried out to evaluate the energetics of different wrapping geometries. The simulated fluorescence spectra in the visible region were obtained by the quantum chemical ZINDO-CI method, using a sampling of structures obtained from the dynamics trajectories. The tested schemes consider polymer chains aligned along the nanotube axis, where chirality has a minimal effect, or forming helical structures, where a preference for high chiral angles is evidenced. Moreover, toluene affects the polymer structure favoring the helical conformation. Simulations show that the most stable hybrid system is the PFO-wrapped (8,6) nanotube, in agreement with the experimentally observed selectivity.

Differential effects of TR ligands on hormone dissociation rates: Evidence for multiple ligand entry/exit pathways

Some nuclear receptor (NR) ligands promote dissociation of radiolabeled bound hormone from the buried ligand binding cavity (LBC) more rapidly than excess unlabeled hormone itself This result was interpreted to mean that challenger ligands bind allosteric sites on the LBD to induce hormone dissociation, and recent findings indicate that ligands bind weakly to multiple sites on the LBD surface. Here we show, that a large fraction of thyroid hormone receptor (TR) ligands promote rapid dissociation (T(1/2) < 2 h) of , radiolabeled T(3) vs. T(3) (T(1/2), approximate to 5-7 h). We cannot discern relationships between this effect and ligand size, activity or affinity for TR beta. One ligand, GC-24, binds the TR LBC and (weakly) to the TR beta-LBD surface that mediates dimer/heterodimer interaction, but we cannot link this interaction to rapid T(3) dissociation. Instead, several lines of evidence suggest that the challenger ligand must interact with the buried LBC to promote rapid T(3) release. Since previous molecular dynamics simulations suggest that TR ligands leave the LBC by several routes, we propose that a subset of challenger ligands binds and stabilizes a partially unfolded intermediate state of TR that arises during T(3) release and that this effect enhances hormone dissociation. (C) 2009 Elsevier Ltd. All rights reserved.

Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

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)

Structural modeling of high-affinity thyroid receptor-ligand complexes

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors` ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.