IBM microcomputer programs that analyze DNA sequences for tRNA genes

A set of four computer programs that search DNA sequence datafiles for transfer RNA genes have been written in IBM (Microsoft)BASIC for the IBM personal computer. These programs locate andplot predicted secondary structures of tRNA genes in the cloverleafconformation. The set of programs are applicable to eukaryotictRNA genes, including those containing intervening sequences,and to prokaryotic and mitochondrial tRNA genes. In addition,two of the programs search up to 150 residues downstream oftRNA gene sequences for possible eukaryotic RNA polymerase IIItermination sites comprised of at least four consecutive T residues.Molecular biologists studying a variety of gene sequences andflanking regions can use these programs to search for the additionalpresence of tRNA genes. Furthermore, investigators studyingtRNA gene structure-to-function relationships would not needto do extensive restriction mapping to locate tRNA gene sequenceswithin their cloned DNA fragments. Received on October 29, 1985; accepted on January 28, 1986     

Automated preparation of DNA sequences for publication.

A computer program which draws DNA sequences is described. A simple method is used which enables the user to highlight or annotate specific parts of a sequence. The sizes of the characters in the sequence to be drawn are specified by the user. In addition, vertical spacing between lines and horizontal spacing between characters can be specified. Sequences can be prepared and high quality output produced on a plotter in a short period of time, making the program advantageous to use over typing, computer printing, or preparation by a graphics department.     

Two computer programs for rapid entry of DNA sequence data

Two computer programs for the IBM personal computer are describedfor rapid and accurate entry of DNA sequence data. The DNA sequencefiles produced can be used directly by the DNA sequence manipulationprograms by R. Staden (the DataBase system), the Universityof Wisconsin Genetics Computer Group, DNASTAR, or D. Mount.The first program, DIGISEQ, utilizes a sonic digitizer for semi-automationof sequence entry. To enter the DNA sequence each band of agel reading is touched by the stylus of the sonic digitizer.DIGISEQ corrects for both changes in lane width and lane curvature.The algorithm is extremely efficient and rarely requires re-entenngthe centers of the lanes. The second program, TYPESEQ, usesonly the keyboard for input. The keyboard is reconfigured toplace nucleotides and ambiguity codes under the fingers of onehand, corresponding to the order of the nucleotides on the geldefined by the user Both programs produce individual tones foreach nucleotide, and certain ambiguity codes. This verifiesinput of the correct nucleotide or ambiguity code, and thuseliminates the need to visually check the screen display duringsequence entry. Received on November 16, 1986; accepted on June 16, 1987     

Spreadsheet-based program for alignment of overlapping DNA sequences.

Molecular biology laboratories frequently face the challenge of aligning small overlapping DNA sequences derived from a long DNA segment. Here, we present a short program that can be used to adapt Excel spreadsheets as a tool for aligning DNA sequences, regardless of their orientation. The program runs on any Windows or Macintosh operating system computer with Excel 97 or Excel 98. The program is available for use as an Excel file, which can be downloaded from the BioTechniques Web site. Upon execution, the program opens a specially designed customized workbook and is capable of identifying overlapping regions between two sequence fragments and displaying the sequence alignment. It also performs a number of specialized functions such as recognition of restriction enzyme cutting sites and CpG island mapping without costly specialized software.     

PERFILS: a program for the quantitative treatment of footprinting data

PERFILS, a computer program written in Borland TurboPascal,performs quantitative analysis of footprinting experiments usingany IBM PC or compatible microcomputer. The program uses theheight of the bands obtained from densitometric scanning offootprinting autoradiographs to calculate a differential cleavageplot. Such a plot displays, on a logarithmic scale, the differenceof susceptibility of a DNA fragment to DNase I, or any othercleaving agent, in the presence of any ligand versus the sequence.PERFILS calculates the fractional cleavage values for controland ligand, giving a table of values for each internucleotidicbond and rendering the differential cleavage plot in only afew seconds.     

Fast computer search for similar DNA sequences.

An extremely fast method of searching a nucleic acid sequence database against a probe sequence is described. The method is based on the detection of deviation from expected number and deviation from random spatial distribution of sub-sequences which are unique within a sequence, and shared between that sequence and the probe. On an IBM 3081 computer, total search of an encoded form of the EMBL nucleic acid sequence database with a 1 kbase probe sequence is completed in a few seconds. Previous best methods for a similar task required a few minutes.     

A tool for aligning very similar DNA sequences

