Toxicology of halogenated aliphatic hydrocarbons: structural and molecular determinants for the disturbance of chromosome segregation and the induction of lipid peroxidation

Ecological and human health risk assessment of total petroleum hydrocarbons in surface water and sediment from Woji Creek in the Niger Delta Estuary of Rivers State, Nigeria

This study was designed to assess total petroleum hydrocarbon (THP) concentrations in the surface water and sediment sampled from Woji Creek and to assess potential ecological and human health risk due to petroleum hydrocarbons along the creek. Physicochemical parameters [pH, temperature (T), electrical conductivity (EC), dissolved oxygen (DO), total dissolved solids (TDS)] were in-situ measured from sediment and surface water; hydrological parameters (width, depth and volume) were used to calculate the flow rate (discharge) at different stations of the creek. Trend of TPH in the surface water samples along the creek were as follows: St4 (3.639 ± 1.121 mg/L) > St3 (2.449 ± 0.623 mg/L) > St1 (1.457 ± 0.244 mg/L) > St2 (1.069 ± 0.228 mg/L) > St5 (1.010 ± 0.120 mg/L) Trend of TPH concentration across the creek was as follows: St1 - 8.758 ± 0.697 mg/kg > St3 - 7.675 ± 0.541 mg/kg > St5 - 5.515 ± 0.401 mg/kg > St4 - 5.075 ± 0.363 mg/kg > St2 - 3.162 ± 0.307 mg/kg. Diagnostic indices indicate that the hydrocarbon in the creek was from petrogenic source. Estimation of ecological risk indicated risk in the surface water but not in the sediment. However, human health risk assessment indicated no risk due to human ingestion of the sediment or surface water.

A descriptor-based analysis to highlight the mechanistic rationale of mutagenicity

Cancer is a main concern for human health and there is a need of alternative methodologies to rapidly screen large quantitative of compounds that may represent a toxicological risk. Here a statistical analyses is performed on a benchmark database of experimental Ames data to identify chemical descriptors discriminating mutagens and non-mutagens. A total of 53 activating and deactivating modulators are identified, that flagged respectively a percentage of mutagen and non-mutagen up to 87%. Modulators are further combined to form synergistic cross-terms, accounting for the effect that combined properties may have on the final toxicity. Exclusion rules are defined as exception to the modulators. Synergistic cross-terms and exclusion rules improve the enrichment of mutagens/non-mutagens with respect of the original abundance in the dataset to values higher than 95%. The external predictivity of modulators and cross-terms reach balanced accuracy up to 0.775 that is analogous to other mutagenicity models from the literature, confirming the suitability of the rules to real-life screening of chemicals. Modulators are discussed for their mechanistic link to mutagenicity. This analysis confirms the key role of some properties (polarizability, shape, mass, presence of reactive functional groups or unsaturated planar systems) as driving elements for the initiation of the mutagenicity.

Exacerbation of substrate toxicity by IPTG in Escherichia coli BL21(DE3) carrying a synthetic metabolic pathway

BACKGROUND

Heterologous expression systems based on promoters inducible with isopropyl-β-D-1-thiogalactopyranoside (IPTG), e.g., Escherichia coli BL21(DE3) and cognate LacI(Q)/P(lacUV5)-T7 vectors, are commonly used for production of recombinant proteins and metabolic pathways. The applicability of such cell factories is limited by the complex physiological burden imposed by overexpression of the exogenous genes during a bioprocess. This burden originates from a combination of stresses that may include competition for the expression machinery, side-reactions due to the activity of the recombinant proteins, or the toxicity of their substrates, products and intermediates. However, the physiological impact of IPTG-induced conditional expression on the recombinant host under such harsh conditions is often overlooked.

RESULTS

The physiological responses to IPTG of the E. coli BL21(DE3) strain and three different recombinants carrying a synthetic metabolic pathway for biodegradation of the toxic anthropogenic pollutant 1,2,3-trichloropropane (TCP) were investigated using plating, flow cytometry, and electron microscopy. Collected data revealed unexpected negative synergistic effect of inducer of the expression system and toxic substrate resulting in pronounced physiological stress. Replacing IPTG with the natural sugar effector lactose greatly reduced such stress, demonstrating that the effect was due to the original inducer's chemical properties.

