Epidemiologic studies have revealed a complex association between human genetic variance and cancer risk. Quantitative biological modeling based on experimental data can play a critical role in interpreting the effect of genetic variation on biochemical pathways relevant to cancer development and progression. Defects in human DNA base excision repair (BER) proteins can reduce cellular tolerance to oxidative DNA base damage caused by endogenous and exogenous sources, such as exposure to toxins and ionizing radiation. If not repaired, DNA base damage leads to cell dysfunction and mutagenesis, consequently leading to cancer, disease, and aging. Population screens have identified numerous single-nucleotide polymorphism variants in many BER proteins and some have been purified and found to exhibit mild kinetic defects. Epidemiologic studies have led to conflicting conclusions on the association between single-nucleotide polymorphism variants in BER proteins and cancer risk. Using experimental data for cellular concentration and the kinetics of normal and variant BER proteins, we apply a previously developed and tested human BER pathway model to (i) estimate the effect of mild variants on BER of abasic sites and 8-oxoguanine, a prominent oxidative DNA base modification, (ii) identify ranges of variation associated with substantial BER capacity loss, and (iii) reveal nonintuitive consequences of multiple simultaneous variants. Our findings support previous work suggesting that mild BER variants have a minimal effect on pathway capacity whereas more severe defects and simultaneous variation in several BER proteins can lead to inefficient repair and potentially deleterious consequences of cellular damage. (Cancer Epidemiol Biomarkers Prev 2006;15(5):1000–8)

DNA is continuously being damaged by endogenous sources, such as oxidative base modification resulting from reactive oxygen species attack (reviewed in ref. 1). Reactive oxygen species are formed as metabolites during normal cellular respiration, with increased levels of reactive oxygen species, as well as increased DNA base damage, associated with ischemia (2, 3), hyperoxia in human lung cells (4), human artherosclerotic plaques (5), and chronic inflammation in ulcerative colitis patients (6). Reactive oxygen species are also formed as secondary damaging agents as a result of exposure to toxins and ionizing radiation (e.g., sunlight, environmental IR). DNA base damage can lead to genotoxicity, replication stalling, apoptosis, and mutagenesis, with physiologic consequences including disease, aging, and cancer, as reviewed in ref. 7. Modified bases not only affect transcriptional integrity and replication but also DNA-protein binding (8, 9). Consequently, oxidative DNA base damage is normally removed and repaired by the efficient base excision repair (BER) pathway system of prokaryotic and eukaryotic cells (reviewed in refs. 1, 10, 11).

Population screens have identified extensive genetic variation, in particular single-nucleotide polymorphisms, altering the primary amino acid sequence of DNA BER proteins (e.g., ref. 12). This genetic variation has been found in some instances to diminish the stability or efficiency of the encoded repair protein, reducing its kinetic effectiveness (13-15). Furthermore, sequence variants with potential kinetic effect have been predicted computationally (16). Splice variants of repair genes have also been identified in cancerous cells, some of which may lead to a dominant-negative protein form with reduced repair efficiency (17, 18).

Nonfunctional DNA repair alleles, such as BRCA1/2 for breast cancer (19) and MSH2 and MLH1 for hereditary nonpolyposis colorectal cancer (20), are “cancer genes” strongly linked to increased cancer risk. However, cancer genes account for only ∼5% of known cancer cases (21, 22). Consequently, reduced function repair protein variants may be exposure-dependent susceptibility alleles responsible for the majority of observed increased cancer risk associated with family history (23) as well as sporadic cancer incidence. Epidemiologic evidence for a correlation of BER gene variants with cancer risk is abundant, albeit conflicting, as reviewed in refs. 24, 25. Nonetheless, extracts from cells of head and neck cancer patients (26) and smokers with lung cancer have been found to exhibit reduced DNA repair efficiency (27) associated with the presence of BER protein variants and defective BER protein activity, respectively.

The BER pathway system includes proteins with multiple enzymatic activities, cooperativity, and compensatory subpathways. This complexity makes it difficult to intuitively estimate how variation in a single repair protein will change the overall ability of a cell to tolerate DNA base damage. Moreover, there is a high likelihood that individuals will possess multiple variations, which in combination can further increase or reduce pathway capacity. Consequently, we apply a mathematical model of human BER that integrates biological knowledge of enzyme mechanisms and biochemical data on enzyme kinetics and protein concentration obtained from the literature. This model has been used to interpret published data from in vitro pathway reconstitution and cell extracts to evaluate mechanistic hypotheses of enzyme cooperativity and coordination and predict the relative significance of the BER subpathways (28). Notably, the model recently predicted that under normal conditions, background oxidative DNA base damage level is at the low end of the widely varying published measurements (i.e., up to hundreds of lesions per cell) and that the background level is relatively stable to small changes in enzyme kinetics (29). In the context of the normally robust operation of BER, we applied the model based on published experimental data to (i) estimate the potential effect if any of the previously characterized BER variant proteins, (ii) identify the magnitude of changes in individual protein kinetics and concentration required to observe a substantial difference in overall BER kinetics and capacity, and (iii) use the model to predict the effect of simultaneous variation in multiple repair proteins on total BER kinetics.

BER Pathway Model Structure

The comprehensive BER pathway model is based on current biological knowledge, which is illustrated in Fig. 1. The BER pathway in the nucleus is initiated by a lesion-specific DNA glycosylase, which recognizes and removes damaged bases (30). For example, 8-oxoguanine (8-oxoG; y1 in Fig. 1) is removed by 8-oxoG DNA glycosylase (Ogg1; e1). Some DNA glycosylases (including Ogg1) also exhibit an apurinic/apyrimidinic (AP) lyase activity, incising 3′ to the abasic (AP) residue, leading to a subpathway (pathway A) in which the phosphodiesterase activity of human AP endonuclease (Ape1; e2) removes the 3′ block (y10; refs. 31, 32). In most cases, “short-patch” BER (pathway B) occurs (33), where Ape1 incises 5′ to the lesion (34). The resultant 5′-deoxyribose phosphate (y4) group is removed by the deoxyribose phosphate lyase activity of DNA polymerase β (Polβ; e3), which then executes single-base replacement at the gapped site (y5 or y11). In our model, the Polβ-catalyzed activities may occur in either order (thus the fork in pathway B). An alternative “long-patch” BER pathway (pathway C) exists in which gap-filling of two to seven nucleotides is done by Polδ (e4) or possibly Polε (excluded from our model due to a lack of kinetic data), in complex with proliferating cell nuclear antigen, followed by flap excision by Fen1 (e5) and subsequent ligation (35). The final step in all subpathways is done by DNA ligase 1 (Lig1; e6) or Lig3 (e7) stabilized in complex with Xrcc1 (36). Lig1 and Lig3/Xrcc1 activities seem to be largely interchangeable (37), although a few recent results indicate that Lig1 may be more significant for long-patch repair (38).

Figure 1.

Schematic of the human BER pathway model, described in the main text, showing substrates and intermediates (y), enzymes (e), and reactions (v). (*, abasic sites can be formed directly, y30, or as intermediates of a bifunctional DNA glycosylase, e.g. Ogg1, y2.).

Figure 1.

Schematic of the human BER pathway model, described in the main text, showing substrates and intermediates (y), enzymes (e), and reactions (v). (*, abasic sites can be formed directly, y30, or as intermediates of a bifunctional DNA glycosylase, e.g. Ogg1, y2.).

Close modal

Model Assumptions

In general, reactions in the model are modeled using Michaelis-Menten enzyme kinetics, with kcat and KM variables defining a reaction velocity (vi in Fig. 1). This implicitly assumes homogeneity and deterministic, continuous reaction kinetics, which are justified by the high concentration of repair proteins as discussed in ref. 28. If Michaelis-Menten assumptions hold for a biochemical system, then kcat is interpreted as the catalytic turnover rate of the enzyme (its “activity”) and KM is a function of the affinity of the enzyme for the substrate.

