Nucleotide sequence of the bacteriophage T4 gene 57 and a deduced amino acid sequence.

A 693 basepair cloned fragment of bacteriophage T4 DNA, which supports specifically growth of T4 amber mutants in gene 57, has been sequenced. A polypeptide can be deduced from this sequence, that is either 54 or 60 amino acids long depending which of two AUG codons, 18 nucleotides apart, are used for initiation. The size of this deduced polypeptide is compatible with the size of a single polypeptide (based on polyacrylamide gel electrophoresis) synthesized in vivo in E. coli under the direction of the cloned T4 DNA fragment.     

mRNA sequence and deduced amino acid sequence of the mumps virus small hydrophobic protein gene.

The mRNA of a putative small hydrophobic protein (SH) of mumps virus was identified in mumps virus-infected Vero cells, and its complete nucleotide sequence was determined by sequencing the genomic RNA and cDNA clones and partial sequencing of mRNA. The SH mRNA is 310 nucleotides long excluding the poly(A) and contains a single open reading frame encoding a protein of 57 amino acids with a calculated molecular weight of 6,719. The predicted protein is highly hydrophobic and contains a stretch of 25 hydrophobic amino acids near the amino terminus which could act as a membrane anchor region. There is no homology between the putative SH protein of mumps virus and the SH protein of simian virus 5, even though the SH genes are located in the same locus in the corresponding genome. One interesting observation is that the hydrophobic domain of simian virus 5 SH protein is at the carboxyl terminus, whereas that of mumps virus putative SH protein is near the amino terminus.     

Prediction of protein secondary structure from amino acid sequence

The conformational parametersP _k for each amino acid species (j=1–20) of sequential peptides in proteins are presented as the product ofP _i,k , wherei is the number of the sequential residues in thekth conformational state (k=-helix,-sheet,-turn, or unordered structure). Since the average parameter for ann-residue segment is related to the average probability of finding the segment in the kth state, it becomes a geometric mean of (P _k )_av=(P _i,k ) ^1/n with amino acid residuei increasing from 1 ton. We then used ln(P_k)_av to convert a multiplicative process to a summation, i.e., ln(P _k ) _av =(1/n)P _i,k (i=1 ton) for ease of operation. However, this is unlike the popular Chou-Fasman algorithm, which has the flaw of using the arithmetic mean for relative probabilities. The Chou-Fasman algorithm happens to be close to our calculations in many cases mainly because the difference between theirP _k and our InP _k is nearly constant for about one-half of the 20 amino acids. When stronger conformation formers and breakers exist, the difference become larger and the prediction at the N- and C-terminal-helix or-sheet could differ. If the average conformational parameters of the overlapping segments of any two states are too close for a unique solution, our calculations could lead to a different prediction.     

Primary structure of rat pulmonary surfactant protein D. cDNA and deduced amino acid sequence.

Surfactant protein D (SP-D) is a carbohydrate-binding glycoprotein containing a collagen-like domain that is synthesized by alveolar type II epithelial cells. The complete primary structure of rat SP-D has been determined by sequencing of a cloned cDNA. The protein consists of three regions: an NH2-terminal segment of 25 amino acids, a collagen-like domain consisting of 59 Gly-X-Y repeats, and a COOH-terminal carbohydrate recognition domain of 153 amino acids. There are 6 cysteine residues present in rat SP-D: 2 in the NH2-terminal noncollagenous segment and 4 in the COOH-terminal carbohydrate-binding domain. The collagenous domain contains one possible N-glycosylation site. The protein is preceded by a cleaved, NH2-terminal signal peptide. SP-D shares considerable homology with the C-type mammalian lectins. Hybridization analysis demonstrates that rat SP-D is encoded by a 1.3-kilobase mRNA which is abundant in lung and highly enriched in alveolar type II cells. Extensive homology exists between rat SP-D and bovine conglutinin.     

Nucleotide and deduced amino acid sequence of stp: the bacteriophage T4 anticodon nuclease gene

