The use of non-tumor data in cancer risk assessment: reflections on butadiene, vinyl chloride, and benzene

Rationale and Roadmap for Developing Panels of Hotspot Cancer Driver Gene Mutations as Biomarkers of Cancer Risk

Cancer driver mutations (CDMs) are necessary and causal for carcinogenesis and have advantages as reporters of carcinogenic risk. However, little progress has been made toward developing measurements of CDMs as biomarkers for use in cancer risk assessment. Impediments for using a CDM-based metric to inform cancer risk include the complexity and stochastic nature of carcinogenesis, technical difficulty in quantifying low-frequency CDMs, and lack of established relationships between cancer driver mutant fractions and tumor incidence. Through literature review and database analyses, this review identifies the most promising targets to investigate as biomarkers of cancer risk. Mutational hotspots were discerned within the 20 most mutated genes across the 10 deadliest cancers. Forty genes were identified that encompass 108 mutational hotspot codons overrepresented in the COSMIC database; 424 different mutations within these hotspot codons account for approximately 63,000 tumors and their prevalence across tumor types is described. The review summarizes literature on the prevalence of CDMs in normal tissues and suggests such mutations are direct and indirect substrates for chemical carcinogenesis, which occurs in a spatially stochastic manner. Evidence that hotspot CDMs (hCDMs) frequently occur as tumor subpopulations is presented, indicating COSMIC data may underestimate mutation prevalence. Analyses of online databases show that genes containing hCDMs are enriched in functions related to intercellular communication. In its totality, the review provides a roadmap for the development of tissue-specific, CDM-based biomarkers of carcinogenic potential, comprised of batteries of hCDMs and can be measured by error-correct next-generation sequencing. Environ. Mol. Mutagen. 61:152-175, 2020. Published 2019. This article is a U.S. Government work and is in the public domain in the USA. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Major groove orientation of the (2S)-N(6)-(2-hydroxy-3-buten-1-yl)-2'-deoxyadenosine DNA adduct induced by 1,2-epoxy-3-butene

1,3-Butadiene (BD) is an environmental and occupational toxicant classified as a human carcinogen. It is oxidized by cytochrome P450 monooxygenases to 1,2-epoxy-3-butene (EB), which alkylates DNA. BD exposures lead to large numbers of mutations at A:T base pairs even though alkylation of guanines is more prevalent, suggesting that one or more adenine adducts of BD play a role in BD-mediated genotoxicity. However, the etiology of BD-mediated genotoxicity at adenine remains poorly understood. EB alkylates the N(6) exocyclic nitrogen of adenine to form N(6)-(hydroxy-3-buten-1-yl)-2'-dA ((2S)-N(6)-HB-dA) adducts ( Tretyakova , N. , Lin , Y. , Sangaiah , R. , Upton , P. B. , and Swenberg , J. A. ( 1997 ) Carcinogenesis 18 , 137 - 147 ). The structure of the (2S)-N(6)-HB-dA adduct has been determined in the 5'-d(C(1)G(2)G(3)A(4)C(5)Y(6)A(7)G(8)A(9)A(10)G(11))-3':5'-d(C(12)T(13)T(14)C(15)T(16)T(17)G(18)T(19) C(20)C(21)G(22))-3' duplex [Y = (2S)-N(6)-HB-dA] containing codon 61 (underlined) of the human N-ras protooncogene, from NMR spectroscopy. The (2S)-N(6)-HB-dA adduct was positioned in the major groove, such that the butadiene moiety was oriented in the 3' direction. At the Cα carbon, the methylene protons of the modified nucleobase Y(6) faced the 5' direction, which placed the Cβ carbon in the 3' direction. The Cβ hydroxyl group faced toward the solvent, as did carbons Cγ and Cδ. The Cβ hydroxyl group did not form hydrogen bonds with either T(16) O(4) or T(17) O(4). The (2S)-N(6)-HB-dA nucleoside maintained the anti conformation about the glycosyl bond, and the modified base retained Watson-Crick base pairing with the complementary base (T(17)). The adduct perturbed stacking interactions at base pairs C(5):G(18), Y(6):T(17), and A(7):T(16) such that the Y(6) base did not stack with its 5' neighbor C(5), but it did with its 3' neighbor A(7). The complementary thymine T(17) stacked well with both 5' and 3' neighbors T(16) and G(18). The presence of the (2S)-N(6)-HB-dA resulted in a 5 °C reduction in the Tm of the duplex, which is attributed to less favorable stacking interactions and adduct accommodation in the major groove.

Increased mitochondrial DNA copy number in occupations associated with low-dose benzene exposure

BACKGROUND

Benzene is an established leukemogen at high exposure levels. Although low-level benzene exposure is widespread and may induce oxidative damage, no mechanistic biomarkers are available to detect biological dysfunction at low doses.

OBJECTIVES

Our goals were to determine in a large multicenter cross-sectional study whether low-level benzene is associated with increased blood mitochondrial DNA copy number (mtDNAcn, a biological oxidative response to mitochondrial DNA damage and dysfunction) and to explore potential links between mtDNAcn and leukemia-related epigenetic markers.

METHODS

We measured blood relative mtDNAcn by real-time polymerase chain reaction in 341 individuals selected from various occupational groups with low-level benzene exposures (> 100 times lower than the Occupational Safety and Health Administration/European Union standards) and 178 referents from three Italian cities (Genoa, Milan, Cagliari).

RESULTS