To account for experimental measurement of 5- to 10-fold reduced diffusion velocities in the crowded cell nucleus (39), KM of all reactions are increased 10-fold (due to the reduction in collision probability). The model includes experimentally measured protein cooperativity by kcat increases as in previous modeling (28, 29), as well as including coordination of consecutive single enzyme activities (such as Polβ and the dual-function DNA glycosylases, e.g., Ogg1) by modeling the second step as a first-order reaction, as described in ref. 28. We also assume no product inhibition for individuals under normal conditions in the pathway, as suggested by observations (32, 40). We do not, however, include the potential role of Xrcc1 as a mediator of a multiprotein coordination complex (41), as the contribution of Xrcc1 to such a complex forming and any consequent enhancement of BER kinetics are unclear (42, 43). In addition, many of the potential coordinative effects are already in the model. Table 1 shows kinetic variables used in the model, including rate enhancements due to cooperativity.

Table 1.

Normal BER protein kinetics

Reactions (vi)kcat (s−1)*KM (nmol/L)*ReferencesRate enhancementsReferences
0.0052 121.5 (14, 15, 32, 109-111) 10× decreased KM (Ape1); 8× increased kcat (Ape1; posttranslational modifications) (32, 112) 
3.2028 34.7 (32, 113-120) — — 
3, 6 0.817 210 (121-124) — — 
4, 5 0.075 500 (125) 6× increased kcat (Ape1) (126, 127) 
7, 10, 14, 16-18 0.0213 56.7 (121, 128, 129) 4× increased kcat (PCNA, RpA) (129, 130) 
8, 15 0.6 100 (131-133) Rate in presence of PCNA — 
0.1335 39 (124, 134) 2× increased kcat (Ape1) (135) 
11 0.000887 7.2 (32) 2.5× increased kcat (Ape1; posttranslational modification) (112, 136) 
12 0.0608 169 (114, 137, 138) — — 
Reactions (vi)kcat (s−1)*KM (nmol/L)*ReferencesRate enhancementsReferences
0.0052 121.5 (14, 15, 32, 109-111) 10× decreased KM (Ape1); 8× increased kcat (Ape1; posttranslational modifications) (32, 112) 
3.2028 34.7 (32, 113-120) — — 
3, 6 0.817 210 (121-124) — — 
4, 5 0.075 500 (125) 6× increased kcat (Ape1) (126, 127) 
7, 10, 14, 16-18 0.0213 56.7 (121, 128, 129) 4× increased kcat (PCNA, RpA) (129, 130) 
8, 15 0.6 100 (131-133) Rate in presence of PCNA — 
0.1335 39 (124, 134) 2× increased kcat (Ape1) (135) 
11 0.000887 7.2 (32) 2.5× increased kcat (Ape1; posttranslational modification) (112, 136) 
12 0.0608 169 (114, 137, 138) — — 

Abbreviations: PCNA, proliferating cell nuclear antigen; RpA, replication protein A.

*

Where multiple sources are used in the literature, the median is used to avoid excessive influence by outliers.

Assuming pathway coordination, this reaction is first order with a rate of kcat.

Lig3 kinetics assumed equal to those for Lig1 kinetics (37).

Model Protein Concentrations

Table 2 shows BER protein concentrations determined from published data calculated by assuming the number of proteins measured per cell is homogeneously distributed in a nucleus with a diameter of 5 μm, as described in refs. 28, 29. Because of the limited data currently available, these protein levels were drawn from various cell lines. Where multiple measurements were available, a conservative assumption was used because some are from tumor cell lines, which have been measured to have higher protein expression levels than normal cells (44). Furthermore, Ape1 (45, 46) and DNA glycosylases, such as Ogg1 (47), may be distributed in both the nucleus and other cell compartments (i.e., cytosol and mitochondria). In our analysis, we neglect the contribution to overall cellular DNA damage of the oxidative damage in mitochondria, which is repaired equally or more efficiently than nuclear DNA and does not occur at levels measurable in whole-cell assays (48-50). Consequently, protein concentrations in Table 2 do not include Ogg1 and Ape1 localized to mitochondria (47, 51, 52). We have found in our analysis of the sensitivity of our results to protein concentration that it closely follows that of sensitivity to kcat, which is published here (results not shown; similar results were also found in our previous work on estimating steady-state BER capacity; ref. 29). This is expected from the mathematical form of the Michaelis-Menten equations used to model most of the BER pathway.

Table 2.

BER protein concentrations

ProteinConcentration (nmol/L)References
Ogg1 406 (51, 111) 
Ape1 2,000 (51, 137, 139) 
Polβ 419 (44, 140) 
Polδ 600 
Fen1 450 
Lig1 254 (141) 
Lig3 254  
ProteinConcentration (nmol/L)References
Ogg1 406 (51, 111) 
Ape1 2,000 (51, 137, 139) 
Polβ 419 (44, 140) 
Polδ 600 
Fen1 450 
Lig1 254 (141) 
Lig3 254  
*

In absence of data, assumed based on concentration of other nuclear proteins.

Assumed equal to Lig1 (37).

Solution of Model Equations

The mathematical structure of the model is identical to that in ref. 29. Differential equations were solved using the “ode15s” stiff solver of MATLAB R14 using variables from Tables 1 and 2, adjusted by the analysis described below. Steady-state reactant levels were the equilibrium point found by simulating the differential equations on a time scale of 24 hours and then setting equations for repair intermediate concentrations equal to zero and using the nonlinear least squares MATLAB routine “fsolve” to solve the resulting algebraic system (results not shown).

Simulation of 8-OxoG Levels as an Estimate for Oxidative DNA Base Repair Capacity

We focus on the repair of 8-oxoG initiated by Ogg1, which includes the full set of pathway steps in Fig. 1. This is motivated by the large amount of data from experimental sources for the kinetics of wild-type and polymorphic forms of Ogg1 (28), as well as cellular and physiologic levels of 8-oxoG (e.g., ref. 53); also reviewed in ref. 29. 8-OxoG is frequently used as a marker for oxidative DNA damage associated with disease consequences (54). As further data accumulate for the kinetics and cellular abundance of other DNA glycosylases (55), we can extend our analysis in the future using the approach described herein. In our studies, we assume that the contribution of NEIL1/2 enzymes recently found to have 8-oxoG DNA glycosylase activity (56) is quantitatively insignificant as compared with Ogg1; recent studies have suggested that their biological role may be primarily for specialized cases, such as damage within bubble DNA structures (57, 58).

We account for the production of other base lesions in our modeling by including an additional damage load of abasic (AP) sites, comprising the products of other DNA glycosylases, as well as spontaneous hydrolysis and base loss, the majority of which result in apurinic sites (59). Overall, we predict the combined numbers of original lesions and repair intermediates, which are a significant component of observed lesions (60) and may be highly genotoxic and mutagenic themselves (61, 62).

Estimated Effect of Mild Protein Variants Found in the Population on BER Kinetics and Capacity

If a protein variant is found to affect the kinetics of a BER reaction, then its effect on the overall pathway is estimated by modifying the appropriate kcat and KM variables for that reaction in the comprehensive BER system model. Published data are available for the kinetics of some BER proteins coded by sequence variants identified in population screens. Ape1 variants were characterized with AP endonuclease efficiencies from 35% to 110% of wild type (13). These are calculated based on the amount of repair of an initial lesion concentration over a fixed period of time, and they may be interpreted as being equivalent to the catalytic rate (kcat) of Ape1. The commonly found Ogg1 polymorphism Ser326Cys has a reported kcat 63% of wild type (15) and other Ogg1 variants were found to have a kcat 44% (Arg154His) and 85% (Arg36Gln) of wild type (14). In the case of Polβ, one quantitatively characterized allelic variant is a dominant-negative inhibitor of wild-type activity (17). Based on the measured equilibrium binding affinity (KD) of the variant with the Polβ substrate and potential cell concentration, the effective KM for Polβ gap-filling reaction is increased 3- to 5-fold (calculated based on competitive inhibition in the Michaelis-Menten model). The estimated effect of protein variants in these ranges of variation are shown in Fig. 2, where we have simulated the effect of Ape1 with 10% of normal kcat (the highest amount at which any difference could be observed), Ogg1 with 50% of normal kcat, and Polβ with a 5-fold increase in KM.

Figure 2.

