A Spectral-SAR Model for the Anionic-Cationic Interaction in Ionic Liquids: Application to Vibrio fischeri Ecotoxicity

Within the recently launched the spectral-structure activity relationship (S-SAR) analysis, the vectorial anionic-cationic model of a generic ionic liquid is proposed, along with the associated algebraic correlation factor in terms of the measured and predicted activity norms. The reliability of the present scheme is tested by assessing the Hansch factors, i.e. lipophylicity, polarizability and total energy, to predict the ecotoxicity endpoints of wide types of ionic liquids with ammonium, pyridinium, phosphonium, choline and imidazolium cations on the aquatic bacteria Vibrio fischeri. The results, while confirming the cationic dominant influence when only lipophylicity is considered, demonstrate that the anionic effect dominates all other more specific interactions. It was also proved that the S-SAR vectorial model predicts considerably higher activity for the ionic liquids than for its anionic and cationic subsystems separately, in all considered cases. Moreover, through applying the least norm-correlation path principle, the complete toxicological hierarchies are presented, unfolding the ecological rules of combined cationic and anionic influences in ionic liquid toxicity.


Introduction
Since their emergence a decade ago, ionic liquids (IL) have had a constantly growing influence on organic, bio-and green chemistry, due to the unique physico-chemical properties manifested by their typical salt structure: a heterocyclic nitrogen-containing organic cation (in general) and an inorganic or organic anion [1], with melting points below 100 ºC and no vapor pressure [2]. 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 [3]. This risk makes it necessary to perform further eco-toxicological studies of IL on various species, in order to improve the "design rules" for synthesized IL with minimal toxicity to environment integrated organisms.
Ionic liquids display variable stability in terms of moisture and solubility in water, polar and nonpolar organic solvents [4]. Various values of ionic liquid hydrophobicity and polarity may be tailored [3] with the help of nucleoside chemistry [5] according to the main principles of green chemistry [6,7]: 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 breaking down into inoffensive degradation products. Most of the ionic liquids with imidazolium, phosphonium, pyridinium and ammonium that were tested were resistant to ready biodegradation [8]. Their toxicity to microorganisms can limit biodegradation, while their toxicity to humans and others organisms is obviously significant. The examination regarding the biodegradation of surfactant compounds focuses on a close resemblance between many quaternary ammonium compounds as well as surfactants based around an imidazolium core. The factors that improved the biodegradation of surfactants have successfully been applied to ionic liquids. For instance, bis(trifluoromethylsulfonyl)imide (TFMSi) and PF6 ionic liquids containing an ester in the side chain exhibit the same hydrophobic character as 1-n-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([BMIM][TFMSi]) and 1-n-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]). The enzymatic hydrolysis step, which initiates a pathway to further breakdown products, improves the biodegradation. Therefore, compound stability and toxicity are the factors biodegradability depends on. The effect of the counter-ion was not noticeable in biodegradability even though modifications of the anion led to changes in physical and chemical properties. Still, the introduction of an organic anion clearly improves the extent of ultimate biodegradation [9].
From the point of view of reactivity, IL generally do not coordinate to metal complexes, enzymes and different organic substrates [10]; however, they are usually the major component of the mixtures having pre-organized structures with the aid of many hydrogen bonds (structural directionality) in contrast to classical salts in which the compounds are mostly formed with the aid of ionic bonds (charge-ordering structures) [11,12].
On the other way, the recycling ability of IL, especially dialkyl-imidazolium based ionic liquids (the most studied until now), is based on their lack of solubility in some key organic solvents (e.g. diethyl ether) and in water -for the special case of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]). Water/soluble ionic liquids are more difficult to recycle, since their water immiscible complements, the secondary products, cannot be easily removed [13].
The detailed examination of relative energies and structural interactions (like ion position, H-bonding and anion conformational variability) in gas-phase ion-pairs has emphasized the way these quantities can be used to build up a picture of the local structure and interactions occurring in ionic liquids. For instance, while 1-n-butyl-3-methylimidazolium chloride ([BMIM][Cl]) forms a highly connected liquid with relatively strong interactions -at one extreme, and 1-n-butyl-3methylimidazolium bis(trifluoromethylsulfonyl)imide ([BMIM][TFMSi]) forms a low connectivity network of weakly linked ions -at the other extreme, [BMIM][BF4] lying between these two extremes forms a weak but more regular network. Melting points and viscosity are partly dependent on 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 [14]. Hunt et al proved that the hydrogen bond is primarily ionic with a moderate covalent character [15]. The fact that the Coulombic attraction is the dominant stabilization force was demonstrated by the analysis of the charge distribution, molecular orbitals and electron density of the dimer complex 1-n-butyl-3-methylimidazolium chloride ([BMIM][Cl]): the interactions governing the top conformers are very different from those in which the Cl anion remains in plane (were Cl anion interacts with the π manifold of orbitals); the ion-pair LUMO is the cation anti-bonding LUMO and thus electron acceptance is not favorable [15]. The effect of the chloride anion on rotation of the butyl chain is investigated and found to lower some rotational barriers while enhancing others [16]. Since the Coulombic forces are significant in an ionic liquid, charge distribution or point charges on the constituent ions are likely to be more significant than in liquids made up of neutral charge molecules [17] and deserve a more detailed study.
In this respect, the costs of all approaches for sustainable product design can be reduced using SAR and QSAR methods [18]. It has already been proved that the anti-microbial activity of quaternary ammonium chlorides is lipophilicity-dependent [19]. While the 1-octanol-water partition coefficient could be seen only as the first approximation for compound lipophylicity, bioaccumulation and toxicity in fish, as well as sorption to soil and sediments assumes that lipophylicity is the main factor of anti-microbial activity [20]. 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 lipophylicity, 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 to predict ionic liquid toxicities concerns the anionic-cationic interaction superimposed on the anionic and cationic subsystems containing ionic liquids. There are basically two complementary ways of attaining this goal. One may address the search of special rules for assessing the anionic-cationic structural separately from the individual anionic and cationic ones, and then generating the QSAR models. Yet, because the cationic and anionic effects on liquid toxicity are merely separately studied at the moment, the appropriate strategy would be to firstly derive the anionic and cationic QSARs and only then to move on to a QSAR of the ionic liquid viewed as an anionic-cationic interaction. The current paper shows, for the first time, how the latter procedure may be implemented by means of the vectorial approach of the spectral-SAR analysis. The illustration of the S-SAR-IL model presented is performed by studying the aquatic bacteria Vibrio fischeri toxicity against a list of twenty-two ionic liquids, appropriately chosen so that they should contain a wide variety of heads, side chains, and anions. This way, the present methodology and results may be extended over a wide range of organisms towards designing specific eco-toxicological ionic liquid batteries.

