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.

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.

Computer-aided design of a novel ligand for retinoic acid receptor in cancer chemotherapy

The isotypes of RAR and RXR are retinoic acid and retinoid X acid receptors, respectively, whose ligand-binding domain contains the ligand-dependent activation function, with distinct pharmacological targets for retinoids, involved in the treatment of various cancers and skin diseases. Due to the major challenge which cancer treatment and cure still imposes after many decades to the international scientific community, there is actually considerable interest in new ligands with increased bioactivity. We have focused on the retinoid acid receptor, which is considered an interesting target for drug design. In this work, we carried out density functional geometry optimizations, and different docking procedures. We performed screening in a large database (hundreds of thousands of molecules which we optimized at the AM1 level) yielding a set of potential bioactive ligands. A new ligand was selected and optimized at the B3LYP/6-31G* level. A flexible docking program was used to investigate the interactions between the receptor and the new ligand. The result of this work is compared with several crystallographic ligands of RAR. Our theoretically more bioactive new-ligand indicates stronger and more hydrogen bonds as well as hydrophobic interactions with the receptor. (c) 2005 Wiley Periodicals, Inc.

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.

Singular Value Decomposition and Ligand Binding Analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Three-dimensional electron microscopy of myriapod hemocyanin, planorbid snail acetylcholine-binding protein and cyanobacterial Vipp1

Three-dimensional electron microscopy (3-D EM) provides a framework for the analysis of large protein quaternary structures. The advantage over the generally higher resolving meth- od of X-ray crystallography is the embedding of the proteins in their physiological environ- ment. However, results of the two methods can be combined to obtain superior structural information. In this work, three different protein types – (i) Myriapod hemocyanin, (ii) vesi- cle-inducing protein in plastids 1 (Vipp1) and (iii) acetylcholine-binding protein (AChBP) – were structurally analyzed by 2-D and 3-D EM and, where possible, functionally interpreted.rnMyriapod hemocyanins have been previously shown to be 6x6-meric assemblies that, in case of Scutigera coleoptrata hemocyanin (ScoHc), show two 3x6-mer planes whith a stag- gering angle of approximately 60°. Here, previously observed structural differences between oxy- and deoxy-ScoHc could be substantiated. A 4° rotation between hexamers of two dif- ferent 3x6-mer planes was measured, which originates at the most central inter-hexamer in- terface. Further information about allosteric behaviour in myriapod hemocyanin was gained by analyzing Polydesmus angustus hemocyanin (PanHc), which shows a stable 3x6-mer and divergent histidine patterns in the inter-hexamer interfaces when compared to ScoHc. Both findings would conclusively explain the very different oxygen binding properties of chilopod and diplopod hemocyanin.rnVipp1 is a protein found in cyanobacteria and higher plants which is essential for thyla- koid membrane function and forms highly variable ring-shaped structures. In the course of this study, the first 3-D analysis of Vipp1 was conducted and yielded reconstructions of six differently sized Vipp1 rings from negatively stained images at resolutions between 20 to 30 Å. Furthermore, mutational analyses identified specific N-terminal amino acids that are essential for ring formation. On the basis of these analyses and previously published results, a hypothetical model of the Vipp1 tertiary and quaternary structure was generated.rnAChBP is a water-soluble protein in the hemolymph of mollusks. It is a structural and functional homologue of the ligand-binding domain of nicotinic acetylcholine receptors. For the freshwater snail Biomphalaria glabrata, we previously described two types of AChBP (BgAChBP1 and BgAChBP2). In this work, a 6 Å 3-D reconstruction of native BgAChBP is presented, which shows a dodecahedral assembly that is unprecedented for an AChBP. Single particle analysis of recombinantely expressed BgAChBP types led to preliminary results show- ing a dodecahedral assembly of BgAChBP1 and a dipentameric assembly of BgAChBP2. This indicates divergent biological functions of the two types.

Purification, cDNA structure and biological significance of a single insulin-like growth factor-binding domain protein (SIBD-1) identified in the hemocytes of the spider Cupiennius salei

