Combined Toxic Effects of Polar and Nonpolar Chemicals on Human Hepatocytes (HepG2) Cells by Quantitative Property-Activity Relationship Modeling.

We determined the toxicity of mixtures of ethyl acetate (EA), isopropyl alcohol (IPA), methyl ethyl ketone (MEK), toluene (TOL) and xylene (XYL) with half-maximal effective concentration (EC50) values obtained using human hepatocytes cells. According to these data, quantitative property-activity relationships (QPAR) models were successfully proposed to predict the toxicity of mixtures by multiple linear regressions (MLR). The leave-one-out cross validation method was used to find the best subsets of descriptors in the learning methods. Significant differences in physico-chemical properties such as boiling point (BP), specific gravity (SG), Reid vapor pressure (rVP) and flash point (FP) were observed between the single substances and the mixtures. The EC50 of the mixture of EA and IPA was significantly lower than that of contained TOL and XYL. The mixture toxicity was related to the mixing ratio of MEK, TOL and XYL (MLR equation EC50 = 3.3081 - 2.5018 × TOL - 3.2595 × XYL - 12.6596 × MEK × XYL), as well as to BP, SG, VP and FP (MLR equation EC50 = 1.3424 + 6.2250 × FP - 7.1198 × SG × FP - 0.03013 × rVP × FP). These results suggest that QPAR-based models could accurately predict the toxicity of polar and nonpolar mixtures used in rotogravure printing industries.

INTRODUCTION

According to the 2014 Environmental Report (1), over 120,000 different chemicals are dispersed worldwide, and over 2,000 new chemical variations are developed and commercialized each year. More than 40,000 chemicals are distributed in Korea, with over 300 new chemicals introduced to the domestic market annually. Thus, the number of chemicals and their distribution are predicted to increase. The primary reason for the development of novel chemicals and their increased distribution is their unique physical/chemical properties, which allow their utilization in various fields. Numerous chemicals are used in Korean industry. In particular, paint manufacturing and printing businesses use high amounts of diverse chemicals. Thus, the risk of health hazards for the employees in these facilities is proportional to the number and amount of chemicals being used. Gravure printing facilities are known to use various chemicals, including benzene (2). These resulted in the regulation of highly toxic chemical use, including benzene. When chemicals exposure at a gravure printing business in Korea was analyzed, the main components in the ink were toluene (CAS No., 108-88-3, TOL), methyl ethyl ketone (CAS No., 78-93-3, MEK), ethyl acetate (CAS No., 141-78-6, EA), isopropyl alcohol (CAS No., 67-63-0, IPA), and xylene (CAS No., 1330-20-7, XYL). Additionally, some inks contained 2-butanol (CAS No., 78-92-2), cyclohexane (CAS No., 110-82-7), and 2-ethoxyethanol (CAS No., 110-80-5) (3). Other studies have also reported that workers in the gravure printing business are exposed to a variety of chemicals, including aromatic hydrocarbons (ethylbenzene, styrene, TOL, XYL, cyclohexane), aliphatic hydrocarbons, ketone bodies, and alcohols. Furthermore, the estimated mixture concentration is higher in gravure printing than in screen or offset printing, indicating that workers in gravure printing businesses are likely exposed to higher chemical concentrations (4). The possibility of health hazards increases with exposure to higher concentrations of a single chemical. However, in the case of mixtures, the effect of a single chemical may be altered owing to interactions between chemicals. Kim et al. (5) reported that dimethylformamide (DMF), TOL, and MEK, which are the major components used in synthetic leather manufacturing, significantly differ in their physical/chemical properties as a single substance and a mixture. Thus, exposure levels can vary depending on their handling form (single or mixture). Interestingly, a previous study observed noticeable hemocyte necrosis after injection of a DMF/TOL mixture, but not with a DMF/MEK mixture (6). This phenomenon can be explained as synergistic hemotoxicity between DMF and TOL, and antagonistic effects between DMF and MEK (7). Importantly, human toxicity prediction studies using mixed substances are ongoing; however, they are insufficient compared to single substance studies. This is most likely because mixed substance studies require more time and have higher costs, since numerous combinations must be used. Furthermore, the number of human and animal studies is restricted owing to research ethics requirements. Therefore, researchers are pursuing studies complementing the aforementioned limitations and producing reliable results. Quantitative property-activity relationships (QPAR) has been used in this field to find key properties (ex. geometric solids, electronic properties, and physical/chemical properties, etc.) without actual experiments or biological information to predict human toxicity, efficacy, and drug reactivity through statistical modeling (8). These studies are primarily performed in the fields of medicine/pharmacology, environmental science, toxicology, and biochemistry, and are gradually expanding to other fields (9,10). Importantly, although the majority of employees handling manufacturing chemicals are exposed to mixed substances, human toxicity studies using mixed substances are limited. Thus, the use of QPAR studies could overcome this challenge. The goal of this study is to compare and analyze the physical/chemical properties and half-maximal effective concentration (EC50) of mixed chemicals handled in gravure printing facilities. We further sought to predict the human toxicity of these mixed chemicals through QPAR analysis.