The Spectral-SAR Ionic Liquid (S-SAR-IL) Model
Although ionic liquid structure is currently defined through three types of substructures, the head group (the positively charged moiety) with the side chain, R1, R2, etc., which are substituents on heads group, and the anion [7], the present approach views IL as consisting mainly of its cationic (head and side chain) and anionic subsystems in mutual interaction. This way, the two sub-systems containing ionic liquid can be modeled through their vectorial activities that sum up into the predicted activity of the ionic liquid as a whole. This picture is further sustained by the possibility of employment of the Spectral-SAR method to the present purposes.

S-SAR Concepts
Without going into details [18], if one has to solve the correlation between a set of biological activities of N-compounds with the set of M-structural properties of each of them, a spectral algebraic algorithm can be applied based on the vectorial view of the descriptors and measured activity combined with the Gram-Schmidt orthogonalization recipe and on the scalar product rule giving out a real number from any two arbitrary N-dimensional vectors Let's also note that when considering the N-biological activities as well as their respective predictor variables grouped into vectors, as displayed in Table 1, the unity vector 1 ...
once it is expanded along its first column, and where With these ingredients the QSAR equation (1) becomes operational, with all coefficients effectively worked out. However, besides the effectiveness of the S-SAR methodology in reproducing the oldfashioned multi-linear QSAR analysis, one of its advances concerns the possibility of introducing norms associated with either predicted (computed) or experimental (measured) activities, thus gaining the possibility of the unique assignment of a number to a specific type of correlation, i.e. performing a sort of final quantification of the models. Nevertheless, the activity norm given in Eq. In fact this new correlation factor definition compares the vectorial lengths of the predicted activity against the measured one, thus being an indicator of the extent to which certain computed property or activity approaches the "dimension" of the observed quantity. However, it was also shown that practically the algebraic correlation factor (8) furnishes higher values than its statistical counterpart definition (7), in a systematical manner [18], thus making it the ideal tool for the present attempt to consider the ionic liquid activity taken from its cationic and anionic sub-system ones.
Moreover, with the help of both spectral norms and correlation factors we can introduce the so called the least path principle in terms of paths between the tested models (or endpoints) [18]: thus providing a practical tool for deciding the dominant hierarchies along the possible ones with the important consequence of picturing the mechanistic and almost temporal evolution of structural causes that trigger the observed effects. This methodology was successfully applied in ecotoxicology [18] and for designing the behavior of the species interactions within a test battery [21], promising to furnish the frame also in ionic liquids analysis.