CONCLUSIONS

IPTG is not an innocuous inducer; instead, it exacerbates the toxicity of haloalkane substrate and causes appreciable damage to the E. coli BL21(DE3) host, which is already bearing a metabolic burden due to its content of plasmids carrying the genes of the synthetic metabolic pathway. The concentration of IPTG can be effectively tuned to mitigate this negative effect. Importantly, we show that induction with lactose, the natural inducer of P lac , dramatically lightens the burden without reducing the efficiency of the synthetic TCP degradation pathway. This suggests that lactose may be a better inducer than IPTG for the expression of heterologous pathways in E. coli BL21(DE3).

Evaluating genotoxicity data to identify a mode of action and its application in estimating cancer risk at low doses: A case study involving carbon tetrachloride

In the new USEPA cancer risk assessment guidelines, mode of action (MoA) information, combined with a determination of whether or not a chemical is mutagenic, plays an important role in determining whether a linear or nonlinear approach should be used to estimate cancer risks at low doses. In this article, carbon tetrachloride (CT) is used as an example to illustrate how mixed genotoxicity data can be evaluated and used to identify a likely MoA. CT is essentially negative in inducing gene mutations in Salmonella, but is consistently positive in inducing recombination and aneuploidy in fungi. Negative or equivocal results were seen in most in vitro and in vivo studies in mammals, including mutation studies in transgenic mice. However, DNA adducts, primarily those derived from oxidation- and lipid-peroxidation-derived products as well as DNA double-strand breaks and micronucleated cells, have been seen repeatedly in the liver of CT-treated animals. On the basis of the weight of evidence, CT should not be considered a directly mutagenic agent. Mutagenic as well as other genotoxic effects, as they occur, will most likely be generated through indirect mechanisms resulting from oxidative and lipid peroxidative damage and/or damage occurring during necrosis or apoptosis. As key events in this process are expected to occur in a nonlinear fashion, the expected relationship between CT dose and carcinogenic response in the liver is likely to be nonlinear with a steep dose response. This conclusion is consistent with rodent cancer bioassay results in which steep nonlinear dose responses have been seen.

Quantitative structure-activity relationships for toxicity and genotoxicity of halogenated aliphatic compounds: wing spot test of Drosophila melanogaster

Halogenated aliphatic compounds were evaluated for toxic and genotoxic effects in the somatic mutation and recombination test employing Drosophila melanogaster. The tested chemicals included chlorinated, brominated and iodinated; mono-, di- and tri-substituted; saturated and unsaturated alkanes: 1,2-dibromoethane, 1-bromo-2-chloroethane, 1-iodopropane, 2,3-dichloropropene, 3-bromo-1-propene, epibromohydrin, 2-iodobutane, 3-chloro-2-methylpropene, 1,2,3-trichloropropane, 1,2-dichloroethane, 1,2-dichlorobutane, 1-chloro-2-methylpropane, 1,3-dichloropropane, 1,2-dichloropropane, 2-chloroethymethylether, 1-bromo-2-methylpropane and 1-chloropentane. N-methyl-N-nitrosourea served as the positive and distilled water as the negative control. The set of chemicals for the toxicological testing was selected by the use of statistical experiment design. Group of unsaturated aliphatic hydrocarbons were generally more toxic than saturated analogues. The genotoxic effect was observed with 14 compounds in the wing spot test, while 3 substances did not show any genotoxicity by using the wing spot test at 50% lethal concentration. The highest number of wing spots was observed in genotoxicity assay with 1-bromo-2-chloroethane, 1,2-dichloroethane, 1,2-dibromoethane and 1-iodopropane. Nucleophilic superdelocalizability calculated by quantum mechanics appears to be a good parameter for prediction of both toxicity and genotoxicity effects of halogenated aliphatic compounds.

Sub-chronic exposure to 1,1-dichloropropene induces frameshift mutations in lambda transgenic medaka

