Study of the Relationships between the Structure, Lipophilicity and Biological Activity of Some Thiazolyl-carbonyl-thiosemicarbazides and Thiazolyl-azoles

Lipophilicity, as one of the most important physicochemical parameters of bioactive molecules, was investigated for twenty-two thiazolyl-carbonyl-thiosemicarbazides and thiazolyl-azoles. The determination was carried out by reversed-phase thin-layer chromatography, using a binary isopropanol-water mobile phase. Chromatographically obtained lipophilicity parameters were correlated with calculated log P and log D and with some biological parameters, determined in order to evaluate the anti-inflammatory and antioxidant potential of the investigated compounds, by using principal component analysis (PCA). The PCA grouped the compounds based on the nature of their substituents (X, R and Y), indicating that their nature, electronic effects and molar volumes influence the lipophilicity parameters and their anti-inflammatory and antioxidant effects. Also, the results of the PCA analysis applied on all the experimental and computed parameters show that the best anti-inflammatory and antioxidant compounds were correlated with medium values of the lipophilicity parameters. On the other hand, the knowledge of the grouping patterns of the tested variables allows the reduction of the number of parameters, determined in order to establish the biological activity.


Introduction
Lipophilicity is an important physicochemical property of bioactive compounds affecting their biological activity with a determinant role in the transport of compounds through biological membranes and in the formation of the ligand-receptor complex.The lipophilic character of an active molecule, defined as the ability of a compound to penetrate through hydrophobic barriers in order to get from the delivery point to the site of action [1,2], is usually quantitatively characterized as the logarithmic forms of the n-octanol-water partition coefficient P ow (LogP ow ) and n-octanol-water distribution coefficient D ow (LogD ow )-experimentally determined [3] or calculated by using a series of mathematical models [4].Due to experimental limitations of the classical "shake-flask" method [5], the most widely used techniques for the measurement of the lipophilic properties of different chemical molecules are nowadays the chromatographic techniques in reversed-phase systems (RP-TLC and RP-HPLC) [6,7].Some of the parameters resulting from a RP-TLC may be associated to the lipophilic character of the analytes [8,9].The retention parameters (R M ) are calculated by means of well-known Batte-Smith and Westall Equation ( 1): R M " log p1{R F ´1q Generally, the R M values determined by RP-TLC are linearly dependent on the concentration of the organic modifier (C) in the mobile phase Equation ( 2): where b and R M0 are, respectively, the slope and the intercept of Equation (2).Extrapolation of the R M value to pure water, leads to determination of R M0 , which can be considered as an estimation of the partitioning of compounds between nonpolar stationary phase and the aqueous system, and hence a way to estimate the lipophilicity of the compounds (Soczewinski-Wachtmeister model) [9,10].Additionally, it has been stated that not only the R M0 , but also the slope, b, as a characteristic of the specific hydrophobic surface area of the compounds, can be used as an estimation of lipophilicity [6,11,12].The advantages of TLC methods consist in the small quantity of compounds needed for the determination, and the rapidity and simplicity of the method.Also, they are inexpensive, rapid and easy to perform [13].
Principal Component Analysis (PCA) is a chemometric tool designed to transform a set of original variables into new uncorrelated variables (axes), which are called principal components.The new variables are linear combinations of the original variables and the new axes lie along the directions of maximum variance.PCA provides an objective way of finding indices of this type, so that the variation in the data can be accounted for, as concisely as possible.By applying PCA to the matrix formed by the retardation factors (R F ) of all compounds in all mobile phases, it is possible to obtain a new lipophilicity scale where the linear combinations of retention indices (the scores) corresponding to the first component (PC1) appears to be a new parameter capable to quantitatively assess the lipophilic character of the compounds [6,[14][15][16].
Based on our experience in the field of heterocyclic compounds [17][18][19], we synthesized and evaluated in vivo, in term of anti-inflammatory and antioxidant activity, a new series of thiazolyl-carbonyl-thiosemicarbazides and hybrid thiazolyl-1,3,4-oxadiazoles, thiazolyl-1,3,4triazoles and thiazolyl-1,3,4-triazoles [20,21].The results demonstrated that the new thiazole compounds exhibit anti-inflammatory effects, lowering the acute-phase response of bone marrow and the oxidative stress.Based on the virtual screening and on the obtained results, the activity may be due to their capacity to reduce NO synthesis by blocking the binding of L-arginine at the active site of inducible nitric oxide synthase (iNOS) [21].
The goal of the present study was to determine the relationships between the computed lipophilicity coefficients and RP-TLC retention parameters of new thiazolyl-carbonyl-thiosemicar bazides and thiazolyl-azoles with the anti-inflammatory/antioxidant potential and their structural parameters obtained by molecular modeling calculations, applying the PCA method.
Molecules 2015, 20, page-page 3 according to a previously described procedure [20], were investigated.The structures of the compounds used in this study are presented in Figure 1.All components of the mobile phases used were of analytical grade purity.

Lipophilicity Evaluation of the Studied Compounds
In all instances, it was found that the retention of the investigated compounds increased with increasing concentration of the organic modifier in the mobile phase.The change of the retardation factor with increasing the organic solvent content in the binary mobile phase, correlates well with the polarity, the substituents (phenyl, allyl, methyl, hydrogen and bromine).The RM values obtained from Equation (1) decrease linearly as the concentration of the organic modifier in the mobile phase increases, and therefore they might be used to assess lipophilicity.The profiles of the RM values present regular changes as the organic modifier content increases, indicating that the same mechanism of lipophilic interaction is dominant (Figure 2).Nevertheless, minor changes of these systematic regularities can be observed with compounds Th-11-17, most probably due to their more rigid molecular structures compared to compounds Th-1-8.
The results of regression analyses using Equation ( 2) (intercept values (RM0), slopes (b) and correlation coefficient (r) are presented in Table 1.These values afford a quantitative estimation for the distribution of the investigated compounds between a non-polar phase (a chemically bound reversed stationary phase) that represents the biological membranes, and a polar phase that represents the extracellular aqueous environment.As expected, it can be observed that there is a linear dependence between the values of RM0 (lipophilic parameter) and the values of b (specific hydrophobic surface) for the majority of these compounds.This linear dependence shows that thiazolyl derivatives might form a homologous series of compounds, as has been suggested by some authors [22].Once more, it can be observed that compound Th-1-8 has a mean correlation coefficient of 0.979 ± 0.005, while compounds Th-9-17 and Th-18-22 have 0.927 ± 0.034 and 0.925 ± 0.018, respectively.In addition,

Lipophilicity Evaluation of the Studied Compounds
In all instances, it was found that the retention of the investigated compounds increased with increasing concentration of the organic modifier in the mobile phase.The change of the retardation factor with increasing the organic solvent content in the binary mobile phase, correlates well with the polarity, the substituents (phenyl, allyl, methyl, hydrogen and bromine).The R M values obtained from Equation (1) decrease linearly as the concentration of the organic modifier in the mobile phase increases, and therefore they might be used to assess lipophilicity.The profiles of the R M values present regular changes as the organic modifier content increases, indicating that the same mechanism of lipophilic interaction is dominant (Figure 2).Nevertheless, minor changes of these systematic regularities can be observed with compounds Th-11-17, most probably due to their more rigid molecular structures compared to compounds Th-1-8.
The results of regression analyses using Equation (2) (intercept values (R M0 ), slopes (b) and correlation coefficient (r) are presented in Table 1.These values afford a quantitative estimation for the distribution of the investigated compounds between a non-polar phase (a chemically bound reversed stationary phase) that represents the biological membranes, and a polar phase that represents the extracellular aqueous environment.As expected, it can be observed that there is a linear dependence between the values of R M0 (lipophilic parameter) and the values of b (specific hydrophobic surface) for the majority of these compounds.This linear dependence shows that thiazolyl derivatives might form a homologous series of compounds, as has been suggested by some authors [22].Once more, it can be observed that compound Th-1-8 has a mean correlation coefficient of 0.979 ˘0.005, while compounds Th-9-17 and Th-18-22 have 0.927 ˘0.034 and 0.925 ˘0.018, respectively.In addition, similar significant differences (p < 0.0003, 12.56 < f < 25.34, ANOVA test) could be observed between these three groups, for R M0 and b values (Table 1).Based on these observations and together with structural aspects (Table 1), the studied compounds have been divided in three groups: compounds Th-1-8 form group A, while compounds Th-9-17 and Th-18-22 form group B and group C, respectively.This classification will be used in further observations and discussions (vide infra).
Molecules 2015, 20, page-page 4 similar significant differences (p < 0.0003, 12.56 < f < 25.34, ANOVA test) could be observed between these three groups, for RM0 and b values (Table 1).Based on these observations and together with structural aspects (Table 1), the studied compounds have been divided in three groups: compounds Th-1-8 form group A, while compounds Th-9-17 and Th-18-22 form group B and group C, respectively.This classification will be used in further observations and discussions (vide infra).Besides the well-known R M0 values which are the experimentally-determined indices for the lipophilicity of the studied compounds, PC1 values, which are obtained by applying PCA to the matrix of the R F values for all mobile phases for all compounds, can be successfully used for quantifying the lipophilicity and may bring complementary information to the R M0 [14,15].The variation of these two parameters for the studied compounds, along with the computed LogP and LogD values, is presented in Figure 3.It can be observed that the R M0 and PC1 profiles are well correlated with the profiles of the computed LogP and LogD values, all describing similar patterns.Regarding the experimental lipophilicity parameters, R M0 and PC1, the compounds present expected variations of the lipophilicity, as the functional groups are modified (both R and X groups).
Molecules 2015, 20, page-page 5 Besides the well-known RM0 values which are the experimentally-determined indices for the lipophilicity of the studied compounds, PC1 values, which are obtained by applying PCA to the matrix of the RF values for all mobile phases for all compounds, can be successfully used for quantifying the lipophilicity and may bring complementary information to the RM0 [14,15].The variation of these two parameters for the studied compounds, along with the computed LogP and LogD values, is presented in Figure 3.It can be observed that the RM0 and PC1 profiles are well correlated with the profiles of the computed LogP and LogD values, all describing similar patterns.Regarding the experimental lipophilicity parameters, RM0 and PC1, the compounds present expected variations of the lipophilicity, as the functional groups are modified (both R and X groups).Analyzing the values obtained for PC1, it was noticed that the cyclisation of thiosemicarbazides Th-1-8 to the corresponding thiazolyl-triazoles Th-9-17 increased the lipophilicity, probably as result of a decrease of the polar character of the former, due to the -CO-NH-NH-CS-NH-fragment.Concerning the nature of the Y substituent, it can be observed that the presence of an electron withdrawing substituent (phenyl) increases the lipophilicity, while electron donating groups decrease it (Th-11 > Th-14 > Th-17) (Figure 3).Regarding the X substituent and the whole molecules, as expected in these homologous series, the lipophilicity increases with the molecular volume, the bromo substituted derivatives being the most lipophilic compounds.Also, the N-substituted thiazolyl-triazoles (Th-9-17) are more lipophilic then the corresponding thiadiazoles (Th-21-22) and oxadiazoles (Th-18-20), the molecular volume decreasing in the order: group B > group C (S > O) (Figure 3).
According to their structure (Figure 1), the studied compounds can be involved in acid-base equilibria, and are hence ionizable, so that the analysis of the distribution coefficient (LogD) at several physiological relevant pH levels is adequate.Figure 3 also presents the profiles of the LogD values at four pH levels (stomach (1.5), duodenum (5), jejuno-ileum (6.5) and blood serum (7.4) pH levels).The compounds from group A (Th-1-8) have slightly basic behavior; the LogD values increase as the pH increases, with a more pronounced effect, thus a much wider interval of distribution, for the first three compounds-which can be explained by the electron withdrawing inductive effect of the phenyl group.Mainly due to the presence of the thiol group, the compounds from groups B and C (Th-9-22) Analyzing the values obtained for PC1, it was noticed that the cyclisation of thiosemicarbazides Th-1-8 to the corresponding thiazolyl-triazoles Th-9-17 increased the lipophilicity, probably as result of a decrease of the polar character of the former, due to the -CO-NH-NH-CS-NH-fragment.Concerning the nature of the Y substituent, it can be observed that the presence of an electron withdrawing substituent (phenyl) increases the lipophilicity, while electron donating groups decrease it (Th-11 > Th-14 > Th-17) (Figure 3).Regarding the X substituent and the whole molecules, as expected in these homologous series, the lipophilicity increases with the molecular volume, the bromo substituted derivatives being the most lipophilic compounds.Also, the N-substituted thiazolyl-triazoles (Th-9-17) are more lipophilic then the corresponding thiadiazoles (Th-21-22) and oxadiazoles (Th-18-20), the molecular volume decreasing in the order: group B > group C (S > O) (Figure 3).
According to their structure (Figure 1), the studied compounds can be involved in acid-base equilibria, and are hence ionizable, so that the analysis of the distribution coefficient (LogD) at several physiological relevant pH levels is adequate.Figure 3 also presents the profiles of the LogD values at four pH levels (stomach (1.5), duodenum (5), jejuno-ileum (6.5) and blood serum (7.4) pH levels).The compounds from group A (Th-1-8) have slightly basic behavior; the LogD values increase as the pH increases, with a more pronounced effect, thus a much wider interval of distribution, for the first three compounds-which can be explained by the electron withdrawing inductive effect of the phenyl group.Mainly due to the presence of the thiol group, the compounds from groups B and C (Th-9-22) present rather an acidic behavior, which is more evident for compounds Th-9-11 and in group C (Th- [18][19][20][21][22].This may suggest that compounds Th-1-3 are better absorbed in intestine than in stomach, while compounds Th-9-11 and Th-18-22 are better absorbed in stomach and less in intestine, due to their higher uncharged population of molecules in the first mentioned compartments.The other compounds (Th-4-9, Th-12-17) are to a lesser extent affected by pH variation at the studied interval (i.e., 1.5-7.4).
Applying PCA on both experimentally and computed lipophilicity indices on all studied compounds, 2D plot of the loadings (circle correlation) showed a regular increase of the correlation between LogD values and experimentally determined lipophilicity indices (PC1, R M0 ) as the pH increases (Figure 4 and Figure S1).However, since the compounds behaved differently according to their group, a PCA application on two separated groups: A (basic compounds) and B and C (acidic compounds) revealed interesting results (Figure 5 and Figure S2).As expected, in both groups, PC1 and R M0 positively correlate with each other (r = 0.804) and these indices present a strong negative correlation with b (r = ´0.965with R M0 and r = ´0.621with PC1) since it is the characteristic of the specific hydrophobic surface area [6,11,12].LogP and LogD values afford a relatively moderate correlation with the experimentally determined indices, most probably due to the limitation of estimation of computed indices for complex compounds, since none of the available methods can take into consideration all the effects of molecular conformation.However, as can be observed from Figure 5, LogD at acidic pH in the case of group A and LogD at basic pH in the case of groups B and C are more correlated with R M0 and PC1 values.
The scatterplot of the scores of the first three principal components after applying PCA on all the previously discussed lipophilicity indices reveals an obvious clustering of the compounds, according to their properties and structural aspects, in very good agreement with the grouping performed in Table 1 (Figure 6).Group A and group B overlap to some extent but group C forms a distinct cluster despite the high scattering within it, most probably due to numerous changes caused by the presence of the heteroatom (Y substituent, Table 1).Compound Th-17 appeared as an outlier due to its distinct higher R M0 value.
Molecules 2015, 20, page-page 6 present rather an acidic behavior, which is more evident for compounds Th-9-11 and in group C (Th- [18][19][20][21][22].This may suggest that compounds Th-1-3 are better absorbed in intestine than in stomach, while compounds Th-9-11 and Th-18-22 are better absorbed in stomach and less in intestine, due to their higher uncharged population of molecules in the first mentioned compartments.The other compounds (Th-4-9, Th-12-17) are to a lesser extent affected by pH variation at the studied interval (i.e., 1.5-7.4).
Applying PCA on both experimentally and computed lipophilicity indices on all studied compounds, 2D plot of the loadings (circle correlation) showed a regular increase of the correlation between LogD values and experimentally determined lipophilicity indices (PC1, RM0) as the pH increases (Figure 4 and Figure S1).However, since the compounds behaved differently according to their group, a PCA application on two separated groups: A (basic compounds) and B and C (acidic compounds) revealed interesting results (Figure 5 and S2).As expected, in both groups, PC1 and RM0 positively correlate with each other (r = 0.804) and these indices present a strong negative correlation with b (r = −0.965with RM0 and r = −0.621with PC1) since it is the characteristic of the specific hydrophobic surface area [6,11,12].LogP and LogD values afford a relatively moderate correlation with the experimentally determined indices, most probably due to the limitation of estimation of computed indices for complex compounds, since none of the available methods can take into consideration all the effects of molecular conformation.However, as can be observed from Figure 5, LogD at acidic pH in the case of group A and LogD at basic pH in the case of groups B and C are more correlated with RM0 and PC1 values.
The scatterplot of the scores of the first three principal components after applying PCA on all the previously discussed lipophilicity indices reveals an obvious clustering of the compounds, according to their properties and structural aspects, in very good agreement with the grouping performed in Table 1 (Figure 6).Group A and group B overlap to some extent but group C forms a distinct cluster despite the high scattering within it, most probably due to numerous changes caused by the presence of the heteroatom (Y substituent, Table 1).Compound Th-17 appeared as an outlier due to its distinct higher RM0 value.

Antioxidant and Anti-Inflammatory Effects of the Studied Compounds
The new thiazole compounds have anti-inflammatory effects by lowering the acute-phase response of bone marrow (L and PMN), phagocytic capacity (PI and PA) and oxidative stress [21].The best anti-inflammatory and antioxidant effects were found for thiazolyl-carbonyl-thiosemicarbazides Th-1-8 (Group A) and for the oxadiazoles Th-18-20 and thiadiazole Th-21 (Group C).Also, the virtual screening of thiazole derivatives revealed that all twenty-two compounds bind the active site of iNOS [21].Based on the virtual screening and on the obtained results, the activity may be due to their capacity to reduce NO synthesis by blocking the active site of iNOS [21].

Antioxidant and Anti-Inflammatory Effects of the Studied Compounds
The new thiazole compounds have anti-inflammatory effects by lowering the acute-phase response of bone marrow (L and PMN), phagocytic capacity (PI and PA) and oxidative stress [21].The best anti-inflammatory and antioxidant effects were found for thiazolyl-carbonyl-thiosemicarbazides Th-1-8 (Group A) and for the oxadiazoles Th-18-20 and thiadiazole Th-21 (Group C).Also, the virtual screening of thiazole derivatives revealed that all twenty-two compounds bind the active site of iNOS [21].Based on the virtual screening and on the obtained results, the activity may be due to their capacity to reduce NO synthesis by blocking the active site of iNOS [21].

Antioxidant and Anti-Inflammatory Effects of the Studied Compounds
The new thiazole compounds have anti-inflammatory effects by lowering the acute-phase response of bone marrow (L and PMN), phagocytic capacity (PI and PA) and oxidative stress [21].The best anti-inflammatory and antioxidant effects were found for thiazolyl-carbonylthiosemicarbazides Th-1-8 (Group A) and for the oxadiazoles Th-18-20 and thiadiazole Th-21 (Group C).Also, the virtual screening of thiazole derivatives revealed that all twenty-two compounds bind the active site of iNOS [21].Based on the virtual screening and on the obtained results, the activity may be due to their capacity to reduce NO synthesis by blocking the active site of iNOS [21].
In order to a model with less variables, we applied PCA on determined biological parameters (L, PMN, PA, PI, NO, TOS, TAR, OSI) (Table 2).After applying PCA only on these biological activity parameters, an acceptable model was obtained, the first three principal components containing 75.87% of total variance (Figure S3).However, it is interesting that the 3D scatterplot of the first three principal components revealed a pattern of clustering of the studied compounds, similar to the one of their lipophilicity (Figure 7).The three groups are well and distinctly grouped, with compound Th-22 as an exception.Here, groups A and B are better separated than in case of the lipophilicity (Figure 6), indicating their more distinct behavior in case of their biological activities.
Molecules 2015, 20, page-page 8 In order to obtain a model with less variables, we applied PCA on determined biological parameters (L, PMN, PA, PI, NO, TOS, TAR, OSI) (Table 2).After applying PCA only on these biological activity parameters, an acceptable model was obtained, the first three principal components containing 75.87% of total variance (Figure S3).However, it is interesting that the 3D scatterplot of the first three principal components revealed a pattern of clustering of the studied compounds, similar to the one of their lipophilicity (Figure 7).The three groups are well and distinctly grouped, with compound Th-22 as an exception.Here, groups A and B are better separated than in case of the lipophilicity (Figure 6), indicating their more distinct behavior in case of their biological activities.

Comparative Evaluation Both Lipophilicity Indices and Antioxidant and Anti-Inflammatory Activities Evaluation for the Studied Compounds
A similar clustering of the compounds to the one seen in the case of their lipophilicity indicates that there are some possible correlations between their lipophilicity and their biological activity.In this regard, PCA was applied on all the studied variables, describing lipophilicity and the anti-inflammatory and antioxidant activities.After applying PCA, an acceptable model was obtained-the first three principal components containing 74.96% of total variance (Figure S4).Analyzing the loadings of the variables and the matrix correlation of the original variables, we observed that NO parameter is positively correlated with the lipophilicity, and TAR activities are negatively correlated with the lipophilicity indices (Figure 8, Table 3).Molecules 2015, 20, page-page 9

Comparative Evaluation of Both Lipophilicity Indices and Antioxidant and Anti-Inflammatory Activities Evaluation for the Studied Compounds
A similar clustering of the compounds to the one seen in the case of their lipophilicity indicates that there are some possible correlations between their lipophilicity and their biological activity.In this regard, PCA was applied on all the studied variables, describing lipophilicity and the anti-inflammatory and antioxidant activities.After applying PCA, an acceptable model was obtained-the first three principal components containing 74.96% of total variance (Figure S4).Analyzing the loadings of the variables and the matrix correlation of the original variables, we observed that NO parameter is positively correlated with the lipophilicity, and TAR activities are negatively correlated with the lipophilicity indices (Figure 8, Table 3).Besides this, PA, TOS and OSI have statistically no correlation, while PMN and L have low correlation with any lipophilicity indices.However, the two previously mentioned variables (TAR and NO) have distinctly higher correlation with all lipophilicity indices.Activation and stimulation of phagocytes in the inflammatory process is dependent on membrane-derived lipid mediators (e.g., arachidonic acid derivatives) [23,24].This may explain the positive correlation between the reactive nitrogen intermediates (RNIs) (i.e., NO), produced by phagocytes, and lipophilicity (Figures 8 and 9).On the other hand, as a consequence of oxidative stress increase by reactive oxygen intermediates (ROIs) and RNIs, the antioxidant defense capacity is reduced.This justifies the negative correlation between TAR and lipophilicity (Figures 8 and 9).
Molecules 2015, 20, page-page and NO) have distinctly higher correlation with all lipophilicity indices.Activation and stimulation of phagocytes in the inflammatory process is dependent on membrane-derived lipid mediators (e.g., arachidonic acid derivatives) [23,24].This may explain the positive correlation between the reactive nitrogen intermediates (RNIs) (i.e., NO), produced by phagocytes, and lipophilicity (Figures 8 and 9).On the other hand, as a consequence of oxidative stress increase by reactive oxygen intermediates (ROIs) and RNIs, the antioxidant defense capacity is reduced.This justifies the negative correlation between TAR and lipophilicity (Figures 8 and 9).The 3D scatterplot of the scores for the first three principal components reported in Figure 10 reveals an even better and more compact clustering for the studied compounds forming the three distinct groups: the thiosemicarbazides Th-1-8 from group A with the best anti-inflammatory and antioxidant effects [21], the triazoles Th-9-17 from group B and the oxadiazoles and thiadiazoles Th-18-22 from group C. The knowledge of the grouping patterns of the tested variables allows for the reduction of experimental determined parameters, the variables which are close to each other in PCA results describing similar properties of the substances.For this reason, it is not necessary to measure and evaluate all variables to achieve the same characterization-thus reducing the amount of experimental work required [25].The 3D scatterplot of the scores for the first three principal components reported in Figure 10 reveals an even better and more compact clustering for the studied compounds forming the three distinct groups: the thiosemicarbazides Th-1-8 from group A with the best anti-inflammatory and antioxidant effects [21], the triazoles Th-9-17 from group B and the oxadiazoles and thiadiazoles Th-18-22 from group C.
Molecules 2015, 20, page-page and NO) have distinctly higher correlation with all lipophilicity indices.Activation and stimulation of phagocytes in the inflammatory process is dependent on membrane-derived lipid mediators (e.g., arachidonic acid derivatives) [23,24].This may explain the positive correlation between the reactive nitrogen intermediates (RNIs) (i.e., NO), produced by phagocytes, and lipophilicity (Figures 8 and 9).On the other hand, as a consequence of oxidative stress increase by reactive oxygen intermediates (ROIs) and RNIs, the antioxidant defense capacity is reduced.This justifies the negative correlation between TAR and lipophilicity (Figures 8 and 9).The 3D scatterplot of the scores for the first three principal components reported in Figure 10 reveals an even better and more compact clustering for the studied compounds forming the three distinct groups: the thiosemicarbazides Th-1-8 from group A with the best anti-inflammatory and antioxidant effects [21], the triazoles Th-9-17 from group B and the oxadiazoles and thiadiazoles Th-18-22 from group C. The knowledge of the grouping patterns of the tested variables allows for the reduction of experimental determined parameters, the variables which are close to each other in PCA results describing similar properties of the substances.For this reason, it is not necessary to measure and evaluate all variables to achieve the same characterization-thus reducing the amount of experimental work required [25].The knowledge of grouping patterns of the tested variables allows for the reduction of experimental determined parameters, the variables which are close to each other in PCA results describing similar properties of the substances.For this reason, it is not necessary to measure and evaluate all variables to achieve the same characterization-thus reducing the amount of experimental work required [25].
The pattern of the profiles of both NO and TAR vs. PC1 and R M0 respectively, clearly sustain the positive and negative correlation (Figure 9).For the first case (NO vs. PC1), the A and C groups present a distinctly higher similarity of the pattern than the B group (Figure 9 left), and the NO synthesis increases with the lipophilicity in agreement with the previous observations.Thus, best anti-inflammatory compounds Th-1-8 have the medium values for the lipophilicity parameters.For the second case, the negative correlation profiles between TAR and R M0 (Figure 9 right) is salient for group A and for some compounds from groups B and C. The negative correlation might be explained by the fact that the increase in hydrophobicity of a given compound may imply an increase in redox potential, thus decreasing its antioxidant ability [26].

Chromatographic Procedure
Chromatography was performed on 20 ˆ20 cm RP-18F 254s TLC precoated silica plates (Merck; Darmstadt, Germany).Solutions (1 mg/mL) of the tested compounds were prepared in iso-propanol, and 3 µL aliquots were spotted in duplicate on the plates by hand, 10 mm from the bottom edge and 20 mm apart.The mobile phases were composed of the iso-propanol-water binary mixtures, with a varying content of organic modifier between 45% and 65% (v/v) in 5% increments, as the study compounds differed considerably in their retention.Chromatography was performed in a normal developing chamber at room temperature, the developing distance being 10 cm.The chromatography chamber was saturated with the mobile phase for 30 min before use.After the development (30-60 min), the plates were air dried at room temperature and examined under a UV lamp (λ = 254 nm), and the R F (retardation factor) values were measured manually by a digital caliper.The experiments were performed in triplicate.

Prediction of Partitioning and Distribution Coefficients
The lipophilicity effects of compounds were predicted as log P and log D using the partitioning module of MarvinSketch 5.5.0.1, 2011 (ChemAxon Kft., Budapest, Hungary, http://www.chemaxon.com).Calculations for log D were estimated for the implicit reference values: log D 1.5 (log D calculated at the physiological pH of stomach), log D 5.0 (log D calculated at the physiological pH of duodenum), log D 6.5 (log D calculated at the physiological pH of jejuno-ileum) and log D 7.4 (log D calculated at the physiological pH of blood serum).Both calculations for log P and log D were performed using four different methods: a calculation algorithm developed by ChemAxon Kft. based on the publication of Viswanadhan and coworkers [27], one based on method developed by Klopman and coworkers [28], another one using the log P data from PHYSPROP© database and an equally weighted calculation algorithm combining the three methods.Supplementary, the software was set to take in account, for each set of calculations/algorithm, the major tautomeric forms of each compound.

Biological Activity Measurements
The synthesized thiazolyl-carbonyl-thiosemicarbazides and thiazolyl-azoles Th-1-Th-22 were previously evaluated in vivo for their anti-inflammatory activity in a turpentine oil-induced inflammation model [21].Their anti-inflammatory activity was assessed by evaluating the acute-phase response of bone marrow by leukocytes count (L) and by differential leucocyte count, expressed as a percentage (PMN) and by phagocytes' capacity (PI%-phagocytosis index,

Figure 2 .
Figure 2. Profiles of RM values at all mobile binary phase composition (% of organic modifier) for all studied compounds.

Figure 2 .
Figure 2. Profiles of R M values at all mobile binary phase composition (% of organic modifier) for all studied compounds.

Figure 3 .
Figure 3. Relationship between computed LogP and LogD (calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) values profiles and experimental lipophilicity indices RM0 and PC1.The L or R information in the brackets refers to either left or right disposition of the values on the y axis.

Figure 3 .
Figure 3. Relationship between computed LogP and LogD (calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) values profiles and experimental lipophilicity indices R M0 and PC1.The L or R information in the brackets refers to either left or right disposition of the values on the y axis.

Figure 4 .
Figure 4. Loadings (circle correlation) for the experimental lipophilicity indices (RM0, PC1, b) and computed (logD and logP-calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) lipophilicity indices for all the studied compounds.

Figure 4 .
Figure 4. Loadings (circle correlation) for the experimental lipophilicity indices (R M0 , PC1, b) and computed (logD and logP-calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) lipophilicity indices for all the studied compounds.

Figure 5 .
Figure 5. Loadings (circle correlation) for the experimental lipophilicity indices (RM0, PC1, b) and computed (LogD and LogP-calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) lipophilicity indices for compounds in group A (left) and for compounds in groups B and C (right).

Figure 6 .
Figure 6.3D plot of the scores for first three principal components after applying PCA on all the discussed parameters, describing lipophilicity for all the studied compounds.

Figure 5 .
Figure 5. Loadings (circle correlation) for the experimental lipophilicity indices (R M0 , PC1, b) and computed (LogD and LogP-calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) lipophilicity indices for compounds in group A (left) and for compounds in groups B and C (right).

Molecules 2015 , 7 Figure 5 .
Figure 5. Loadings (circle correlation) for the experimental lipophilicity indices (RM0, PC1, b) and computed (LogD and LogP-calculated by the equally weighted calculation algorithm combining the three methods described in the Experimental Section) lipophilicity indices for compounds in group A (left) and for compounds in groups B and C (right).

Figure 6 .
Figure 6.3D plot of the scores for first three principal components after applying PCA on all the discussed parameters, describing lipophilicity for all the studied compounds.

Figure 6 .
Figure 6.3D plot of the scores for first three principal components after applying PCA on all the discussed parameters, describing lipophilicity for all the studied compounds.

Figure 7 .
Figure 7. 3D scatterplot of the scores for the first three principal components after applying PCA on all the discussed parameters, describing biological activity for all the studied compounds.

Figure 7 .
Figure 7. 3D scatterplot of the scores for the first three principal components after applying PCA on all the discussed parameters, describing biological activity for all the studied compounds.

Figure 9 .
Figure 9. Profiles of NO vs. PC1 (left panel) and TAR vs. RM0 (right panel) for all studied compounds.

Figure 10 .
Figure 10.3D scatterplot of the scores for the first three principal components after applying PCA on all the discussed parameters describing biological activity and lipophilicity for all studied compounds.

Figure 9 .
Figure 9. Profiles of NO vs. PC1 (left panel) and TAR vs. R M0 (right panel) for all studied compounds.

Figure 9 .
Figure 9. Profiles of NO vs. PC1 (left panel) and TAR vs. RM0 (right panel) for all studied compounds.

Figure 10 .
Figure 10.3D scatterplot of the scores for the first three principal components after applying PCA on all the discussed parameters describing biological activity and lipophilicity for all studied compounds.

Figure 10 .
Figure 10.3D scatterplot of the scores for the first three principal components after applying PCA on all the discussed parameters describing biological activity and lipophilicity for all studied compounds.

Table 1 .
Retention parameters obtained from equation R M = R M0 + bC.

Table 2 .
Effect of thiazolyl-carbonyl-thiosemicarbazides and thiazolyl-azoles on the determined biological parameters for anti-inflammatory and antioxidant activities.

Table 2 .
Effect of thiazolyl-carbonyl-thiosemicarbazides and thiazolyl-azoles on the determined biological parameters for anti-inflammatory and antioxidant activities.

Table 3 .
Correlation matrix of the lipophilicity indices and antioxidant and anti-inflammatory effects for all the studied compounds.Values greater than 0.5 are in bold face characters.

Table 3 .
Correlation matrix of the lipophilicity indices and antioxidant and anti-inflammatory effects for all the studied compounds.Values greater than 0.5 are in bold face characters.