Results: We have produced a computer program, named sim3, thatsolves the following computational problem. Two DNA sequencesare given, where the shorter sequence is very similar to somecontiguous region of the longer sequence. Sim3 determines sucha similar region of the longer sequence, and then computes anoptimal set of single-nucleotide changes (i.e. insertions, deletionsor substitutions) that will convert the shorter sequence tothat region. Thus, the alignment scoring scheme is designedto model sequencing errors, rather than evolutionary processes.The program can align a 100 kb sequence to a 1 megabase sequencein a few seconds on a workstation, provided that there are veryfew differences between the shorter sequence and some regionin the longer sequence. The program has been used to assemblesequence data for the Genomes Division at the National Centerfor Biotechnology Information. Availability: A version of sim3 for UNIX machines can be obtainedby anonymous ftp from ncbi. nlm. nih. gov, in the pub/sim3 directory. Contact: For portable versions for Macs and PCs, contact zjing@sunset.nlm. nih. gov.     

A computer program that simulates cloning of DNA

Easy Cloner is a computer program that manipulates DNA sequences as in cloning experiments and produces maps of the resulting plasmids. The program runs in the graphics mode of an IBM PC or compatible computer and is operated by using a mouse to point to the required actions. The program is available in the public domain.     

A Macintosh computer program for designing DNA sequences that code for specific peptides and proteins

A computer program (PINCERS) is described for use in the design of synthetic genes and mixed-probe DNA sequences. A protein sequence is reverse translated with generation of synonymous codons at each position producing a degenerate sequence. In order to locate potential restriction enzyme sites, the degenerate sequence is searched with a library of restriction enzymes for sites that utilize any combination of synonymous codons. These sites are indicated in a map so that they may be incorporated into the synthetic gene sequence. The program allows the user to select the appropriate codon usage table for the organism of interest and then to set a threshold usage frequency below which codons are not generated. PINCERS may also be used to assist in planning the synthesis of mixed-probe DNA sequences for cross-hybridization experiments. It can identify regions of specified length with the protein sequence that have the least overall degeneracy, thereby minimizing the number of probes to be synthesized and, therefore, maximizing the concentration of a given probe sequence.     

HYLAS: program for generating H curves (abstract three-dimensional representations of long DNA sequences)

The computer program HYLAS generates from a standard DNA lettersequence a three-dimensional space curve (H curve) which embodiesthe entire information content of the original nucleotide sequence.The program can display H curves either as two-dimensional (frontand side view) projections or as stereo-pair images. The curvescan be marked at specific nucleotide locations, annotated, rotatedfor observation from any viewing angle, and manipulated forconvenient side-by-side comparisons. Unlike the cumbersome lettersequences, H curves can be drastically condensed in size withoutlosing their ability to reflect the global nucleotide-distributionpattern of the entire DNA sequence. Often, biologically importantloci can be visually identified on the H curves. HYLAS is writtenin FORTRAN with separate mainframe (IBM- VM/CMS) and microcomputer(MS-DOS) versions. It uses the Tektronix-TCS library of graphicsubroutines. Received on October 24, 1988; accepted on July 15, 1989     

The relationship between base composition and codon usage in bacterial genes and its use for the simple and reliable identification of protein-coding sequences

Bacterial genes that code for proteins appear to possess a codon usage characteristic of their overall base composition. This results in different but predictable non-random distributions of nucleotides within codons, permitting the recognition of protein-coding sequences in a wide range of bacterial species. The nature of this distribution depends on the base composition of the coding sequence. The position-specific differences are especially conspicuous in genes of extreme G + C content, allowing the particularly reliable prediction of the reading frame and coding strand of experimentally determined DNA sequences. This fmding has been exploited to identify the coding sequence of the viomycin phosphotransferase (vph) gene of Streptomyces vinaceus. An easily applied computer program (“Frame”) has been written to carry out and display such analyses.     

Whole genome sequence-enabled prediction of sequences performed for random PCR products of Escherichia coli

The sequence of an unknown PCR product generated by random (and conventional) PCR could be determined without sequencing when it is provided with the template DNA sequence. Theoretically, this was based on formerly established ideas which assert that the amount of random PCR product mainly depends on the stability of the primer-binding structures and that the dynamic solution structure of DNA is essentially governed by the Watson–Crick base pairing. However, it has not been clear whether this holds true for larger genomes of mega- to gigabase size, beside the λ phage genome (of 50 kb) used previously, nor has it been ascertained to uniquely specify the sequence of a random PCR product. Here, we jointly use two computer programs together with experimental data from Genome Profiling (i.e. TGGE analysis of random PCR products). The first procedure carried out by a newly remodeled computer program (PCRAna-A1) was shown to be competent to calculate a set of random PCR products from Escherichia coli genome DNA (4.7 Mb). The other procedure performed with another program (Poland-H) played a critical role in determining the final candidate sequence by theoretically offering the initial melting temperature and the melting pattern of unspecified candidate sequences. The success attained here not only proved our method to be useful for sequence prediction but also confirmed the above-mentioned ideas as rational. We believe that this is the first case to computer-utilize a genome sequence as a whole.     

