1 Bioinformatics Laboratory, International Institute of Cell and Molecular Biology, ul. ks. Trojdena 4, 02109 Warsaw and 2 BioInfoBank, ul. Limanowskiego 24A, 60744 Poznan, Poland
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Keywords: bioinformatics/homology modeling/molecular evolution/mRNA capping/protein structure/RNA methyltransferase
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
While the `minimal' capping apparatus from yeast and vaccinia has been extensively characterized, little is known about the cap 1 and cap 2 MTases except for the vaccinia virus protein VP39, which acts as a cap 1 MTase and a non-catalytic, smaller subunit of the heterodimeric poly(A) polymerase. Structural and functional studies on VP39 have highlighted many aspects of this protein's 2'-O-MTase function (Hodel et al., 1996, 1998
, 1999
). All the cap MTases characterized so far utilize S-adenosylmethionine (AdoMet) as the methyl group donor and belong to a family of proteins with a common, albeit weakly conserved, signature of the cofactor-binding region (Kagan and Clarke, 1994
; Cheng and Blumenthal, 1999
). Therefore, one may envisage identification of other classes of cap MTases and phylogenetic classification based on their degree of similarity to the VP39 protein. However, sequence similarity searches failed to identify any homologs of VP39 outside pox viruses, presumably due to the extreme divergence of subfamilies of various viral and cellular cap 1 MTases. On the other hand, extensive sequence analysis of proteins from positive-strand RNA viruses carried out by Koonin and co-workers allowed them to delineate a putative MTase domain, tentatively associated with the activity required for cap formation (Koonin et al., 1992
; Rozanov et al., 1992
; Koonin, 1993
). Nevertheless, it was unclear which of the cap methylations are carried out by the predicted domain, since its sequence did not show pronounced similarity to any other cap MTase subfamilies.
Recently, the structure of another 2'-O-MTase, namely the 23S rRNA:U2552 specific RrmJ MTase from Escherichia coli, has been solved at 1.5 Å resolution, revealing the three-dimensional fold and architecture of the active site common with VP39, despite the fact that these proteins lack overall sequence similarity (Bugl et al., 2000). Also, in the mammalian reovirus
2 protein structure solved at 3.6 Å resolution, two AdoMet-dependent MTase-like domains have been delineated (Reinisch et al., 2000
). Remarkably, owing to the lack of similarity to other proteins in the sequence databases, the authors could not unambiguously determine which of these domains is the cap 0 and which is the cap 1 MTase and, according to our analysis, their tentative assignment should be reversed (Bujnicki and Rychlewski, 2001
). Problems with identification of the sequence elements specific for 2'-O-ribose MTases are a good illustration of the degree of divergence in this protein family.
The non-segmented, negative-sense, single-stranded RNA viruses [order Mononegavirales (MNV)] comprise many human and animal pathogens of significant epidemiological importance, including respiratory syncytial virus, measles, mumps, rabies, parainfluenza, vesicular stomatitis and Marburg and Ebola viruses, and several plant pathogens (abbreviations of the names of viruses analyzed in this work are given in Table I). Complete nucleotide sequences have been determined for more than 20 MNV, revealing five common genes maintained in highly similar order (NPMGL), varied in some cases by insertions (Pringle and Easton, 1997
; Conzelmann, 1998
). One of the bona fide homologous components of the viral ribonucleoprotein core is the large (L) protein, which functions as the RNA-dependent RNA polymerase (Tordo et al., 1992
). Other activities attributed to the L protein include mRNA capping, cap 0 and cap 1 methylation, poly(A) polymerase and protein kinase.
|
In vesicular stomatitis virus (VSV), the monomethylated G(5')ppp(5')Am structure, generated by the cap 1 MTase, is a preferred substrate for the cap 0 (m7G) MTase (Testa and Banerjee, 1977; Hammond and Lesnaw, 1987
). In this respect, the VSV system is distinct from the vaccinia and reovirus systems, in which the cap 0 structure is necessary for cap 1 methylation to occur (Furuichi and Shatkin, 2000
). It has been shown that the VSV mRNAs lacking the cap 0 structures are poor templates for protein synthesis in vitro (Testa and Banerjee, 1977
; Horikami and Moyer, 1982
). Therefore, inhibitors of either cap 0 or cap 1 MTase activities of the L protein may be used as antiviral drugs, but to date the regions involved in catalysis of the methyl transfer have not been mapped on to the primary structure of the L protein, which hampers knowledge-based drug design. Hence, in silico structure prediction of the domains responsible for the L protein activities is of obvious importance.
In this paper, we report the identification of a putative MTase domain in the L protein. Based on extensive bioinformatics analysis, including iterative database searches, structure prediction and molecular modeling, we demonstrate that this domain shares key features with known 2'-O-ribose MTases. We analyze sequence and structural similarities to other 2'-O-ribose MTases, propose possible roles for conserved residues and discuss evolutionary relationships among the representative members of the family of the putative cap 1 MTases of MNV.
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The non-redundant (nr) database at NCBI was extensively searched with the PSI-BLAST algorithm (Altschul et al., 1997), using the sequence of RrmJ, VP39 and subsequently some newly retrieved representative sequences as queries. The expectation (e)-value cutoff value was varied in the range 1061012, depending on the visual inspection of the alignments reported by the program for different queries. Low-complexity sequence regions were left unmasked. Full-length protein sequence alignments were reconstructed using ClustalX (Thompson et al., 1997
) based of the degapped PSI-BLAST output processed using BIB-VIEW (http://bioinfo.pl/bibview.pl). Secondary and tertiary structure predictions were carried out via the protein structure prediction MetaServer-Pcons interface (http://bioinfo.pl/meta/) (Bujnicki et al., 2001
; Lundstrom et al., 2001
) using 150400 amino acid fragments of the L protein as queries.
Molecular modeling
Homology modeling was carried out following a modified version of the `multiple models' approach (Pawlowski et al., 1997) using MODELLER (Sali and Blundell, 1993
) to generate several alternative preliminary models based on threading-derived pairwise targettemplate alignments and PROMODII (Guex and Peitsch, 1997
) to merge the best-scored fragments of preliminary models into the final structure. The preliminary models were obtained using unrefined pairwise alignments reported by PSI-BLAST (Altschul et al., 1997
), FFAS (Rychlewski et al., 2000
), 3DPSSM (Kelley et al., 2000
), BIOINBGU (Fischer, 2000
) and GenThreader (Jones, 1999
) and merged using secondary structure prediction results and energy evaluation to resolve ambiguities. The structure of an insertion 20 aa long was predicted using the ab initio protein folding server I-SITES/ROSETTA (Simons et al., 1997
) and manually inserted into the homology-modeled core. Energy minimization was carried out using GROMOS96 (Scott et al., 1999
) until all inconsistencies in geometry were rectified and all the short contacts were relieved. The stereochemical and energetic properties of modeling intermediates and of the final model were evaluated using WHATCHECK (Hooft et al., 1996
) and VERIFY3D (Eisenberg et al., 1997
). Semi-automated and manual manipulations with protein structures and sequencestructure alignments were conducted using SWISS-PDB VIEWER (Guex and Peitsch, 1997
).
Phylogenetic analysis
Phylogenetic inference was carried out using the conserved regions of the refined sequencestructure alignment of catalytic domains of the predicted viral cap 1 MTase based on the neighbor-joining method (Saitou and Nei, 1987). A corrected distance matrix was calculated from sequences according to the JTT model (Jones et al., 1992
). Bootstrapping analysis was performed, generating 100 replicates of the sequence alignment. The majority-rule consensus tree was visualized using TREEVIEW (Page, 1996
).
![]() |
Results and discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The evolutionary relationships among various nucleic acid MTase families have been studied (Gustafsson et al., 1996), but for a long time there was no indication of any structural or evolutionary relatedness between various families of 2'-O-ribose MTases. As a part of a larger project, aiming at identification and classification of novel RNA MTases among the uncharacterized or putative proteins in sequence databases, we conducted exhaustive database searches using sequences and structures of bona fide 2'-O-ribose MTases as queries. VP39 and its close homologs turned out to be a poor query, since even with the relaxed cutoff (e-value <0.1) all PSI-BLAST (Altschul et al., 1997
) searches converged after several iterations yielding no sequences except the 2'-O-ribose MTases from poxviruses. Conversely, when the sequence of E.coli RrmJ was used as an initial probe with a default cutoff (e-value = 103), the search `exploded', ultimately reporting similarities to many AdoMet-dependent MTase families, including enzymes known to modify proteins, various small molecules and bases in RNA and DNA (data not shown). Therefore, we tested alternative cutoff values, monitoring the addition of MTases exhibiting conservation of both the AdoMet-binding region and the ribose-binding residues common to VP39 and RrmJ (Bugl et al., 2000
).
The stringent cutoff of 5x106 turned out to be optimal for the search for putative 2'-O-ribose MTases initiated with the RrmJ sequence, since all hits reported retained the residues implicated in ribose binding in structurally characterized members of the family. The second iteration revealed significant similarity (e-value = 9x107) of RrmJ and its homologs to the putative methyltransferase domain of the large non-structural protein of ssRNA positive-strand viruses (Rozanov et al., 1992; Koonin, 1993
; Koonin and Dolja, 1993
). Remarkably, whereas Koonin and co-workers could not predict the function of cap 0 or cap 1 MTase (or both) based on their analysis of amino acid conservation patterns, our study of recently solved structures strongly suggests that the MTase domain that they identified shares the key features of the 2'-O-ribose MTase and therefore most likely functions as the cap 1 MTase. Iterating the PSI-BLAST search resulted in the accumulation of numerous closely related sequences of the putative cap 1 MTase of positive-strand viruses and members of the TlyA family (Aravind and Koonin, 1999
). Further, in the third iteration a similarity of 2'-O-MTases to a fragment of the RNA polymerase (L protein) sequence from the mumps virus, which belongs to the genus Rubulavirus of ssRNA negative-strand viruses, has been reported with a score of 4x106 and in further iterations more related sequences from rubulaviruses and paramyxoviruses were retrieved with e-values as high as 1013. We have also detected a new putative 2'-O-MTase YgdE in E.coli and related bacteria (Bujnicki and Rychlewski, 2000
).
Visual analysis of the PSI-BLAST output saved as multiple sequence alignment revealed that all sequence fragments of the L protein of negative-strand viruses retained the presumptive ribose-binding C-terminal subdomain with three conserved side chains, but lacked the N-terminal AdoMet-binding subdomain with the characteristic glycine-rich motif I. Consequently, the fourth residue typical of 2'-O-MTases, namely invariant Lys from motif X, which in RrmJ is most proximal to the N-terminus, could not be observed in the initial alignment. To test the possibility that the L protein of viruses from the family Paramyxoviridae and possibly from other MNV includes an intact domain similar to 2'-O-MTases, we followed two distinct strategies. Firstly, we carried out a series of reciprocal PSI-BLAST searches with default cutoff values, using the reported fragments of the L protein with varying extensions at both termini to localize the previously missed conserved Lys residue and motif I. Secondly, we used the same sequence fragments to carry out sequence-to-structure threading in order to determine if the domain under consideration is indeed similar to the structurally characterized MTases. Both the sequence- and structure-based strategy turned out to be successful. Querying PSI-BLAST with the fragment of the mumps virus L protein spanning residues 17642010 resulted in a multiple sequence alignment encompassing both of the previously identified C-termini of the catalytic domain, and also two additional sequence blocks, one of which contained an invariant Lys residue and the other closely resembled the GxGxG pattern typical for the conserved AdoMet-binding loop (Kagan and Clarke, 1994; Fauman et al., 1999
). On the other hand, threading algorithms evaluated the MTase fold as the most likely candidate for the structure of the query sequence. For instance, 3D-PSSM reported similarity of the above-mentioned region to the RrmJ structure with a highly significant score of 0.00271 (the detailed results of the search reported by all servers are available online at http://bioinfo.pl/meta/target.pl?id=2357). The results of the threading analysis carried out for other sequences in the data set of potential MTase domains extracted from the L proteins confirmed their high propensity to assume the MTase fold. Further, carrying out PSI-BLAST searches for the same sequence fragments resulted in the retrieval of a bulk of sequences from MNV in the earliest iterations, followed by eukaryotic and prokaryotic homologs of the RrmJ MTase and the putative cap 1 MTases from positive-strand viruses (data not shown).
Interestingly, in the sixth iteration of the search initiated with the PI2HT putative MTase domain a hit to the grass carp reovirus (GCRV) VP1 protein has been reported. This region in GCRV VP1 corresponds to a part of the MTase I domain in the recently solved crystal structure of the human reovirus 2 protein. Extending the profile-to-profile alignment of sequences from the MNV to the pair of reoviral proteins using FFAS (Rychlewski et al., 2000
) resulted in a perfect match of the four invariant residues (Figure 1
). No significant overall sequence similarity or generally conserved residues was observed when we attempted to align the 2'-O-MTase family to the MTase II domain. This result supports our independent prediction, based on analysis of protein structures alone, that the reoviral MTase I domain is more similar to 2'-O-ribose cap 1 MTases, while the MTase II domain most likely functions as the m7G-specific cap 2 MTase (Bujnicki and Rychlewski, 2001
|
In addition to the conserved motifs common for the majority of MTases, the putative cap 1 MTase domain in the L protein contains a variable region between motifs II and III. This region is absent from other 2'-O-ribose MTases. Interestingly, proteins from the subfamily of DNA:m6A MTases possess a variable region between motifs II and III, which forms an autonomous domain (the so-called `target recognition domain', TRD) implicated in recognition of the specific sequence in the DNA (Tran et al., 1998
). The chemotaxis receptor methyltransferase CheR also possesses an additional small domain in the same region; this domain is dissimilar to the TRD of
-m6A MTases and is involved in specific interactions with the methylated receptor (Djordjevic and Stock, 1998
). However, the `variable' domains of
-m6A and CheR MTases are of relatively constant length and exhibit conservation of the key hydrophobic residues, suggesting that a common fold is retained within each family, while the length of the insertion present in the predicted cap 1 MTase domain is extremely variable. Moreover, we could detect only one Pro residue that may be conserved in this subfamily, which suggests that this part of the protein is structurally variable.
Molecular modeling
It is known that detection of distant homologs either by sequence searches or by threading does not necessarily translate into correct alignments, from which sequencestructure function relationships could be accurately inferred (Smith et al., 1997). The independent verification of threading, refinement of sequence alignments with known structures and estimation of reliability of certain regions in the alignments can be addressed by modeling. For that reason and to gain insight into the molecular basis of intriguing similarities in the active site and into considerable differences in the sequence of the `variable' regions, the structure of the putative 2'-O-ribose MTase of Ebola virus was predicted by homology modeling, using the coordinates of RrmJ as the template. We did not use the structures of the vaccinia virus and reovirus cap 1 MTases, since their sequences exhibited much lower similarity to the whole range of potential target sequences from MNV and we observed large discrepancies between alternative threading-based alignments including these structures. Besides, we could not identify a common cap-binding motif in these two MTases or a region at the N- or C-terminus of the predicted MTase domain in the L protein that would exhibit similarity to the cap-binding site of either cap 1 MTase of known structure. Therefore, only the conserved core spanning the AdoMet- and ribose-binding sites could be predicted with confidence. Still, the modeling was not trivial, since RrmJ and the viral proteins shared low sequence identity and therefore we resorted to a modified version of the `multiple models approach' of Pawlowski et al. (Pawlowski et al., 1997
) (see Materials and methods) to ensure that the number of possibly misaligned sequence segments between the target and the template was reduced to the minimum. The preliminary models showed very good agreement in the predicted core elements; there were only few relative shifts of sequence segments, almost exclusively in the peripheral structures, which very strongly supports the presented fold-recognition results (data not shown).
The final averaged and optimized model passed all the tests implemented in the stereochemistry-evaluating WHATCHECK suite (Hooft et al., 1996) and in the VERIFY3D program, which uses contact potentials to assess whether the modeled amino acid residues occur in the environment typical for globular proteins with hydrophobic core and solvent-exposed surface (Eisenberg et al., 1997
). It is worth emphasizing that the stereochemistry of even a plain wrong model can be refined to an acceptable degree; however, the calculation of energy based on observed contacts would still indicate that the polypeptide chain is misfolded. Moreover, the reasonable energies are rarely observed for misfolded structures. Thus, the scores reported for our model by WHATCHECK (Z-score 4.1) and VERIFY3D (average score 0.3, no regions scored lower than 0) suggest that both its three-dimensional fold and the conformation of individual residues are reasonable.
Sequencestructurefunction relationships
The modeled putative cap 1 MTase domain of the Ebola virus L protein, (aa 18082009) resembles its modeling template, structure of the 2'-O-ribose RrmJ MTase (Figure 2). All insertions in the sequence of the viral MTase are localized in the connectors between the secondary structure elements of RrmJ and are modeled as loops; these modifications do not interfere with the active site of the MTase. The only major difference between the two structures is the insertion in the viral protein that forms an elaboration of the edge of the common ß-sheet, distant from the cofactor- and ribose-binding sites. It is rather unlikely that it is involved in recognition of the substrate; however, in the absence of experimental data it is not possible to propose any function for this region other than interactions with other proteins, possibly other domains of the L protein.
|
Phylogenetic analysis
In spite of the low degree of sequence similarity among putative cap 1 MTases from individual families of MNV, the present study shows that they all originate from a common ancestral enzyme. To evaluate the evolutionary relationships between the MTase domains and verify if they are consistent with relationships between the polymerase domain and the established taxonomy of MNV (http://www.ncbi.nlm.nih.gov/ICTV/overview/negssrna.html) (van Regenmortel et al., 2000), phylogenetic trees were inferred from the alignment as described in Materials and methods. A consensus tree inferred using the neighbor-joining method with default parameters is shown in Figure 3.
|
Conclusions
Using iterative searches of sequence databases and sequence-to-structure threading, we have detected the presence of a putative 2'-O-ribose (cap 1) MTase domain in the L protein of MNV and predicted the key residues involved in binding and the methyl transfer reaction. The homology of the putative cap 1 MTase domain of MNV to other 2'-O-ribose MTases, which served as a basis for comparative molecular modeling and phylogenetic analysis, is not prima facie evident and is also not recognizable by standard, non-iterated algorithms for pairwise sequence comparison. Therefore, the alignment presented in this work will be a good starting point for the creation of specific sequence profiles used for identification of other 2'-O-ribose MTases in genomic data and the structural model may be useful in development of antiviral drugs. Our results will help in understanding better the biology of MNV and their relationship to other viruses. Further biochemical and structural studies of the L proteins from MNV, including site-directed mutagenesis and X-ray crystallography, may be used in clarifying the structural constraints imposed by the function of these enzymes and their divergence from the common ancestor shared with many ribose MTases and other AdoMet-dependent enzymes. In addition, our phylogenetic analysis revealed the unusual relationship of the cap 1 MTase domain of pneumoviruses to other paramyxoviruses, which suggests that the evolution of the L proteins should be studied based on comparison of a series of phylogenies inferred from alignments of all individual domains.
![]() |
Notes |
---|
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Aravind,L. and Koonin,E.V. (1999) J. Mol. Evol., 48, 291302.[ISI][Medline]
Banerjee,A.K. (1980) Microbiol. Rev., 44, 175205.[ISI]
Bugl,H., Fauman,E.B., Staker,B.L., Zheng,F., Kushner,S.R., Saper,M.A., Bardwell,J.C. and Jakob,U. (2000) Mol. Cell, 6, 349360.[ISI][Medline]
Bujnicki,J.M. and Rychlewski,L. (2000) Acta Microbiol. Pol., 49, 253260.[Medline]
Bujnicki,J.M. and Rychlewski,L. (2001) Genome Biol., 2,38.138.6
Bujnicki,J.M., Elofsson,A., Fischer,D. and Rychlewski,L. (2001) Bioinformatics, 17, 750751.
Cheng,X. and Blumenthal,R.M. 1999. S-Adenosylmethionine-dependent Methyltransferases: Structures and Functions. World Scientific, Singapore.
Conzelmann,K.K. (1998) Annu. Rev. Genet., 32, 123162.[CrossRef][ISI][Medline]
Djordjevic,S. and Stock,A.M. (1998) Nature Struct. Biol., 5, 446450.[ISI][Medline]
Eisenberg,D., Luthy,R. and Bowie,J.U. (1997) Methods Enzymol., 277, 396404.[ISI][Medline]
Fauman,E.B., Blumenthal,R.M. and Cheng,X. (1999) In Cheng,X. and Blumenthal,R.M. (eds), S-Adenosylmethionine-dependent Methyltransferases: Structures and Functions. World Scientific, Singapore, pp. 138.
Fischer,D. (2000) Pac. Symp. Biocomput., 119130.
Fischer,U. and Luhrmann,R. (1990) Science, 249, 786790.[ISI][Medline]
Flamand,A. and Bishop,D.H. (1973) J. Virol., 12, 12381252.[ISI][Medline]
Furuichi,Y. and Shatkin,A.J. (2000) Adv. Virus Res., 55, 135184.[ISI][Medline]
Guex,N. and Peitsch,M.C. (1997) Electrophoresis, 18, 27142723.[ISI][Medline]
Gustafsson,C., Reid,R., Greene,P.J. and Santi,D.V. (1996) Nucleic Acids Res., 24, 37563762.
Hammond,D.C. and Lesnaw,J.A. (1987) Virology, 159, 229236.[ISI][Medline]
Hodel,A.E., Gershon,P.D., Shi,X. and Quiocho,F.A. (1996) Cell 85,247256.
Hodel,A.E., Gershon,P.D. and Quiocho,F.A. (1998) Mol. Cell, 1, 443447.[ISI][Medline]
Hodel,A.E., Quiocho,F.A. and Gershon,P.D. (1999) In Cheng,X. and Blumenthal,R.M. (eds), S-Adenosylmethionine-dependent Methyltransferases: Structures and Functions. World Scientific, Singapore, pp. 255282.
Hooft,R.W., Vriend,G., Sander,C. and Abola,E.E. (1996) Nature, 381, 272[ISI][Medline]
Horikami,S.R. and Moyer,S.A. (1982) Proc. Natl Acad. Sci. USA, 79, 76947698.[Abstract]
Jones,D.T. (1999) J. Mol. Biol., 287, 797815.[CrossRef][ISI][Medline]
Jones,D.T., Taylor,W.R. and Thornton,J.M. (1992) Comput. Appl. Biosci., 8, 275282.[Abstract]
Kagan,R.M. and Clarke,S. (1994) Arch. Biochem. Biophys., 310, 417427.[CrossRef][ISI][Medline]
Kelley,L.A., McCallum,C.M. and Sternberg,M.J. (2000) J. Mol. Biol., 299, 501522.[CrossRef]
Koonin,E.V. (1993) J. Gen. Virol., 74, 733740.[Abstract]
Koonin,E.V. and Dolja,V.V. (1993) Crit. Rev. Biochem. Mol. Biol., 28, 375430.[Abstract]
Koonin,E.V., Gorbalenya,A.E., Purdy,M.A., Rozanov,M.N., Reyes,G.R. and Bradley,D.W. (1992) Proc. Natl Acad. Sci. USA, 89, 82598263.[Abstract]
Langberg,S.R. and Moss,B. (1981) J Biol. Chem., 256, 1005410060.
Lundstrom,J., Rychlewski,L., Bujnicki,J.M. and Elofsson,A. (2001) Protein Sci., 10, 23542362.
McClure,M.A. and Perrault,J. (1989) Virology, 172, 391397.[ISI][Medline]
Page,R.D. (1996) Comput. Appl. Biosci., 12, 357358.[Medline]
Pawlowski,K., Jaroszewski,L., Bierzynski,A. and Godzik,A. (1997) Pac. Symp. Biocomput.,328339.
Pringle,C.R. and Easton,A.J. (1997) Semin. Virol., 8, 4957.[CrossRef]
Reinisch,K.M., Nibert,M.L. and Harrison,S.C. (2000) Nature, 404, 960967.[CrossRef][ISI][Medline]
Rozanov,M.N., Koonin,E.V. and Gorbalenya,A.E. (1992) J. Gen. Virol., 73, 21292134.[Abstract]
Rychlewski,L., Jaroszewski,L., Li,W. and Godzik,A. (2000) Protein Sci., 9, 232241.[Abstract]
Saitou,N. and Nei,M. (1987) Mol. Biol. Evol., 4, 406425.[Abstract]
Sali,A. and Blundell,T.L. (1993) J. Mol. Biol., 234, 779815.[CrossRef][ISI][Medline]
Scott,W.R.P. et al. (1999) J. Phys. Chem., 103, 35963607.[ISI]
Simons,K.T., Kooperberg,C., Huang,E. and Baker,D. (1997) J. Mol. Biol., 268, 209225.[CrossRef][ISI][Medline]
Smith,T.F., Lo,C.L., Bienkowska,J., Gaitatzes,C., Rogers,R.G.J. and Lathrop,R. (1997) J. Comput. Biol., 4, 217225.[ISI][Medline]
Stec,D.S., Hill,M.G. and Collins,P.L. (1991) Virology, 183, 273287.[ISI][Medline]
Testa,D. and Banerjee,A.K. (1977) J. Virol., 24, 786793.[ISI][Medline]
Thompson,J.D., Gibson,T.J., Plewniak,F., Jeanmougin,F. and Higgins,D.G. (1997) Nucleic Acids Res., 25, 48764882.
Tordo,N., de Haan,P., Goldbach,R. and Poch,O. (1992) Virology, 3, 341357.[CrossRef]
Tran,P.H., Korszun,Z.R., Cerritelli,S., Springhorn,S.S. and Lacks,S.A. (1998) Structure, 6, 15631575.[ISI][Medline]
van Regenmortel,M.H.V. et al. (2000) Virus Taxonomy: the Classification and Nomenclature of Viruses. The Seventh Report of the International Committee on Taxonomy of Viruses. Academic Press, San Diego.
Varani,G. (1997) Structure, 5, 855858.[ISI][Medline]
Zanotto,P.M., Gibbs,M.J., Gould,E.A. and Holmes,E.C. (1996) J. Virol., 70, 60836096.[Abstract]
Received March 27, 2001; revised October 5, 2001; accepted November 11, 2001.