Evolution and taxonomy of positive-strand RNA viruses: implications of comparative analysis of amino acid sequences

Despite the rapid mutational change that is typical of positive-strand RNA viruses, enzymes mediating the replication and expression of virus genomes contain arrays of conserved sequence motifs. Proteins with such motifs include RNA-dependent RNA polymerase, putative RNA helicase, chymotrypsin-like and papain-like proteases, and methyltransferases. The genes for these proteins form partially conserved modules in large subsets of viruses. A concept of the virus genome as a relatively evolutionarily stable "core" of housekeeping genes accompanied by a much more flexible "shell" consisting mostly of genes coding for virion components and various accessory proteins is discussed. Shuffling of the "shell" genes including genome reorganization and recombination between remote groups of viruses is considered to be one of the major factors of virus evolution. Multiple alignments for the conserved viral proteins were constructed and used to generate the respective phylogenetic trees. Based primarily on the tentative phylogeny for the RNA-dependent RNA polymerase, which is the only universally conserved protein of positive-strand RNA viruses, three large classes of viruses, each consisting of distinct smaller divisions, were delineated. A strong correlation was observed between this grouping and the tentative phylogenies for the other conserved proteins as well as the arrangement of genes encoding these proteins in the virus genome. A comparable correlation with the polymerase phylogeny was not found for genes encoding virion components or for genome expression strategies. It is surmised that several types of arrangement of the "shell" genes as well as basic mechanisms of expression could have evolved independently in different evolutionary lineages. The grouping revealed by phylogenetic analysis may provide the basis for revision of virus classification, and phylogenetic taxonomy of positive-strand RNA viruses is outlined. Some of the phylogenetically derived divisions of positive-strand RNA viruses also include double-stranded RNA viruses, indicating that in certain cases the type of genome nucleic acid may not be a reliable taxonomic criterion for viruses. Hypothetical evolutionary scenarios for positive-strand RNA viruses are proposed. It is hypothesized that all positive-strand RNA viruses and some related double-stranded RNA viruses could have evolved from a common ancestor virus that contained genes for RNA-dependent RNA polymerase, a chymotrypsin-related protease that also functioned as the capsid protein, and possibly an RNA helicase.     

Designing amino acid sequences to fold with good hydrophobic cores

We present two methods for designing amino acid sequences of proteins that will fold to have good hydrophobic cores. Given the coordinates of the desired target protein or polymer structure, the methods generate sequences of hydrophobic (H) and polar (P) monomers that are intended to fold to these structures. One method designs hydrophobic inside, polar outside; the other minimizes an energy function in a sequence evolution process. The sequences generated by these methods agree at the level of 60-80% of the sequence positions in 20 proteins in the Protein Data Bank. A major challenge in protein design is to create sequences that can fold uniquely, i.e. to a single conformation rather than to many. While an earlier lattice-based sequence evolution method was shown not to design unique folders, our method generates unique folders in lattice model tests. These methods may also be useful in designing other types of foldable polymer not based on amino acids.     

Representation of amino acid sequences as two-dimensional point patterns

An algorithm for the representation of amino acid sequences as two-dimensional point patterns (2-D plot) is described. The algorithm is based on chaos game representation (CGR) for DNA sequences and was extended for amino acid sequences. The 2-D plot depicts the sequentiality of amino acids and the amino acid composition of a protein. Changes in a protein sequence as insertion, deletion and repeats of amino acids are characterized by specific geometrical properties and changes in the 2-D plots. The 2-D plot may be considered as a two-dimensional "fingerprint" of a protein. The properties of the algorithm are explained by user-defined amino acid sequences. As an example the 2-D plots of two selected heart proteins are generated. The sequences of these proteins are obtained from the protein sequence database SWISS-PROT.     

Analysis of amino acid sequences of Rhodobacter sphaeroides glutamine synthetase

A comparative analysis of the amino acid sequence of glutamine synthetase (GS) of the photosynthetic purple bacterium Rhodobacter sphaeroides revealed that the enzyme is typical for first type procaryotic GSs and structurally resembles GSs of enteric bacteria. The data obtained indicate that the complex phenotype of purple bacterial mutants at the glnA gene coding for GS may be conditioned by specific regulation of nitrogen metabolism in bacterial cells rather than by structural-and-functional peculiarities of GS.     

Assigning amino acid sequences to 3-dimensional protein folds