A new method for DNA surface-property display and its application to E. coli promoter sequences

A personal computer program to visualize and compare the property of double-stranded DNA surface has been developed. Comparison of the surface property between Watson-Crick base-pairs in B-form DNA has elucidated that the base-pair replacement between a "degenerated base-pairs" conserves the pattern of potential hydrogen-bonding sites in both major and minor grooves. The idea of the "degenerated base-pairs" was applied for the problem of the base-sequence variation from the consensus sequence in the -35 region of E. coli promoter. The sequence variation is found to have tendency to occur among the degenerated base-pairs.     

DNA保守序列识别算法的并行化和MPI集群环境构建

DNA序列中保守序列的识别需要较大的计算量。开发了一个转录因子结合位点识别的并行算法,能够从多条DNA序列中识别指定长度的序列模式。算法使用概率模型进行序列模式保守性的度量,利用迭代过程实现保守序列的搜索。使用C编程结合MPI消息传递模型开发了相应的程序,并在Windows平台下构建了一个3节点的集群环境,利用20个长度均为200的序列数据集进行测试,实现了模体识别工作,结果表明并行算法使模体识别的效率得到提高。     

Investigating extended regulatory regions of genomic DNA sequences.

MOTIVATION: Despite the growing volume of data on primary nucleotide sequences, the regulatory regions remain a major puzzle with regard to their function. Numerous recognising programs considering a diversity of properties of regulatory regions have been developed. The system proposed here allows the specific contextual, conformational and physico-chemical properties to be revealed based on analysis of extended DNA regions. RESULTS: The Internet-accessible computer system RegScan, designed to analyse the extended regulatory regions of eukaryotic genes, has been developed. The computer system comprises the following software: (i) programs for classification dividing a set of promoters into TATA-containing and TATA-less promoters and promoters with and without CpG islands; (ii) programs for constructing (a) nucleotide frequency profiles, (b) sequence complexity profiles and (c) profiles of conformational and physico-chemical properties; (iii) the program for constructing the sets of degenerate oligonucleotide motifs of a specified length; and (iv) the program searching for and visualising repeats in nucleotide sequences. The system has allowed us to demonstrate the following characteristic patterns of vertebrate promoter regions: the TATA box region is flanked by regions with an increased G+C content and increased bending stiffness, the TATA box content is asymmetric and promoter regions are saturated with both direct and inverted repeats. AVAILABILITY: The computer system RegScan is available via the Internet at http://www.mgs.bionet.nsc. ru/Systems/RegScan, http://www.cbil.upenn.edu/mgs/systems/r egscan/.     

A computer aided thermodynamic approach for predicting the formation of Z-DNA in naturally occurring sequences.

The ease with which a particular DNA segment adopts the left-handed Z-conformation depends largely on the sequence and on the degree of negative supercoiling to which it is subjected. We describe a computer program (Z-hunt) that is designed to search long sequences of naturally occurring DNA and retrieve those nucleotide combinations of up to 24 bp in length which show a strong propensity for Z-DNA formation. Incorporated into Z-hunt is a statistical mechanical model based on empirically determined energetic parameters for the B to Z transition accumulated to date. The Z-forming potential of a sequence is assessed by ranking its behavior as a function of negative superhelicity relative to the behavior of similar sized randomly generated nucleotide sequences assembled from over 80,000 combinations. The program makes it possible to compare directly the Z-forming potential of sequences with different base compositions and different sequence lengths. Using Z-hunt, we have analyzed the DNA sequences of the bacteriophage phi X174, plasmid pBR322, the animal virus SV40 and the replicative form of the eukaryotic adenovirus-2. The results are compared with those previously obtained by others from experiments designed to locate Z-DNA forming regions in these sequences using probes which show specificity for the left-handed DNA conformation.     

SILMUT: a computer program for the identification of regions suitable for silent mutagenesis to introduce restriction enzyme recognition sequences.

We describe a set of IBM-compatible computer programs designed to selectively identify the potential sites for silent mutagenesis within a target DNA sequence. This program is based on a novel strategy of identifying amino acid motifs compatible with each restriction site (BioTechniques 12:382-384, 1991). The programs can be used to identify the suitability for the introduction of any 6-base nucleic acid sequences, such as restriction enzyme sites in cassette mutagenesis strategies. The Table program generates a table of multiple amino acid motifs for each restriction enzyme, obtained by translating each unique recognition sequence in all three reading frames. The Silmut program, which utilizes the features of Table, will further identify the presence of a match between any amino acid motif of each restriction enzyme and the input target sequence. Minor manipulations of the data base files will enable the individual researcher to identify the potential for introduction of any 6-base sequences by silent mutagenesis.