MATERIALS AND METHODS

Targeted chemicals and mixing composition

We examined the component ratio of chemicals and mixed substances related to gravure printing facilities through data collection. Data were divided into experimental groups, as shown in Table 1.

Table 1

Classification of experimental groups and mixing ratio of chemicals

Experimental groups	Mixing ratio (vol/vol)
Single
Ethyl acetate (EA)
Isopropyl alcohol (IPA)
Methyl ethyl ketone (MEK)
Toluene (TOL)
Xylene (XYL)
Mixtures
MEK + EA (G1)	9 : 1
MEK + IPA (G2)	9 : 1
MEK + EA + IPA (G3)	8 : 1 : 1
EA + IPA + MEK + TOL (G4)	3 : 1 : 1 : 5
EA + IPA + MEK + XYL (G5)	3 : 1 : 1 : 5
IPA + MEK + TOL (G6)	1 : 1 : 5
IPA + MEK + XYL (G7)	1 : 1 : 5
EA + IPA + MEK + TOL + XYL (G8)	3 : 1 : 1 : 2.5 : 2.5

Letters in parenthesis represents the groups (G1–G8).

Ambient monitoring in workplace

Based on the collected data, the chemicals in the target facilities were confirmed using a manual sampling device (Organic Vapor Monitor 3500, 3M, St. Paul, MN, USA) and active sampling device (LFS-113, Gilian, St. Petersburg, FL, USA) attached to a charcoal tube (Cat. No. 226-01, 226-81, SKC Ltd., Eighty Four, PA, USA). Samples collected by a manual sampling device were subjected to desorption for 1 hr with 2 mL CS2 (1% butanol). Samples from the charcoal tube collected by the active sampling device were subjected to desorption with 1 mL CS2 (1% butanol). Samples were analyzed according to the official test method 1501 (11) and 2500 (12), as recommended by National Institute of Occupational Safety and Health (NIOSH).

Physico-chemical properties of chemicals

Data on physic-chemical properties of each EA, IPA, MEK, TOL and XYL were obtained from material safety data sheet (MSDS) provided by the Korea Occupational Safety and Health Agency. Physico-chemical properties such as boiling point (BP) (13), specific gravity (SG) (14), Reid vapor pressure (rVP) (15) and flash point (FP) (16) in mixtures were measured.

Experimental determination of half-maximal effective concentrations (EC50)