Cupiennius salei single insulin-like growth factor-binding domain protein (SIBD-1), which exhibits an IGFBP N-terminal domain-like profile, was identified in the hemocytes of the spider C. salei. SIBD-1 was purified by RP-HPLC and the sequence determined by a combination of Edman degradation and 5'-3'- RACE PCR. The peptide (8676.08 Da) is composed of 78 amino acids, contains six intrachain disulphide bridges and carries a modified Thr residue at position 2. SIBD-1 mRNA expression was detected by quantitative real-time PCR mainly in hemocytes, but also in the subesophageal nerve mass and muscle. After infection, the SIBD-1 content in the hemocytes decreases and, simultaneously, the temporal SIBD-1 expression seems to be down-regulated. Two further peptides, SIBD-2 and IGFBP-rP1, also exhibiting IGFBP N-terminal domain variants with unknown functions, were identified on cDNA level in spider hemocytes and venom glands. We conclude that SIBD-1 may play an important role in the immune system of spiders.

Crystal structures of progressive Ca2+ binding states of the Ca2+ sensor Ca2+ binding domain 1 (CBD1) from the CALX Na+/Ca2+ exchanger reveal incremental conformational transitions.

Na(+)/Ca(2+) exchangers (NCX) constitute a major Ca(2+) export system that facilitates the re-establishment of cytosolic Ca(2+) levels in many tissues. Ca(2+) interactions at its Ca(2+) binding domains (CBD1 and CBD2) are essential for the allosteric regulation of Na(+)/Ca(2+) exchange activity. The structure of the Ca(2+)-bound form of CBD1, the primary Ca(2+) sensor from canine NCX1, but not the Ca(2+)-free form, has been reported, although the molecular mechanism of Ca(2+) regulation remains unclear. Here, we report crystal structures for three distinct Ca(2+) binding states of CBD1 from CALX, a Na(+)/Ca(2+) exchanger found in Drosophila sensory neurons. The fully Ca(2+)-bound CALX-CBD1 structure shows that four Ca(2+) atoms bind at identical Ca(2+) binding sites as those found in NCX1 and that the partial Ca(2+) occupancy and apoform structures exhibit progressive conformational transitions, indicating incremental regulation of CALX exchange by successive Ca(2+) binding at CBD1. The structures also predict that the primary Ca(2+) pair plays the main role in triggering functional conformational changes. Confirming this prediction, mutagenesis of Glu(455), which coordinates the primary Ca(2+) pair, produces dramatic reductions of the regulatory Ca(2+) affinity for exchange current, whereas mutagenesis of Glu(520), which coordinates the secondary Ca(2+) pair, has much smaller effects. Furthermore, our structures indicate that Ca(2+) binding only enhances the stability of the Ca(2+) binding site of CBD1 near the hinge region while the overall structure of CBD1 remains largely unaffected, implying that the Ca(2+) regulatory function of CBD1, and possibly that for the entire NCX family, is mediated through domain interactions between CBD1 and the adjacent CBD2 at this hinge.

Cytochrome P-450 reductase FAD binding domain: Cloning, enzymatic activity and flavin binding characterization

The Mixed Function Oxidase System metabolizes a wide range of biochemicals including drugs, pesticides and steroids. Cytochrome P450 reductase is a key enzymatic component of this system, supplying reducing equivalents from NADPH to cytochrome P450. The electrons are shuttled through reductase via two flavin moieties: FAD and FMN. Although the exact mechanism of flavins action is not known, the enzymatic features of reductase greatly depleted of either FMN of FAD have been characterized. Additionally, flavin location within reductase has been proposed by homology and chemical modification studies. This study seeks to extend the flavin depletion analysis in a more controlled system by eliminating the proposed FMN binding domain with recombinant DNA techniques and biochemical analysis. Two P450 reductase cDNA clones containing only the FMN and NADPH binding domain were isolated, expressed and the protein products purified and analysed. This study confirms the proposed FAD binding site, role of FAD in electron shuttling pathway and provides new methods to study the FAD binding domain. ^

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.