Spectral-Structure Activity Relationship (Spectral-SAR) Assessment of Ionic Liquids’ in Silico Ecotoxicity

The use of green solvents, including ionic liquids (IL), in the synthesis of new materials is currently highlighted worldwide [1-5]. The green high thermal stabilities, and can remain in the liquid state over a wide range of temperatures (from -50°C to 200°C). Furthermore, ILs exhibit low toxicity because of their exceptionally low volatility. This property is why ILs could be used as “green” alternatives to volatile organic solvents in many processes [6,7].


Introduction
The use of green solvents, including ionic liquids (IL), in the synthesis of new materials is currently highlighted worldwide [1][2][3][4][5]. The green high thermal stabilities, and can remain in the liquid state over a wide range of temperatures (from -50°C to 200°C). Furthermore, ILs exhibit low toxicity because of their exceptionally low volatility. This property is why ILs could be used as "green" alternatives to volatile organic solvents in many processes [6,7]. Accordingly, ILs have been used as a reaction medium for inorganic materials, which demonstrates the pre-organized structure of the ILs as a template for generating porous nanomaterials [8,9]. It is well-known that the properties and structures of mesoporous materials vary by the choice of synthesis conditions, such as temperature, reagents, base concentrations, and the nature of the organic groups [10]. By organic groups, one means the organic chains and structures that affect the functionality of the precursor; for example, where there is Si-C bond, depending on the silica to surfactant ratio, different pore topologies (hexagonal, cubic and lamellar) are present, and when only ionic liquid was used as the surfactant, mesoporous silica with an irregular shape may be obtained [11]. However, when a small amount of cetyl trimethylammonium bromide (CTAB) was added to the reaction, mesoporous silica with a spherical shape was obtained. When a mixture of ionic liquids and CTAB was used as the surfactant templates, larger mesopores were formed in the silica spheres [12]; however, ILs can generally dramatically decrease the polarity and dielectric constant of water and increase the solubility of surfactants when used as co-solvents [13,14]. However, the interaction between CTAB and IL has also been observed, which causes the IL to behave as a co-surfactant [14]. Starting from these new concepts, ILs with different cations (e.g., imidazolium, pyridinium and ammonium) and anions (e.g., hexafluorophosphate, tetrafluoro-borate, octyl sulfate and bromide) and the ordinary CTAB and the mixture of CTAB and ILs are used in the synthesis of mesoporous silica.
Furthermore, ILs provide an interesting method for controlling the chemical contents of the collapsing bubble through the reduction of the solvent vapor pressure [15] because the rates of sonochemical reactions can be increased by decreasing the vapor pressure of the solvent or by choosing a less volatile one [16]. ILs are usually more viscous and denser than other organic solvents, and producing cavitations should be more difficult under such conditions with large cohesive forces [16,17].
Studies of the thermal decomposition of ILs confined in nanopores compared to bulky pores are also currently on the forefront of nanomaterials research. Many physicochemical properties of ILs have been observed to change when confined in nanopores, which results from the interaction of the IL cations/anions with the pore walls, which could be weak, as in silica nanopores [18,19]. These studies influence the next step, i.e., the surface functionalization of mesoporous silica materials through the covalent bonding of organic groups, which is achieved through the comparative co-condensation and post-grafting of amino and alkyl functionalized silica alkoxide precursors [20]. Importantly, the comparative study is necessary to optimize the synthesis methods for nanomaterials because, even if the post-grafting method produces wellordered functionalized mesostructured materials, it often produces non-uniformly distributed organic groups [20]. Furthermore, the co-condensation synthetic method for mesoporous materials involves a one-step procedure and allows better control of the loading and distribution of the organic groups [20], as was detailed in the recent review on the advances in (bio)responsive nanomaterials [21]. Supplementary information on the effects of ILs on the texture of gels can be obtained from the pore size distribution. This aspect is interesting because the ionic liquid concentration changes the pore volume due to the low volatility of the IL, which reduces the surface tension associated with the pore collapse [22][23][24][25].
The next step in the ecotoxicological evaluation of the synthetized mesoporous silica is to use the logistic enzyme kinetic and Spectral-SAR methods [26][27][28][29][30] with acetylcholine esterase as a working enzyme. Acetylcholine is a notably good marker for monitoring the evolution in ionic liquids, due to the similar chemical structures of many ionic liquids cations and acetylcholine and especially because this enzyme does not follow the standard Michaelis-Menten kinetics mechanism. The logistic enzyme kinetics [31][32][33][34][35] is an analytical substitute for the Michaelis-Menten mechanism, and it is based on some probabilistic considerations that better accommodate the QSSA (quasi-steady state approximation) conditions, while providing the analytical transformation of the W-Lambert function into a logarithmic function; this so-called logistic transformation was successfully evaluated for different cases of enzyme kinetics [31][32][33][34][35][36][37][38][39]. The analytical information contained in the logistic enzymatic temporal curves could be further combined with a QSAR algorithm, such as the Spectral-SAR, for designing molecular mechanisms of the ecotoxicity; this procedure has already been applied for bare ILs [27,28,38,40]. Eventually, this procedure is intended to be used for some toxicity tests [41,42] or even for obtaining the efficiency of the material as a drug carrier in biological media, if the proper conditions are satisfied.
Regarding the action/toxicity mechanism of ionic liquids (M.O.A.), possible mechanisms of toxic actions are identified through membrane disruption because they have structural similarity with detergents, pesticides and antibiotics [43][44][45][46] that induce polar narcosis, due to their interfacial properties, and may cause membrane-bond protein disruption [47][48][49]. As a possible mode of action, disrupted membranes and hydrophobic molecules have a greater ability to accumulate at this interface [50]. The low K ow values of imidazolium ionic liquids indicate the low permeability of these ionic liquids [47]. However, other mechanisms arise from acethylcolinesterase inhibition [51], relate to a common cellular structure or process [45], involve structural DNA damage [47,52], or are mediated by certain bacteria that could potentially break down imidazolium into different metabolites that could be inhibited by the ionic liquid itself [43,44]. Nevertheless, knowing if the chemicals could reach the target site in the organism and by which mechanism is of paramount importance, and at minimum, a conceptual extrapolation from in vitro to in vivo studies is required [51] because other factors, such as biodegradation and bioaccumulation, are necessary before conclusions can be drawn [49]. The dose-response approach for estimating the lethal effects of toxicants on organisms is being criticized because of the lack of real ecological meaning; furthermore, regulatory norms have been established around LC50 values that can be compared between different toxicants and organisms [47,53]. In this respect, it was observed that the side chain toxic mode of action of ionic liquids is not specific to species, and a similar relationship has been observed [47]. However, the systematic variation of, for instance, R1 and R2 at identical head groups and anions from the molecular to individual organism level lead to the conclusion in most published data that the shorter the chain lengths of the side chains, the lower the cytotoxicity (higher EC50 values) [50,54]. The electronic portion of the bond factor extends along a chain that has no more than 5 alkane carbons, which results in a decrease in the overall hydrophobicity; therefore, chains longer than 5-6 carbon atoms will have a considerably greater permeability through the cell membrane (for example, [DMIM][BF4]), see ref. [55]. Therefore, chemical transformations of the side chains of ionic liquids may reduce their toxicity because the metabolites are less toxic compared to their parent chemicals [54]. With respect to the cationic effects, see Figure 1, the imidazolium ring plays a major role due its structure as a delocalized aromatic system with a high electron acceptor potential; therefore, the nitrogen atoms are not capable of forming any hydrogen bonds, which results in a rigid and sterically hindered system. The elongation of the R2 residue causes a continuous increase of flexibility, which implies more conformational freedom [56]. The chemical nature of a head group also influences the biological activity [54]. The fragmental hydrophobicity of each carbon connected to a quaternary amine combines a geometric bond factor that applies to the neutral solute together with a negative electronic bond factor, which decreases in magnitude with the square of the distance from the central nitrogen atom [55,57]. Complementarily, the anion could play a central role as a technicophore because it exhibits a high potential for changing the materials technological properties (solubility, viscosity); however, there are reports that the anion presents no significant effects [44,45,55]. The first study reporting the influence of the anion was by the Jastorff group [54], where it was reported that the anion in the ion pair can partially decompose, as was previously observed [55,58,59]. Finally, at the computational level, the detailed examination of the relative energies and structural interactions (such as ion position, H-bonding and variability in the anion conformation) in gas-phase ion-pairs has emphasized how these quantities can be used to construct a picture of the local structure and interactions that occur in ionic liquids. For instance, [BMIM][Cl] (1-n-butyl-3-methylimidazolium chloride) forms a highly connected liquid with relatively strong interactions, [BMIM][TFMSi] (1-n-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) forms a low connectivity network of weakly linked ions, whereas [BMIM][BF4] lies between these extremes and forms a weak but more regular network. The melting points and viscosity are somewhat dependent on the local interactions between an ion and other molecules in the first solvation shell. For imidazolium based cations, nine sites of interactions are preferred by the anion [60][61][62][63]. Hunt et al. demonstrated that the hydrogen bond is primarily ionic with a moderate covalent character. The fact that the Coulombic attraction is the dominant stabilization force was demonstrated by analyzing the charge distribution, molecular orbitals and electron density of the dimer complex [BMIM][Cl] [63]. The interactions governing the top conformers are notably different compared to those where the Cl anion remains in-plane (where the Cl anion interacts with the π-manifold of the orbitals); the LUMO of the ion-pair is the cation antibonding orbital; therefore, electron acceptance is not favorable [63]. The effect of the chloride anion on the rotation of the butyl chain was consequently investigated and observed to lower certain rotational barriers while enhancing others [61]. The Coulombic forces in an ionic liquid, and hence the point charges or charge distribution on the constituent ions are anticipated to be more important than in liquids composed of neutrally charged molecules [60]; these ions allow a type of semi-classical analysis for the toxicity of ILs, which is well-adapted for quantitative structure-activity studies (QSAR), to which the present review belongs. Nevertheless, even the well-regarded QSAR methodology has to be refined and adapted or extended to include the chemical electrostatic interactions present in the IL and of their transduction into various organism triggering (or not) toxicity. These aspects will be reviewed and exemplified in the following.