In each city, benzene-exposed participants showed higher mtDNAcn than referents: mtDNAcn was 0.90 relative units in Genoa bus drivers and 0.75 in referents (p = 0.019); 0.90 in Milan gas station attendants, 1.10 in police officers, and 0.75 in referents (p-trend = 0.008); 1.63 in Cagliari petrochemical plant workers, 1.25 in referents close to the plant, and 0.90 in referents farther from the plant (p-trend = 0.046). Using covariate-adjusted regression models, we estimated that an interquartile range increase in personal airborne benzene was associated with percent increases in mtDNAcn equal to 10.5% in Genoa (p = 0.014), 8.2% (p = 0.008) in Milan, 7.5% in Cagliari (p = 0.22), and 10.3% in all cities combined (p < 0.001). Using methylation data available for the Milan participants, we found that mtDNAcn was associated with LINE-1 hypomethylation (-2.41%; p = 0.007) and p15 hypermethylation (+15.95%, p = 0.008).

CONCLUSIONS

Blood MtDNAcn was increased in persons exposed to low benzene levels, potentially reflecting mitochondrial DNA damage and dysfunction.

Evaluating chromosomal damage in workers exposed to hexavalent chromium and the modulating role of polymorphisms of DNA repair genes

PURPOSE

Welders have been chronically exposed to hexavalent chromium with potential consequences on chromosomal integrity. Our study is focused on the extent of any such chromosomal aberrations with respect to chromium levels in the blood of welders as well as on the tentative modulating role of polymorphisms in DNA repair genes XPD Lys751Gln, XPG Asn114His, XPC Lys939Gln, hOGG1 Ser326Cys and XRCC1 Arg399Gln on chromosomal damage.

METHODS

The study was conducted on 144 individuals consisting of 73 welders exposed to chromium for 10.2 ± 1.67 years and 71 control individuals without known exposures. Chromosomal aberrations, their chromatid-type and chromosome-type aberrations were detected by conventional cytogenetic analysis. XPD, XPG, XPC, hOGG1 and XRCC1 gene polymorphisms were assayed for by Taqman SNP genotyping assay ("Assay-by-Demand") using Real-Time allelic discrimination on AB 7500 equipment. Chromium concentration in the blood was determined by atomic absorption spectrophotometry.

RESULTS

The level of chromium in the blood of welders ranged between 0.032 and 0.182 μmol l(-1) and was significantly higher than that in controls (0.07 ± 0.04 μmol l(-1) vs. 0.03 ± 0.007 μmol l(-1)). Parameters of chromosomal damage were similar in both the exposed and the control individuals (1.89% vs. 1.70% for total chromosomal aberrations, 0.97% vs. 0.88% for chromosome-type and 0.92% vs. 0.80% for chromatid-type, respectively). Chromatid-type of aberrations positively correlated with the level of chromium in the blood (r = 0.28; P = 0.02). Significantly higher total chromosomal aberrations were detected in individuals with homozygous variant polymorphism in XRCC1 Arg399Gln gene as compared to those with heterozygous and homozygous wild-type genotypes (2.20, 1.89 and 1.48%, respectively; P = 0.01). A similar tendency was found for chromatid-type aberrations (1.30% for homozygous variant genotype bearers, 0.94% for those with heterozygous genotype and 0.75% for carriers of homozygous wild-type genotype, respectively; P = 0.04).

CONCLUSIONS

Although no apparent increase in chromosomal damage was recorded in chromium-exposed welders in comparison with controls, genetic make-up in DNA repair genes may increase susceptibility toward adverse effect of chromium.

Development and application of an LC-MS/MS method for the detection of the vinyl chloride-induced DNA adduct N(2),3-ethenoguanine in tissues of adult and weanling rats following exposure to [(13)C(2)]-VC

In the 1970s, exposure to vinyl chloride (VC) was shown to cause liver angiosarcoma in VC workers. We have developed a new LC-MS/MS method for analyzing the promutagenic DNA adduct N(2),3-ethenoguanine (εG) and have applied this to DNA from tissues of both adult and weanling rats exposed to 1100 ppm [(13)C(2)]-VC for 5 days or 1100 ppm VC for 1 day. This assay utilizes neutral thermal hydrolysis and an HPLC cleanup prior to quantitation by LC-MS/MS. The number of endogenous and exogenous εG adducts in DNA from tissues of adult rats exposed to [(13)C(2)]-VC for 5 days was 4.1 ± 2.8 adducts/10(8) guanine of endogenous and 19.0 ± 4.9 adducts/10(8) guanine of exogenous εG in the liver, 8.4 ± 2.8 adducts/10(8) guanine of endogenous and 7.4 ± 0.5 adducts/10(8) guanine of exogenous εG in the lung, and 5.9 ± 3.3 adducts/10(8) guanine of endogenous and 5.7 ± 2.1 adducts/10(8) guanine of exogenous εG in the kidney (n = 4). Additionally, the data from weanling rats demonstrated higher numbers of exogenous εG, with ∼4-fold higher amounts in the liver DNA of weanlings (75.9 ± 17.9 adducts/10(8) guanine) in comparison to adult rats and ∼2-fold higher amounts in the lung (15.8 ± 3.6 adducts/10(8) guanine) and kidney (12.9 ± 0.4 adducts/10(8) guanine) (n = 8). The use of stable isotope labeled VC permitted accurate estimates of the half-life of εG for the first time by comparing [(13)C(2)]-εG in adult rats with identically exposed animals euthanized 2, 4, or 8 weeks later. The half-life of εG was found to be 150 days in the liver and lung and 75 days in the kidney, suggesting little or no active repair of this promutagenic adduct.

1,3-Butadiene: III. Assessing carcinogenic modes of action

