Related Phages Research Articles

Residues at the carboxy terminus of T4 DNA polymerase are important determinants for interaction with the polymerase accessory proteins.

Three T4 DNA polymerase accessory proteins (44P/62P and 45P) stimulate the polymerase (pol) activity and the 3'-5' exonuclease (exo) activity of T4 DNA polymerase (43P) on long, double-stranded DNA substrates. The 44P/62P "clamp loader" facilitates the binding of 45P, the "sliding clamp", to DNA that is primed for replication. Using a series of truncated 43P mutants, we identified a region at the extreme carboxy terminus of the DNA polymerase that is required for its interaction with accessory proteins. Truncation mutants of 43P lacking the carboxy-terminal 3, 6, or 11 residues retained full pol and exo activity on short synthetic primer-templates. However, the ability of the accessory proteins to enhance these activities on long double-stranded DNA templates was drastically reduced, and the extent of the reduction in activity was greater as more residues were deleted. One of the truncation mutants (N881), which had 17 residues removed from the carboxy terminus, showed reduced binding affinity and diminished pol activity but enhanced exo activity upon incubation with a small primer-template. The exo activity of the N881 mutant, on short, single-stranded DNA was unchanged, however, compared to the wild-type enzyme. These results are consistent with inferences drawn from the crystal structure of a DNA polymerase from a related T-even phage, RB69, where the carboxy-terminal 12 residues (equivalent to the 11 residues of 43P from phage T4) protrude from the thumb domain and are free to interact with complementary surfaces of the accessory proteins. The structural integrity of the thumb region in the N881 mutant is probably perturbed and could account for its reduced binding affinity and pol activity when incubated with short, double-stranded DNA substrates.

Temperate Myxococcus xanthus phage Mx8 encodes a DNA adenine methylase, Mox.

Temperate bacteriophage Mx8 of Myxococcus xanthus encapsidates terminally repetitious DNA, packaged as circular permutations of its 49-kbp genome. During both lytic and lysogenic development, Mx8 expresses a nonessential DNA methylase, Mox, which modifies adenine residues in occurrences of XhoI and PstI recognition sites, CTCGAG and CTGCAG, respectively, on both phage DNA and the host chromosome. The mox gene is necessary for methylase activity in vivo, because an amber mutation in the mox gene abolishes activity. The mox gene is the only phage gene required for methylase activity in vivo, because ectopic expression of mox as part of the M. xanthus mglBA operon results in partial methylation of the host chromosome. The predicted amino acid sequence of Mox is related most closely to that of the methylase involved in the cell cycle control of Caulobacter crescentus. We speculate that Mox acts to protect Mx8 phage DNA against restriction upon infection of a subset of natural M. xanthus hosts. One natural isolate of M. xanthus, the lysogenic source of related phage Mx81, produces a restriction endonuclease with the cleavage specificity of endonuclease BstBI.

Kinetic analysis of the endonuclease activity of phage lambda terminase: assembly of a catalytically competent nicking complex is rate-limiting.

The terminase enzyme from bacteriophage lambda is responsible for excision of a single genome from a concatameric DNA precursor and its insertion into an empty viral procapsid. The enzyme possesses a site-specific endonuclease activity which is responsible for excision of the viral genome and the formation of the 12 base-pair single-stranded "sticky" ends of mature lambda DNA. We have previously reported a kinetic analysis of the endonuclease activity of lambda terminase which showed an enzyme concentration-dependent change in the kinetic time course of the reaction [Tomka, M. A., & Catalano, C. E. (1993b) J. Biol. Chem. 268, 3056-3065]. We presented a model which suggested that the rate-limiting step in the nuclease reaction was the assembly of a catalytically competent prenicking complex. Here, we provide additional evidence for a slow assembly step in the nuclease reaction and demonstrate that the observed rate is affected by protein concentration, but not by the length of the DNA substrate. Consistent with our model, preincubation of terminase with DNA also yields an observable fast phase of the reaction, but only when large (> or = 3 kb) DNA substrates are used. Finally, we present data which demonstrate that phage lambda terminase can efficiently utilize DNA from the closely related phage phi21 as an endonuclease substrate and that the enzyme binds efficiently to the cosB region of both phage genomes. The implications of these results with respect to the assembly of a catalytically competent nucleoprotein complex required to initiate genome packaging are discussed.

