BIGPROBE: a computer program that predicts the sequence of long oligonucleotide probes with high reliability

We have written a computer program, BIGPROBE, which facilitates the design of long nucleic acid probes from the partial or complete amino acid sequence of a protein. BIGPROBE relies upon information on codon usage, intercodon dinucleotide frequency, and potential probe self-complementarity. We have examined the accuracy with which the program predicts coding sequences using sample human and rat genes and probe lengths of 30-60 nucleotides. Rat probe sequences selected by BIGPROBE using either codon usage or dinucleotide frequency data alone averaged 86-92% homology with the known exons of the corresponding gene sequences. Predictive accuracy with rat gene probes could be improved to 89-94%, depending upon probe length, by applying codon usage and dinucleotide frequency data in combination. Similar accuracy was achieved for human genes.     

Nucleic acid and protein sequence databases

Nucleic acid and protein sequences contain a wealth of informationof interest to molecular biologists. The advent of molecularsequence databases provides a unique opportunity for the computeranalysis of all available sequences. Sequence databases servetwo main functions: (i) to facilitate comparisons with newlydetermined sequences, and (ii) to act as a source of data forthe generation and testing of hypotheses concerning molecularsequence organisation and evolution. The large amounts of sequencedata now becoming available require that algorithms for databasesearching be fast and efficient and considerable progress isbeing made in this area.     

Secondary structure-based profiles: use of structure-conserving scoring tables in searching protein sequence databases for structural similarities

The profile method, for detecting distantly related proteins by sequence comparison, has been extended to incorporate secondary structure information from known X-ray structures. The sequence of a known structure is aligned to sequences of other members of a given folding class. From the known structure, the secondary structure (alpha-helix, beta-strand or "other") is assigned to each position of the aligned sequences. As in the standard profile method, a position-dependent scoring table, termed a profile, is calculated from the aligned sequences. However, rather than using the standard Dayhoff mutation table in calculating the profile, we use distinct amino acid mutation tables for residues in alpha-helices, beta-strands or other secondary structures to calculate the profile. In addition, we also distinguish between internal and external residues. With this new secondary structure-based profile method, we created a profile for eight-stranded, antiparallel beta barrels of the insecticyanin folding class. It is based on the sequences of retinol-binding protein, insecticyanin and beta-lactoglobulin. Scanning the sequence database with this profile, it was possible to detect the sequence of avidin. The structure of streptavidin is known, and it appears to be distantly related to the antiparallel beta barrels. Also detected is the sequence of complement component C8, which we therefore predict to be a member of this folding class.     

PATMAT: a searching and extraction program for sequence, pattern and block queries and databases

A program has been developed that provides molecular biologistswith multiple tools for searching databases, yet uses a verysimple interface. PATMATcan use protein or (translated) DNAsequences, patterns or blocks of aligned proteins as queriesof databases consisting of amino acid or nucleotide sequences,patterns or blocks. The ability to search databases of blocksby ‘on-the-fly’ conversion to scoring matrices providesa new tool for detection and evaluation of distant relationships.PATMAT uses a pull-down, menu-driven interface to carry outits multiple searching, extraction and viewing functions. Eachquery or database type is recognized, reported, and the appropriatesearch carried out, with matches and alignments reported inwindows as they occur. Any of the high scoring matches can beexported to a file, viewed and recalled as a query using onlya few keystrokes or mouse selections. Searches of multiple databasefiles are carried out by user selection within a window. PATMATruns under DOS; the searching engine also runs under UNIX.     

3MOTIF: visualizing conserved protein sequence motifs in the protein structure database

SUMMARY: 3MOTIF is a web application that visually maps conserved sequence motifs onto three-dimensional protein structures in the Protein Data Bank (PDB; Berman et al., Nucleic Acids Res., 28, 235-242, 2000). Important properties of motifs such as conservation strength and solvent accessible surface area at each position are visually represented on the structure using a variety of color shading schemes. Users can manipulate the displayed motifs using the freely available Chime plugin. AVAILABILITY: http://motif.stanford.edu/3motif/     

Distribution and phasing of sequence motifs that facilitate CRISPR adaptation

1. Download : Download high-res image (184KB)
2. Download : Download full-size image

     