S-SAR for Ionic Liquids
Basically, since we may consider the ionic liquid as composed by the anionic and cationic subsystems, we can employ the spectral vectorial space of anionic and cationic activities to form the ionic liquid one throughout the vectorial resultant as displayed also in Figure 1. The vectorial summation in (11) may also be seen as the interference of the anionic and cationic activities in the space of ionic liquid so that their mutual interaction is included. Actually, searching for the predicted norm of the vectorial ionic liquid activity one can notice the interference effects between the anionic and cationic systems through the correlation angle between the associated vectors, cosθ AC . Nevertheless, its practical definition can be achieved, since the anionic-cationic norm (13) is rewritten employing the vectorial relation (11) to the scalar product between the anionic and cationic vectors Now, from the two equivalent expressions (13) and (14) Basically, the degree of anionic-cationic activity interaction in ionic liquid is fixed by the value of the angle (15) for each envisaged QSAR model or endpoint.
Finally, aiming to obtain the ionic liquid working QSAR equation from the cationic and anionic counterparts one simply needs to expand the norm and the scalar product in (14) as the square sum of components, according to the norm and scalar product definitions (12) and (2) With these all the required concepts and analytical tools are given for that any ionic liquid anioniccationic interaction and resulted activity analysis on certain species or organism be performed, while concrete example will in next be exposed.

The Working System
Usually, the classical assessment of new industrial chemicals is inflexible, being directed by regulations and standardized procedures, while biological test systems are very expensive and the costs increase with the number of new chemicals released. So, a flexible test strategy is needed in correlating the chemical compounds with the systems they act upon. Although qualitative (eco) toxicological algorithms for selection of test systems have been proposed in terms of identification of individual effects of different head-groups (with identical side chains and anions) and different anions (with identical cations) [7], a quantitative theoretical prediction algorithm for tested compounds in biological systems is still need be proposed due the lack of readily accessible ecotoxicological data [8].
The principles of green chemistry say that one has to consider the whole process (life cycle analysis) rather than individual components of reaction (single issue sustainability), based on the fact that the acute toxicity measurements do not provide a complete characterization of the full impact of a substance release into environment, but are only part of the environmental impact assessment [6,8].
In this context, ionic liquids with cations like pyridinium, imidazolium and pyrrolidinium have already been nominated in the United States National Toxicology Program (NTP) for toxicological testing based upon their potential of new solvents but also due they ability to enter in aquatic system [22]; if an accidentally discharge of ionic liquids into water occur, many of them being water soluble, they may be an environmental risk to aquatic plants and animals [23].
For this reason the current application will develop the S-SAR complete algorithm up to the mechanistic prediction of the correlated structural causes and ecotoxicological effects for the ionic liquids of Figure 2, containing ammonium, pyridinium, phosphonium, choline, and imidazolium cations, on aquatic bacteria Vibrio fischeri [24]. Whereas the Vibrio fisheri species was previously found to be one of the most resistant species to ordinary phenol compounds toxicity [21], the present ionic liquids show quite a wide structural variety to furnish useful information of their environmental action. The cationic and anionic structural properties, the lipophylicity, the electronic polarizability and the total energy, where computed with the HyperChem computational environment [25] and displayed in Table 2 together with the reported measured activity of their containing ionic liquids, respectively.
The data in Table 2 are suitable arranged so that the S-SAR analysis is performed successively at the cationic and anionic level and then at the ionic liquid levels based on equations (4)- (10) and (15) when accounting for different combinations between the lipophylicity (hydrophobic character), polarizability (electronic character) and total energy (steric character) factors, respectively. The spectral hierarchy of the predicted activities with respect to different concerned endpoints would lead to the mechanistic scheme according with the cationic and anionic sides as well as the overall ionic liquids influences the environmental (Vibrio fischeri) toxicity.

