Actions of Agonists, Fipronil and Ivermectin on the Predominant In Vivo Splice and Edit Variant (RDLbd, I/V) of the Drosophila GABA Receptor Expressed in Xenopus laevis Oocytes

Ionotropic GABA receptors are the targets for several classes of insecticides. One of the most widely-studied insect GABA receptors is RDL (resistance to dieldrin), originally isolated from Drosophila melanogaster. RDL undergoes alternative splicing and RNA editing, which influence the potency of GABA. Most work has focussed on minority isoforms. Here, we report the first characterisation of the predominant native splice variant and RNA edit, combining functional characterisation with molecular modelling of the agonist-binding region. The relative order of agonist potency is GABA> muscimol> TACA> β-alanine. The I/V edit does not alter the potency of GABA compared to RDLbd. Docking calculations suggest that these agonists bind and activate RDLbdI/V through a similar binding mode. TACA and β-alanine are predicted to bind with lower affinity than GABA, potentially explaining their lower potency, whereas the lower potency of muscimol and isoguvacine cannot be explained structurally from the docking calculations. The A301S (resistance to dieldrin) mutation reduced the potency of antagonists picrotoxin, fipronil and pyrafluprole but the I/V edit had no measurable effect. Ivermectin suppressed responses to GABA of RDLbdI/V, RDLbd and RDLbdI/VA301S. The dieldrin resistant variant also showed reduced sensitivity to Ivermectin. This study of a highly abundant insect GABA receptor isoform will help the design of new insecticides

A model for sequential evolution of ligands by exponential enrichment (SELEX) data

A Systematic Evolution of Ligands by EXponential enrichment (SELEX) experiment begins in round one with a random pool of oligonucleotides in equilibrium solution with a target. Over a few rounds, oligonucleotides having a high affinity for the target are selected. Data from a high throughput SELEX experiment consists of lists of thousands of oligonucleotides sampled after each round. Thus far, SELEX experiments have been very good at suggesting the highest affinity oligonucleotide, but modeling lower affinity recognition site variants has been difficult. Furthermore, an alignment step has always been used prior to analyzing SELEX data. We present a novel model, based on a biochemical parametrization of SELEX, which allows us to use data from all rounds to estimate the affinities of the oligonucleotides. Most notably, our model also aligns the oligonucleotides. We use our model to analyze a SELEX experiment containing double stranded DNA oligonucleotides and the transcription factor Bicoid as the target. Our SELEX model outperformed other published methods for predicting putative binding sites for Bicoid as indicated by the results of an in-vivo ChIP-chip experiment.Comment: Published in at http://dx.doi.org/10.1214/12-AOAS537 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Large-Scale Turnover of Functional Transcription Factor Binding Sites in Drosophila

The gain and loss of functional transcription factor binding sites has been proposed as a major source of evolutionary change in cis-regulatory DNA and gene expression. We have developed an evolutionary model to study binding-site turnover that uses multiple sequence alignments to assess the evolutionary constraint on individual binding sites, and to map gain and loss events along a phylogenetic tree. We apply this model to study the evolutionary dynamics of binding sites of the Drosophila melanogaster transcription factor Zeste, using genome-wide in vivo (ChIP–chip) binding data to identify functional Zeste binding sites, and the genome sequences of D. melanogaster, D. simulans, D. erecta, and D. yakuba to study their evolution. We estimate that more than 5% of functional Zeste binding sites in D. melanogaster were gained along the D. melanogaster lineage or lost along one of the other lineages. We find that Zeste-bound regions have a reduced rate of binding-site loss and an increased rate of binding-site gain relative to flanking sequences. Finally, we show that binding-site gains and losses are asymmetrically distributed with respect to D. melanogaster, consistent with lineage-specific acquisition and loss of Zeste-responsive regulatory elements

Crystal structures of a GABAA-receptor chimera reveal new endogenous neurosteroid-binding sites.

γ-Aminobutyric acid receptors (GABAARs) are vital for controlling excitability in the brain. This is emphasized by the numerous neuropsychiatric disorders that result from receptor dysfunction. A critical component of most native GABAARs is the α subunit. Its transmembrane domain is the target for many modulators, including endogenous brain neurosteroids that impact anxiety, stress and depression, and for therapeutic drugs, such as general anesthetics. Understanding the basis for the modulation of GABAAR function requires high-resolution structures. Here we present the first atomic structures of a GABAAR chimera at 2.8-Å resolution, including those bound with potentiating and inhibitory neurosteroids. These structures define new allosteric binding sites for these modulators that are associated with the α-subunit transmembrane domain. Our findings will enable the exploitation of neurosteroids for therapeutic drug design to regulate GABAARs in neurological disorders

JGromacs: A Java Package for Analyzing Protein Simulations

UNLABELLED: In this paper, we introduce JGromacs, a Java API (Application Programming Interface) that facilitates the development of cross-platform data analysis applications for Molecular Dynamics (MD) simulations. The API supports parsing and writing file formats applied by GROMACS (GROningen MAchine for Chemical Simulations), one of the most widely used MD simulation packages. JGromacs builds on the strengths of object-oriented programming in Java by providing a multilevel object-oriented representation of simulation data to integrate and interconvert sequence, structure, and dynamics information. The easy-to-learn, easy-to-use, and easy-to-extend framework is intended to simplify and accelerate the implementation and development of complex data analysis algorithms. Furthermore, a basic analysis toolkit is included in the package. The programmer is also provided with simple tools (e.g., XML-based configuration) to create applications with a user interface resembling the command-line interface of GROMACS applications. AVAILABILITY: JGromacs and detailed documentation is freely available from http://sbcb.bioch.ox.ac.uk/jgromacs under a GPLv3 license

