PrefaceThe Viracept (nelfinavir)—Ethyl methanesulfonate case: A threshold risk assessment for human exposure to a genotoxic drug contamination?☆
Introduction
Following reports by patients of a strange smell of Viracept® 250 mg tablets, analysis revealed that Roche's protease inhibitor nelfinavir mesylate produced between March 2007 and May 2007 was contaminated with up to 2300 ppm ethyl methanesulfonate (EMS), as a result of a production accident (Gerber and Toelle, 2009). This led to a global recall of the drug on June 6th 2007 and initiated an extensive testing strategy in close coordination with regulatory agencies (Müller and Singer, 2009) in order to establish a comprehensive risk assessment for an estimated daily EMS dose of 0.055 mg/kg/d in patients taking 10 tablets daily over 3 months (Müller et al., 2009).
EMS is a genotoxic carcinogen, inducing DNA damage and mutation by DNA ethylation. Although EMS has widely been used as model mutagen, the database regarding carcinogenic potency is poor and information on genotoxicity limited to high dose levels (Gocke et al., 2009b). An initial cancer risk assessment had therefore to be based on very limited data and the default procedure for extrapolation to low dose, i.e., a linear dose–response relationship. A maximum cancer risk based on lifetime exposure was estimated to be on the order of 10−3, a 3-month exposure was considered to be associated with less than one additional tumor in 10,000 exposed patients (Gocke et al., 2009c).
In view of the absence of reliable toxicity data, a large testing program was set up. This included a standard 4-week toxicity study in Wistar rats that showed a NOAEL of 20 mg/kg per day (Pfister and Eichinger-Chapelon, 2009). A major effort was made for in vivo genotoxicity studies, including the mouse bone marrow micronucleus test (MNT) and the MutaMouse assay for lacZ mutant induction. A wide dose range was covered and periods of exposure were up to 28 days (Gocke et al., 2009a). The rationale for these efforts was a report on threshold-like dose responses for a number of endpoints of genotoxicity by EMS in cultured cells (Doak et al., 2007). Ethylation of the N-terminal valine of mouse globin was determined in parallel (Gocke et al., 2009a). The data showed pronounced sublinearity in the dose–response (slope increasing with dose) for EMS-induced genotoxicity, while globin ethylation increased approximately linearly with dose. Doses up to 25 mg/kg per day (for lacZ mutant frequency) and up to 80 mg/kg per day (for micronucleus induction) did not produce any deviation from background measures in mouse bone marrow, while the slope above this putative threshold dose was prominent. Statistical assessment of the data (Gocke and Wall, 2009) indicated that the null hypothesis of linearity could be rejected with high significance against a hockey stick model (Lutz and Lutz, 2009). The distinct threshold-like appearance of the dose–response curve can be explained by saturation of the alkylguanine DNA alkyltransferase repair system above a certain threshold dose. This occurs whenever the rate of adduct formation exceeds the rate of instantaneous repair.
Under these conditions, the concentration of EMS at the DNA rather than the AUC becomes the critical dose metrics. In order to estimate cmax and AUC in humans for the estimated dose of about 0.05 mg/kg per day, a comprehensive characterization of the kinetics of EMS in vitro and in vivo was necessary (Lavé et al., 2009a), including translation of the data to exposed patients by appropriate modeling (Lavé et al., 2009b). This all done, a comprehensive human risk assessment was performed, based on the assumption of a threshold and using margins of exposure (Müller et al., 2009).
Contamination of drug mesylates by EMS is a general pharmaceutical problem, and a daily dose of 1.5 μg EMS per patient is considered the threshold of toxicological concern (TTC). In view of the threshold-type risk assessment described here, a permitted daily exposure (PDE) for EMS of 100 μg per person is proposed (Müller and Gocke, 2009). This exposure would still include a product safety factor of 12,000.
The last contribution to this special issue discusses impact and provides an outlook by external experts in the field (Walker et al., 2009). The authors recommend biologically based models in lieu of default models also for risk assessments for DNA-reactive compounds.
Section snippets
Types of “thresholds”
In the review process for the manuscripts in this issue, it became evident that lack of definition of the term threshold stirred up some controversy. For a good introduction, I recommend reading of (Slob, 1999). Fig. 1 shows three shapes of dose–response curves that could all be considered thresholds. They differ with respect to the initial slope of the curve above dose zero, being insignificantly positive (A), zero (B), or negative (C).
Type A is plausible biologically (Lutz, 1998) and is an
Conflict of interest
The author was a consultant to F. Hoffmann-La Roche Ltd. in connection with the Viracept-EMS incident.
Acknowledgements
Thanks go to the reviewers and to my secretary Susanne Schraut for perfect help with handling submission and review process for the manuscripts contributed to this special issue.
References (27)
Mathematical aspects of dose–response studies in carcinogenesis—the concept of thresholds
Oncology
(1976)Efficiency of the adaptive response of Escherichia coli to alkylating agents
Nature
(1980)- et al.
Hormesis: a generalizable and unifying hypothesis
Crit. Rev. Toxicol.
(2001) - et al.
Population variability in biological adaptive responses to DNA damage and the shapes of carcinogen dose–response curves
Toxicol. Appl. Pharmacol.
(2005) - et al.
Nonmonotonic dose–response relationships: mechanistic basis, kinetic modeling, and implications for risk assessment
Toxicol. Sci.
(2004) - et al.
Mechanistic influences for mutation induction curves after exposure to DNA-reactive carcinogens
Cancer Res.
(2007) - et al.
What happened: the chemistry side of the incident with EMS contamination in Viracept tablets
Toxicol. Lett.
(2009) - et al.
MNT and Muta™mouse studies to define the in vivo dose response relations of the genotoxicity of EMS and ENU
Toxicol. Lett.
(2009) - et al.
Literature review on the genotoxicity, reproductive toxicity, and carcinogenicity of ethyl methanesulfonate
Toxicol. Lett.
(2009) - et al.
EMS in Viracept—initial (‘traditional’) assessment of risk to patients based on linear dose response relations
Toxicol. Lett.
(2009)
In vivo genotoxicity of EMS: statistical assessment of the dose response curves
Toxicol. Lett.
MGMT: key node in the battle against genotoxicity, carcinogenicity and apoptosis induced by alkylating agents
DNA Repair (Amst.)
In vivo and in vitro characterization of ethyl methanesulfonate pharmacokinetics in animals and in human
Toxicol. Lett.
Cited by (18)
Assessment of mechanisms driving non-linear dose-response relationships in genotoxicity testing
2015, Mutation Research - Reviews in Mutation ResearchCitation Excerpt :In the majority of the reviewed literature, single effect levels, i.e. NO[G]EL or LO[G]EL are experimentally derived and have been used to describe non-linear MoA. Only recently, more emphasis has been given to mathematical modeling and dose–response curves by using different approaches (e.g. hockeystick model, BMD; reviewed by [14,125,193,194]). Modeling of data more precisely describes the shape of the dose–response curve and can be used to extrapolate to lower dose levels than experimentally assessed.
Interpretation of the margin of exposure for genotoxic carcinogens - Elicitation of expert knowledge about the form of the dose response curve at human relevant exposures
2013, Food and Chemical ToxicologyCitation Excerpt :The majority view expressed during the workshop was that, in reality, both mechanisms (i.e. (a) the cancer process is stochastic and (b) the cancer process is a reflection of underlying inter-individual susceptibility) are likely to be in operation, but that there is little or no information on the balance between them. Consequently, it would be pragmatic to start by assuming linearity and consider various models for departures from linearity depending on the particular chemical concerned, e.g. saturation of clearance, saturation of metabolic activation, induction of repair, saturation of repair, induction of cytotoxicity, cell division slow down, age effects, etc. (Gocke and Müller, 2009; Lutz, 2009). However, although experts may have considered a number of different philosophical perspectives, they appeared to agree that species differences, sensitive subpopulation/s, and most particularly the mode of action (and consequently extrapolation between chemicals) are amongst the most significant factors influencing the true form of the dose–response curve for genotoxic carcinogens.
The assessment of impurities for genotoxic potential and subsequent control in drug substance and drug product
2013, Journal of Pharmaceutical SciencesMutagenicity of N-ethyl-N-nitrosourea, N-methyl-N-nitrosourea, methyl methanesulfonate and ethyl methanesulfonate in the developing Syrian hamster fetus
2010, Mutation Research - Genetic Toxicology and Environmental MutagenesisEmpagliflozin: Validation of Stability-Indicating LC Method and in silico Toxicity Studies
2024, Journal of Chromatographic Science
- ☆
For the Special Issue of Toxicology Letters entitled “Assessment of human toxicological risk of Viracept patients accidentally exposed to ethyl methanesulfonate (EMS) based on preclinical investigations with EMS and ethylnitrosourea”.