Genetics, Vol. 162, 1003-1018, November 2002, Copyright © 2002

Dissecting the Fidelity of Bacteriophage RB69 DNA Polymerase: Site-Specific Modulation of Fidelity by Polymerase Accessory Proteins

Anna Bebenek1,a,c, Geraldine T. Carvera, Holly Kloos Dressman2,a, Farid A. Kadyrova, Joseph K. Hasemanb, Vasiliy Petrovd, William H. Konigsberge, Jim D. Karamd, and John W. Drakea
a Laboratory of Molecular Genetics, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, North Carolina 27709-2233,
b Biometry Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, North Carolina 27709-2233,
c Institute of Biochemistry and Biophysics, Polish Academy of Science, 02-106 Warsaw, Poland,
d Department of Biochemistry, Tulane University Health Sciences Center, New Orleans, Louisiana 70112-2699
e Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut 06510-3219

Corresponding author: John W. Drake, Room E-344, National Institute of Environmental Health Sciences, P.O. Box 12233, Research Triangle Park, NC 27709-2233., drake{at}niehs.nih.gov (E-mail)

Communicating editor: S. LOVETT


*  ABSTRACT
*TOP
*ABSTRACT
*MATERIALS AND METHODS
*RESULTS
*DISCUSSION
*LITERATURE CITED

Bacteriophage RB69 encodes a replicative B-family DNA polymerase (RB69 gp43) with an associated proofreading 3' exonuclease. Crystal structures have been determined for this enzyme with and without DNA substrates. We previously described the mutation rates and kinds of mutations produced in vivo by the wild-type (Pol+ Exo+) enzyme, an exonuclease-deficient mutator variant (Pol+ Exo-), mutator variants with substitutions at Tyr567 in the polymerase active site (PolM Exo+), and the double mutator PolM Exo-. Comparing the mutational spectra of the Pol+ Exo- and Pol+ Exo+ enzymes revealed the patterns and efficiencies of proofreading, while Tyr567 was identified as an important determinant of base-selection fidelity. Here, we sought to determine how well the fidelities of the same enzymes are reflected in vitro. Compared to their behavior in vivo, the three mutator polymerases exhibited modestly higher mutation rates in vitro and their mutational predilections were also somewhat different. Although the RB69 gp43 accessory proteins exerted little or no effect on total mutation rates in vitro, they strongly affected mutation rates at many specific sites, increasing some rates and decreasing others.


THE DNA polymerase (gp43, encoded by gene 43) of the T4-related phages is responsible for the replication of the ~170-kb phage genome (KARAM and KONIGSBERG 2000 Down). Lacking DNA mismatch repair, these phages achieve the high fidelity of replication, ~2 x 10-8/bp, through the combined polymerase (Pol) and exonuclease (Exo) activities of gp43 (DRESSMAN et al. 1997 Down; BEBENEK et al. 2001 Down). Like T4 gp43 and the gp43-like DNA polymerases of some archaeons, gp43 of the T4-related coliphage RB69 is a member of the B-family (Pol{alpha}-like) polymerases, which include the eukaryotic replicative polymerases -{alpha}, -{delta}, and -{epsilon} (BRAITHWAITE and ITO 1993 Down). The crystal structure of RB69 gp43 shows five distinct domains: N (N terminal), palm (seat of the Pol catalytic center), fingers, thumb, and Exo (seat of the exonuclease catalytic activity; WANG et al. 1997 Down). The recently reported structures of three B-family polymerases from archaea show marked similarities to RB69 gp43 (HOPFNER et al. 1999 Down; RODRIGUEZ et al. 2000 Down; HASHIMOTO et al. 2001 Down). Two additional RB69 gp43 crystallographic structures are now available: a complex with DNA with the primer terminus occupying the Exo site (SHAMOO and STEITZ 1999 Down) and a ternary complex with partially double-stranded template-primer DNA and an incoming dNTP (FRANKLIN et al. 2001 Down). The three structures provide a basis for testing structure-function relationships in the fidelity of this polymerase.

As with other DNA polymerases, the Pol catalytic site of RB69 gp43 is composed of residues in the palm domain, thumb domain (which binds primer-duplex DNA), and fingers domain (which contains most of the side chains that bind the incoming nucleotide; FRANKLIN et al. 2001 Down). In unliganded gp43, the enzyme displays an open conformation. After binding a primer-template and a nucleotide, the resulting ternary complex displays a rotated fingers domain and changes in the conformation of the template DNA. This liganded state is referred to as the closed conformation (FRANKLIN et al. 2001 Down). Nearly all crystallographic structures for polymerases in such ternary complexes show the same rotation of the fingers domain toward the polymerase catalytic center (reviewed in KUNKEL and BEBENEK 2000 Down). When the fingers close down completely, a tight binding pocket is formed around the nascent base pair. The geometry of this pocket closely accommodates only a base pair with proper geometry and thus discriminates against mispairs. Therefore, although hydrogen bonding contributes strongly to base selection, geometry alone can provide reasonably correct nucleotide selection (GOODMAN 1997 Down; MORAN et al. 1997 Down). In addition, the polymerase is envisioned to test the geometry of the nascent base pair through specific hydrogen bonding between the enzyme and the pyrimidine O2 or purine N3 atoms in the minor groove of duplex DNA. Chemistry to produce primer extension depends strongly on the hydrogen-bonding potential of the 3'-terminal nucleotide (MORALES and KOOL 2000 Down) and primer extension may provide the signal for transferring the primer terminus to the Exo catalytic site for proofreading (JOHNSON 1993 Down). Both structural and biochemical studies suggest that the switch from primer extension to proofreading involves a conformational transition in the enzyme from a closed to an open (or Exo) mode (DOUBLIE et al. 1998 Down; SHAMOO and STEITZ 1999 Down).

Tyr567 in the crystal structure of closed RB69 gp43 is positioned near the Pol catalytic site at the end of helix P of the fingers domain. This tyrosine is conserved in all B-family polymerases and is a part of the highly conserved region III, KXXXNXYG (WANG et al. 1997 Down). The closed RB69 gp43 ternary structure suggests that Tyr567 makes a water-mediated hydrogen bond with a minor-groove acceptor on the template base of the terminal primer-template base pair. Thus, Tyr567 is proposed to be a part of the Pol binding pocket and to be an important residue for sensing proper geometry in a newly formed base pair (FRANKLIN et al. 2001 Down). We recently reported that substitution of Tyr567 with any of several amino acid residues having smaller side chains dramatically increases base-pair substitution errors in vivo but has little effect on frameshift mutagenesis (BEBENEK et al. 2001 Down). Kinetic measurements in vitro showed that the Y567A substitution (PolM enzyme) increases the rate of formation of some mispairs, but at the same time reduces some rates of mispair extension.