Figure 2:
The ionic liquids acting on Vibrio fischeri species studied in this work.  Figure 2 of those toxic activities A= Log(EC 50 ) on Vibrio fischeri were considered [24], with the marked values taken from [4], along structural parameters LogP, POL (Å 3 ), and E TOT (kcal/mol) as accounting for the hydrophobicity, electronic (polarizability) and steric (total energy at optimized 3D geometry) effects, computed with the help of HyperChem program [25], for each cation and anion containing ionic liquid, respectively.

Results and Discussion
As earlier mentioned, the first step in our analysis consists in deriving the cationic and anionic QSARs that link the structural lipophilic-electronic-steric parameters of the ionic liquids of Figure 2 and Table 2 with the observed activities of the whole containing ionic liquids upon the Vibrio fischeri species. The results are presented in Tables 3 and 4 for the cationic and anionic subsystems for all main combinations, i.e. generating the uni-modes Ia-to-Ic when only one structural parameter is correlated, the two-modes IIa-to-IIc when two combined structural factors are taken into account and for the three-mode correlation III with all structural factors involved, respectively. For each such mode of action, the associated endpoint norm, the statistic and algebraic correlation factors were reported, computed with the equations (6), (7), and (8), respectively. Table 3: Spectral structure activity relationships (S-SAR) of the ionic liquids of Figure 2 against Vibrio fischeri toxicity, and the associated computed spectral norms with ║|Y EXP >║=4.41537, statistic and algebraic correlation factors, computed upon the relations (6), (7), and (8), throughout the possible correlation models considered from the cationic data in Table 2, respectively. As a general observation, in all cases there was recorded a systematic increase of the correlation factor when computed in spectral space, i.e. using the algebraic definition, as compared with the standard statistical values. We would like to take this opportunity to advocate the use of the algebraic definition instead of the statistical one since the first one has the physical meaning of the "length of action" respecting the old-fashioned dispersion analysis. Nevertheless, dispersion being a consequence of the appropriateness of the fit, the vectorial norm or the "length of action" accounts merely for the degree with which a certain model approaches the observed, or measured or manifested (chemicalbiological) interaction. In this respect, it is also worth noting that both cationic and anionic predicted activities poorly resemble the experimentally expected activities with a smooth increase on the anionic influence, for all computed models except two cases based on the lipophylicity correlation (Ia and IIa).