Quantitative principles of cis-translational control by general mRNA sequence features in eukaryotes.

BackgroundGeneral translational cis-elements are present in the mRNAs of all genes and affect the recruitment, assembly, and progress of preinitiation complexes and the ribosome under many physiological states. These elements include mRNA folding, upstream open reading frames, specific nucleotides flanking the initiating AUG codon, protein coding sequence length, and codon usage. The quantitative contributions of these sequence features and how and why they coordinate to control translation rates are not well understood.ResultsHere, we show that these sequence features specify 42-81% of the variance in translation rates in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Arabidopsis thaliana, Mus musculus, and Homo sapiens. We establish that control by RNA secondary structure is chiefly mediated by highly folded 25-60 nucleotide segments within mRNA 5' regions, that changes in tri-nucleotide frequencies between highly and poorly translated 5' regions are correlated between all species, and that control by distinct biochemical processes is extensively correlated as is regulation by a single process acting in different parts of the same mRNA.ConclusionsOur work shows that general features control a much larger fraction of the variance in translation rates than previously realized. We provide a more detailed and accurate understanding of the aspects of RNA structure that directs translation in diverse eukaryotes. In addition, we note that the strongly correlated regulation between and within cis-control features will cause more even densities of translational complexes along each mRNA and therefore more efficient use of the translation machinery by the cell

Survey of large protein complexes D. vulgaris reveals great structural diversity

An unbiased survey has been made of the stable, most abundant multi-protein complexes in Desulfovibrio vulgaris Hildenborough (DvH) that are larger than Mr {approx} 400 k. The quaternary structures for 8 of the 16 complexes purified during this work were determined by single-particle reconstruction of negatively stained specimens, a success rate {approx}10 times greater than that of previous 'proteomic' screens. In addition, the subunit compositions and stoichiometries of the remaining complexes were determined by biochemical methods. Our data show that the structures of only two of these large complexes, out of the 13 in this set that have recognizable functions, can be modeled with confidence based on the structures of known homologs. These results indicate that there is significantly greater variability in the way that homologous prokaryotic macromolecular complexes are assembled than has generally been appreciated. As a consequence, we suggest that relying solely on previously determined quaternary structures for homologous proteins may not be sufficient to properly understand their role in another cell of interest

Subcellular Localization of Proteins in the Anaerobic Sulfate Reducer Desulfovibrio vulgaris via SNAP-tag Labeling and Photoconversion

Systems Biology studies the temporal and spatial 3D distribution of macromolecular complexes with the aim that such knowledge will allow more accurate modeling of biological function and will allow mathematical prediction of cellular behavior. However, in order to accomplish accurate modeling precise knowledge of spatial 3D organization and distribution inside cells is necessary. And while a number of macromolecular complexes may be identified by its 3D structure and molecular characteristics alone, the overwhelming number of proteins will need to be localized using a reporter tag. GFP and its derivatives (XFPs) have been traditionally employed for subcelllar localization using photoconversion approaches, but this approach cannot be taken for obligate anaerobic bacteria, where the intolerance towards oxygen prevents XFP approaches. As part of the GTL-funded PCAP project (now ENIGMA) genetic tools have been developed for the anaerobe sulfate reducer Desulfovibrio vulgaris that allow the high-throughput generation of tagged-protein mutant strains, with a focus on the commercially available SNAP-tag cell system (New England Biolabs, Ipswich, MA), which is based on a modified O6-alkylguanine-DNA alkyltransferase (AGT) tag, that has a dead-end reaction with a modified O6-benzylguanine (BG) derivative and has been shown to function under anaerobic conditions. After initial challenges with respect to variability, robustness and specificity of the labeling signal we have optimized the labeling. Over the last year, as a result of the optimized labeling protocol, we now obtain robust labeling of 20 out of 31 SNAP strains. Labeling for 13 strains were confirmed at least five times. We have also successfully performed photoconversion on 5 of these 13 strains, with distinct labeling patterns for different strains. For example, DsrC robustly localizes to the periplasmic portion of the inner membrane, where as a DNA-binding protein localizes to the center of the cell, where the chromosome is located. Two other proteins - Thiosulfate reductase and ATP binding protein were found to be cytoplasmically distributed, whereas a molybdenum transporter was found to locate to the cell periphery. We judge labeling outcome by (1) SDS gel electrophoresis, followed by direct fluorescence imaging of the gel to address specificity of labeling/confirm expected molecular weight, and subsequent Coomassie analysis to ensure comparable protein levels (2) fluorescence intensity of culture by plate reader for statistical sampling (after adjustment for respective cell numbers) and (3) fluorescence microscopy for addressing cell-to-cell signal variation and potential localization patterns. All three assays were usually found to be consistent with one another. While we have been able to improve the efficacy of photoconversion by drastically reducing (eliminating) non-specific binding with our altered labeling protocol, we are currently working on reducing non-specific photoconversion reaction arising occasionally in non-labeled cells. In addition, we have confirmed the presence of SNAP tagged constructs in three recently cloned E.coli strains under promotor control, and are in the process of utilizing them for evaluating the sensitivity of the photoconversion protocol. Fluorescent Activated Cell Sorting was successfully applied to labeled E.coli cells containing SNAP tagged AtpA protein. Different batches of sorted cells, representing low and high labeling intensity, were re-grown and re-labeled and displayed a labeling efficiency similar to the starter culture, supporting the notion that cell-to-cell differences in labeling reflect difference in protein expression, rather then genetic differences