QSAR/Spectral-SAR Modeling of the Structure-EcoToxicity of Ionic Liquids
It is an already established fact that the costs of all approaches for sustainable product design can be reduced using the SAR and QSAR methods [56]. For instance, the anti-microbial activity of quaternary ammonium chlorides is dependent on the lipophilicity, and the 1-octanol-water partition coefficient, K ow , could provide a good approximation for the lipophilicity of the compound. The relationship between the chemical structure and anti-microbial activity of several new choline-like quaternary ammonium derivatives was analyzed using the QSAR method, and all of the studied compounds were active against the microorganisms. The statistically significant correlation of the obtained QSAR equations confirms that the lipophilicity is the primary factor that governs anti-microbial activity [43]. The octanolwater partition coefficient has been correlated with bioaccumulation and toxicity in fish and to sorption in soil and sediments. The K ow values for the BMIM cation range from 0.003 to 1.1, depending on the choice of anion, and the K ow value increases with increasing length of the alkyl chain on the cation. Replacing the H atom between the 2 nitrogen atoms in the ring with a methyl group has almost no effect on the K ow value. Because all of the measured K ow values in the analyzed ionic liquids, i.e., imidazolium ionic liquids with butyl, ethyl, dimethyl, hexyl and octyl side chains with different anions, such as tetrafluoroborate, nitrate, bis(trifluoromethylsulfonyl)imide), are less than 15, these ionic liquids will not accumulate or concentrate in the environment. It was also observed that the K ow values for ionic liquids with the TFMSi (bis(trifluoromethylsulfonyl)imide) anion (a notably hydrophobic anion) are dependent on concentration, even at the diluted concentrations studied [64].
However, QSPRs (Quantitative Structure Properties Relationships) methods could be used to correlate and predict the toxicity of ionic liquids. With at most 4 molecular descriptors, log EC50 and log EC5 data are reproduced with a R 2 (statistical Pearson correlation coefficient) of 0.78-0.88 using tests with Vibrio fischeri and Daphnia magna. These methods calculated the electronic, spatial, structural, thermodynamic and topological descriptors for both the cation and anion separately, whereas the geometries were also optimized separately. Concerning Vibrio fischeri, the QSPR equations that used cationic (first equation) and cationic + anionic (second equations) descriptors fit the training set well, with R 2 =0.887 and R 2 =0.782, respectively. However, the majority of the training set compounds were imidazolium and pyridinium based ionic liquids. Then, Couling and co-workers [49] used the equation to predict the toxicity of new compounds, which contained ammonium and phosphonium cations, but the toxicity of these new additional compounds was not predicted with the same accuracy, which led to the conclusion that the equation is only good for compounds of the same class. The best accuracy for a new predicted compound (using the first QSPR equation) was for [BdMAPy] [Br] (1-n-butyl-4-dimethylaminopyridinium bromide) because it contains an aromatic pyridinium cation, and it also contains an amino group, unlike all of the other training set compounds. Concerning Daphnia magna, the QSPR equation fits the training data set well (R 2 =0.862), and the descriptors are similar to those provided for Vibrio fischeri, which suggests that there may be similar indicators of toxicity observed in many different species. The equation was used to predict the toxicity of another compound, and contrary to the case of Vibrio fischeri, the correlation was good, R 2 =0.775; however, the neglected compounds were more similar to those of the training set compounds [49]. Such results further motivated the use of the actual Spectral-SAR method for these two paradigmatic species response upon IL action because it is an important step in constructing a consistent ecotoxicological test battery. The reliability of the Spectral-SAR/QSAR studies primarily relies on the lack of readily accessible ecotoxicological data [65]. Note that because the principles of green chemistry state that one should consider the entire process (life cycle analysis), rather than the individual components of a reaction (single issue sustainability), the acute toxicity measurements do not fully characterize the full impact of the release of a substance into the environment, but they are only part of the environmental impact assessment [65,66]. In particular, ionic liquids with cations, such as pyridinium, imidazolium and pyrrolidinium, have been nominalized by the United States National Toxicology Program (NTP) for toxicological testing because of their potential use as new solvents and potential to enter the aquatic system [67]: if an accidental discharge of ionic liquids into water were to occur, there may be an environmental risk to aquatic plants and animals because many of the ILs are water-soluble [48]. The recorded averaged descriptions of the acute toxicity based on the LC50 (mg/mL) are very highly toxic (<10 -4 ), highly toxic (10 -4 -10 -3 ), moderately toxic (10 -3 -10 -2 ), slightly toxic (10 -2 -10 -1 ) and not acutely toxic (>10 -1 ) [55,68].

