The FunFOLD2 server for the prediction of protein-ligand interactions

The FunFOLD2 server is a new independent server that integrates our novel protein–ligand binding site and quality assessment protocols for the prediction of protein function (FN) from sequence via structure. Our guiding principles were, first, to provide a simple unified resource to make our function prediction software easily accessible to all via a simple web interface and, second, to produce integrated output for predictions that can be easily interpreted. The server provides a clean web interface so that results can be viewed on a single page and interpreted by non-experts at a glance. The output for the prediction is an image of the top predicted tertiary structure annotated to indicate putative ligand-binding site residues. The results page also includes a list of the most likely binding site residues and the types of predicted ligands and their frequencies in similar structures. The protein–ligand interactions can also be interactively visualized in 3D using the Jmol plug-in. The raw machine readable data are provided for developers, which comply with the Critical Assessment of Techniques for Protein Structure Prediction data standards for FN predictions. The FunFOLD2 webserver is freely available to all at the following web site: http://www.reading.ac.uk/bioinf/FunFOLD/FunFOLD_form_2_0.html.

Melatonin Protects CD4(+) T Cells from Activation-Induced Cell Death by Blocking NFAT-Mediated CD95 Ligand Upregulation

Over the past 20 y, the hormone melatonin was found to be produced in extrapineal sites, including cells of the immune system. Despite the increasing data regarding the biological effects of melatonin on the regulation of the immune system, the effect of this molecule on T cell survival remains largely unknown. Activation-induced cell death plays a critical role in the maintenance of the homeostasis of the immune system by eliminating self-reactive or chronically stimulated T cells. Because activated T cells not only synthesize melatonin but also respond to it, we investigated whether melatonin could modulate activation-induced cell death. We found that melatonin protects human and murine CD4(+) T cells from apoptosis by inhibiting CD95 ligand mRNA and protein upregulation in response to TCR/CD3 stimulation. This inhibition is a result of the interference with calmodulin/calcineurin activation of NFAT that prevents the translocation of NFAT to the nucleus. Accordingly, melatonin has no effect on T cells transfected with a constitutively active form of NFAT capable of migrating to the nucleus and transactivating target genes in the absence of calcineurin activity. Our results revealed a novel biochemical pathway that regulates the expression of CD95 ligand and potentially other downstream targets of NFAT activation. The Journal of Immunology, 2010, 184: 3487-3494.

Crystallographic structure and substrate-binding interactions of the molybdate-binding protein of the phytopathogen Xanthomonas axonopodis pv. citri

In Xanthomonas axonopodis pv. citri (Xac or X citri), the modA gene codes for a periplasmic protein (ModA) that is capable of binding molybdate and tungstate as part of the ABC-type transporter required for the uptake of micronutrients. In this study, we report the crystallographic structure of the Xac ModA protein with bound molybdate. The Xac ModA structure is similar to orthologs with known three-dimensional structures and consists of two nearly symmetrical domains separated by a hinge region where the oxyanion-binding site lies. Phylogenetic analysis of different ModA orthologs based on sequence alignments revealed three groups of molybdate-binding proteins: bacterial phytopathogens, enterobacteria and soil bacteria. Even though the ModA orthologs are segregated into different groups, the ligand-binding hydrogen bonds are mostly conserved, except for Archaeglobus fulgidus ModA. A detailed discussion of hydrophobic interactions in the active site is presented and two new residues, Ala(38) and Ser(151), are shown to be part of the ligand-binding pocket. (c) 2007 Elsevier B.V All rights reserved.

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.

Ligand induced interaction of thyroid hormone receptor beta with its coregulators

Thyroid hormones exert most of their physiological effects through two thyroid hormone receptor (TR) subtypes, TR alpha and TR beta, which associate with many transcriptional coregulators to mediate activation or repression of target genes. The search for selective TR beta ligands has been stimulated by the finding that several pharmacological actions mediated by TR beta might be beneficial in medical conditions such as obesity, hypercholesterolemia and diabetes. Here, we present a new methodology which employs surface plasmon resonance to investigate the interactions between TR beta ligand binding domain (LBD) complexes and peptides derived from the nuclear receptor interaction motifs of two of its coregulators, SRC2 and DAX1. The effect of several TR beta ligands, including the TR beta selective agonist GC-I and the TR beta selective antagonist NH-3, were investigated. We also determined the kinetic rate constants for the interaction of TR beta-T3 with both coregulators, and accessed the thermodynamic parameters for the interaction with DAX1. Our findings Suggest that flexibility plays an important role in the interaction between the receptor and its coregulators. and point out important aspects of experimental design that should be addressed when using TR beta LBD and its agonists. Furthermore, the methodology described here may be useful for the identification of new TR beta ligands. (C) 2008 Elsevier Ltd. All rights reserved.

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.

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.

Structural requirement for PPAR gamma binding revealed by a meta analysis of holo-crystal structures

PPAR gamma is a ligand regulated transcriptional factor that modulates the transcription of several genes involved in fat and sugar metabolism. Due to its easy bacterial expression and crystallization, several crystal structures of holo-PPAR gamma have been reported and deposited in the Protein Data Bank. Here, we investigated the three-dimensional electrostatic properties of 55 PPAR gamma ligands and used this information for clustering them through principal component analysis. We found out that, according to their electrostatic potential, these ligands can be separated in three groups, with different binding features. We also observed that non-selective and selective ligands show different 3D electrostatic properties and are separated in different clusters. The relevance of this analysis for the development of new binders is discussed. (C) 2010 Elsevier Masson SAS. All rights reserved.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

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.

