Comparative accuracy of methods for protein sequence similarity search

MOTIVATION: Searching a protein sequence database for homologs is a powerful tool for discovering the structure and function of a sequence. Two new methods for searching sequence databases have recently been described: Probabilistic Smith-Waterman (PSW), which is based on Hidden Markov models for a single sequence using a standard scoring matrix, and a new version of BLAST (WU-BLAST2), which uses Sum statistics for gapped alignments. RESULTS: This paper compares and contrasts the effectiveness of these methods with three older methods (Smith- Waterman: SSEARCH, FASTA and BLASTP). The analysis indicates that the new methods are useful, and often offer improved accuracy. These tools are compared using a curated (by Bill Pearson) version of the annotated portion of PIR 39. Three different statistical criteria are utilized: equivalence number, minimum errors and the receiver operating characteristic. For complete-length protein query sequences from large families, PSW's accuracy is superior to that of the other methods, but its accuracy is poor when used with partial-length query sequences. False negatives are twice as common as false positives irrespective of the search methods if a family-specific threshold score that minimizes the total number of errors (i.e. the most favorable threshold score possible) is used. Thus, sensitivity, not selectivity, is the major problem. Among the analyzed methods using default parameters, the best accuracy was obtained from SSEARCH and PSW for complete-length proteins, and the two BLAST programs, plus SSEARCH, for partial-length proteins.      

Large acceleration of alpha-chymotrypsin-catalyzed dipeptide formation by 18-crown-6 in organic solvents

The effects of 18-crown-6 on the synthesis of peptides catalyzed by alpha-chymotrypsin are reported. Lyophilization of the enzyme in the presence of 50 equivalents of 18-crown-6 results in a 425-fold enhanced activity when the reaction between the 2-chloroethylester of N-acetyl-L-phenylalanine and L-phenylalaninamide is carried out in acetonitrile. Addition of crown ether renders the dipeptide synthesis in nonaqueous solvents catalyzed by alpha-chymotrypsin possible on a preparative scale. The acceleration is observed in different solvents and for various peptide precursors. Copyright 1998 John Wiley & Sons, Inc.     

cAMP and in vitro inotropic actions of secretin and VIP in rat papillary muscle.

Secretin and VIP stimulate cardiac adenylyl cyclase activity and exert a positive inotropic action in several mammalian species. This study examined positive inotropic activity and cAMP levels in rat papillary muscle. Isoproterenol and secretin increased contractions by 150+/-31% and 129+/-27%, respectively. VIP increased contraction by 30+/-21% only at 10 microM. Isoproterenol significantly increased cAMP levels by 82%, whereas increases by secretin (58%) and VIP (56%) were not significant. These results are consistent with reports that secretin and VIP stimulate cardiac adenylyl cyclase in the rat, but suggest that cAMP tissue levels cannot totally explain the positive inotropic responses to secretin and VIP.     

Actin filaments elongate from their membrane-associated ends

In limulus sperm an actin filament bundle 55 mum in length extends from the acrosomal vacuole membrane through a canal in the nucleus and then coils in a regular fashion around the base of the nucleus. The bundle expands systematically from 15 filaments near the acrosomal vacuole to 85 filaments at the basal end. Thin sections of sperm fixed during stages in spermatid maturation reveal that the filament bundle begins to assemble on dense material attached to the acrosomal vacuole membrane. In micrographs fo these early stages in maturation, short bundles are seen extending posteriorly from the dense material. The significance is that these short, developing bundles have about 85 filaments, suggesting that the 85-filament end of the bundle is assembled first. By using filament bundles isolated and incubated in vitro with G actin from muscle, we can determine the end “preferred” for addition of actin monomers during polymerization. The end that would be associated with the acrosomal vacuole membrane, a membrane destined to be continuous with the plasma membrane, is preferred about 10 times over the other, thicker end. Decoration of the newly polymerized portions of the filament bundle with subfragment 1 of myosin reveals that the arrowheads point away from the acrosomal vacuole membrane, as is true of other actin filament bundles attached to membranes. From these observations we conclude that the bundle is nucleated from the dense material associated with the acrosomal vacuole and that monomers are added to the membrane-associated end. As monomers are added at the dense material, the thick first-made end of the filament bundle is pushed down through the nucleus where, upon reaching the base of the nucleus, it coils up. Tapering is brought about by the capping of the peripheral filaments in the bundle.     