1,3-Butadiene (BD) is a multisite carcinogen in laboratory rodents following lifetime exposure, with greater potency in the mouse than the rat, and is associated with an increase in leukemia mortality in highly exposed workers. Species differences in the formation of reactive metabolites underlie observed species differences in sensitivity to the carcinogenic effects of BD. The modes of action (MOAs) for human leukemia and rodent tumors are both likely related to mutagenic potencies of one or more of these metabolites. However, differences in the nature of genotoxic lesions associated with human leukemia and rodent tumors, along with their implications for risk assessment, require that they be discussed separately. The MOAs for BD are assessed in this review using the modified Hill criteria and human relevance framework. Key events in MOAs for human and rodent cancers are identified, along with important species differences and sources of nonlinearity for each event that can affect extrapolations made from high- to low-dose exposures. Because occupational exposures to BD have also included co-exposures to styrene and dimethyldithiocarbamide (DMDTC), potential interactions with BD carcinogenicity are also discussed. The MOAs for BD carcinogenesis will be used to guide key decisions made in the quantitative cancer dose-response assessment.

Creating context for the use of DNA adduct data in cancer risk assessment: I. Data organization

The assessment of human cancer risk from chemical exposure requires the integration of diverse types of data. Such data involve effects at the cell and tissue levels. This report focuses on the specific utility of one type of data, namely DNA adducts. Emphasis is placed on the appreciation that such DNA adduct data cannot be used in isolation in the risk assessment process but must be used in an integrated fashion with other information. As emerging technologies provide even more sensitive quantitative measurements of DNA adducts, integration that establishes links between DNA adducts and accepted outcome measures becomes critical for risk assessment. The present report proposes an organizational approach for the assessment of DNA adduct data (e.g., type of adduct, frequency, persistence, type of repair process) in concert with other relevant data, such as dosimetry, toxicity, mutagenicity, genotoxicity, and tumor incidence, to inform characterization of the mode of action. DNA adducts are considered biomarkers of exposure, whereas gene mutations and chromosomal alterations are often biomarkers of early biological effects and also can be bioindicators of the carcinogenic process.

Advances in understanding benzene health effects and susceptibility

Benzene is a ubiquitous chemical in our environment that causes acute leukemia and probably other hematological cancers. Evidence for an association with childhood leukemia is growing. Exposure to benzene can lead to multiple alterations that contribute to the leukemogenic process, indicating a multimodal mechanism of action. Research is needed to elucidate the different roles of multiple metabolites in benzene toxicity and the pathways that lead to their formation. Studies to date have identified a number of polymorphisms in candidate genes that confer susceptibility to benzene hematotoxicity. However, a genome-wide study is needed to truly assess the role of genetic variation in susceptibility. Benzene affects the blood-forming system at low levels of occupational exposure, and there is no evidence of a threshold. There is probably no safe level of exposure to benzene, and all exposures constitute some risk in a linear, if not supralinear, and additive fashion.

Stochasticity in physiologically based kinetics models: implications for cancer risk assessment

In case of low-dose exposure to a substance, its concentration in cells is likely to be stochastic. Assessing the consequences of this stochasticity in toxicological risk assessment requires the coupling of macroscopic dynamics models describing whole-body kinetics with microscopic tools designed to simulate stochasticity. In this article, we propose an approach to approximate stochastic cell concentration of butadiene in the cells of diverse organs. We adapted the dynamics equations of a physiologically based pharmacokinetic (PBPK) model and used a stochastic simulator for the system of equations that we derived. We then coupled kinetics simulations with a deterministic hockey stick model of carcinogenicity. Stochasticity induced substantial modifications relative to dose-response curve, compared with the deterministic situation. In particular, there was nonlinearity in the response and the stochastic apparent threshold was lower than the deterministic one. The approach that we developed could easily be extended to other biological studies to assess the influence of stochasticity at macroscopic scale for compound dynamics at the cell level.

Vinyl Chloride—A Classical Industrial Toxicant of New Interest

The carcinogenicity of vinyl chloride in humans was recognized in 1974 based on observations of hepatic angiosarcomas in highly exposed workers. A multiplicity of endpoints has been demonstrated. The primary target organ, the liver, displays differential susceptibilities of hepatocytes and sinusoidal cells, which are modified by factors of age and dose. There is consistency in organotropism between experimental animals and humans. Vinyl chloride is a pluripotent carcinogen, predominantly directed toward hepatic endothelial (sinusoidal) cells, and second toward the parenchymal cells of the liver. The similarity of results between experimental animals and humans is a solid basis of an amalgamation of experimental and epidemiological risk estimates. Vinyl chloride requires metabolic activation for carcinogenicity and mutagenicity, and toxicokinetics are a key to interpret the dose response. Practically the entire initial metabolism of vinyl chloride is oxidative. At higher exposure concentrations this is nonlinear, and metabolic saturation of metabolism in rats is reached at about 250 ppm. This is consistent with the plateau of hepatic angiosarcoma incidence in rat bioassays. Physiologically based pharmacokinetic/toxicokinetic (PBPK) models have been developed and successfully applied within the frame of human cancer risk assessments. The major DNA adduct induced by vinyl chloride (approximately 98% of total adducts in rats), 7-(2-oxoethyl)guanine, is almost devoid of promutagenic activity. The clearly promutagenic "etheno" adducts N2,3-ethenoguanine and 3,N4-ethenocytosine each represent approximately 1% of the vinyl chloride DNA adducts in rats, and 1,N6-ethenoadenine is found at even lower concentrations. Etheno adducts appear to have a long persistence and are repaired by glycosylases. Vinyl chloride represents a human carcinogen for which a series of mechanistic events connects exposure with the carcinogenic outcome. These include (1) metabolic activation (to form chloroethylene oxide), (2) DNA binding of the reactive metabolite (to exocyclic etheno adducts), (3) promutagenicity of these adducts, and (4) effects of such mutations on protooncogenes/tumor suppressor genes at the gene and gene product levels. In rat hepatocytes, a further event is a biomarker response. Cancer prestages (enzyme-altered foci), as quantitative biomarkers, provide a tool to study dose response even within low dose ranges where a carcinogenic risk cannot be seen in cancer bioassays directly. Such biomarker responses support a linear nonthreshold extrapolation for low-dose assessment of carcinogenic risks due to vinyl chloride. Published risk estimates based on different sets of data (animal experiments, epidemiological studies) appear basically consistent, and on this basis an angiosarcoma risk of approximately 3 x 10(-4) has been deduced by extrapolation, for exposure to 1 ppm vinyl chloride over an entire human working lifetime. An important point that should be considered in regulatory standard settings is the presence of a physiological background of those etheno DNA adducts, which are also produced by vinyl chloride. Likely reasons for this background are oxidative stress and lipid peroxidation. In essence, fundamentals of the hepatocarcinogenicity of vinyl chloride appear now well established, providing a solid scientific basis for regulatory activities.

