Computer programs for the analysis and the management of DNA sequences

Prolyl isomerases do not catalyze isomerization of non-prolyl peptide bonds

Prolyl isomerases accelerate the cis <--> trans isomerization of prolyl peptide bonds during protein folding and probably also in folded proteins. We asked whether this catalytic function is in fact restricted to prolyl bonds or whether the isomerizations of 'normal' non-prolyl peptide bonds are catalyzed as well. By using the P39A variant of ribonuclease T1 as a substrate we find that the trans --> cis isomerization of the Tyr38-Ala39 bond in the refolding of this protein is not catalyzed by prolyl isomerases of the cyclophilin, FKBP and parvulin families. These enzymes are neither able to catalyze amide bond isomerizations in the proline-free model peptide Ala-Ala-Tyr-Ala-Ala.

X-HUSAR, an X-based graphical interface for the analysis of genomic sequences

Management and analysis of nucleotide and protein sequence and structure data constitute a traditional area of bioinformatics. Since the analytical programs are frequently developed by researchers, rather than software engineers, they tend to suffer from idiosyncratic and non-ergonomic man-machine interfaces. We report on HUSAR, our 140+ collection of third-party, as well as in-house developed or adapted, sequence manipulation and analysis tools, well integrated into the UNIX operating system environment and accessible via consistent menu-aware interface. Most of the HUSAR programs can be completely specified by UNIX command-line options; they can thus be run in batches or combined into pipes. Adding such a program into the HUSAR environment is almost a 'plug-and-play' exercise. HUSAR has been recently complemented with a graphical client interface, X-HUSAR, to support users on UNIX platforms with X11 windowing systems. The whole X-HUSAR interface is based on a single generic program, COMLIGEN, and a number of specific configuration files. COMLIGEN interprets those files and renders appropriate windows, menus, and other interactive elements, which help the end user in selecting application programs and specifying their options. Efforts of extending both HUSAR and X-HUSAR are roughly linear to the size of the collection.

The humoral response to heat shock protein 70 in human and murine Schistosomiasis mansoni

We have expressed two cDNA sequences encoding 121 and 230 amino acids of the C-terminus of the Schistosoma mansoni Hsp 70 in Escherichia coli. The products were synthesized as polypeptides fused to the RNA polymerase of bacteriophage MS2, and their reactivities were tested in ELISAs, using sera from human and murine infections. Anti-Hsp70 antibodies were detected in a significant number of individuals suffering from chronic schistosomiasis mansoni, but not in patients with known recent infections. This, together with the finding that antibodies directed at S. mansoni-specific Hsp70 determinants during the course of infection of experimental mice were not detectable until 5-6 weeks post-infection, suggests that the protein may be a useful marker for distinguishing late and early infections. The diagnostic specificity of Hsp70 was evaluated with sera from humans infected with different schistosome species and other parasitic diseases. While some subjects infected with S. haematobium produced antibodies which recognized the S. mansoni Hsp70, no such antibodies were generated in S. japonicum infected individuals. However, cross-reactive antibodies were elicited in donors with other parasitic diseases such as filariasis and malaria. The absence of antibodies in early infection and the observed cross-reactivities led us to conclude that Hsp70 will be of limited value in the diagnosis of schistosomiasis.

Primary structures of Sm31/32 diagnostic proteins of Schistosoma mansoni and their identification as proteases

We have constructed cDNA clones containing the complete nucleotide sequences coding for two highly antigenic Schistosoma mansoni adult worm proteins, Sm31 and Sm32. The predicted amino acid sequence of Sm31 shows significant homology to mouse, rat and human cathepsin B. The nucleotide sequence of Sm32 is identical to that reported by others for S. mansoni "haemoglobinase'. The different nucleotide sequences demonstrate the existence of two different proteolytic enzymes, both of which are synthesised in the form of precursor molecules. Structural homology of the schistosome cathepsin B to the mammalian ones indicates that the mature protein is processed from a propeptide. The calculated molecular weight of haemoglobinase of 47,000 suggests that post-translational processing is also involved in generating an active protease.

Subtyping of European foot-and-mouth disease virus strains by nucleotide sequence determination