Mode Vectors Cationic S-SAR
Apparently, this behaviour disagrees with the previously reported studies in which the anionic effect was only marginal in cytotoxicity [1,4,26], but in some special cases of certain ionic liquids tests [7]. Table 4: Spectral structure activity relationships (S-SAR) of the ionic liquids of Figure 2 against Vibrio fischeri toxicity, and the associated computed spectral norms with ║|Y EXP >║=4.41537, statistic and algebraic correlation factors, computed upon the relations (6), (7), and (8), through all possible correlation models considered from the anionic data in Table 2, respectively. This situation was based on the observation that, for instance, since imidazolium ring in cations is a delocalized aromatic system with high electron acceptor potential, the nitrogen atoms are not capable to form any hydrogen bonds and the result is a very rigid and sterically inflexible system, while the elongation of R2 residue in side chain leads to a continuous increase of flexibility implying more conformational freedom [19]. The reason for this actions relies on the fragmental hydrophobicity of each carbon connected to a quaternary amine which, combining the geometric bond factor (that applies to the neutral solute) with a negative electronic bond factor, decreases its magnitude with the square of the distance from the central nitrogen atom [26,27]. Table 5: Spectral structure activity relationships (S-SAR) of the ionic liquids of Figure  2 against Vibrio fischeri toxicity, and the associated computed spectral norms with ║|Y EXP >║=4.41537, statistic and algebraic correlation factors, computed upon the relations (6), (7), and (8), through all possible correlation models considered for the cationicanionic interaction (11), respectively. Moreover, the systematic variation of R1, R2, etc., at identical head groups and anions, in all published data, from the molecular to individual organism level, leads to the conclusion that the shorter the chain lengths of side chains the lower the cytotoxicity (higher EC50 values) [7,9]. The electronic portion of the bond factor extends along a chain of no more than 5 alkane carbons causes a decrease in overall hydrophobicity so the chains longer than 5-6 atom carbons will have a greater permeability through the cell membrane (see for example [DMIM][BF4]) [26]. This, probably because the chemical transformation of the side chains of ionic liquids may reduce toxicity as far as the metabolites are less toxic compared to their parent chemicals [7]. Table 6: The variation of the cosines of the anion-cationic correlation angle in vectorial space, according with the Eq. (15), for all considered modes of action of the ionic liquids of Figure 2 with the cationic and anionic subsystems S-SAR activities given in Tables 3 and 4 Next, aiming to combine the two somewhat separate effects of cations and anions in the ionic liquid activities, the spectral-SAR results are given in the Table 5, showing the working ionic liquid QSAR equation as well the actual vectorial norm, statistic and algebraic correlation factors for each mode of action envisaged so far.

Mode Ionic Liquid S-SAR
The data in Table 5 clearly demonstrate that the ionic liquid S-SAR models always predict higher norms in endpoint activities thus providing the considerable increase in the algebraic correlation factors with respect to the cationic and anionic subsystem effects. It is very interesting to see that the statistical values not only furnish lower values than the algebraic outputs for the ionic liquids, but often lie even below the corresponding statistical values of the cationic and anionic subsystems. This situation gives us a chance to establish the limits of using the dispersion based correlation factor definition since it does not properly reproduced the addition effect of the two interacting subsystems as the anionic-cationic interaction in the ionic liquid structures. On the other hand, the mixture effect of the cationic and anionic vectorial actions in (eco)toxicological studies is well established as far as the single substances of a mixture acts in similar and close quantified manner in all considered modes.
The fact that this is the present case can be firstly visualized by the close inspection of Tables 3 and  4 where the spectral norms of the predicted activities feature close values in relatively narrow domain of actions. Moreover, the reinforcement of this idea comes from the data of Table 6 in which the angle of interactions between the cationic and anionic subsystems are computed, as given by the formula (15), for all considered endpoints with almost constant values around the 0.600 cosines of the angle, while the only higher fluctuation deviation appears in the Ia and IIa cases -the same previously evidenced for the cationic dominancy over the anion effects.
Overall, once the electronic and steric mechanisms have also been considered, the anion could play a central role as technicophore because it exhibits a high potential for change technological properties (solubility, viscosity) or due its peculiarity to partially decompose itself in the ion pair interaction [26,28], leading to its dominant effect in the chemical-biological engaged activity.
With these consideration we can safely assume that the Spectral-SAR in general and the algebraic correlation factors in particular are especially suited for modelling the (eco)toxicological activities of ionic liquids from its anionic and cationic component effects. Table 7: Synopsis of the statistic and algebraic values of paths connecting the S-SAR models of Table 5, in the norm-correlation spectral-space of Figure 3, for the ionic liquids of Figure 2 against Vibrio fischeri toxicity. The primary, secondary and tertiarythe so called alpha (α), beta (β) and gamma (γ) paths, are indicated according to the least path principle in spectral norm-correlation space with the statistic and algebraic variants of the correlation factors used, respectively. The final part of discussion is devoted for picturing a mechanistically mode of action for cationic, anionic and of their summed effects in containing ionic liquids on the considered Vibrio fischeri species. That is, the minimum path procedure among all possible ways connecting endpoints from each category of models (i.e. with one, two or three factors dependency) is to be considered. The path lengths were computed employing the equation (10) to all cationic, anionic and ionic liquids data of Tables 3, 4, and 5, and the results are displayed in the Table 7, respectively.