Molecular Divergence of Lysozymes and α-Lactalbumin

The vast number of proteins that sustain the currently living organisms have been generated from a relatively small number of ancestral genes that has involved a variety of processes. Lysozyme is an ancient protein whose origin goes back an estimated 400 to 600 million years. This protein was originally a bacteriolytic defensive agent and has been adapted to serve a digestive function on at least two occasions, separated by nearly 40 million years. The origins of the related goose type and T4 phage lysozyme that are distinct from the more common C type are obscure. They share no discernable amino acid sequence identity and yet they possess common secondary and tertiary structures. Lysozyme C gene also gave rise, after gene duplication 300 to 400 million years ago, to a gene that currently codes for α-lactalbumin, a protein expressed only in the lactating mammary gland of all but a few species of mammals. It is required for the synthesis of lactose, the sugar secreted in milk. α-Lactalbumin shares only 40% identity in amino acid sequence with lysozyme C, but it has a closer spatial structure and gene organization. Although structurally similar, functionally they are quite distinct. Specific amino acid substitutions in α-lactalbumin account for the loss of the enzyme activity of lysozyme and the acquisition of the features necessary for its role in lactose synthesis. Evolutionary implications are as yet unclear but are being unraveled in many laboratories.

Gene organization and transcription of a late-expressed region of a Lactococcus lactis phage.

The lactococcal phage bIL41 belongs to the small isometric-headed phages of the 936 quasi-species and is resistant to the abortive infection determined by abiB. A 10.2-kb segment from this phage, in which late transcription is initiated, has been sequenced. Thirteen open reading frames (ORFs) organized in one transcriptional unit have been identified. The location of two of them and the structural features of the proteins they code for are evocative of terminase subunits. Five other ORFs specify proteins which are highly homologous to structural proteins from the closely related phage F4-1. By comparing the phage bIL41 sequence with partial sequences available for four related phages, we were able to deduce a chimerical phage map covering the middle- and a large part of the late-expressed regions. Phages from this quasi-species differ by the insertion or deletion of either 1 to about 400 bp in noncoding regions or an entire ORF. Transcription was initiated 9 min after infection at a promoter with a -10 but no -35 consensus sequence. Synthesis of a phage activator protein was needed for initiation of transcription. A large 16-kb transcript covering all of the late-expressed region of the genome was synthesized. This transcript gave rise to smaller units. One of these units most probably resulted from a RNase E processing.

Activation of Replication Origins in ϕ29-related Phages Requires the Recognition of Initiation Proteins to Specific Nucleoprotein Complexes

Protein p6 of Bacillus subtilis phage phi29 activates the initiation of viral DNA replication by forming a multimeric nucleoprotein complex at the origins of replication, located at both ends of the linear genome. This activation requires a precise positioning of the protein p6 array with respect to the initiation site. To investigate this activation mechanism, we have purified the phi29 protein p6 counterparts from the related phages Nf and GA-1 and analyzed the formation of complexes with DNA. In the homologous protein p6-DNA complexes the phi29 and Nf protein arrays showed an identical positioning, different than that of the GA-1 protein array. In contrast, in the heterologous complexes the protein showed a different arrangement except in the case of the Nf protein-phi29 DNA complex. We have also purified the proteins involved in the initiation of replication (terminal protein and DNA polymerase) from phages Nf and GA-1 and measured the ability of the different p6 proteins to activate homologous and heterologous replication origins. The results obtained indicate that the activation requires not only the formation of a specific nucleoprotein complex but also its specific recognition by the proteins involved in the initiation of DNA replication.

