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.

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.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: Implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis. (C) 2010 Elsevier Inc. All rights reserved.

Whydration effects on DNA double helix stability modulates ligand binding to natural DNA in response to changes in water activity

In this work we present evidence that water molecules are actively involved on the control of binding affinity and binding site discrimination of a drug to natural DNA. In a previous study, the effect of water activity (a(w)) on the energetic parameters of actinomycin-D intercalation to natural DNA was determined using the osmotic stress method (39). This earlier study has shown evidence that water molecules act as an allosteric regulator of ligand binding to DNA via the effect of water activity on the long-range stability of the DNA secondary structure. In this work we have carried out DNA circularization experiments using the plasmid pUC18 in the absence of drugs and in the presence of different neutral solutes to evaluate the contribution of water activity to the energetics of DNA helix unwinding. The contribution of water to these independent reactions were made explicit by the description of how the changes in the free energy of ligand binding to DNA and in the free energy associated with DNA helix torsional deformation are linked to a(w) via changes in structural hydration. Taken together, the results of these studies reveal an extensive linkage between ligand binding affinity and site binding discrimination, and long range helix conformational changes and DNA hydration, This is strong evidence that water molecules work as a classical allosteric regulator of ligand binding to the DNA via its contribution to the stability of the double helix secondary structure, suggesting a possible mechanism by which the biochemical machinery of DNA processing takes advantage of the low activity of water into the cellular milieu.

Sequence, evolution and ligand binding properties of mammalian Duffy antigen/receptor for chemokines

The Duffy antigen/receptor for chemokine, DARC, acts as a widely expressed promiscuous chemokine receptor and as the erythrocyte receptor for Plasmodium vivax. To gain insight into the evolution and structure/function relations of DARC, we analyzed the binding of anti-human Fy monoclonal antibodies (mAbs) and human chemokines to red blood cells (RBCs) from 11 nonhuman primates and two nonprimate mammals, and we elucidated the structures of the DARC genes from gorilla, gibbon, baboon, marmoset, tamarin, night monkey and cattle. CXCL-8 and CCL-5 chemokine binding analysis indicated that the promiscuous binding profile characteristic of DARC is conserved across species. Among three mAbs that detected the Fy6 epitope by flow cytometric analysis of human and chimpanzee RBCs, only one reacted with night monkey and squirrel monkey. Only chimpanzee RBCs bound a significant amount of the anti-Fy3 mAb. Fy3 was also poorly detected on RBCs from gorilla, baboon and rhesus monkey, but not from new world monkeys. Alignment of DARC homologous sequences allowed us to construct a phylogenetic tree in which all branchings were in accordance with current knowledge of primate phylogeny. Although DARC was expected to be under strong internal and external selection pressure, in order to maintain chemokine binding and avoid Plasmodium vivax binding, respectively, our present study did not provide arguments in favor of a selection pressure on the extracellular domains involved in ligand specificity. The amino acid variability of DARC-like polypeptides was found to be well correlated with the hydrophylicity indexes, with the highest divergence on the amino-terminal extracellular domain. Analysis of the deduced amino acid sequences highlighted the conservation of some amino acid residues, which should prove to be critical for the structural and functional properties of DARC.

Coupling Between GABA(A)-R Ligand-Binding Activity and Membrane Organization in beta-Cyclodextrin-Treated Synaptosomal Membranes from Bovine Brain Cortex: New Insights from EPR Experiments

