Assessing computational tools for the discovery of transcription factor binding sites

The prediction of regulatory elements is a problem where computational methods offer great hope. Over the past few years, numerous tools have become available for this task. The purpose of the current assessment is twofold: to provide some guidance to users regarding the accuracy of currently available tools in various settings, and to provide a benchmark of data sets for assessing future tools.

Herbicide-binding Sites Revealed in the Structure of Plant Acetohydroxyacid Synthase

The sulfonylureas and imidazolinones are potent commercial herbicide families. They are among the most popular choices for farmers worldwide, because they are nontoxic to animals and highly selective. These herbicides inhibit branched-chain amino acid biosynthesis in plants by targeting acetohydroxyacid synthase (AHAS, EC 2.2.1.6). This report describes the 3D structure of Arabidopsis thaliana AHAS in complex with five sulfonylureas (to 2.5 angstrom resolution) and with the imidazolinone, imazaquin (IQ; 2.8 angstrom). Neither class of molecule has a structure that mimics the substrates for the enzyme, but both inhibit by blocking a channel through which access to the active site is gained. The sulfonylureas approach within 5 angstrom of the catalytic center, which is the C2 atom of the cofactor thiamin diphosphate, whereas IQ is at least 7 angstrom from this atom. Ten of the amino acid residues that bind the sulfonylureas also bind IQ. Six additional residues interact only with the sulfonylureas, whereas there are two residues that bind IQ but not the sulfonylureas. Thus, the two classes of inhibitor occupy partially overlapping sites but adopt different modes of binding. The increasing emergence of resistant weeds due to the appearance of mutations that interfere with the inhibition of AHAS is now a worldwide problem. The structures described here provide a rational molecular basis for understanding these mutations, thus allowing more sophisticated AHAS inhibitors to be developed. There is no previously described structure for any plant protein in complex with a commercial herbicide.

Sulfadoxine resistance in Plasmodium vivax is associated with a specific amino acid in dihydropteroate synthase at the putative sulfadoxine-binding site

Sulfadoxine is predominantly used in combination with pyrimethamine, commonly known as Fansidar, for the treatment of Plasmodium falciparum. This combination is usually less effective against Plasmodium vivax, probably due to the innate refractoriness of parasites to the sulfadoxine component. To investigate this mechanism of resistance by P. vivax to sulfadoxine, we cloned and sequenced the P. vivax dhps (pvdhps) gene. The protein sequence was determined, and three-dimensional homology models of dihydropteroate synthase (DHPS) from P. vivax as well as P. falciparum were created. The docking of sulfadoxine to the two DHPS models allowed us to compare contact residues in the putative sulfadoxine-binding site in both species. The predicted sulfadoxine-binding sites between the species differ by one residue, V585 in P. vivax, equivalent to A613 in P. falciparum. V585 in P. vivax is predicted by energy minimization to cause a reduction in binding of sulfadoxine to DHPS in P. vivax compared to P. falciparum. Sequencing dhps genes from a limited set of geographically different P. vivax isolates revealed that V585 was present in all of the samples, suggesting that V585 may be responsible for innate resistance of P. vivax to sulfadoxine. Additionally, amino acid mutations were observed in some P. vivax isolates in positions known to cause resistance in P. falciparum, suggesting that, as in P. falciparum, these mutations are responsible for acquired increases in resistance of P. vivax to sulfadoxine.

Determination of alpha-conotoxin binding modes on neuronal nicotinic acetylcholine receptors

alpha-Conotoxins, from cone snails, and alpha-neurotoxins, from snakes, are competitive inhibitors of nicotinic acetylcholine receptors (nAChRs) that have overlapping binding sites in the ACh binding pocket. These disulphide-rich peptides are used extensively as tools to localize and pharmacologically characterize specific nAChRs subtypes. Recently, a homology model based on the high-resolution structure of an ACh binding protein (AChBP) allowed the three-fingered alpha-neurotoxins to be docked onto the alpha7 nAChR. To investigate if alpha-conotoxins interact with the nAChR in a similar manner, we built homology models of human alpha7 and alpha3beta2 nAChRs, and performed docking simulations of alpha-conotoxins ImI, PnIB, PnIA and MII using the program GOLD. Docking revealed that alpha-conotoxins have a different mode of interaction compared with alpha-neurotoxins, with surprisingly few nAChR residues in common between their overlapping binding sites. These docking experiments show that Imi and PnIB bind to the ACh binding pocket via a small cavity located above the beta9/beta10 hairpin of the (+)alpha7 nAChR subunit. Interestingly, PnIB, PnIA and MII were found to bind in a similar location on alpha7 or alpha3beta2 receptors mostly through hydrophobic interactions, while ImI bound further from the ACh binding pocket, mostly through electrostatic interactions. These findings, which distinguish alpha-conotoxin and alpha-neurotoxin binding modes, have implications for the rational design of selective nAChR antagonists. Copyright (C) 2004 John Wiley Sons, Ltd.

Evolutionary turnover of mammalian transcription start sites

Alignments of homologous genomic sequences are widely used to identify functional genetic elements and study their evolution. Most studies tacitly equate homology of functional elements with sequence homology. This assumption is violated by the phenomenon of turnover, in which functionally equivalent elements reside at locations that are nonorthologous at the sequence level. Turnover has been demonstrated previously for transcription-factor-binding sites. Here, we show that transcription start sites of equivalent genes do not always reside at equivalent locations in the human and mouse genomes. We also identify two types of partial turnover, illustrating evolutionary pathways that could lead to complete turnover. These findings suggest that the signals encoding transcription start sites are highly flexible and evolvable, and have cautionary implications for the use of sequence-level conservation to detect gene regulatory elements.

Molecular determinants of ginkgolide binding in the glycine receptor pore

Ginkgolides are potent blockers of the glycine receptor Cl- channel (GlyR) pore. We sought to identify their binding sites by comparing the effects of ginkgolides A, B and C and bilobalide on alpha 1, alpha 2, alpha 1 beta and alpha 2 beta GlyRs. Bilobalide sensitivity was drastically reduced by incorporation of the beta subunit. In contrast, the sensitivities to ginkgolides B and C were enhanced by beta subunit expression. However, ginkgolide A sensitivity was increased in the alpha 2 beta GlyR relative to the alpha 2 GlyR but not in the alpha 1 beta GlyR relative to the alpha 1 GlyR. We hypothesised that the subunit-specific differences were mediated by residue differences at the second transmembrane domain 2' and 6' pore-lining positions. The increased ginkgolide A sensitivity of the alpha 2 beta GlyR was transferred to the alpha 1 beta GlyR by the G2'A (alpha 1 to alpha 2 subunit) substitution. In addition, the alpha 1 subunit T6'F mutation abolished inhibition by all ginkgolides. As the ginkgolides share closely related structures, their molecular interactions with pore-lining residues were amenable to mutant cycle analysis. This identified an interaction between the variable R2 position of the ginkgolides and the 2' residues of both alpha 1 and beta subunits. These findings provide strong evidence for ginkgolides binding at the 2' pore-lining position.