With the advent of genome sequencing projects, the amino acid sequences of thousands of proteins are determined every year. Each of these protein sequences must be identified with its function and its 3-dimensional structure for us to gain a full understanding of the molecular biology of organisms. To meet this challenge, new methods are being developed for fold recognition, the computational assignment of newly determined amino acid sequences to 3-dimensional protein structures. These methods start with a library of known 3-dimensional target protein structures. The new probe sequence is then aligned to each target protein structure in the library and the compatibility of the sequence for that structure is scored. If a target structure is found to have a significantly high compatibility score, it is assumed that the probe sequence folds in much the same way as the target structure. The fundamental assumptions of this approach are that many different sequences fold in similar ways and there is a relatively high probability that a new sequence possesses a previously observed fold. We review various approaches to fold recognition and break down the process into its main steps: creation of a library of target folds; representation of the folds; alignment of the probe sequence to a target fold using a sequence-to-structure compatibility scoring function; and assessment of significance of compatibility. We emphasize that even though this new field of fold recognition has made rapid progress, technical problems remain to be solved in most of the steps. Standard benchmarks may help identify the problem steps and find solutions to the problems.     

The comparative amino acid sequences, substrate specificities and gene or cDNA nucleotide sequences of some prokaryote and eukaryote amidinotransferases: implications for evolution

The amino acid sequences of the amidinotransferases and the nucleotide sequences of their genes or cDNA from four Streptomyces species (seven genes) and from the kidneys of rat, pig, human and human pancreas were compared. The overall amino acid and nucleotide sequences of the prokaryotes and eukaryotes were very similar and further, three regions were identified that were highly identical. Evidence is presented that there is virtually zero chance that the overall and high identity regions of the amino acid sequence similarities and the overall nucleotide sequence similarities between Streptomyces and mammals represent random match. Both rat and lamprey amidinotransferases were able to use inosamine phosphate, the amidine group acceptor of Streptomyces. We have concluded that the structure and function of the amidinotransferases and their genes has been highly conserved through evolution from prokaryotes to eukaryotes. The evolution has occurred with: (1) a high degree of retention of nucleotide and amino acid sequences; (2) a high degree of retention of the primitive Streptomyces guanine + cytosine (G + C) third codon position composition in certain high identity regions of the eukaryote cDNA; (3) a decrease in the specificities for the amidine group acceptors; and (4) most of the mutations silent in the regions suggested to code for active sites in the enzymes.     

A novel derivatization method with 5-bromonicotinic acid N-hydroxysuccinimide for determination of the amino acid sequences of peptides

We have developed a novel method that effectively identifies the N-terminal product ions produced in the tandem mass spectrometry (MS/MS) analysis of peptides done in conjunction with the specific derivatization of the N-terminal amino group using 5-bromonicotinic acid N-hydroxysuccinimide ester (BrNA-NHS). Electrospray ionization with low-energy collision-induced dissociation (CID) MS/MS clearly differentiated the N-terminal product ions labeled with the 5-bromonicotinyl group from other ions, on the basis of the appearance of CID peaks with a doublet pattern characteristically separated by 2 mass units produced by the equal natural abundances of 79Br and 81Br. The tracing of a series of these bromine-containing product ions allows the easy amino acid sequencing of peptides. Using Gln-Arg-Leu-Gln-Ser-Asn-Gln-Leu-Lys as the test peptide, we found that within 30 minutes at pH 6.5 and 37 degrees C its alpha-amino group was completely acylated with BrNA-NHS (peptide: BrNA-NHS = 1:40; mol/mol). The epsilon-amino group of the C-terminal lysine residue was less likely to be acylated under these conditions, being only partly modified (about 20%). This suggests the possibility of keeping the epsilon-amino group free from acylation. The method was successfully applied to the determination of the amino acid sequences of peptides from porcine kidney aminoacylase I produced by digestion with lysyl endopeptidase and with Staphylococus aureus V8 protease.     

A new method of inference of ancestral nucleotide and amino acid sequences

A statistical method was developed for reconstructing the nucleotide or amino acid sequences of extinct ancestors, given the phylogeny and sequences of the extant species. A model of nucleotide or amino acid substitution was employed to analyze data of the present-day sequences, and maximum likelihood estimates of parameters such as branch lengths were used to compare the posterior probabilities of assignments of character states (nucleotides or amino acids) to interior nodes of the tree; the assignment having the highest probability was the best reconstruction at the site. The lysozyme c sequences of six mammals were analyzed by using the likelihood and parsimony methods. The new likelihood-based method was found to be superior to the parsimony method. The probability that the amino acids for all interior nodes at a site reconstructed by the new method are correct was calculated to be 0.91, 0.86, and 0.73 for all, variable, and parsimony-informative sites, respectively, whereas the corresponding probabilities for the parsimony method were 0.84, 0.76, and 0.51, respectively. The probability that an amino acid in an ancestral sequence is correctly reconstructed by the likelihood analysis ranged from 91.3 to 98.7% for the four ancestral sequences.     