Short synthetic CDR-peptides forming the antibody combining site of the monoclonal antibody against RNA bacteriophage fr neutralize the phage activity

The construction of a mouse hybridoma FR52 secreting neutralizing monoclonal antibody specific for RNA bacteriophages fr, MS2 and GA is reported. The genes encoding the variable domains of the monoclonal antibody FR52 heavy and light chains were cloned and sequenced and the corresponding complementarity determining region (CDR) peptides were chemically synthesized. The CDR-peptides were tested for their ability to neutralize the activity of RNA phage fr and related RNA phages MS2 and GA. The CDR-derived peptides H2, L2 and L3 interacted with the fr phage particles and neutralized fr phage activity. Two of these peptides--H2 and L3 also had the ability to neutralize partly the activity of related bacteriophage MS2, but L1 and especially L3 neutralize the activity of the RNA phage GA. These results provide an excellent system for further antibody-antigen interaction studies and raise the possibility that simple CDR-peptides may serve as a new class of anti-viral molecules.

Defective site-specific integration elements are present in the genome of virulent bacteriophage LL-H of Lactobacillus delbrueckii

The phage attachment site, attP, and the integrase-encoding gene, int, are sufficient to promote site-specific integration of the temperate phage mv4 genome into the chromosome of the Lactobacillus delbrueckii host (L. Dupont, B. Boizet-Bonhoure, M. Coddeville, F. Auvray, and P. Ritzenthaler, J. Bacteriol. 177:586--595, 1995). The mv4 genome region containing these elements was compared at the nucleotide and amino acid levels with that of the closely related virulent phage LL-H. Complex DNA rearrangements were identified; a truncated integrase gene and two sites homologous to the mv4 attP site were detected in the genome of the virulent phage LL-H. These observations suggest that the two phages derive from a common temperate ancestor.

The crystal structure of bacteriophage Qβ at 3.5 å resolution

Background: The capsid protein subunits of small RNA bacteriophages form a T=3 particle upon assembly and RNA encapsidation. Dimers of the capsid protein repress translation of the replicase gene product by binding to the ribosome binding site and this interaction is believed to initiate RNA encapsidation. We have determined the crystal structure of phage Qβ with the aim of clarifying which factors are the most important for particle assembly and RNA interaction in the small phages. Results The crystal structure of bacteriophage Qβ determined at 3.5 å resolution shows that the capsid is stabilized by disulfide bonds on each side of the flexible loops that are situated around the fivefold and quasi-sixfold axes. As in other small RNA phages, the protein capsid is constructed from subunits which associate into dimers. A contiguous ten-stranded antiparallel β sheet facing the RNA is formed in the dimer. The disulfide bonds lock the constituent dimers of the capsid covalently in the T=3 lattice. Conclusion The unusual stability of the Qβ particle is due to the tight dimer interactions and the disulfide bonds linking each dimer covalently to the rest of the capsid. A comparison with the structure of the related phage MS2 shows that although the fold of the Qβ coat protein is very similar, the details of the protein–protein interactions are completely different. The most conserved region of the protein is at the surface, which, in MS2, is involved in RNA binding.

A Genomic Region of Lactococcal Temperate Bacteriophage TP901-1 Encoding Major Virion Proteins

Two major structural proteins, MHP (major head protein) and MTP (major tail protein), from the lactococcal temperate phage TP901-1 were sequenced at their amino acid termini, and derived degenerate oligonucleotides were used to locate the corresponding genes in the phage genome. This genomic region was sequenced. The sequence characterized includes a total of 11 open reading frames (ORFs) showing an operon structure. Upstream of each ORF, except ORFb2and ORFx,potential ribosome-binding sites were found, suggesting independent translation. However, coupled translation is suggested for ORFxand as a possibility for ORFb3and ORFc2,which have ribosome-binding sites located more distant from their start codons. ORFb2may be translationally fused withmhpat a low frequency. Themhpandmtpgenes are transcribed as a 3.7-kb mRNA with at least six additional ORFs. The organization of the genomic region analyzed resembles that of other distantly related phages, providing possible roles for the uncharacterized ORFs.