Correlations between GABA(A) receptor (GABA(A)-R) activity and molecular organization of synaptosomal membranes (SM) were studied along the protocol for cholesterol (Cho) extraction with beta-cyclodextrin (beta-CD). The mere pre-incubation (PI) at 37A degrees C accompanying the beta-CD treatment was an underlying source of perturbations increasing [H-3]-FNZ maximal binding (70%) and K (d) (38%), plus a stiffening of SMs' hydrocarbon core region. The latter was inferred from an increased compressibility modulus (K) of SM-derived Langmuir films, a blue-shifted DPH fluorescence emission spectrum and the hysteresis in DPH fluorescence anisotropy (A (DPH)) in SMs submitted to a heating-cooling cycle (4-37-4A degrees C) with A (DPH,heating) < A (DPH,cooling). Compared with PI samples, the beta-CD treatment reduced B (max) by 5% which correlated with a 45%-decrement in the relative Cho content of SM, a decrease in K and in the order parameter in the EPR spectrum of a lipid spin probe labeled at C5 (5-SASL), and significantly increased A (TMA-DPH). PI, but not beta-CD treatment, could affect the binding affinity. EPR spectra of 5-SASL complexes with beta-CD-, SM-partitioned, and free in solution showed that, contrary to what is usually assumed, beta-CD is not completely eliminated from the system through centrifugation washings. It was concluded that beta-CD treatment involves effects of at least three different types of events affecting membrane organization: (a) effect of PI on membrane annealing, (b) effect of residual beta-CD on SM organization, and (c) Cho depletion. Consequently, molecular stiffness increases within the membrane core and decreases near the polar head groups, leading to a net increase in GABA(A)-R density, relative to untreated samples.

A fragment-based approach for ligand binding affinity and selectivity for the liver X receptor beta

Selective modulation of liver X receptor beta (LXR beta) has been recognized as an important approach to prevent or reverse the atherosclerotic process. In the present work, we have developed robust conformation-independent fragment-based quantitative structure-activity and structure-selectivity relationship models for a series of quinolines and cinnolines as potent modulators of the two LXR sub-types. The generated models were then used to predict the potency of an external test set and the predicted values were in good agreement with the experimental results, indicating the potential of the models for untested compounds. The final 2D molecular recognition patterns obtained were integrated to 3D structure-based molecular modeling studies to provide useful insights into the chemical and structural determinants for increased LXR beta binding affinity and selectivity. (C) 2011 Elsevier Inc. All rights reserved.

Structural Insights into Human Peroxisome Proliferator Activated Receptor Delta (PPAR-Delta) Selective Ligand Binding

Peroxisome proliferator activated receptors (PPARs delta, alpha and gamma) are closely related transcription factors that exert distinct effects on fatty acid and glucose metabolism, cardiac disease, inflammatory response and other processes. Several groups developed PPAR subtype specific modulators to trigger desirable effects of particular PPARs without harmful side effects associated with activation of other subtypes. Presently, however, many compounds that bind to one of the PPARs cross-react with others and rational strategies to obtain highly selective PPAR modulators are far from clear. GW0742 is a synthetic ligand that binds PPAR delta more than 300-fold more tightly than PPAR alpha or PPAR gamma but the structural basis of PPAR delta: GW0742 interactions and reasons for strong selectivity are not clear. Here we report the crystal structure of the PPAR delta:GW0742 complex. Comparisons of the PPAR delta:GW0742 complex with published structures of PPARs in complex with alpha and gamma selective agonists and pan agonists suggests that two residues (Val312 and Ile328) in the buried hormone binding pocket play special roles in PPAR delta selective binding and experimental and computational analysis of effects of mutations in these residues confirms this and suggests that bulky substituents that line the PPAR alpha and gamma ligand binding pockets as structural barriers for GW0742 binding. This analysis suggests general strategies for selective PPAR delta ligand design.

The Ras inhibitors caveolin-1 and docking protein 1 activate peroxisome proliferator-activated receptor ? through spatial relocalization at helix 7 of its ligand-binding domain

Peroxisome proliferator-activated receptor ? (PPAR?) is a transcription factor that promotes differentiation and cell survival in the stomach. PPAR? upregulates and interacts with caveolin-1 (Cav1), a scaffold protein of Ras/mitogen-activated protein kinases (MAPKs). The cytoplasmic-to-nuclear localization of PPAR? is altered in gastric cancer (GC) patients, suggesting a so-far-unknown role for Cav1 in spatial regulation of PPAR? signaling. We show here that loss of Cav1 accelerated proliferation of normal stomach and GC cells in vitro and in vivo. Downregulation of Cav1 increased Ras/MAPK-dependent phosphorylation of serine 84 in PPAR? and enhanced nuclear translocation and ligand-independent transcription of PPAR? target genes. In contrast, Cav1 overexpression sequestered PPAR? in the cytosol through interaction of the Cav1 scaffolding domain (CSD) with a conserved hydrophobic motif in helix 7 of PPAR?'s ligand-binding domain. Cav1 cooperated with the endogenous Ras/MAPK inhibitor docking protein 1 (Dok1) to promote the ligand-dependent transcriptional activity of PPAR? and to inhibit cell proliferation. Ligand-activated PPAR? also reduced tumor growth and upregulated the Ras/MAPK inhibitors Cav1 and Dok1 in a murine model of GC. These results suggest a novel mechanism of PPAR? regulation by which Ras/MAPK inhibitors act as scaffold proteins that sequester and sensitize PPAR? to ligands, limiting proliferation of gastric epithelial cells.

