Computer building and folding of fictitious transfer-RNA sequences

In order to evaluate the common occurrence with which polynucleotides may adopt the cloverleaf configuration, 1150 random sequences were computer built and folded into their most stable secondary structure. Various constraints modulated the generation of the sequences: i) the base-pairing pattern, ii) the nucleotide composition, iii) the presence of assigned bases (modified or not) at certain sites, and iv) the chain length. In many cases, artificial tRNAs appear to require a more complex organization than a cloverleaf pairing scheme to achieve, as do natural molecules, the corresponding secondary structure. Moreover, the preferred foldings of sequences from 50 to 90 nucleotide long without an imposed pairing pattern usually contain two rather than three hairpin-loops. Implications concerning the emergence and the evolution of the protein-synthesis apparatus are discussed.     

Analysis of highly repeated DNA sequences of rat with EcoR1 endonuclease

Cleavage of rat liver nuclear DNA with EcolR1 restriction endonuclease yields 14 discrete fragments ranging from 2300 to 93 base pairs in length, representing approx. 10.5% of the rat genome. Fragments of 1500, 180, and 93 base pairs are reiterated over 100 000 times; fragments of 2300, 880, 290, and 200 base pairs are reiterated over 20 000 times; the remaining fragments are present in over 1000 copies per genome. When compared to whole rate DNA, 11 were 1-5% richer in A . T base pairs and five were 1.5-2.5 times more methylated. From the criteria of the banding patterns in complete and incomplete digests, base composition and extent of methylation, none of these fragments appeared to be generated as oligomers of a basic shorter repeat. The reassociation of EcoR1 fragments was monitored on hydroxyapatite and by S1 nuclease treatment in order to assess band reiteration frequency and the possibility of interpersion or short internal repeats. The renaturation of the four smallest EcoR1 fragments gave no indication of short internal repeats from hyperpolymer formation nor interpersion with lower frequency sequences by size reduction after S1 nuclease treatment. Anomalous renaturation of several large fragments was observed, possibly due to internal repeats.     

Primary structure of the A chain of human complement-classical-pathway enzyme C1r. N-terminal sequences and alignment of autolytic fragments and CNBr-cleavage peptides.

Activated human complement-classical-pathway enzyme C1r has previously been shown to undergo autolytic cleavages occurring in the A chain [Arlaud, Villiers, Chesne & Colomb (1980) Biochim. Biophys. Acta 616, 116-129]. Chemical analysis of the autolytic products confirms that the A chain undergoes two major cleavages, generating three fragments, which have now been isolated and characterized. The N-terminal alpha fragment (approx. 210 residues long) has a blocked N-terminus, as does the whole A chain, whereas N-terminal sequences of fragments beta and gamma (approx. 66 and 176 residues long respectively) do not, and their N-terminal sequences were determined. Fragments alpha, beta and gamma, which are not interconnected by disulphide bridges, are located in this order within C1r A chain. Fragment gamma is disulphide-linked to the B chain of C1r, which is C-terminal in the single polypeptide chain of precursor C1r. CNBr cleavage of C1r A chain yields seven major peptides, CN1b, CN4a, CN2a, CN1a, CN3, CN4b and CN2b, which were positioned in that order, on the basis of N-terminal sequences of the methionine-containing peptides generated from tryptic cleavage of the succinylated (3-carboxypropionylated) C1r A chain. About 60% of the sequence of C1r A chain (440-460 residues long) was determined, including the complete sequence of the C-terminal 95 residues. This region shows homology with the corresponding parts of plasminogen and chymotrypsinogen and, more surprisingly, with the alpha 1 chain of human haptoglobin 1-1, a serine proteinase homologue.     

Mean free energy topology for nucleotide sequences of varying composition based on secondary structure calculations.

The mean free energy generated from the secondary structure of RNA sequences of varying length and composition has been studied by way of probability theory. The expected boundaries or maximal and minimal values of a given distribution are explored and a method for estimating error as a function of the number of shuffled sequences is also examined. For typical nucleotide sequences found in biologically active organisms, the mean free energy, free energy distributions and errors appear to be scalable in terms of a fixed set of algorithm-dependent parameters and the nucleotide composition of the particular sequence under evaluation. In addition, a general semi-analytical formula for predicting the mean free energy is proposed which, at least to first-order approximation, can be used to rapidly predict the mean free energy of any sequence length and composition of RNA. The general methodology appears to be algorithm independent. The results are expected to provide a reference point for certain types of analysis related to structure of RNA or DNA sequences and to assist in measuring the somewhat related matter of complexity in algorithm development. Some related applications are discussed.     