In the same study we examined the fidelity of an RB69 gp43 whose Exo activity was inactivated by two alanine substitutions at critical metal-coordinating Asp residues. The Exo- mutant exhibited the expected strongly increased mutation rate and promoted several kinds of mutations, including base-pair substitutions and additions and deletions of single base pairs. The PolM and the Exo- mutators each increased mutation rates in vivo by ~500-fold. We also observed that the PolM Exo- double mutator increased average mutation rates by ~2300-fold (BEBENEK et al. 2001 Down). This result sharply contradicts the expectations of one simple model in which the Pol and Exo fidelity contributions operate independently, where the outcome would be a multiplicative rather than an additive increase. This observation and other evidence indicate that the Pol and Exo functions of RB69 gp43 interact with each other, but the mechanism of this interaction remains unexplained.

In the studies reported here, we explored the extent to which the fidelities of the same RB69 gp43 mutant enzymes measured in vitro mirror those determined previously in vivo. We used the well-described M13mp2 lacZ{alpha} in vitro system (BEBENEK and KUNKEL 1995 Down). This system monitors all 12 possible base mispairs and many kinds of base addition and deletion mutations in a variety of nucleotide sequence contexts and has been used to characterize the fidelities of numerous wild-type and mutant DNA polymerases (KUNKEL and BEBENEK 2000 Down). Our study had two major aims. First, we sought to determine how well the rates and patterns of mutagenesis ascertained in vitro reflect the same parameters ascertained in vivo. This important comparison, which has been performed only infrequently with other polymerases, assesses the advantages and drawbacks of the two types of approach to measuring mutation rates and specificities. Second, we sought to measure the effects of the polymerase accessory proteins upon fidelity. Such measurements have been performed a few times before but only for the overall rate of mutation in a reporter gene or for rates of reversion at a small number of sites, rather than for rates at a large number of specific DNA sites.

Gp43 is poorly processive (DAS and FUJIMURA 1979 Down; NEWPORT et al. 1980 Down) except in the presence of the phage-induced sliding-clamp (a gp45 trimer), clamp loader (a 4:1 complex of gp44:gp62), and the ssDNA-binding protein (SSB, gp32; YOUNG et al. 1992 Down; NOSSAL 1994 Down). In contrast to our studies in vivo, the in vitro approach allowed us to compare the replication fidelities of different RB69 gp43 mutants acting either alone or in the presence of gp32, gp44/gp62, and gp45.

Here we report that the fidelities of the wild-type and mutator RB69 gp43's measured in vitro only partially resemble their fidelities measured in vivo. Error rates are generally somewhat higher in vitro than in vivo. Both the Pol+ Exo- and PolM Exo+ enzymes displayed mutator activities in vitro and, as in vivo, the PolM Exo- enzyme displayed only modestly higher mutator activity than that of either of its components. The Pol+ Exo- polymerase produced diverse base substitutions, single-base deletions, and large deletions. The PolM Exo+ polymerase produced mostly base substitutions, but, unlike in vivo, it also produced a substantial number of single-base deletions. The double-mutator PolM Exo- polymerase produced mostly base substitutions, but of types not predicted by its component mutators acting alone. In addition to producing single-base deletions, PolM Exo- also produced single-base additions, again unlike either of its component mutators. While the accessory proteins exhibited little or no impact on total mutation rates, they substantially increased mutation rates at some sites while decreasing rates at other sites. In addition to offering potential insights into factors that may determine the fidelity of DNA synthesis by gp43, these results indicate that fidelity analyses conducted in vitro may not agree fully with analyses conducted in vivo.


*  MATERIALS AND METHODS
*TOP
*ABSTRACT
*MATERIALS AND METHODS
*RESULTS
*DISCUSSION
*LITERATURE CITED

RB69 DNA polymerases and accessory proteins:
The production and purification of wild-type RB69 gp43 (Pol+ Exo+), D222A/D327A gp43 (Pol+ Exo-), Y567A gp43 (PolM Exo+), and Y567A/D222A/D327A gp43 (PolM Exo-) have been described (BEBENEK et al. 2001 Down). The RB69 accessory proteins gp45 and gp44/62 were prepared as described (YEH 1998 Down; YEH et al. 1998 Down). T4 gp32, which can support RB69 DNA replication (BORJAC 1997 Down), was prepared as described by KADYROV and DRAKE 2001 Down.

Gap-filling synthesis and forward-mutation assay:
Gap-filling reactions were performed with phage M13mp2 lacZ{alpha} gapped substrates prepared as described by BEBENEK and KUNKEL 1995 Down. Reactions with polymerase alone (25 µl) contained 1.6 nM DNA substrate, 25 mM Tris acetate, 10 mM magnesium acetate, 150 mM potassium acetate, 2 mM dithiothreitol, and 40 nM gp43. Reactions were started by adding all four dNTP's to final concentrations of 1 mM each. The mixtures were incubated for 15 min at 37° and were terminated by adding ethylenediamine tetraacetic acid (EDTA) to 15 mM.

Reactions with polymerase plus accessory proteins were conducted in the same mixture supplemented with 1 mM ATP, 6 nM gp44/gp62 complex, 80 nM gp45 (as a monomer) and 5.6–11.2 µM gp32. Reaction mixtures were preincubated for 30 sec at 37° to allow the accessory proteins to assemble on the gapped DNA substrate before polymerase was added to 20 nM. This mixture was incubated at 37° for another 30 sec to allow the polymerase to bind to the nascent replication complex and primer extension was then started by adding all four dNTP's to final concentrations of 1 mM each. After 10 min at 37° reactions were terminated by adding EDTA to 15 mM. A 20-µl sample from each reaction was analyzed on a 0.8% agarose gel to verify the extent of gap filling. The remaining 5 µl was diluted with 100 µl of distilled water and centrifuged twice through Microcon 30 columns, the DNA being recovered in a volume of 10–50 µl. Small portions (1–5 µl) of this DNA were subsequently electroporated into competent Escherichia coli MC1061 cells that were then plated to form M13 plaques. All plaques were counted and putative mutant plaques were scored as white or light blue against a background of dark-blue nonmutant plaques. Apparent M13 lacZ{alpha} mutants were picked and replated to verify their phenotypes and single isolated mutant plaques were resuspended in 100 µl of water. A representative set of mutants from each collection was sequenced to obtain mutational spectra.

Mutation-frequency and spectral data were usually compiled from several reactions, frequently with several transfections per reaction. Mainly because of small sample sizes, mutation frequencies usually varied by a fewfold (typically twofold or less) among samples for a given polymerase. Where appropriate, we call attention to mutation-frequency differences that we judge to be not strongly reliable due to underlying variation or small numbers of mutants scored.

