High-throughput sequence (HTS) analysis of combinatorial selection populations accelerates lead discovery and optimization and offers dynamic insight into selection processes. An underlying principle is that selection enriches high-fitness sequences as a fraction of the population, whereas low-fitness sequences are depleted. HTS analysis readily provides the requisite numerical information by tracking the evolutionary trajectory of individual sequences in response to selection pressures. Unlike genomic data, for which a number of software solutions exist, user-friendly tools are not readily available for the combinatorial selections field, leading many users to create custom software. FASTAptamer was designed to address the sequence-level analysis needs of the field. The open source FASTAptamer toolkit counts, normalizes and ranks read counts in a FASTQ file, compares populations for sequence distribution, generates clusters of sequence families, calculates fold-enrichment of sequences throughout the course of a selection and searches for degenerate sequence motifs. While originally designed for aptamer selections, FASTAptamer can be applied to any selection strategy that can utilize next-generation DNA sequencing, such as ribozyme or deoxyribozyme selections, in vivo mutagenesis and various surface display technologies (peptide, antibody fragment, mRNA, etc.). FASTAptamer software, sample data and a user's guide are available for download at http://burkelab.missouri.edu/fastaptamer.html.

No abstract

Chlorophyll (Chl) is often viewed as having preceded bacteriochlorophyll (BChl) as the primary photoreceptor pigment in early photosynthetic systems because synthesis of Chl requires one fewer enzymatic reduction than does synthesis of BChl. We have conducted statistical DNA sequence analyses of the two reductases involved in Chl and BChl synthesis, protochlorophyllide reductase and chlorin reductase. Both are three-subunit enzymes in which each subunit from one reductase shares significant amino acid identity with a subunit of the other, indicating that the two enzymes are derived from a common three-subunit ancestral reductase. The "chlorophyll iron protein" subunits, encoded by the bchL and bchX genes in the purple bacterium Rhodobacter capsulatus, also share amino acid sequence identity with the nitrogenase iron protein, encoded by niff. When nitrogenase iron proteins are used as outgroups, the chlorophyll iron protein tree is rooted on the chlorin reductase lineage. This rooting suggests that the last common ancestor of all extant photosynthetic eubacteria contained BChl, not Chl, in its reaction center, and implies that Chl-containing reaction centers were a late invention unique to the cyanobacteria/chloroplast lineage.

In pre-steady-state, fast-quench kinetic analysis, the tertiary-stabilized hammerhead ribozyme ''RzB'' cleaves its substrate RNA with maximal measured k obs values of ; ;3000 min À1 in 1 mM Mn 2+ and ; ;780 min À1 in 1 mM Mg 2+ at 37°C (pH 7.4). Apparent pKa for the catalytic general base is ; ;7.8-8.5, independent of the corresponding metal hydrate pKa, suggesting potential involvement of a nucleobase as general base as suggested previously from nucleobase substitution studies. The pH-rate profile is bell-shaped for Cd 2+ , for which the general catalytic acid has a pKa of 7.3 6 0.1. Simulations of the pH-rate relation suggest a pKa for the general catalytic acid to be ; ;9.5 in Mn 2+ and >9.5 in Mg 2+ . The acid pKa's follow the trend in the pKa of the hydrated metal ions but are displaced by ; ;1-2 pH units in the presence of Cd 2+ and Mn 2+ . One possible explanation for this trend is direct metal ion coordination with a nucleobase, which then acts as general acid.

Systematic evolution of ligands through exponential enrichment (SELEX) is a well-established method for generating nucleic acid populations that are enriched for specified functions. High-throughput sequencing (HTS) enhances the power of comparative sequence analysis to reveal details of how RNAs within these populations recognize their targets. We used HTS analysis to evaluate RNA populations selected to bind type I human immunodeficiency virus reverse transcriptase (RT). The populations are enriched in RNAs of independent lineages that converge on shared motifs and in clusters of RNAs with nearly identical sequences that share common ancestry. Both of these features informed inferences of the secondary structures of enriched RNAs, their minimal structural requirements and their stabilities in RT-aptamer complexes. Monitoring population dynamics in response to increasing selection pressure revealed RNA inhibitors of RT that are more potent than the previously identified pseudoknots. Improved potency was observed for inhibition of both purified RT in enzymatic assays and viral replication in cell-based assays. Structural and functional details of converged motifs that are obscured by simple consensus descriptions are also revealed by the HTS analysis. The approach presented here can readily be generalized for the efficient and systematic post-SELEX development of aptamers for down-stream applications.