A novel assembly pipeline and functional annotations for targeted sequencing: A case study on the globally threatened Margaritiferidae (Bivalvia: Unionida)

The proliferation of genomic sequencing approaches has significantly impacted the field of phylogenetics. Target capture approaches provide a cost-effective, fast and easily applied strategy for phylogenetic inference of non-model organisms. However, several existing target capture processing pipelines are incapable of incorporating whole genome sequencing (WGS). Here, we develop a new pipeline for capture and de novo assembly of the targeted regions using whole genome re-sequencing reads. This new pipeline captured targeted loci accurately, and given its unbiased nature, can be used with any target capture probe set. Moreover, due to its low computational demand, this new pipeline may be ideal for users with limited resources and when high-coverage sequencing outputs are required. We demonstrate the utility of our approach by incorporating WGS data into the first comprehensive phylogenomic reconstruction of the freshwater mussel family Margaritiferidae. We also provide a catalogue of well-curated functional annotations of these previously uncharacterized freshwater mussel-specific target regions, representing a complementary tool for scrutinizing phylogenetic inferences while expanding future applications of the probe set.     

P1 and NR1 Plasmid Replication during the Cell Cycle of Escherichia coli

Replication patterns of the miniP1 plasmid pZC176, the miniNR1 plasmid pRR933, and the high-copy miniNR1 derivative pRR942 were examined during the Escherichia coli cell division cycle and compared to the cycle-specific replication pattern of a minichromosome and the cycle nonspecific pattern of pBR322. In E. coli cells growing with doubling times of 40 and 60 min, the miniP1 plasmid was found to replicate with a slight periodicity during the division cycle. The periodicity was not nearly as pronounced as that of the minichromosome, was not affected by the presence of a minichromosome, and was not evident in cells growing more rapidly with a doubling time of 25 min. Both miniNR1 plasmids, pRR933 and pRR942, replicated with patterns indistinguishable from that of pBR322 and clearly different from that of the minichromosome. It is concluded that both P1 and NR1 plasmids can replicate at all stages of the cell cycle but that P1 displays a slight periodicity in replication probability in the cycle of slower growing cells. This periodicity does not appear to be coupled to a specific age in the cycle, but could be associated with the achievement of a specific cell mass per plasmid. During temperature shifts of a dnaC(Ts) mutant, the miniP1 plasmid and pBR322 replicated with similar patterns that differed from that of the minichromosome, but were consistent with a brief eclipse between rounds of replication.     

Mountaintop island age determines species richness of boreal mammals in the American Southwest

Models that describe the mechanisms responsible for insular patterns of species richness include the equilibrium theory of island biogeography and the nonequilibrium vicariance model. The relative importance of dispersal or vicariance in structuring insular distribution patterns can be inferred from these models. Predictions of the alternative models were tested for boreal mammals in the American Southwest. Age of mountaintop islands of boreal habitat was determined by constructing a geographic cladogram based on characteristics of intervening valley barriers. Other independent variables included area and isolation of mountaintop islands. Island age was the most important predictor of species richness. In contrast with previous studies of species richness patterns in this system, these results supported the nonequilibrium vicariance model, which indicates that vicariance has been the primary determinant of species distribution patterns in this system.     

Partially glucose-capped oligosaccharides are found on the hemoglobins of the deep-sea tube worm Riftia pachyptila

We report here the structural determination of N-linked oligosaccharides found on extracellular hemoglobins of the hydrothermal vent tube worm Riftia pachyptila. Structures were elucidated by a combination of electrospray ionization tandem mass spectrometry, matrix- assisted laser desorption/ionization mass spectrometry, normal-phase high performance liquid chromatography, and exoglycosidase digestion. The sugar chains were found to consist mainly of high-mannose-type glycans with some structures partially capped by one or two terminal glucose residues. The present study represents the first report of the occurrence of glucose capping of N-linked carbohydrates in a secreted glycoprotein of a metazoan. Previously, glucose capping has only been described for a membrane-bound surface glycoprotein from the unicellular parasite Leishmania mexicana amazonensis.