DNA sequencing:
Prior to sequencing a mutant, the entire lacZ{alpha} gene was amplified by the polymerase chain reaction (PCR) using two 20-mer primers for PCR purchased from Research Genetics (Birmingham, AL). The forward primer was 5'-TAAGGGATTTTGCCGATTTC, the reverse primer was 5'-CAGTTTGAGGGGACGACGAC, and the length of the final product was 642 bp. The PCR consisted of 30 cycles of 1 min at 94°, 1 min at 58°, 1 min at 72°, with a final extension time of 10 min at 72° using Taq large-fragment polymerase (Display System Biotech TAQFL from PGC Scientifics). PCR products were purified with the QIAGEN (Valencia, CA) PCR purification kit. Sequencing was performed using the same primers and an ABI Prizm 377 automatic sequencer and dRhodamine terminator cycle sequencing kits (PE Applied Biosystems, Foster City, CA). Each mutation was identified by sequencing in both directions. When mutational spectra were constructed, only mutations known from historical experience to produce a mutant phenotype were included. Although silent mutations appearing singly are exceedingly rare, they are more common in mutants that also bear a phenotypically detectable mutation. The detectable base substitutions comprise those described by BEBENEK and KUNKEL 1995 Down as well as three others discovered by us (retaining the criterion that such a mutation must be observed as a single event at least twice): G -> C at 178, T -> C at 183, and C -> A at 189 using the numbering system of Fig 1.






View larger version (104K):
In this window
In a new window
Download PPT slide
 
Figure 1. Mutational spectra for RB69 DNA polymerases with and without accessory proteins. The 5' -> 3' sequence of the viral template strand of the lacZ{alpha} sequence in M13mp2 is shown from position -84 through +190 where +1 is the first transcribed base, the vertical line separates - and + bases, and the initiation codon is coded by ATG at 39–41; the last seven bases are not shown because they generated no mutations in this study. The direction of DNA synthesis is leftward from 190. +AP indicates the presence of accessory proteins and -AP indicates their absence during the reaction. Capital letters indicate base substitutions. The deletion of a single base is indicated by {Delta} and the addition of a single base by {nabla} except that the addition of a base unlike a template base is indicated by an arrow pointing between the two bases between which the insertion occurred and with the inserted base indicated below the arrow. The deletion or addition of two or more bases is indicated by - or + before the underlined capital letters. Larger deletions are described in Table 3 and complex mutations in Table 4. Significant asymmetries between the spectra with and without accessory proteins are indicated by ** for P < 0.05 or * for P < 0.10 and by displaying the more frequent kind of mutation in boldface. (A) Pol+ Exo+. (B) Pol+ Exo-. (C) PolM Exo+. (D) PolM Exo-. Note that the mutation shown as -TCG below the first 10 bases in the second line of template sequence in this spectrum could result from deletions of either TCG or CGT from the fugal sequence TCGT.

Statistical methods:
To determine whether the distribution of mutations across the 274 lacZ{alpha} bases displayed in the mutational spectra was similar in reactions with or without the polymerase accessory proteins, a chi-square statistic was calculated on the basis of the observed distribution of the mutations. However, because of the low expected mutation frequencies at each site, this statistic was not compared to standard chi-square tables. Instead, using the approach suggested by PIEGORSCH and BAILER 1994 Down, a Monte Carlo simulation study was performed, conditional upon the observed totals of mutations from reactions with or without accessory proteins as well as upon the total number of mutations at each site. The overall significance of the difference between reactions with or without accessory proteins with regard to the observed pattern of mutation frequencies was determined by comparing the observed chi-square statistic to chi-square statistics from 5000 simulated data sets. Pairwise comparisons of mutation frequencies at each site were made by two-sided Fisher's exact tests.


*  RESULTS
*TOP
*ABSTRACT
*MATERIALS AND METHODS
*RESULTS
*DISCUSSION
*LITERATURE CITED

Fidelities of RB69 DNA polymerase mutants:
Having characterized the fidelities of several variant DNA polymerases in vivo (BEBENEK et al. 2001 Down), we proceeded to conduct parallel analyses in vitro by measuring the accuracy of DNA replication during the filling of a 407-nucleotide (nt) single-stranded gap in bacteriophage M13mp2 RF DNA. The gap contained the mutation-reporter lacZ{alpha} as a template for DNA synthesis (BEBENEK and KUNKEL 1995 Down). This forward-mutation assay allows the detection of many base substitutions, most single-base additions and deletions, and many more complex mutations. Mutants are detected as light-blue or colorless M13 plaques against a background of dark-blue nonmutant plaques. The gap was filled by either the wild-type or a mutant RB69 gp43 in reactions either with or without the gp43 accessory proteins. T4 gp43 is poorly processive in the absence of the accessory proteins (DAS and FUJIMURA 1979 Down; NEWPORT et al. 1980 Down). Because RB69 gp43 differs from its T4 counterpart at ~39% of its residues (WANG et al. 1995 Down), we compared the processivities of the two proteins in vitro. Both enzymes displayed identical limited processivity, the predominant products being about 20 nt long but with some products of 100–140 nt (data not shown).

For each of the four gp43 constructs studied here, Table 1 describes mutation rates, Table 2 describes the kinds of mutations produced, and Fig 1 shows the mutational spectra. In the spectra, certain sites are highlighted because the accessory proteins either enhanced or inhibited mutagenesis at those sites; these are necessarily mutation-rich sites where significant differences between numbers of mutations produced with and without accessory proteins can be recognized. Error rates for specific classes of mutations can be uniquely well estimated in the lacZ{alpha} system because previous studies have established which changes at each site yield detectable mutations (BEBENEK and KUNKEL 1995 Down). Thus, error rates for each specific kind of mutation can be normalized to the number of informative targets, and these rates are provided in Table 5 and Table 6.


 
View this table:
In this window
In a new window

 
Table 1. Mutation frequencies (MF) and rates (µ) for wild-type and mutant RB69 DNA polymerases in vitro and in vivo


 
View this table:
In this window
In a new window

 
Table 2. Kinds of mutations produced in vitro by wild-type and mutant RB69 DNA polymerases without and with accessory proteins


 
View this table:
In this window
In a new window

 
Table 3. Characteristics of deletion mutations produced by wild-type and mutant RB69 DNA polymerases in vitro without and with accessory proteins


 
View this table:
In this window
In a new window

 
Table 4. Characteristics of complex mutations produced by wild-type and mutant RB69 DNA polymerases in vitro without and with accessory proteins


 
View this table:
In this window
In a new window

 
Table 5. Base-substitution rates (µb) for mutant RB69 DNA polymerases without and with accessory proteins


 
View this table:
In this window
In a new window

 
Table 6. Base addition/deletion rates for mutant RB69 DNA polymerases without and with accessory proteins

Below, we show results suggesting that the mutations in the Pol+ Exo+ (wild-type gp43) reaction products were produced during the growth of the M13 phage vector and the processing of the M13 DNA template and not by the RB69 polymerase itself. This background contributes modestly to mutations generated by the Pol+ Exo- polymerase, slightly to those produced by the PolM Exo+ polymerase, and insignificantly to those produced by the PolM Exo- polymerase. Significant contributions will be noted at the appropriate places. When we estimated such contributions, we pooled the Pol+ Exo+ results with and without accessory proteins because most of the mutations preceded the reactions in vitro.