HepG2 cells (human hepatocytes) were acquired from the Korean Cell Line Bank. The cells were cultivated in DMEM (10% FBS, 100 unit/mL penicillin, and 100 μg/mL streptomycin) and MEM badges in a 5% CO2 atmosphere at 37°C. They were seeded onto 96-well plates (Corning Inc., Corning, New York, USA) at a concentration of 5 × 104 cells/well and were mixed with single and chemical mixtures after 24 hrs. The cells were transferred into MEM culture medium to be stabilized for 24 hrs at a concentration of 5 × 104 cells/well in a 96-well plate. Samples were grouped as shown in Table 1 and treated with the appropriate solutions for 3, 6, 12, and 24 hrs. After culture, the medium was removed and the cells were diluted by 1/10 with the CCK-8 assay kit (Woongbee, Seoul, Korea) in DMEM culture medium, and then incubated for 1.5 hrs at 37°C. EC50 values were determined by measuring the absorption of formazan at 450 nm (7).

Statistical validation of the model

The goodness-of-fit with the physico-chemical properties or mixing ratio of EA, IPA, MEK, TOL and XYL, as well as the EC50 was carried out using SAS JMP PRO software (ver. 12.0, SAS Institute Inc., Cary, NC, USA). The descriptors were subjected to stepwise selection of multiple linear regression analysis which resulted in descriptors model with least squares fit (R2) and leave-one-out cross-validated R2 (LOO Q2). The performance of the QPAR regression models has been evaluated by R2, RMSE, and Q2 as follows.

The QPAR regression model was analyzed by measuring the degree of relevance between the properties used and the estimated EC50 using the linear regression equation (1) below with a training set divided into the pre-set number of folds following a cross-validation test.

Kruskal-Wallis test was performed to assess the differences in physico-chemical properties among the experimental groups, and the QPAR regression model was cross-validated. All results were presented as mean percentages with standard deviation.

RESULTS

Identification of chemical substances

Chemical substances handled in the gravure printing facilities were confirmed by ambient monitoring (Table 2). Handling substances differed depending on the printing product; however, the ambient monitoring results showed that EA, MEK, IPA, TOL, and XYL are primarily used.

Table 2

Identification and ambient levels of chemicals

Companies	Mean concentrations of ambient chemicals (ppm)
	EA	IPA	MEK	Toluene	Xylene
Rotogravure paint manufacturing companies (2 companies)
A	3.4 ± 5.1	0.00	0.00	12.7 ± 16.1	12.6 ± 19.1
B	20.6 ± 19.9	3.2 ± 2.6	20.6 ± 33.1	26.5 ± 32.2	21.5 ± 21.8
Rotogravure printing companies (3 companies)
C	210.4 ± 51.1	0.00	81.9 ± 24.7	77.6 ± 28.7	0.00
D	218.1 ± 291.9	20.9 ± 51.1	219.4 ± 307.4	38.9 ± 55.7	0.7 ± 2.1
E	140.8 ± 52.8	4.7 ± 2.2	178.7 ± 69.1	0.00	3.7 ± 2.4
Total	132.1 ± 170.5	7.2 ± 26.5	125.8 ± 179.1	26.8 ± 41.5	6.6 ± 13.6
TWA	400	200	200	50	100

EA, ethyl acetate; IPA, isopropyl alcohol; MEK, methyl ethyl ketone; TWA, time weighted average.

Composition of physico-chemical properties by experimental groups

The BP, SG, rVP and FP were measured (Table 3). The BP of TOL and XYL alone were 110.0°C and 139.2°C respectively. These significantly decreased when mixed with EA, MEK, and IPA (p < 0.01). The SG of EA alone (0.90 g/mL) significantly decreased when mixed with other substances (p < 0.01). The rVP of TOL and XYL significantly increased when mixed with EA, IPA, and MEK, as compared to the substance alone (p < 0.01). The FP of EA, which has a low FP value (−4.0°C), was significantly decreased (p < 0.01) when mixed with MEK (−9.0°C), whereas the FP of MEK was significantly increased when mixed with other substances.

Table 3

Results of the physico-chemical properties