Genetic evidence for an activator required for induction of colicin-like bacteriocin 28b production in Serratia marcescens by DNA-damaging agents.

Bacteriocin 28b production is induced by mitomycin in wild-type Serratia marcescens 2170 but not in Escherichia coli harboring the bacteriocin 28b structural gene (bss). Studies with a bss-lacZ transcriptional fusion showed that mitomycin increased the level of bss gene transcription in S. marcescens but not in the E. coli background. A S. marcescens Tn5 insertion mutant was obtained (S. marcescens 2170 reg::Tn5) whose bacteriocin 28b production and bss gene transcription were not increased by mitomycin treatment. Cloning and DNA sequencing of the mutated region showed that the Tn5 insertion was flanked by an SOS box sequence and three genes that are probably cotranscribed (regA, regB, and regC). These three genes had homology to phage holins, phage lysozymes, and the Ogr transcriptional activator of P2 and related bacteriophages, respectively. Recombinant plasmid containing this wild-type DNA region complemented the reg::Tn5 regulatory mutant. A transcriptional fusion between a 157-bp DNA fragment, containing the apparent SOS box upstream of the regA gene, and the cat gene showed increased chloramphenicol acetyltransferase activity upon mitomycin treatment. Upstream of the bss gene, a sequence similar to the consensus sequence proposed to bind Ogr protein was found, but no sequence similar to an SOS box was detected. Our results suggest that transcriptional induction of bacteriocin 28b upon mitomycin treatment is mediated by the regC gene whose own transcription would be LexA dependent.

Secondary Structure Model for the First Three Domains of Qβ RNA. Control of A-protein synthesis

We present a secondary structure model for the first 860 nucleotides of Qβ RNA. The model is supported by phylogenetic comparison, nuclease S 1structure probing and computer prediction using energy minimization and a Monte Carlo approach. To provide the necessary data for the comparative analysis we have sequenced the single-stranded RNA coliphages MX1, M11 and NL95. Together with the known sequences of Qβ and SP, this yields five sequences with sufficient sequence diversity to be useful for the analysis. The part of the Qβ genome examined contains the 60 nucleotide 5′ untranslated region and the first 800 nucleotides of the maturation protein gene. The RNA adopts a highly ordered structure in which all hairpins are held in place by a network of long-distance interactions, which form three-way and four-way junctions. Only the 5′-terminal hairpin is unrestrained, while connected by a few single-stranded nucleotides to the body of the RNA. The start region of the A-protein gene, which is part of the network of long-distance interactions, is base-paired to three non-contiguous downstream sequences. As a result, translation is expected to be progressively quenched when the length of the nascent chains increases. This feature explains the previous observation that A-protein synthesis on Qβ RNA can start only on short nascent strands. Translational control of the A protein in the distantly related phage MS2 was recently shown to be controlled by the kinetics of RNA folding. This basic difference and its possible biological purpose can be explained by the different RNA folding pathways in Qβ and MS2. Interestingly, due to the presence of G·U pairs, structure prediction for the minus strand differs in some aspects from that for the plus strand. More specifically, there is a minus-strand specific, long-distance interaction bordering the minus-strand equivalent of the 5′-terminal hairpin. This interaction extends at the expense of the lower part of the terminal helix, thereby exposing the terminal C residues at which replication starts. This long-distance interaction, which was recently shown to be required for minus-strand replication, is strongly supported by our comparative data.

Crystal Structures of MS2 Capsids with Mutations in the Subunit FG Loop