The VP1-coding regions of foot-and-mouth disease virus strains from 18 recent European outbreaks and of 9 strains isolated more than 20 years ago and used in part as vaccines were determined by direct cDNA sequencing. Comparison of the sequences revealed that most of the isolated outbreak viruses are closely related to the vaccine strains used. Isolates from the Italian epizootic of 1984 to 1985 correspond, for example, to the vaccine strain A5 Parma 62; the outbreak in 1984 in Bernbeuren, Federal Republic of Germany, was induced by A5 Allier 60; outbreaks in 1982 in Funen, Denmark, and in Murchin, German Democratic Republic, were caused by O1 Lausanne 65. Viruses isolated during the 1983 Iberian epizootic show a close relationship to the vaccine strain A5 Allier 60 but were probably derived from another not yet identified vaccine strain from Spain. Only two minor outbreaks in the Federal Republic of Germany, A Aachen in 1976 and O Wuppertal in 1982, did not correspond to the classical European strains but were obviously introduced from outside. We suggest that nucleotide sequence analysis should be used as a standard method of diagnosis, because when compared with other techniques it more clearly reveals the origin and course of epizootics and offers the possibility of preventing further outbreaks.

Nonepidermal members of the keratin multigene family: cDNA sequences and in situ localization of the mRNAs

A keratin set which consists of a type I 47kd and a type II 57kd protein occurs as a major constituent of the keratin patterns of various internal stratified epithelia of the mouse. We have isolated specific cDNA clones of the two complementary keratin subunits from a cDNA library constructed with polyA+RNA of mouse tongue epithelium and present the complete nucleotide and deduced amino acid sequences of the 57kd protein and about 75% of the corresponding data of the 47kd protein. The comparison of the sequence data with those of known epidermal keratin mRNAs coding for the two types of keratin proteins reveals a fundamentally identical and type-specific organization of the mRNAs into both highly conserved and variable domains. In order to avoid cross-reactions with other members of the keratin multigene family, appropriately taylored 35S-labeled cDNA probes comprising the low and non-homologous 3' coding and noncoding domains of the mRNAs were used for in situ hybridization to tissue sections. The localization and distribution of the corresponding transcripts indicates a strongly compartmentalized keratin expression in mouse tongue epithelium.

A peptide to DNA conversion program