Algebraic QSAR: Spectral-SAR [69,70]
The key concept in the SAR discussion concerns the algebraic consideration of biological activity and the structural parameters in Table 1. As a consequence, we may further employ this feature to quantify the basic SAR through an orthogonal space. The idea is to transform the columns of structural data in Table 1 into an abstract orthogonal space, where all of the predictor variables are independent, solve the SAR problem in orthogonal space and subsequently compare the result to the initial data using a coordinate transformation. Because QSAR models aim to develop correlations between the molecular structures of interest and the measured (or otherwise evaluated) activity, it naturally appears that the structure part of the problem is accommodated within quantum theory and its formalisms.

Activity
Structural predictor variables In fact, there are a few quantum characteristics that we are using within the present approach [29]: • Any molecular structural state (dynamical because it undergoes interactions with the species and organisms) may be represented by a | ket state vector in an abstract Hilbert space, which follows the bra | ket Dirac formalism [71]; such states are to be represented using any reliable molecular index, or specifically, in our study, by the hydrophobicity | LogP , polarizability | POL , and the total optimized energy | E tot . These parameters are only the so-called Hansch parameters, which are usually employed to account for the diffusion, electrostatic and steric effects for molecules acting, for instance, within the cells of organisms, respectively; • The (quantum) superposition principle, which ensures that the summation of molecular states map onto another resulting molecular state, which is interpreted here as the bio-, eco-or toxico-logical activity, e.g., resents the "noise" activity (present even when all other influences are absent); • The orthogonalization feature of quantum states, which is a crucial condition where the superimposed molecular states generate other molecular states (here quantified as the activity of the molecular ligand -linking receptor); analytically, the orthogonalization condition is represented by the bra | ket scalar product of two envisaged states (molecular indices). If bra | ket = 0 , then the states are said to be orthogonal, and the molecular descriptors are independent; therefore, they are suitable to be added as contributing states in the resulting activity and as molecular indices in the activity correlation.
As such, the analytical procedure is decomposed into three fundamental steps.