Chromosomal aberrations in tire plant workers and interaction with polymorphisms of biotransformation and DNA repair genes

We evaluated chromosomal aberrations in lymphocytes of 177 workers exposed to xenobiotics in a tire plant and in 172 controls, in relation to their genetic background. Nine polymorphisms in genes encoding biotransformation enzymes and nine polymorphisms in genes involved in main DNA repair pathways were investigated for possible modulation of chromosomal damage. Chromosomal aberration frequencies were the highest among exposed smokers and the lowest in non-smoking unexposed individuals (2.5+/-1.8% vs. 1.7+/-1.2%, respectively). The differences between groups (ANOVA) were borderline significant (F=2.6, P=0.055). Chromosomal aberrations were higher in subjects with GSTT1-null (2.4+/-1.7%) than in those with GSTT1-plus genotype (1.8+/-1.4%; F=7.2, P=0.008). Considering individual groups, this association was significant in smoking exposed workers (F=4.4, P=0.040). Individuals with low activity EPHX1 genotype exhibited significantly higher chromosomal aberrations (2.3+/-1.6%) in comparison with those bearing medium (1.7+/-1.2%) and high activity genotype (1.5+/-1.2%; F=4.7, P=0.010). Both chromatid- and chromosome-type aberration frequencies were mainly affected by exposure and smoking status. Binary logistic regression analysis revealed that frequencies of chromatid-type aberrations were modulated by NBS1 Glu185Gln (OR 4.26, 95%CI 1.38-13.14, P=0.012), and to a moderate extent, by XPD Lys751Gln (OR 0.16, 95%CI 0.02-1.25, P=0.081) polymorphisms. Chromosome-type aberrations were lowest in individuals bearing the EPHX1 genotype conferring the high activity (OR 0.38, 95%CI 0.15-0.98, P=0.045). Present results show that exposed individuals in the tire production, who smoke, exhibit higher chromosomal aberrations frequencies, and the extent of chromosomal damage may additionally be modified by relevant polymorphisms.

HPRT mutations in lymphocytes from 1,3-butadiene-exposed workers in China

BACKGROUND

1,3-Butadiene (BD) is an important industrial chemical and an environmental and occupational pollutant. The carcinogenicity of BD in rodents has been proved, but its carcinogenic and mutagenic molecular mechanism(s) are not fully elucidated in humans.

OBJECTIVES

In the present study, we compared the mutation frequencies and exon deletions of BD-exposed workers with that of control subjects in China to identify the characteristic mutations associated with BD exposure in the human HPRT (hypoxanthine-guanine-phosphoribosyltransferase) gene.

METHODS

Seventy-four workers exposed to BD via inhalation and 157 matched controls were evaluated in Nanjing, China. Molecular analysis of HPRT mutant T lymphocytes from BD-exposed workers and nonexposed control subjects was conducted to identify changes in the structure of the HPRT gene. A total of 783 HPRT mutants were analyzed by multiplex polymerase chain reaction, in which 368 HPRT mutants were isolated from BD-exposed workers and 415 mutants from control subjects.

RESULTS

The BD-exposed workers showed a higher mutation frequency (18.2 +/- 9.4 x 10(-6)) than the control subjects (12.7 +/- 7.3 x 10(-6)), but the difference was not significant (p > 0.05). The frequency of exon deletions in BD-exposed workers (27.4%) was significantly higher than that in control subjects (12.5%) (p < 0.05), which mainly included multiplex exon deletions (2-8 exons).

CONCLUSIONS

The results of the present study suggest that BD should increase the frequency of large deletions of HPRT gene in human lymphocytes This change confirms and supports the previous findings in BD-exposed workers.

Regulatory assessment of chemicals within OECD member countries, EU and in Russia

The chemical risk assessment is essesntial part of new chemical legislation registration, evaluation, and authorization of chemicals (REACH). The article presents a review of chemical legislation policies in the European Union (EU) and in Russia, and changes in chemicals regulations to meet the requirement of REACH. The risk assessment paradigm, toxicological parameters, safe limits and classification criteria used by different agencies and authorities in different countries are reported. Our investigation also focuses on comparison of chemical risk assessment criteria used in OECD member countries and in Russia. Tendencies in harmonization in accordance with the globally harmonized system of classification and labeling of chemicals (GHS) are discussed.

Propylene oxide: genotoxicity profile of a rodent nasal carcinogen