Fidelity of the Pol+ Exo+ RB69 DNA polymerase:
The wild-type RB69 and T4 DNA polymerases were previously shown to be so accurate in copying lacZ{alpha} that their mutation frequencies (3–5 x 10-4) were at or slightly below the historical background frequency (5–7 x 10-4) for unfilled template DNA (DRESSMAN et al. 1997 Down). This background frequency simply reflects the intrinsic average mutation rate per base pair for phage M13 replication (DRAKE 1991 Down), but it can hinder analyses of the fidelities of very accurate polymerases. A similar result was encountered here: In reactions either with or without accessory proteins, the mutation frequency was similar to or indistinguishable from the historical background (Table 1). We do not consider the apparent difference between the mutation frequency without (~8 x 10-4) and with (~4 x 10-4) accessory proteins to be reliable because of the small numbers of mutants scored. The kinds of mutations observed (Table 2) and their spectra (Fig 1A) are indistinguishable in the presence or absence of accessory proteins, as expected from mutations that already existed in the template DNA prior to the reactions conducted in vitro. Among the 43 sequenced mutants, 26 were C -> T transitions, 5 were single-base deletions, and 5 were larger deletions, a distribution similar to one described previously (KUNKEL 1985 Down).

Fidelity of the Pol+ Exo- RB69 DNA polymerase:
A pair of alanine substitutions at two catalytic aspartic-acid residues (D222A/D327A) in the exonuclease domain of RB69 gp43 completely inactivates the exonuclease function (ABDUS SATTAR et al. 1996 Down; FRANKLIN et al. 2001 Down). This Exo- defect had little effect (<10% reduction) on DNA synthesis and phage production in vivo (BEBENEK et al. 2001 Down).

The exonuclease-deficient RB69 gp43 is a mutator polymerase in vitro, exhibiting a lacZ{alpha} mutation rate of 35 x 10-4 that is only a little larger (1.6-fold) than the rI forward mutation rate in vivo of 22 x 10-4 (Table 1). The accessory proteins do not significantly affect the total mutation rate. The Pol+ Exo+ background mutation frequency contributes ~23% of the mutations observed in the two Pol+ Exo- collections.

The mutations produced by the Pol+ Exo- polymerase in vivo are about 50% transitions, 25% transversions, and 25% additions or deletions of single base pairs (BEBENEK et al. 2001 Down). This polymerase produces roughly similar proportions of mutations in vitro and, in addition, some larger deletions. The specific kinds of base substitutions produced in vivo and in vitro are also similar (Table 2, Table 5, and Table 6) except that single-base additions were favored in vivo whereas single-base deletions predominate in vitro.

The only significant effect of the background contribution is to reduce the number of polymerase-induced C -> T transitions about twofold (from 45 and 58 to 24 and 32 without and with accessory proteins, respectively) and to reduce the numbers of transitions by about one-third (from 67 and 94 to 45 and 67 among similar numbers of mutations). The background contributed most of the C -> T mutations at positions -32, 75, 108, and 166.

The accessory proteins altered the distributions of mutations among types in several ways. When numbers of mutations are normalized to relative mutation rates with and without accessory proteins (Table 1) and to sample sizes (Table 2), the accessory proteins modestly increased the frequencies of transitions (by 12%) and transversions (by 32%). They decreased the overall frequency of single-base deletions by 36%, but increased the frequency most (by 69%) at single-base targets while increasing the frequencies modestly (P > 0.1) at a four-C run and strongly (by 540%) at a five-C run. They reduced the frequency of larger deletions by 86%.

The two spectra display four specific sites (at positions -36, -32, 82, and 132–136) in which the accessory proteins significantly altered the numbers of mutations (Fig 1B). At two of these sites the mutation rate is higher in reactions with accessory proteins and at two it is lower. (When this spectrum was adjusted by subtracting the calculated contributions of the Pol+ Exo+ background, the resulting spectrum displayed a similar pattern.) Thus, although the accessory proteins had no effect on the gross mutation rate for this polymerase, they substantially altered mutational specificity in class- and site-specific ways.

Fidelity of the PolM Exo+ RB69 DNA polymerase:
The Y567A substitution at the Pol site is strongly mutagenic in vivo, generating transitions almost exclusively (BEBENEK et al. 2001 Down). The PolM Exo+ enzyme is also a mutator in vitro (Table 1) and is sufficiently strong that the Pol+ Exo+ background contributes only ~10% of the mutations.

The PolM Exo+ mutation rate in vitro (100 x 10-4) is about five times higher than the corresponding rate in vivo (21 x 10-4). However, as discussed later, this difference is likely to have been augmented by our use of dNTP concentrations high enough to partially inhibit proofreading. The mutational propensities of the PolM Exo+ gp43 in vitro differ in three ways from those observed in vivo (Table 2). First, the mutator activity of this enzyme is about threefold stronger than that of the Pol+ Exo- enzyme in vitro, whereas the two mutators were identically strong in vivo. Second, a substantial minority of the mutations produced in vitro are single-base deletions, in contrast to almost none in vivo. Third, the favored transition is G·C -> A·T in vivo but T -> C in vitro. Larger deletions were absent from the PolM Exo+ spectrum both in vivo and in vitro.

The only significant impact of subtracting the spontaneous background from the values in Table 2 is to reduce the numbers of C -> T transitions. These are reduced from 17 to 7 without accessory proteins and from 8 to 0 with accessory proteins.

The accessory proteins appear to slightly increase (by 27%) the overall PolM Exo+ mutation frequency (Table 1), but this increase is within the variability of the individual measurements pooled to produce the values in Table 1. The accessory proteins also seem to have little effect on the proportions of mutations of various kinds (Table 2). Nevertheless, they significantly modify site-specific mutation rates at five positions (-58, 71, 112, 147, and 149), where rates were increased at three positions and decreased at two (Fig 1C). The spontaneous background has almost no effect on these two mutational spectra.

Fidelity of the PolM Exo- RB69 DNA polymerase:
When examined in vivo (BEBENEK et al. 2001 Down), this double-mutator gp43 was only a little stronger (an estimated 2300-fold increase over the Pol+ Exo+ mutation rate) than its component PolM and Exo- mutators (each with ~500-fold increases). This enzyme exhibited a significantly higher mutation rate in vitro in the presence of accessory proteins (~400 x 10-4) than in their absence (~260 x 10-4; Table 1) and both of these values exceed the range of mutation rates (74–120 x 10-4) that we estimated using the rI system in vivo. The Pol+ Exo+ background does not contribute significantly to mutations produced by the PolM Exo- enzyme. Because of the high lacZ{alpha} background mutation frequency contributed by growth in M13, the Pol+ Exo+ mutation frequency could not be determined and mutator factors therefore cannot be expressed as fold increases over the Pol+ Exo+ value. However, as observed previously in vivo, the PolM Exo- mutation rate in vitro is only a fewfold higher than the sum of its component single mutations regardless of the extent to which the PolM Exo+ rate may have been overestimated because of the high concentration of dNTP's used here (discussed below).

As previously observed in vivo, the PolM Exo- enzyme is a base-substitution mutator that produces many more transitions than transversions. While its frequencies of A -> G transitions are low, as also seen with both single mutators, the frequencies of the other three transitions are not well predicted by their frequencies in the single mutators.