I.
Given a set of N molecules that are being examined for their biological activity, they produce through the considered M structural indicators all of the input information (the states) that may be vectorially expressed by the columns in Table 1 and correlated through the following equation: added to account for the noise term. For equation (1) to represent a reliable model of the given activities, the assumed molecular states (indices) should constitute an orthogonal set; having this constraint is a quantum mechanical fundamental, as described above. However, unlike other important studies that have addressed this problem [72,73], the present use of the Spectral-SAR method assumes that the prediction error vector in eq. (1) arises from the predicted activity's being orthogonal because it cannot consider the input data otherwise being not known a priori any correlation is made. In this manner, it follows from eqs. (1) and (3) that the prediction error vector has to be orthogonal on all other descriptor states of the predicted activity.
In other words, conditions (3) and (4) confirm the form of (1) in the sense that the prediction vector and the prediction activity |Y PRED (with all of its sub-intended states | X i=0,M ) belong to disjointed (thus orthogonal) Hilbert spaces; that is, one can state that the Hilbert space of the observed activity |Y OBS may be decomposed into a predicted and an error independent Hilbert sub-spaces of states. Therefore, within the Spectral-SAR procedure, the first step in the orthogonalization procedure orthogonalizes the predicted activity to its prediction error, whereas the remaining orthogonalization algorithm does not search for optimizing the minimization of errors, but it searches for the optimum method for producing the ideal correlation between |Y PRED and the given descriptors | X i=0,M .

II.
Next, the Gram-Schmidt orthogonalization algorithm is applied by constructing the orthogonal set of descriptors using the consecrated iteration [74][75][76] 0 0 providing the orthogonal correlation

