Cancer Research CTRC-AACR San Antonio Breast Cancer Symposium Tumor Immunology: New Perspectives
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wright, B. E.
Right arrow Articles by Reschke, D. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wright, B. E.
Right arrow Articles by Reschke, D. K.
[Cancer Research 62, 5641-5644, October 15, 2002]
© 2002 American Association for Cancer Research

Advances in Brief

Hypermutable Bases in the p53 Cancer Gene Are at Vulnerable Positions in DNA Secondary Structures1

Barbara E. Wright2, Jacqueline M. Reimers, Karen H. Schmidt and Dennis K. Reschke

Division of Biological Sciences University of Montana, Missoula, Montana 59812


   ABSTRACT
Top
 ABSTRACT
Introduction
 Results and Discussion
 REFERENCES
 
A DNA folding analysis indicates that the most hypermutable bases in exons 5, 7, and 8 of the p53 tumor suppressor gene are located immediately next to stems in stable DNA stem-loop structures. On the basis of the highest negative energy (-{Delta}G) value of the structures containing each mutable bases and on the extent to which each base is unpaired during transcription, their relative mutabilities are calculated using a new computer algorithm. These predicted mutation frequencies correlate well with those observed in 14,000 human cancers (R2 = 0.76), whereas there is no such correlation (R2 = 0.0005) for nearby control bases. The correlation of hypermutable base frequencies with -{Delta}G values is poor (R2 = 0.19), indicating that the extent to which a base is unpaired during transcription is a significant contribution to predicting mutation frequencies.


   Introduction
Top
 ABSTRACT
 Introduction
Results and Discussion
 REFERENCES
 
Chromosomal DNA is negatively supercoiled, and supercoiling creates and stabilizes SLSs3 that are sequence determined, i.e., they form as a consequence of inverted complementary sequences in ssDNA (1 , 2) . These SLSs contain regions of ssDNA with unpaired and mispaired bases vulnerable to mutation (3 , 4) . SLSs not only expose specific bases to mutation but also increase the time of their exposure by periodically blocking the progress of transcribing RNA polymerase complexes during transcript elongation. Such pausing is also known to be sequence specific and to promote mutagenesis (5 , 6) . Transcription enhances supercoiling, which is highly localized (7, 8, 9) . The direct effects of various kinds of stress on supercoiling have also been shown to be quite specific within the genome and to affect primarily those genes related to the stress involved. Examples are starvation (10) ,4 reactive oxygen species (11) , metal toxicity (12) , virulence (13) , and osmotic shock (14) . Transcription-enhanced mutations would therefore occur preferentially at pause sites in SLSs of supercoiled DNA formed in the wake of the advancing transcription bubble. If specific types of stress directly target particular areas of the genome for enhanced supercoiling, the formation of secondary structures will be facilitated; if stress first activates transcription, this will in turn cause localized supercoiling. By either scenario, the resulting increase in concentration of SLSs results in localized hypermutation. A mechanism that targets increased mutation rates to genes that must mutate to overcome prevailing conditions of stress has likely been selected during evolution (15) . Direct evidence linking increased rates of transcription with localized hypermutation in unpaired bases of SLSs has been obtained recently, in Escherichia coli.4 Here, evidence is presented suggesting that the hypermutable bases in p53 are also uniquely vulnerable to mutation because of their location and position in high negative energy SLSs.


   Results and Discussion
Top
 ABSTRACT
 Introduction
 Results and Discussion
REFERENCES
 