DNABIND: an interactive microcomputer program searching for nucleotide sequences that may code for conserved DNA-binding protein motifs

This paper presents a simple program for interactive searchingfor nucleotide sequences that may code for the helix—turn—helix,zinc finger or leucine zipper motifs in proteins. The helix—turn—helixmotifs are predicted using the recently published method ofDodd and Egan, while zinc fingers and leucine zippers are searchedfor by our original methods. DNABIND is shown to detect allfour known helix—turn—helix motifs in bacteriophagelambda genes and both zinc fingers of the adrl gene of yeast.     

RSDB: representative protein sequence databases have high information content

MOTIVATION: Biological sequence databases are highly redundant for two main reasons: 1. various databanks keep redundant sequences with many identical and nearly identical sequences 2. natural sequences often have high sequence identities due to gene duplication. We wanted to know how many sequences can be removed before the databases start losing homology information. Can a database of sequences with mutual sequence identity of 50% or less provide us with the same amount of biological information as the original full database? RESULTS: Comparisons of nine representative sequence databases (RSDB) derived from full protein databanks showed that the information content of sequence databases is not linearly proportional to its size. An RSDB reduced to mutual sequence identity of around 50% (RSDB50) was equivalent to the original full database in terms of the effectiveness of homology searching. It was a third of the full database size which resulted in a six times faster iterative profile searching. The RSDBs are produced at different granularity for efficient homology searching. AVAILABILITY: All the RSDB files generated and the full analysis results are available through internet: ftp://ftp.ebi.ac. uk/pub/contrib/jong/RSDB/http://cyrah.e bi.ac.uk:1111/Proj/Bio/RSDB     

PROMOT: a FORTRAN program to scan protein sequences against a library of known motifs

Information about the three-dimensional structure or functionof a newly determined protein sequence can be obtained if theprotein is found to contain a characterized motif or patternof residues. Recently a database (PROSITE) has been establishedthat contains 337 known motifs encoded as a list of allowedresidue types at specific positions along the sequence. PROMOTis a FORTRAN computer program that takes a protein sequenceand examines if it contains any of the motifs in PROSITE. Theprogram also extends the definitions of patterns beyond thoseused in PROSITE to provide a simple, yet flexible, method toscan either a PROSITE or a user-defined pattern against a proteinsequence database. Received on October 17, 1990; accepted on November 15, 1990     

ADEPT: a computer program for prediction of protein antigenic determinants

ADEPT, a program that can be used for prediction of proteinantigenic determinants from the amino acid sequence alone forthe cases of humoral and cellular immune response, is described.Most methods presently utilized for this purpose are implementedin the program along with some original parameters for the caseof humoral immune response. There is also a possibility to combineseveral methods for the case of humoral immune response, whichprovides a means to create new predictive methods. ADEPT alsoincludes a literature-derived database of proteins in the SWISS-PROTstandard with experimentally determined antigenic determinants,so the predictive ability of new methods can be assessed withinthe program and corresponding statistical information is calculated.ADEPT may be useful for solving vamious tasks involving thedelineation of antigenic regions on proteins or characterizationof peptide fragments, arising, for examnple, in protein engineering(e.g. selection of protein portions most suitable for insertionof peptide sequences that are antigenically active or possessingother desired characteristics), and in particular for constructionof vaccines containing B- and/or T-cell epitopes, as well asfor other problems of medico-biological research.     

BBReader: A computer program for the combined use of the BioMagResBank and PDB databases

A computer program (BBReader) was developed which performs an inverse search in theBioMagResBank database. Given (cross) peak positions of a protein, the program searchesfor atoms with matching chemical shifts and suggests possible assignments for user-specifiedhomo- and heteronuclear one- to three-dimensional COSY- and NOESY-type experiments.It can handle 1H, 13C and 15N spectra. Distance information from PDB files can be utilizedfor filtering possible NOESY cross peak assignments.     

3MATRIX and 3MOTIF: a protein structure visualization system for conserved sequence motifs