Experimental groups	Physico-chemical properties
	BP (°C)	SG (g/mL)	rVP (kPa)	FP (°C)
Ethyl acetate (EA)	77.0	0.90	22.5	−4.0
Isopropyl alcohol (IPA)	83.0	0.79	12.1	11.7
Methyl ethyl ketone (MEK)	80.0	0.80	21.6	−9.0
Toluene (TOL)	111.0	0.86	7.0	4.0
Xylene (XYL)	139.2	0.87	2.3	29.0
MEK + EA (G1)	78.7 ± 0.1	0.8181 ± 0.0008	21.23 ± 0.21	−7.37 ± 0.29
MEK + IPA (G2)	78.5 ± 0.2	0.8015 ± 0.0001	21.60 ± 0.36	−6.85 ± 0.29
MEK + EA + IPA (G3)	77.6 ± 0.4	0.8110 ± 0.0001	22.47 ± 0.38	−6.53 ± 0.29
EA + IPA + MEK + TOL (G4)	83.9 ± 0.3	0.8615 ± 0.0001	18.00 ± 0.44	−2.15 ± 0.50
EA + IPA + MEK + XYL (G5)	88.1 ± 0.3	0.8613 ± 0.0001	14.63 ± 0.12	0.43 ± 0.01
IPA + MEK+TOL (G6)	87.9 ± 0.1	0.8466 ± 0.0001	15.93 ± 0.15	−0.57 ± 0.50
IPA + MEK + XYL (G7)	94.4 ± 0.2	0.8454 ± 0.0001	11.77 ± 0.12	6.31 ± 0.44
EA + IPA + MEK + TOL + XYL (G8)	86.1 ± 0.2	0.8610 ± 0.0001	16.20 ± 0.20	−1.01 ± 0.25
χ2 (Kuskal Wallis test)	22.375	22.683	22.426	22.406
P value	0.01	0.01	0.01	0.01

BP, boiling point; SG, specific gravity; rVP, Reid vapor pressure; FP, flash point.

Letters in parenthesis represents the groups (G1–G8).

EC50 values

The EC50 values for each experimental group were measured using HepG2 cells (Fig. 1). The EC50 values in the single substance treated groups were 4.427 μL/100 μL for EA, 4.341 μL/100 μL for IPA, 2.505 μL/100 μL for MEK, respectively, the EC50 values were lower in mixed substances 0.792 μL/100 μL for TOL, and 0.146 μL/100 μL for XYL than single substances, and the differences between the experimental groups were dependent on the EC50 values of the single substances.

Fig. 1

EC50 (HepG2) values. EA, ethyl acetate; IPA, isopropyl alcohol; MEK, methyl ethyl ketone; G1, MEK + EA; G2, MEK + IPA; G3, MEK + EA + IPA; G4, EA + IPA + MEK + Toluene; G5, EA + IPA + MEK + Xylene; G6, IPA + MEK + Toluene; G7, IPA + MEK + Xylene; G8, EA + IPA + MEK + Toluene + Xylene.

Regression model for predicting EC50 in HepG2 cells

Regression model by chemical mixing ratios

To predict human toxicity according to the mixing ratio of EA, IPA, MEK, TOL, and XYL, optimized regression models were obtained using a linear regression (Table 4). When we compared the goodness of fit (R2) for the training data through regression models by chemical mixing ratios, mixed substances containing TOL and XYL showed goodness of fit between 0.7385 and 0.9114. To determine the descriptors that affect the model optimization, leave-one-out cross-validation was performed using the leave-many-out (LMO) method. Our results indicated that R2 ranged from 0.3670 to 0.7197 for mixed substances containing TOL and XYL. The multiple linear regression (MLR) models by mixing ratios of TOL + XYL (Q2 = 0.6006) and MEK + TOL + XYL (Q2 = 0.6053) showed Q2 over 0.5, which is generally considered an acceptance level, but models containing EA were Q2 < 0.5, which is a standard for validity verification. Furthermore, models containing MEK * XYL as TOL + XYL + MEK * XYL (Q2 = 0.7197, y = 3.3081 − 2.5018 * TOL − 3.2595 * XYL − 12.6596 * MEK * XYL) and MEK + TOL + XYL + MEK * XYL (Q2 = 0.6678, y = 4.1971 − 1.5544 * MEK − 3.5487 * TOL − 4.2526 * XYL − 14.2728 * MEK * XYL) gave cross-validated Q2 value of more than 0.6 and showed that the EC50 value decreased with increasing TOL, XYL and MEK * XYL concentrations. Model containing TOL + XYL + MEK * XYL showed good predictive ability owing to highest cross-validated Q2 value.