Crystallographic studies on the binding of isonicotinyl-NAD adduct to wild-type and isoniazid resistant 2-trans-enoyl-ACP (CoA) reductase from Mycobacterium tuberculosis

The resumption of tuberculosis led to an increased need to understand the molecular mechanisms of drug action and drug resistance, which should provide significant insight into the development of newer compounds. Isoniazid (INH), the most prescribed drug to treat TB, inhibits an NADH-dependent enoyl-acyl carrier protein reductase (InhA) that provides precursors of mycolic acids, which are components of the mycobacterial cell wall. InhA is the major target of the mode of action of isoniazid. INH is a pro-drug that needs activation to form the inhibitory INH-NAD adduct. Missense mutations in the inhA structural gene have been identified in clinical isolates of Mycobacterium tuberculosis resistant to INH. To understand the mechanism of resistance to INH, we have solved the structure of two InhA mutants (121V and S94A), identified in INH-resistant clinical isolates, and compare them to INH-sensitive WT InhA structure in complex with the INH-NAD adduct. We also solved the structure of unliganded INH-resistant S94A protein, which is the first report on apo form of InhA. The salient features of these structures are discussed and should provide structural information to improve our understanding of the mechanism of action of, and resistance to, INH in M. tuberculosis. The unliganded structure of InhA allows identification of conformational changes upon ligand binding and should help structure-based drug design of more potent antimycobacterial agents. (c) 2007 Elsevier B.V. All rights reserved.

Water regulation of actinomycin-D binding to DNA: the interplay among drug affinity, DNA long-range conformation, and hydration

Actiaomycin-D (actD) binds to natural DNA at two different classes of binding sites, weak and strong. The affinity for these sites is highly dependent on DNA se(sequence and solution conditions, and the interaction appears to be purely entropic driven Although the entropic character of this reaction has been attributed to the release of water molecules upon drug to DNA complex formation, the mechanism by which hydration regulates actD binding and discrimination between different classes of binding sites on natural DNA is still unknown. In this work, we investigate the role of hydration on this reaction using the osmotic stress method. We skew that the decrease of solution water activity, due to the addition of sucrose, glycerol ethylene glycol, and betaine, favors drug binding to the strong binding sites on DNA by increasing both the apparent binding affinity Delta G, and the number of DNA base pairs apparently occupied by the bound drug n(bp/actD). These binding parameters vary linearly with the logarithm of the molar fraction of water in solution log(X-w), which indicates the contribution of water binding to the energetic of the reaction. It is demonstrated that the hydration change measured upon binding increases proportionally to the apparent size of the binding site n(bp/uctD). This indicates that n(bp/actD) measured from the Scatchard plod is a measure of the size of the DNA molecule changing conformation due to ligand binding. We also find that the contribution of DNA deformation, gauged by n(bp/act) to the total free energy of binding Delta G, is given by Delta G = Delta G(local) + n(bp/actD) x delta G(DNA), where Delta G(local), = -8020 +/- 51 cal/mol of actD bound and delta G(DNa) = -24.1 +/- 1.7cal/mol of base pair at 25 degrees C. We interpret Delta G(local), as the energetic contribution due to the direct interactions of actD with the actual tetranucleotide binding site, and it n(bp/actB) X delta G(DNA) as that due to change inconformation, induced by binding, of it n(bp/actD) DNA base pairs flanking the local site. This interpretation is supported by the agreement found between the value of delta G(DNA) and the torsional free energy change measured independently. We conclude suggesting an allosteric model for ligand binding to DNA, such that the increase in binding affinity is achieved by increasing the relaxation of the unfavorable free energy of binding storage at the local site through a larger number of DNA base pairs. The new aspect on this model is that the size of the complex is not fixed but determined by solutions conditions, such as water activity, which modulate the energetic barrier to change helix conformation. These results may suggest that long-range allosteric transitions of duplex DNA are involved in the inhibition of RNA synthesis by actD, and more generally, in the regulation of transcription. (C) 2000 John Wiley & Sons, Inc.

Novel aryl beta-aminocarbonyl derivatives as inhibitors of Trypanosoma cruzi trypanothione reductase: binding mode revised by docking and GRIND2-based 3D-QSAR procedures

Trypanothione reductase has long been investigated as a promising target for chemotherapeutic intervention in Chagas disease, since it is an enzyme of a unique metabolic pathway that is exclusively present in the pathogen but not in the human host, which has the analog Glutathione reductase. In spite of the present data-set includes a small number of compounds, a combined use of flexible docking, pharmacophore perception, ligand binding site prediction, and Grid-Independent Descriptors GRIND2-based 3D-Quantitative Structure-Activity Relationships (QSAR) procedures allowed us to rationalize the different biological activities of a series of 11 aryl beta-aminocarbonyl derivatives, which are inhibitors of Trypanosoma cruzi trypanothione reductase (TcTR). Three QSAR models were built and validated using different alignments, which are based on docking with the TcTR crystal structure, pharmacophore, and molecular interaction fields. The high statistical significance of the models thus obtained assures the robustness of this second generation of GRIND descriptors here used, which were able to detect the most important residues of such enzyme for binding the aryl beta-aminocarbonyl derivatives, besides to rationalize distances among them. Finally, a revised binding mode has been proposed for our inhibitors and independently supported by the different methodologies here used, allowing further optimization of the lead compounds with such combined structure- and ligand-based approaches in the fight against the Chagas disease.

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.