Sensitivity of 8-oxoG BER capacity and kinetics to mild kinetic defects in Ogg1 (50% of normal kcat), Ape1 (10% of kcat), and Polβ (5-fold increased KM). A. Total number of 8-oxoG, abasic, and repair intermediate bases per cell given increasing 8-oxoG rates of formation (assuming simultaneous production of an additional 20,000 AP sites/cell/d). B. Time required to clear an instantaneous pulse of 100,000 AP sites/cell and an increasing number of 8-oxoG lesions/cell and restore damage to the background steady-state level for continuous formation of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d. Simulated data for the pathway with normal kinetics and with a 10% kcat Ape1 variant are nearly overlapping.

Figure 2.

Sensitivity of 8-oxoG BER capacity and kinetics to mild kinetic defects in Ogg1 (50% of normal kcat), Ape1 (10% of kcat), and Polβ (5-fold increased KM). A. Total number of 8-oxoG, abasic, and repair intermediate bases per cell given increasing 8-oxoG rates of formation (assuming simultaneous production of an additional 20,000 AP sites/cell/d). B. Time required to clear an instantaneous pulse of 100,000 AP sites/cell and an increasing number of 8-oxoG lesions/cell and restore damage to the background steady-state level for continuous formation of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d. Simulated data for the pathway with normal kinetics and with a 10% kcat Ape1 variant are nearly overlapping.

Close modal

Figure 2 illustrates two possible kinds of effect on the BER pathway: the capacity of BER to tolerate a constant rate of base lesion formation (8-oxoG and AP site formation combined) and the kinetics of BER repairing an instantaneous pulse of oxidative DNA damage (also consisting of both 8-oxoG and AP site lesions). Such a pulse simulates acute exposure (e.g., a dose of radiation therapy, increased metabolism during intense exercise, or acute inflammation during an infection). We consider repair of transient damage pulses in the context of continuing background base lesion formation. Base lesion formation rates (2,000 8-oxoG/cell/day and 20,000 AP sites/cell/day) are based on experimental data for normal background rates of formation (59, 63, 64). They are adjusted upwards to reflect (i) increased 8-oxoG corresponding to exercise and other variations in metabolism (65) and (ii) AP sites produced by the action of additional DNA glycosylases processing other base lesions. Damage “clearance” is defined as the time at which the background steady-state level of lesions corresponding to the constant damage formation rate is restored. Notably, as Fig. 2 shows, 8-oxoG repair rate is more sensitive to small changes in protein kinetics than the steady-state repair capacity.

Table 3 summarizes the effect on BER capacity (steady-state number of lesions per cell on exposure to the formation of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d) and BER kinetics (time to clear 10,000 8-oxoG and 100,000 AP sites) of each protein in the pathway with 10% and 50% of their normal (Table 1) kcat. Results are given for each of the activities of multiple-function proteins separately. Data are provided as a percentage increase in steady-state damage and clearance time with respect to normal protein kinetics. Table entries of “0%” represent a change <0.01% (results not shown). Notably, because of the dominance of pathway B for our model of 8-oxoG repair, assuming normal Polβ function, pathway C proteins, Polδ and Fen1, will have negligible effect on kinetics, as was found in our previous modeling analysis (28).

Table 3.

Predicted effect of protein kinetic variants on background damage levels and time to clear lesions (% of wild type)

ProteinActivity (vi)10% kcat
50% kcat
DamageTimeDamageTime
Ogg1 Base excision 29% 341% 3% 9% 
Ape1 5′-incision 2% 1% 0% 0% 
Polβ 3, 6 Gap-filling 13% 7% 4% 2% 
Polβ 4, 5 5′-Deoxyribose phosphate lyase 32% 56% 4% 1% 
Lig1* 7, 10, 14, 16-18  53% 90% 21% 34% 
Polδ 8, 15  1% 0% 0% 0% 
Fen1  0% 0% 0% 0% 
Ogg1 11 AP lyase 0% 0% 0% 0% 
Ape1 12 3′-Deoxyribose phosphate lyase 0% 0% 0% 0% 
ProteinActivity (vi)10% kcat
50% kcat
DamageTimeDamageTime
Ogg1 Base excision 29% 341% 3% 9% 
Ape1 5′-incision 2% 1% 0% 0% 
Polβ 3, 6 Gap-filling 13% 7% 4% 2% 
Polβ 4, 5 5′-Deoxyribose phosphate lyase 32% 56% 4% 1% 
Lig1* 7, 10, 14, 16-18  53% 90% 21% 34% 
Polδ 8, 15  1% 0% 0% 0% 
Fen1  0% 0% 0% 0% 
Ogg1 11 AP lyase 0% 0% 0% 0% 
Ape1 12 3′-Deoxyribose phosphate lyase 0% 0% 0% 0% 

NOTE: Steady-state number of lesions per cell is based on 2,000 8-oxoG and 20,000 AP sites/cell/d; time to clear is calculated for 10,000 8-oxoG and 100,000 AP sites/cell.

*

Based on model assumptions, Lig3 will have the same results.

Sensitivity Analysis of BER Pathway Capacity to Hypothetical Protein Variants

In Fig. 2, curves showing an increased lesion burden and clearance time as a result of increased damaged load show a quantitatively small difference for mild kinetic variants of BER proteins. This suggests that the BER pathway is apparently robust to mild kinetic variation. Accordingly, we explored a wider range of protein kinetic variables to identify levels at which variants may have a qualitatively observable effect on BER pathway capacity, as shown in Fig. 3 for kcat and KM (modified relative to the normal values in Table 1; not shown are variants for BER reactions that were found to have no effect on steady-state damaged base levels, e.g., Fen1 kcat and KM). Whereas the analysis in Fig. 3 is specifically for steady-state damage levels given physiologically relevant levels of 8-oxoG and AP site formation, it qualitatively reflects the general sensitivity of the BER pathway to its individual components. Regions on the plots in Fig. 3 that show a rapidly changing curve represent variable values for which BER capacity is highly sensitive and less robust to increases in damage level. Regions for which no curve is plotted at all represent variable ranges of BER pathway kinetics that preclude a steady-state equilibrium point. For these hypothetical kinetic variants, the BER pathway breaks down, with lesions and repair intermediates accumulating uncontrollably. For normal BER reaction kinetics, this breaking point is ∼125,000 8-oxoG/cell/d (29), a level that seems to be physiologically unrealistic. By comparison, in cell culture, 8-oxoG formation was measured at 0.01/106 bases/Gy (66), and 8-oxoG levels on the order of 100/106 have been measured in naked DNA exposed to highly oxidative agents or >5 Gy radiation (67). However, for some modeled changes in the kcat and KM of Ogg1, Ape1, Polβ, and Lig1, it can be on the order of 104 8-oxoG/cell/d (data not shown) or less (e.g., where the curves are not shown in Fig. 3). Such variants will result in a BER pathway incapable of tolerating even normal damage formation rates (i.e., 103 8-oxoG/cell/d and 104 AP sites/cell/d) without accumulating unrepaired lesions.

Figure 3.