Pre-existing host tRNAs are reprocessed during bacteriophage T4 infection of certain Escherichia coli strains. In this pathway, tRNALys is cleaved 5' to the wobble base by anticodon nuclease and is later restored in polynucleotide kinase and RNA ligase reactions. Anticodon nuclease depends on prr, a locus found only in host strains that restrict T4 mutants lacking polynucleotide kinase and RNA ligase; and on stp, the T4 suppressor of prr restriction. stp was cloned and the nucleotide sequences of its wild-type and mutant alleles determined. Their comparison defined an stp open reading frame of 29 codons at 162.8 to 9 kb of T4 DNA (1 kb = 10(3) base-pairs). We suggest that stp encodes a subunit of anticodon nuclease, perhaps one that harbors the catalytic site; while additional subunits, such as a putative prr gene product, impart protein folding environment and tRNA substrate recognition.     

Branched-chain amino acid aminotransferase of Escherichia coli: nucleotide sequence of the ilvE gene and the deduced amino acid sequence

The ilvE gene of the Escherichia coli K-12 ilvGEDA operon, which encodes branched-chain amino acid aminotransferase [EC 2.6.1.42], was cloned. The nucleotide sequence of 1.5 kilobase pairs containing the gene was determined. The coding region of the ilvE gene contained 927 nucleotide residues and could encode 309 amino acid residues. The predicted molecular weight, amino acid composition and the sequence of the N-terminal 15 residues agreed with the enzyme data reported previously (Lee-Peng, F.-C., et al. (1979) J. Bacteriol. 139, 339-345). From the deduced amino acid sequence, the secondary structure was predicted.     

Analysis of protein function and its prediction from amino acid sequence