1,1-Dichloropropene (1,1-DCP) is a contaminant present in both ground and surface waters used as sources for drinking water. Structural similarity to several compounds with known mutagenicity and carcinogenicity, and recent demonstration of mutagenicity in vitro, suggest this compound may be similarly mutagenic in vivo. A transgenic fish model, the lamda transgenic medaka, was used to evaluate the potential mutagenicity of this contaminant in vivo following sub-chronic exposure for 6 weeks. Mutant frequencies of the cII target gene (MF) increased six-fold in the livers of fish exposed to the lowest 1,1-DCP exposure concentration (0.44 mg/L, MF = 18.4 x 10(-5), and increased with each treatment, culminating in a 32-fold induction in fish from the highest 1,1-DCP treatment (16.60 mg/L, MF = 96.3 x 10(-5). Mutations recovered from treated fish showed a distinctive mutational spectrum comprised predominantly of +1 frameshift mutations, induced 166-fold above that of untreated animals. The majority of frameshifts were +1 insertions at thiamine and adenine. These results represent the first evidence of mutagenicity of 1,1-DCP in vivo, and of the highly characteristic spectrum of induced mutations dominated by +1 frameshift mutations. Based upon results from previous in vitro studies, the similar role of glutathione S-transferase (GSTT1-1) in the activation of 1,1-DCP to a mutagen in vivo is also suggested. This study further illustrates the utility of the lamda transgenic medaka as a model for identifying and characterizing potential genetic health risks associated with chemical exposures in the environment.

Proteomic Maps−Toxicity Relationship of Halocarbons Studied with Similarity Index and Genetic Algorithm

In this work we analyzed proteomic maps obtained from hepatocytes, which were treated with 14 halocarbons. A similarity index was introduced as a robust measure of similarity between two maps or between two selections of spots within the maps. A searching algorithm was used to identify the spots that may play an important role in toxicity mechanism. The highest correlation coefficients obtained between the similarity index and biological parameter were larger than 0.9.

Combining Chemodescriptors and Biodescriptors in Quantitative Structure−Activity Relationship Modeling

In view of the wide distribution of halocarbons in our world, their toxicity is a public health concern. Previous work has shown that various measures of toxicity can be predicted with standard molecular descriptors. In our work, biodescriptors of the effect of halocarbons on the liver were obtained by exposing hepatocytes to 14 halocarbons and a control and by producing two-dimensional electrophoresis gels to assess the expressed proteome. The resulting spot abundances provide additional biological information that might be used in toxicity prediction. QSAR models were fitted via ridge regression to predict eight dependent toxicity measures: d37, arr, EC50MTT, EC50LDH, EC20SH, LECLP, LECROS, and LECCAT. Three predictor sets were used for each-the chemodescriptors alone, the biodescriptors alone, and the combined set of both chemo- and biodescriptors. The results differed somewhat from one dependent to another, but overall it was shown that better results could be obtained by using both chemo- and biodescriptors in the model than by using either chemo- or biodescriptors alone. The library of compounds used was small and quite homogeneous, so our immediate conclusions are correspondingly limited in scope, but we believe the underlying methodologies have broad applicability at the interface of chemical and biological descriptors.

Simple and alpha,beta-unsaturated aldehydes: correct prediction of genotoxic activity through structure-activity relationship models

Aldehydes are widespread environmental and industrial compounds, able to stimulate a range of adverse health effects (e.g., general toxicity, allergenic reactions, mutagenicity, and carcinogenicity). We have previously presented quantitative structure-activity relationships (QSARs) for the genotoxicity of simple and alpha,beta-unsaturated aliphatic aldehydes. In this study, we show that the QSAR models are able to correctly predict--based only on the knowledge of the chemical structure--the genotoxicity of other aldehydes, not considered in the development phase of the models. This adds confidence to the reliability of our QSAR models as tools for the theoretical assessment of the genotoxic hazard posed by aldehydes. The analysis of SOS Chromotest induction ability and the ease of formation of DNA adducts by the aldehydes provided further mechanistic insights.

Prediction of cellular toxicity of halocarbons from computed chemodescriptors: a hierarchical QSAR approach

A hierarchical quantitative structure-activity relationship (HiQSAR) approach was used to estimate toxicity and genetic toxicity for a set of 55 halocarbons using computed chemodescriptors. The descriptors consisted of topostructural (TS), topochemical (TC), geometrical, semiempirical (AM1) quantum chemical, and ab initio (STO-3G, 6-31G(d), 6-311G, 6-311G(d), and aug-cc-pVTZ) quantum chemical indices. For the two toxicity endpoints investigated, ARR and D(37), the TC indices gave the best cross-validated R(2) values. The 3-D indices also performed either as well as or slightly superior to the TC indices. For the four categories of quantum chemical indices used for the development of predictive models, the AM1 parameters gave the worst performance, and the most advanced ab initio (B3LYP/aug-CC-pVTZ) parameters gave the best results when used alone. This was also the case when the quantum chemical indices were used in the hierarchical QSAR approach for both of the toxicity endpoints, ARR and D(37). The models resulting from HiQSAR are of sufficiently good quality to estimate toxicity of halocarbons from structure.

Quantum descriptors for predictive toxicology of halogenated aliphatic hydrocarbons

In order to improve Quantitative Structure-Activity Relationships (QSARs) for halogenated aliphatics (HA) and to better understand the biophysical mechanism of toxic response to these ubiquitous chemicals, we employ improved quantum-mechanical descriptors to account for HA electrophilicity. We demonstrate that, unlike the lowest unoccupied molecular orbital energy, ELUMO, which was previously used as a descriptor, the electron affinity can be systematically improved by application of higher levels of theory. We also show that employing the reciprocal of ELUMO, which is more consistent with frontier molecular orbital (FMO) theory, improves the correlations with in vitro toxicity data. We offer explanations based on FMO theory for a result from our previous work, in which the LUMO energies of HA anions correlated surprisingly well with in vitro toxicity data. Additional descriptors are also suggested and interpreted in terms of the accepted biophysical mechanism of toxic response to HAs and new QSARs are derived for various chemical categories that compose the data set employed. These alternate descriptors provide important insight and could benefit other classes of compounds where the biophysical mechanism of toxic response involves dissociative attachment.

The importance of hydrophobicity and electrophilicity descriptors in mechanistically-based QSARs for toxicological endpoints

Quantitative structure-activity relationship (QSAR) analysis of four toxicological data sets is described. The toxicological data include three data sets retrieved from the literature (the toxic and metabolic effects of 23 aliphatic alcohols on the perfused rat liver; the toxicity of 21 pyridines to mice; the lethality of 55 halogenated hydrocarbons to the mould Aspergillus nidulans). In addition, the toxicity of 13 mono- and di-substituted nitrobenzenes in a 15 min assay using the alga Chlorella vulgaris was analysed. QSARs were developed successfully using descriptors to describe uptake in the organism (i.e. hydrophobicity as quantified by the logarithm of the octanol-water partition coefficient, log P) and reactivity at the site of action (i.e. electrophilicity as quantified by the energy of the lowest unoccupied molecular orbital, E(LUMO)). A further parameter describing molecular branching as also required to model the data for the aliphatic alcohols. The results demonstrate that mechanistically based QSARs can be developed for these diverse endpoints which are, in terms of statistical quality as good as, if not better, than QSARs based on less mechanistically interpretable descriptors.

QSAR modeling of oxidative stress in vitro following hepatocyte exposures to halogenated methanes

Volatile halogenated aliphatic compounds are among those chemicals that can cause oxidative stress in vitro and in vivo. Relationships can be identified between the potential of these chemicals to elicit certain biological responses and their specific chemical descriptors, such as molecular orbital energies (LUMO) or partition coefficients (logP). A quantitative structure-activity relationship (QSAR) model has not been reported previously for the potential of a series of brominated and chlorinated methanes to induce oxidative stress in primary rat hepatocytes. By utilizing a novel in vitro methodology to expose cultures of rat primary hepatocytes to volatile chemicals, biological responses were assessed from exposures of hepatocytes to individual halogenated methanes. Indicators of lipid peroxidation, reactive oxygen species and cytotoxicity were measured. For the 10 brominated and chlorinated methanes tested, semi-empirical molecular orbital methods were used to calculate the physical/chemical descriptors used in the QSAR models. These models were used to explain the relative potential for a given halogenated methane to induce markers of oxidative stress or related damage in vitro. The results showed that certain descriptors, such as the molecular orbital energies, bond lengths, and lipophilicity are quantitatively correlated with induction of indicators for oxidative stress and cytotoxicity by halogenated methanes in primary rat hepatocytes.

Structures and electron attachment properties of halomethanes (CX(n)Y(m), X=H, F; Y=Cl, Br, I; n=0,4; m=4-n)

Theoretical studies of structures of neutral molecules and their anions as well as dissociative electron attachment properties are presented for the halomethanes of general formula CX(n) Y(m); X=H, F; Y=Cl, Br, I; n=0,4; m=4-n. The dissociative electron attachment seems to be the primary process resulting in toxicity of these species. The halomethane anions containing hydrogens are formed as radical-anion adducts. When H is replaced by F, these species become true sigma* radicals. The electron affinities are computed using a variety of computational techniques including the four-order Møller-Plesset (MP4) technique that included 250 basis functions. It is challenging to compare the computed results with experiment due to dearth of experimental data and uncertainties in the existing experimental data. Thus in certain cases larger differences are found between the computed and experimental values.

Improved QSARs for predictive toxicology of halogenated hydrocarbons

In our continuing efforts to provide a predictive toxicology capability, we seek to improve QSARs (quantitative structure-activity relationships) for chemicals of interest. Currently, although semi-empirical molecular orbital methods are hardly the state of the art for studying small molecules, AM1 calculations appear to be the method of choice when calculating quantum-chemical descriptors. However, with the advent of modern computational capabilities and the development of fast algorithms, ab initio molecular orbital and first principles density functional methods can be expeditiously applied in current QSAR studies. We present a study on halogenated alkanes to assess whether more accurate quantum methods result in QSARs that correlate better with experimental data. Furthermore, improved QSARs can also be obtained through development of new descriptors with explicit physical interpretations that should lead to better understanding of the mechanisms involved in the toxic response. We show that descriptors calculated from chemical intermediates may be useful in future QSARs.

Structural and mechanistic bases for the induction of mitotic chromosomal loss and duplication ('malsegregation') in the yeast Saccharomyces cerevisiae: relevance to human carcinogenesis and developmental toxicology

MultiCASE has the ability to automatically determine the structural features responsible for the biological activity of chemicals. In the present study, 93 chemicals tested for their ability to induce chromosomal 'malsegregation' in the yeast Saccharomyces cerevisiae were analyzed. This 'malsegregation' mimics molecular events that occur during human development and carcinogenesis resulting in an effective loss of one chromosome of an autosomal pair and duplication of the homologue. Structural features associated with the ability to induce such chromosome loss and duplication were identified and compared with those obtained from examination of other toxicological data bases. The most significant structural similarities were identified between the induction of chromosomal malsegregation and several toxicological phenomena such as cellular toxicity, induction of sister chromatid exchanges in vitro and rodent developmental toxicity. Very significant structural similarities were also found with systemic toxicity, induction of micronuclei in vivo and human developmental toxicity. Less significant structural overlaps were found between yeast malsegregation and rodent carcinogenicity, DNA reactivity and mutagenicity, and the induction of chromosome aberrations in vitro and sister chromatid exchanges in vivo. These overlaps may indicate mechanistic similarities between the induction of chromosomal malsegregation and other toxicological phenomena. The predictivity of the SAR model derived from the present data base is relatively low, however. This may be merely a reflection of the small size and composition of the data base, however, further analyses suggest that it reflects primarily the multiple mechanisms responsible for the induction of chromosomal malsegregation in yeast and the complexity of the phenomenon.