III.
Remarkably, while the studies dedicated to the orthogonal problem usually stop at this stage, the Spectral-SAR method uses this stage to provide the solution for the original searched correlation, eq. (1). This process can be adequately achieved by rearranging eqs. (6) and (8) such that the system of all descriptors in Table 1 can be written in terms of orthogonal descriptors ...
If the determinant of eq. (11) is expanded on its first column, and the result rearranged so that to have |Y PRED on the left side and the rest of states/indicators on the right side, the searched QSAR solution of the initial problem of eq. (1) is obtained as the Spectral-SAR vectorial expansion (from where the "spectral" name is also justified) with the error vector already absorbed in the orthogonalization procedure. In fact, the Spectral-SAR procedure uses double conversion passages: one forward from the given problem of eq. (1) to the orthogonal one of eq. (8) where the error vector is orthogonally "dissolved"; and the reverse one, back from the orthogonal to the real descriptors throughout the system (10), which provides the determinant (11) to be expanded as the QSAR solution.
The result is that the QSAR/Spectral-SAR equation is now directly delivered by the determinant (11) and not through matrices products, as in the statistical Pearson approach, while directly providing the Spectral-SAR correlation equation and not only the parameters of multi-variate correlation [77][78][79][80][81][82][83]. Furthermore, the Spectral-SAR algorithm is also invariant to the order of descriptors that are chosen in orthogonalization procedure, which provides equivalent determinants only with rearranged lines; this is a matter that was not previously achieved by other orthogonalization techniques [84][85][86][87].
However, it is worth being convinced by comparing the present Spectral-SAR method with the standard statistical one by specializing the general problem (1) to the linear case and to determine whether it furnishes the linear regression parameters given by the consecrated least squares analysis through the Spectra-SAR equation (11) [81][82][83]. In this respect, we actually deal with the particular equation which is immediately rearranged to ( ) such that identifying the actual linear coefficients When evaluating expressions (15) and (16) within the Spectral-SAR algorithm, there are instructions to identify only the relevant variables from Table 1 using the convenient notation Other working tools are the zero-th and the first orthogonal vectors, which are accordingly considered and computed as (5) with (2) and respectively by.
with the help of the coefficient (18) specialized from the general definition (7).
In the same manner, the other specific Spectral-SAR coefficients from the general orthogonal recipe (9) are now computed as the zero-th order contribution for linear regression (19) whereas the first order contribution precisely recovers the linear slope [69,70] as prescribed by the correspondence of (15). Additionally, its companion -the noise term coefficient of eq. (14) -may be now directly evaluated with (16) which successfully regains the noise term that is otherwise consecrated by means of a variational statistical (optimization of errors' squares summation) procedure.
With this description, it is clear that the SPECTRAL-SAR methodology not only recovers the standard statistical QSAR least square correlation results but also generalizes it analytically to a great extent toward a better assessment of mechanistically ordering and influences in practical eco-and bio-logical applications.

From in Cerebro to in Silico Principles of Ecotoxicity
When utilizing the analytical model of QSAR/Spectral-SAR for environmental interactions, one should consider the framework of the principles both at the general and applied levels.
As such, regarding the general principles or green chemistry and engineering, they are provided in Table 2 to ensure that they can be readily compared [88]. Note that the fundamental principles constitute the background or the general framework or desiderate that is eventually supported by the associate engineering principle; the green engineering principle is primarily based on "minimizing" or "maximizing" the time, space, energy and costs, and it is constituted either in an economical enterprise and an extension of the main principles of nature or in terms of optimizing mass-energy and time-space. From this point forward, the most basic physical and chemical principles are observed as those acting on each process, system, or state to be created, maintained or modified [89][90][91].
However, while restraining the analysis to the specific interactions between chemical structures and biological species, the Organization for Economic Co-operation and Development (OECD) advanced a set of standard principles for the validation and for regulatory purposes of the (quantitative) structure-activity relationship models [92][93][94][95]: Within this context, the present QSAR-Spectra-SAR (QSAR-SSAR) approach "responds" to these OECD-QSAR principles by the present Spectral-SAR ecotoxicological principles' realization [96], with special reference to ionic liquids (Principle 3): • ECOTOX-SSAR Principle 1 is assured by: the "length" of the predicted/measured (eco)biological action follows the self-scalar product rule of the computed endpoint activity ( ) • ECOTOX-SSAR Principle 2 is assured by: the "orthogonality" of assumed molecular factors that correlate with eco-and bio-effects is assured by the spectral decomposition of the associate activity respecting them (see eq. (8)), and the orthogonal-real space transformation given by the Spectral-determinant, eq. (11), giving in fact the searched structure-activity relationship model; • ECOTOX-SSAR Principle 3 is assured by: the method of considering the structural parameters and the activities with which they should be correlated, for specific (target) class of molecules; for instance, the consecrated Hansch quantitative structure activity relationships (QSARs) generally prescribes the activity expansion under the generic minimal but meaningful form: which provides sufficient information about the transport, electronic affinity and specific interaction at the molecular level, respectively; whereas the hydrophobicity index, LogP, describes, at best, the quality of molecular transport through cellular membranes. For the electronic and steric contributions, many structural parameters may be considered [97,98]; among them, the polarizability (POL) measures the inductive electronic effect that reflects the long range or van der Waals bonding, whereas for the steric component, the total energy (E TOT ) is assumed to be the representative index because it is calculated at the optimum molecular geometry at which the stereo-specificity is included. These parameters have been demonstrated to be quite reliable in modeling the ecotoxicological interactions [26,29,30], and they will also be used in the present IL applications. To this aim, when information about the eco-biological influence of ionic liquids (IL) is desired, the particular anionic-cationic structure has to be properly considered because almost all structural information about ionic liquids is based on the superposition of the separate anionic and cationic contributions. In this situation, two different additive models for modeling anionic-cationic interaction can be considered. The first model is based on the vectorial summation of the produced anionic and cationic biological effects. In other words, this so-called |1+> model is constructed from the superposition of the anionic (subscripted with A) and cationic (subscripted with C) activities, and can be formally represented as [27,30]: with Hansch combinations , , Practically, with the |1+> model, the SPECTRAL-SAR procedure is separately performed for the anionic and cationic subsystems, and it is subsequently summed in the resulting IL-activity. The second SPECTRAL-SAR model can be advanced here when the additive stage is considered at the incipient stage of the SPECTRAL-SAR operator (1) such that the considered Hansch factors, for instance, are first combined to produce the anionic-cationic (subscripted with AC) indices that are further used to produce the spectral mechanistic map of the concerned interaction, producing the so-called |0+> model [28,30]: with the Hansch specification of the spectral vectors: The open issue addresses whether the |0+> & |1+> states leave with the same results or what aspects of the SPECTRAL-SAR operator (1) might differ in the IL ecotoxicity, a matter that is solved by computing the so-called ionic liquid internal angle between the anion-cationic activity vectors, Figure 2, with y iA , y iC , i = 1, N components following the prescription [27,28,30, • ECOTOX-SSAR Principle 4 is assured by the "intensity" of the chemical-eco-bio-interaction that is determined by the ratio of the expected to measured activity norms [26,29,69] as a counterpart of the classical statistical correlation factor [69,98,99] ( ) • ECOTOX-SSAR Principle 5 is assured by the "selection" of the manifested chemical-eco-(bio-)binding that parallels the minimum distances of paths [26,96,97] [ ] connecting all possible endpoints in the norm-correlation hyperspace In this manner, the "validation" of the obtained mechanistic picture is achieved by requiring that the influential minimum paths are numbered by the cardinal of the input structural fac-tors set such that, excepting that the final endpoint that is always considered as the final evolution target, all other endpoints are activated one time and one time only.
The present chapter reviews the presented Spectral-SAR-IL models as applied to studying the ecotoxicity of the aquatic species Vibrio fischeri and Daphnia magna by the tested ionic liquids, which were appropriately chosen such that they contained a wide variety of heads, side chains, and anions.

Two Cases of Ionic Liquids Ecotoxicity
Long and short term tests are used to explore the toxicity of the ionic liquids [100][101][102].  Table 3. The series of ionic liquids whose toxic activities A= Log(EC 50 ) on Vibrio fischeri were considered [49], with the marked values being taken from [45] along with the structural parameters LogP, POL (Å 3 ), and E TOT (kcal/mol) that account for the hydrophobicity, electronic (polarizability) and steric (total energy at optimized 3D geometry) effects, computed with the help of the HyperChem program [103], for each cation and anion containing ionic liquid, respectively [27].

NAME
The results from using the Ames test for mutagenity with the Salmonella typhyimurium species indicated that none of the imidazolium, pyridinium and quaternary ammonium ionic liquids caused mutations. However, the Ames test is a short-term test, and it cannot be possible to fully predict the carcinogenicity in animals [46]. The antimicrobial activity (against strains of Grampositive and -negative bacteria and fungi) increase with increasing alkyl chain lengths in pyridinium, imidazolium and quaternary ammonium salts [43]. However, because bacteria have a short generation time, they are an ideal starting point for investigating the structure-activity relationship in ionic liquids and serve as a basis for further toxicity tests to higher organisms and more complex systems; a few examples are listed below [45,67]: (1-n-butyl-3-methylimidayolium bromide) was observed to be considerably less toxic than [BMPy] [Br] (1-n-butyl-3-methylpyridinium bromide); • Cation toxicity increases with the increase in the hydrophobicity of the molecules, from Bpy (1-n-butyl-3-pyridinium) to BMPy (1-n-butyl-3-methylpyridinium) to BdMPy (1-nbutyl-3-dimethylpyridinium); • Butyl substituted ionic liquids have considerably lower logP values and are more water soluble; • In comparison with some commonly used industrial solvents (acetone, methanol, ethyl acetate, and dichloromethane), the octyl and hexyl substituted ionic liquids are more toxic.

• [BMIM][Cl] (1-n-butyl-3-methylimidayolium bis(trifluoromethylsulfonyl)imide chloride) and [BPy][Cl]
(1-n-butyl-3-pyridinium chloride) may have similar activity to the pesticide parquet.   Table 4. Spectral structure activity relationships (S-SAR) of the ionic liquids in Table 3 against their toxicity to Vibrio fischeri, and the associated computed spectral norms with ║|Y EXP >║=4.41537, statistical and algebraic correlation factors, computed upon the relations (22), (31), and (32), throughout the possible correlation models considered from the anionic, cationic and composed ionic liquids data from  Table 5. The variation of the cosines of the anion-cationic correlation angle in vectorial space based on Eq. (30) for all considered modes of action of the ionic liquids in Table 3 with the cationic and anionic subsystems S-SAR predicted activities, respectively [27].  Table 6. Synopsis of the statistic and algebraic values of paths connecting the S-SAR models of Table 4 in the normcorrelation spectral-space of eq. (34) for the ionic liquids in Table 3 against the toxicity to Vibrio fischeri. The primary, secondary and tertiary -the so called alpha (α), beta (β) and gamma (γ) paths -are indicated according to the least path principle in spectral norm-correlation space with the statistical and algebraic variants of the correlation factors used, respectively [27].

Path
The Microtox Acute Toxicity Test is often used to determine the toxicity of single compounds in sediment contamination studies and for monitoring industrial effluents in environmental water quality surveys [49,100]. Before releasing the ionic liquids into the environment, the antimicrobial properties using Vibrio fischeri in the Microtox method were thoroughly examined [45]. Equally, by other tests with Vibrio fisheri, an increase in the toxicity corresponding to an increase in the chain length and an increase in the number of alkyl groups substituted on the cation ring was also demonstrated [45,101].
Salts used for anion substitution, such as sodium bromide and sodium dicyanoamide, were less toxic to Vibrio fisheri than the compounds used for synthesizing the cation. For instance, 3-methyl pyridine and 1-methylimidazole may be used as starting compounds. The addition of a butyl chain to the C1 (carbon atom) of the pyridinium and imidazolium cation slightly increased the toxicity of the ionic liquid. Furthermore, the addition of a hexyl or octyl chain will also increase the toxicity. The same observation holds for quaternary ammonium salts. All of these observations suggest that the effects are related to the lipophilicity of the cation, thus explaining why many of the imidazolium and pyridinium ionic liquids were more tox-ic than acetonitrile, acetone and methanol. In contrast, monatomic anions (such as bromide, chloride) are predicted to be less toxic than large anions that contain regions of positive charge, e.g., TFMSi (bis(trifluoromethylsulfonyl)imide) [49]. Figure 3. The spectral representation of the chemical-biological interaction paths across the S-SAR to the modeled endpoints of the Vibrio fischeri, according to the least (shortest) path rule within the spectral norm-correlation space applied to the data in Table 6 for the cationic, anionic and resulting ionic liquid norms of Table 4 for the statistical and algebraic versions of correlation factors from up to down and left to right, respectively. The primary-alpha, secondarybeta and tertiary-gamma path hierarchies of Table 6 are indicated by decreasing thicknesses of the connecting lines [27].
At the organism level, which are Vibrio fisheri and Daphnia magna in this review, the toxicity generally increases with the cation type as ammonium<pyridinium<imidazolium<triazoli-um<tetrazolium. The choline ionic liquids (with a negatively charged oxygen atom) and the quaternary ammonium ionic liquids with short chain lengths were relatively non-toxic to Vibrio fisheri. Furthermore, the methylation of the aromatic ring of the cation reduces the toxicity to Vibrio fisheri and Daphnia magna [49]. The toxicity increases with the number of nitrogen atoms, whereas the anions play a secondary role in the toxicity (even if the presence of positively charged atoms on the anion is predicted to slightly increase the toxicity) [49].
Note that Daphnia is an important link between the microbial and higher trophic levels [49,102]. By acting on Daphnia magna, ionic liquids with longer alkyl chain substituents have toxicities comparable to phenol, whereas those with shorter substituents (for example, [BMIM][Br]) were more toxic to Daphnia magna than benzene and methanol [47,49].
The lethal concentration to D. magna (also called the acute toxicity) was observed to be considerably lower for ILs that employed imidazolium as cation than salts with Na + ; therefore, the toxicity was related to the imidazolium cation and not to the anions. With respect to the life history of the freshwater crustacean, salts with a sodium cation (NaPF 6 sodium hexafluorophosphate and NaBF 4 , sodium tetrafluoroborate) affect the reproduction of D. magna at high concentrations. Sub-lethal effects on the life-history traits occurred when the concentrations of the ionic liquids were an order of magnitude lower than those for acute effects. The clonally variation was a possible reason why the average brood size in controls with those for ionic liquids. This was the first study on aquatic eukaryotes [47]. The EC50 values on the toxicity to a single species are necessary for further studies at the community and ecosystems levels [47]. Overall, qualitatively, the toxic effects were not observed to be notably different for the imidazolium and pyridinium ionic liquids, although the toxicity was greater with an increase in the length of the side chain. The K ow values for the IL were lower than those of chemicals that bioaccumulate in the tissue of organisms. Quantitative analysis by the presented QSAR/Spectral-SAR algorithm follows.  Table 7. The actions of the studied ionic liquids on the Daphnia magna species with the toxic activities Aexp= Log(EC50) [49], while the marked values were taken from [47] along with the structural parameters LogP, POL, and ETOT to account for the hydrophobicity, electronic (polarizability) and steric (total energy at optimized 3D geometry) effects, computed with the HyperChem program [103], for each cation and anion fragment, and for the anioniccationic |0+> composed state by means of equations (27)- (29), respectively [28].

Ionic Liquid Compound A exp LogP POL [Å 3 ] E TOT [kcal/mol]
The application of S-SAR-IL to Vibrio Fischeri conforms with the |1+> model presented above for the data presented in Table 3. Accordingly, the Spectral-SAR models for anions and cations are reported in Table 4, and the internal angle for all of the possible models are listed in Table 5, which confirms that model |1+>, with the results at the global level of ionic liquids being shown on the last column of Table 4.
Picturing a mechanistically mode of action for ionic liquids containing cations and anions and of their summed effects on the considered Vibrio fischeri species remains as the final purpose of the QSAR method, OECD-QSAR normative, and the Spectral-SAR algorithm, as extensively presented. In this regard, the spectral paths analysis is presented in Table 6, which is based on eqs. (33) and (34) with the cationic, anionic and ionic liquids data of Table 4. The minimum path procedure assumes the identification of the minimum paths (equal in numbers with the number of structural parameters considered, thus alpha, beta, and gamma for the present case), in an "ergodic manner"; this means that the overall path is firstly identi-fied; if two or more paths appear with equal lengths, the one to be chosen is the one that contains the "minimum first movement", which is that path that links the closest first two norms that belong to successive norms; therefore, the second overall minimum path is identified such that it does not contain the models "touched" by the alpha path, thereby excluding the common final endpoint (model III here). In this manner, the remaining beta and gamma paths are also selected. Such analysis is performed for both statistical and algebraic norms; see Table 6. In this manner, interesting results are obtained, allowing conceptual interpretation, such as the following [27]: • While in cationic case, the statistical and algebraic paths do not coincide (e.g., what the alpha Ia-IIc-III path in statistic differs than the alpha Ia-IIa-III path in algebraic views), in the anionic case, they are identically predicted (except the fact that in the algebraic case, the shortened paths are registered), which together provide a mixed behavior for the resulting ionic liquid (only the beta path Ic-IIc-III is overlapping between statistical and algebraic views); • While in the cationic and anionic subsystems, the path hierarchies are reversed as α→β→γ and γ→β→α. In the resulting ionic liquid, the mixture effect is again observed because the succession α→γ→β against the successions of starting endpoints Ia→Ib→Ic, respectively; • The cationic alpha path is started on the lipophylicity causes (Ia), which is the same as for the containing ionic liquid. A different situation arises for the alpha anionic path that begins with the steric influence (Ic); in this way, the previously noted dominance of the cationic influence when correlated with lipophilicity and the observed anionic influence related with steric effects are theoretically confirmed in this picture; • Figure 3 clearly illustrates the fact that while anionic and cationic activity tendencies are somewhat complementary, they do not cancel each other in the ionic liquid that contains them but add up to attain the overall observed toxicity, in the spectral norm -correlation factor space. Furthermore, other useful data can be extracted from Figure 3 concerning the major path of structural causes in manifested toxicological action, as revealed below; • While the algebraic paths are systematically lower that the corresponding statistical ones for cationic and anionic subsystems, in the ionic liquid case, the situation is reversed; the interpretation is that it also confirms that the chemical-biological ionic liquid dispersive (not specific) actions in environment are merely through its subsystem components than from itself as a whole; • The ecotoxicological paths in the cationic and anionic subsystems are summed up in the paths of the corresponding ionic liquids in a nontrivial manner: the anionic gamma path effect is marginal over the cationic alpha path • the cationic and anionic beta paths decay into the gamma ionic liquid path when joined such that recording a sort of reciprocal cancellation of their effects • the anionic alpha path effect is reinforced over the cationic gamma path averaging both at the beta path level of the resulting ionic liquid  Table 8. Spectral structure activity relationships (SPECTRAL-SAR) of the ionic liquids toxicity of Table 7 against the Daphnia magna species, and the associated computed spectral norms, with ║|YEXP>║=9.59481, statistic and algebraic correlation factors [26,30], throughout the possible correlation models considered from the anionic, cationic, and ionic liquid |1+> and |0+> states, respectively [28].  Table 9. The values of the cosines of the anion-cationic vectorial angles [27,30] for all considered modes of action of Table 8 indicating the |0+> states of the considered ionic liquids [28].  Table 10. Synopsis of the statistical and algebraic values of the paths connecting the SPECTRAL-SAR models of Table  8, in the norm-correlation spectral-space, for Daphnia magna species against the ionic liquids toxicity of Table 7; the primary, secondary and tertiary -the so called alpha (α), beta (β) and gamma (γ) paths, are indicated according to the "selection" and "validation" principles in norm-correlation spectral space when the statistic and algebraic variants of the correlation factors are respectively used [28].

Mode
Spectral-Structure Activity Relationship (Spectral-SAR) Assessment of Ionic Liquids' in Silico Ecotoxicity http://dx.doi.org/10.5772/51657  Table 7 on Daphnia magna species: the alpha path (α) initiates on the polarizability (Ib) anionic-cationic interaction (in the left-bottom hypersurface region), being followed by the beta path (β), which originates on the steric (Ic) anionic-cationic interaction (in the left-top hypersurface region hypersurface region), and successively by the gamma path (γ) based on the hydrophobic (Ia) anionic-cationic interaction (in the extreme left-top hypersurface region) of the norm-correlation spectral space of Table 8 with the decaying order of the thickness of the connecting arrows, respectively [28].
When turning to the study of the action of ionic liquids on the Daphnia species, the pool of molecules in Table 7 are employed for assessing the SAR-Ionic Liquid Ecotoxicological |0+> Model for IL-Daphnia chemical-biological interaction. This fact is confirmed by employing eqs. (27)- (29) for the cationic and anionic data of Table 7 and later computing the S-SAR determinants (11) associated with all models and combinations presented in Table 8, which leads to the internal angle computations and the results in Table 9 and assures the application of the model |0+> is in accordance with the prescription given by eq. (30). Within the Spectral-SAR algorithm and allied eco-toxicological principles, the results allow the specific conclusions [28]: • From the toxicological actions of Table 8, it can be observed that both anionic and cationic fragments have important contributions to the "length" and "intensity" of the ionic liquids ecotoxicity through the computed spectral norms and algebraic correlation factors, respectively, which are close to the experimental one, i.e., to 9.59481; • In all cases, the mode of action where all three Hansch factors were considered (mode III with LogP+POL+E TOT ) records the best norm and correlations that are the closest description of the ionic liquids-Daphnia magna chemical-biological interaction; • The cationic influence is observed with the dominant contribution over the anionic effects in ecotoxicity, in all considered Hansch modes of action; • The statistical correlation factors always yield smaller values than the corresponding algebraically ones, see Table 8; • There are recorded imaginary statistical correlations of the computed endpoints |Y AC-Mode > 1+ that indicate certain limitations of its use for activity modeling in ecotoxicology; for these cases, the algebraically outputs provide almost the sum of the anionic and cationic length and intensity endpoint activity. This result can be phenomenologically explained by the so-called "resonance effect" when the angles between the anionic and cationic endpoint vectors are almost zero, as clearly evidenced by the cosine values of Table 9; • Within the |0+> model, all of the lengths and intensities of the endpoints |Y AC-Mode > 0+ behave as an average of the anionic and cationic ecotoxicological effects with a smooth increase over the individual cationic effects for the modes Ib (POL), Ic (E TOT ), IIa (LogP +POL), and IIb (LogP+E TOT ); however, further selection for the binding mechanism is performed by identifying the minimum analysis of the Spectral-paths, Table 10.
The Spectral path analysis is unfolded in the same manner as previously used for the |1+> models, i.e., by "ergodic" selection of the models per paths, with the ecotoxicological results in Table 10 and correspondingly interpreted as follows [28]: • the additive parametric and endpoint models, |0+> and |1+>, provide the same hierarchies of the paths for the chemical-biological actions; • the statistical imaginary correlation values for the ionic liquids |1+> are avoided from the mechanistic principle and do not belong to any selected path in Table 10; • the dominant cationic effects can also be noted here at the least paths level because the nature of the cationic mechanism is preserved to the ionic liquids nature according with the spectral path equations: • the results of all SPECTRAL-SAR ecotoxicological principles applied to ionic liquids-Daphnia magna case of chemical-eco-biological interaction can be unitarily presented in the Figure 4, where the spectral hypersurface was generated by the 3D interpolation of all lengths (norms) for all the endpoint modes of Table 8, for all cationic, anionic, |0+> and | 1+> states of ionic liquids of Table 7. The alpha dominant paths are easily identified according to Table 10, as originating in the Ib, i.e., on POLarizability or van der Waals mo-lecular mode of action, while the beta and gamma ones starts with the steric (Ic: E TOT ) and hydrophobic (Ia: LogP) specific chemical-biological binding, respectively.
Overall, one may assess the sets of eqs. (35)- (37) and (38)- (40) as specific for ionic liquids action over biological marine species within the additive and parametric models |1+> and |0+>, respectively, that should be further confirmed or extended by future studies with other structural parameters (beyond Hansch descriptors, i.e., by topological and quantum molecular factors) and/or with other species on similar congeneric ILs and working algebraic chemical-biological interaction models.

Conclusions
Since their emergence a decade ago, ionic liquids (ILs) have had a constantly increasing influence on organic, bio-and green chemistry, due to their unique physico-chemical properties manifested by their typical salt structure: a heterocyclic nitrogen-containing organic cation (in general) and an inorganic or organic anion [43] with melting points below 100 ºC and no vapor pressure [47]. The latter property leads to the practical replacement of conventional volatile organic compounds (VOCs) from the point of view of atmospheric emissions, though they do present the serious drawback that a small amount of IL could enter the environment through groundwater [104]. This risk makes it necessary to perform further ecotoxicological studies of IL on various species to improve the "design rules" for synthesized ILs with minimal toxicity to the environmental integrated organisms. Ionic liquids display variable stability in terms of moisture and solubility in water, polar and nonpolar organic solvents [45]. Various values of ionic liquid hydrophobicity and polarity may be tailored [104] with the help of nucleoside chemistry [105] according to the main principles of green chemistry [54,66]: the new chemicals must be designed to preserve effectiveness of function while reducing toxicity and not persisting in the environment at the end of their usage but rather breaking down into inoffensive degradation products.
In this respect, the costs of all approaches for sustainable product design can be reduced using the SAR and QSAR methods [26][27][28][29][30]. While the 1-octanol-water partition coefficient could be observed only as the first approximation for compound lipophilicity, bioaccumulation and toxicity in fish, and sorption to soil and sediments, it assumes that lipophilicity is the main factor of anti-microbial activity [44,56]. Nevertheless, aiming at a deeper understanding of the specific mechanistic description of IL eco-toxicity, it is worth considering that the ionic liquid properties are more comprehensively quantified through lipophilicity, polarizability and total energy as a unitarily complex of factors in developing appropriate structure-activity relationship (SAR) studies. However, the main problem in assessing the viable QSAR studies for predicting ionic liquid toxicities concerns the anionic-cationic interaction superimposed on the anionic and cationic subsystems containing ionic liquids. There are two main complementary ways of attaining this goal [27,28,30]. One may address the search of special rules for assessing the anionic-cationic structural separately from the individual anionic and cationic ones and later generating the QSAR models -the so-called mod-el |1+>. Otherwise, when the so-called internal angle of ionic liquids sub-systems (cation and anion), eq. (30), is high enough that a sort of "resonance effects" would appear between them, that outcome should be avoided by considering structural parameter composition (superposition, as orthogonal states in quantum mechanics) such that it constitutes the way of parameter depending the predicted endpoints, that is, the so-called causal model |0+>. Together, these properties may be unitarily expressed by the operatorial equation with the Spectral-SAR operator, consistently defined by the successive rules: possible predicted end-points), thus providing the appropriate mechanistic picture of the envisaged ecotoxicity [26,30]. In this way, the presented S-SAR model appears to provide a uniform picture of the anioniccationic interaction in ionic liquids as conciliating the anionic and cationic effects observed to date. However, further studies on different species with diverse computational schemes and parameters are required, especially those that employ the quantum and algebraic features of the general operatorial equation derived in this study (41). Furthermore, by involving the enzyme-substrate modeling for receptor-effector interactions for biological activity driving chemical reactivity, as was recently reported [107], this study conceptually assessed a definitive theory of ionic liquid inter-and intra-mode of action such that it can be used for being integrated over a wide range of organisms toward designing specific eco-toxicological batteries [97], see Figure 5, for the ionic liquid's chemical-biological interactions.