Table 4

QPAR model according to chemical mixing ratio by linear regression equation

Solvent groups	Training		LOO CV	Linear regression equation
	R2	RMSE	Q2
MEK, XYL	0.4487	1.1305	0.0300	y = 2.6437 + 0.0293 * MEK − 2.9631 * XYL
TOL, XYL	0.7385	0.7786	0.6006	y = 3.2074 − 2.4922 * TOL − 3.6719 * XYL
EA, TOL, XYL	0.7980	0.7214	0.4725	y = 2.9383 + 1.2337 * EA − 2.2903 * TOL − 3.4700 * XYL
MEK, TOL, XYL	0.8502	0.6211	0.6053	y = 4.0349 − 1.4681 * MEK − 3.4798 * TOL − 4.6595 * XYL
TOL, XYL, MEK * XYL	0.7870	0.7408	0.7197	y = 3.3081 − 2.5018 * TOL − 3.2595 * XYL − 12.6596 * MEK * XYL
MEK, TOL, XYL, MEK * XYL	0.9114	0.5006	0.6678	y = 4.1971 − 1.5544 * MEK − 3.5487 * TOL − 4.2526 * XYL − 14.2728 * MEK * XYL
EA, MEK, TOL, XYL	0.8560	0.6459	0.3670	y = 3.8154 + 0.4561 * EA − 1.2552 * MEK − 3.2620 * TOL − 4.4416 * XYL

EA, ethyl acetate; IPA, isopropyl alcohol; MEK, methyl ethyl ketone; TOL, toluene; XYL, xylene; LOO CV, leave-one-out cross-validation; R2 and Q2, coefficient; RMSE, root mean square error.

Physico-chemical regression model

Using the physico-chemical properties BP, SG, rVP, and FP as descriptors for the linear regression equation, an optimized regression model was obtained (Table 5). Comparing the goodness of fit (R2) of the training data in each experimental group by cross-validation showed that, when more than two descriptors were selected, the goodness of fit (R2) was 0.5434–0.9623. Goodness of fit was observed to be high when rVP and FP was combined compared in other physic-chemical property. We conducted the leave-one-out cross validation to determine the descriptors that affect the model optimization. As a results, the models with rVP * FP R2 > 0.9 and Q2 > 0.5. Model containing FP + SG * FP, rVP * FP showed good predictive ability owing to highest cross-validated Q2 value and revealed that EC50 value decreased with increasing SG and rVP values, and with decreasing FP values.

Table 5

QPAR model according to physico-chemical properties by linear regression equation

Physico-chemical properties	Training		LOO CV	Linear regression equation
	R2	RMSE	Q2