Approximately 14% of the mutations produced by the PolM Exo- enzyme in vitro are composed of losses or additions of single bases, compared with about 2% in vivo. The accessory proteins produced at most modest effects on frequencies of base substitutions, but perhaps diminished frequencies of single-base additions and of larger deletions (Table 2). On the other hand, their effects on site-specific mutation rates are striking. The accessory proteins decreased mutation rates at positions -68, -66, 87, and 89 and increased mutation rates at positions 118, 139, and 183 (Fig 1D). Because the overall mutation rate is increased ~1.5-fold by the accessory proteins, it is likely that increases occurred at many other sites but were not detectable because of the small number of mutations at most sites.


*  DISCUSSION
*TOP
*ABSTRACT
*MATERIALS AND METHODS
*RESULTS
*DISCUSSION
*LITERATURE CITED

Mutational targets:
The targets we used in vivo and in vitro to conduct these RB69 gp43 fidelity studies display similarities as well as differences. The T4 rI target used in vivo contains 291 translated bases while the lacZ{alpha} target used in vitro contains 281 bases, of which 197 are transcribed and 159 are translated. The mRNA strand of T4 rI has the base composition A = 91, T = 94, G = 59, and C = 47 and an A·T content of 63.5%, which is close to the 64.7% calculated for the whole T4 genome (http://www.ncbi.nlm.nih.gov:80/cgi-bin/Entrez/getff?gi=15081&db=Genome). The lacZ{alpha} mRNA strand has the base composition A = 68, T = 68, G = 62, and C = 83, or 48.4% A·T, which is close to the 49.2% A·T for the E. coli genome from which it originated (BLATTNER et al. 1997 Down) but with a modest excess of C over G. Both targets contain numerous single-base runs of various lengths that are prone to frameshift mutagenesis, and both targets contain numerous adjacent and distant repeats that are prone to deletion mutagenesis. The main advantages of the rI target are that it reflects biochemical realities in vivo and that, for the wild-type polymerase, it has the low average spontaneous mutation rate per base pair expected for a T4 sequence, ~2 x 10-8 (DRAKE 1991 Down). The main advantages of the lacZ{alpha} target are that it differentiates between the mutabilities of each base in a base pair and that a historical record of detectable base substitutions is available.

Mutation rates in vivo and in vitro:
Table 1 shows that the mutation rates of the three mutator RB69 gp43's are somewhat higher in vitro than in vivo, by 1.6-fold for the Pol+ Exo- enzyme, by 5-fold for the PolM Exo+ enzyme, and by ~4-fold for the Pol+ Exo- enzyme in the presence of accessory proteins. The Pol+ Exo- difference is small and probably insignificant. The PolM Exo- difference is uncertain because of the indirect method used to estimate the rate in vivo and the disparate values obtained with and without accessory proteins in vitro, but is likely to be significantly higher in vitro. The PolM Exo+ difference may have been overestimated because we used dNTP's at 1 mM, which is high enough to inhibit proofreading activity to some extent. Proofreading is generally inhibited by high dNTP concentrations (such as 1 mM vs. 10 µM) in DNA polymerases with 3'-exonucleases such as T4 gp43 (CLAYTON et al. 1979 Down), E. coli Pol I (KUNKEL et al. 1981 Down; BEBENEK et al. 1990 Down), and E. coli Pol III holoenzyme (FERSHT 1979 Down; PHAM et al. 1998 Down). However, MATHEWS and SINHA 1982 Down estimated the effective dTTP concentration at the compartmentalized replication fork to be ~220 µM, ~5-fold lower than the concentration used in this study. More importantly, the kinds of mutations produced by the PolM Exo+ enzyme are not at all what might be expected from a strong contribution by a phenotypically Exo- condition. In one of several examples, in the absence of accessory proteins, the PolM Exo+ enzyme produced no deletions of >=2 bases, whereas 23(147/159) {cong} 21 were expected from the results with the Pol+ Exo- enzyme. Thus, the PolM Exo+ difference in vitro would probably decrease but not disappear at lower dNTP concentrations. In summary, the mutation rates observed with the RB69 gp43 mutators are only modestly higher in vitro than in vivo.

These modest differences in mutation rates cannot be adequately explained by low polymerase processivity, because including the accessory proteins markedly affected rates at numerous specific sites. We offer five possible explanations for these observations.

  1. Additional undiscovered replication accessory proteins might act in vivo to increase fidelity. This conjecture is mildly inconsistent with the striking lack of mutator mutations in genes other than those encoding the polymerase and its accessory proteins (DRAKE and RIPLEY 1994 Down). Note, however, that the functions of roughly half of the {cong}280 T4 genes remain undescribed (KUTTER et al. 1994 Down).

  2. Conditions in vitro might be inadequate to support optimal fidelity. For instance, the cationic composition of our reaction mixture is dominated by acetate, whereas it is likely to be glutamate in vivo. In addition, concentrations of anionic and cationic species and even replication proteins in vivo are not well described for the putative metabolic compartments in which T4 (and presumably RB69) DNA replication occurs (MATHEWS and ALLEN 1983 Down; GREENBERG et al. 1994 Down).

  3. As suggested by PHAM et al. 1998 Down, a replication fork with coupled synthesis of the leading and lagging strand might be intrinsically more accurate than half a fork. Gap filling in vitro certainly differs significantly from coupled, helicase-driven leading-strand synthesis and primer-dependent lagging-strand synthesis in vivo.

  4. A DNA polymerase might replicate a noncognate template (with which it did not coevolve) less accurately than a cognate template.

  5. The lacZ{alpha} reporter sequence may, for unknown reasons, be intrinsically more mutable than the rI reporter sequence.

Two other comparisons have been made of forward-mutation systems in vivo vs. in vitro. In one, FUJII et al. 1999 Down studied the E. coli rpsL gene moved to a plasmid. Whereas the mutation frequency of this ~500-base reporter gene was similar to that of lacZ{alpha} when on M13, it was far lower (~0.007 x 10-4) when plasmid borne, but close to that expected for the mutation rate of E. coli, ~5 x 10-10/bp (DRAKE 1991 Down). Rolling-circle replication by an in vitro mixture containing the Pol III holoenzyme generated an rpsL mutation frequency of ~2 x 10-4, a 300-fold higher frequency that must reflect at the least the lack of DNA mismatch repair in vitro. The predominant mutations produced in vitro (in the absence of mismatch repair) were additions and deletions of base pairs, whereas base-pair substitutions predominated in vivo. The base-substitution mutational spectra were very different in vivo and in vitro on the same reporter sequence, but the frameshift mutational spectra were rather similar. In the other comparison, the E. coli lacI gene was the reporter sequence (SCHAAPER and DUNN 1987 Down). In this case, the relevant rate in vivo is in a genetic background defective in DNA mismatch repair, in which case fidelity is determined primarily by the holoenzyme. Applying the method of DRAKE 1991 Down yields a mutation rate in vivo of 6.6 x 10-6. The rate in vitro using the Pol III holoenzyme and combining the results for 10 and 50 µM dNTP's was ~900 x 10-6 based on a small and uncertain difference between the mutation frequency with and without enzyme (PHAM et al. 1998 Down). Even higher rates were obtained in vitro when proofreading was absent because only the {alpha} (polymerase) subunit of Pol III was used (MO and SCHAAPER 1996 Down) or because synthesis was conducted with the holoenzyme using a proofreading-impaired mutant together with 1 mM dNTP's, which also impairs proofreading (PHAM et al. 1998 Down); the rate with {alpha} alone was 16,600 x 10-6 and with proofreading-impaired holoenzyme was 18,700 x 10-6, or ~20-fold higher than that with the proofreading-proficient holoenzyme. While the ratio of base substitutions to frameshift mutations was similar in vivo and with the wild-type holoenzyme in vitro at low dNTP concentrations, both of the conditions in vitro that impaired proofreading enhanced frameshift mutagenesis by about an order of magnitude more than they enhanced base-substitution mutagenesis. Thus, in both the rpsL and the lacI studies, mutation rates were far higher in vitro than in vivo, in contrast to the small differences seen in this study.

Mutation specificities:
The base-pair substitution hotspots produced by all four of the RB69 gp43 constructs in vivo tended to occur at certain specific G·C-rich 6-mers and especially at GG/CC dimers (BEBENEK et al. 2001 Down). This is a particularly striking tendency because of the low G·C content of the T4 genome. Like most of the E. coli genome, the lacZ{alpha} target is roughly half (G + C) but with some bias toward C. In vitro, the three spectra generated by mutator gp43's included 47 different sites with four or more base substitutions. These sites were taken to be much more mutable than average. (In this analysis, the Pol+ Exo- spectrum was first adjusted for the contribution from the Pol+ Exo+ background.) Unlike the hotspots produced in vivo, these 47 sites revealed no bias for either a G or a C to occupy either or both neighboring bases ({chi}2 = 3.4, P {cong} 0.4). However, these 47 sites comprised 0 A, 22 T, 15 G, and 10 C residues. This distribution is very different from a random sample of lacZ{alpha} bases ({chi}2 = 24.5, P < 0.0001): The number of hotspots at a G or a C is close to the expectations from random sampling, while the number at an A is much lower and the number at a T is much higher. We examined the distribution of T residues neighboring these hotspots but observed little bias ({chi}2 = 4.65, P {cong} 0.2). Thus, lacZ{alpha} hotspots generated by the mutator RB69 gp43's studied here occur preferentially at T sites, but not preferentially at TT or TTT sites. WANG et al. 2002 Down have recently described a process by which E. coli Pol II produces an excess of mutations at A·T-rich sites, but it is unclear whether this applies to the present system because the number of hotspots at A sites is relatively low while the number at T sites is relatively higher.

While the sequence determinants of base-substitution hotspots are usually either unknown or poorly understood, one model has been particularly helpful. In this model, forward misalignment of the primer terminus by one base, followed by correct incorporation, followed by realignment generates a base mispair that is then a candidate for extension (KUNKEL and SONI 1988 Down; BOOSALIS et al. 1989 Down; reviewed in BEBENEK and KUNKEL 2000 Down). Several of our base-substitution hotspots (particularly in PolM Exo+ and PolM Exo- backgrounds) contained mutations of the form XYY -> XXY (such as CTT -> CCT), which can be a signature of base substitution initiated by transient misalignment. However, we did not perform the tests that are required to establish the operation of this mechanism.

The sequence determinants of base addition and/or deletion hotspots are somewhat better understood. The first insight was the slippage model of STREISINGER et al. 1966 Down, which applies to events occurring within short homonucleotide runs or repeated short sequences. In this model, a primer base or short repeat misaligns backward or forward while retaining conventional base pairing, and primer extension then locks in the addition or deletion. Most frameshift hotspots in many spectra are of this type. In the second model, which is not confined to sequence repeats, the usual disinclination to extend a mispair can favor realignment of a mispaired primer-terminal base to achieve correct pairing provided the adjacent base is its complement, producing a readily extended primer and generating a single-base addition or deletion (with deletions strongly favored) (KUNKEL and SONI 1988 Down; BEBENEK and KUNKEL 1990 Down; reviewed in BEBENEK and KUNKEL 2000 Down). In the third model, active-site misalignment (which is also not confined to sequence repeats), the template base misaligns in a way that allows an incoming dNTP to pair with the next template base; if chemistry occurs and is followed by primer extension, a single-base deletion results (KUNKEL 1986 Down). This model is supported by both kinetic and structural studies (EFRATI et al. 1997 Down; HASHIM et al. 1997 Down; LING et al. 2001 Down). In this study, we observed numerous deletions of a single base in a nonrun environment (particularly of G and particularly in the Pol+ Exo- background). The high frequency of single-base deletions at positions 120–122 (CTT) and 136–139 (CTTT) in the PolM Exo+ background might have been promoted by the misinsertion of G opposite the 5' T. Again, however, we did not perform tests to establish the mode of frameshift mutagenesis.

We detected a total of 25 deletions of from 7 to 436 bases (Table 3). Of these, 18 appeared once, 1 appeared twice, and 1, a 317-base deletion, appeared five times and was a deletion frequently observed in this system (KUNKEL 1985 Down). Most of these deletions arose between direct repeats of 2–9 bases. These repeats are G·C-rich, 75% on average vs. 52% for lacZ{alpha}, an observation suggesting that distantly misaligned replication intermediates are more stable the higher the melting temperature of the misaligned segment. The 317-base deletion may have been promoted not only by the direct CCCGC repeat but also by a contribution from a short palindromic repeat (KUNKEL and SONI 1988 Down). Except for four deletions in the wild-type spectrum that may not have been generated in vitro, all the rest were generated by an Exo- polymerase. This observation suggests that large deletions generated in vitro by RB69 gp43 either are produced by the polymerase and are well proofread or are produced when the primer terminus partitions to the defective Exo site but returns to a distant site more frequently than from an Exo+ site. In the second mechanism, accurate realignment of a primer terminus would require a normal Exo site, because of either the energetics of cleavage or the correct primer orientation. Perhaps the most interesting observation concerning the deletions of >=2 bases is that their frequency was sharply reduced by the accessory proteins. For a short sequence to misalign with a distant repeat, several bases of the primer terminus must separate from the template strand and must then diffuse sufficiently to discover a distant unpaired complement. It would not be surprising if such a process were inhibited by a processivity factor and/or by gp32.

We observed four complex mutations of two types (Table 4). In one type, two mutations associated with a C-run were generated by a Pol+ Exo- polymerase and correspond to the replacement of two adjacent bases by one different base. These mutations could have arisen by any of several possible coupled mispair-and-slippage errors. In the second type, two mutations generated by a PolM Exo+ polymerase represent complicated palindrome-associated events, first recognized by GLICKMAN and RIPLEY 1984 Down, that may involve strand-switching templating errors during primer extension.

Certain kinds of mutations were produced more frequently in vitro than in vivo. These included single-base deletions by all three polymerases, larger deletions by the Pol+ Exo- polymerase (although most of these were prevented by the accessory proteins), and transitions at A and T sites relative to those at G and C sites by the PolM Exo+ polymerase (at least in part because of the hotspot biases described above). In contrast, the Pol+ Exo- polymerase produced many fewer single-base additions in vitro than in vivo with or without accessory proteins. The higher mutation rates observed in vitro vs. in vivo cannot be explained simply by the additional classes of mutations sometimes produced in vitro. The litany of possible reasons for these differences is the same as that invoked above for the differences in mutation rates, but here pertaining to specific mutational pathways.

The PolM-Exo- interaction:
The double mutator displayed two unanticipated traits in vivo: Its mutational specificity was not fully predicted by either of its two component mutators, and its mutator strength was only about twofold greater than the sum of the mutator strengths of its components. The same two traits are displayed in vitro: Both G -> A transitions and single-base additions are more frequent than those with either of the single mutators, and the strength of the double mutator is again only about twofold greater than the sum of the mutator strengths of its components. In one simple view of replication fidelity, proofreading removes a constant fraction of nascent mutations of various types. This view predicts that the double mutator should be far stronger than observed, even though this cannot be specifically predicted because we lack a Pol+ Exo+ value. The modest increase in the strength of the double mutator does not seem to be due to the efficient extension of mispaired primer termini because the PolM enzyme extends mispairs at least as poorly as does the Pol+ enzyme (BEBENEK et al. 2001 Down), but several other explanations are tenable: (i) It might reflect an effect of the Y567 Pol mutation on the partitioning of the primer terminus to the Exo site; (ii) it might result from an effect of the Exo deficiency on fidelity at the Pol site; (iii) if this gp43 usually dissociated from a mispair with subsequent loss of the DNA substrate, many mutations might be lost to our measurement system; or (iv–v) conversely, if this gp43 usually stalled irreversibly at a mispair, at either the Pol or the Exo binding sites, many mutations could be lost.

The impact of accessory proteins on polymerase fidelity:
The effect of the accessory proteins on the fidelity of DNA replication in vitro has been the subject of several studies:

  1. The T7 DNA polymerase processivity factor is a tightly bound host protein, thioredoxin. The T7 holoenzyme is too accurate to produce a reliable signal in the lacZ{alpha} system, but an Exo- holoenzyme did produce a signal (KUNKEL et al. 1994 Down). In reversion tests, the presence of thioredoxin decreased base-substitution frequencies by about threefold and decreased -1 and -2 base deletions substantially in nonrun contexts, but strongly increased additions of one or more bases within runs.

  2. As with T7 DNA polymerase, the mitochondrial DNA polymerase-{gamma} accessory subunit increased fidelity for some base substitutions on synthetic oligonucleotide primer templates (JOHNSON and JOHNSON 2001 Down). A similar result was obtained by LONGLEY et al. 2001 Down, who showed that fidelity increased, not during mispair formation but during mispair extension. In extensive tests with the lacZ{alpha} system, the overall fidelity of wild-type Pol{gamma} was unaffected by the accessory subunit, although the subunit promoted the specific template·primer mispair T·G with the wild-type enzyme. The overall fidelity of an Exo- Pol{gamma} was decreased threefold by the accessory subunit, mainly by increases in the formation of the specific mispair A·A and in single-base deletions in nonruns or in runs of size 2. In runs of size 4 or 5, the accessory subunit considerably decreased the rate of single-base deletions.

  3. The fidelity of mammalian DNA polymerase-{delta} on synthetic oligonucleotide primer templates was decreased considerably by the inclusion of its processivity clamp, proliferating cell nuclear antigen (MOZZHERIN et al. 1996 Down).

  4. When the fidelity of the E. coli Pol III {alpha} (polymerase) subunit alone (MO and SCHAAPER 1996 Down) was compared with that of a holoenzyme that was defective in proofreading because of the combination of a genetic defect in the proofreading subunit (mutD5) and a high concentration (1 mM) of dNTP's (PHAM et al. 1998 Down), there was no difference in the overall mutation rate. A small decrease in the relatively high frameshift rate (primarily single-base deletions) was offset by a tripling of the base-substitution rate.

T4 gp43 removed preformed terminal mismatches more efficiently when gp32, gp45, and gp44/gp62 were present (SINHA 1987 Down). However, we cannot compare these results with those reported here because the high lacZ{alpha} background in our system conceals the specificity of the Pol+ Exo+ enzyme. When T4 gp43 infidelity was measured as dNTP -> dNMP turnover using homopolymer primers and templates, the accessory proteins improved fidelity (TOPAL and SINHA 1983 Down).

The study of the effects of accessory proteins on fidelity that is most relevant to the current report used a T4 Pol+ Exo- gp43 and the lacZ{alpha} system (KROUTIL et al. 1998 Down). Adding the gp45 clamp increased the forward-mutant frequency from 24 to 47 per 104, in apparent contrast to our observation of an insignificant decrease from 29 to 25 per 104 (Table 1). The frequency of large deletions was decreased ~10-fold, a factor similar to that observed by us (Table 2). Adding the gp45 clamp increased the TGA revertant frequency from 4 to 10 per 105, an increase distributed irregularly over several of the possible mispairs. However, the wild-type precursor (TGG) to their TGA reporter at bases 87–89 is particularly prone to accessory-protein effects in the RB69 system (Fig 1B, Fig C, and Fig D), exhibiting both increases and decreases. Adding the gp45 clamp had no significant effect on the frequency of single-base deletions in a TTTTTTT run. Kinetic parameters for misincorporation were similar with and without the accessory proteins. A procedural difference should be noted between the T4 and the RB69 studies. In the T4 studies, most of the comparisons were between reactions with and without gp45, the other accessory proteins being present throughout, whereas in the RB69 studies tabulated above the comparison was between gp43 alone and gp43 with all four accessory proteins. However, this procedural difference seems unlikely to be profound, because we also measured mutant frequencies (but did not sequence mutants) in reactions lacking only gp45. For comparison with the values in Table 1, the uncorrected mutant frequencies for the four polymerases (in the same order) were 4.4, 20, 71, and 275 per 104. In each case, the frequency with gp32 and gp44/gp62 present was very similar with and without gp45.

Including the SSB, the processivity clamp and the clamp-loading proteins in our reactions had little or no effect upon gross mutation rates (Table 1). Although we did not measure the impact of the RB69 accessory proteins on processivity, other studies have examined this matter with the related T4 proteins, both in the context of leading-strand and lagging-strand synthesis (KADYROV and DRAKE 2001 Down and references therein) and during gap filling, all gaps being filled in 15 sec in their presence and none in their absence (KROUTIL et al. 1998 Down). In addition, the dramatic impact of the RB69 accessory proteins on site-specific mutation rates argues strongly for their participation in polymerase fidelity. There were only a few convincing effects on rates of specific types of base substitutions, additions or deletions (Table 2, Table 5, and Table 6), such as modestly increased frequencies of T -> C transitions and sharply decreased frequencies of larger deletions in the two Exo- mutators. However, the mutational spectra (Fig 1) revealed numerous specific sites where rates were affected. The striking aspect of these effects is that site-specific rates were sometimes increased and sometimes decreased, with no obvious preference for the direction of the bias.

Computer simulations first established the statistical significance of the differences in the overall distributions of mutation frequencies with and without accessory proteins for each of the three RB69 gp43's, but did not identify the particular sites at which these differences might be most notable. (As expected and noted previously, there were no significant differences with the wild-type gp43.) Subsequent Fisher's exact tests then identified specific hotspots showing evidence of such differences. Because of the large number of sites examined, we expected occasional false positives, that is, P values below 0.10 or 0.05 that can be attributed to chance alone. Sites having few or no mutations do not contribute false positives because even the most extreme distribution of mutations possible at such sites would not be statistically significant. For each of the four polymerases studied here, between 21 and 27 sites could have produced a false-positive outcome on the basis of the frequency of mutations observed at those sites. Calculations based on the distribution of the Fisher's exact test statistic revealed that for each polymerase, we would expect to find approximately one instance of P < 0.05 and two of P < 0.10 by chance alone. This compares with the two to five instances of P < 0.05 and four to seven of P < 0.10 actually observed for each of the variant polymerases. We conclude that while it is possible that a few of the site-specific differences in mutation frequency associated with the presence or absence of accessory proteins may be chance occurrences, the majority appear to reflect real effects upon both base substitutions and base deletions. Many other such biases probably occur but are not detected with our sample sizes.

Although our evaluation of each site was based on all mutations detected at that site regardless of type, in most instances a single type of mutation accounted for the statistically significant bias. The biases occur for a variety of base substitutions (at least the transitions G -> A, T -> C, and C -> T and at least the transversions T -> A and G -> T) and base deletions (at least -T, -C, and -G). For most kinds of mutations for which three or more asymmetries were recorded, the biases were in both directions. Of the 16 different biased sites we observed, none occurred more than once among the three mutator spectra. There was a clear tendency for a biased site to be flanked by a pyrimidine (P {cong} 0.01) but C (12) and T (9) were represented among the flanking pyrimidines in nearly the same proportion (83:68) as in the lacZ{alpha} target. However, structural information about the nature of the interactions among DNA, polymerase, clamp, and SSB are far too sparse to provide an explanation for the bias toward flanking pyrimidines.

In summary, the RB69 DNA-replication accessory proteins have at most small overall effects on mutation rates in vitro but strongly inhibit or promote mutation at specific sites. (This is another instance in which conclusions drawn from small targets, as in most reversion tests, may be misleading and should be confirmed in a forward-mutation test.) Evolution has optimized the interactions among the proteins of DNA metabolism to achieve an optimal mutation rate perhaps driven by the deleterious effects of most mutations and the costs that must be paid to keep mutation rates low (DRAKE 1991 Down). The huge variety of local DNA sequences that interact with the replication proteins at and near the primer terminus results in high mutation rates at some sites and low rates at others. Our results may offer an opportunity to probe the specific roles of accessory proteins in causing such variability.


*  FOOTNOTES

1 Present address: Institute of Biochemistry and Biophysics, Polish Academy of Science, 02-106 Warsaw, Poland. Back
2 Present address: Department of Genetics, Box 3054, Duke University Medical Center, Durham, NC 27710. Back


*  ACKNOWLEDGMENTS

We thank Tom Kunkel and Roel Schaaper for their critical readings of this article.

Manuscript received May 9, 2002; Accepted for publication August 6, 2002.


*  LITERATURE CITED
*TOP
*ABSTRACT
*MATERIALS AND METHODS
*RESULTS
*DISCUSSION
*LITERATURE CITED

ABDUS SATTAR, K. M. A., T.-C. LIN, C. JONES, and W. H. KONIGSBERG, 1996  Functional consequences and exonuclease kinetic parameters of point mutations in bacteriophage T4 DNA polymerase. Biochemistry 35:16621-16629.[Medline]

BEBENEK, A., H. K. DRESSMAN, G. T. CARVER, S. N. NG, and V. PETROV et al., 2001  Interacting fidelity defects in the replicative DNA polymerase of bacteriophage RB69. J. Biol. Chem. 276:10387-10397.[Abstract/Free Full Text]

BEBENEK, K. and T. A. KUNKEL, 1990  Frameshift errors initiated by nucleotide misincorporation. Proc. Natl. Acad. Sci. USA 87:4946-4950.[Abstract/Free Full Text]

BEBENEK, K. and T. A. KUNKEL, 1995  Analyzing the fidelity of DNA polymerases. Methods Enzymol. 262:217-232.[Medline]

BEBENEK, K. and T. A. KUNKEL, 2000  Streisinger revisited: DNA synthesis errors mediated by substrate misaligments. Cold Spring Harbor Symp. Quant. Biol. 65:81-91.[Medline]

BEBENEK, K., C. M. JOYCE, M. P. FITZGERALD, and T. A. KUNKEL, 1990  The fidelity of DNA synthesis catalyzed by derivatives of Escherichia coli DNA polymerase I. J. Biol. Chem. 265:13878-13887.[Abstract/Free Full Text]

BLATTNER, F. R., G. PLUNKETT, III, C. A. BLOCK, N. T. PERNA, and V. BURLAND et al., 1997  The complete genome sequence of Escherichia coli K-12. Science 277:1453-1474.[Abstract/Free Full Text]

BOOSALIS, M. S., D. W. MOSBAUGH, R. HAMATAKE, A. SUGINO, and T. A. KUNKEL et al., 1989  Kinetic analysis of base substitution mutagenesis by transient misaligment of DNA and by miscoding. J. Biol. Chem. 264:11360-11366.[Abstract/Free Full Text]

BORJAC, J., 1997 RNA determinants of specificity of two phylogenetically related single-stranded DNA binding proteins. Ph.D. Thesis, Tulane University, New Orleans.

BRAITHWAITE, D. K. and J. ITO, 1993  Compilation, alignment, and phylogenetic relationships of DNA polymerases. Nucleic Acids Res. 21:787-802.[Free Full Text]

CLAYTON, L. K., M. F. GOODMAN, E. W. BRANSCOMB, and D. J. GALAS, 1979  Error induction and correction by mutant and wild type T4 DNA polymerase. J. Biol. Chem. 254:1902-1912.[Abstract/Free Full Text]

DAS, S. K. and R. K. FUJIMURA, 1979  Processiveness of DNA polymerases, a comparative study using a simple procedure. J. Biol. Chem. 254:1227-1232.[Abstract/Free Full Text]

DOUBLIÉ, S., S. TABOR, A. M. LONG, C. C. RICHARDSON, and T. ELLENBERGER, 1998  Crystal structure of a bacteriophage T7 DNA replication complex at 2.2 Å resolution. Nature 391:251-258.[Medli