Understanding protein function is one of the keys to understanding life at the molecular level. It is also important in the context of human disease because many conditions arise as a consequence of alterations of protein function. The recent availability of relatively inexpensive sequencing technology has resulted in thousands of complete or partially sequenced genomes with millions of functionally uncharacterized proteins. Such a large volume of data, combined with the lack of high-throughput experimental assays to functionally annotate proteins, attributes to the growing importance of automated function prediction. Here, we study proteins annotated by Gene Ontology (GO) terms and estimate the accuracy of functional transfer from protein sequence only. We find that the transfer of GO terms by pairwise sequence alignments is only moderately accurate, showing a surprisingly small influence of sequence identity (SID) in a broad range (30-100%). We developed and evaluated a new predictor of protein function, functional annotator (FANN), from amino acid sequence. The predictor exploits a multioutput neural network framework which is well suited to simultaneously modeling dependencies between functional terms. Experiments provide evidence that FANN-GO (predictor of GO terms; available from http://www.informatics.indiana.edu/predrag) outperforms standard methods such as transfer by global or local SID as well as GOtcha, a method that incorporates the structure of GO.     

The formation of protein secondary structure. Its connection with amino acid sequence

Statistical analysis of the occurrence of tetrapeptides in 35 globular proteins was performed. It was found that the amino acids along the polypeptide chain are close to being randomly distributed and that the same tetrapeptide segments exist in different types of secondary structure. Therefore, a new method was proposed for locating 'microdomains' in protein interiors. Amino acid replacements in the hydrophobic core of six proteins were analyzed. The results show that the locations of amino acids belonging to defined microdomains are extremely conserved. It is suggested that the structures found may play a role as nucleation centers in protein folding.     

Combining the GOR V algorithm with evolutionary information for protein secondary structure prediction from amino acid sequence

We have modified and improved the GOR algorithm for the protein secondary structure prediction by using the evolutionary information provided by multiple sequence alignments, adding triplet statistics, and optimizing various parameters. We have expanded the database used to include the 513 non-redundant domains collected recently by Cuff and Barton (Proteins 1999;34:508-519; Proteins 2000;40:502-511). We have introduced a variable size window that allowed us to include sequences as short as 20-30 residues. A significant improvement over the previous versions of GOR algorithm was obtained by combining the PSI-BLAST multiple sequence alignments with the GOR method. The new algorithm will form the basis for the future GOR V release on an online prediction server. The average accuracy of the prediction of secondary structure with multiple sequence alignment and full jack-knife procedure was 73.5%. The accuracy of the prediction increases to 74.2% by limiting the prediction to 375 (of 513) sequences having at least 50 PSI-BLAST alignments. The average accuracy of the prediction of the new improved program without using multiple sequence alignments was 67.5%. This is approximately a 3% improvement over the preceding GOR IV algorithm (Garnier J, Gibrat JF, Robson B. Methods Enzymol 1996;266:540-553; Kloczkowski A, Ting K-L, Jernigan RL, Garnier J. Polymer 2002;43:441-449). We have discussed alternatives to the segment overlap (Sov) coefficient proposed by Zemla et al. (Proteins 1999;34:220-223).     

The primary structure of phosphoenolpyruvate carboxylase of Escherichia coli. Nucleotide sequence of the ppc gene and deduced amino acid sequence

The nucleotide sequence of the ppc gene, the structural gene for phosphoenolpyruvate carboxylase [EC 4.1.1.31], of Escherichia coli K-12 was determined. The gene codes for a polypeptide comprising 883 amino acid residues with a calculated molecular weight of 99,061. The amino acid sequence deduced from the nucleotide sequence was entirely consistent with the protein chemical data obtained with the purified enzyme, including the NH2- and COOH-terminal sequences and amino acid composition. The coding region is preceded by two putative ribosome binding sites, and is followed closely by a good representative of rho-independent terminator. The codon usage in the ppc gene suggests a moderate expression of the gene. The secondary structure of the enzyme was predicted from the deduced amino acid sequence.     

Nucleotide sequence of the murine prothymosin alpha cDNA and its deduced primary and secondary protein structure

Searching for proliferation-related and cell cycle phase-specific genes we detected a full-length cDNA for the murine prothymosin alpha mRNA which was sequenced on the DNA level. The amino acid sequence deduced from the nucleotide sequence shows a high degree of positional identities with prothymosin alpha from man and rat. However, the minor differences in the primary structures largely influence predictions for the secondary structures of prothymosin alpha from different species. These differences in the secondary structure could explain the differences of activity of prothymosin alpha from different origin in immuno-protection assays.     

Putative amino acid sequence of chick calcium-binding protein deduced from a complementary DNA sequence.

Two DNA fragments coding for chick CaBP have been isolated and sequenced. cDNA was prepared from enriched intestinal mRNA and cloned in pUC12. The recombinant clones were screened by differential hybridisation with 32P-cDNA probes synthesized from vitamin D replete and deficient chick intestinal mRNA. Two clones had outstanding affinity with the +D probe. Hybrid-arrested and hybrid-selected translation systems showed that both clones hybridised to mRNA coding for immunoprecipitable CaBP. The mRNA for CaBP has a 100 bp G,C rich sequence before a 786 bp coding region followed by 1250 nucleotides 3' untranslated region. Nucleotides coding for the Ca-binding sites show a high degree of homology for Ca-binding sites in chick calmodulin and rat intestinal CaBP. The amino acid sequence specified by the longest open reading frame contains five Ca-binding sites but is too large for the native CaBP; post-translational modification must therefore occur.     

Rotavirus neutralizing protein VP7: antigenic determinants investigated by sequence analysis and peptide synthesis.

The rotavirus neutralizing antigen, VP7, is a 37,000-molecular-weight glycoprotein which is a major component of the outer shell of the virion. The amino acid sequence of VP7 for strain S2 (human serotype 2) and Nebraska calf diarrhea virus (bovine serotype) has been inferred from the nucleic acid sequence of cloned copies of genomic segment nine. Comparison of the amino acid sequences of these two VP7 proteins with those already determined for other rotavirus strains reveals extensive sequence conservation between serotypes with clusters of amino acid differences sited predominantly in hydrophilic domains of the protein. Six peptides have been synthesized that span the hydrophilic regions of the molecule. Antisera to these peptides both recognize the respective homologous peptides in a solid-phase radioimmunoassay and bind to denatured VP7 in a Western blot. However, none of the antisera either recognize virus or exhibit significant neutralizing activity, indicating that these peptide sequences are not available on the surface of the virus.