Quantitative Structure-Activity Relationship Study of Some Antipsychotics by Multiple Linear Regressions ()
1. Introduction
Chromatography is a powerful technique for the measurement of physicochemical properties that are used as parameters to correlate biological activities with structures in quantitative structure-activity relationship (QSAR) studies. QSAR modeling pertains to the construction of predictive models of biological activities as a function of structural and molecular information of compound library. These studies have had a tremendous impact in the fields of drug design, toxicology, and environmental monitoring. Typical molecular descriptors that are correlated to electronic properties, hydrophobicity, steric effect and topology, can be determined empirically through experimentation or theoretically via computational chemistry. Computationally determined property parameters have become crucial in identifying potential drug candidates [1] -[4] .
Psychosis is a syndrome, which is part of a group of very serious mental disorders in which some loss of contact with reality has occurred. The symptoms of psychosis may be lessened by antipsychotic medications. Antipsychotics affect many neurotransmitter systems. Common causes of chronic psychosis in the elderly include: dementia, depression, delirium, Parkinson’s disease, manic depressive illness, and schizophrenia [5] -[7] . Schizophrenia is the overwhelming mental disorder characterized by severe distortions of reality and disturbances in perception, intellectual performance, behaviour and motor activities. Various typical antipsychotics like chlorpromazine, haloperidol have been introduced which showed improvement in positive symptoms of schizophrenia by blocking dopaminergic transmission in the brain [8] . It is very important to control the antipsychotics distribution to the place of action. The lipophilicity has a significant impact on the absorption, distribution, metabolism, and excretion of compounds (ADME properties). Antipsychotics are targeting the central nervous system, so they must have certain lipophilicity to be able to pass the blood brain barrier by P-glycoprotein [9] [10] . Also, the gastrointestinal resorption and distribution of drugs through the bloodstream, by means of albumin, is dependent on lipophilicity [11] . The general rule is that the more lipid soluble a molecule or drug is, the more readily it will tend to enter the brain tissue [10] .
In preclinical trials the determination of lipophilicity of potential drug (through either experimental measurement or prediction) is one of the first selection criteria. The classical procedure is the determination of lipophilicity in terms of log P (the partition coefficient), which is descriptor of the differential partitioning of a neutral compound between two immiscible solvents, usually octan-1-ol and water [12] . Instead of this traditional shakeflask method, reversed-phase thin layer chromatography (RP TLC) is frequently used to estimate the lipophilicity of observed compounds due to the fact that the same basic intermolecular interactions determine the behavior of compound in both biological and chromatographic environment [13] -[17] . For this purpose the intercepts of the linear relationships between the logarithm of retention constants and the volume fraction of the organic modifier in a binary mobile phase obtained in RP TLC experiments are the most suitable. The relation is given by equation:
(1)
where j stands for the concentration of the organic component in the mobile phase and m is the slope, which indicates the rate at which the solubility of the solute in the mobile phase increases with changes in its composition. Lipophilicity, measured as
, represents the relative affinity of different compounds for the non-aqueous environment in the biological system.
is the value of RM in pure water, and therefore it could reflect the dependence of the hydrophobic properties of investigated compound on its structure.
Multivariate regression models in chemistry and other sciences quantitatively relate a response (dependent) variable y to a block of predictor variables X, in the form of a mathematical equation y = f(X), where the predictors can be determined experimentally or computationally. Among the best known of such quantitative-X-y relationships is quantitative structure-activity relationships, in which y is a biological response and any of the predictors, designated as descriptors, may account for a microscopic (i.e., determined by molecular structure) or a macroscopic property. QSAR has become important in medicinal chemistry, pharmacy, toxicology and environmental science, because it deals with bioactive substances such as drugs and toxicants [18] .
The objectives of this work were to investigate the retention behavior of this kind of compound using reversed-phase thin-layer chromatography, to determine their lipophilicity, to investigate the effect of different mobile-phase modifiers on the retention as well as on the lipophilicity parameter, and to estimate the quantitative structure-activity relationship using molecular descriptors.
2. Experimental
2.1. Chromatography
Chromatographic investigations were carried out by horizontal thin layer chromatography on silica gel RP 18 plates, 10 × 10 cm (Merck, Darmstadt, Germany) using a Camag horizontal HPTLC development chamber in the tank configuration. Standard solutions (5 mg/mL) of the compounds were prepared in chloroform. The plates were spotted with 1.0 μL aliquots of freshly prepared solutions of the corresponding compound. Before development, the spotted plates were equilibrated for 30 min in a chromatographic chamber saturated with mobile phase vapor. All solvents used throughout the present study were of analytical-grade purity. The applied mobile phase was a mixture of different mobile-phase modifiers (acetone, tetrahydrofuran, methanol) and water. The concentration of organic modifier in the mobile phase ranged from 50% to 80% (v/v) in 5% increments. The investigated compounds were chromatographed simultaneously. After development, the spots were colored by Dragendorff reagent. RF values were determined as averages from three independent measurements. All measurements were carried out at ambient temperature (22˚C ± 2˚C).
2.2. Calculations
All structures were drawn with the HyperChem Professional software (version 7.0, Hybercube, Gainseville, FL, USA). In order to obtain molecular descriptors; the geometry optimization of molecules was performed by the molecular mechanics MM+ force field method. The single point calculation was done with the semi-empirical quantum chemical method AM1. Polak-Ribiere algorithm with the convergence limit set at 0.1 kcal/mol was used during modeling process. The HyperChem was used to calculate Total energy (ETotal), Dipole moment (µ), Energy of the highest occupied molecular orbital (EHOMO), Energy of the lowest unoccupied molecular orbital (ELUMO), Surface area (SA), Hydration energy (Ehydratation), Refractivity and Polarizability. Molar volume, Molar depth, Hydrophilic lipophilic balance (HLB), Solubility parameter (Sp), Hansen dispersion, Hansen polarity, Hydrogen bond acceptor (HBAcc), Hydrogen bond donor (HBD), Parachor, Kappa 2, Total e state, Water solubility, Connectivity indices, Valence indices, logP Crippen, Surface tension in water (STW), Hydrogen bond number (HBN), Hydrophilic surface area (HSA) and Polar surface area (PSA) were calculated by MMP Plus [19] . ChemDraw 8.0 was used to calculate logP. Statistical calculations, variable selection routine and multiple linear regression analysis (MLR) were performed by NCSS 2004 software package [20] . We carried out multivariate variable selection to select descriptors for MLR. The reduced collection of descriptors was used as the input for final MLR analysis. The quality criteria of the fit in MLR analysis were squared correlation coefficients (r2), cross-validated coefficient (
), the mean square error (MSE) and Fischer significance value (F). The prediction performance was validated using a “leave-one-out” cross validation method. The cross-validated
values reflect the overall predictive ability of the model defined as (SSY-PRESS)/SSY. PRESS is predicted residual error sum of squares and SSY is the sum of the squared deviations of the dependent variable values from their mean. The significance level of the performed calculations was above 95%. RM values were calculated by use of the Bate-Smith and Westall equation [21] .
3. Results and Discussion
Structures of the investigated antipsychotics are shown in Figure 1 while their names, IUPAC names and molecular formula are listed in Table 1.
The retention parameters (RF and RM) of observed antipsychotics were determined at several compositions of the three different binary solvent systems composed of organic modifier and water: methanol-water, acetonewater and tetrahydrofuran-water. The investigated compounds, regardless on differences in their structure, showed some regularity in chromatographic behavior under applied chromatographic condition. Hydrophobic interactions dominated in reversed-phase thin layer chromatography. The most nonpolar compound has the highest retention in all applied mobile-phases. Also, regular retention behavior was observed, i.e. retention decreased regularly with increasing concentration of organic modifier in the mobile phase.
The linear relationship between the RM values and the concentration of the organic modifier in the mobile phase (j), expressed by Equation (1).
is lipophilicity parameter and m is the slope and represents the hydrophobic surface area. Lipophilicity, measured as
, represents the relative affinity of different compounds for the non-aqueous environment in the biological system.
is the value of RM in pure water, and therefore it could reflect the dependence of the hydrophobic properties of investigated compound on its structure.
The slope (m) and intercept
values, and the statistical data (correlation coefficient (r) and standard deviation (s)) for each binary system are listed in Table 2.
The chromatographically obtained lipophilicity parameter,
, is in agreement with structures of the investigated antipsychotics. The results from Table 2 show difference between the
of investigated compounds. For all investigated compounds, the
values are highest in methanol, which has the lowest elution strength among all the organic modifiers applied. Taking into account the observed retention, it can be concluded that compounds 7 and 4 exhibit stronger retention compared to 1 and 3. Retention of the compounds 7 and 4 is increased due to delocalization of electrons of substituents along the rings. Based on RM values (Table 2) the tested compounds are in accordance with the increasing lipophilicity (for the methanol and acetone as organic modifiers) described as follows: 1 < 3 < 2 < 5 < 6 < 4 < 7. This sequence is different when tetrahydrofuran applied as organic modifier. Also, compound 1 has a higher mobility and lower
values, in all applied mobile
Figure 1. Structures of the investigated antipsychotics.
Table 1. IUPAC names, and molecular formula of investigated antypsichotics.
Table 2. Lipophilicity parameter and regression data.
phases, as most polar compound. Lipophilicity parameter of compound 1 is lowest, when as organic modifier was applied tetrahydrofuran, which is in accordance with polarity of this organic solvent (4.0) in comparison with methanol and acetone (5.1) [22] . High correlation was obtained between the intercept
and the slopes, m, values (Table 3).
The m values are specific hydrophobic surface of the observed compound while the
is lipophilicity parameter. High correlation coefficients indicate that the substances investigated could be regarded as a homologous series [23] .
Also, the chromatographically determined lipophilicity parameter values with different organic modifier were correlated. The equations of these linear relationships are:



These differences are most probably a consequence of the different chemical natures of the three organic modifier (such as polarity as well as elution power) and indicates that the values basicaly reflect the some molecular propeties of the solute in the mobile phases used.
In order to established quantitative-structure models, we used the calculated molecular descriptors as a quantitative measure of individual structural features. Due to the collinearity problem in MLR analysis, the collinear descriptors (r > 0.9) can be removed before model development, by applying the heuristic method [24] . A correlation check for the descriptors was performed. Established QSAR models for different organic modifiers in mobile phase are shown in Table 4.
This is in accordance with experimentally obtained retention parameters. Based on statistically parameters we checked the robustness of established QSAR models. The statistically best model was obtained when organic modifier was acetone (model 1, Table 4). It is well known that the smaller PRESS is, the better the predictability of the model [25] . If PRESS is smaller than SSY the model predict is better than chance and can be considered to be statistically significant. In a reasonable QSAR model, PRESS/SSY should be smaller than 0.4 (model 1 and 2, Table 4). Obtained value of PRESS/SSY ratio is smaller than 0.1 indicates an excellent established model (model 1, Table 4) [26] . Also, QSAR model 1 shows a high degree of correlation between experimentally obtained and predicted parameters of retention, RM (Figure 2).
Established model 1 gives us the possibility to gain insight into factors responsible for the retention as well as lipophilicity of the investigated set of the compounds. The results obtained indicate that the most relevant descriptors influencing lipophilicity parameters are: Energy of the highest occupied molecular orbital, Energy of
Table 3. Correlation coefficients for correlation between the intercept
and the slopes (m) values.
Table 4. Regression QSAR models obtained for three different mobile phase modifiers.
Sp-Solubility parameter; EHOMO-Energy of the highest occupied molecular orbital; ELUMO-Energy of the lowest unoccupied molecular orbital; HSAHydrophilic surface area; r2-correlation coefficient; F-Fischer significance value; MSE-mean square error;
-cross-validated coefficient; PRESSpredicted residual error sum of squares; SSY-sum of the squared deviations.
Figure 2. Predicted
estimated by MLR versus experimental
(model 1).
the lowest unoccupied molecular orbital, and Solubility parameter. These descriptors can be accounted for the structural features responsible for chromatographic behavior as well as chromatographically determined lipophilicity parameter of investigated compounds. Lipophilicity parameter of investigated compounds increases with increasing EHOMO and ELUMO descriptors and decrease with increasing of solubility parameter. The order of significance of the descriptors is: ELUMO > EHOMO > Sp. It is obviously that ELUMO has the highest significance impact in comparison with EHOMO and Sp. The ELUMO is the property of electronic structure and represents the electron affinity of a molecule or its reactivity as an electrophile. Good electrophiles are those were the ELUMO is “low-lying” [27] . The presence of ELUMO in this model suggests that the higher this energy is, the weaker are the interactions of the compounds with mobile phases, and also higher values of
were observed. This is in accordance with chromatographic behavior of these compounds. A higher EHOMO suggests higher affinity of the molecule to react as a nucleophile. Therefore, the compounds with high values of EHOMO shows the highest retention under applied reversed-phase chromatographic conditions as well as highest
as can be seen from established model 1 (Table 4). Solubility parameter indicate that the stronger the intermolecular interactions between molecules and mobile phase are, the analytes are less retained on the stationary phase and lower RM as well as
are obtained.
Therefore, from Table 4 it is evident that the same molecular descriptors describe lipophilicity in models 1 and 2 regardless on applied organic modifiers. This is probably a consequence of the different eluotropic strength of applied organic modifiers. Model 2 is less statistically significant than model 1, probably due to ability of molecules of methanol to form monolayer at the surface of stationary phase [28] . The influence of applied organic modifier is obvious when tetrahydrofuran is in the mobile phase (model 3, Table 4). In this model the opposite influence of Solubility parameter was obtained and Hydrophilic surface area occurs as relevant molecular descriptor. This is in accordance with mechanisms in reversed-phase thin-layer chromatography considering properties of applied organic modifiers in eluotropic series based on polarity parameters [22] .
4. Conclusion
Quantitative structure-activity relationship studies have been performed to correlate the molecular characteristics of observed compounds with their retention as well as with their chromatographically determinated lipophilicity parameters. RP TLC proved to be a reliable and accurate method of describing the lipophilic nature of observed antipsychotics. Obtained results are in agreement with polarity of applied organic modifiers as well as structure of investigated compounds. In addition, we used molecular descriptors to establish QSAR models for all applied mobile phases. Established QSAR model for acetone as organic modifier (model 1) is excellent. For this model PRESS/SSY value is smaller than 0.1. This study provides good information about pharmacologically important physico-chemical parameters of observed antipsychotics relevant to variations in molecular lipophilicity and chromatographic behavior. A very advantageous feature of established models is that it allows us to understand chromatographic behavior of novel, not yet synthesized compounds, solely from their structural descriptors, and to estimate lipophilicity for similar compounds. These QSAR models could help multi-target novel antipsychotic compound design.
Acknowledgements
This work was performed within the framework of the research project No 172017 supported by the Ministry of Education, Science and Technological development of Serbia.