BP	0.4811	1.0457	0.2215	y = 7.0812 − 0.05513 * BP
rVP	0.4185	1.1070	0.2654	y = −0.1222 + 0.1419 * rVP
BP, FP	0.5915	0.9731	0.4664	y = 11.0604 − 0.1006 * BP + 0.08954 * FP
rVP, FP	0.5434	1.0288	0.4570	y = −2.8460 + 0.3052 * rVP + 0.1116 * FP
BP, rVP, FP	0.6658	0.9277	0.1462	y = 5.3946 − 0.06978 * BP + 0.1732 * rVP + 0.1393 * FP
BP * SG, BP * FP	0.5950	0.9687	0.3622	y = 10.4646 − 0.1127 * BP * SG = 0.0007994 * BP * FP
BP, SG * FP	0.5910	0.9734	0.4680	y = 11.3045 − 0.1034 * BP + 0.1104 * SG * FP
FP, SG * FP, rVP * FP	0.9479	0.3664	0.5932	y = 1.3424 + 6.2250 * FP − 7.1198 * SG * FP − 0.03013 * rVP * FP
SG, FP, SG * FP, rVP * FP	0.9623	0.3303	0.5381	y = 5.7976 − 5.2754 * SG + 6.1237 * FP − 7.0023 * SG * FP − 0.02938 * rVP * FP

BP, boiling point; FP, flash point; rVP, Reid vapor pressure; SG, specific gravity; LOO CV, leave-one-out cross-validation; R2 and Q2, coefficient; RMSE, root mean square error.

DISCUSSION

