The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

/pdf/clustal-w-improving-the-sensitivity-of-progressive-multiple-2na2s9sq2h.pdf

Clustal w: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

Some simple formulae were obtained which enable us to estimate evolutionary distances in terms of the number of nucleotide substitutions (and, also, the evolutionary rates when the divergence times are known). In comparing a pair of nucleotide sequences, we distinguish two types of differences; if homologous sites are occupied by different nucleotide bases but both are purines or both pyrimidines, the difference is called type I (or “transition” type), while, if one of the two is a purine and the other is a pyrimidine, the difference is called type II (or “transversion” type). Letting P and Q be respectively the fractions of nucleotide sites showing type I and type II differences between two sequences compared, then the evolutionary distance per site is K = — (1/2) ln {(1 — 2P — Q) }. The evolutionary rate per year is then given by k = K/(2T), where T is the time since the divergence of the two sequences. If only the third codon positions are compared, the synonymous component of the evolutionary base substitutions per site is estimated by K'S = — (1/2) ln (1 — 2P — Q). Also, formulae for standard errors were obtained. Some examples were worked out using reported globin sequences to show that synonymous substitutions occur at much higher rates than amino acid-altering substitutions in evolution.

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences.

The University of Wisconsin Genetics Computer Group (UWGCG) has been organized to develop computational tools for the analysis and publication of biological sequence data. A group of programs that will interact with each other has been developed for the Digital Equipment Corporation VAX computer using the VMS operating system. The programs available and the conditions for transfer are described.

/pdf/a-comprehensive-set-of-sequence-analysis-programs-for-the-luj6w12mfm.pdf

A comprehensive set of sequence analysis programs for the VAX

The abbreviated name,‘mfold web server’,describes a number of closely related software applications available on the World Wide Web (WWW) for the prediction of the secondary structure of single stranded nucleic acids. The objective of this web server is to provide easy access to RNA and DNA folding and hybridization software to the scientific community at large. By making use of universally available web GUIs (Graphical User Interfaces),the server circumvents the problem of portability of this software. Detailed output,in the form of structure plots with or without reliability information,single strand frequency plots and ‘energy dot plots’, are available for the folding of single sequences. A variety of ‘bulk’ servers give less information,but in a shorter time and for up to hundreds of sequences at once. The portal for the mfold web server is http://www.bioinfo.rpi.edu/applications/ mfold. This URL will be referred to as ‘MFOLDROOT’.

/pdf/mfold-web-server-for-nucleic-acid-folding-and-hybridization-3pzoe1h05m.pdf

Mfold web server for nucleic acid folding and hybridization prediction

https://njc.rockefeller.edu/pdf4/Class2-LeeAmbros1993.pdf

The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14

The structure and evolution of the human β-globin gene family

A 'rescuable' plasmid containing globin gene sequences allowing recombination with homologous chromosomal sequences has enabled us to produce, score and clone mammalian cells with the plasmid integrated into the human beta-globin locus. The planned modification was achieved in about one per thousand transformed cells whether or not the target gene was expressed.

Insertion of DNA sequences into the human chromosomal beta-globin locus by homologous recombination.

Argiris Efstratiadis Department of Biological Chemistry Harvard Medical School Boston, Massachusetts 02115 James W. Posakony, Tom Maniatis, Richard M. Lawn* and Catherine O’Connell+ Division of Biology California Institute of Technology Pasadena, California 91125 Richard A. Spritz, Jon K. DeRiel,# Bernard G. Forget and Sherman M. Weissman Departments of Genetics and Internal Medicine Yale University School of Medicine New Haven, Connecticut 06510 Jerry L. Slightom, Ann E. Blechl and Oliver Smithies Laboratory of Genetics University of Wisconsin Madison, Wisconsin 53706 Francisco E. Baralle, Carol C. Shoulders and Nicholas J. ProudfootQ MRC Laboratory of Molecular Biology Hills Road Cambridge CB2 2QH, England Summary We present the results of a detailed comparison of the primary structure of human p-like globin genes and their flanking sequences. Among the se- quences located 5’ to these genes are two highly conserved regions which include the sequences ATA and CCAAT located 31 2 1 and 77 + 10 bp, respectively, 5’ to the mRNA capping site. Similar sequences are found in the corresponding locations in most other eucaryotic structural genes. Calcula- tion of the divergence times of individual @like globin gene pairs provides the first description of the evolutionary relationships within a gene family based entirely on direct nucleotide sequence com- parisons. In addition, the evolutionary relationship of the embryonic e-globin gene to the other human P-like globin genes is defined for the first time. Finally, we describe a model for the involvement of short direct repeat sequences in the generation of deletions in the noncoding and coding regions of B-like globin genes during evolution.

https://www.cell.com/cell/pdf/0092-8674(80)90429-8.pdf

Oliver Smithies

Papers

A comprehensive set of sequence analysis programs for the VAX

The structure and evolution of the human β-globin gene family

Insertion of DNA sequences into the human chromosomal beta-globin locus by homologous recombination.

The Structure and Evolution of the Human /?-Globin Gene Family

Human fetal gγ- and Aγ-globin genes: Complete nucleotide sequences suggest that DNA can be exchanged between these duplicated genes