Complete amino acid sequences and phosphorylation sites, determined by Edman degradation and mass spectrometry, of rat parotid destrin- and cofilin-like proteins

Beta-adrenergic or cholinergic stimulation of the rat parotid gland was earlier shown to induce dephosphorylation of endogenous destrin- and cofilin-like proteins, which are phosphorylated in resting cells at Ser residues probably present near the N-terminals. The primary structures and phosphorylation sites were determined here. The rat destrin-like protein had a sequence 95% identical to the cDNA-derived sequence of porcine destrin. The rat cofilin-like protein was 98% identical to that of porcine cofilin. Each protein lacked the initiator Met and began with an acetylalanine residue followed by a Ser residue. The N-terminal peptides generated with endoproteinase Asp-N were isolated; they were each phosphorylated at Ser-2. Earlier work had shown that partial cleavage of the phosphorylated destrin- and cofilin-like proteins with cyanogen bromide provides unphosphorylated 16.7- and 18.3-kDa fragments, respectively. It was here confirmed that they contained all the Ser residues other than those present in the N-terminal peptides. From these observations, it was now concluded that the destrin- and cofilin-like proteins are rat parotid destrin and cofilin (non-muscle type), respectively, and that each protein is phosphorylated exclusively at Ser-2 in resting cells and dephosphorylated at this site in response to beta-adrenergic or cholinergic stimulation.     

The evolution of proteins from random amino acid sequences. I. Evidence from the lengthwise distribution of amino acids in modern protein sequences

We examine in this paper one of the expected consequences of the hypothesis that modern proteins evolved from random heteropeptide sequences. Specifically, we investigate the lengthwise distributions of amino acids in a set of 1,789 protein sequences with little sequence identify using the run test statistic (ro) of Mood (1940, Ann. Math. Stat. 11, 367-392). The probability density of ro for a collection of random sequences has mean = 0 and variance = 1 [the N(0,1) distribution] and can be used to measure the tendency of amino acids of a given type to cluster together in a sequence relative to that of a random sequence. We implement the run test using binary representations of protein sequences in which the amino acids of interest are assigned a value of 1 and all others a value of 0. We consider individual amino acids and sets of various combinations of them based upon hydrophobicity (4 sets), charge (3 sets), volume (4 sets), and secondary structure propensity (3 sets). We find that any sequence chosen randomly has a 90% or greater chance of having a lengthwise distribution of amino acids that is indistinguishable from the random expectation regardless of amino acid type. We regard this as strong support for the random-origin hypothesis. However, we do observe significant deviations from the random expectation as might be expected after billions years of evolution. Two important global trends are found: (1) Amino acids with a strong alpha-helix propensity show a strong tendency to cluster whereas those with beta-sheet or reverse-turn propensity do not. (2) Clustered rather than evenly distributed patterns tend to be preferred by the individual amino acids and this is particularly so for methionine. Finally, we consider the problem of reconciling the random nature of protein sequences with structurally meaningful periodic "patterns" that can be detected by sliding-window, autocorrelation, and Fourier analyses. Two examples, rhodopsin and bacteriorhodopsin, show that such patterns are a natural feature of random sequences.     

Muscimol potentiation of acidic amino acid release from cerebellar synaptosomes is chloride dependent

In previous studies we have shown that the depolarization-induced release of preaccumulated acidic amino acids and newly synthesized glutamate from cerebellar synaptosomal preparations is potentiated by gamma-aminobutyric acid (GABA) agonists through a GABAergic presynaptic mechanism. Here we report a systematic analysis of the ionic requirements of the potentiating effect of muscimol on the high K+-evoked release of D-[3H]aspartate. Our studies show that: Ca2+, Na+, and Mg2+ are not required for muscimol to exert its effect; a depolarizing concentration of K+ is a necessary, but not sufficient, condition to observe the presynaptic effect in question; and a minimal Cl- concentration (50-70 mM) is also required. A possible model based on these findings is proposed.     

Interaction of a bZip oligopeptide model with oligodeoxyribonucleotides modelling DNA binding sites. The effect of flanking sequences

A leucine zipper (bZip) binding peptide BP1 was constructed based on the DNA binding sequence of the GCN4 protein, slightly modified to make it more similar to the sequence of other bZip proteins (Jun) with related DNA binding specificity. Self-complementary DNA hexadecanucleotides containing ATF/CRE, AP-1 and C/EPB target sites were used to study peptide-DNA complex formation. Conformation changes in both components that occur on complex formation were studied by circular dichroism (CD) spectroscopy. The results show that the amount of alpha-helix formed in the peptide strongly depends not only on the target site present, but also on the type of the sequence flanking the ATF/CRE target site. Highest amount of the alpha-helix induced in the peptide was observed when homopurine homopyrimidine flanking sequences were present, whereas the presence of alternating sequences, especially of the CA/TG type, showed considerably lower effects. The change in DNA conformation on complex formation was generally small, but also depended on the type of the flanking sequence. It appears that the sequences flanking the target site can considerably modify the ability of the target sequence to bind specifically the bZip peptide, probably by slightly varying the overall DNA conformation.     