Chemical substance development and production is useful in many fields. However, employees who produce these chemicals have a higher risk for health problems, due to chemical exposure (17). Because most of these employees are exposed to mixed compounds (6), it is difficult to elucidate the chemical substances that directly affect health, as well as the interactions between these chemicals. Therefore, we evaluated the chemical substances in mixed compounds and determined which physico-chemical properties have the greatest effect on predicting human toxicity after mixed compound exposure, using single and mixed compounds and EC50 value as descriptors. When we compared the physico-chemical properties and EC50 values of single and mixed compounds, we found that human toxicity due to mixed compounds was dependent on the physico-chemical properties of the single substances. When TOL and XYL were mixed with EA, IPA, and MEK, the BP and FP were lower, whereas the rVP was higher (Table 3). Furthermore, the EC50 values of chemicals with lower toxicity (EA and EPA) decreased (toxicity increased) when they were mixed with other highly toxic chemicals, similar to the changes in physico-chemical properties. Furthermore, the EC50 value of TOL and XYL, which have relatively high toxicity, increased (toxicity decreased) when they were mixed with low toxicity chemicals (Fig. 1). The individual substances composing mixed compounds undergo selective interaction, due to differences in physico-chemical properties and structures. As a result, differences in substance metabolism and toxicity are observed (18). Many studies are being conducted on the absorption and metabolic interactions of mixed compounds in the body; (19,20) however, studies determining toxicity and its mechanism are lacking (21,22). Liira et al. (23) reported that MEK inhibited the metabolism of XYL. Additionally, a study by Tardif et al. (24) reported that mixed compounds with low TOL and XYL concentrations had no effect on substance metabolism at high concentrations. Interestingly, TOL and XYL have metabolic interactions, and inhibited the metabolism of the other. Freundt et al. (25) stated that, after exposure to combined TOL and EA, EA increased the metabolism of TOL, and that this increase was associated with the exposure concentration. Thus, metabolic interactions in mixed compounds are associated with physico-chemical properties and structural similarities, and are dependent on the exposure level (19). Xenobiotics can induce toxicity following absorption by altering pharmacokinetic processes, including absorption, distribution, metabolism, and excretion, depending on their chemical structure and physico-chemical properties. Toxicity can occur at any step in these processes. In contrast, mixed compounds undergo more complicated metabolic processes than single substances, due to the individual pharmacokinetics of the single substances and the pharmacokinetic/pharmacodynamic interactions between individual substances. Human toxicity occurs largely owing to functional inhibition or decreased interaction between metabolic intermediates and adducts, which are generated during the metabolic process, and homeostatic regulatory factors. This is because interactions between components in mixed compounds either increases or inhibits metabolism, thereby antagonizing or synergistically enhancing toxicity (7). When we compared the cell survival rate and EC50 values to examine toxicity between single substances (EA, IPA, MEK, TOL, and XYL) and mixed compounds, there were many differences between the two groups. The cell survival rate was lower with mixed compounds than single substances, and mixing TOL and XYL with EA, IPA, and MEK significantly decreased the cell survival rate. These seem to occur due to changes in physico-chemical properties following interactions between single substances comprising the mixed compound (5). Furthermore, the cell survival rate was positive correlated with the EC50 value in this study. EC50 values for EA and IPA alone were 4.427 μL/100 μL and 4.341 μL/100 μL, respectively; however, they decreased slightly when mixed with MEK and 4-fold when mixed with TOL and XYL. Interestingly, XYL showed a marked decrease, as compared to TOL. When the toxicity level is compared with the time weight average (TWA) standard of the Ministry of Employment and Labor, TOL toxicity is higher (TOL is 50 ppm and XYL is 100 ppm) and LD50 value is lower (2,600 mg/kg rat) than XYL (3,500 mg/kg rat), indicating that TOL is more toxic. However, Croute et al. (18) reported that TOL cytotoxicity is higher than that of benzene, owing to the side chain attached to the benzene ring and the effect of the lipophilic properties. The EC50 value of XYL is lower than that of TOL in this study, likely due to the structural properties of the chemicals. In this case, it is presumed to be the effect of the methyl group on the benzene ring. Thus, when the EC50 and LD50 values are used as biological indices to evaluate environmental and human toxicity by chemical substances, evaluation properties should be considered. The toxicity of mixed substances must be predicted using physico-chemical properties, due to the time, cost, and ethical considerations required for experimental research. Thus, in this study, we predicted the toxicity with a QPAR liner regression model using physico-chemical properties and EC50 values as descriptors. As a result, the Q2 value for the TOL + XYL and MEK + TOL + XYL model were 0.6006 and 0.6053, respectively, exceeding the standard for goodness of fit (Q2 = 0.5). The linear regression analysis showed that TOL and XYL reduced the EC50 value among the components of the mixed compounds. In particular, models containing MEK * XYL as TOL + XYL + MEK * XYL (Q2 = 0.7197, y = 3.3081 − 2.5018 * TOL − 3.2595 * XYL − 12.6596 * MEK * XYL) and MEK + TOL + XYL + MEK * XYL (Q2 = 0.6678, y = 4.1971 − 1.5544 * MEK − 3.5487 * TOL − 4.2526 * XYL − 14.2728 * MEK * XYL) showed goodness of fit than other solvent models. Furthermore, the results verifying the model optimization according to physico-chemical properties showed that Q2 values for FP, SG × FP and rVP × FP, and a combination of SG, FP, SG × FP and rVP × FP were 0.5932 and 0.5381, respectively, exceeding the standard for the goodness of fit (Q2 = 0.5). When we compared the experimental value and EC50 prediction value for descriptors affecting their optimization, the linear regression equation (y = 5.7976 − 5.2754 × SG + 6.1237 × FP −7.0023 × SG × FP × 0.02938 × rVP × FP) showed that the EC50 decreased as combinations of BP × SG, BP + FP, SG × FP and rVP × FP increased. These results are most likely due to differences in the toxicity expression, which results from changes in physico-chemical properties via interactions between components in the mixed compounds.

In this study, we used physico-chemical properties and EC50 values as descriptors for QPAR modeling to predict the toxicity of mixed compounds containing EA, IPA, MEK, TOL, XYL, which are used in gravure printing. Our results showed that physico-chemical properties and the EC50 values of mixed compounds were dependent on the physico-chemical properties of single substances. Toxicity predictions from the linear regression equation using experimental values and prediction EC50 values showed that mixtures containing MEK, TOL and XYL displayed a decreased EC50, but in mixtures containing MEK * XYL dramatically decrease EC50 value. Furthermore, either decrease of SG and rVP or increase of FP for solvent mixture increase EC50.