Amino acid sequence studies on the alpha chain of human fibrinogen. Overlapping sequences providing the complete sequence

The complete amino acid sequence of the alpha chain of human fibrinogen has been determined. It contains 610 amino acid residues and has a calculated molecular weight of 66,124. The chain has 10 methionines, and fragmentation with cyanogen bromide yields 11 peptides [Doolittle, R.F., Cassman, K.G., Cottrell, B.A., Friezner, S.J., Hucko, J.T., & Takagi, T. (1977) Biochemistry 16, 1703]. The arrangement of the 11 fragments was determined by the isolation of peptide overlaps from plasmic and staphylococcal protease digests of fibrinogen and/or alpha chains. In addition, certain of the cyanogen bromide fragments, preliminary reports of whose sequences have appeared previously, have been reexamined in order to resolve several discrepancies. The alpha chain is homologous with the beta and gamma chains of fibrinogen, although a large repetitive segment of unusual composition is absent from the latter two chains. The existence of this unusual segment divides the sequence of the alpha chain into three zones of about 200 residues each that are readily distinguishable on the basis of amino acid composition alone.     

The flexibility of alternating dA-dT sequences

The flexibility of alternating poly (dA-dT) has been investigated by the technique of transient electric dichroism. Rotational relaxation times, which are very sensitive to changes in the end-to-end length of flexible polymers, are determined from the field free dichroism decay curves of four, well defined fragments of poly (dA-dT) ranging in size from 136 to 270 base pairs. Persistence lengths, calculated from the results of Hagerman and Zimm (Biopolymers (1981) 29, 1481-1502), are in the range 200-250 A. This makes alternating dA-dT sequences about twice as flexible as naturally occurring, "random" sequence DNA. Considering a bend around a nucleosome, for example, this difference in persistence length translates to an energy difference between poly (dA-dT) and random sequence DNA of 0.17 kT/base pair or 1 kcal per 10 base pair stretch. This energy difference is sufficiently large to suggest that dA-dT sequences could serve as markers in DNA packaging, for example, at sites where DNA must tightly bend to accommodate structures.     

An Eulerian path approach to global multiple alignment for DNA sequences.

With the rapid increase in the dataset of genome sequences, the multiple sequence alignment problem is increasingly important and frequently involves the alignment of a large number of sequences. Many heuristic algorithms have been proposed to improve the speed of computation and the quality of alignment. We introduce a novel approach that is fundamentally different from all currently available methods. Our motivation comes from the Eulerian method for fragment assembly in DNA sequencing that transforms all DNA fragments into a de Bruijn graph and then reduces sequence assembly to a Eulerian path problem. The paper focuses on global multiple alignment of DNA sequences, where entire sequences are aligned into one configuration. Our main result is an algorithm with almost linear computational speed with respect to the total size (number of letters) of sequences to be aligned. Five hundred simulated sequences (averaging 500 bases per sequence and as low as 70% pairwise identity) have been aligned within three minutes on a personal computer, and the quality of alignment is satisfactory. As a result, accurate and simultaneous alignment of thousands of long sequences within a reasonable amount of time becomes possible. Data from an Arabidopsis sequencing project is used to demonstrate the performance.     

Analysis of rat repetitive DNA sequences.

Parameters of repetitive sequence organization have been measured in the rat genome. Experiments using melting, hydroxylapatite binding, and single strand specific nuclease digestion have been used to measure the number, length, and arrangement of repeated DNA sequences. Renaturation and melting or S1 nuclease digestion of 1.0 kbp DNA fragment show about 20% of rat DNA sequences are 3000-fold repeated. Renatured duplexes from 4.0 kbp DNA fragments display two repetitive size fractions after nuclease digestion. About 60% of the repeated sequences are 0.2-0.4 kbp long while the remainder are longer than 1.5 kbp. The arrangement of the repeated sequences has been measured by hydroxylapatite fractionation of DNA fragments of varying lengths bearing a repeated sequence. Repeated DNA sequences are interspersed among 2.5 kbp long nonrepeated sequences throughout more than 70% of the rat genome. There are approximately 350 different 3000-fold short repeated sequences in the rat interspersed among 600,000 nonrepeated DNA sequences.     

Identification of proteins in sequence databases from amino acid composition data.

Having obtained the amino acid composition of a protein, chemists and molecular biologists may wish to identify the protein from this data alone. In general such data will have errors associated with them and the length of the protein may be known only approximately or not at all. In this paper a method is described which enables searching of protein sequence databases for sequences or fragments of sequences which have a composition similar to the one being sought. Such searches are generally quite discriminating as shown by the examples provided. This method has been implemented as part of the computer program Scrutineer and is being freely distributed. It is simple to use.