Sensitivity over a broad range of kcat (A) and KM (B) of enzymatic activities in BER (Ogg1 excision of 8-oxoG, Ape1 nicking at the 5′ end, Polβ gap-filling and 5′-deoxyribose phosphate (dRp) lyase reactions, Polδ gap-filling, and the ligase activity of Lig1. Changes in kinetic variables are fractions and multiples of the values for those reactions in Table 1 (with all KM increased 10-fold as described in the text). The steady-state total number of 8-oxoG, AP sites, and repair intermediates is calculated based on a constant total damage formation rate of 2,000 8-oxoG lesions/cell/d and 20,000 AP sites/cell/d. The plot is truncated at the point where a continued KM increase results in no physically possible steady state.

Figure 3.

Sensitivity over a broad range of kcat (A) and KM (B) of enzymatic activities in BER (Ogg1 excision of 8-oxoG, Ape1 nicking at the 5′ end, Polβ gap-filling and 5′-deoxyribose phosphate (dRp) lyase reactions, Polδ gap-filling, and the ligase activity of Lig1. Changes in kinetic variables are fractions and multiples of the values for those reactions in Table 1 (with all KM increased 10-fold as described in the text). The steady-state total number of 8-oxoG, AP sites, and repair intermediates is calculated based on a constant total damage formation rate of 2,000 8-oxoG lesions/cell/d and 20,000 AP sites/cell/d. The plot is truncated at the point where a continued KM increase results in no physically possible steady state.

Close modal

Figure 3 shows the sensitivity of repair capacity to changes in only Lig1 kinetic efficiency. The concentration and kinetics of Lig3 were assumed to be equivalent to those of Lig1. Consequently, simulated predictions for Lig3 are identical to those for Lig1, showing redundancy in the DNA ligation step of the repair pathway. This agrees partially with experimental results in ref. 37, which do suggest that in some cases one DNA ligase may be preferred to another. Due to the redundant role of Lig1 and Lig3 in our model, simulated data predict that BER capacity is sensitive to reduction in Lig1 kinetics to a certain point; after which, further reduction no longer affects the pathway, as shown in Fig. 3.

Observation of Unexpected Synergies for Hypothetical Multiple Protein Variants in the BER Pathway

In some cases, the complexity of individual protein kinetics within the context of the whole BER pathway leads to unexpected synergistic variation. To illustrate this, we estimate the consequence of simultaneous variation of Ape1 5′-incision activity and Polβ gap-filling and 5′-deoxyribose phosphate lyase catalytic activities (kcat) on overall BER dynamics for the repair of a transient pulse of DNA damage, assuming that the Polβ rate enhancement by Ape1 is unaffected. In the steady-state analysis (not shown), there is only a 3.75 to 4 lesions/cell difference (including repair intermediates) between Polβ variants with 200-fold reduced kcat and 10-fold reduced kcat (assuming constant rates of formation of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d), with the exception of the situation with a 1,000-fold Ape1 kcat reduction, in which case the number of steady-state lesions differs by 27 per cell. This contrasts with the findings shown in Fig. 4 for transient repair dynamics. The time to fully clear a pulse of 8-oxoG and abasic lesions and repair intermediates for the Polβ-inefficient pathway is at its maximum when Ape1 kcat is 70% of normal. As Fig. 4B shows, this point of respective Ape1 and Polβ activity modification is a peak (darkest region). Repair efficiency improves (more lightly shaded regions of Fig. 4B) as Ape1 kinetics either increase or decrease, and the overall system shows less sensitivity to reduced Polβ kinetics.

Figure 4.

Simulation of the effect of varying Ape1 and Polβ kcat simultaneously on the time required to clear 10,000 8-oxoG sites and 100,000 AP sites, with a constant background damage formation rate of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d. A. Clearance time for varying Ape1 kcat (multiple of normal variable value in Table 1) for fractions of normal Polβ kcat. B. Contour plot of clearance time for simultaneous changes (log ratio to normal kcat) in Ape1 and Polβ; darker regions correspond to higher clearance times.

Figure 4.

Simulation of the effect of varying Ape1 and Polβ kcat simultaneously on the time required to clear 10,000 8-oxoG sites and 100,000 AP sites, with a constant background damage formation rate of 2,000 8-oxoG/cell/d and 20,000 AP sites/cell/d. A. Clearance time for varying Ape1 kcat (multiple of normal variable value in Table 1) for fractions of normal Polβ kcat. B. Contour plot of clearance time for simultaneous changes (log ratio to normal kcat) in Ape1 and Polβ; darker regions correspond to higher clearance times.

Close modal

Based on a study of the estimated levels of BER intermediates during the simulation, it seems that at that particular critical level of Ape1, there is a shifting of repair to pathways A and C—which are slower than pathway B—and sites being repaired through pathway A are particularly sensitive to Polβ kinetics following the first steps. Consequently, there is a persistence of repair intermediates and prolongation of the time to reach the normal steady-state background. As this model analysis shows, the complexity of multiple subpathways in the BER pathway can be a mechanism for unexpected synergies between protein variants, which are not merely additive results of reduced kinetic rates.

There have been numerous investigations of potential correlations between BER gene variants (e.g., single-nucleotide polymorphisms, truncations, and frameshifts) with levels of biomarkers, including DNA damage and mitotic delay, prevalence in tumor cells, and increased cancer susceptibility. To date, many of these studies have been inconclusive or contradictory. For instance, the commonly found Ser326Cys variant of Ogg1 has been found to be correlated epidemiologically with increased cancer risk, increased measurements of cellular DNA base damage, and reduced base repair kinetics in several studies and assays (68-75), whereas null or statistically insignificant correlations have been found in others (refs. 69, 76-79); also a recent extensive review in ref. 80). Applying a previously developed and tested comprehensive BER pathway model integrating reaction mechanisms and variables obtained from published experimental data, we estimate that Ogg1 base excision kcat reduced by 50%, comparable to that of the Ogg1 Ser326Cys polymorphism (15), results in a decrease of ∼5% to 10% in BER capacity, as defined by the steady-state number of lesions per cell for a wide range of constant damage formation rates (as shown in Fig. 2). The transient kinetics of repairing a pulse of damage are somewhat more sensitive to mild kinetic variants but the effects are still minimal.

Overall, our modeling supports a complex role for genetics in tolerating exposure to genotoxic stresses in determining DNA base damage tolerance and cancer risk. Under normal and increased damage loads and rates (simulated based on experimental measurements), results favor the null or marginal experimental findings for the effect of Ogg1 Ser326Cys on cellular tolerance of oxidative damage and cancer susceptibility. Furthermore, Ape1 population variants characterized in ref. 13 are estimated by our modeling to have a negligible effect on BER kinetics and capacity. These results are supported by a null finding for cancer susceptibility with the Ape1 variant Asp146Glu in pooled studies (75).

Our analysis (Fig. 3) reveals that further increases in individual repair protein kinetics have a minor effect on overall BER pathway kinetics, as suggested by the long evolutionary history of BER, conserved to a relatively high degree from E. coli through human cells (7). Of all its steps, the BER pathway is most sensitive to decreases in Ogg1 base excision, Ape1 5′-incision of the abasic site, and Polβ gap-filling and 5′-deoxyribose phosphate lyase reaction kinetics, with limited compensated sensitivity to DNA ligase activity. The high sensitivity to Ogg1 kinetics is supported by evidence that allelic loss (which reduces Ogg1 concentration and reaction velocity) is prevalent in tumor cells with reduced 8-oxoG repair activity (69, 76, 81, 82). Furthermore, whereas there is little evidence for Polβ variants in cells of normal tissues in the human population, several variants with diminished activity have been found in tumor cell lines and are associated with reduced overall BER pathway capacity and a Polβ mutator phenotype, in accord with estimates presented here (17, 18, 83-88). This suggests a potentially significant role for inherited Polβ variants in increasing cellular susceptibility to DNA base damage and cancer risk, suggesting targets for further epidemiologic studies.

Whereas the assumption of redundant Lig1 and Lig3 activities results in a modeling prediction of some compensation for the case of variants, there is still notable sensitivity to kinetics that would be exacerbated if activities are not completely overlapping. A Lig1 Arg771Trp variant has been characterized in an immunodeficient patient, with cell cultures showing hypersensitivity to DNA damaging agents and impaired BER activity with 3% to 5% Lig1 kinetics (89), and a noncoding variant of Lig3 has been associated with increased esophageal cancer risk (90). In addition, Lig3 activity requires Xrcc1, and there are reports associating Xrcc1 variants with increased cancer susceptibility (70, 90-94). However, other studies indicate null or statistically marginal results for Xrcc1 variants, highlighting the complexity of variant protein roles in the BER pathway (95-97).

Whereas the association of most proteins involved in oxidative DNA base damage repair with cancer remains a controversial question, there is recent evidence of a direct link between protein sequence of MYH (the DNA glycosylase that removes a mismatched adenine across from the 8-oxoG lesion) and colorectal cancer (refs. 98-100; reviewed in ref. 101). There is also some evidence that MYH variants associated with increased cancer risk and polyposis have reduced catalytic efficiency but published studies have been on human MYH sequence variants generated by site-directed mutagenesis not found in the population (102) or mutants of the E. coli homologue corresponding to the population variants (103). As further enzymatic data are obtained for MYH activity, our comprehensive BER system model can be extended to include MYH along with Ogg1 and other DNA glycosylases, allowing for estimation of variant effect on the whole pathway.

In the post-genome era, there has been tremendous investment in the search for genetic variants in the population related to increased cancer risk and chemoresistance, particularly single-nucleotide polymorphisms (104). The relative prevalence of single-nucleotide polymorphisms in the population suggests that many are low-penetrance alleles, which will reduce the efficiency of biological pathways without eliminating their function entirely. As this study shows, biochemically modeling the effect of variants reveals complex behavior and pathway sensitivity and can provide insight to interpret incomplete and contradictory results obtained from both epidemiology and laboratory experiments. Progress in understanding how complex genotypes can predict the response to endogenous background and exogenous environmental exposure to DNA damaging stresses will allow an individualized, predictive approach in the interpretation of the medical significance of DNA damage quantification (105), identification of chemoprevention strategies for individuals with known susceptibility factors (106, 107), and unveiling of mechanisms of chemoresistance (108) to therapeutics. In general, the analysis presented here underscores the critical role quantitative biological system modeling and experiments will play in the long-term goal of the post-genomic era, effectively interpreting complex genetic studies and translating their findings to benefit human health.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1
Evans MD, Dizdaroglu M, Cooke MS. Oxidative DNA damage and disease: induction repair and significance.
Mutat Res
2004
;
567
:
1
–61.
2
Lin LH, Cao S, Yu L, Cui J, Hamilton WJ, Liu PK. Up-regulation of base excision repair activity for 8-ydroxy-2′-deoxyguanosine in the mouse brain after forebrain ishemia-reperfusion.
J Neurochem
2000
;
74
:
1098
–105.
3
Lan J, Li W, Zhang F, et al. Inducible repair of oxidative DNA lesions in the rat brain after transient focal ischemia and reperfusion.
J Cereb Blood Flow Metab
2003
;
23
:
1324
–39.
4
Wu M, He YH, Kobune M, Xu Y, Kelley MR, Martin WJ II. Protection of human lung cells against hyperoxia using the DNA base excision repair genes hOgg1 and Fpg.
Am J Respir Crit Care Med
2002
;
166
:
192
–9.
5
Martinet W, Knaapen MWM, De Meyer GRY, Herman AG, Kocks MM. Elevated levels of oxidative DNA damage and DNA repair enzymes in human artherosclerotic plaques.
Circulation
2002
;
106
:
927
–32.
6
Hofseth LJ, Khan MA, Ambrose M, et al. The adaptive imbalance in base excision-repair enzymes generates microsatellite instability in chronic inflammation.
J Clin Invest
2003
;
112
:
1887
–94.
7
Coates PJ, Lorimore SA, Wright EG. Cell and tissue responses to genotoxic stress.
J Pathol
2005
;
205
:
221
–35.
8
Parsian AJ, Funk MC, Tao TY, Hunt CR. The effect of DNA damage on the formation of protein/DNA complexes.
Mutat Res
2002
;
501
:
105
–13.
9
Opresko PL, Fan J, Danzy S, Wilson DM III, Bohr VA. Oxidative damage in telomeric DNA disrupts recognition by TRF1 and TRF2.
Nucleic Acids Res
2005
;
33
:
1230
–9.
10
Fortini P, Pascucci B, Parlanti E, D'Errico M, Simonelli V, Dogliotti E. The base excision repair: mechanisms and its relevance for cancer susceptibility.
Biochimie
2003
;
85
:
1053
–71.
11
Barnes DE, Lindahl T. Repair and genetic consequences of endogenous DNA base damage in mammalian cells.
Annu Rev Genet
2004
;
38
:
445
–76.
12
Mohrenweiser HW, Xi T, Vazquez-Matias J, Jones IM. Identification of 127 amino acid substitution variants in screening 37 dna repair genes in humans.
Cancer Epidemiol Biomarkers Prev
2002
;
11
:
1054
–64.
13
Hadi MZ, Coleman MA, Fidelis K, Mohrenweiser HW, Wilson DM III. Functional characterization of Ape1 variants identified in the human population.
Nucleic Acids Res
2000
;
28
:
3871
–9.
14
Dherin C, Radicella JP, Dizdaroglu M, Boiteux S. Excision of oxidatively damaged DNA bases by the human α-hOgg1 protein and the polymorphic α-hOgg1 (Ser326Cys) protein which is frequently found in human populations.
Nucleic Acids Res
1999
;
27
:
4001
–7.
15
Audebert M, Radicella JP, Dizdaroglu M. Effect of single mutations in the OGG1 gene found in human tumors on the substrate specificity of the Ogg1 protein.
Nucleic Acids Res
2000
;
28
:
2672
–8.
16
Xi T, Jones IM, Mohrenweiser HW. Many amino acid substitution variants identified in DNA repair genes during human population screenings are predicted to impact protein function.
Genomics
2004
;
83
:
970
–9.
17
Bhattacharyya N, Banerjee S. A novel role of XRCC1 in the functions of a DNA polymerase β variant.
Biochemistry
2001
;
40
:
9005
–13.
18
Thompson TE, Rogan PK, Risinger JI, Taylor JA. Splice variants but not mutations of DNA polymerase β are common in bladder cancer.
Cancer Res
2002
;
62
:
3251
–6.
19
Ormiston W. Hereditary breast cancer.
Eur J Cancer Care (Engl)
1996
;
5
:
13
–20.
20
Peltomaki P, Vasen HF. Mutations predisposing to heriditary nonpolyposis colorectal cancer: database and results of a collaborative study. The International Collaborative Group on Hereditary Nonpolyposis Colorectal Cancer.
Gastroenterology
1997
;
113
:
1146
–58.
21
Madigan MP, Ziegler RG, Benichou J, Byrne C, Hoover RN. Proportion of breast cancer cases in the United States explained by well-established risk factors.
J Natl Cancer Inst
1995
;
87
:
1681
–5.
22
Newman B, Mu H, Butler LM, Millikan RC, Moorman PG, King MC. Frequency of breast cancer attributable to BRCA1 in a population-based series of American women.
JAMA
1998
;
279
:
915
–21.
23
Mohrenweiser HW, Jones IM. Variation in DNA repair is a factor in cancer susceptibility: a paradigm for the promises and perils of individual and population risk estimation?
Mutat Res
1998
;
400
:
15
–24.
24
Goode EL, Ulrich CM, Potter JD. Polymorphisms in DNA repair genes and association with cancer risk.
Cancer Epidemiol Biomarkers Prev
2002
;
11
:
1513
–30.
25
Mohrenweiser HW, Wilson DM III, Jones IM. Challenges and complexities in estimating both the impact and the disease risk associated with the extensive genetic variation in human DNA repair genes.
Mutat Res
2003
;
526
:
93
–125.
26
Sturgis EM, Castillo EJ, Li L, et al. Polymorphisms of DNA repair gene XRCC1 in squamous cell carcinoma of the head and neck.
Carcinogenesis
1999
;
20
:
2125
–9.
27
Paz-Elizur T, Krupsky M, Blumenstein S, Elinger D, Schectman E, Livneh Z. DNA repair activity for oxidative damage and risk of lung cancer.
J Natl Cancer Inst
2003
;
95
:
1312
–9.
28
Sokhansanj BA, Rodrigue GR, Fitch JP, Wilson DM III. A quantitative model of human DNA base excision repair. I. Mechanistic insights.
Nucleic Acids Res
2002
;
30
:
1817
–25.
29
Sokhansanj BA, Wilson DM III. Oxidative DNA damage background estimated by a system model of base excision repair.
Free Radic Biol Med
2004
;
37
:
422
–7.
30
Scharer OD, Jiricny J. Recent progress in the biology, chemistry and structural biology of DNA glycosylases.
Bioessays
2001
;
23
:
270
–81.
31
Pascucci B, Maga G, Hubscher U, et al. Reconstitution of the base excision repair pathway for 7,8-dihydro-8-oxoguanine with purified human proteins.
Nucleic Acids Res
2002
;
30
:
2124
–30.
32
Hill JW, Hazra TK, Izumi T, Mitra S. Stimulation of human 8-oxoguanine-DNA glycosylase by AP-endonuclease: potential coordination of the initial steps in base excision repair.
Nucleic Acids Res
2001
;
29
:
430
–8.
33
Dianov G, Bischoff C, Piotrowski J, Bohr VA. Repair pathways for processing of 8-oxoguanine in DNA by mammalian cell extracts.
J Biol Chem
1998
;
273
:
33811
–6.
34
Wilson DM III, Barsky D. The major human abasic endonuclease Ape1: formation, consequences, and repair of abasic lesions in DNA.
Mutat Res
2001
;
484
:
283
–307.
35
Matsumoto Y, Kim K, Hurwitz J, et al. Reconstitution of proliferating cell nuclear antigen-independent repair of apurinic/apyrimidinic sites with purified human proteins.
J Biol Chem
1999
;
274
:
33703
–8.
36
Tomkinson AE, Mackey ZB. Structure and function of mammalian DNA ligases.
Mutat Res
1998
;
407
:
1
–9.
37
Ho ELY, Satho MS. Repair of single-strand DNA interruptions by redundant pathways and its implication in cellular sensitivity to DNA-damaging agents.
Nucleic Acids Res
2003
;
31
:
7032
–40.
38
Sleeth KM, Robson RL, Dianov GL. Exchangeability of mammalian DNA ligases between base excision repair pathways.
Biochemistry
2004
;
43
:
12924
–30.
39
Phair RD, Misteli T. High mobility of proteins in the mammalian cell nucleus.
Nature
2000
;
404
:
604
–9.
40
Wilson SH, Kunkel TA. Passing the baton in base excision repair.
Nat Struct Biol
2000
;
7
:
176
–8.
41
Marsin S, Vidal AE, Sossou M, et al. Role of XRCC1 in the coordination and stimulation of oxidative DNA damage repair initiated by the DNA glycosylase hOGG1.
J Biol Chem
2003
;
278
:
44068
–74.
42
Parsons JL, Dianova II, Allinson SL, Dianov GL. DNA polymerase β promotes recruitment of DNA ligase IIIα-XRCC1 to sites of base excision repair.
Biochemistry
2005
;
44
:
10613
–9.
43
Wong HK, Kim D, Hogue BA, McNeill DR, Wilson DM III. DNA damage levels and biochemical repair capacities associated with XRCC1 deficiency.
Biochemistry
2005
;
44
:
14335
–43.
44
Srivastava DK, Husain I, Arteaga CL, Wilson SH. DNA polymerase β expression differences in selected human tumors and cell lines.
Carcinogenesis
1999
;
20
:
1049
–54.
45
Evans AR, Limp-Foster M, Kelley MR. Going APE over ref-1.
Mutat Res
2000
;
461
:
83
–108.
46
Pines A, Perrone L, Bivi N, et al. Activation of APE1/Ref-1 is dependent on reactive oxygen species generated after purinergic receptor stimulation by ATP.
Nucleic Acids Res
2005
;
33
:
4379
–94.
47
Nakabeppu Y. Regulation of intracellular localization of human MTH1, OGG1, and MYH proteins for repair of oxidative DNA damage.
Prog Nucleic Acid Res Mol Biol
2001
;
68
:
75
–94.
48
Hudson EK, Hogue B, de Souza-Pinto NC, et al. Age-associated change in mitochondrial DNA damage.
Free Radic Res
1998
;
29
:
573
–9.
49
Bohr VA, Stevnsner T, de Souza-Pinto NC. Mitochondrial DNA repair of oxidative damage in mammalian cells.
Gene
2002
;
286
:
127
–34.
50
Thorslund T, Sunesen M, Bohr VA, Stevnsner T. Repair of 8-oxoG is slower in endogenous nuclear genes than in mitochondrial DNA and is without strand bias.
DNA Repair
2002
;
1
:
261
–73.
51
Cappelli E, Hazra T, Hill JW, Slupphaug G, Bogliolo M, Frosina G. Rates of base excision repair are not solely dependent on levels of initiating enzymes.
Carcinogenesis
2001
;
22
:
387
–93.
52
Nishioka K, Ohtsubo T, Oda H, et al. Expression and differential intracellular localization of two major forms of human 8-oxoguanine DNA glycosylase encoded by alternatively spliced OGG1 mRNAs.
Mol Biol Cell
1999
;
10
:
1637
–52.
53
ESCODD, Gedik CM, Collins A. Establishing the background level of base oxidation in human lymphocyte DNA: results of an interlaboratory validation study.
FASEB J
2005
;
19
:
82
–4.
54
Olinski R, Gackowski D, Rozalski R, Foksinski M, Bialkowski K. Oxidative DNA damage in cancer patients: a cause or a consequence of the disease development?
Mutat Res
2003
;
531
:
177
–90.
55
Dizdaroglu M. Base-excision repair of oxidative DNA damage by DNA glycosylases.
Mutat Res
2005
;
591
:
45
–59. doi:10.1016/j.mrfmmm.2005.01.033.
56
Hazra TK, Izumi T, Boldogh I, et al. Identification and characterization of a human DNA glycosylase for repair of modified bases in oxidatively damaged DNA.
Proc Natl Acad Sci U S A
2002
;
99
:
3523
–8.
57
Dou H, Mitra S, Hazra TK. Repair of oxidized bases in DNA bubble structures by human DNA glycosylases NEIL1 and NEIL2.
J Biol Chem
2003
;
278
:
49679
–84.
58
Parsons JL, Zharkov DO, Dianov GL. NEIL1 excises 3′ end proximal oxidative DNA lesions resistant to cleavage by NTH1 and OGG1.
Nucleic Acids Res
2005
;
33
:
4849
–56.
59
Lindahl T. Instability and decay of the primary structure of DNA.
Nature
1993
;
362
:
709
–15.
60
Nakamura J, La DK, Swenberg JA. 5′-nicked apurinic/apyrimidinic sites are resistant to β-elimination by β-polymerase and are persistent in human cultured cells after oxidative stress.
J Biol Chem
2000
;
275
:
5323
–8.
61
Simonelli V, Narciso L, Dogliotti E, Fortini P. Base excision repair intermediates are mutagenic in mammalian cells.
Nucleic Acids Res
2005
;
33
:
4404
–11.
62
Sobol RW, Kartalou M, Almeida KH, et al. Base excision repair intermediates induce p53-independent cytotoic and genotoxic responses.
J Biol Chem
2003
;
278
:
39951
–9.
63
Lindahl T, Nyberg B. Rate of depurination of native deoxyribonucleic acid.
Biochemistry
1972
;
11
:
3610
–8.
64
Nakamura J, Walker VE, Upton PB, Chiang SY, Kow YW, Swenberg JA. Highly sensitive apurinic/apyrimidinic site assay can detect spontaneous chemically induced depurination under physiological conditions.
Cancer Res
1998
;
58
:
222
–5.
65
Tsai K, Hsu TG, Hsu KM, et al. Oxidative DNA damage in human peripheral leukocytes induced by massive aerobic exercise.
Free Radic Biol Med
2001
;
31
:
1465
–72.
66
Pouget JP, Douki T, Richard MJ, Cadet J. DNA damage induced in cells by γ and UVA radiation as measured by HPLC/GC-MS and HPLC-EC and Comet assay.
Chem Res Toxicol
2000
;
13
:
541
–9.
67
Kennedy LJ, Moore K, Jr., Caulfield JL, Tannenbaum SR, Dedon PC. Quantitation of 8-oxoguanine and strand breaks produced by four oxidizing agents.
Chem Res Toxicol
1997
;
10
:
386
–92.
68
Xu J, Zheng SL, Turner A, et al. Associations between hOGG1 sequence variants and prostate cancer susceptibility.
Cancer Res
2002
;
62
:
2253
–7.
69
Kohno T, Shinmura K, Tosaka M, et al. Genetic polymorphisms and alternative splicing of the hOGG1 gene, that is involved in the repair of 8-hydroxyguanine in damaged DNA.
Oncogene
1998
;
25
:
3219
–25.
70
Mateuca R, Aka PV, De Boeck M, Hauspie R, Kirsch-Volders M, Lison D. Influence of hOGG1, XRCC1 and XRCC3 genotypes on biomarkers of genotoxicity in workers exposed to cobalt or hard metal dusts.
Toxicol Lett
2005
;
156
:
277
–88.
71
Park J, Chen L, Tockman MS, Elahi A, Lazarus P. The human 8-oxogunaine DNA N-glycosylase 1 (hOGG1) DNA repair enzyme and its association with lung cancer risk.
Pharmacogenetics
2004
;
14
:
103
–9.
72
Trzeciak AR, Nyaga SG, Jaruga P, Lohani A, Dizdaroglu M, Evans MK. Cellular repair of oxidatively induced DNA base lesions is defective in prostate cancer cell lines, PC-3 and DU-145.
Carcinogenesis
2004
;
25
:
1359
–70.
73
Aka P, Mateuca R, Buchet JP, Thierens H, Kirsch-Volders M. Are genetic polymorphisms in OGG1, XRCC1 and XRCC3 genes predictive for the DNA strand break repair phenotype and genotoxicity in workers exposed to low dose ionising radiation?
Mutat Res
2004
;
556
:
169
–81.
74
Peng T, Shen HM, Liu ZM, et al. Oxidative DNA damage in peripheral leukocytes and its association with expression and polymorphisms of hOGG1: a study of adolescents in a high risk region for hepatocellular carcinoma in China.
World J Gastroenterol
2003
;
9
:
2186
–93.
75
Hung RJ, Hall J, Brennan P, Boffetta P. Genetic polymorphisms in the base excision repair pathway and cancer risk: a HuGE review.
Am J Epidemiol
2005
;
162
:
1
–18.
76
Wikman H, Risch A, Klimek F, et al. hOGG1 polymorphism and loss of heterozygosity (LOH): significance for lung cancer susceptibility in a caucasian population.
Int J Cancer
2000
;
88
:
932
–7.
77
Kim SR, Matsui K, Yamada M, et al. Suppression of chemically induced and spontaneously occurring oxidative mutagenesis by three alleles of human OGG1 gene encoding 8-hydroxyguanine DNA glycosylase.
Mutat Res
2004
;
554
:
365
–74.
78
Janssen K, Schlink K, Gotte W, Hippler B, Kaina B, Oesch F. DNA repair activity of 8-oxoguanine DNA glycosylase 1 (OGG1) in human lymphocytes is not dependent on genetic polymorphism Ser326/Cys326.
Mutat Res
2001
;
486
:
207
–16.
79
Hansen R, Saebo M, Skjelbred CF, et al. GPX Pro198Leu and OGG1 Ser326Cys polymorphisms and risk of development of colorectal adenomas and colorectal cancer.
Cancer Lett
2005
;
229
:
85
–91.
80
Weiss JM, Goode EL, Ladiges WC, Ulrich CM. Polymorphic variation in hOGG1 and risk of cancer: a review of the functional and epidemiological literature.
Mol Carcinog
2005
;
42
:
127
–41.
81
Hardie LJ, Briggs JA, Davidson LA, et al. The effect of hOGG1 and glutathione peroxidase I genotypes and 3p chromosomal loss on 8-hydroxydeoxyguanosine levels in lung cancer.
Carcinogenesis
2000
;
21
:
167
–72.
82
Hagiwara A, Kitajima Y, Sato S, Miyazaki K. Allelic loss of the DNA repair gene OGG1 against oxidative damage in esophageal squamous cell carcinoma.
Oncol Rep
2005
;
13
:
1009
–16.
83
Bhattacharyya N, Banerjee S. A variant of DNA polymerase β acts as a dominant negative mutant.
Proc Natl Acad Sci U S A
1997
;
94
:
10324
–9.
84
Iwanaga A, Ouchida M, Miyazaki K, Hori K, Mukai T. Functional mutation of DNA polymerase β found in human gastric cancer—inability of the base excision repair in vitro.
Mutat Res
1999
;
435
:
121
–8.
85
Vens C, Dahman-Mooren E, Verwijs-Janssen M, et al. The role of DNA pollymerase β in determining sensitivity to ionizing radiation in human tumor cells.
Nucleic Acids Res
2002
;
30
:
2995
–3004.
86
Husain I, Morton BS, Beard WA, et al. Specific Inhibition of DNA polymerase β by its 14 kDa domain: role of single- and double-stranded DNA binding and 5′-phosphate recognition.
Nucleic Acids Res
1995
;
23
:
1597
–603.
87
Lang T, Maitra M, Starcevic D, Li SX, Sweasy JB. A DNA polymerase β mutant from colon cancer cells induces mutations.
Proc Natl Acad Sci U S A
2004
;
101
:
6074
–9.
88
Starcevic D, Dalal S, Sweasy JB. Is there a link between DNA polymerase β and cancer?
Cell Cycle
2004
;
3
:
998
–1001.
89
Prigent C, Satoh MS, Daly G, Barnes DE, Lindahl T. Aberrant DNA repair and DNA replication due to an inherited enzymatic defect in human DNA ligase I.
Mol Cell Biol
1994
;
14
:
310
–7.
90
Hao B, Wang H, Zhou K, et al. Identification of genetic variants in base excision repair pathway and their associations with risk of esophageal squamous cell carcinoma.
Cancer Res
2004
;
64
:
4378
–84.
91
Joseph T, Kusumakumary P, Chacko P, Abraham A, Pillai MR. DNA repair gene XRCC1 polymorphisms in childhood aculte lymphoblastic leukemia.
Cancer Lett
2005
;
217
:
17
–24.
92
Hu JJ, Smith TR, Miller MS, Lohman K, Case LD. Genetic regulation of ionizing radiation sensitivity and breast cancer risk.
Environ Mol Mutagen
2002
;
39
:
208
–15.
93
Ito H, Matsuo K, Hamajima N, et al. Gene-environment interactions between the smoking habit and polymorphisms in the DNA repair genes, APE1 Asp148Glu and XRCC1 Arg399Glu, in Japanese cancer risk.
Carcinogenesis
2004
;
25
:
1395
–401.
94
Hong YC, Lee KH, Kim WC, et al. Polymorphisms of XRCC1 gene, alcohol consumption and colorectal cancer.
Int J Cancer
2005
;
116
:
428
–32.
95
Matsuo K, Hamajima N, Suzuki R, et al. Lack of association between DNA base excision repair gene XRCC1 Arg399Gln polymorphism and risk of malignant lymphoma in Japan.
Cancer Genet Cytogenet
2004
;
149
:
77
–89.
96
Huang WY, Olshan AF, Schwartz SM, et al. Selected genetic polymorphisms in MGMT, XRCC1, XPD, and XRCC3 and risk of head and neck cancer: a pooled analysis.
Cancer Epidemiol Biomarkers Prev
2005
;
14
:
1747
–53.
97
Hu Z, Ma H, Chen F, Wei Q, Shen H. XRCC1 polymorphisms and cancer risk: a meta-analysis of 38 case-control studies.
Cancer Epidemiol Biomarkers Prev
2005
;
33
:
1810
–8.
98
Al-Tassan N, Chmiel NH, Maynard J, et al. Inherited variants of MYH associated with somatic G:C to T:A mutations in colorectal tumors.
Nat Genet
2002
;
30
:
227
–32.
99
Croitoru ME, Cleary SP, Di Nicola N, et al. Association between biallelic and monoallelic germline MYH gene mutations and colorectal cancer risk.
J Natl Cancer Inst
2004
;
96
:
1631
–4.
100
Sieber OM, Lipton L, Crabtree M, et al. Multiple colorectal adenomas, classic adenomatous polyposis, and germ-line mutations in MYH.
N Engl J Med
2003
;
348
:
791
–9.
101
Sampson JR, Jones S, Dolwani S, Cheadle JP. MutYH (MYH) and colorectal cancer.
Biochem Soc Trans
2005
;
33
:
679
–83.
102
Bai H, Jones S, Guan X, Wilson TM, Sampson JR, Cheadle JP. Functional characterization of two human MutY homolog (hMYH) missense mutations (R227W and V232F) that lie within the putative hMSH6 binding domain and are associated with hMYH polyposis.
Nucleic Acids Res
2005
;
33
:
597
–604.
103
Chmiel NH, Livingston AL, David SS. Insight into the functional consequences of inherited variants of the hMYH adenine glycosylase associated with colorectal cancer: complementation assays with hMYH variants and pre-steady-state kinetics of the corresponding mutated E. coli enzymes.
J Mol Biol
2003
;
327
:
431
–43.
104
Erichsen HC, Chanock SJ. SNPs in cancer research and treatment.
Br J Cancer
2004
;
90
:
747
–51.
105
Sharma RA, Farmer PB. Biological relevance of adduct detection to the chemoprevention of cancer.
Clin Cancer Res
2004
;
10
:
4901
–12.
106
Yu HJA, Lin KM, Ota DM, Lynch HT. Herediatry nonpolyposis colorectal cancer: preventative management.
Cancer Treat Rev
2003
;
29
:
461
–70.
107
Garcea G, Dennison AR, Steward WP, Berry DP. Chemoprevention of gastrointestinal malignancies.
ANZ J Surg
2003
;
73
:
680
–6.
108
Lenning PE. Genes causing inherited cancer as beacons to identify the mechanisms of chemoresistance.
Trends Mol Med
2004
;
10
:
113
–8.
109
Leipold MD, Workman H, Muller JG, Burrows CJ, David SS. Recognition and removal of oxidized guanines in duplex DNA by the base excision repair enzymes hOGG1, yOGG1, and yOGG2.
Biochemistry
2003
;
43
:
11373
–81.
110
Asagoshi K, Yamada T, Terato H, et al. Distinct repair activities of human 7,8-dihydro-8-oxoguanine DNA glycosylase and formamidopyrimidine DNA glycosylase for formamidopyrimidine and 7,8-dihydro-8-oxoguanine.
J Biol Chem
2000
;
275
:
4956
–64.
111
Mokkapati SK, Wiederhold L, Hazra TK, Mitra S. Stimulation of DNA glycosylase activity of OGG1 by NEIL1: functional collaboration between two human DNA glycosylases.
Biochemistry
2004
;
43
:
11596
–604.
112
Hu J, Imam SZ, Hashiguchi K, Souza-Pinto NC, Bohr VA. Phosphorylation of human oxoguanine DNA glycosylase (α-OGG1) modulates its function.
Nucleic Acids Res
2005
;
33
:
3271
–82.
113
Shen JC, Loeb LA. Mutations in the α-8 loop of human APE1 alter binding and cleavage of DNA containing an abasic site.
J Biol Chem
2003
;
278
:
46994
–7001.
114
Suh D, Wilson DM III, Povirk LF. 3′-phosphodiesterase activity of human apurinic/apyrimidinic endonuclease at DNA double-stranded break ends.
Nucleic Acids Res
1997
;
25
:
2495
–500.
115
Xu Y, DeMott MS, Hwang JT, Greenberg MM, Demple B. Action of human apurinic endonuclease (Ape1) on C1′-oxidized deoxyribose damage in DNA.
DNA Repair
2003
;
2
:
175
–85.
116
Strauss PR, Beard WA, Patterson TA, Wilson SH. Substrate binding by human apurinic/apyrimidinic endonuclease indicates a Briggs-Haldane mechanism.
J Biol Chem
1997
;
272
:
1302
–7.
117
Mol CD, Izumi T, Mitra S, Tainer JA. DNA-bound structures and mutants reveal abasic DNA binding by APE1 DNA repair and coordination.
Nature
2000
;
403
:
451
–6.
118
Hang B, Rothwell D, Sagi J, Hickson ID, Singer B. Evidence for a common active site for cleavage of an AP site and the benzene-derived exocyclic adduct, 3, N4-Benzethono-dC, in the major human AP endonuclease.
Biochemistry
1997
;
36
:
15411
–8.
119
Izumi T, Malecki J, Chaudhry MA, et al. Intragenic suppression of an active site mutation in the human apurinic/apyrimidinic endonuclease.
J Mol Biol
1999
;
287
:
47
–57.
120
Rothwell DG, Hang B, Gorman MA, Freemont PS, Singer B, Hickson ID. Substitution of Asp-210 in HAP1 (APE/Ref-1) eliminates endonuclease activity but stabilises substrate binding.
Nucleic Acids Res
2000
;
28
:
2207
–13.
121
Chagovetz AM, Sweasy JB, Preston BD. Increased activity and fidelity of DNA polymerase β on single-nucleotide gapped DNA.
J Biol Chem
1997
;
272
:
27501
–4.
122
Srivastava DK, Vande Berg BJ, Prasad R, et al. Mammalian abasic site base excision repair.
J Biol Chem
1998
;
273
:
21203
–9.
123
Beard WA, Osheroff WP, Prasad R, et al. Enzyme-DNA interactions required for efficient nucleotide incorporation and discrimination in human DNA polymerase β.
J Biol Chem
1996
;
271
:
12141
–4.
124
Liu Y, Beard WA, Shock DD, Prasad R, Hou EW, Wilson SH. DNA polymerase β and flap endonuclease 1 enzymatic specificities sustain DNA synthesis for long patch base excision repair.
J Biol Chem
2004
;
280
:
3665
–74.
125
Prasad R, Beard WA, Strauss PR, Wilson SH. Human DNA polymerase β deoxyribose phosphate lyase.
J Biol Chem
1998
;
273
:
15263
–70.
126
Bennett RAO, Wilson DM III, Wong D, Demple B. Interaction of human apurinic endonuclease and DNA polymerase β in the base excision repair pathway.
Proc Natl Acad Sci U S A
1997
;
94
:
7166
–9.
127
Wong D, Demple B. Modulation of the 5′-deoxyribose-5-phosphate lyase and DNA synthesis activities of mammalian DNA polymerase β by apurinic/apyrimidinic endonuclease I.
J Biol Chem
2004
;
279
:
25268
–75.
128
Teraoka H, Sawai M, Tsukada K. Kinetic studies on the reaction catalyzed by DNA ligase from calf thymus.
Biochim Biophys Acta
1983
;
747
:
117
–22.
129
Ranalli TA, DeMott MS, Bambara RA. Mechanism underlying replication protein A stimulation of DNA ligase I.
J Biol Chem
2002
;
277
:
1719
–27.
130
Tom S, Henricksen LA, Park MS, Bambara RA. DNA ligase I and proliferating cell nuclear antigen form a functional complex.
J Biol Chem
2001
;
276
:
24817
–25.
131
Einolf HJ, Guengerich FP. Kinetic analysis of nucleotide incorporation by mammalian DNA polymerase δ.
J Biol Chem
2000
;
275
:
16316
–22.
132
Einolf HJ, Guengerich FP. Fidelity of nucleotide insertion at 8-oxo-7,8-dihydroguanine by mammalian DNA polymerase β.
J Biol Chem
2001
;
276
:
3764
–71.
133
Mozzherin DJ, Fisher PA. Human DNA polymerase epsilon: enzymologic mechanism and gap-filling synthesis.
Biochemistry
1996
;
35
:
3572
–7.
134
Lee BI, Wilson DM III. The RAD2 domain of human exonuclease 1 exhibits 5′ to 3′ exonuclease and flap structure-specific endonuclease activities.
J Biol Chem
1999
;
274
:
37763
–9.
135
Danova I, Bohr VA, Dianov GL. Interaction of human AP endonculease 1 with flap endonuclease 1 and proliferating cell nuclear antigen involved in long patch repair.
Biochemistry
2001
;
40
:
12639
–44.
136
Saitoh T, Shinmura K, Yamaguchi S, et al. Enhancement of OGG1 protein AP lyase activity by increase of APEX protein.
Mutat Res
2000
;
486
:
31
–40.
137
Wiederhold L, Leppard JB, Kedar P, et al. AP endonuclease-independent DNA base excision repair in human cells.
Mol Cell
2004
;
15
:
209
–20.
138
Winters TA, Henner WD, Russell PS, McCullough A, Jorgensen TJ. 3′-phosphoglycolate from DNA strand-break damage in an oligonucleotide substrate by recombinant human apurinic/apyrimidinic endonuclease 1.
Nucleic Acids Res
1994
;
22
:
1866
–73.
139
Chen DS, Herman T, Demple B. Two distinct human DNA diesterases that hydrolyze 3′-blocking deoxyribose fragments from oxidized DNA.
Nucleic Acids Res
1991
;
19
:
5907
–14.
140
Horton JK, Srivastava DK, Zmudzka BZ, Wilson SH. Strategic down-regulation of DNA polymerase β by antisense RNA sensitizes mammalian cells to specific DNA damaging agents.
Nucleic Acids Res
1995
;
23
:
3810
–5.
141
Tomkinson AE, Lasko DD, Daly G, Lindahl T, Mammalian DNA ligases, Catalytic domain and size of DNA ligase I.
J Biol Chem
1990
;
265
:
12611
–7.