Propylene oxide (PO) is a DNA-reactive genotoxic agent; that is, it reacts with DNA to produce lesions in the genetic material. PO also induces tumors in rodents, although only at high concentrations and at portals of entry. This review of PO's genotoxicity profile is organized according to endpoints measured, that is, nonmutational or mutational endpoints, and as to whether the results were from in vitro or in vivo studies. In addition to results of experimental studies, PO's genotoxicity for humans is assessed by reviewing results of published biomarker studies. The weight of evidence indicates that although it is genotoxic, PO's potency as a DNA-reactive mutagen is weak. Other aspects of PO's overall tissue toxicities are also reviewed, with attention to glutathione (GSH) depletion and its consequences, that is, cell proliferation, death, and necrosis. These toxic tissue responses occur in the same anatomical regions in rodents as do the PO-induced tumors. Furthermore, some of these tissue toxicities can produce effects that may either augment PO's DNA-reactive mutagenicity or be genotoxic in themselves, not dependent on PO's DNA reactivity. Although its DNA reactivity may be a necessary component of PO's overall genotoxicity and rodent carcinogenicity, it is likely not sufficient, and the associated tissue toxicities, which are rate-limiting, also seem to be required. This complex mode of action has implications for estimations of PO's cancer potential in humans, especially at low exposure concentrations.

The effect of benzene on the activity of adenosine deaminase in tissues of rats

Swiss Albino (Rat rattus norvegicus) rats were intraperitoneally injected with a 100 mg kg(-1) dosage of benzene, a toxic and carcinogenic agent widely used for industrial purposes. Changes in the adenosine deaminase (ADA) activity in the liver, kidney and serum of rats were investigated at 0, 2, 4, 8, 16, 32 and 64 h following injection. Serum physiological was administered to each control group. Enzyme activities were measured spectrophotometrically. Our purpose was to further investigations of some diseases caused by benzene, and present evidence of variations in the activity of ADA enzyme effected by benzene. While benzene caused significant inhibitions in ADA activity in the liver at 16 and 32 h and at 0.05 probability level, no significant inhibition or activation occurred at other test periods (hours). ADA activity did not present any significant variation in the kidneys. It was observed that ADA activity displayed similar patterns in the control groups. Comparisons of ADA activities in the two groups showed a statistically significant decrease between 4(th) and 64(th) hours (p< 0.05), demonstrating a direct correlation between benzene and its effects on ADA enzymes.

1,3-Butadiene exposure and cardiovascular disease

This review summarizes the epidemiologic, biochemical and genetic evidence associating occupational, environmental or experimental exposure to 1,3-butadiene (BD) with subsequent development of cardiovascular disease, with the primary focus on atherosclerosis. The potential role of BD in the known atherosclerotic effects of environmental tobacco smoke as well as correlations between polymorphisms in BD phase II enzymes and development of atherosclerosis are presented.