The loop between the F and G β strands (FG loop) of the bacteriophage MS2 coat protein subunit forms inter-subunit contacts around the 5-fold and 3-fold (quasi 6-fold) axes of the T=3 protein shell. In capsids, the loop is found in two very different conformations, one in B subunits, which form the 5-fold contact, and one in A and C subunits, which form the quasi 6-fold contact. One proline residue, Pro78, is strictly conserved in the coat protein of all related bacteriophages, and in the case of MS2 this proline residue is preceded by a cispeptide bond in the B subunit. In order to probe the role of the FG loop in capsid assembly, we have determined the crystal structures of two MS2 capsids, formed by coat proteins with mutations at two positions in the FG loop, P78N or E76D. These mutants show conformational changes in the FG loops that explain the reduced temperature stability of the capsids. The P78N mutant has a normal transpeptide bond at position 78.

Components of multiprotein-RNA complex that controls transcription elongation in Escherichia coli phage lambda

This chapter discusses the components of multiprotein-RNA complex that controls transcription elongation in Escherichia Coli phage λ. Studies of bacteriophage A led to the discovery of several basic mechanisms of transcriptional regulation. One of these is transcriptional antitermination, the process in which genes whose transcription is otherwise blocked by premature termination are expressed through termination suppression. In λ and related phages, genome-specific antiterminators convert RNA polymerase (RNAP) into a termination-resistant form during early phases of transcription elongation. The chapter describes the current understanding of how one such antiterminator, the λ N gene product, works. The chapter reviews the methods of overproduction, isolation, and assay of the N protein as well as several accessory factors that modulate transcription elongation in Escherichia coli in the form of a multiprotein-RNA complex. The availability of N and Nus factors in large scale should now make it possible not only to attempt to determine the structures of the individual protein components, and the various RNA-protein and protein-protein complexes by biophysical methods, but also to examine these interactions by conventional biochemical methods.

Thermodynamics of folding of the RNA pseudoknot of the T4 gene 32 autoregulatory messenger RNA.

Nucleotides U(-67) to C(-40) at the extreme 5' end of the gene 32 mRNA in bacteriophage T4 have been shown to fold into an RNA pseudoknot proposed to be important for translational autoregulation. The thermal denaturation of three in vitro transcribed RNAs corresponding to the pseudoknot region has been investigated as a function of Mg2+ concentration to begin to elucidate the determinants of the structure and stability of this conformation. T4-35 is a 35-nucleotide RNA containing a 5' G followed by the natural T4 sequence starting with the mature 5' end of the mRNA, nucleotides A(-71) to C(-38). A 32-nucleotide RNA, termed T4-32, contains the native sequence form U(-67) to C(40) with 5'GC and 5'CA single-stranded regions appended to the 5' and 3' ends of the core sequence, respectively. T4-28 contains only the 28 core nucleotides, and the predicted closing U(-67)-A(-52) base pair in stem 1 has been replaced with a phylogenetically allowed G(-67)-C(-52) base pair. Ribonuclease mapping of T4-32 and imino proton NMR experiments of T4-35 show that both sequences adopt a pseudoknotted conformation. At pH 6.9 and 50 mM NaCl, T4-35 and T4-32 RNAs are characterized by a single major melting transition over a wide range of [Mg2+] (0-6 mM). The delta H degree of unfolding for T4-35 and T4-32 shows a large dependence on Mg2+ concentration; the maximum delta H degree occurs at about 2.0 mM Mg2+ with further addition of Mg2+ simply increasing the tm. Investigation of the [Mg2+] dependence of the tm suggests that a net of one Mg2+ ion is released upon denaturation of T4-35 and T4-32 RNAs. Over the entire [Mg2+] range, the delta G degree (37 degrees C) for the folding of T4-35 is consistently 1-1.5 kcal mol(-1) more negative than T4-32 due to a higher stabilization enthalpy for the natural sequence molecule. In contrast to this behavior, T4-28 gives consistently higher tm's but less negative enthalpies and is destabilized (at 37 degrees C) by about 0.5-1.5 kcal mol(-1) relative to T4-32 and by about 2-3 kcal mol(-1) relative to T4-35, depending upon cation concentration. (1)H NMR experiments suggest that, even in the presence of 4.0 mM Mg2+, T4-28 RNA does not adopt a stable pseudoknotted conformation. These data show that the stability of the pseudoknot in the gene 32 mRNA encoded by the 28-nucleotide core sequence is significantly influenced by the number and nature of the immediately adjacent "single-stranded" 5' and/or 3' nucleotides appended to the core structure. These findings are discussed within the context of the structural model for the evolutionarily related phage T2 and T6 gene 32 mRNA pseudoknots presented in the following paper [Du, Z., Giedroc, D. P., & Hoffman, D. W. (1996) Biochemistry 35, 4187-4198].