Value
However, in order to identify the shortest paths in each category of endpoint connections, the following rules are applied: the first choice is the overall minimum path in a certain column of Table 7 (i.e. a system with a specific way of correlation factor); if the overall minimum belongs to many equivalent paths (as is the case of cationic with algebraic factor column in Table 7, for instance) the minimum path will be considered the one that links the starting endpoint with the closest endpoint in the sense of norms (as is the norm of IIa mode the closest to the norm of Ia mode in cationic-algebraic column of Table 7, for example); the overall minimum path will set the dominant hierarchy of the mechanistically mode of action towards the experimentally observed activity and will be called the alpha path (α); once the alpha path has been set the next minimum path will be looked for such that the starting endpoint should be different from the one involved in the alpha path (that is, if in the alpha path the starting end point was Ia, the next path to be identified will begin either from the Ib or the Ic mode); the next minimum paths are chosen on the same rules as before and will be called as beta and gamma paths, β and γ, respectively. At the end of this procedure each mode of action is "touched" only once, except for the final endpoint, here III (as a computational substitute of EXP), so that all the methodology being regarded as searching minimum path throughout variation of paths with the fixed final endpoint. Now, the alpha, beta and gamma path can be easily identified in Table 7 and there are accordingly marked.

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 on Table 7 data, for the cationic, anionic and resulted ionic liquid norms of Tables 3-5, for the statistic and algebraic versions of correlation factors, from up to down and left to right, respectively. The primary-alpha, secondary-beta and tertiary-gamma path hierarchies of Table 7 are indicated by decreasing the thickness of the connecting lines.
With the aim of giving the complete results of the analysis in a single shot, Figure 3 is a "spectral" representation of the data, with the norms for each cationic, anionic, and resulted cationic liquid being drawn, linked by the major hierarchical paths against the considered mode of actions toward the observed "length" (norm) of the toxicity activity of the Vibrio fischeri species, for both statistical and algebraic correlation pictures. Figure 7 clearly underlines 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. Besides this, some other useful information can be extracted from Figure 7 concerning the major path of structural causes in manifested toxicological action as revealed bellow.
While in cationic case the statistical and algebraic path does not coincide at all (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 (excepting the fact that in the algebraic case the shortened paths are registered), together providing a mixed behavior for the resulted ionic liquid (only the beta path Ic-IIc-III is overlapping between statistic and algebraic views). While in cationic and anionic subsystems the path hierarchies are reversed, as α→β→γ and γ→β→α, in the resulting ionic liquid again the mixture effect is recorded since the succession α→γ→β, against the successions of starting endpoints Ia→Ib→Ic, respectively. It is worth pointing out here that, indeed, the cationic alpha path is started on the lipophylicity causes (Ia), the same as for the containing ionic liquid, a different situation arising for the alpha anionic path that is beginning with the steric influence (Ic). This way the previously noted dominance of cationic influence when correlated with lipophylicity as well the recorded anionic influence related with steric effects are in this picture theoretically confirmed.
Other interesting features about ecotoxicological paths rely on the fact that, while for cationic and anionic subsystems the algebraic paths are systematically lower that the corresponding statistical ones, in the ionic liquid case the situation is reversed. The interpretation of this fascinating result 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. This result further motivates the design of ionic liquids by tailoring the properties of its containing anionic and cationic components for prevent their hazardous toxicity.
Finally, let us also note that the ecotoxicological paths in cationic and anionic subsystems are summed up in the paths of the corresponding ionic liquids in a not trivial ways: From these S-SAR equations some conceptual, however computationally based, ecotoxicological path rules for ionic liquids can be concluded, namely: the anionic gamma path effect is marginal over the cationic alpha path -equation (19); the cationic and anionic beta paths decay into the gamma ionic liquid path when met together so that recording a sort of reciprocal cancellation of their effects -equation (20); the anionic alpha path effect is reinforcing over the cationic gamma path averaging both at the beta path level of the resulted ionic liquid -equation (21).
This way, the S-SAR model presented seems to provide a unitary picture of the anionic-cationic interaction in ionic liquids as conciliating the anionic and cationic effects observed so far. However, further studies on different species with diverse computation schemes and parameters are required in order to conceptually asses a definitely [a definitive] theory of ionic liquid inter-and intra-mode of action.

Conclusions
Despite the promise to change "the face" of organic chemistry, the ionic liquids are becoming a central issue in green chemistry as well, due to their unique physico-chemical properties. For instance, their high viscosities compared with conventional molecular solvents have impact on lowering the reaction rates while through their higher degree of immersion into ground soil and sediments the ecotoxicity of various species and organisms are recorded. The latter effects are based on their common supposed mechanism through membrane disruption: having structural similarity with detergents, pesticides and antibiotics [4,20,29,30] they induce polar narcosis due to their interfacial properties and may cause membrane-bond protein disruption [2,24]. Ionic liquids modes of action are based on the fact that the disrupt membranes and hydrophobic molecules have a greater ability to accumulate at this interface [9], while other mechanisms may arise from acetylcholinesterase inhibition [31], from common cellular structure or processes [4], from structural DNA damage [23,32]. Eventually, some bacteria could potentially break down imidazolium into different metabolites which follow further inhibition on themselves [1]. In any case, low lipophylicity imidazolium based ionic liquids values indicate low permeability of these ionic liquids [2]. In addition, to the certain role that the lipophylicity has in the assessment of ionic liquid toxicity on various species, before constructing kit-or battery tests involving chemical-biological level of organization, from enzyme to cellular to short and long lived organism [2,4,[7][8][9]23,24,26,[31][32][33][34][35][36][37][38][39][40][41][42][43][44], worth focusing on the quantitative algorithm with the help of which the experimental data will be interpret. This is because, criticized for its lacking in real ecological meaning, the dose-response approach for estimating the lethal effects of toxicants on organisms regulatory norms have been built around LC50 values that can be compared between different toxicants and organisms [2,45], while the statistical correlation factors of the proposed QSARs so far are based on dispersion (i.e. hazardous) action and less accounting on the specific ways of ionic liquids action. In this regard, it is important to know if the tested chemicals could reach the target site into organism in a specifically manner first, so that an extrapolation from in vitro to in vivo is necessarily [31] while other factors like biodegradation and bioaccumulation are necessarily be considered before conclusions are drawn [24,46,47].
Having to consider both specific and hazardous paths of action concerning the environmental risk posed by the newly released ionic liquids the recently introduced Spectral-SAR scheme of quantification of chemical-biological interaction [18] is adapted here to the particular ionic liquid cationic-anionic structural interaction. The present method is based on the previous evaluation of QSARs with specific norms and algebraic correlation factors for cationic and anionic influences followed by their vectorial interference in the spectral norm space. With a short generation lifetime, bacteria is an ideal starting point for ionic liquid structure-activity relationship investigation and serves as a basis for further toxicity tests to higher organisms and more complex systems. This way, the S-SAR-IL model adapted here was applied to Vibrio fischeri species, this species being reported as the most resistant to other conventional environmental releasing chemical compounds [21]. The results are promising, confirming at some level of analysis all previous certified facts, the cationic generally dominance [48], the special anionic effects [49], the ionic pair effect [50], etc., adding in an elegant manner the least (or minimum) path rule in the spectral norm-correlation factor space as an effective tool in which the ecotoxicological rules of the ionic liquids can be derived. Worth, however, noted that when choosing the Hansch factors, although there was already pointed out that the use of free enthalpy instead of total energy provides better statistical correlation factors [51], in the present study the total energy was preferred due to its close connection with the sterical effects through the computational geometrical optimization procedure involved. This way, the minimal comprehensive set of QSAR descriptors, i.e. hydrophobic, electronic, and steric, when predict the ecotoxicity endpoints was assured [21]. In this respect, further studies may address the effect of the length of alkyl chain of the imidazolium-based ionic liquid on the activity.
The present approach leaves room for further investigation at both the conceptual and computational levels of expertise when improving the specificity of analysis in terms of considered structural factors and of their SAR combinations or when generalizing it to the eco-design of multiple species battery, respectively.