Ligand binding induces Cbl-dependent EphB1 receptor degradation through the lysosomal pathway

Eph receptor tyrosine kinases play a critical role in embryonic patterning and angiogenesis. In the adult, they are involved in carcinogenesis and pathological neovascularization. However, the mechanisms underlying their role in tumor formation and metastasis remain to be defined. Here, we demonstrated that stimulation of EphB1 with ephrinB1/Fc led to a marked downregulation of EphB1 protein, a process blocked by the lysosomal inhibitor bafilomycin. Following ephrinB1 stimulation, the ubiquitin ligase Cbl was recruited by EphB1 and then phosphorylated. Both Cbl phosphorylation and EphB1 ubiquitination were blocked by the Src inhibitor PP2. Overexpression of wild-type Cbl, but not of 70Z mutant lacking ligase activity, enhanced EphB1 ubiquitination and degradation. This negative regulation required the tyrosine kinase activity of EphB1 as kinase-dead EphB1-K652R was resistant to Cbl. Glutathione S-transferase binding experiments showed that Cbl bound to EphB1 through its tyrosine kinase-binding domain. In aggregate, we demonstrated that Cbl induces the ubiquitination and lysosomal degradation of activated EphB1, a process requiring EphB1 and Src kinase activity. To our knowledge, this is the first study dissecting the molecular mechanisms leading to EphB1 downregulation, thus paving the way to new means of modulating their angiogenic and tumorigenic properties.

AMPA RECEPTOR ALLOSTERISM: MEASUREMENT OF THE CONFORMATIONAL CHANGES IN THE LIGAND BINDING DOMAIN OF A FUNCTIONAL RECEPTOR

Ionotropic glutamate receptors are important excitatory neurotransmitter receptors in the mammalian central nervous system that have been implicated in a number of neuropathologies such as epilepsy, ischemia, and amyotrophic lateral sclerosis. Glutamate binding to an extracellular ligand binding domain initiates a series of structural changes that leads to the formation of a cation selective transmembrane channel, which consequently closes due to desensitization of the receptor. The crystal structures of the AMPA subtype of the glutamate receptor have been particularly useful in providing initial insight into the conformational changes in the ligand binding domain; however, these structures are limited by crystallographic constraint. To gain a clear picture of how agonist binding is coupled to channel activation and desensitization, it is essential to study changes in the ligand binding domain in a dynamic, physiological state. In this dissertation, a technique called Luminescence Resonance Energy Transfer was used to determine the conformational changes associated with activation and desensitization in a functional AMPA receptor (ÄN*-AMPA) that contains the ligand binding domain and transmembrane segments; ÄN*-AMPA has been modified such that fluorophores can be introduced at specific sites to serve as a readout of cleft closure or to establish intersubunit distances. Previous structural studies of cleft closure of the isolated ligand binding domain in conjunction with functional studies of the full receptor suggest that extent of cleft closure correlates with extent of activation. Here, LRET has been used to show that a similar relationship between cleft closure and activation is observed in the “full length” receptor showing that the isolated ligand binding domain is a good model of the domain in the full length receptor for changes within a subunit. Similar LRET investigations were used to study intersubunit distances specifically to probe conformational changes between subunits within a dimer in the tetrameric receptor. These studies show that the dimer interface is coupled in the open state, and decoupled in the desensitized state, similar to the isolated ligand binding domain crystal structure studies. However, we show that the apo state dimer interface is not pre-formed as in the crystal structure, hence suggesting a mechanism for functional transitions within the receptor based on LRET distances obtained.

RNA binding characteristics and overall topology of the vaccinia virus polyadenylation complex