Enhancer-independent variants of phage Mu transposase: enhancer-specific stimulation of catalytic activity by a partner transposase.

Assembly of the functional tetrameric form of phage Mu transposase (A protein) requires specific interactions between the Mu A monomer and its cognate sequences at the ends of the Mu genome (attL and attR) as well as those internal to it (the enhancer element). We describe here deletion variants of Mu A that show enhancer-independence in the assembly of the strand cleavage complex. These deletions remove the amino-terminal region of Mu A required for its interactions with the enhancer elements. The basal enhancer-independent activity of the variant proteins can be stimulated by a partner variant harboring an intact enhancer-binding domain. By exploiting the identical att-binding, and nonidentical enhancer-binding specificities of Mu A and D108 A (transposase of the Mu related phage D108), we show that the stimulation of activity is enhancer-specific. Taken together, these results suggest that the domain of Mu A that includes the enhancer-interacting region may exert negative as well as positive modulatory effects on the strand cleavage reaction. We discuss the implications of these results in the framework of a recent model for the assembly of shared active sites within the Mu A tetramer.

Virion Positions and Relationships of Lactococcal Temperate Bacteriophage TP901-1 Proteins

The major proteins of phage TP901-1 virion were characterized by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and structural relations were determined using specific antibodies, obtained by affinity purification from a polyclonal serum. A 23-kDa protein was identified as the major tail protein, and a 31-kDa molecule as the major head protein, respectively. Labeling experiments with antibodies against two proteins, with molecular masses of 20 and 19 kDa, indicated that they were baseplate-related components. A 72-kDa protein was found to be part of a neck passage structure, which includes a collar. Evidence for the presence of attached whiskers was also obtained. T7 RNA polymerase-mediated expression of the two major proteins confirmed the gene location of the previously sequenced region of the phage genome. The relation to other lactococcal phages was determined by DNA hybridization and antibody probing, showing that despite low DNA similarity, TP901-1 NPS epitopes were detected in both related and unrelated small isometric-headed phages.

A Domain of the Gene 4 Helicase/Primase of Bacteriophage T7 Required for the Formation of an Active Hexamer

The bacteriophage T7 gene 4 protein, like a number of helicases, is believed to function as a hexamer. The amino acid sequence of the T7 gene 4 protein from residue 475 to 491 is conserved in the homologous proteins of the related phages T3 and SP6. In addition, part of this region is conserved in DNA helicases such as Escherichia coli DnaB protein and phage T4 gp41. Mutations within this region of the T7 gene 4 protein can reduce the ability of the protein to form hexamers. The His475-->Ala and Asp485-->Gly mutant proteins show decreases in nucleotide hydrolysis, single-stranded DNA binding, double-stranded DNA unwinding, and primer synthesis in proportion to their ability to form hexamers. The mutation Arg487-->Ala has little effect on oligomerization, but nucleotide hydrolysis by this mutant protein is inhibited by single-stranded DNA, and it has a higher affinity for dTTP, suggesting that this protein is defective in the protein-protein interactions required for efficient nucleotide hydrolysis and translocation on single-stranded DNA. Gene 4 protein can form hexamers in the absence of a nucleotide, but dTTP increases hexamer formation, as does dTDP, to a lesser extent, demonstrating that the protein self-association affinity is influenced by the nucleotide bound. Together, the data demonstrate that this region of the gene 4 protein is important for the protein-protein contacts necessary for both hexamer formation and the interactions between the subunits of the hexamer required for coordinated nucleotide hydrolysis, translocation on single-stranded DNA, and unwinding of double-stranded DNA. The fact that the gene 4 proteins form dimers, but not monomers, even while hexamer formation is severely diminished by some of the mutations, suggests that the proteins associate in a manner with two separate and distinct protein-protein interfaces.