Single-molecule detection of specific nucleic acid sequences in unamplified genomic DNA

A new technique is described for the rapid detection of specific nucleic acid sequences in unamplified DNA samples. The method consists of using two nucleic acid probes complementary to different sites on a target DNA sequence. The two probes are each labeled with different fluorescent dyes. When mixed with a sample containing the target DNA, the two probes hybridize to their respective binding sites on the same target DNA molecule. The sample is then analyzed by a laser-based ultrasensitive fluorescence system capable of detecting single fluorescent molecules at two different wavelength channels simultaneously. Since the probes are bound to the same target DNA molecule, their signals appear simultaneously. Thus, coincident detection of both dyes provides the necessary specificity to detect an unamplified, single-copy target DNA molecule in a homogeneous assay. If the target is not present, only uncorrelated events originating from free probes will be observed at either channel. Phage lambda DNA in a background of salmon genomic DNA was detected as a two-dye coincident signal at a relative concentration of one lambda molecule per salmon genome. In a control sample, cleavage of the lambda DNA between the two probe binding sites eliminated the coincident signals. In a second experiment, a single-copy transgene was detected in maize. Detection parameters and possible future applications to genetic analysis are discussed.     

Rice embryo globulins: amino-terminal amino acid sequences, cDNA cloning and expression

A globulin fraction prepared from rice embryos contained polypeptides or polypeptide groups of 49 kDa (designated REG1), 46 kDa (designated REG2), about 35 kDa, 32 kDa and 25 kDa. The amino-terminal sequences of REG1 and the major polypeptide in the 35-kDa group were identical, suggesting that the REG1 polypeptide undergoes partial proteolytic processing that removes a carboxy-terminal region. A cDNA clone, designated pcREG2, encoding REG2 was isolated, and its nucleotide sequence was determined. The deduced amino acid sequence of REG2 was found to be 68% identical to that of the maize GLB2 globulin. Reg2 mRNA was present at high levels during embryo development for up to 14 days after flowering (DAF). Lower levels were found 20 DAF when the maturation of embryos was almost completed, and at the dry mature stage. Reg2 mRNA almost disappeared upon imbibition of isolated dry mature embryos but it was re-induced at a low level by further treatment with ABA. The expression of Reg2 was not induced by ABA in suspension-cultured cells, unlike that of Osem, one of the late embryogenesis abundant protein (LEA) genes.     

The hepato-renal syndrome: renal amino acid transport in bile duct ligated rats (DL)--influence of treatment with triiodothyronine or dexamethasone on renal amino acid handling in amino acid loaded rats

The influence of triiodothyronine or dexamethasone on renal amino acid handling was investigated in anaesthetized, bile duct-ligated (DL) adult female rats. 3 days after DL, glomerular filtration rate (GFR) was unchanged whereas urine flow was decreased. Plasma concentrations of 5 out of 16 amino acids were significantly enhanced after DL. On the other hand, the fractional excretion (FE) of 11 out of 16 amino acids was significantly reduced as a sign of improved reabsorption capacity. Bolus injections of leucine (20 mg/100 g b.wt.), glutamine (45 mg/100 g b.wt.), or taurine (45 mg/100 g b.wt.) were followed by a temporary increase in the FE of the administered amino acids as well of the endogenous amino acids which were not administered. This phenomenon was more pronounced in DL than in control rats. Under load conditions, dexamethasone (60 microg/100 g b.wt.) or triiodothyronine (20 microg/100 g b.wt.) treatment for 3 days, i.p. once daily, was followed by a stimulation of renal amino acid reabsorption in DL rats. The increase in fractional amino acid excretion after amino acid load was significantly lower than in untreated rats. This effect was also more pronounced in DL rats.     

Barley beta-glucosidase: expression during seed germination and maturation and partial amino acid sequences

Confusion regarding proper use of the terms rate and risk persists in the literature. This has implications for the proper modeling of prognosis and transition between health states in decision analysis and related techniques. The issue is complicated by the plethora of terms related to rate and risk. Although the suggestion to use the terms force and probability as substitutes for rate and risk has some appeal, the change in terminology by itself is unlikely to solve all the confusion or misuse of terms. This paper clarifies the proper definitions and estimations of rates and risks and suggests critical factors for the decision analyst to remember when using, modeling, or interpreting transition rates and risks.