mRNA 3′ polyadenylation is central to mRNA biogenesis in prokaryotes and eukaryotes, and is implicated in numerous aspects of mRNA metabolism, including efficiency of mRNA export from the nucleus, message stability, and initiation of translation. However, due to the great complexity of the eukaryotic polyadenylation apparatus, the mechanisms of RNA 3 ′ end processing have remained elusive. Although the RNA processing reactions leading to polyadenylated messenger RNA have been studied in many systems, and much progress has been made, a complete understanding of the biochemistry of the poly(A) polymerase enzyme is still lacking. My research uses Vaccinia virus as a model system to gain a better understanding of this complicated polyadenylation process, which consist of RNA binding, catalysis and polymerase translocation. ^ Vaccinia virus replicates in the cytoplasm of its host cell, so it must employ its own poly(A) polymerase (PAP), a heterodimer of two virus encoded proteins, VP55 and VP39. VP55 is the catalytic subunit, adding 30 adenylates to a non-polyadenylated RNA in a rapid processive manner before abruptly changing to a slow, non-processive mode of adenylate addition and dissociating from the RNA. VP39 is the stimulatory subunit. It has no polyadenylation catalytic activity by itself, but when associated with VP55 it facilitates the semi-processive synthesis of tails several hundred adenylates in length. ^ Oligonucleotide selection and competition studies have shown that the heterodimer binds a minimal motif of (rU)2 (N)25 U, the “heterodimer binding motif”, within an oligonucleotide, and its primer selection for polyadenylation is base-type specific. ^ Crosslinking studies using photosensitive uridylate analogs show that within a VP55-VP39-primer ternary complex, VP55 comes into contact with all three required uridylates, while VP39 only contacts the downstream uridylate. Further studies, using a backbone-anchored photosensitive crosslinker show that both PAP subunits are in close proximity to the downstream −10 to −21 region of 50mer model primers containing the heterodimer binding motif. This equal crosslinking to both subunits suggests that the dimerization of VP55 and VP39 creates either a cleft or a channel between the two subunits through which this region of RNA passes. ^ Peptide mapping studies of VP39 covalently crosslinked to the oligonucleotide have identified residue R107 as the amino acid in close proximity to the −10 uridylate. This helps us project a conceptual model onto the known physical surface of this subunit. In the absence of any tertiary structural data for VP55, we have used a series of oligonucleotide selection assays, as well as crosslinking, nucleotide transfer assays, and gel shift assays to gain insight into the requirements for binding, polyadenylation and translocation. Collectively, these data allow us to put together a comprehensive model of the structure and function of the polyadenylation ternary complex consisting of VP39, VP55 and RNA. ^

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.

Structural factors controlling ligand binding to myoglobin: A kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as “kinetic hole burning.” The study of these correlations gives direct insights into the structure–function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Structural changes in hemoglobin during adsorption to solid surfaces: Effects of pH, ionic strength, and ligand binding

We have studied the adsorption of two structurally similar forms of hemoglobin (met-Hb and HbCO) to a hydrophobic self-assembled methyl-terminated thiol monolayer on a gold surface, by using a Quartz Crystal Microbalance (QCM) technique. This technique allows time-resolved simultaneous measurements of changes in frequency (f) (c.f. mass) and energy dissipation (D) (c.f. rigidity/viscoelastic properties) of the QCM during the adsorption process, which makes it possible to investigate the viscoelastic properties of the different protein layers during the adsorption process. Below the isoelectric points of both met-Hb and HbCO, the ΔD vs. Δf graphs displayed two phases with significantly different slopes, which indicates two states of the adsorbed proteins with different visco-elastic properties. The slope of the first phase was smaller than that of the second phase, which indicates that the first phase was associated with binding of a more rigidly attached, presumably denatured protein layer, whereas the second phase was associated with formation of a second layer of more loosely bound proteins. This second layer desorbed, e.g., upon reduction of Fe3+ of adsorbed met-Hb and subsequent binding of carbon monoxide (CO) forming HbCO. Thus, the results suggest that the adsorbed proteins in the second layer were in a native-like state. This information could only be obtained from simultaneous, time-resolved measurements of changes in both D and f, demonstrating that the QCM technique provides unique information about the mechanisms of protein adsorption to solid surfaces.