A modification and extension of the computer program REVCUT (Blumenthal et al, Nucl. Acids Res. 10, 91-101 (1982) is described. The new program searches for restriction endonuclease recognition sites that are not coding DNA sequences of a protein of known aminoacid sequence using bit patterns. The modifications make the program more accurate and extend the range of the restriction endonucleases.

Complete sequence of a bovine type I cytokeratin gene: conserved and variable intron positions in genes of polypeptides of the same cytokeratin subfamily

The complete sequence of a bovine gene encoding an epidermal cytokeratin of mol. wt. 54 500 (No VIb) of the acidic (type I) subfamily is presented, including an extended 5' upstream region. The gene (4377 bp, seven introns) which codes for a representative of the glycine-rich subtype of cytokeratins of this subfamily, is compared with genes coding for: another subtype of type I cytokeratin; a basic (type II) cytokeratin gene; and vimentin, a representative of another intermediate filament (IF) protein class. The positions of the five introns located within the highly homologous alpha-helix-rich rod domain are identical or equivalent, i.e., within the same triplet, in the two cytokeratin I genes. Four of these intron positions are also identical with intron sites in the vimentin gene, and three of these intron positions are identical or similar in the type I and type II cytokeratin subfamilies. On the other hand, the gene organization of both type I cytokeratins differs from that of the type II cytokeratin in the rod region in five intron positions and in the introns located in the carboxy-terminal tail region, with the exception of one position at the rod-tail junction. Remarkably, the two type I cytokeratins also differ from each other in the positions of two introns located at and in the region coding for the hypervariable, carboxy-terminal portion. The introns and the 5' upstream regions of the cytokeratin VIb gene do not display notable sequence homologies with the other IF protein genes, but sequences identical with--or very similar to--certain viral and immunoglobulin enhancers have been identified.(ABSTRACT TRUNCATED AT 250 WORDS)

Amino acid sequence microheterogeneities of basic (type II) cytokeratins of Xenopus laevis epidermis and evolutionary conservativity of helical and non-helical domains

Three clones coding for carboxy-terminal portions of type II cytokeratins have been isolated from a cDNA library constructed from the epidermis of the frog Xenopus laevis. These clones have been identified by hybridization-selection-translation and Northern blot analysis, and contain sequences complementary to mRNAs of similar size that code for three different polypeptides of the Mr 64,000 group, Ia-c, i.e. the only major type II cytokeratins expressed in this tissue. A comparison of the corresponding nucleotide sequences and the amino acid sequences deduced therefrom shows only minor differences in these polypeptides, most of which occur as isolated point mutations. This indicates that coding sequences of the different type II cytokeratin genes in epidermis of Xenopus are very similar, in contrast to the more extended differences of type II cytokeratin genes expressed in mammalian epidermis, which probably reflects a lower degree of evolutionary divergence of members of this protein family in amphibia. A comparison of the Xenopus sequences with those of mammalian type II cytokeratins reveals the same characteristic features, i.e. an alpha-helical domain ending with the familiar consensus sequence T Y R (X Y) L E G E, followed by a non-helical domain Cl enriched in hydroxyamino acids. Both domains are remarkably conserved in sequence between Xenopus and mammals. The following glycine-rich domain (C2) displays similar oligopeptide repeats (mostly of the type G G G M in the frog keratins), and the terminal C3 domain is characterized by a region exceptionally rich in hydroxyamino acids, which immediately precedes a cluster of basic amino acids at the carboxy terminus. Our results show that the typical features of the domain of type II cytokeratins are already established in amphibia and that these homologies are not restricted to the alpha-helical rod of these proteins but, in principle, extend to the other domains located in the so-called hypervariable tail portion. This suggests that the hypervariable regions are not subject to random variability but contain functionally important domains that have been well conserved during evolution.

Comparative sequence analysis of duck and human hepatitis B virus genomes

We have cloned and sequenced an infectious, functionally active genome of a duck hepatitis B virus (DHBV). It is 3,021 base pairs (bp) in length and shows little DNA sequence homology to the genome of human hepatitis B virus (HBV). However, the amino acid sequences of predicted viral gene products are similar between DHBV and HBV, and the genome organization present in DHBV reflects that of HBV. As in the mammalian virus the long minus strand of the DHBV genome encodes three long overlapping reading frames designated as P, S, and C. The fourth open reading frame, termed X, is absent in DHBV. A comparison with a sequence of a second DHBV isolate [Mandart et al, Journal of Virology 49:782-792, 1984] revealed a nucleotide sequence variation of 5.6% and confirmed the presented overall gene organization of DHBV.

Identification of two types of keratin polypeptides within the acidic cytokeratin subfamily I

Cytoskeletal filaments of the alpha-keratin type (cytokeratins) are a characteristic of epithelial cells. In diverse mammals (man, cow and rodents) these cytokeratins consist of a family of approximately 20 polypeptides, which may be divided into the more acidic (I) and the more basic (II) subfamilies. These two subfamilies show only limited amino acid sequence homology. In contrast, nucleic acid hybridization experiments and peptide maps have been interpreted to show that polypeptides of the same subfamily share extended sequence homology. We compare two polypeptides of the acidic cytokeratin subfamily, VIb (Mr 54,000) and VII (Mr 50,000), which are co-expressed in large amounts in bovine epidermal keratinocytes. These two epidermal keratins can be distinguished by specific antibodies and show different patterns of expression among several bovine tissues and cultured cells. In addition, they differ in the stability of their complexes with basic keratin polypeptides and in their tryptic peptide maps. The amino acid sequences deduced from the nucleotide sequences of complementary DNA clones containing the 3' ends of the messenger RNAs for these keratins are compared with each other and with available amino acid sequences of human, murine and amphibian epidermal keratins. Bovine keratins VIb and VII share considerable sequence homology in the alpha-helical portion (68% residues identical) but lack significant homology in the extrahelical portion. Bovine keratin VIb shows, in its alpha-helical region, a pronounced sequence homology (88% identity) to the murine epidermal keratin of Mr 59,000. In addition, the non-helical carboxy-terminal regions of both proteins are glycine-rich and contain a canonic sequence GGGSGYGG, which may be repeated several times. Moreover, their mRNAs present a highly conserved stretch of 236 nucleotides containing, in the murine sequence, the end of the coding and all of the non-coding region (81% identical nucleotides). Bovine keratin VII is considerably different from the murine Mr 59,000 keratin but is almost identical to the human cytokeratin number 14 of Mr 50,000, both in the alpha-helical and in the non-alpha-helical regions of the proteins, and the mRNAs of the human and the bovine keratins also display a high homology in their 3' non-coding ends. The results show that in the same species keratins of the same subfamily can differ considerably, whereas equivalent keratin polypeptides of different species are readily identified by characteristic sequence homologies in the alpha-helical and the non-helical regions as well as in the 3' non-coding portions of their mRNAs.(ABSTRACT TRUNCATED AT 400 WORDS)

Nucleotide sequence and genome organization of foot-and-mouth disease virus

A continuous 7802 nucleotide sequence spanning the 94% of foot and mouth disease virus RNA between the 5'-proximal poly(C) tract and the 3'-terminal poly(A) was obtained from cloned cDNA, and the total size of the RNA genome was corrected to 8450 nucleotides. A long open reading frame was identified within this sequence starting about 1300 bases from the 5' end of the RNA genome and extending to a termination codon 92 bases from its polyadenylated 3' end. The protein sequence of 2332 amino acids deduced from this coding sequence was correlated with the 260 K FMDV polyprotein. Its processing sites and twelve mature viral proteins were inferred from protein data, available for some proteins, a predicted cleavage specificity of an FMDV encoded protease for Glu/Gly(Thr, Ser) linkages, and homologies to related proteins from poliovirus. In addition, a short unlinked reading frame of 92 codons has been identified by sequence homology to the polyprotein initiation signal and by in vitro translation studies.

Structure of the FMDV translation initiation site and of the structural proteins

A cDNA clone of Foot and Mouth Diseases Virus (FMDV), strain C1, has been sequenced. The limits of the structural genes were defined by comparison with the available protein data. We identified two potential translation initiation sites for the viral polyprotein separated by 84 nucleotides. We suggest that these two initiation sites could be used to express two proteins differing only at the N-terminal, P16 and P20a. This model is supported by the fact that antiserum against a bacterially synthesized polypeptide corresponding to the anterior region of the polyprotein precipitates specifically both P16 and P20a. Comparison of the C1 sequence with two other serotypes, O1K and A10 revealed variability in the major immunogenic structural protein, VP1, and also in two other capsid proteins, VP2 and VP3. P16/P20a, VP4, and the N-terminal part of the precursor of the nonstructural genes, P52, are rather conserved between the different FMDV strains.

Cross-reaction of snRNA and an Alu I-like sequence from rat with DNAs from different eucaryotic species

Sequence homologies to rat U1-snRNA and U2-snRNA were investigated in DNA from 23 eucaryotic species (3 lower eucaryotes, 4 plants, and 16 animals) using dot hybridization at various stringency conditions. Cross-hybridization among very distantly related species in e.g. plants-insects or mold-vertebrates is not the rule; there are, however, examples for stronger homologies like Rattus-Dictyostelium. Furthermore, DNA from all 23 species was analysed for sequence homologies with the repetitive DNA sequence B1 (an Alu I family equivalent) from rat. We observed a wide range of homologies covering some plants and insects, up to vertebrates. Hybridization at increasing stringency conditions revealed species with higher degree of homology to the rat B1 sequence: maize, chicken, and hamster.

A tandem repeat gene in a picornavirus

Three closely related genes for the small genome-linked protein (VPg) of picornaviruses have been identified by sequence analysis as a tandem repeat in the genome of Foot and Mouth Disease Virus (FMDV), strain O1K. This unusual structure was also found in the genome of strain C1O, belonging to a different FMDV serotype. Predicted biochemical properties of the three VPg gene products are in excellent agreement with the data from protein analysis of a heterogeneous VPg population from a third FMDV serotype, strain A10 (1). Taken together, these data indicate that the VPgs from all three genes function equally well in vivo. This is the first report of a tandem repeat gene in a viral genome.

Nucleotide sequence and exact localization of the neomycin phosphotransferase gene from transposon Tn5

The nucleotide sequence of 1200 bp from the unique region of transposon Tn5 containing the neomycin phosphotransferase gene (neo) was determined, and the location of the neo gene was identified by deletion mutants in a translational reading frame of 792 bp. The derived gene product, an aminoglycoside 3'-phosphotransferase (APH) II, consists of 264 amino acid residues and has a calculated Mr of 29053. Its amino acid sequence shows sequence homologies to the APH type I enzyme coded for by transposon Tn903 (Oka et al., 1981).