Stress-induced Mutations in E. coli.
Localized supercoiling and hypermutation should occur as a consequence of starvation-induced derepression and transcriptional activation (15) . To test this hypothesis, mutations in derepressed genes of starving E. coli cells were analyzed and compared with mutations in cells not under stress.4 Zuker’s DNA folding computer program (16) was used to reveal the most thermodynamically stable SLSs and their negative-free energies of folding (-{Delta}G) for each nucleotide sequence analyzed. The higher the -{Delta}G value, the more likely a stable SLS will form and result in a pause site. Using Zuker’s program, 30, 40, and 50 nt segments were folded, and the 30 nt folds showed the most specific correlation of maximum peak -{Delta}G values with the presence of each mutant nt in the derepressed genes. A high -{Delta}G per se is due primarily to the stability of the stem(s) and does not necessarily indicate the presence of a mutable base. To predict the mutability of a base, it is also necessary to know the extent to which the base is unpaired in all its foldings during transcription (3 , 4 , 15) . We therefore developed an interfacing program that calculates: (a) the percentage of folds in which each base is unpaired; and (b) a Mutability Index for each base, which is expressed as an absolute value and defined as follows: MI = (percentage of folds in which the base is unpaired) x (highest -{Delta}G of all folds in which it is unpaired; http://biology.dbs.umt.edu/wright/upload/mfg.html).

Thus, MIs indicate the relative mutation frequencies of bases in predicted SLSs. The highest -{Delta}G term gives the most weight to the most stable fold containing the unpaired base. That is, primary importance is given to the structure in which the base is exposed for the longest period of time and therefore has the highest probability of mutation.

Mutable bases in derepressed genes are preferentially located in high -{Delta}G SLSs compared with the unstressed mutable bases. In the mutable bases of transcribing genes, there is a good correlation between relative MIs and mutation frequencies determined experimentally, whereas there is no correlation for mutable bases in unstressed cells.4 The results of these investigations in prokaryotes are consistent with the location of hypermutable bases in high -{Delta}G SLSs of transcribing genes, suggesting that environmental stressors may cause mutations in such genes. The relevance of these findings to higher organisms has now been investigated in p53. A database of sequences for point mutations in the human p53 gene is available (17) for analyzing DNA secondary structure as a possible cause of hypermutation and cancer.

Hypermutable p53 Codons in High -{Delta}G SLSs.
Exons 5–8 in p53 were analyzed with the DNA folding program to determine whether the hypermutable bases are located in high -{Delta}G regions of the gene. The size of the folding window used appeared to be less critical than in prokaryotes because folding sequences between 35 and 50 nts gave essentially the same results. About 100 nts containing hypermutable bases (exons 5, 7, and 8) or not containing hypermutable bases (exon 6) were folded in successive, overlapping 43 nt windows, beginning at each fifth nt (data not shown). The average -{Delta}Gs for the regions, including exons 5, 7, and 8, with the hypermutable bases are 6.2, 5.6, and 4.4, respectively, as compared with 3.0 for exon 6 with no hypermutable bases.

The individual -{Delta}Gs for the SLS containing the hypermutable bases in the first and second positions of codons 175, 245, 248, 249, 273, and 282 are given in Table 1Citation (the third position is silent, except for 249). Also shown in this table are the predicted MIs for each base compared with its observed mutation frequency recorded from 14,000 sequences of the database. As seen in Fig. 1Citation , the correlation between observed and predicted mutation frequencies is significant (R2 = 0.40) using a 35 nt window (C) and highly significant (R2 = 0.76 and 0.63) for 43 nt (A) and 50 nt (D) windows, respectively. These results suggest that, in vivo, the average number of nts folded during transcription is in the range of 40–50 nts. The correlation between MIs and mutation frequencies presumably depends upon both of the variables involved, i.e., the -{Delta}Gs of the relevant SLSs, as well as the extent to which the bases are unpaired during transcription. To test this assumption, mutation frequencies were plotted against -{Delta}G using a 43 nt window (Fig. 1E)Citation . The R2 value of 0.19 is poor compared with 0.76, indicating the contribution of both variables to the predictive value of MIs. The highest -{Delta}G SLS containing the most mutable base in each of the six codons are shown in Fig. 2Citation . In each case, the unpaired base is immediately next to a stem. Proximity to a stem also characterized the more mutable bases in a similar analysis of prokaryotes.4 The last column of Table 1Citation indicates that 9 of 12 hypermutable bases are located immediately next to a stem, compared with 5 of 24 control mutable bases that were arbitrarily chosen based on their proximity to the hypermutable codons (Table 2)Citation . No correlation is seen (R2 = 0.0005) between predicted and observed mutation frequencies in these control mutable bases (Fig. 1B)Citation .


View this table:
[in this window]
[in a new window]
  		Table 1 Comparison of mutation frequencies, MIs, and -{Delta}Gs in p53 hypermutable codons

 


View larger version (15K):
[in this window]
[in a new window]
  		Fig. 1. Linear regression analyses of the data for hypermutable bases (Table 1)Citation and control bases (Table 2)Citation are shown. MIs determined by folding 43 nt windows are plotted against mutation frequencies of hypermutable bases (A) and control bases (B). Windows of 35 nts (C) and 50 nts (D) are also shown for the hypermutable bases. The mutation frequency is also plotted against -{Delta}G alone (E) to contrast with the fit obtained in (A).

 


View larger version (13K):
[in this window]
[in a new window]
  		Fig. 2. SLSs of the most mutable bases in each of the hypermutable codons. Blue numbers and lines indicate codons, red arrows point to the mutable bases, and the 5' end of the SLSs are indicated by a green dot.

 

View this table:
[in this window]
[in a new window]
  		Table 2 Comparison of mutation frequencies, MIs, and -{Delta}G in p53 control codons

 
The Mechanism and Consequences of Mutability.
The process of creating SLSs during transcription is robust. Regardless of probable continuous variations in the length of the nt segments folded, the location of mutable bases in these structures generally coincide with peak -{Delta}G values. In this dynamic and ever-changing scene, the stems, which create SLSs, are the invariant part. Thus, the most reliable location for an unpaired base that evolved to be mutable is immediately next to a stem, as seen in most of the hypermutable bases (Table 1)Citation . Except for one or perhaps two nts next to a stem, ssDNA segments at the beginning or end of a SLS may, in reality, be paired to their normal complement in double-stranded DNA. As this is not recognized by the computer program, MI values may not be valid for such bases. This would be the case for the C of codon 175 and the G of codon 282 (SLS not shown) but would not change the total mutation frequencies of these codons significantly. Ten of the control bases would be affected in this manner.

The data suggest that exposure time of the hypermutable bases during transcription is a critical variable in the complex of events causing human cancers. Proximity to a stem in a high -{Delta}G SLS may help explain why the hypermutable bases in p53 are so mutable and why they are targets for various kinds of environmental stress such as nitric oxide and oxyradicals (18 , 19) . Mutagens are known to stimulate p53 transcription (20 , 21) as well as supercoiling (11 , 12) , and the MI of a base may dictate the extent to which a toxin can affect mutation frequency. Because stressors cannot affect MIs directly, their effects on mutation frequency may also be indirect by enhancing supercoiling and SLSs directly or indirectly via transcriptional activation (10, 11, 12, 13, 14, 15) . By the same token, suppressing such mutations may require inhibition of the supercoiling and/or transcription involved.

Whereas MIs determine relative mutation frequencies, the spectrum of mutations produced is apparently determined by the chemical mechanisms involved. This conclusion is substantiated by the fact that the spectrum of p53 mutations that occurs is similar to that of background E. coli lacI mutations (22) : C mutates to T much more frequently than to A or G, and G mutates to A much more frequently than to C or T. The thermodynamic properties of hydrolytic reactions that occur in nucleic acids under physiological conditions are such that the deamination of C is much greater than A, and depurination of G or A is much greater than depyrimidation of C or T (reviewed in Ref. 23 ). Also, the oxidation potential of a G located 5' to another G (codons 245, 248, and 249) is greater than it is next to a C or T (reviewed in Ref. 24 ). Because of its size, an A is more likely than a C or T to replace a G. The total G and C content of the SLS containing the bases in Table 1Citation is 32 and 30%, respectively (data not shown). In Table 1Citation , 91% of the bases are Gs or Cs compared with 50% of the control bases (Table 2)Citation . In contrast to the hypermutable bases, no Gs and only one C in the control bases are immediately next to a stem (codon 250, Table 2Citation ). The data suggest that while mutagens increase mutation frequency, their specificity for targeting Gs or Cs cancel out in a database of 14,000 mutations such that the spectrum is similar to that of background mutations. The spectrum of mutants produced in p53 or lacI is independent of MI. Whether or not a mutagen is present, the underlying cause of relative mutation frequencies in p53 appears to be the inevitable consequence of folding the sequences that constitute this gene. Most, if not all, mutagens are more reactive on ssDNA, and all of the mechanisms by which they are presumed to cause mutation occur at lesser (but significant) levels in an untreated genome because of the inherent nature of the bond involved. Thus, it is more appropriate to use the term background rather than "spontaneous" mutation (23) .

In conclusion, it appears as though GC richness, thermodynamic potential of depurination, deamination, and oxidation, as well as proximity to a stem conspired during evolution to create the hypermutable bases of the p53 gene. In unicellular organisms, adverse environmental conditions can bestow mutability on those genes regulating the cell’s response to each type of stress, thus minimizing genome-wide genetic damage while creating the most appropriate variants for selection of the fittest in response to stress (15) . In multicellular organisms, however, variants are not independent cells but part of a whole, and selection of the fittest results in abnormal growth. This would be an ironic legacy for mechanisms that were selected to accelerate evolution.


   ACKNOWLEDGMENTS
 
We thank Drs. Mike Minnick, Steve Lodmell, Frank Rosenzweig, Brooke Martin, and Kent Sugden for valuable suggestions and discussions.


   FOOTNOTES
 
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 Supported by NIH Grants R15CA88893 and R55CA99242 and the Stella Duncan Memorial Research Institute. Back

2 To whom requests for reprints should be addressed, at Phone: (406) 243-6676; Fax: (406) 243-4184; E-mail: bewright{at}selway.umt.edu Back

3 The abbreviations used are: SLS, stem-loop structure; ssDNA, single-stranded DNA; nt, nucleotide; MI, mutability index; -{Delta}G, negative energy of folding (kcal/mole). Back

4 B. E. Wright, D. K. Reschke, K. H. Schmidt, J. M. Reimers, and W. Knight. Predicting mutation frequencies in stem-loop structures of derepressed genes: implications for evolution, submitted for publication. Back

Received 7/12/02. Accepted 9/ 3/02.


   REFERENCES
Top
 ABSTRACT
 Introduction
 Results and Discussion
 REFERENCES
 

  1. Lilley D. M. J. The inverted repeat as a recognizable structural feature in supercoiled molecules. Proc. Natl. Acad. Sci. USA, 77: 6468-6472, 1980.[Abstract/Free Full Text]
  2. Panayotatos N., Wells R. D. Cruciform structures in supercoiled DNA. Nature (Lond.), 289: 446-470, 1981.[Medline]
  3. Farabaugh P. J., Schmeisser U., Hofer M., Miller J. H. Genetic studies on the lac repressor. VII. On the molecular nature of spontaneous hot spots in the lacI gene of Escherichia coli. J. Mol. Biol., 126: 847-863, 1978.[Medline]
  4. Ripley L. S., Glickman B. W. Unique self-complementarity of palindromic sequences provides DNA structural intermediates for mutation. Cold Spring Harbor Symp. Quant. Biol., 47: 851-861, 1983.
  5. Bebenek K., Abbotts J., Wilson S. H., Kunkel T. A. Error-prone polymerization by HIV-1 reverse transcriptase. J. Biol. Chem., 268: 10324-10334, 1993.[Abstract/Free Full Text]
  6. Pananicolaou C., Ripley L. S. An in vitro approach to identifying specificity determinants of mutagenesis mediated by DNA misalignments. J. Mol. Biol., 221: 805-821, 1991.[Medline]
  7. Liu L., Wang J. C. Supercoiling and the DNA template during transcription. Proc. Natl. Acad. Sci. USA, 84: 7024-7027, 1987.[Abstract/Free Full Text]
  8. Dayn A., Malkhosyan S., Mirkin S. M. Transcriptionally driven cruciform formation in vivo. Nucleic Acids Res., 20: 5991-5997, 1992.[Abstract/Free Full Text]
  9. Rahmouni A. R., Wells R. D. Direct evidence for the effect of transcription on local DNA supercoiling in vivo. J. Mol. Biol., 223: 131-144, 1992.[Medline]
  10. Wright B. E., Longacre A., Reimers J. M. Hypermutation in derepressed operons of Escherichia coli K12. Proc. Natl. Acad. Sci. USA, 96: 5089-5094, 1999.[Abstract/Free Full Text]
  11. Weinstein-Fischer D., Elgrably-Weiss M., Altuvia S. Escherichia coli response to hydrogen peroxide: a role for DNA supercoiling, Topoisomerase I and Fis. Mol. Microbiol., 35: 1413-1420, 2000.[Medline]
  12. Ansari A. Z., Chael M. L., O’Halloran T. V. Allosteric underwinding of DNA is a critical step in positive control of transcription by Hg-MerR. Nature (Lond.), 355: 87-89, 1992.[Medline]
  13. Dorman C. J. DNA topology and the global control of bacterial gene expression: Implications for the regulation of virulence gene expression. Microbiology, 141: 1271-1280, 1995.[Medline]
  14. Jordi B. J. A. M., Owen-Hughes T. A., Hulton C. S. J., Higgins C. F. DNA twist, flexibility and transcription of the osmoregulated proU promoter of Salmonella typhimurium. EMBO J., 14: 5690-5700, 1995.[Medline]
  15. Wright B. E. A biochemical mechanism for nonrandom mutations and evolution. J. Bacteriol., 182: 2993-3001, 2000.[Free Full Text]
  16. Zuker M., Mathews D. H., Turner D. H. A practical guide Barciszewski J. Clark B. F. C. eds. . RNA Biochemistry and Biotechnology, 11-43, Kluwer Acad. Publ. Dordrecht, the Netherlands 1999.
  17. Oliver, M., Eeles, R., Hollstein, M., Khan, M. A., Harris, C. C., and Hainaut, P. The IARC TP53 database: new online mutation analysis and recommendations to users. Hum. Mutat., in press, 2002.
  18. Wink D. A., Kasprzak K. S., Margos C. M., Elespuru R. K., Misra M., Dumams T. M., Cebula T. A., Koch W. H., Andrews A. W., Allen J. S., Keefer L. K. DNA deaminating ability and genotoxicity of nitric oxide and its progenitors. Science (Wash. DC), 254: 1001-1003, 1991.[Abstract/Free Full Text]
  19. Hussain S. P., Aguilar F., Amstad P., Cerutti P. Oxy-radical induced mutagenesis of hotspot codons 248 and 249 of the human p53 gene. Oncogene, 9: 2277-2281, 1994.[Medline]
  20. Rogel A., Popliker M., Webb C. G., Oren M. p53 cellular tumor antigen: analysis of mRNA levels in normal adult tissues, embryos, and tumors. Mol. Cell. Biol., 5: 2851-2855, 1985.[Abstract/Free Full Text]
  21. Reisman D., Loging W. T. Transcriptional regulation of the p53 tumor suppressor gene. Semin. Cancer Biol., 8: 317-324, 1998.[Medline]
  22. Schaaper R. M., Dunn R. L. Spontaneous mutation in the Escherichia coli lacI gene. Genetics, 129: 317-326, 1991.[Abstract]
  23. Singer B., Kusmierek J. T. Chemical mutagenesis. Annu. Rev. Biochem., 52: 655-693, 1982.
  24. Sugden K. D., Stearns D. M. The role of chromium (V) in the mechanism of chromate-induced oxidative DNA damage and cancer. J. Environ. Pathol. Toxicol. Oncol., 19: 215-230, 2000.[Medline]



This article has been cited by other articles:


Home page
 MicrobiologyHome page
E. Burkala, J. M. Reimers, K. H. Schmidt, N. Davis, P. Wei, and B. E. Wright
Secondary structures as predictors of mutation potential in the lacZ gene of Escherichia coli
Microbiology, July 1, 2007; 153(7): 2180 - 2189.
[Abstract] [Full Text] [PDF]

Home page
 J. Bacteriol.Home page
C. Ross, C. Pybus, M. Pedraza-Reyes, H.-M. Sung, R. E. Yasbin, and E. Robleto
Novel Role of mfd: Effects on Stationary-Phase Mutagenesis in Bacillus subtilis
J. Bacteriol., November 1, 2006; 188(21): 7512 - 7520.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
O. Amosova, R. Coulter, and J. R. Fresco
Self-catalyzed site-specific depurination of guanine residues within gene sequences
PNAS, March 21, 2006; 103(12): 4392 - 4397.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
R. D. Wells, R. Dere, M. L. Hebert, M. Napierala, and L. S. Son
Advances in mechanisms of genetic instability related to hereditary neurological diseases
Nucleic Acids Res., July 8, 2005; 33(12): 3785 - 3798.
[Abstract] [Full Text] [PDF]

Home page
 J. Mol. Diagn.Home page
M. Akbari, M. Dore Hansen, J. Halgunset, F. Skorpen, and H. E. Krokan
Low Copy Number DNA Template Can Render Polymerase Chain Reaction Error Prone in a Sequence-Dependent Manner
J. Mol. Diagn., February 1, 2005; 7(1): 36 - 39.
[Abstract] [Full Text] [PDF]

Home page
 J. Biol. Chem.Home page
A. Bacolla and R. D. Wells
Non-B DNA Conformations, Genomic Rearrangements, and Human Disease
J. Biol. Chem., November 12, 2004; 279(46): 47411 - 47414.
[Full Text] [PDF]

Home page
 MicrobiologyHome page
J. M. Reimers, K. H. Schmidt, A. Longacre, D. K. Reschke, and B. E. Wright
Increased transcription rates correlate with increased reversion rates in leuB and argH Escherichia coli auxotrophs
Microbiology, May 1, 2004; 150(5): 1457 - 1466.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wright, B. E.
Right arrow Articles by Reschke, D. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wright, B. E.
Right arrow Articles by Reschke, D. K.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online