Protocol
Abstract
Background: Chemicals are not required to be tested systematically for their neurotoxic potency, although they may contribute to the development of several neurological diseases. The absence of systematic testing may be partially explained by the current Organisation for Economic Co-operation and Development (OECD) Test Guidelines, which rely on animal experiments that are expensive, laborious, and ethically debatable. Therefore, it is important to understand the risks to exposed workers and the general population exposed to domestic products. In this study, we propose a strategy to test the neurotoxicity of solvents using the commonly used glycol ethers as a case study.
Objective: This study aims to provide a strategy that can be used by regulatory agencies and industries to rank solvents according to their neurotoxicity and demonstrate the use of toxicokinetic modeling to predict air concentrations of solvents that are below the no observed adverse effect concentrations (NOAECs) for human neurotoxicity determined in in vitro assays.
Methods: The proposed strategy focuses on a complex 3D in vitro brain model (BrainSpheres) derived from human-induced pluripotent stem cells (hiPSCs). This model is accompanied by in vivo, in vitro, and in silico models for the blood-brain barrier (BBB) and in vitro models for liver metabolism. The data are integrated into a toxicokinetic model. Internal concentrations predicted using this toxicokinetic model are compared with the results from in vivo human-controlled exposure experiments for model validation. The toxicokinetic model is then used in reverse dosimetry to predict air concentrations, leading to brain concentrations lower than the NOAECs determined in the hiPSC-derived 3D brain model. These predictions will contribute to the protection of exposed workers and the general population with domestic exposures.
Results: The Swiss Centre for Applied Human Toxicology funded the project, commencing in January 2021. The Human Ethics Committee approval was obtained on November 16, 2022. Zebrafish experiments and in vitro methods started in February 2021, whereas recruitment of human volunteers started in 2022 after the COVID-19 pandemic–related restrictions were lifted. We anticipate that we will be able to provide a neurotoxicity testing strategy by 2026 and predicted air concentrations for 6 commonly used propylene glycol ethers based on toxicokinetic models incorporating liver metabolism, BBB leakage parameters, and brain toxicity.
Conclusions: This study will be of great interest to regulatory agencies and chemical industries needing and seeking novel solutions to develop human chemical risk assessments. It will contribute to protecting human health from the deleterious effects of environmental chemicals.
International Registered Report Identifier (IRRID): DERR1-10.2196/50300
doi:10.2196/50300
Keywords
Introduction
Environmental and occupational exposure to chemicals may contribute to the development of several neurological diseases [
, ]. In particular, organic solvents used in industries such as car repair, painting, furniture manufacturing, printing, and cleaning have been associated with several central nervous system (CNS) conditions. These include mild to severe toxic encephalopathy [ ]; deficits in cognitive function [ - ]; and, in some cases, neurodegenerative diseases [ , ]. The diffuse neuropathological effects of acute solvent intoxication reflect neurophysiological abnormalities involving multiple brain regions. With increasingly intense or prolonged exposure, the severity of acute impairment may progress along the spectrum of delirium. Chronic high-level exposure may lead to global cognitive impairment including deficits in memory, attention, energy, and personality, which are well-described forms of dementia [ - ]. Much of the initial work on organic solvent toxicity originated in Scandinavia, where a neurobehavioral syndrome in painters leading to their early retirement was first described [ ]. However, although the neurotoxicity of solvents such as toluene, trichloroethylene, and n-hexane has been recognized, the neurotoxicity of common solvents currently on the market has not been evaluated. Notably, neurotoxicity testing is only required if the chemical is deemed a pesticide; otherwise, all other chemicals are evaluated on a case-by-case basis. The only exception is if the compound structure is suspected to have nervous system targets and no data are available for read-across or when effects on the nervous system are found in single-dose (Organisation for Economic Co-operation and Development [OECD] Test Guidelines [TGs] 402, 403, 420, 423, or 425) or repeated-dose toxicity studies (TG 407 or 408). Because the nervous system effect endpoints considered as triggers (ie, modifications of wet brain weight or basic histopathology, or both) are quite insensitive, high amounts of potentially neurotoxic compounds are available on the market. The European Union Classification, Labelling, and Packaging Regulation does not include a classification for neurotoxicity. Exposure to organic solvents, especially among workers with higher exposure than the general population, can produce neurotoxic effects, depending on the internal dose. Therefore, to efficiently protect the population from possible solvent toxicity, it is important to determine the air concentrations at which neurotoxicity does not occur. To this end, we propose a strategy using a combination of in vivo (zebrafish embryo), in vitro, and in silico tools coupled with controlled in vivo human exposure experiments to assess the neurotoxicity of solvents ( ).The recognized method for the evaluation of the neurotoxic potential of chemicals, the OECD TG 424 (neurotoxicity in rodents), uses complex in vivo tests in rodents, which are laborious, expensive, difficult to apply in a standardized manner, and ethically debatable. Regulators from different agencies worldwide as well as the scientific community are becoming increasingly aware of the limitations of the current toxicity testing paradigm. Animal-based high-dose testing in typically 1 stand-alone guideline test is not always relevant for human exposure scenarios [
]. One of the most challenging aspects of this animal-centric approach is the impossibility of coping with the thousands of chemicals for which data are still lacking. Conducting animal tests is time consuming and expensive. Therefore, they cannot be carried out routinely because of the sheer number of chemicals that are currently on the market and those anticipated to enter it in the coming years [ ]. In addition, there are shortcomings regarding interspecies concordance between different mammalian or rodent species as well as with respect to extrapolation from experimental animals to humans. These ambiguities in results or poor reproducibility performance call into question the relevance of such test methods for human risk assessment [ - ]. All this prompts a move away from animal testing toward a combination of in vitro and in silico approaches that address functional mechanistic endpoints [ , , ].Numerous in vitro models have been proposed for the evaluation of neurotoxicity in the last decades. Monolayer cultures of a single brain cell type are far from representing the human brain in terms of architecture and functionality. Given the sophistication of brain cell-to-cell interactions, some complexity is required to recapitulate human-relevant cellular processes and functions in vitro. However, a good balance must be found between this complexity and the simplicity needed to have robust and reproducible systems that can be applied for chemical screening in a high-throughput manner [
]. The 3D hiPSC–derived brain test system (BrainSpheres) we previously developed [ ] will be used in this study, as it fulfills these requirements.The blood-brain barrier (BBB) protects the brain parenchymal cells from the deleterious effects of xenobiotics. However, some chemicals are able to cross or impair the BBB [
]. The transient or permanent opening of the BBB provides xenobiotics, plasma proteins, and immunoregulatory mediators access to the CNS, where they can induce toxic effects. Therefore, we will implement a predictive model to assess the impact of solvents on BBB based on in vivo (zebrafish embryo), in vitro (human brain microvascular endothelial cells [hCMEC/D3]), and in silico models [ ] to assist in the interpretation of the results obtained in in vitro neurotoxicity testing.Glycol ethers will be used as a case study to evaluate the feasibility of our protocol. Glycol ethers form a wide family of a few dozen solvents with different physicochemical properties making them versatile and usable in a variety of industrial applications ranging from pharmaceuticals and microelectronics to domestic cleaning, personal care, and printing. The 2 main groups of glycol ethers are the E series (ie, ethylene glycol ethers [EGEs]) and the P series (ie, propylene glycol ethers [PGEs]). EGEs and PGEs show differences in their toxicological properties regarding teratogenicity, hemolysis, and testicular atrophy [
], apparently resulting from their distinct production of metabolites [ , ]. EGEs have a primary alcohol group and are oxidized by alcohol dehydrogenase and aldehyde dehydrogenase to form the toxic alkoxyacetic acid. Therefore, PGEs are progressively introduced as a less toxic replacement of EGEs. PGEs are sold as a mixture of 2 isomers, with the bulk having a secondary alcohol (a) group and generally <5% of primary alcohol (b) groups [ ]. The b-isomer is oxidized by alcohol dehydrogenase and aldehyde dehydrogenase to form the toxic alkoxypropionic acid in the body. However, the actual toxicity of PGEs is poorly characterized. Therefore, it is essential for our strategy to study not only the parent compound but also the metabolites. Metabolically competent liver cell models will be used to screen for liver toxicity of solvents and for the production of potential metabolites.Exposure to solvents in human volunteers is fundamental when quantifying possible risks from chemicals because toxicological effects are related to internal exposure, that is, the concentration of a chemical inside the body and its biotransformation. Therefore, in vivo human volunteer studies are necessary to quantify absorption, distribution, metabolism, and excretion (ADME) kinetics in humans under controlled exposure conditions, which are necessary to understand the relationship between solvent concentrations in the air and biological samples, such as blood, urine, or exhaled air. These experiments provide important data on the total absorbed dose after inhaled solvent concentration, absorption rate (ie, from the site of absorption into the bloodstream), biotransformation rate (ie, metabolization of the parent chemical with production of metabolites), and elimination rate (ie, excretion of the parent chemical and metabolites from the body) [
]. These toxicokinetic parameters are representative of the target organ or of the tissue concentrations that may trigger an effect and are, therefore, relevant for understanding chemical risks in humans.Toxicokinetic models quantitatively describe the body’s ADME of a chemical or substance (different terms for this concept are preferred in different fields, including “toxicokinetics,” “pharmacokinetics,” and “biokinetics”) [
]. Furthermore, using a toxicokinetic model in reverse dosimetry, we can predict the solvent air concentrations leading to brain concentrations below the levels found to produce neurotoxic effects in vitro in the BrainSpheres. The toxicokinetic model will incorporate metabolism parameters derived from the in vitro liver system and passage through or toxicity to BBB. Once calibrated, the toxicokinetic model can be used to simulate chronic exposure scenarios to predict cumulative brain concentrations and used in reverse dosimetry to predict air concentrations that will not likely result in brain concentrations associated with toxicity.Methods
Ethical Considerations
Ethics committee approval was obtained from Swiss ethics (Commission cantonale d’éthique de la recherche sur l’être humain) in 2022 for this nonclinical human study (2022-01567). Healthy women and men were recruited as participants in our study. Each participant signed a written informed consent form before inclusion in the study. The participants will be reimbursed for their time and inconvenience according to the Swiss guidelines.
Global Strategy
The choice of solvents to be included in the study will be based on the amount annually placed on the European market and the number of products registered containing known glycol ethers. The selected organic solvents will be applied to various in vitro models to determine their neurotoxicity. They must be amphiphilic to be solubilized in the cell culture media. We established the following solvent selection criteria: (1) used or produced >1 metric ton per year, (2) incorporated in numerous industrial and commercial products, and (3) water solubility. This selection process involves consulting government databases and contacting different industry sectors. We will start the project concomitantly by testing 2 solvents of the P series, propylene glycol methyl ether (PGME), for which we have already developed a toxicokinetic model, and propylene glycol butyl ether. We will test 1 additional solvent from the E series with the in vitro test system, namely, ethylene glycol methyl ether (EGME), which has been banned for use in cosmetic products in Europe [
].The study is organized into 5 work packages (WPs). The information workflow between the WPs is shown in
. All results collected from the abovementioned systems will contribute to refining the toxicokinetic model we previously developed for PGME [ ]. The toxicokinetic parameters of the solvent and the metabolites will then be characterized in human volunteers after exposure to PGE vapors under controlled conditions. These results will be used to calibrate and expand our toxicokinetic model [ ]. The toxicokinetic model will be constructed to predict brain concentrations of selected solvents and, consequently, will include a brain compartment to predict the target organ solvent and metabolite concentrations. BBB parameters such as barrier transport, transport of the solvent once in the brain, and solvent-brain binding will be incorporated. Solvent neurotoxicity may depend on the metabolic modifications of the substances; therefore, we will incorporate the parameters for the parent compound and the metabolites found in the in vitro liver system. The data necessary to build the model will be retrieved from peer-reviewed scientific literature for tissue:blood partition coefficients (PCs) [ ] following the fit-for-purpose dose-response analysis approach. The toxicokinetic model should be able to predict human brain concentration for each of the tested solvents after inhalation exposure, given the air concentration of vapors and duration of exposure. The simulated human brain concentrations will then be compared with the no observed adverse effect concentrations (NOAECs) obtained from the neurotoxicity in vitro system (BrainSpheres). In addition, the toxicokinetic model will be used to predict solvent air concentrations that are unlikely to lead to brain concentration equal to or superior to the brain NOAECs using reverse dosimetry. The specific aims of each WP are summarized in .WP | Name | Specific aims |
WP1 | In vitro neurotoxicity testing |
|
WP2 | In vivo or in vitro or in silico BBBc functionality testing |
|
WP3 | In vitro hepatic metabolism and clearance |
|
WP4 | In vivo volunteer exposure |
|
WP5 | In silico PBTKg modeling |
|
aNOAEC: no observed adverse effect concentration.
bKE: key event.
cBBB: blood-brain barrier.
dhCMEC/D3: human brain microvascular endothelial cells.
eVmax: maximum velocity.
fKm: Michaelis constant.
gPBTK: physiologically based toxicokinetic.
WP1: In Vitro Neurotoxicity Testing
Testing strategies are needed to evaluate the neurotoxicity of chemicals in a more cost-effective, efficient, and ethical manner. Participating in an international effort, we developed a 3D human-induced pluripotent stem cells (hiPSC)–derived brain model containing several subtypes of neurons, astrocytes, and oligodendrocytes [
]. This system allows the cells to reach a high level of differentiation and cellular maturation, exemplified by the presence of functional synapses and compact myelin. The presence of myelin is important for this project because solvents more easily target lipid-rich structures [ ]. This 3D human brain model has already proven its usefulness for neurotoxicity testing [ - ]. We hypothesized that glycol ethers are neurotoxic. Therefore, we propose to take advantage of our hiPSC-derived BrainSpheres model to study the neurotoxicity induced by uncharacterized glycol ethers present on the market, which will be compared with the neurotoxicity of well-characterized solvents known to induce human encephalopathy.Solvents are data-poor substances. It was originally hypothesized that they exert their toxic effects largely through nonspecific physicochemical effects that modulate membrane fluidity and perturb the hydrophobic force regulating macromolecular interactions [
]. However, recent evidence supports the view that solvents interact with lipophilic areas on protein receptors [ , ]. They have also been shown to induce lipid peroxidation, leading to mitochondrial dysfunction, failure of electron transport, and energy production [ , ]. In this study, omics (eg, proteomics, metabolomics, and lipidomics) technology will be used to decipher the mechanisms of glycol ether neurotoxicity and to identify potential biomarkers of toxic effects.Primary 3D hiPSC-derived brain cell cultures will be prepared and maintained as previously described [
]. This model contains neurons that form synapses, astrocytes, and oligodendrocytes myelinating the axons ( ). Cytotoxicity will be determined by a resazurin assay after repeated exposure (7 d) to the selected glycol ethers (parent compounds and metabolites). Gene expression for cell type–specific genes, markers of synapses and myelin, and markers of cell stress will be quantified by quantitative reverse transcription polymerase chain reaction at concentrations of solvents under half-maximal effective concentration (EC50) for cytotoxicity. Immunostaining will be performed to assess the effects of solvents on synapses, myelin, and astrocyte reaction, and immunofluorescence will be quantified. NOAECs ( ) will be determined for all tested endpoints, as previously shown for gene expression [ ]. Brain cell cultures will also be exposed to the metabolites of PGME, propylene glycol butyl ether, and EGME and to the metabolites of the newly selected uncharacterized solvents, potentially produced by liver metabolism (WP3).To establish the in vitro distribution kinetics of selected solvents necessary for toxicokinetic modeling, 3D brain cell cultures and medium will be collected 3, 6, 24, and 48 hours after the first exposure and after the last exposure of the repeated treatment. The solvent and its main metabolites (if relevant) will be quantified to establish a time course of disappearance from the medium and appearance in the cells as well as to assess the potential accumulation for the entire period of exposure. The fraction bound to culture plates’ plastic will be quantified after desorption. In silico modeling of glycol ethers in vitro distribution kinetics will then be developed. This model will be able to predict the change in cell-associated concentrations of solvents in BrainSpheres with time, as previously shown for amiodarone [
].WP2: In Vivo, In Vitro, and In Silico BBB Functionality Testing
We previously established the zebrafish as a predictive vertebrate screening model to study the systemic circulation and tissue distribution of particulate drug carriers [
, ]. At 72 hours postfertilization, zebrafish embryos have a functional CNS and, presumably, a fully functional BBB. Anatomical structures such as the vascular endothelium can be visualized using transgenic fish lines expressing fluorescent proteins ( ). Defined exposure of the zebrafish can be achieved by the simple addition of solvents to the fish incubation medium within a closed container. Other advantages of the model include the possibility of studying BBB functionality under physiological conditions in vivo and the high throughput. A well-known in vitro model for the human brain endothelium, hCMEC/D3 cell line ( ) showing the formation of tight junctions and the expression of most transporters and receptors of the in vivo BBB [ ], cultured in a transwell system, will also be used. Furthermore, extrapolation of in silico, in vitro, and zebrafish data to higher vertebrates seems feasible [ ].Zebrafish larvae are frequently used in developmental biology or toxicological studies. However, in this study, we will use zebrafish larvae exclusively to study BBB integrity and functionality [
]. Fluorescently labeled reference compounds ( ) will be intravenously injected into the Duct of Cuvier, as markers of paracellular permeability (eg, fluorescein isothiocyanate dextran 70 or fluorescently labeled liposomes), substrates of drug export transporters (eg, rhodamine-123 as P-glycoprotein substrate), or nutrient transporters (eg, fluorescently labeled transferrin as a marker for receptor-mediated transcytosis). PGME will be the first reference compound to be tested because of its high water miscibility. To precisely assess exposure, analytical methods (gas chromatography-tandem mass spectrometry [GC-MS/MS]) will be used to determine the concentrations of solvents and their metabolites in zebrafish medium, in the headspace of closed incubation vessels and tissue samples (ie, zebrafish homogenates). Circulation, tissue distribution, and brain uptake of the reference compounds will be monitored by confocal laser scanning microscopy (live imaging of anesthetized fish embryos for up to 24 hours). The concentration-dependent toxicity of solvents or their metabolites will be monitored based on the viability and malformations of embryos. The integrity of vasculature will be visualized in transgenic zebrafish kdrl: enhanced green fluorescent protein embryos. The metabolic capacity of the zebrafish will be determined by quantifying potential metabolites (determined in WP3; ) in zebrafish tissue homogenates. The concentration-dependent toxicity to BBB and the coefficient of permeation of solvents will additionally be evaluated in the hCMEC/D3 cell line cultured in a transwell system.Finally, quantitative estimates of passive cellular uptake and BBB permeability of solvents and their metabolites will be provided based on computational modeling using physicochemical molecular descriptors according to the methods we previously established [
, ]. These methods provide very high throughput, allowing the screening of web-based chemical libraries.WP3: In Vitro Hepatic Metabolism and Clearance
Because the liver is the main organ responsible for metabolism and a large contributor to compound clearance, we will implement a system suitable for predicting the hepatic metabolism of solvents. In recent years, 3D liver cell models have been proposed as an alternative to less physiological 2D cell monolayers, and their applications have progressed substantially [
]. They are widely used for the assessment of hepatotoxicity [ - ]. An advantage of spheroids is that they overcome the limitation of rapid decline of drug-metabolizing enzyme activities in primary human hepatocyte suspension culture and cell lysates [ ], such as microsomes and liver S9 fractions. In this study, we will use 3D liver cultures of the well-characterized human HepaRG cell line ( ), which represent a promising model to evaluate hepatotoxicity and hepatic metabolism [ ].Determining the appropriate experimental test system (eg, cell plate, incubation time, and exposure concentration) will be an essential part of the development of the 3D model. Moreover, analytical methods (GC-MS/MS and liquid chromatography-tandem mass spectrometry) to detect and quantify the solvents and the metabolites formed must be developed to calculate hepatic metabolism and clearance. Then, proof of metabolic competence and maintenance of the 3D HepaRG cells will be carried out by assessing the metabolism of known P450 substrates. The presence and secretion of albumin as a specific hepatocyte marker will be assessed using immunostaining and enzyme-linked immunosorbent assay. Solvent- and metabolite-induced cytotoxicity will be assessed after 48 hours and 7 days (repeated) of exposure to determine the nontoxic concentration range for subsequent experiments. The metabolic abilities of the 3D HepaRG model and 3D primary human hepatocytes will be compared. Furthermore, the clearance data obtained from the 3D HepaRG model will be compared with the short-term clearance measured in the human liver cell lysate (S9 fractions). In addition, Michaelis-Menten-Kinetic parameters (Vmax and Km) for the formation of the metabolites will be derived using the S9 fractions. These data will be used to build a physiologically based toxicokinetic (PBTK) model (
).WP4: In Vivo Volunteer Exposure
Human biomonitoring refers to monitoring exposure-related health risks by analyzing biological samples, usually blood and urine samples [
]. The biomonitoring limit values (BMLVs) are set to protect human populations against the potential toxic effects of chemical substances. These limit values account for all routes through which a chemical can enter the body. These are most often the inhalation and skin routes in occupational and environmental settings. Kinetic studies that provide absorption, biotransformation, and elimination rates as well as the absorption and elimination half-lives of the parent compound and its metabolites are necessary to set BMLVs. The apparent urinary elimination half-lives of the parent compound and its metabolites will later be used to develop a biomonitoring method. Sample collection time is crucial and is determined by the apparent elimination half-life of the chemical. Blood concentrations will be used to calibrate the air:blood PC for the toxicokinetic models.We will recruit 4 participants for 2 of the selected solvents. All participants must meet the following criteria: they should be healthy individuals who do not smoke or use contraceptive hormones, do not consume alcohol, be aged between 18 and 65 years, have normal red blood cells and hemoglobin concentrations, maintain a BMI between 18 and 25, and should not be working with glycol ethers. Pregnant and breastfeeding women will be excluded from this study. Participants will be recruited using flyers and announcements distributed at the teaching hospital, university websites, and bulletin boards. All participants will sign a written informed consent form before being included in the study.
The participants will be exposed to a single glycol ether for 4 hours under controlled conditions in an exposure chamber (12 m3). PGE concentrations will be set at or below the Swiss occupational exposure level (OEL) if one exists. In the absence of an OEL, we will rely on existing OELs for other propylene glycols. The parent compound (free and conjugated) and the oxidative metabolites (free and conjugated) of the selected glycol ethers will be monitored in blood, urine, and exhaled air samples. These are noninvasive methods used for human participants, and the results will be used in WP5 to estimate brain concentrations. All compounds will be quantified using capillary gas (parent compound in blood, urine, and exhaled air) or liquid (metabolites in blood and urine) chromatograms with tandem mass spectroscopy detection.
WP5: In Silico PBTK Modeling
PBTK models can be used to estimate human brain concentrations. The risk of neurotoxic effects can be estimated by comparing the predicted solvent-brain concentrations with the NOAEC obtained from the in vitro models. Mathematical models such as PBTK models can be used to predict the ADME of a chemical and its metabolites. In these PBTK models, the body is represented by 1 or more compartments. Each compartment represents 1 or more tissues that are kinetically homogeneous, that is, that have similar perfusion rates and an assumed similar substance solubility. PBTK models are described by a set of parameters that define the compartments and a set of mass balance differential equations for each compartment.
A previously developed toxicokinetic model for PGME with metabolism that is assumed to follow Michaelis-Menten kinetics calibrated for different age groups serves as the basis for our development [
, , ]. We aim to modify this previously developed toxicokinetic model to include a separate compartment for the brain using BBB flux rates obtained from in vivo, in vitro, and in silico models (WP2). In addition, we will implement PC obtained from empirical human experiments (WP4) and metabolic parameters assessed in a hepatocyte assay (WP3). The toxicokinetic models will be able to simulate not only acute but also chronic exposures; therefore, both short-term and long-term exposures can be explored in silico. We will develop the toxicokinetic model into a physiologically based pharmacokinetic model based on the existing inhalation-only toxicokinetic model originally developed for PGME [ ] and build it in the Berkeley Madonna software or equivalent. We will model the brain as a single compartment with direct contact with the blood flow and where organic solvent uptake will be assumed to be diffusion limited, which is in line with other physiologically based pharmacokinetic models [ ]. Values for physiological parameters (volume of vascular brain, as fraction of brain volume [FVvb], volume of extravascular brain, as fraction of brain volume [FVevb], volume fraction of brain tissue [FVB; as percent of body weight], BBB surface [Sh] in cm2, fraction of cardiac output in brain at rest [BFbrainrest]/cardiac output in brain at light work [BFbrain]) required to build the TK model are from the scientific literature. Depending on the substance, values of chemical-specific parameters such as the pulmonary retention (Rpulm), central:air PC (Pca), blood:air PC (Pba), and brain tissue:vascular brain PC (Pevb_vb) are either taken from the literature or estimated in silico. Since the partitioning of organic compounds between human tissue homogenate and blood is a function of water and lipid content of tissues and the n-octanol:water PC (Kow), PCs are estimated in silico based on LogKow. Kinetic coefficients needed for each organic solvent included in this study will be found in WP3 for liver metabolism (Michaelis-Menten parameters [Vmax and Km]), WP2 for BBB uptake (BBB permeability-surface area product [PS]). The fraction unbound in blood (Fu_blood) will be estimated based on the fraction unbound in plasma (Fu_plasma) and the blood-to-plasma ratio (Rb), and the fraction unbound in brain (Fu_brain) will be considered when modeling each solvent as only the free fraction is able to distribute to different tissues and is biologically active. The model will be calibrated by comparing the predicted and actual urinary organic solvent concentrations obtained from the controlled human experiments (WP4). Both the free and total organic solvent concentrations (free+conjugated) will be obtained for calibration.Results
With this project, we expect to provide a strategy to rank uncharacterized solvents and their potential liver-formed metabolites, according to their potential neurotoxicity, and in comparison with the banned EGME. More importantly, a series of PBTK simulations will be conducted to predict occupational exposure, assuming 8 hours of exposure per day, 5 days per week, physical activity for 12 hours per day, and rest for the remaining 12 hours. We will use the PBTK model in reverse dosimetry to estimate air concentrations that do not produce brain concentrations determined as neurotoxic in the hiPSC-derived 3D brain model. We will recommend that authorities setting occupational exposure and public health limits consult these values. Keeping the exposure below the brain effect level should ultimately increase the protection of exposed workers and the general population with domestic exposures.
We also anticipate gaining insights into the mechanisms of action of solvents of the glycol ether family. We will elucidate the possible toxic endpoints in the brain, liver, and zebrafish models. Furthermore, we will be able to establish how toxicity is related to the compounds’ lipophilicity and metabolites.
Discussion
Overall, our strategy combining multiple, fit-for-purpose 3D advanced cell culture systems; zebrafish larvae; biomarker analysis; human ADME experiments; and in silico prediction is expected to contribute to the improvement of human risk assessment. Although we identified some risks we could encounter during the project, we are confident that our already determined mitigation measures will be able to overcome potential pitfalls.
Determination of the passage of solvent through the BBB may be challenging; hence, we are applying 3 different complementary methods: in vivo zebrafish larvae, in vitro human cells (hCMEC/D3), and in silico models. We are also considering and assessing the effects of the hepatic metabolites of the solvents on human BBB cells. With this experimental strategy, issues regarding the potential direct effect of solvents on cell membranes, the relatively low miscibility of solvents with water, and the physiological differences between zebrafish and humans (eg, metabolism and route of expected) should be overcome.
We have extensive experience in recruiting human volunteers for controlled human exposure sessions in the exposure chamber. Sometimes, recruitment takes longer than anticipated, and if that is the case, we will extend the timeline to not compromise the size of the study. New analytical chemical methods will need to be determined, which is time consuming. However, we will use a laboratory with extensive experience in analyzing PGME in urine and blood samples. This will also have to be accommodated with a delay in the timeline.
Future developments, not included in this study, are a strategy extended to include developmental neurotoxicity by determining other endpoints, such as proliferation and neurite outgrowth, after exposure of BrainSpheres to solvents at earlier developmental stages and by adding an in vitro test system to take into account the passage of solvents through the placental barrier [
]. We might also consider combining zebrafish embryo behavioral assays (eg, spontaneous tail coiling) with the BrainSpheres model as readouts for developmental neurotoxicity. Finally, the PBTK model could be adapted to determine solvent air concentrations that are unlikely to cause neurotoxic effects in fetuses or pregnant women.Acknowledgments
This work was partially funded by the Swiss Centre for Human Applied Toxicology (SCAHT), Basel, Switzerland.
Data Availability
Data will be available from the corresponding author on reasonable request.
Authors' Contributions
NBH, LSD, JH, and MGZ conceptualized the paper. NBH, LSD, JH, and MGZ wrote the original draft. All authors reviewed and edited the paper. LH, DP, HP, RDP, and SW visualized the paper. NBH, LSD, JH, and MGZ administered the project. NBH, LSD, JH, and MGZ acquired the funding.
Conflicts of Interest
None declared.
References
- Brown RC, Lockwood AH, Sonawane BR. Neurodegenerative diseases: an overview of environmental risk factors. Environ Health Perspect. Sep 2005;113(9):1250-1256. [FREE Full text] [CrossRef] [Medline]
- Pearce N, Kromhout H. Neurodegenerative disease: the next occupational disease epidemic? Occup Environ Med. Sep 2014;71(9):594-595. [CrossRef] [Medline]
- van Valen E, Wekking E, van Hout M, van der Laan G, Hageman G, van Dijk F, et al. Chronic solvent-induced encephalopathy: course and prognostic factors of neuropsychological functioning. Int Arch Occup Environ Health. Oct 2018;91(7):843-858. [FREE Full text] [CrossRef] [Medline]
- Firestone JA, Gospe SMJ. Organic solvents. In: Dobbs MR, editor. Clinical Neurotoxicology: Syndromes, Substances, Environments. Amsterdam, The Netherlands. Elsevier Health Sciences; 2009.
- Bates MN, Reed BR, Liu S, Eisen EA, Hammond SK. Solvent exposure and cognitive function in automotive technicians. Neurotoxicology. Dec 2016;57:22-30. [CrossRef] [Medline]
- Berr C, Vercambre MN, Bonenfant S, Manoux AS, Zins M, Goldberg M. Occupational exposure to solvents and cognitive performance in the GAZEL cohort: preliminary results. Dement Geriatr Cogn Disord. 2010;30(1):12-19. [FREE Full text] [CrossRef] [Medline]
- Meyer-Baron M, Blaszkewicz M, Henke H, Knapp G, Muttray A, Schäper M, et al. The impact of solvent mixtures on neurobehavioral performance: conclusions from epidemiological data. Neurotoxicology. May 2008;29(3):349-360. [CrossRef] [Medline]
- Hachim M, Elemam NM, Maghazachi AA. The beneficial and debilitating effects of environmental and microbial toxins, drugs, organic solvents and heavy metals on the onset and progression of multiple sclerosis. Toxins (Basel). Mar 05, 2019;11(3):147. [FREE Full text] [CrossRef] [Medline]
- Liu M, Shin EJ, Dang DK, Jin CH, Lee PH, Jeong JH, et al. Trichloroethylene and Parkinson's disease: risk assessment. Mol Neurobiol. Jul 2018;55(7):6201-6214. [CrossRef] [Medline]
- Sabbath EL, Gutierrez LA, Okechukwu CA, Singh-Manoux A, Amieva H, Goldberg M, et al. Time may not fully attenuate solvent-associated cognitive deficits in highly exposed workers. Neurology. May 13, 2014;82(19):1716-1723. [FREE Full text] [CrossRef] [Medline]
- Golbabaei F, Dehghani F, Saatchi M, Zakerian SA. Evaluation of occupational exposure to different levels of mixed organic solvents and cognitive function in the painting unit of an automotive industry. Health Promot Perspect. Oct 27, 2018;8(4):296-302. [FREE Full text] [CrossRef] [Medline]
- Filley CM, Heaton RK, Rosenberg NL. White matter dementia in chronic toluene abuse. Neurology. Mar 1990;40(3 Pt 1):532-534. [CrossRef] [Medline]
- Filley CM, Halliday W, Kleinschmidt-DeMasters BK. The effects of toluene on the central nervous system. J Neuropathol Exp Neurol. Jan 2004;63(1):1-12. [CrossRef] [Medline]
- Gregersen P, Mikkelsen S, Klausen H, Døssing M, Nielsen H, Thygesen P. [Chronic cerebral syndrome in painters. Dementia due to inhalation or is it of cryptogenic origin?]. Ugeskr Laeger. Jul 03, 1978;140(27):1638-1644. [Medline]
- Rovida C, Alépée N, Api AM, Basketter DA, Bois FY, Caloni F, et al. Integrated testing strategies (ITS) for safety assessment. ALTEX. 2015;32(1):25-40. [FREE Full text] [CrossRef] [Medline]
- Schmeisser S, Miccoli A, von Bergen M, Berggren E, Braeuning A, Busch W, et al. New approach methodologies in human regulatory toxicology - not if, but how and when!. Environ Int. Aug 2023;178:108082. [FREE Full text] [CrossRef] [Medline]
- Karmaus AL, Mansouri K, To KT, Blake B, Fitzpatrick J, Strickland J, et al. Evaluation of variability across rat acute oral systemic toxicity studies. Toxicol Sci. Jun 28, 2022;188(1):34-47. [FREE Full text] [CrossRef] [Medline]
- Kleinstreuer NC, Ceger PC, Allen DG, Strickland J, Chang X, Hamm JT, et al. A curated database of rodent uterotrophic bioactivity. Environ Health Perspect. May 2016;124(5):556-562. [FREE Full text] [CrossRef] [Medline]
- Luechtefeld T, Maertens A, Russo DP, Rovida C, Zhu H, Hartung T. Global analysis of publicly available safety data for 9,801 substances registered under REACH from 2008-2014. ALTEX. 2016;33(2):95-109. [FREE Full text] [CrossRef] [Medline]
- Ly Pham L, Watford S, Pradeep P, Martin MT, Thomas R, Judson R, et al. Variability in studies: defining the upper limit of performance for predictions of systemic effect levels. Comput Toxicol. Aug 01, 2020;15(August 2020):1-100126. [FREE Full text] [CrossRef] [Medline]
- Wang B, Gray G. Concordance of noncarcinogenic endpoints in rodent chemical bioassays. Risk Anal. Jun 2015;35(6):1154-1166. [CrossRef] [Medline]
- Hartung T. Opinion versus evidence for the need to move away from animal testing. ALTEX. 2017;34(2):193-200. [FREE Full text] [CrossRef] [Medline]
- Hogberg HT, Smirnova L. The future of 3D brain cultures in developmental neurotoxicity testing. Front Toxicol. Jan 27, 2022;4:808620. [FREE Full text] [CrossRef] [Medline]
- Pamies D, Barreras P, Block K, Makri G, Kumar A, Wiersma D, et al. A human brain microphysiological system derived from induced pluripotent stem cells to study neurological diseases and toxicity. ALTEX. 2017;34(3):362-376. [FREE Full text] [CrossRef] [Medline]
- Wardill HR, Mander KA, Van Sebille YZ, Gibson RJ, Logan RM, Bowen JM, et al. Cytokine-mediated blood brain barrier disruption as a conduit for cancer/chemotherapy-associated neurotoxicity and cognitive dysfunction. Int J Cancer. Dec 15, 2016;139(12):2635-2645. [FREE Full text] [CrossRef] [Medline]
- Bendels S, Kansy M, Wagner B, Huwyler J. In silico prediction of brain and CSF permeation of small molecules using PLS regression models. Eur J Med Chem. Aug 2008;43(8):1581-1592. [CrossRef] [Medline]
- Gupta PK. Principles and basic concepts of toxicokinetics. In: Fundamentals of Toxicology: Essential Concepts and Applications. Cambridge, MA. Academic Press; 2016.
- Breen M, Ring CL, Kreutz A, Goldsmith MR, Wambaugh JF. High-throughput PBTK models for in vitro to in vivo extrapolation. Expert Opin Drug Metab Toxicol. Aug 2021;17(8):903-921. [FREE Full text] [CrossRef] [Medline]
- Substance regulatory obligations. European Chemicals Agency. URL: https://echa.europa.eu/legislation-obligation/-/obligations/100.003.377 [accessed 2023-12-19]
- Hopf NB, Vernez D, Berthet A, Charriere N, Arnoux C, Tomicic C. Effect of age on toxicokinetics among human volunteers exposed to propylene glycol methyl ether (PGME). Toxicol Lett. May 20, 2012;211(1):77-84. [FREE Full text] [CrossRef] [Medline]
- Reale E, Sandstrom J, Culot M, Hechon J, Wellens S, Heymans M, et al. Predicting human neurotoxicity of propylene glycol methyl ether (PGME) by implementing in vitro neurotoxicity results into toxicokinetic modelling. Sci Total Environ. Aug 15, 2023;886:163767. [FREE Full text] [CrossRef] [Medline]
- Filley CM, Kleinschmidt-DeMasters BK. Toxic leukoencephalopathy. N Engl J Med. Aug 09, 2001;345(6):425-432. [CrossRef] [Medline]
- Pamies D, Block K, Lau P, Gribaldo L, Pardo CA, Barreras P, et al. Rotenone exerts developmental neurotoxicity in a human brain spheroid model. Toxicol Appl Pharmacol. Sep 01, 2018;354:101-114. [FREE Full text] [CrossRef] [Medline]
- Zhong X, Harris G, Smirnova L, Zufferey V, Sá RC, Baldino Russo F, et al. Antidepressant paroxetine exerts developmental neurotoxicity in an iPSC-derived 3D human brain model. Front Cell Neurosci. Feb 21, 2020;14:25. [FREE Full text] [CrossRef] [Medline]
- Nunes C, Singh P, Mazidi Z, Murphy C, Bourguignon A, Wellens S, et al. An in vitro strategy using multiple human induced pluripotent stem cell-derived models to assess the toxicity of chemicals: a case study on paraquat. Toxicol In Vitro. Jun 2022;81:105333. [FREE Full text] [CrossRef] [Medline]
- Chesnut M, Paschoud H, Repond C, Smirnova L, Hartung T, Zurich MG, et al. Human IPSC-derived model to study myelin disruption. Int J Mol Sci. Aug 31, 2021;22(17):9473. [FREE Full text] [CrossRef] [Medline]
- Goldstein DB. The effects of drugs on membrane fluidity. Annu Rev Pharmacol Toxicol. 1984;24:43-64. [CrossRef] [Medline]
- Bowen SE, Batis JC, Paez-Martinez N, Cruz SL. The last decade of solvent research in animal models of abuse: mechanistic and behavioral studies. Neurotoxicol Teratol. 2006;28(6):636-647. [CrossRef] [Medline]
- Balster RL. Neural basis of inhalant abuse. Drug Alcohol Depend. 1998;51(1-2):207-214. [CrossRef] [Medline]
- Fukumori R, Takarada T, Nakazato R, Fujikawa K, Kou M, Hinoi E, et al. Selective inhibition by ethanol of mitochondrial calcium influx mediated by uncoupling protein-2 in relation to N-methyl-D-aspartate cytotoxicity in cultured neurons. PLoS One. Jul 16, 2013;8(7):e69718. [FREE Full text] [CrossRef] [Medline]
- Revilla AS, Pestana CR, Pardo-Andreu GL, Santos AC, Uyemura SA, Gonzales ME, et al. Potential toxicity of toluene and xylene evoked by mitochondrial uncoupling. Toxicol In Vitro. Aug 2007;21(5):782-788. [CrossRef] [Medline]
- Zurich MG, Stanzel S, Kopp-Schneider A, Prieto P, Honegger P. Evaluation of aggregating brain cell cultures for the detection of acute organ-specific toxicity. Toxicol In Vitro. Jun 2013;27(4):1416-1424. [CrossRef] [Medline]
- Nunes C, Proença S, Ambrosini G, Pamies D, Thomas A, Kramer NI, et al. Integrating distribution kinetics and toxicodynamics to assess repeat dose neurotoxicity using human BrainSpheres: a case study on amiodarone. Front Pharmacol. Sep 6, 2023;14:1248882. [FREE Full text] [CrossRef] [Medline]
- Sieber S, Grossen P, Bussmann J, Campbell F, Kros A, Witzigmann D, et al. Zebrafish as a preclinical in vivo screening model for nanomedicines. Adv Drug Deliv Rev. 2019;151-152:152-168. [CrossRef] [Medline]
- Rizzo LY, Golombek SK, Mertens ME, Pan Y, Laaf D, Broda J, et al. In vivo nanotoxicity testing using the zebrafish embryo assay. J Mater Chem B. Jun 10, 2013;1:10.1039/C3TB20528B. [FREE Full text] [CrossRef] [Medline]
- Weksler B, Romero IA, Couraud PO. The hCMEC/D3 cell line as a model of the human blood brain barrier. Fluids Barriers CNS. Mar 26, 2013;10(1):16. [FREE Full text] [CrossRef] [Medline]
- Sieber S, Grossen P, Detampel P, Siegfried S, Witzigmann D, Huwyler J. Zebrafish as an early stage screening tool to study the systemic circulation of nanoparticulate drug delivery systems in vivo. J Control Release. Oct 28, 2017;264:180-191. [CrossRef] [Medline]
- Quiñonez-Silvero C, Hübner K, Herzog W. Development of the brain vasculature and the blood-brain barrier in zebrafish. Dev Biol. Jan 15, 2020;457(2):181-190. [FREE Full text] [CrossRef] [Medline]
- Suenderhauf C, Hammann F, Huwyler J. Computational prediction of blood-brain barrier permeability using decision tree induction. Molecules. Aug 31, 2012;17(9):10429-10445. [FREE Full text] [CrossRef] [Medline]
- Duval K, Grover H, Han LH, Mou Y, Pegoraro AF, Fredberg J, et al. Modeling physiological events in 2D vs. 3D cell culture. Physiology (Bethesda). Jul 2017;32(4):266-277. [FREE Full text] [CrossRef] [Medline]
- Prestigiacomo V, Weston A, Messner S, Lampart F, Suter-Dick L. Pro-fibrotic compounds induce stellate cell activation, ECM-remodelling and Nrf2 activation in a human 3D-multicellular model of liver fibrosis. PLoS One. Jun 30, 2017;12(6):e0179995. [FREE Full text] [CrossRef] [Medline]
- Prestigiacomo V, Suter-Dick L. Nrf2 protects stellate cells from Smad-dependent cell activation. PLoS One. Jul 20, 2018;13(7):e0201044. [FREE Full text] [CrossRef] [Medline]
- Riede J, Wollmann BM, Molden E, Ingelman-Sundberg M. Primary human hepatocyte spheroids as an in vitro tool for investigating drug compounds with low clearance. Drug Metab Dispos. Jun 01, 2021:DMD-AR-2020-000340. [FREE Full text] [CrossRef] [Medline]
- Ingelman-Sundberg M, Lauschke VM. 3D human liver spheroids for translational pharmacology and toxicology. Basic Clin Pharmacol Toxicol. Jan 2022;130 Suppl 1:5-15. [FREE Full text] [CrossRef] [Medline]
- Proctor WR, Foster AJ, Vogt J, Summers C, Middleton B, Pilling MA, et al. Utility of spherical human liver microtissues for prediction of clinical drug-induced liver injury. Arch Toxicol. Aug 2017;91(8):2849-2863. [FREE Full text] [CrossRef] [Medline]
- Preiss LC, Lauschke VM, Georgi K, Petersson C. Multi-well array culture of primary human hepatocyte spheroids for clearance extrapolation of slowly metabolized compounds. AAPS J. Mar 11, 2022;24(2):41. [CrossRef] [Medline]
- Ramaiahgari SC, Waidyanatha S, Dixon D, DeVito MJ, Paules RS, Ferguson SS. From the cover: three-dimensional (3D) HepaRG spheroid model with physiologically relevant xenobiotic metabolism competence and hepatocyte functionality for liver toxicity screening. Toxicol Sci. Sep 01, 2017;159(1):124-136. [FREE Full text] [CrossRef] [Medline]
- Costa C, Teixeira JP. Biomonitoring. In: Wexler P, editor. Encyclopedia of Toxicology (Third Edition). Cambridge, MA. Academic Press; 2014.
- Tomicic C, Berode M, Oppliger A, Castella V, Leyvraz F, Praz-Christinaz SM, et al. Sex differences in urinary levels of several biological indicators of exposure: a human volunteer study. Toxicol Lett. May 10, 2011;202(3):218-225. [FREE Full text] [CrossRef] [Medline]
- Jongeneelen F, ten Berge W. Simulation of urinary excretion of 1-hydroxypyrene in various scenarios of exposure to polycyclic aromatic hydrocarbons with a generic, cross-chemical predictive PBTK-model. Int Arch Occup Environ Health. Aug 2012;85(6):689-702. [CrossRef] [Medline]
- Cherubini M, Erickson S, Haase K. Modelling the human placental interface in vitro-a review. Micromachines (Basel). Jul 27, 2021;12(8):884. [FREE Full text] [CrossRef] [Medline]
Abbreviations
ADME: absorption, distribution, metabolism, and excretion |
BBB: blood-brain barrier |
BFbrain: cardiac output in brain at light work |
BFbrainrest: cardiac output in brain at rest |
BMLV: biomonitoring limit value |
CNS: central nervous system |
EGE: ethylene glycol ether |
EGME: ethylene glycol methyl ether |
Fu_blood: fraction unbound in blood |
Fu_brain: fraction unbound in brain |
Fu_plasma: fraction unbound in plasma |
FVB: volume fraction of brain tissue |
FVevb: volume of extravascular brain, as fraction of brain volume |
FVvb: volume of vascular brain, as fraction of brain volume |
GC-MS/MS: gas chromatography-tandem mass spectrometry |
hCMEC/D3: human brain microvascular endothelial cells |
hiPSC: human induced pluripotent stem cell |
NOAEC: no observed adverse effect concentration |
OECD: Organisation for Economic Co-operation and Development |
OEL: occupational exposure level |
Pba: blood:air partition coefficient |
PBTK: physiologically based toxicokinetic |
PC: partition coefficient |
Pca: central:air partition coefficient |
Pevb_vb: brain tissue:vascular brain partition coefficient |
PGE: propylene glycol ether |
PGME: propylene glycol methyl ether |
PS: blood-brain barrier permeability-surface area product |
Rb: blood-to-plasma ratio |
Rpulm: pulmonary retention |
Sh: blood-brain barrier surface |
TG: Test Guideline |
Vmax: maximum velocity |
WP: work package |
Edited by A Mavragani; submitted 26.06.23; peer-reviewed by T Dong, J Parkin Kullmann; comments to author 18.09.23; revised version received 10.10.23; accepted 13.10.23; published 18.01.24.
Copyright©Nancy B Hopf, Laura Suter-Dick, Jörg Huwyler, Myriam Borgatta, Lucie Hegg, David Pamies, Hélène Paschoud, Ramya Deepthi Puligilla, Elena Reale, Sophie Werner, Marie-Gabrielle Zurich. Originally published in JMIR Research Protocols (https://www.researchprotocols.org), 18.01.2024.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Research Protocols, is properly cited. The complete bibliographic information, a link to the original publication on https://www.researchprotocols.org, as well as this copyright and license information must be included.