Evolution of T4-related phages

Much progress has been made in understanding T-even phage biology in the last 50 years. We now know the entire sequence of T4, encoding nearly 300 genes, only 69 of which have been shown to be essential under standard laboratory conditions; no specific function is yet known for about 140 of them. The origin of most phage genes is unclear, and only 42 genes in T4 have significant similarity to anything currently included in GenBank. Comparative analysis of related phages is now being used to gain insight into both the evolutionary origins and interrelationships of these phage genes, and the functions of their protein products. The genomes of phages isolated from Tbilisi hospitals, Long Island sewage plants, the Denver zoo, and Khabarovsk show basic similarity. However, these phages show substantial insertions and deletions in a number of regions relative to each other, and closer investigation of specific sequences often reveals much more complex relationships. There are only a few cases in T4-related phages in which there is evidence for evolution through DNA duplication. These include the fibrous products of genes 12, 34, and 37; head proteins gp23 and gp24; and the Alt enzyme and its downstream neighbors. T4 also contains 13 apparent relatives of group I and group II intron homing endonucleases. Distal portions of the tail fibers of various T-even phages contain segments closely related to tail-fiber regions of other DNA coliphages, such as Mu, P1, P2, and lambda. Horizontal gene transfer clearly emerges as a major factor in the evolution of at least the tail-fiber regions, where site-specific recombination probably is involved in the exchange of host-range determinants.

Secondary Structure Model of the Last Two Domains of Single-stranded RNA Phage Qβ

We have determined the nucleotide sequence of three positive single-stranded RNA coliphages and have used this information, together with the known sequences of the related phages Qβ and SP, to construct a secondary structure model for the two distal domains of Qβ RNA. The 3′ terminal domain, which is about 100 nucleotides long, contains most of the 3′ untranslated region and folds into four short, regular hairpins. The adjacent 3′ replicase domain contains about 1100 nucleotides. Hairpins in this protein-coding domain are much longer and more irregular than in the 3′ untranslated region. Both domains are defined by long-distance interactions. The secondary structure is not a collection of hairpin structures connected by single-stranded regions. Rather, the RNA stretches between the stem-loop structures are all involved in an extensive array of long-distance interactions that contract the molecule to a rigid structure in which all hairpins are predicted to have a fixed position with respect to each other. A general feature of the model is that helices tend to be organized in four-way junctions with little or no unpaired nucleotides between them. As a result, there is a potential for coaxial stacking of adjacent stems. The essential features of the model are supported by the SInuclease cleavage pattern.Viral RNA sequences are strongly constrained by their coding function. As a result, structural evolution by simple base-pair substitutionis not always possible, as this usually requires the juxtaposition of the codon wobble positions in stems. Rather, we often observe co-ordinate base substitutions that maintain the stem, but tend to change the position at which bulges or internal loops are found. Structures that differ in this way are apparently equally fit. Also, the relative position of hairpin loops can shift several nucleotides through an alignment based on maximal sequence identity i.e. amino acid homology. The fact that these structural irregularities do not occur at the 3′ untranslated region suggests indeed that the coding function of the RNA constrains the secondary structure. Hairpins with the stable tetraloop motif GNRA and UNCG or their complement are over-represented. This suggests their involvement in segregation of plus and minus strand. The genome of the coliphages contains a well-defined high affinity binding site for the coat protein, which serves to suppress replicase translation and also acts as a nucleation point in capsid formation. Close to the 3′ end we find an additional conserved helix that meets the described consensus criteria for coat-protein binding.