A new DNA sequence assembly program.

We describe the Genome Assembly Program (GAP), a new program for DNA sequence assembly. The program is suitable for large and small projects, a variety of strategies and can handle data from a range of sequencing instruments. It retains the useful components of our previous work, but includes many novel ideas and methods. Many of these methods have been made possible by the program's completely new, and highly interactive, graphical user interface. The program provides many visual clues to the current state of a sequencing project and allows users to interact in intuitive and graphical ways with their data. The program has tools to display and manipulate the various types of data that help to solve and check difficult assemblies, particularly those in repetitive genomes. We have introduced the following new displays: the Contig Selector, the Contig Comparator, the Template Display, the Restriction Enzyme Map and the Stop Codon Map. We have also made it possible to have any number of Contig Editors and Contig Joining Editors running simultaneously even on the same contig. The program also includes a new 'Directed Assembly' algorithm and routines for automatically detecting unfinished segments of sequence, to which it suggests experimental solutions.     

A DNA sequence handling program

A computer program that aids in recording, editing, and analysis of the base sequences of DNA and RNA is presented. A tape containing copies of the program and the user manual for it are available at cost.     

TRACTS: A program to map oligopurine.oligopyrimidine and other binary DNA tracts

A program to map the locations and frequencies of DNA tracts composed of only two bases ('Binary DNA') is described. The program, TRACTS (URL http://bioportal.weizmann.ac.il/tracts/tracts.html and/or http://bip.weizmann.ac.il/miwbin/servers/tracts) is of interest because long tracts composed of only two bases are highly over-represented in most genomes. In eukaryotes, oligopurine.oligopyrimidine tracts ('R.Y tracts') are found in the highest excess. In prokaryotes, W tracts predominate (A,T 'rich'). A pre-program, ANEX, parses database annotation files of GenBank and EMBL, to produce a convenient one-line list of every gene (exon, intron) in a genome. The main unit lists and analyzes tracts of the three possible binary pairs (R.Y, K.M and S;W). As an example, the results of R.Y tract mapping of mammalian gene p53 is described.     

A computer program for translating DNA sequences into protein.

This paper describes a comprehensive program for translating one or two DNA sequences into amino acid sequences. Written in FORTRAN, it was designed for maximum flexibility of use and easy maintenance, modification and portability. It has full comments throughout.     

MULTAN: a program to align multiple DNA sequences.

I describe a computer program which can align a large number of nucleic acid sequences with one another. The program uses an heuristic, iterative algorithm which has been tested extensively, and is found to produce useful alignments of a variety of sequence families. The algorithm is fast enough to be practical for the analysis of large number of sequences, and is implemented in a program which contains a variety of other functions to facilitate the analysis of the aligned result.     

NEBcutter: A program to cleave DNA with restriction enzymes

NEBcutter, version 1.0, is a program available via a web server (http://tools.neb.com/NEBcutter) that will accept an input DNA sequence and produce a comprehensive report of the restriction enzymes that will cleave the sequence. It produces a variety of outputs including restriction enzyme maps, theoretical digests and links into the restriction enzyme database, REBASE (http://www.neb.com/rebase). Importantly, its table of recognition sites is updated daily from REBASE and it marks all sites that are potentially affected by DNA methylation (Dam, Dcm, etc.). Many options exist to choose the enzymes used for digestion, including all known specificities, subsets of those that are commercially available or sets of enzymes that produce compatible termini.     

A program for selecting DNA fragments to detect mutations by denaturing gel electrophoresis methods.

A computer program was developed to automate the selection of DNA fragments for detecting mutations within a long DNA sequence by denaturing gel electrophoresis methods. The program, MELTSCAN, scans through a user specified DNA sequence calculating the melting behavior of overlapping DNA fragments covering the sequence. Melting characteristics of the fragments are analyzed to determine the best fragment for detecting mutations at each base pair position in the sequence. The calculation also determines the optimal fragment for detecting mutations within a user specified mutational hot spot region. The program is built around the statistical mechanical model of the DNA melting transition. The optimal fragment for a given position is selected using the criteria that its melting curve has at least two steps, the base pair position is in the fragment's lowest melting domain, and the melting domain has the smallest number of base pairs among fragments that meet the first two criteria. The program predicted fragments for detecting mutations in the cDNA and genomic DNA of the human p53 gene.     

A flexible new computer program for handling DNA sequence data.

A compact new computer program for handling nucleic acid sequence data is presented. It consists of a number of different subsets, which may be used according to a given code system. The program is designed for the determination of restriction enzyme and other recognition sites in correlation with translation patterns, and allows tabulation of codon frequencies and protein molecular weights within specified gene boundaries. The program is especially designed for detection of overlapping genes. The language, is FORTRAN and thus the program may be used on small computers; it may also be used without any prior computer experience. Copies are available on request.     

Tandem repeats finder: a program to analyze DNA sequences.

A tandem repeat in DNA is two or more contiguous, approximate copies of a pattern of nucleotides. Tandem repeats have been shown to cause human disease, may play a variety of regulatory and evolutionary roles and are important laboratory and analytic tools. Extensive knowledge about pattern size, copy number, mutational history, etc. for tandem repeats has been limited by the inability to easily detect them in genomic sequence data. In this paper, we present a new algorithm for finding tandem repeats which works without the need to specify either the pattern or pattern size. We model tandem repeats by percent identity and frequency of indels between adjacent pattern copies and use statistically based recognition criteria. We demonstrate the algorithm's speed and its ability to detect tandem repeats that have undergone extensive mutational change by analyzing four sequences: the human frataxin gene, the human beta T cellreceptor locus sequence and two yeast chromosomes. These sequences range in size from 3 kb up to 700 kb. A World Wide Web server interface atc3.biomath.mssm.edu/trf.html has been established for automated use of the program.     

A large database DNA sequence handling program with generalized searching specifications.

The program described allows for the creation and manipulation of files of DNA sequence data up to very great lengths. The program uses its own paging system to load segments of the sequence into a small internal buffer so that the program does not have excessive memory requirements. The program offers a menu of functions to the user, and has been written to be forgiving of user errors. A code for the generalised specification of bases as a series of groups (i.e. A or T, Purine, etc.) has been devised and can be used in search specifications or in sequence files. Versions of the program have been developed to run with special efficiency under DIGITAL's RT11 operating system or to run under systems with a suitable implementation of FORTRAN VI.     

Methylation blockage and other improvements to a comprehensive DNA analysis program.

A comprehensive DNA analysis computer program was described in the second special issue of Nucleic Acids Research on the applications of computers to research on nucleic acids by Stone and Potter (1). Criteria used in designing the program were user friendliness, ability to handle large DNA sequences, low storage requirement, migratability to other computers and comprehensive analysis capability. The program has been used extensively in an industrial-research environment. This paper talks about improvements to that program. These improvements include testing for methylation blockage of restriction enzyme recognition sites, homology analysis, RNA folding analysis, integration of a large DNA database (GenBank), a site specific mutagenesis analysis, a protein database and protein searching programs. The original design of the DNA analysis program using a command executive from which any analytical programs can be called, has proven to be extremely versatile in integrating both developed and outside programs to the file management system employed.     

Information transfer from DNA to peptide nucleic acids by template-directed syntheses.

Peptide nucleic acids (PNAs) are analogs of nucleic acids in which the ribose-phosphate backbone is replaced by a backbone held together by amide bonds. PNAs are interesting as models of alternative genetic systems because they form potentially informational base paired helical structures. Oligocytidylates have been shown to act as templates for formation of longer oligomers of G from PNA G2 dimers. In this paper we show that information can be transferred from DNA to PNA. DNA C4T2C4 is an efficient template for synthesis of PNA G4A2G4 using G2 and A2 units as substrates. The corresponding synthesis of PNA G4C2G4 on DNA C4G2C4 is less efficient. Incorporation of PNA T2 into PNA products on DNA C4A2C4 is the least efficient of the three reactions. These results, obtained using PNA dimers as substrates, parallel those obtained using monomeric activated nucleotides.     

Cooperative, non-specific binding of a zinc finger peptide to DNA.

The DNA binding and structural properties of Xfin-31 (Lee, M.S., Gippert, G.P., Soman, K.V., Case, D.A. and Wright, P.E., 1989, Science 245, 635-637), a twenty five amino acid zinc finger peptide, in the reduced, oxidized and zinc complex forms, as well as the fourteen residue helical segment of the zinc finger (residues 12-25) have been compared using affinity coelectrophoresis (ACE) and circular dichroism (CD) spectroscopy. The zinc complex and oxidized peptides bind cooperatively to DNA although the cooperativity factor, omega, is more than 15-fold greater for the zinc complex. The reduced peptide in the absence of zinc and the helical segment do not bind cooperatively (omega = 1). Hence, the binding constant for singly contiguous sites (K omega) ranges over 100-fold for the various peptides even though the intrinsic binding constants (K) are similar. An increase in binding order and affinity for the other forms of Xfin-31 is correlated with an increasing similarity of the CD spectrum to that of the Xfin-31 zinc complex. The surprising DNA binding activity of the oxidized peptide may result from hydrophobic interactions between the amino-terminal loop formed by the Cys3-Cys6 disulfide bond and conserved hydrophobic residues in the carboxyl-terminal segment. Xfin-31 may be a particularly useful model for studying several poorly understood aspects of cooperative, non-specific DNA binding since it is small, has a stable, well-defined structure, and structures of zinc fingers bound to DNA have been determined.     

Proteolytic conversion of rat thymus DNA polymerase alpha to a more active form.

The relationship between two DNA polymerase alpha species from mammalian tissues has been resolved with the isolation of a protease from rat thymus which converts the larger alpha polymerase (7.3S) to a smaller (5.4S) size. The proteolytic activity is present only in the chromatin fraction and the limited proteolysis is accompanied by an increase in activity of the DNA polymerase, possibly consistent with a biological control function for this phenomenon.     

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 program for the identification of tRNA-like structures in DNA sequence data.

A computer algorithm has been developed which identifies tRNA genes and tRNA-like structures in DNA sequences. The program searches the sequence string for specific base positions that correspond to the invariant and semi-invariant bases found in tRNAs. The tRNA nature of the sequence is confirmed by the presence of complementary base pairing at the tRNA's calculated 5' and 3' ends (which in situ constitutes the amino-acyl stem region). The program achieves greater than 96% accuracy when run against known tRNA sequences in the Genbank database. The program is modular and is readily modified to allow searching either a file or database. The program is written in "C" and operates on a D.E.C. Vax 750. The utility of the algorithm is demonstrated by the identification of a distinctive tRNA structure in an intron of a published bovine hemoglobin gene.     

New peptide conjugates with 9-aminoacridine: synthesis and binding to DNA.

New peptides-9-aminoacridine conjugates with an ethylene diamine linker-have been synthesized (both solution and solid phase methods were used) and their interactions with DNA have been studied. The affinity of H-Phe-Gln-Gly-Ile(2)-NHCH(2)CH(2)NH-Acr conjugate and of its extended analogue containing 6-aminohexanoic acid to DNA were lower than that of a standard H-Gly-NHCH(2)CH(2)NH-Acr conjugate. The results fit well into our concept of peptide conjugates with lowered binding activity to DNA, which could be capable of unlimited extravascular distribution. Moreover, new structures could be potentially useful as the mild tuners of DNA interaction with strong bis-acridine binders.     

A peptide released from plasma fibrin stabilzing factor in the conversion to the active enzyme by thrombin

A peptide, which was released accompanying with the activation of bovine plasma fibrin stabilizing factor (FSF) by thrombin, was isolated and characterized. The peptide consisted of Asp4, Thr₃, Ser₄, Glu₄, Pro₅, Gly₄, Ala₄, Val₂, Ile₁, Leu₂, Phe₁, and Arg₃. The content of proline was highest in all of these amino acids. The carboxyl-terminal residue of the peptide was identified as arginine. However, no N-terminal amino acid reactive with phenyl$\text{iso}$thiocyanate and dansyl chloride could be determined. Edman degradation on the inactive FSF showed glutamic acid or glutamine as one N-terminal residue. After the activation of FSF by thrombin, glycine was identified as a second N-terminal residue, in addition to glutamic acid (glutamine).These results indicate that the transformation of FSF to the active enzyme by thrombin involves proteolysis of an arginyl-glycyl bond located in the N-terminal region of one of the subunits of the proenzyme.