Abstract
Introduction
A comparative toxicological assessment of particulates is being an important subject experimental study within the framework of particle toxicology as a research area in its own right including toxicology of chemically similar metal-containing particles falling within different size ranges, for instance, nanometric particles compared to micrometric ones. At the same time, however, toxic effects of particles of
That is why, the effects of various NP combinations were the focus of a range of comparative experimental and mathematical modeling studies carried out by our research team, including the nanotoxicology of Ni and Mn 5,6 ; Cu, Zn, Pb 7 ; Al, Ti and Si, 8 and Pb and Cd 9 combined in various proportions. Pikula et al., 2020 10 experimented with CdS and ZnS combination.
Alongside this line of research, our aspiration has been to advance the general theory of combined toxicity ( 11 -16 and many others). In these studies, various acute, subchronic or chronic intoxications were modeled by exposing rats to repeated administration of dissolved salts or of nanoparticles of different toxic elements in different combinations. The health effects thus induced were assessed by a wide range of functional, biochemical and morphometric criteria, and the results could be summarized as follows: “effect additivity” and “dose additivity” generally recognized as 2 opposite paradigms of combined toxicology theory are, in fact, interchangeable. Therefore we assumed that these paradigms depend on which particular method is chosen for describing experimental data mathematically rather than on essential differences in fundamental mechanisms of combined toxic action. Moreover, we maintain that the main types of combined toxic action, i.e. additivity, superadditivity (synergism), subadditivity (hidden antagonism) and opposite action (explicit antagonism) are not inherent invariable characteristics of a specific mixture. Indeed, even under exposure to one and the same simple binary combination, these types are revealed depending on the toxic effect considered, the level of this effect, the ratio of the doses involved, and on the presence, if any, of a third, concomitantly acting toxic agent.
Different Cd-Pb combinations are most often present in workplace and ambient air in technologies associated with copper smelting and refining. Both lead and cadmium are toxic for the organism in practically all chemical forms and, thus, are of special research interest. One of the first combinations we assessed for toxicity in a subchronic experiment on rats was one of Cd and Pb in their ionic forms. 17 The data from that experiment were then processed with the help of the Response Surface Method (RSM), 18 which garnered additional support to the above theoretical postulates. This experimental study had been prompted by the results of an epidemiological one 19 that assessed the combined nephrotoxic effect of these 2 metals.
Some time ago we decided to assess experimentally the combined vaso-cardiotoxic effects of lead and cadmium considering the part they both play in human cardiovascular disease. All experimental toxicological studies published thus far involved water solutions of lead and cadmium salts while epidemiological studies considered exposure to lead and/or cadmium without specifying their physical form. In the foregoing, the occupational and environmental hazards of metal-containing nanoparticles for humans were mentioned as important in practical terms, and so the toxic properties of lead and, to a lesser extent, cadmium have been considered in a number of experimental studies, particularly by our team. 7,9,20,21 However, little attention (if any) has been given just to the vaso-cardiotoxicity of such nanoparticles.
As far as we know, the
Data analysis with the help of a response surface linear model with a cross-term provided fresh support to the above postulate that one and the same pair of toxic agents may display a diversity of combined action types.
The present paper focuses on the
Materials and Methods
For the experimental modeling of different nano-intoxications, sufficiently stable suspensions of nanoparticles were produced by laser ablation of respective 99.99% pure targets (in this research, of lead sulfide or cadmium sulfide) under a layer of deionized water. This technique is regularly employed by our team and has been described in sufficient detail more than once (e.g. Katsnelson et al., 1 ). In this way, we have obtained spherical NPs of cadmium and lead sulfides with diameter 37 ± 5 nm and 24 ± 4 nm, respectively)
The
For exposure to nanoparticles, the cells were seeded in a 96-well plate (TPP Techno Plastic Products AG, Trasadingen, Switzerland), 100,000 cells per well in 100 µL nutrient medium. These plates were kept for 24 hours before adding NPs at 37ºC under an atmosphere of 5% CO2 in a DMEM medium with L-glutamine, 1 g/L glucose, 10% fetal bovine serum (FBS) and 0.5% gentamicin antibiotic. The final NP concentrations in the incubation medium (mcg/mL) with concurrently exposed samples were as follows: CdS-NP 3.75-5-7.5; PbS-NP 7.5-11-15; CdS-NP 3.75 + PbS-NP 7.5; CdS-NP 7.5 + PbS-NP 15; CdS-NP 5 + PbS-NP 15; CdS-NP 7.5 + PbS-NP 7.5; CdS-NP 5 + PbS-NP 11 µg/mL. All exposure variants were performed in quadruplicates. Cell incubation in the NP-containing medium lasted 24 hours.
CellTiter-Glo Luminescent Cell Viability Assay (Promega) was used according to the manufacturer’s instruction. The plates were equilibrated to room temperature for approx. 30 minutes and the reconstituted CellTiter-Glo Reagent was added 1:1 to the amount of cell culture medium present in each well. The plates were shaken for 2 min and incubated for additional 10 min at RT before transferring well contents to a luminescence compatible 96-well plate, and luminescence was read on LM-01 T with Kilia software (Immunotech, Beckman Coulter Company, Praha, Czech Republic). The measurement results were expressed in relative luminescence units (RLU).
ET-1 was determined quantitatively in cell lysate using a diagnostic test system by Biomedica (Vienna, Austria). The method is based on enzyme-linked immunosorbent assay. The results were measured in the Epoch Microplate Spectrophotometer (BioTek, USA) at a wavelength of 450 nm and expressed in fmol/ml per 100,000 cells. Cells lysate was obtained by freezing-thawing cell destruction.
For cell size determination, cells were removed from the plate, put on a glass slide and imaged in an optical microscope, 3D Cell Explorer (Nanolive, Switzerland). ImageJ 1.48v image processing program (Wayne Rasband, National Institutes of Health, USA) was used to measure cell and cell-nucleus area, reading the results in µm2.
Data Mathematical Treatment and Analysis
The experimental data demonstrated an explicitly non-monotonic dose-response dependence. Moreover, within the dose ranges used, this dependence is bi- or tri-phasic for some of the indices. This required approximating it by more complex expressions than the usual models of hormesis (see (Panov et al. 2020 9 and articles cited therein). Since polyphase functions describe a complex and multifactorial biological phenomenon, it would be unlikely to expect that such a dose-response dependence could be adequately expressed by simple functions. For instance, in 9 this dependence is represented by power functions depending on 8 parameters.
In the present paper, a different approximation approach was employed, in which the complexity of the dose-response relationship is hidden within the approximation functions while the form of such approximation itself is very simple: it is a linear combination of the chosen functions. In particular, the template used for approximating the multiphase dose-response relationship was given by the sum:
where
Note that the proposed form of dose-response approximation is effective exactly where a polyphase hormesis is present. On the contrary, for classical hormesis or a usual monotonic dose-response relationship it would be sufficient to apply simple analytical expressions as used below.
As well as in our previous works on combined toxicity, 5,7,8,11,13,18,12,22 the mathematical modeling of Cd-Pb combined action was performed by the Response Surface Methodology (RSM).
The regression equation describing the response surface Y = Y (x1, x2) in our case is:
where Y is a toxicity index, x1 and x2 are the doses of the agents participating in the combination (cadmium and lead, respectively). This equation was constructed by fitting the coefficients
Results and Discussion
Like in the study mentioned in the Introduction, 9 performed with the same NP species on an HL-1 culture of cardiomyocytes, the dose-response function for the culture of endothelial cells often displayed a non-monotonic character. The reader is referred to the above publication for the references to works from which we derived the mathematical models used as such or modified by us for satisfactory approximation of this kind of dose-response relationships. In the same recent publication, we discussed specially the understanding of such non-monotonic dose-response functions within the hormesis paradigm in the broad sense of the term. In the current paper, therefore, we do not repeat this discussion, and confine ourselves to graphic illustration of the relevant results and their model approximation, as well as brief comments relating to the meaning of specific responses.
Figures 1 and 2 display a rather similar pattern of diminishing culture viability as the dosage of CdS-NP or PbS-NP is increased, although in different dose ranges. Note, however, that in our previous experiment on HL-1 cardiomyocytesCCЫЛКА for this integral index of cytotoxicity, a similar type of dose-response dependence was observed for PbS-NP only, whereas for CdS-NP it was monotonic and was described well, for instance, by the hyperbolic function.

Reduction in ATP-dependent luminescence in endothelium cell culture under exposure to various doses of CdS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mcg/mL; the ordinate plots the values of the Relative Light Units (RLU) index in the exposed or control groups. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is

Reduction in ATP-dependent luminescence in endothelium cell culture under exposure to various doses of PbS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mg/mL; the ordinate plots the values of the Relative Light Units (RLU) index in the exposed or control groups. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is
As seen from Figures 3 and 4, a dependence of the endothelial cell area on the dose of the same factors is, on the contrary, triphasic for CdS-NP and monotonically hyperbolic for PbS-NP.

Changes in endothelium cell area under exposure to various doses of CdS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mg/mL; the ordinate plots the cell area in mcm2 in the exposed or control group. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is

Changes in endothelium cell area under exposure to various doses of PbS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mcg/mL; the ordinate plots the cell area in mcm2 in the exposed or control group. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is
It is to be emphasized that even when cell size increased somewhat with increasing concentration of CdS-NP in the second phase of the dose-response relationship, it still did not differ significantly from the control value. Statistically highly significant was only the impact of the CdS-NP dose that induced an opposite effect, while this hypotrophic effect was characteristic of all PbS-NP doses tested. Note that in our animal experiments
Taking into consideration results of an animal experiment presented and discussed by Sutunkova et al.,
24
of special interest is the question whether or not the

Changes in endothelin-1 concentration in cell lysate under exposure to various doses of CdS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mcg/mL; the ordinate plots the endothelin-1 concentration, fmol/mL in the exposed or control group. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is

Changes in endothelin-1 concentration in cell lysate under exposure to various doses of PbS-NP. The abscissa plots the acting concentrations of nanoparticles in the incubation medium, mcg/mL; the ordinate plots the endothelin concentration, fmol/mL in the exposed or control group. The dots indicate the mean values with the standard error of the mean. The equation for the approximating curve is
As can be seen from the isobolograms constructed on the basis of the RSM model (Figure 7), the combined toxicity of CdS-NP and PbS-NP is unidirectional and, at the same time, additive for ATP-dependent luminescence reduction but contra-directional for the other 2 effects considered above. Such variability of combined binary action typology, depending, among other things, on which effect it is estimated by, pertains to the important postulates of the theory developed by us as was emphasized in the Introduction. In addition to the ample support for this postulate obtained in

Isobolograms characterizing the combined toxic action of CdS-NP and PbS-NP on an endothelial cell culture as estimated by its effects on (A) ATP-dependent luminescence, (B) cell area and (C) endothelin-1 concentration. Numbers at the axes are respective NP concentrations in mcg/mL; numbers at the isoboles are the values of the effect to which they correspond.
At the same time, these studies (Minigalieva et al., 2017 6 in particular) showed that a specific combined toxicity type characteristic of a particulate’s effect on a cell culture and of a seemingly toxicodynamically-related index may mismatch a similar index on organ/system level. The same may be noted when comparing the result shown in Figure 7 with the effect of combined lead-cadmium subchronic intoxication on ET-1 level in blood serum 24 ).
These inconsistencies are hardly surprising considering the fact that the outcome of an effect of just 2 toxic agents on any index characterizing the status of a whole animal organism is mediated by a considerably greater number of primary and feedback mechanisms than the outcome of an impact on a cell pulled out from this mechanistic framework (not to mention micro-environmental differences).
Therefore, as early as when discussing the comparative results of assessing the NiO-NP + Mn3O4-NP combination on several cell cultures and on rats
6
we emphasized that we consider
Conclusions
It may be assumed that the described by us separately adverse vascular-regulatory effects of animal intoxications with lead and cadmium are mediated by various toxicodynamic mechanisms operating at system and organism levels rather than being connected with the direct cytotoxic action of lead- and cadmium-containing nanoparticles on cells controlling the production of endothelin-1. This assumption is suggested by the fact that in the present experiments in which an established line of arterial endothelial cells was subjected to such direct action, the inhibition of ET-1 production was observed in a certain range of doses for cadmium-containing NPs only, and not for lead-containing ones at all. At the same time, the dose-dependent adverse effect on the viability and size of these cells was more definite for the action of the lead-containing nanoparticles.
It has been confirmed once more that the combined binary action typology is essentially variable depending on which effect it is estimated by. However, a specific combined toxicity type characteristic of a particular effect on a cell culture and of a related index of