Computational methods such as sequence alignment and motif construction are useful in grouping related proteins into families, as well as helping to annotate new proteins of unknown function. These methods identify conserved amino acids in protein sequences, but cannot determine the specific functional or structural roles of conserved amino acids without additional study. In this work, we present 3MATRIX (http://3matrix.stanford.edu) and 3MOTIF (http://3motif.stanford.edu), a web-based sequence motif visualization system that displays sequence motif information in its appropriate three-dimensional (3D) context. This system is flexible in that users can enter sequences, keywords, structures or sequence motifs to generate visualizations. In 3MOTIF, users can search using discrete sequence motifs such as PROSITE patterns, eMOTIFs, or any other regular expression-like motif. Similarly, 3MATRIX accepts an eMATRIX position-specific scoring matrix, or will convert a multiple sequence alignment block into an eMATRIX for visualization. Each query motif is used to search the protein structure database for matches, in which the motif is then visually highlighted in three dimensions. Important properties of motifs such as sequence conservation and solvent accessible surface area are also displayed in the visualizations, using carefully chosen color shading schemes.     

Definition of protein kinase sequence motifs that trigger high affinity binding of Hsp90 and Cdc37

Hsp90 cooperates with its co-chaperone Cdc37 to provide obligatory support to numerous protein kinases involved in the regulation of cellular signal transduction pathways. In this report, the crystal structure of the Src family tyrosine kinase Lck was used to guide the creation of kinase constructs to determine features recognized by Hsp90 and its "kinase-specific" co-chaperone Cdc37. Two parameters were assayed: the ability and extent to which the constructs bound to Hsp90 and Cdc37, and the ability of the constructs to trigger salt-resistant high affinity complexes with Hsp90 and Cdc37 independent of the presence of molybdate. Although Hsp90 interacted with both the N-terminal and C-terminal lobes (NL and CL, respectively) of the catalytic domains of the kinases, the lobes themselves were not sufficient to trigger the high affinity binding of Hsp90. Only constructs containing a complete N- or C-terminal lobe and part of the adjacent lobe bound to Hsp90 and Cdc37 in salt-stable complexes independent of molybdate. The two minimum constructs that bound Hsp90 and Cdc37 contained the alpha-C-helix and the beta4- and beta5-strands of the NL through to end of the CL and the NL through to the alpha-E-helix and the amino acids that cap the helix. Cdc37 interacted with only the NL and minimally required the alpha-C-helix and beta4- and beta5-strands of this lobe of Lck. The results indicate that the high affinity binding activity of Hsp90 is triggered through its interaction with adjacent subdomain structures of kinase catalytic domains. Furthermore, the alpha-C-helix and part of its adjoining loop connection to the beta4-strand appear to be the primary determinants recognized by Cdc37.     

Nup155 is a novel nuclear pore complex protein that contains neither repetitive sequence motifs nor reacts with WGA

We have molecularly cloned and sequenced a rat liver nuclear pore complex (NPC) protein of calculated molecular mass of 155 kD. Consistent with recently proposed nomenclature this protein is termed nucleoporin 155, or nup155. Unlike other nups that have so far been molecularly cloned and sequenced, nup155 does not contain repetitive sequence domains. It does not show similarity to the sequences of other proteins, including any nups, so far compiled in the data bases. Like other vertebrate nups which have been characterized nup155 possesses abundant (46 in total) consensus sites for various kinases. By immunoelectron microscopy, nup155 is associated with both the nucleoplasmic and the cytoplasmic aspect of the NPC and is therefore possibly a component of the symmetrically arranged NPC substructures. In mitotic cells, nup155 assumes a diffuse cytoplasmic distribution. Nup155 is among the integral of 30 proteins that were extracted from rat liver nuclear envelopes by 2.0 M urea/1.0 mM EDTA, separated from WGA-reactive proteins by WGA-Sepharose and further subfractionated by SDS-hydroxylapatite. These proteins are potential candidates for being nups.     

RSITE: a computer program to predict the recognition sequence of a restriction enzyme.

A computer program (RSITE) was developed which predicts the recognition sequence of a restriction endonuclease. The sizes of fragments experimentally determined on cleavage of a DNA of known sequence were input. Possible recognition sequences producing fragments of sizes matching those determined empirically were printed out. The program faithfully predicted the specificity of restriction enzymes of known recognition sequence and also determined the recognition sequence of a new restriction enzyme from Haemophilus influenzae GU (HinGU II).