Structure of the 1,4-Bis(2'-deoxyadenosin-N(6)-yl)-2S,3S-butanediol intrastrand DNA cross-link arising from butadiene diepoxide in the human N-ras codon 61 sequence

The 1,4-bis(2'-deoxyadenosin-N(6)-yl)-2S,3S-butanediol intrastrand DNA cross-link arises from the bis-alkylation of tandem N(6)-dA sites in DNA by R,R-butadiene diepoxide (BDO(2)). The oligodeoxynucleotide 5'-d(C(1)G(2)G(3)A(4)C(5)X(6)Y(7)G(8)A(9)A(10)G(11))-3'.5'-d(C(12)T(13)T(14)C(15)T(16)T(17)G(18)T(19)C(20)C(21)G(22))-3' contains the BDO(2) cross-link between the second and third adenines of the codon 61 sequence (underlined) of the human N-ras protooncogene and is named the (S,S)-BD-(61-2,3) cross-link (X,Y = cross-linked adenines). NMR analysis reveals that the cross-link is oriented in the major groove of duplex DNA. Watson-Crick base pairing is perturbed at base pair X(6).T(17), whereas base pairing is intact at base pair Y(7).T(16). The cross-link appears to exist in two conformations, in rapid exchange on the NMR time scale. In the first conformation, the beta-OH is predicted to form a hydrogen bond with T(16) O(4), whereas in the second, the beta-OH is predicted to form a hydrogen bond with T(17) O(4). In contrast to the (R,R)-BD-(61-2,3) cross-link in the same sequence (Merritt, W. K., Nechev, L. V., Scholdberg, T. A., Dean, S. M., Kiehna, S. E., Chang, J. C., Harris, T. M., Harris, C. M., Lloyd, R. S., and Stone, M. P. (2005) Biochemistry 44, 10081-10092), the anti-conformation of the two hydroxyl groups at C(beta) and C(gamma) with respect to the C(beta)-C(gamma) bond results in a decreased twist between base pairs X(6).T(17) and Y(7).T(16), and an approximate 10 degrees bending of the duplex. These conformational differences may account for the differential mutagenicity of the (S,S)- and (R,R)-BD-(61-2,3) cross-links and suggest that stereochemistry plays a role in modulating biological responses to these cross-links (Kanuri, M., Nechev, L. V., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580).

Molecular epidemiological studies in 1,3-butadiene exposed Czech workers: female-male comparisons

Results of a recent molecular epidemiological study of 1,3-butadiene (BD) exposed Czech workers, conducted to compare female to male responses, have confirmed and extended the findings of a previously reported males only study (HEI Research Report 116, 2003). The initial study found that urine concentrations of the metabolites 1,2-dihydroxy-4-(acetyl) butane (M1) and 1-dihydroxy-2-(N-acetylcysteinyl)-3-butene (M2) and blood concentrations of the hemoglobin adducts N-[2-hydroxy-3-butenyl] valine (HB-Val) and N-[2,3,4-trihydroxy-butyl] valine (THB-Val) constitute excellent biomarkers of exposure, both being highly correlated with BD exposure levels, and that GST genotypes modulate at least one metabolic pathway, but that irreversible genotoxic effects such as chromosome aberrations and HPRT gene mutations are neither associated with BD exposure levels nor with worker genotypes (GST [glutathione-S-transferase]-M1, GSTT1, CYP2E1 (5' promoter), CYP2E1 (intron 6), EH [epoxide hydrolase] 113, EH139, ADH [alcohol dehydrogenase]2 and ADH3). The no observed adverse effect level (NOAEL) for chromosome aberrations and HPRT mutations was 1.794 mg/m(3) (0.812 ppm)--the mean exposure level for the highest exposed worker group in this initial study. The second Czech study, reported here, initiated in 2003, included 26 female control workers, 23 female BD exposed workers, 25 male control workers and 30 male BD exposed workers (some repeats from the first study). Multiple external exposure measurements (10 full 8-h shift measures by personal monitoring per worker) over a 4-month period before biological sample collections showed that BD workplace levels were lower than in the first study. Mean 8-h TWA exposure levels were 0.008 mg/m(3) (0.0035 ppm) and 0.397 mg/m(3) (0.180 ppm) for female controls and exposed, respectively, but with individual single 8-h TWA values up to 9.793 mg/m(3) (4.45 ppm) in the exposed group. Mean male 8-h TWA exposure levels were 0.007 mg/m(3) (0.0032 ppm) and 0.808 mg/m(3) (0.370 ppm) for controls and exposed, respectively; however, the individual single 8-h TWA values up to 12.583 mg/m(3) (5.72 ppm) in the exposed group. While the urine metabolite concentrations for both M1 and M2 were elevated in exposed compared to control females, the differences were not significant, possibly due to the relatively low BD exposure levels. For males, with greater BD exposures, the concentrations of both metabolites were significantly elevated in urine from exposed compared to control workers. As in the first study, urine metabolite excretion patterns in both sexes revealed conjugation to be the minor detoxification pathway (yielding the M2 metabolite) but both M1 and M2 concentration values were lower in males in this second study compared to their concentrations in the first, reflecting the lower external exposures of males in this second study compared to the first. Of note, females showed lower concentrations of both M1 and M2 metabolites in the urine per unit of BD exposure than did males while exhibiting the same M1/(M1+M2) ratio, reflecting the same relative utilization of the hydrolytic (producing M1) and the conjugation (producing M2) detoxification pathways as males. Assays for the N,N-(2,3-dihydroxy-1,4-butadyl) valine (pyr-Val) hemoglobin (Hb) adduct, which is specific for the highly genotoxic 1,2,3,4-diepoxybutane (DEB) metabolite of BD, have been conducted on blood samples from all participants in this second Czech study. Any adduct that may have been present was below the limits of quantitation (LOQ) for this assay for all samples, indicating that production of this important BD metabolite in humans is below levels produced in both mice and rats exposed to as little as 1.0 ppm BD by inhalation (J.A. Swenberg, M.G. Bird, R.J. Lewis, Future directions in butadiene risk assessment, Chem. Biol. Int. (2006), this issue). Results of assays for the HB-Val and THB-Val hemoglobin adducts are pending. HPRT mutations, determined by cloning assays, and multiple measures of chromosome level changes (sister-chromatid exchanges [SCE], aberrations determined by conventional methods and FISH) again showed no associations with BD exposures, confirming the findings of the initial study that these irreversible genotoxic changes do not arise in humans occupationally exposed to low levels of BD. Except for lower production of both urine metabolites in females, no female-male differences in response to BD exposures were detected in this study. As in the initial study, there were no significant genotype associations with the irreversible genotoxic endpoints. However, as in the first, differences in the metabolic detoxification of BD as reflected in relative amounts of the M1 and M2 urinary metabolites were associated with genotypes, this time both GST and EH.

The 5'-GNC site for DNA interstrand cross-linking is conserved for diepoxybutane stereoisomers

The bifunctional alkylating agent 1,2,3,4-diepoxybutane forms interstrand DNA-DNA cross-links between the N7 positions of deoxyguanosine residues on opposite strands of the duplex. For racemic diepoxybutane, these cross-links predominate within 5'-GNC/3'CNG sequences, where N is any nucleotide. We used denaturing polyacrylamide gel electrophoresis (dPAGE) to examine the role of stereochemistry in the cross-linking reaction, subjecting a restriction fragment to cross-linking with S,S-DEB, R,R-DEB, or meso-DEB. DNA cross-links generated by each isomer were isolated by dPAGE, and the sites of cross-linking were identified by sequencing gel analysis of DNA fragments generated by hot piperidine cleavage. We found that the 5'-GNC consensus sequence of racemic DEB is conserved, but the efficiencies of cross-linking vary, with S,S- > R,R- > meso-DEB. These results help explain the observed differences between the biological activities of DEB stereoisomers.

Structure of the 1,4-bis(2'-deoxyadenosin-N6-yl)-2R,3R-butanediol cross-link arising from alkylation of the human N-ras codon 61 by butadiene diepoxide

The solution structure of the 1,4-bis(2'-deoxyadenosin-N(6)-yl)-2R,3R-butanediol cross-link arising from N(6)-dA alkylation of nearest-neighbor adenines by butadiene diepoxide (BDO(2)) was determined in the oligodeoxynucleotide 5'-d(CGGACXYGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained codon 61 (underlined) of the human N-ras protooncogene. The cross-link was accommodated in the major groove of duplex DNA. At the 5'-side of the cross-link there was a break in Watson-Crick base pairing at base pair X(6).T(17), whereas at the 3'-side of the cross-link at base pair Y(7).T(16), base pairing was intact. Molecular dynamics calculations carried out using a simulated annealing protocol, and restrained by a combination of 338 interproton distance restraints obtained from (1)H NOESY data and 151 torsion angle restraints obtained from (1)H and (31)P COSY data, yielded ensembles of structures with good convergence. Helicoidal analysis indicated an increase in base pair opening at base pair X(6).T(17), accompanied by a shift in the phosphodiester backbone torsion angle beta P5'-O5'-C5'-C4' at nucleotide X(6). The rMD calculations predicted that the DNA helix was not significantly bent by the presence of the four-carbon cross-link. This was corroborated by gel mobility assays of multimers containing nonhydroxylated four-carbon N(6),N(6)-dA cross-links, which did not predict DNA bending. The rMD calculations suggested the presence of hydrogen bonding between the hydroxyl group located on the beta-carbon of the four-carbon cross-link and T(17) O(4), which perhaps stabilized the base pair opening at X(6).T(17) and protected the T(17) imino proton from solvent exchange. The opening of base pair X(6).T(17) altered base stacking patterns at the cross-link site and induced slight unwinding of the DNA duplex. The structural data are interpreted in terms of biochemical data suggesting that this cross-link is bypassed by a variety of DNA polymerases, yet is significantly mutagenic [Kanuri, M., Nechev, L. V., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580].

Dual roles of glycosyl torsion angle conformation and stereochemical configuration in butadiene oxide-derived N1 beta-hydroxyalkyl deoxyinosine adducts: a structural perspective

The solution structure of the N1-[1-hydroxy-3-buten-2(R)-yl]-2'-deoxyinosine adduct arising from the alkylation of adenine N1 by butadiene epoxide (BDO), followed by deamination to deoxyinosine, was determined in the oligodeoxynucleotide 5'-d(CGGACXAGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained the BDO adduct at the second position of codon 61 of the human N-ras protooncogene (underlined) and was named the ras61 R-N1-BDO-(61,2) adduct. 1H NMR revealed a weak C5 H1' to X6 H8 nuclear Overhauser effects (NOE), followed by an intense X6 H8 to X6 H1' NOE. Simultaneously, the X6 H8 to X6 H3' NOE was weak. The resonances arising from the T16 and T17 imino protons were not observed. 1H NOEs between the butadiene moiety and the DNA positioned the adduct in the major groove. Structural refinement based upon a total of 394 NOE-derived distance restraints and 151 torsion angle restraints yielded a structure in which the modified deoxyinosine was in the syn conformation about the glycosyl bond, with a glycosyl bond angle of 83 degrees , and T17, the complementary nucleotide, was stacked into the helix but not hydrogen bonded with the adducted inosine. The refined structure provides a plausible hypothesis as to why these N1 deoxyinosine adducts strongly code for the incorporation of dCTP during trans lesion DNA replication, irrespective of stereochemistry, both in Escherichia coli [Rodriguez, D. A., Kowalczyk, A., Ward, J. B. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2001) Environ. Mol. Mutagen. 38, 292-296] and in mammalian cells [Kanuri, M., Nechev, L. N., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580]. Rotation of the N1 deoxyinosine adduct into the syn conformation may facilitate incorporation of dCTP via Hoogsteen type templating with deoxyinosine, generating A to G mutations. However, conformational differences between the R- and the S-N1-BDO-(61,2) adducts, involving the positioning of the butenyl moiety in the major groove of DNA, suggest that adduct stereochemistry plays a secondary role in modulating the biological response to these adducts.

Use of mode of action in risk assessment: past, present, and future

The evolution of chemical risk assessment has been marked by a steadily increasing expectation for the use of chemical-specific dosimetric and mechanistic information to tailor the risk assessment approach. The information to be used can range from the broad physical properties of the chemical to detailed information on the mechanism by which it causes a particular toxic outcome, and the risk assessment decisions effected can in turn range from how to define equivalent exposures across species to whether a particular animal outcome is relevant to a human health assessment. A concept that has proven useful in support of these considerations is the "mode of action," a term coined by the USEPA in their new guidelines for carcinogen risk assessment. This paper describes the increasing use of mode-of-action considerations in risk assessment, beginning with early examples involving quantitative dosimetry on the one hand, and qualitative relevance on the other, which foreshadowed the current interest in mode of action. It then describes more recent developments regarding the use of the mode-of-action concept for the selection of a low-dose extrapolation approach, for harmonization of cancer and noncancer risk assessment approaches, and for cross-chemical evaluations. Finally, examples of recent controversies associated with the use of mode-of-action information in risk assessment are provided to demonstrate the challenges that must be overcome to assure the continued viability of the mode-of-action approach.

Analysis of Diepoxide-Specific Cyclic N-Terminal Globin Adducts in Mice and Rats after Inhalation Exposure to 1,3-Butadiene

1,3-Butadiene is an important industrial chemical used in the production of synthetic rubber and is also found in gasoline and combustion products. It is a multispecies, multisite carcinogen in rodents, with mice being the most sensitive species. 1,3-Butadiene is metabolized to several epoxides that form DNA and protein adducts. Previous analysis of 1,2,3-trihydroxybutyl-valine globin adducts suggested that most adducts resulted from 3-butene-1,2-diol metabolism to 3,4-epoxy-1,2-butanediol, rather than from 1,2;3,4-diepoxybutane. To specifically examine metabolism of 1,3-butadiene to 1,2;3,4-diepoxybutane, the formation of the 1,2;3,4-diepoxybutane-specific adduct N,N-(2,3-dihydroxy-1,4-butadiyl)-valine was evaluated in mice treated with 3, 62.5, or 1250 ppm 1,3-butadiene for 10 days and rats exposed to 3 or 62.5 ppm 1,3-butadiene for 10 days, or to 1000 ppm 1,3-butadiene for 90 days, using a newly developed immunoaffinity liquid chromatography tandem mass spectrometry assay. In addition, 2-hydroxy-3-butenyl-valine and 1,2,3-trihydroxybutyl-valine adducts were determined. The analyses of several adducts derived from 1,3-butadiene metabolites provided new insight into species and exposure differences in 1,3-butadiene metabolism. Mice formed much higher amounts of N,N-(2,3-dihydroxy-1,4-butadiyl)-valine than rats. The formation of 2-hydroxy-3-butenyl-valine and N,N-(2,3-dihydroxy-1,4-butadiyl)-valine was similar in mice exposed to 3 or 62.5 ppm 1,3-butadiene, whereas 2-hydroxy-3-butenyl-valine was 3-fold higher at 1250 ppm. In both species, 1,2,3-trihydroxybutyl-valine adducts were much higher than 2-hydroxy-3-butenyl-valine and N,N-(2,3-dihydroxy-1,4-butadiyl)-valine. Together, these data show that 1,3-butadiene is primarily metabolized via the 3-butene-1,2-diol pathway, but that mice are much more efficient at forming 1,2;3,4-diepoxybutane than rats, particularly at low exposures. This assay should also be readily adaptable to molecular epidemiology studies on 1,3-butadiene-exposed workers.

Carcinogenicity categorization of chemicals-new aspects to be considered in a European perspective

Existing systems of classification of carcinogens are a matter of discussion, world-wide. There is agreement that it should be distinguished between genotoxic and non-genotoxic chemicals. The risk assessment approach used for non-genotoxic chemicals is similar among different regulatory bodies: insertion of an uncertainty (safety) factor permits the derivation of permissible exposure levels at which no relevant human cancer risks are anticipated. For genotoxic carcinogens, case studies of chemicals point to a whole array of possibilities. Positive data of chromosomal effects only, in the absence of mutagenicity, may support the characterization of a compound that produces carcinogenic effects only at high, toxic doses. Non-DNA-reactive genotoxins, such as topoisomerase inhibitors or inhibitors of the spindle apparatus are considered in this respect. In such cases, arguments are in favour of the existence of "practical" thresholds. Taking existing concepts together, it is proposed to basically distinguish between "perfect" and "practical" thresholds. There is a wide consensus that for non-DNA-reactive genotoxins such as aneugens (aneuploidy, chromosome loss, non-disjunction) thresholds should be defined. It is being discussed as to whether the identification of possible threshold effects should also include other mechanisms of genotoxicity, in addition to aneugenic effects. Specific mechanisms of clastogenicity have been repeatedly addressed as also having thresholds, such as topoisomerase II poisons or mechanisms based on reactive oxygen. Oxidative stress as an important mechanism is triggered by exposure to exogenous factors such as ultraviolet (UV) and ionizing radiation, anoxia and hyperoxia, and by chemicals producing reactive oxygen species. The idea is receiving increased support that reactive oxygen species (ROS)-mediated processes of carcinogenesis have practical thresholds. Since reactive oxygen species are genotoxic in principle, questions arise whether chemicals that increase ROS production will superimpose to an endogenously produced background level of DNA lesions, related to mechanisms that may result in non-linear dose-effect relationships. The existence of "endogenous" DNA adducts has been generally accepted, and possible regulatory implications of the presence of endogenous carcinogens have been discussed. It is now becoming evident that a diversity of methods of carcinogenic risk extrapolation to low doses must be considered, dependent on the mode of action. Although there is an increasing international awareness of these developments, the system of classification of carcinogens of the European Union still remains static. This should be changed, as the philosophy of separation of a strictly sequential "hazard assessment" and "risk assessment" appears out-of-date.

Diepoxybutane induces caspase and p53-mediated apoptosis in human lymphoblasts

Diepoxybutane (DEB) is the most potent metabolite of the environmental chemical 1,3-butadiene (BD), which is prevalent in petrochemical industrial areas. BD is a known mutagen and human carcinogen, and possesses multiorgan systems toxicity that includes bone marrow depletion, spleen, and thymus atrophy. Toxic effects of BD are mediated through its epoxy metabolites. In working towards elucidating the cellular and molecular mechanisms of BD toxicity, we investigated the ability of DEB to induce apoptosis in human lymphoblasts. DEB induced a concentration and exposure time-dependent apoptosis, which accounted for the DEB-induced loss of cell viability observed in TK6 lymphoblasts. The DEB-induced apoptosis was inhibited by inhibitors of caspases 3 and 9. The role of p53 in mediating the DEB-induced apoptosis was also investigated. DEB induced elevated p53 levels in direct correlation to the extent of DEB-induced apoptosis, as the concentration of DEB increased up to 5 microM. The extent of DEB-induced apoptosis was dramatically higher in TK6 lymphoblasts as compared to the genetically paired p53-deficient NH32 lymphoblasts under the same experimental conditions. Our results confirm and extend observations on the occurrence of apoptosis in DEB exposed cells, and demonstrate for the first time the elevation of p53 levels in human lymphoblasts in response to DEB exposure. In addition, our results demonstrate for the first time that DEB-induced apoptosis is mediated by caspases 3 and 9, as well as the p53 protein. It is possible that DEB-induced apoptosis may explain BD-induced bone marrow depletion, spleen and thymus atrophy in BD-exposed animals.