Improved assay-dependent searching of nucleic acid sequence databases

Nucleic acid-based biochemical assays are crucial to modern biology. Key applications, such as detection of bacterial, viral and fungal pathogens, require detailed knowledge of assay sensitivity and specificity to obtain reliable results. Improved methods to predict assay performance are needed for exploiting the exponentially growing amount of DNA sequence data and for reducing the experimental effort required to develop robust detection assays. Toward this goal, we present an algorithm for the calculation of sequence similarity based on DNA thermodynamics. In our approach, search queries consist of one to three oligonucleotide sequences representing either a hybridization probe, a pair of Padlock probes or a pair of PCR primers with an optional TaqMan™ probe (i.e. in silico or ‘virtual’ PCR). Matches are reported if the query and target satisfy both the thermodynamics of the assay (binding at a specified hybridization temperature and/or change in free energy) and the relevant biological constraints (assay sequences binding to the correct target duplex strands in the required orientations). The sensitivity and specificity of our method is evaluated by comparing predicted to known sequence tagged sites in the human genome. Free energy is shown to be a more sensitive and specific match criterion than hybridization temperature

FAD binding, cobinamide binding and active site communication in the corrin reductase (CobR)

Adenosylcobalamin, the coenzyme form of vitamin B12, is one Nature's most complex coenzyme whose de novo biogenesis proceeds along either an anaerobic or aerobic metabolic pathway. The aerobic synthesis involves reduction of the centrally chelated cobalt metal ion of the corrin ring from Co(II) to Co(I) before adenosylation can take place. A corrin reductase (CobR) enzyme has been identified as the likely agent to catalyse this reduction of the metal ion. Herein, we reveal how Brucella melitensis CobR binds its coenzyme FAD (flavin dinucleotide) and we also show that the enzyme can bind a corrin substrate consistent with its role in reduction of the cobalt of the corrin ring. Stopped-flow kinetics and EPR reveal a mechanistic asymmetry in CobR dimer that provides a potential link between the two electron reduction by NADH to the single electron reduction of Co(II) to Co(I)

MICA: desktop software for comprehensive searching of DNA databases

BACKGROUND: Molecular biologists work with DNA databases that often include entire genomes. A common requirement is to search a DNA database to find exact matches for a nondegenerate or partially degenerate query. The software programs available for such purposes are normally designed to run on remote servers, but an appealing alternative is to work with DNA databases stored on local computers. We describe a desktop software program termed MICA (K-Mer Indexing with Compact Arrays) that allows large DNA databases to be searched efficiently using very little memory. RESULTS: MICA rapidly indexes a DNA database. On a Macintosh G5 computer, the complete human genome could be indexed in about 5 minutes. The indexing algorithm recognizes all 15 characters of the DNA alphabet and fully captures the information in any DNA sequence, yet for a typical sequence of length L, the index occupies only about 2L bytes. The index can be searched to return a complete list of exact matches for a nondegenerate or partially degenerate query of any length. A typical search of a long DNA sequence involves reading only a small fraction of the index into memory. As a result, searches are fast even when the available RAM is limited. CONCLUSION: MICA is suitable as a search engine for desktop DNA analysis software

Genomorama: genome visualization and analysis

<p>Abstract</p> <p>Background</p> <p>The ability to visualize genomic features and design experimental assays that can target specific regions of a genome is essential for modern biology. To assist in these tasks, we present Genomorama, a software program for interactively displaying multiple genomes and identifying potential DNA hybridization sites for assay design.</p> <p>Results</p> <p>Useful features of Genomorama include genome search by DNA hybridization (probe binding and PCR amplification), efficient multi-scale display and manipulation of multiple genomes, support for many genome file types and the ability to search for and retrieve data from the National Center for Biotechnology Information (NCBI) Entrez server.</p> <p>Conclusion</p> <p>Genomorama provides an efficient computational platform for visualizing and analyzing multiple genomes.</p

Characterisation of PduS, the pdu Metabolosome Corrin Reductase, and Evidence of Substructural Organisation within the Bacterial Microcompartment

PduS is a corrin reductase and is required for the reactivation of the cobalamin-dependent diol dehydratase. It is one component encoded within the large propanediol utilisation (pdu) operon, which is responsible for the catabolism of 1,2-propanediol within a self-assembled proteinaceous bacterial microcompartment. The enzyme is responsible for the reactivation of the cobalamin coenzyme required by the diol dehydratase. The gene for the cobalamin reductase from Citrobacter freundii (pduS) has been cloned to allow the protein to be overproduced recombinantly in E. coli with an N-terminal His-tag. Purified recombinant PduS is shown to be a flavoprotein with a non-covalently bound FMN that also contains two coupled [4Fe-4S] centres. It is an NADH-dependent flavin reductase that is able to mediate the one-electron reductions of cob(III)alamin to cob(II)alamin and cob(II)alamin to cob(I)alamin. The [4Fe-4S] centres are labile to oxygen and their presence affects the midpoint redox potential of flavin. Evidence is presented that PduS is able to bind cobalamin, which is inconsistent with the view that PduS is merely a flavin reductase. PduS is also shown to interact with one of the shell proteins of the metabolosome, PduT, which is also thought to contain an [Fe-S] cluster. PduS is shown to act as a corrin reductase and its interaction with a shell protein could allow for electron passage out of the bacterial microcompartment

CSAR Benchmark Exercise of 2010: Selection of the Protein–Ligand Complexes

ABSTRACT: A major goal in drug design is the improvement of computational methods for docking and scoring. The Community Structure Activity Resource (CSAR) aims to collect available data from industry and academia which may be used for this purpose (www.csardock.org). Also, CSAR is charged with organizing community-wide exercises based on the collected data. The first of these exercises was aimed to gauge the overall state of docking and scoring, using a large and diverse data set of protein ligand complexes. Participants were asked to calculate the affinity of the complexes as provided and then recalculate with changes which may improve their specific method. This first data set was selected from existing PDB entries which had binding data (Kd or Ki) in Binding MOAD, augmented with entries from PDBbind. The final data set contains 343 diverse protein ligand complexes and spans 14 pKd. Sixteen proteins have three or more complexes in the data set, from which a user could start an inspection of congeneric series. Inherent experimental error limits the possible correlation between scores and measured affinity; R 2 is limited to ∼0.9 when fitting to the data set without over parametrizing. R 2 is limited to ∼0.8 when scoring the data set with a method trained on outside data. The details of how the data set was initially selected, and the process by which it matured t